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.

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.

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.

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.

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.

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.

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.

Analysis of Normal Human Mammary Epigenomes Reveals Cell-Specific Active Enhancer States and Associated Transcription Factor Networks.

The normal adult human mammary gland is a continuous bilayered epithelial system. Bipotent and myoepithelial progenitors are prominent and unique components of the outer (basal) layer. The inner (luminal) layer includes both luminal-restricted progenitors and a phenotypically separable fraction that lacks progenitor activity. We now report an epigenomic comparison of these three subsets with one another, with their associated stromal cells, and with three immortalized, non-tumorigenic human mammary cell lines. Each genome-wide analysis contains profiles for six histone marks, methylated DNA, and RNA transcripts. Analysis of these datasets shows that each cell type has unique features, primarily within genomic regulatory regions, and that the cell lines group together. Analyses of the promoter and enhancer profiles place the luminal progenitors in between the basal cells and the non-progenitor luminal subset. Integrative analysis reveals networks of subset-specific transcription factors.

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.

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.

Genome-Wide Profiles of Extra-cranial Malignant Rhabdoid Tumors Reveal Heterogeneity and Dysregulated Developmental Pathways.

Malignant rhabdoid tumors (MRTs) are rare lethal tumors of childhood that most commonly occur in the kidney and brain. MRTs are driven by SMARCB1 loss, but the molecular consequences of SMARCB1 loss in extra-cranial tumors have not been comprehensively described and genomic resources for analyses of extra-cranial MRT are limited. To provide such data, we used whole-genome sequencing, whole-genome bisulfite sequencing, whole transcriptome (RNA-seq) and microRNA sequencing (miRNA-seq), and histone modification profiling to characterize extra-cranial MRTs. Our analyses revealed gene expression and methylation subgroups and focused on dysregulated pathways, including those involved in neural crest development.

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.

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.

A pan-cancer analysis of alternative splicing events reveals novel tumor-associated splice variants of matriptase.

High-throughput transcriptome sequencing allows identification of cancer-related changes that occur at the stages of transcription, pre-messenger RNA (mRNA), and splicing. In the current study, we devised a pipeline to predict novel alternative splicing (AS) variants from high-throughput transcriptome sequencing data and applied it to large sets of tumor transcriptomes from The Cancer Genome Atlas (TCGA). We identified two novel tumor-associated splice variants of matriptase, a known cancer-associated gene, in the transcriptome data from epithelial-derived tumors but not normal tissue. Most notably, these variants were found in 69% of lung squamous cell carcinoma (LUSC) samples studied. We confirmed the expression of matriptase AS transcripts using quantitative reverse transcription PCR (qRT-PCR) in an orthogonal panel of tumor tissues and cell lines. Furthermore, flow cytometric analysis confirmed surface expression of matriptase splice variants in chinese hamster ovary (CHO) cells transiently transfected with cDNA encoding the novel transcripts. Our findings further implicate matriptase in contributing to oncogenic processes and suggest potential novel therapeutic uses for matriptase splice variants.

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.

Effects of conjugated linoleic acids on prostaglandin secretion by bovine endometrial epithelial cells in vitro.

OBJECTIVE: To determine the effects of 2 conjugated linoleic acid (CLA) isomers (cis-9, trans-11 and trans-10, cis-12) on synthesis of prostaglandin (PG) E(2) and F(2α) and expression of prostaglandin H synthase-2 (PGHS-2) of adult and fetal bovine endometrial epithelial cells in vitro. SAMPLE: Primary cultures of endometrial epithelial cells obtained from 4 adult cows and 4 fetal bovine carcasses. PROCEDURES: Cells were exposed to 0, 50, 100, or 200µM cis-9, trans-11 or trans-10, cis-12 CLA isomers for 24 hours. Culture media collected before and after 6 hours of stimulation of cells with phorbol 12-myristate 13-acetate were assayed to detect PGE(2) and PGF(2α) via ELISA. After stimulation, cells were collected for western blot analysis to quantify PGHS-2. RESULTS: Concentrations of PGF(2α) and PGE(2) were significantly lower in culture media of adult and fetal endometrial epithelial cells exposed to any concentration of either CLA than they were in media of cells not exposed to CLAs. The trans-10, cis-12 CLA isomer seemed to decrease PG production more markedly than did the cis-9, trans-11 CLA isomer. Most concentrations of both CLAs significantly reduced culture media PGE(2):PGF(2α) concentration ratios of cells. Exposure of cells to CLAs did not affect expression of PGHS-2 protein. CONCLUSIONS AND CLINICAL RELEVANCE: Results of this study indicated CLAs significantly decreased PGF(2α) and PGE(2) concentrations and PGE(2):PGF(2α) concentration ratios for cultures of adult and fetal endometrial epithelial cells with no apparent effect on PGHS-2 expression. Similar effects in cows could have effects on maternal recognition of pregnancy and immune function.

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.

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.