Integrative analysis of 111 reference human epigenomes.

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

Point Mutations in Exon 1B of APC Reveal Gastric Adenocarcinoma and Proximal Polyposis of the Stomach as a Familial Adenomatous Polyposis Variant.

Gastric adenocarcinoma and proximal polyposis of the stomach (GAPPS) is an autosomal-dominant cancer-predisposition syndrome with a significant risk of gastric, but not colorectal, adenocarcinoma. We mapped the gene to 5q22 and found loss of the wild-type allele on 5q in fundic gland polyps from affected individuals. Whole-exome and -genome sequencing failed to find causal mutations but, through Sanger sequencing, we identified point mutations in APC promoter 1B that co-segregated with disease in all six families. The mutations reduced binding of the YY1 transcription factor and impaired activity of the APC promoter 1B in luciferase assays. Analysis of blood and saliva from carriers showed allelic imbalance of APC, suggesting that these mutations lead to decreased allele-specific expression in vivo. Similar mutations in APC promoter 1B occur in rare families with familial adenomatous polyposis (FAP). Promoter 1A is methylated in GAPPS and sporadic FGPs and in normal stomach, which suggests that 1B transcripts are more important than 1A in gastric mucosa. This might explain why all known GAPPS-affected families carry promoter 1B point mutations but only rare FAP-affected families carry similar mutations, the colonic cells usually being protected by the expression of the 1A isoform. Gastric polyposis and cancer have been previously described in some FAP-affected individuals with large deletions around promoter 1B. Our finding that GAPPS is caused by point mutations in the same promoter suggests that families with mutations affecting the promoter 1B are at risk of gastric adenocarcinoma, regardless of whether or not colorectal polyps are present.

DMRFusion: A differentially methylated region detection tool based on the ranked fusion method.

DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Computational analysis of differentially methylated regions (DMRs) could explore the underlying reasons of methylation. DMRFusion is presented as a useful tool for comprehensive DNA methylation analysis of DMRs on methylation sequencing data. This tool is designed base on the integration of several ranking methods; Information gain, Between versus within Class scatter ratio, Fisher ratio, Z-score and Welch's t-test. In this study, DMRFusion on reduced representation bisulfite sequencing (RRBS) data in chronic lymphocytic leukemia cancer displayed 30 nominated regions and CpG sites with a maximum methylation difference detected in the hypermethylation DMRs. We realized that DMRFusion is able to process methylation sequencing data in an efficient and accurate manner and to provide annotation and visualization for DMRs with high fold difference score (p-value and FDR<0.05 and type I error: 0.04).

Micro-ribonucleic acid-155 is a direct target of Meis1, but not a driver in acute myeloid leukemia.

Micro-ribonucleic acid-155 (miR-155) is one of the first described oncogenic miRNAs. Although multiple direct targets of miR-155 have been identified, it is not clear how it contributes to the pathogenesis of acute myeloid leukemia. We found miR-155 to be a direct target of Meis1 in murine Hoxa9/Meis1 induced acute myeloid leukemia. The additional overexpression of miR-155 accelerated the formation of acute myeloid leukemia in Hoxa9 as well as in Hoxa9/Meis1 cells in vivo However, in the absence or following the removal of miR-155, leukemia onset and progression were unaffected. Although miR-155 accelerated growth and homing in addition to impairing differentiation, our data underscore the pathophysiological relevance of miR-155 as an accelerator rather than a driver of leukemogenesis. This further highlights the complexity of the oncogenic program of Meis1 to compensate for the loss of a potent oncogene such as miR-155. These findings are highly relevant to current and developing approaches for targeting miR-155 in acute myeloid leukemia.

The clonal and mutational evolution spectrum of primary triple-negative breast cancers.

Primary triple-negative breast cancers (TNBCs), a tumour type defined by lack of oestrogen receptor, progesterone receptor and ERBB2 gene amplification, represent approximately 16% of all breast cancers. Here we show in 104 TNBC cases that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some having only a handful of coding somatic aberrations in a few pathways, whereas others contain hundreds of coding somatic mutations. High-throughput RNA sequencing (RNA-seq) revealed that only approximately 36% of mutations are expressed. Using deep re-sequencing measurements of allelic abundance for 2,414 somatic mutations, we determine for the first time-to our knowledge-in an epithelial tumour subtype, the relative abundance of clonal frequencies among cases representative of the population. We show that TNBCs vary widely in their clonal frequencies at the time of diagnosis, with the basal subtype of TNBC showing more variation than non-basal TNBC. Although p53 (also known as TP53), PIK3CA and PTEN somatic mutations seem to be clonally dominant compared to other genes, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal, cell shape and motility proteins occurred at lower clonal frequencies, suggesting that they occurred later during tumour progression. Taken together, our results show that understanding the biology and therapeutic responses of patients with TNBC will require the determination of individual tumour clonal genotypes.

Novel candidate genes may be possible predisposing factors revealed by whole exome sequencing in familial esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma is one of the deadliest of all the cancers. Its metastatic properties portend poor prognosis and high rate of recurrence. A more advanced method to identify new molecular biomarkers predicting disease prognosis can be whole exome sequencing. Here, we report the most effective genetic variants of the Notch signaling pathway in esophageal squamous cell carcinoma susceptibility by whole exome sequencing. We analyzed nine probands in unrelated familial esophageal squamous cell carcinoma pedigrees to identify candidate genes. Genomic DNA was extracted and whole exome sequencing performed to generate information about genetic variants in the coding regions. Bioinformatics software applications were utilized to exploit statistical algorithms to demonstrate protein structure and variants conservation. Polymorphic regions were excluded by false-positive investigations. Gene-gene interactions were analyzed for Notch signaling pathway candidates. We identified novel and damaging variants of the Notch signaling pathway through extensive pathway-oriented filtering and functional predictions, which led to the study of 27 candidate novel mutations in all nine patients. Detection of the trinucleotide repeat containing 6B gene mutation (a slice site alteration) in five of the nine probands, but not in any of the healthy samples, suggested that it may be a susceptibility factor for familial esophageal squamous cell carcinoma. Noticeably, 8 of 27 novel candidate gene mutations (e.g. epidermal growth factor, signal transducer and activator of transcription 3, MET) act in a cascade leading to cell survival and proliferation. Our results suggest that the trinucleotide repeat containing 6B mutation may be a candidate predisposing gene in esophageal squamous cell carcinoma. In addition, some of the Notch signaling pathway genetic mutations may act as key contributors to esophageal squamous cell carcinoma.

Recurrent somatic DICER1 mutations in nonepithelial ovarian cancers.

BACKGROUND: Germline truncating mutations in DICER1, an endoribonuclease in the RNase III family that is essential for processing microRNAs, have been observed in families with the pleuropulmonary blastoma-family tumor and dysplasia syndrome. Mutation carriers are at risk for nonepithelial ovarian tumors, notably sex cord-stromal tumors. METHODS: We sequenced the whole transcriptomes or exomes of 14 nonepithelial ovarian tumors and noted closely clustered mutations in the region of DICER1 encoding the RNase IIIb domain of DICER1 in four samples. We then sequenced this region of DICER1 in additional ovarian tumors and in certain other tumors and queried the effect of the mutations on the enzymatic activity of DICER1 using in vitro RNA cleavage assays. RESULTS: DICER1 mutations in the RNase IIIb domain were found in 30 of 102 nonepithelial ovarian tumors (29%), predominantly in Sertoli-Leydig cell tumors (26 of 43, or 60%), including 4 tumors with additional germline DICER1 mutations. These mutations were restricted to codons encoding metal-binding sites within the RNase IIIb catalytic centers, which are critical for microRNA interaction and cleavage, and were somatic in all 16 samples in which germline DNA was available for testing. We also detected mutations in 1 of 14 nonseminomatous testicular germ-cell tumors, in 2 of 5 embryonal rhabdomyosarcomas, and in 1 of 266 epithelial ovarian and endometrial carcinomas. The mutant DICER1 proteins had reduced RNase IIIb activity but retained RNase IIIa activity. CONCLUSIONS: Somatic missense mutations affecting the RNase IIIb domain of DICER1 are common in nonepithelial ovarian tumors. These mutations do not obliterate DICER1 function but alter it in specific cell types, a novel mechanism through which perturbation of microRNA processing may be oncogenic. (Funded by the Terry Fox Research Institute and others.).

ALEA: a toolbox for allele-specific epigenomics analysis.

The assessment of expression and epigenomic status using sequencing based methods provides an unprecedented opportunity to identify and correlate allelic differences with epigenomic status. We present ALEA, a computational toolbox for allele-specific epigenomics analysis, which incorporates allelic variation data within existing resources, allowing for the identification of significant associations between epigenetic modifications and specific allelic variants in human and mouse cells. ALEA provides a customizable pipeline of command line tools for allele-specific analysis of next-generation sequencing data (ChIP-seq, RNA-seq, etc.) that takes the raw sequencing data and produces separate allelic tracks ready to be viewed on genome browsers. The pipeline has been validated using human and hybrid mouse ChIP-seq and RNA-seq data. AVAILABILITY: The package, test data and usage instructions are available online at http://www.bcgsc.ca/platform/bioinfo/software/alea CONTACT: : mkarimi1@interchange.ubc.ca or sjones@bcgsc.ca Supplementary information: Supplementary data are available at Bioinformatics online.

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.).

deFuse: an algorithm for gene fusion discovery in tumor RNA-Seq data.

Gene fusions created by somatic genomic rearrangements are known to play an important role in the onset and development of some cancers, such as lymphomas and sarcomas. RNA-Seq (whole transcriptome shotgun sequencing) is proving to be a useful tool for the discovery of novel gene fusions in cancer transcriptomes. However, algorithmic methods for the discovery of gene fusions using RNA-Seq data remain underdeveloped. We have developed deFuse, a novel computational method for fusion discovery in tumor RNA-Seq data. Unlike existing methods that use only unique best-hit alignments and consider only fusion boundaries at the ends of known exons, deFuse considers all alignments and all possible locations for fusion boundaries. As a result, deFuse is able to identify fusion sequences with demonstrably better sensitivity than previous approaches. To increase the specificity of our approach, we curated a list of 60 true positive and 61 true negative fusion sequences (as confirmed by RT-PCR), and have trained an adaboost classifier on 11 novel features of the sequence data. The resulting classifier has an estimated value of 0.91 for the area under the ROC curve. We have used deFuse to discover gene fusions in 40 ovarian tumor samples, one ovarian cancer cell line, and three sarcoma samples. We report herein the first gene fusions discovered in ovarian cancer. We conclude that gene fusions are not infrequent events in ovarian cancer and that these events have the potential to substantially alter the expression patterns of the genes involved; gene fusions should therefore be considered in efforts to comprehensively characterize the mutational profiles of ovarian cancer transcriptomes.

Using next-generation sequencing for the diagnosis of rare disorders: a family with retinitis pigmentosa and skeletal abnormalities.

Linkage analysis with subsequent candidate gene sequencing is typically used to diagnose novel inherited syndromes. It is now possible to expedite diagnosis through the sequencing of all coding regions of the genome (the exome) or full genomes. We sequenced the exomes of four members of a family presenting with spondylo-epiphyseal dysplasia and retinitis pigmentosa and identified a six-base-pair (6-bp) deletion in GNPTG, the gene implicated in mucolipidosis type IIIγ. The diagnosis was confirmed by biochemical studies and both broadens the mucolipidosis type III phenotype and demonstrates the clinical utility of next-generation sequencing to diagnose rare genetic diseases.

Subtype-specific mutation of PPP2R1A in endometrial and ovarian carcinomas.

PPP2R1A mutations have recently been described in 3/42 (7%) of clear cell carcinomas of the ovary. PPP2R1A encodes the α-isoform of the scaffolding subunit of the serine/threonine protein phosphatase 2A (PP2A) holoenzyme. This putative tumour suppressor complex is involved in growth and survival pathways. Through targeted sequencing of PPP2R1A, we identified somatic missense mutations in 40.8% (20/49) of high-grade serous endometrial tumours, and 5.0% (3/60) of endometrial endometrioid carcinomas. Mutations were also identified in ovarian tumours at lower frequencies: 12.2% (5/41) of endometrioid and 4.1% (2/49) of clear cell carcinomas. No mutations were found in 50 high-grade and 12 low-grade serous carcinomas. Amino acid residues affected by these mutations are highly conserved across species and are involved in direct interactions with regulatory B-subunits of the PP2A holoenzyme. PPP2R1A mutations in endometrial high-grade serous carcinomas are a frequent and potentially targetable feature of this disease. The finding of frequent PPP2R1A mutations in high-grade serous carcinoma of the endometrium but not in high-grade serous carcinoma of the ovary provides clear genetic evidence that these are distinct diseases.

Effect of mastitis during the first lactation on production and reproduction performance of Holstein cows.

The aim of this study was to evaluate the effect of postpartum mastitis between first calving and subsequent conception on production and reproduction performance as well as culling of Holstein cows. A data set of 9,183 first lactation cows was used. Results showed that the first cumulative 100 days' milk production and the milk yield standardized to 305 days were affected by the interval from calving to first mastitis (P < 0.05). Cows with one episode of mastitis produced more milk than those with repeated episodes of mastitis (P < 0.01). Increase in the number of mastitis episodes and also decrease in interval between first calving and mastitis increased services per conception (P < 0.001). Mastitis episode and the interval between calving and first mastitis had no apparent impact on the calving to conception interval (P > 0.05). Calving year, calving difficulty score, and cumulative first 60 days milk production had significant impacts on mastitis risk (P < 0.05). The interval from calving to the first incidence of mastitis decreased over the period studied (P < 0.001). Productive life tended to be decreased due to mastitis (P = 0.07). Survival analysis showed a significant difference between the lengths of productive life for cows with different intervals from calving to first mastitis (P < 0.01). The results demonstrated that clinical mastitis between first calving and conception reduced production and reproduction performance with an increase in chance of culling.

Epithelial-Mesenchymal Plasticity Induced by Discontinuous Exposure to TGFß1 Promotes Tumour Growth.

Transitions between epithelial and mesenchymal cellular states (EMT/MET) contribute to cancer progression. We hypothesize that EMT followed by MET promotes cell population heterogeneity, favouring tumour growth. We developed an EMT model by on and off exposure of epithelial EpH4 cells (E-cells) to TGFß1 that mimics phenotypic EMT (M-cells) and MET. We aimed at understanding whether phenotypic MET is accompanied by molecular and functional reversion back to epithelia by using RNA sequencing, immunofluorescence (IF), proliferation, wound healing, focus formation and mamosphere formation assays as well as cell xenografts in nude mice. Phenotypic reverted epithelial cells (RE-cells) obtained after MET induction presented epithelial morphologies and proliferation rates resembling E cells. However, the RE transcriptomic profile and IF staining of epithelial and mesenchymal markers revealed a uniquely heterogeneous mixture of cell subpopulations with a high self-renewal ability. RE cell heterogeneity was stably maintained for long periods after TGFß1 removal both in vitro and in large tumours derived from the nude mice. Overall, we show that phenotypic reverted epithelial cells (RE cells) do not return to the molecular and functional epithelial state and present mesenchymal features related to aggressiveness and cellular heterogeneity that favour tumour growth in vivo. This work strengthens epithelial cell reprogramming and cellular heterogeneity fostered by inflammatory cues as a tumour growth-promoting factor in vivo.

Single-cell analysis of RORα tracer mouse lung reveals ILC progenitors and effector ILC2 subsets.

Lung group 2 innate lymphoid cells (ILC2s) drive allergic inflammation and promote tissue repair. ILC2 development is dependent on the transcription factor retinoic acid receptor-related orphan receptor (RORα), which is also expressed in common ILC progenitors. To elucidate the developmental pathways of lung ILC2s, we generated RORα lineage tracer mice and performed single-cell RNA sequencing, flow cytometry, and functional analyses. In adult mouse lungs, we found an IL-18Rα+ST2- population different from conventional IL-18Rα-ST2+ ILC2s. The former was GATA-3intTcf7EGFP+Kit+, produced few cytokines, and differentiated into multiple ILC lineages in vivo and in vitro. In neonatal mouse lungs, three ILC populations were identified, namely an ILC progenitor population similar to that in adult lungs and two distinct effector ILC2 subsets that differentially produced type 2 cytokines and amphiregulin. Lung ILC progenitors might actively contribute to ILC-poiesis in neonatal and inflamed adult lungs. In addition, neonatal lung ILC2s include distinct proinflammatory and tissue-repairing subsets.

Epigenomic programming in early fetal brain development.

Aim: To provide a comprehensive understanding of gene regulatory networks in the developing human brain and a foundation for interpreting pathogenic deregulation. Materials & methods: We generated reference epigenomes and transcriptomes of dissected brain regions and primary neural progenitor cells (NPCs) derived from cortical and ganglionic eminence tissues of four normal human fetuses. Results: Integration of these data across developmental stages revealed a directional increase in active regulatory states, transcription factor activities and gene transcription with developmental stage. Consistent with differences in their biology, NPCs derived from cortical and ganglionic eminence regions contained common, region specific, and gestational week specific regulatory states. Conclusion: We provide a high-resolution regulatory network for NPCs from different brain regions as a comprehensive reference for future studies.

MicroRNA-708 is a novel regulator of the Hoxa9 program in myeloid cells.

MicroRNAs (miRNAs) are commonly deregulated in acute myeloid leukemia (AML), affecting critical genes not only through direct targeting, but also through modulation of downstream effectors. Homeobox (Hox) genes balance self-renewal, proliferation, cell death, and differentiation in many tissues and aberrant Hox gene expression can create a predisposition to leukemogenesis in hematopoietic cells. However, possible linkages between the regulatory pathways of Hox genes and miRNAs are not yet fully resolved. We identified miR-708 to be upregulated in Hoxa9/Meis1 AML inducing cell lines as well as in AML patients. We further showed Meis1 directly targeting miR-708 and modulating its expression through epigenetic transcriptional regulation. CRISPR/Cas9 mediated knockout of miR-708 in Hoxa9/Meis1 cells delayed disease onset in vivo, demonstrating for the first time a pro-leukemic contribution of miR-708 in this context. Overexpression of miR-708 however strongly impeded Hoxa9 mediated transformation and homing capacity in vivo through modulation of adhesion factors and induction of myeloid differentiation. Taken together, we reveal miR-708, a putative tumor suppressor miRNA and direct target of Meis1, as a potent antagonist of the Hoxa9 phenotype but an effector of transformation in Hoxa9/Meis1. This unexpected finding highlights the yet unexplored role of miRNAs as indirect regulators of the Hox program during normal and aberrant hematopoiesis.

Familial Esophageal Squamous Cell Carcinoma with damaging rare/germline mutations in KCNJ12/KCNJ18 and GPRIN2 genes.

In Iran, esophageal cancer is the fourth common cancers in women and sixth common cancers in men. Here we evaluated the importance of familial risk factors and the role of genetic predisposition in Esophageal Squamous Cell Carcinoma (ESCC) using Whole-Exome Sequencing (WES). Germline damaging mutations were identified in WES data from 9 probands of 9 unrelated ESCC pedigrees. Mutations were confirmed with Sanger sequencing and evaluated amplification-refractory mutation system-Polymerase Chain Reaction (ARMS-PCR) in 50 non-related ethnically matched samples and in complete genomics database. Sixteen candidate variants were detected in ESCC 9 probands. Four of these 16 variants were rare damaging mutations including novel mutations in KCNJ12/KCNJ18, and GPRIN2 genes. This WES study in Iranian patients with ESCC, provides insight into the identification of novel germline mutations in familial ESCC. Our data suggest an association between specific mutations and increased risk of ESCC.

Characterization of the human thyroid epigenome.

The thyroid gland, necessary for normal human growth and development, functions as an essential regulator of metabolism by the production and secretion of appropriate levels of thyroid hormone. However, assessment of abnormal thyroid function may be challenging suggesting a more fundamental understanding of normal function is needed. One way to characterize normal gland function is to study the epigenome and resulting transcriptome within its constituent cells. This study generates the first published reference epigenomes for human thyroid from four individuals using ChIP-seq and RNA-seq. We profiled six histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K36me3, H3K9me3, H3K27me3), identified chromatin states using a hidden Markov model, produced a novel quantitative metric for model selection and established epigenomic maps of 19 chromatin states. We found that epigenetic features characterizing promoters and transcription elongation tend to be more consistent than regions characterizing enhancers or Polycomb-repressed regions and that epigenetically active genes consistent across all epigenomes tend to have higher expression than those not marked as epigenetically active in all epigenomes. We also identified a set of 18 genes epigenetically active and consistently expressed in the thyroid that are likely highly relevant to thyroid function. Altogether, these epigenomes represent a powerful resource to develop a deeper understanding of the underlying molecular biology of thyroid function and provide contextual information of thyroid and human epigenomic data for comparison and integration into future studies.

Nucleosome Density ChIP-Seq Identifies Distinct Chromatin Modification Signatures Associated with MNase Accessibility.

Nucleosome position, density, and post-translational modification are widely accepted components of mechanisms regulating DNA transcription but still incompletely understood. We present a modified native ChIP-seq method combined with an analytical framework that allows MNase accessibility to be integrated with histone modification profiles. Application of this methodology to the primitive (CD34+) subset of normal human cord blood cells enabled genomic regions enriched in one versus two nucleosomes marked by histone 3 lysine 4 trimethylation (H3K4me3) and/or histone 3 lysine 27 trimethylation (H3K27me3) to be associated with their transcriptional and DNA methylation states. From this analysis, we defined four classes of promoter-specific profiles and demonstrated that a majority of bivalent marked promoters are heterogeneously marked at a single-cell level in this primitive cell type. Interestingly, extension of this approach to human embryonic stem cells revealed an altered relationship between chromatin modification state and nucleosome content at promoters, suggesting developmental stage-specific organization of histone methylation states.