Microarrays as a functional approach to the transcriptome

Knowing a cell’s transcriptome is a fundamental requisite in order to analyze its response to the environment. Microarrays have supposed a revolution on this field as they are able to yield an overview of gene expression at any environmental condition on a genome-wide scale. This technique consists in the hybridisation of a nucleic acid sample, previously marked, with a probe (which might be made up of cDNA, oligonucleotides or PCR products) anchored to a solid surface (made of glass, plastic, silicon...) giving as a result a dot grid which reveals, after image analysis, which genes are being expressed. Nevertheless, this only can be achieved if information on the species genome has been generated. Different kinds of expression microarrays exist attending to the probe’s nature and the method used in its synthesis. In this poster two of these will be treated: Spotted Microarrays, for which the probe is synthesised prior to its fixation to the array and allow the analysis of two targets simultaneously. They can be easily customized, but lack high reproducibility and sensitivity. Oligonucleotide Microarrays, which are characterized by the direct printing of the probe on the array. In this case the probes consist on, invariably, oligonucleotides that are complementary to a small fraction of the gene it is representing at the microarray. Their application is somewhat restricted. This fact, however, makes them more reproducible. Currently, the approach towards the transcriptome studies from the Next Generation Sequencing technologies offers a large volume of information in a short amount of time needing less previous information on the target organism than that needed by microarrays, but their expensive price limits their use. The versatility of the latter, together with their reduced costs in comparison to other techniques, makes them an interesting resource in applications that may need less complexity.

Complete genome of a dengue virus serotype 4 strain from Amazonas, Brazil

Dengue virus (DENV) infections represent a significant concern for public health worldwide, being considered as the most prevalent arthropod-borne virus regarding the number of reported cases. In this study, we report the complete genome sequencing of a DENV serotype 4 isolate, genotype II, obtained in the city of Manaus, directly from the serum sample, applying Ion Torrent sequencing technology. The use of a massive sequencing technology allowed the detection of two variable sites, one in the coding region for the viral envelope protein and the other in the nonstructural 1 coding region within viral populations.

Differential expression in “rogue” paramutation in peas (Pisum sativum L.)

Dissertação de Mestrado, Engenharia Biológica, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2014

Studies of MHC class I antigen presentation & the origins of the immunopeptidome

La présentation d'antigène par les molécules d'histocompatibilité majeure de classe I (CMHI) permet au système immunitaire adaptatif de détecter et éliminer les agents pathogènes intracellulaires et des cellules anormales. La surveillance immunitaire est effectuée par les lymphocytes T CD8 qui interagissent avec le répertoire de peptides associés au CMHI présentés à la surface de toutes cellules nucléées. Les principaux gènes humains de CMHI, HLA-A et HLA-B, sont très polymorphes et par conséquent montrent des différences dans la présentation des antigènes. Nous avons étudié les différences qualitatives et quantitatives dans l'expression et la liaison peptidique de plusieurs allotypes HLA. Utilisant la technique de cytométrie de flux quantitative nous avons établi une hiérarchie d'expression pour les quatre HLA-A, B allotypes enquête. Nos résultats sont compatibles avec une corrélation inverse entre l'expression allotypique et la diversité des peptides bien que d'autres études soient nécessaires pour consolider cette hypothèse. Les origines mondiales du répertoire de peptides associés au CMHI restent une question centrale à la fois fondamentalement et dans la recherche de cibles immunothérapeutiques. Utilisant des techniques protéogénomiques, nous avons identifié et analysé 25,172 peptides CMHI isolées à partir des lymphocytes B de 18 personnes qui exprime collectivement 27 allotypes HLA-A,B. Alors que 58% des gènes ont été la source de 1-64 peptides CMHI par gène, 42% des gènes ne sont pas représentés dans l'immunopeptidome. Dans l'ensemble, l’immunopeptidome présenté par 27 allotypes HLA-A,B ne couvrent que 17% des séquences exomiques exprimées dans les cellules des sujets. Nous avons identifié plusieurs caractéristiques des transcrits et des protéines qui améliorent la production des peptides CMHI. Avec ces données, nous avons construit un modèle de régression logistique qui prédit avec une grande précision si un gène de notre ensemble de données ou à partir d'ensembles de données indépendants génèrerait des peptides CMHI. Nos résultats montrent la sélection préférentielle des peptides CMHI à partir d'un répertoire limité de produits de gènes avec des caractéristiques distinctes. L'idée que le système immunitaire peut surveiller des peptides CMHI couvrant seulement une fraction du génome codant des protéines a des implications profondes dans l'auto-immunité et l'immunologie du cancer.

POPLAR GENE EXPRESSION DATA ANALYSIS PIPELINES

Analyzing large-scale gene expression data is a labor-intensive and time-consuming process. To make data analysis easier, we developed a set of pipelines for rapid processing and analysis poplar gene expression data for knowledge discovery. Of all pipelines developed, differentially expressed genes (DEGs) pipeline is the one designed to identify biologically important genes that are differentially expressed in one of multiple time points for conditions. Pathway analysis pipeline was designed to identify the differentially expression metabolic pathways. Protein domain enrichment pipeline can identify the enriched protein domains present in the DEGs. Finally, Gene Ontology (GO) enrichment analysis pipeline was developed to identify the enriched GO terms in the DEGs. Our pipeline tools can analyze both microarray gene data and high-throughput gene data. These two types of data are obtained by two different technologies. A microarray technology is to measure gene expression levels via microarray chips, a collection of microscopic DNA spots attached to a solid (glass) surface, whereas high throughput sequencing, also called as the next-generation sequencing, is a new technology to measure gene expression levels by directly sequencing mRNAs, and obtaining each mRNA’s copy numbers in cells or tissues. We also developed a web portal (http://sys.bio.mtu.edu/) to make all pipelines available to public to facilitate users to analyze their gene expression data. In addition to the analyses mentioned above, it can also perform GO hierarchy analysis, i.e. construct GO trees using a list of GO terms as an input.

Studies of MHC class I antigen presentation & the origins of the immunopeptidome

La présentation d'antigène par les molécules d'histocompatibilité majeure de classe I (CMHI) permet au système immunitaire adaptatif de détecter et éliminer les agents pathogènes intracellulaires et des cellules anormales. La surveillance immunitaire est effectuée par les lymphocytes T CD8 qui interagissent avec le répertoire de peptides associés au CMHI présentés à la surface de toutes cellules nucléées. Les principaux gènes humains de CMHI, HLA-A et HLA-B, sont très polymorphes et par conséquent montrent des différences dans la présentation des antigènes. Nous avons étudié les différences qualitatives et quantitatives dans l'expression et la liaison peptidique de plusieurs allotypes HLA. Utilisant la technique de cytométrie de flux quantitative nous avons établi une hiérarchie d'expression pour les quatre HLA-A, B allotypes enquête. Nos résultats sont compatibles avec une corrélation inverse entre l'expression allotypique et la diversité des peptides bien que d'autres études soient nécessaires pour consolider cette hypothèse. Les origines mondiales du répertoire de peptides associés au CMHI restent une question centrale à la fois fondamentalement et dans la recherche de cibles immunothérapeutiques. Utilisant des techniques protéogénomiques, nous avons identifié et analysé 25,172 peptides CMHI isolées à partir des lymphocytes B de 18 personnes qui exprime collectivement 27 allotypes HLA-A,B. Alors que 58% des gènes ont été la source de 1-64 peptides CMHI par gène, 42% des gènes ne sont pas représentés dans l'immunopeptidome. Dans l'ensemble, l’immunopeptidome présenté par 27 allotypes HLA-A,B ne couvrent que 17% des séquences exomiques exprimées dans les cellules des sujets. Nous avons identifié plusieurs caractéristiques des transcrits et des protéines qui améliorent la production des peptides CMHI. Avec ces données, nous avons construit un modèle de régression logistique qui prédit avec une grande précision si un gène de notre ensemble de données ou à partir d'ensembles de données indépendants génèrerait des peptides CMHI. Nos résultats montrent la sélection préférentielle des peptides CMHI à partir d'un répertoire limité de produits de gènes avec des caractéristiques distinctes. L'idée que le système immunitaire peut surveiller des peptides CMHI couvrant seulement une fraction du génome codant des protéines a des implications profondes dans l'auto-immunité et l'immunologie du cancer.

Cartographie des cassures bicaténaires du remodelage chromatinien du spermatide et développement des outils techniques associés.

Résumé : La phase haploïde de la spermatogenèse (spermiogenèse) est caractérisée par une modification importante de la structure de la chromatine et un changement de la topologie de l’ADN du spermatide. Les mécanismes par lesquels ce changement se produit ainsi que les protéines impliquées ne sont pas encore complètement élucidés. Mes travaux ont permis d’établir la présence de cassures bicaténaires transitoires pendant ce remodelage par l’essai des comètes et l’électrophorèse en champ pulsé. En procédant à des immunofluorescences sur coupes de tissus et en utilisant un extrait nucléaire hautement actif, la présence de topoisomérases ainsi que de marqueurs de systèmes de réparation a été confirmée. Les protéines de réparation identifiées font partie de systèmes sujets à l’erreur, donc cette refonte structurale de la chromatine pourrait être génétiquement instable et expliquer le biais paternel observé pour les mutations de novo dans de récentes études impliquant des criblages à haut débit. Une technique permettant l’immunocapture spécifique des cassures bicaténaires a été développée et appliquée sur des spermatides murins représentant différentes étapes de différenciation. Les résultats de séquençage à haut débit ont montré que les cassures bicaténaires (hotspots) de la spermiogenèse se produisent en majorité dans l’ADN intergénique, notamment dans les séquences LINE1, l’ADN satellite et les répétions simples. Les hotspots contiennent aussi des motifs de liaisons des protéines des familles FOX et PRDM, dont les fonctions sont entre autres de lier et remodeler localement la chromatine condensée. Aussi, le motif de liaison de la protéine BRCA1 se trouve enrichi dans les hotspots de cassures bicaténaires. Celle-ci agit entre autres dans la réparation de l’ADN par jonction terminale non-homologue (NHEJ) et dans la réparation des adduits ADN-topoisomérase. De façon remarquable, le motif de reconnaissance de la protéine SPO11, impliquée dans la formation des cassures méiotiques, a été enrichi dans les hotspots, ce qui suggère que la machinerie méiotique serait aussi utilisée pendant la spermiogenèse pour la formation des cassures. Enfin, bien que les hotspots se localisent plutôt dans les séquences intergéniques, les gènes ciblés sont impliqués dans le développement du cerveau et des neurones. Ces résultats sont en accord avec l’origine majoritairement paternelle observée des mutations de novo associées aux troubles du spectre de l’autisme et de la schizophrénie et leur augmentation avec l’âge du père. Puisque les processus du remodelage de la chromatine des spermatides sont conservés dans l’évolution, ces résultats suggèrent que le remodelage de la chromatine de la spermiogenèse représente un mécanisme additionnel contribuant à la formation de mutations de novo, expliquant le biais paternel observé pour certains types de mutations.

The effect of the HLA B27 allele on the immune response to acute HCV in HIV infected patients

In mono-infected individuals, the HLA-B27 allele is strongly associated with spontaneous clearance of HCV in association with a strong CD8+ response targeted against a single epitope within the HCV RNA-dependent RNA polymerase (NS5B). We studied variation across the whole HCV genome and T cell responses over time in a rare cohort of HLA-B27+ patients with acute HCV and HIV co-infection, the majority of whom progressed to chronicity. We used next generation sequencing to detect changes within and outwith the immuno-dominant HLA-B27 restricted HCV-specific CD8+ T cell epitope NS5B2841-2849 (ARMILMTHF) during evolving progression of early HCV infection. Within the Acute HCV UK cohort, 10 patients carried the HLA B27 allele. Of these, 3/8 patients (37.5%) with HIV infection and 2/2 (100%) without HIV spontaneously cleared HCV (p=0.44). Sequential samples from nine HLA-B27+ patients (2 with monoinfection and 7 with HIV co-infection) were available for analysis (four spontaneous clearers and five evolving progressors). Mutations identified using NGS were assessed using a replicon genotype 1a system to evaluate viral fitness. Multiple mutations within the HLA-B27 restricted NS5B2841-2849 epitope were associated with progression to chroncity whereas patients who cleared the HCV infection spontaneously had no or only one mutation at this site (p=0.03). A triple NS5B2841-2849 mutant observed during progression to chronicity was associated with restored replication when compared to wild-type virus while single or double mutants were significantly associated with impaired replication (p=0.0495). T cell responses measured in these patients using ELISpot and flow cytometry. HLA-B27+ patients had significantly higher IFN-γ responses than patients who were HLA-B27- (p=0.0014). Those who progressed to chronicity had lower IFN-γ responses than those who cleared HCV (p=0.0011). Mono-infected patients had higher IFN-γ responses compared to co-infected patients (p=0.0015). HIV co-infection is associated with a lower likelihood of spontaneous clearance of HCV in HLA B27+ patients and this is associated with impaired T cell function in this group.

The impact of denture related disease on the oral microbiome of denture wearers

Advances in healthcare over the last 100 years has resulted in an ever increasing elderly population. This presents greater challenges for adequate systemic and oral healthcare delivery. With increasing age there is a natural decline in oral health, leading to the loss of teeth and ultimately for some having to wear denture prosthesis. It is currently estimated that approximately one fifth of the UK and US populations have some form of removable prosthesis. The microbiology of denture induced mucosal inflammation is a pivotal factor to consider in denture care management, similar to many other oral diseases of microbial influence, such as caries, gingivitis and periodontitis. Dentures support the growth of microbial biofilms, structures commonly known as denture plaque. Microbiologically, denture stomatitis (DS) is a disease primarily considered to be of yeast aetiology, with the literature disproportionately focussed on Candida spp. However, the denture surface is capable of carrying up to 1011 microbes per milligram, the majority of which are bacteria. Thus it is apparent that denture plaque is more diverse than we assume. There is a fundamental gap in our understanding of the bacterial composition of denture plaque and the role that they may play in denture related disease such as DS. This is categorised as inflammation of the oral mucosa, a disease affecting around half of all denture wearers. It has been proposed that bacteria and fungi interact on the denture surface and that these polymicrobial interactions lead to synergism and increased DS pathogenesis. Therefore, understanding the denture microbiome composition is the key step to beginning to understand disease pathogenesis, and ultimately help improve treatments and identify novel targets for therapeutic and preventative strategies. A group of 131 patients were included within this study in which they provided samples from their dentures, palatal mucosa, saliva and dental plaque. Microbes residing on the denture surface were quantified using standard Miles and Misra culture technique which investigated the presence of Candida, aerobes and anaerobes. These clinical samples also underwent next generation sequencing using the Miseq Illumina platform to give a more global representation of the microbes present at each of these sites in the oral cavity of these denture wearers. This data was then used to compare the composition and diversity of denture, mucosal and dental plaque between one another, as well as between healthy and diseased individuals. Additional comparisons included denture type and the presence or absence of natural teeth. Furthermore, microbiome data was used to assess differences between patients with varying levels of oral hygiene. The host response to the denture microbiome was investigated by screening the patients saliva for the presence and quantification of a range of antimicrobial peptides that are associated with the oral cavity. Based on the microbiome data an in vitro biofilm model was developed that reflected the composition of denture plaque. These biofilms were then used to assess quantitative and compositional changes over time and in response to denture cleansing treatments. Finally, the systemic implications of denture plaque were assessed by screening denture plaque samples for the presence of nine well known respiratory pathogens using quantitative PCR. The results from this study have shown that the bacterial microbiome composition of denture wearers is not consistent throughout the mouth and varies depending on sample site. Moreover, the presence of natural dentition has a significant impact on the microbiome composition. As for healthy and diseased patients the data suggests that compositional changes responsible for disease progression are occurring at the mucosa, and that dentures may in fact be a reservoir for these microbes. In terms of denture hygiene practices, sleeping with a denture in situ was found to be a common occurrence. Furthermore, significant shifts in denture microbiome composition were found in these individuals when compared to the denture microbiome of those that removed their denture at night. As for the host response, some antimicrobial peptides were found to be significantly reduced in the absence of natural dentition, indicating that the oral immune response is gradually impaired with the loss of teeth. This study also identified potentially serious systemic implications in terms of respiratory infection, as 64.6% of patients carried respiratory pathogens on their denture. In conclusion, this is the first study to provide a detailed understanding of the oral microbiome of denture wearers, and has provided evidence that DS development is more complex than simply a candidal infection. Both fungal and bacterial kingdoms clearly play a role in defining the progression of DS. The biofilm model created in this study demonstrated its potential as a platform to test novel actives. Future use of this model will aid in greater understanding of host: biofilm interactions. Such findings are applicable to oral health and beyond, and may help to identify novel therapeutic targets for the treatment of DS and other biofilm associated diseases.

Transcriptome analysis of wheat under drought.

2013

Detection and molecular characterization of Grapevine yellow speckle viroid 1 isolates infecting grapevines in Brazil.

Abstract Presently, Hop stunt viroid(HSVd) and Citrus exocortis viroid (CEVd) are the only viroids reported to infect grapevines (Vitis spp.) in Brazil, among the seven viroid species already reported infecting this host in other countries. All grapevine viroid diseases are graft-transmissible and can induce losses especiallywhenassociatedwithviruses.Theaimofthisworkwas to confirm infection by Grapevine yellow speckle viroid 1(GYSVd-1) in grapevine samples exhibiting yellow speckle symptoms in the leaves and in asymptomatic samples sequenced by next generation sequencing (NGS). The occurrence of this viroid in Brazil was further investigated in a second study. Total RNAs and dsRNAs were extracted from five symptomatic plants and 16 asymptomatic samples, respectively. Specific primers were used for RT-PCR and amplified DNA fragments were cloned and sequenced by the Sanger method. Eleven complete nucleotide sequences of GYSVd-1 isolates (366 ?367 nt) were obtained from NGS and from RT-PCR amplicons. Comparisons showed high identities (95.9 ?100 %) among ten isolates and an identity of 87.2 ?90.4 % with a divergent isolate (RM-BR). Phylogenetic analyses placed GYSVd-1 isolates in four clusters (types 1, 2, 3 and 4). All GYSVd-1 infections were confirmed by conventional RT-PCR and RT-qPCR using specific oligonucleo-tides and a labeled probe. This is the first report and molecular characterization of GYSVd-1 infecting grapevines in Brazil, and our survey indicates that this viroid could be widespread in the major grape producing regions of Brazil. Keywords GYSVd-1 . Incidence . Next generation sequencing. Secondary structure. Vine.

Addressing tumor heterogeneity in esophageal adenocarcinoma through different molecular approaches

Several studies have shown epidemiologic, clinical, immune-histochemical and molecular differences among esophageal adenocarcinomas (EAC). Since pathogenesis and biology of this tumor are far to be well defined, our study aimed to examine intra- and inter-tumor heterogeneity and to solve crucial controversies through different molecular approaches. Target sequencing was performed for sorted cancer subpopulations from formalin embedded material obtained from 38 EACs, not treated with neoadjuvant therapy. 35 out 38 cases carried at least one somatic mutation, not present in the corresponding sorted stromal cells. 73.7% of cases carried mutations in TP53 and 10.5% in CDKN2A. Mutations in other genes occurred at lower frequency, including HNF1A, not previously associated with EAC. Sorting allowed us to isolate clones with different mutational loads and/or additional copy number amplifications, confirming the high intra-tumor heterogeneity of these cancers. In our cohort TP53 gene abnormalities correlated with a better survival (P = 0.028); conversely, loss of SMAD4 protein expression was associated with a higher recurrence rate (P = 0.015). Shifting the focus on the epigenetic characterization of EAC, miR-221 and miR-483-3p resulted upregulated from the MicroRNA Array card analysis and confirmed with further testing. The up-regulation of both miRNAs correlated with clinical outcomes, in particular with a reduced cancer-specific survival (miR483-3p P=0.0293; miR221 P=0.0059). In vitro analyses demonstrated an increase for miR-483-3p (fold-change=2.7) that appear to be inversely correlated with SMAD4 expression in FLO-1 cell-line. In conclusion, selective sorting allowed to define the real mutation status and to isolate different cancer subclones. MiRNA expression analysis revealed a significant up-regulation of miR-221 and miR-483-3p, which correlated with worst prognosis, implying that they can be considered oncogenic factors in EAC. Therefore, cell sorting technologies, coupled with next generation sequencing, and the analysis of microRNA profiles seem to be promising strategies to guide treatment and help classify cancer prognosis.

Integrated molecular analysis of endometrial carcinoma: translation of biomarker profiles into the clinical practice

The Cancer Genome Atlas (TCGA) collaborative project identified four distinct prognostic groups of endometrial carcinoma (EC) based on molecular alterations: (i) the ultramutated subtype that encompassed POLE mutated (POLE) cases; (ii) the hypermutated subtype, characterized by MisMatch Repair deficiency (MMRd); (iii) the copy-number high subtype, with p53 abnormal/mutated features (p53abn); (iv) the copy-number low subtype, known as No Specific Molecular Profile (NSMP). Although the prognostic value of TCGA molecular classification, NSMP tumors present a wide variability in molecular alterations and biological aggressiveness. This study aims to investigate the impact of ARID1A and CTNNB1/β-catenin alterations by targeted Next-generation sequencing (NGS) and immunohistochemistry (IHC) in a consecutive series of 125 molecularly classified ECs. NGS and IHC were used to assign surrogate TCGA groups and to identify molecular alterations of multiple target genes including POLE, PTEN, ARID1A, CTNNB1, TP53. Associations with clinicopathologic parameters, molecular subtypes, and outcomes identified NSMP category as the most heterogeneous group in terms of clinicopathologic features and outcome. Integration of surrogate TCGA molecular classification with ARID1A and β-catenin analysis showed NSMP cases with ARID1A mutation characterized by the worst outcome with early recurrence, while NSMP tumors with ARID1A wild-type and β-catenin alteration had indolent clinicopathologic features and no recurrence. This study indicates how the identification of ARID1A and β-catenin alterations in EC represents a simple and effective way to characterize NSMP tumor aggressiveness and metastatic potential.

Computational resources for structural and metagenomic characterization of functions in bacteria and bacterial communities. Antimicrobial resistance, pathways for xenobiotic degradation and relationship between gut microbiome and ageing.

Prokaryotic organisms are one of the most successful forms of life, they are present in all known ecosystems. The deluge diversity of bacteria reflects their ability to colonise every environment. Also, human beings host trillions of microorganisms in their body districts, including skin, mucosae, and gut. This symbiosis is active for all other terrestrial and marine animals, as well as plants. With the term holobiont we refer, with a single word, to the systems including both the host and its symbiotic microbial species. The coevolution of bacteria within their ecological niches reflects the adaptation of both host and guest species, and it is shaped by complex interactions that are pivotal for determining the host state. Nowadays, thanks to the current sequencing technologies, Next Generation Sequencing, we have unprecedented tools for investigating the bacterial life by studying the prokaryotic genome sequences. NGS revolution has been sustained by the advancements in computational performance, in terms of speed, storage capacity, algorithm development and hardware costs decreasing following the Moore’s Law. Bioinformaticians and computational biologists design and implement ad hoc tools able to analyse high-throughput data and extract valuable biological information. Metagenomics requires the integration of life and computational sciences and it is uncovering the deluge diversity of the bacterial world. The present thesis work focuses mainly on the analysis of prokaryotic genomes under different aspects. Being supervised by two groups at the University of Bologna, the Biocomputing group and the group of Microbial Ecology of Health, I investigated three different topics: i) antimicrobial resistance, particularly with respect to missense point mutations involved in the resistant phenotype, ii) bacterial mechanisms involved in xenobiotic degradation via the computational analysis of metagenomic samples, and iii) the variation of the human gut microbiota through ageing, in elderly and longevous individuals.

Tumor heterogeneity in Esophageal Adenocarcinoma: a genomic approach

INTRODUCTION: Esophageal adenocarcinoma (EAC) is a severe malignancy in terms of prognosis and mortality rate. Because its great genetic heterogeneity, disputes regarding classification, prevention and treatments are still unsolved. AIM: We investigated intra- and inter-EAC heterogeneity by defining EAC’s somatic mutational profile and the role of candidate microRNAs, to correlate the molecular profile of tumors to clinical outcomes and to identify biomarkers for classification. METHODS: 38 EAC cases were analyzed via high-throughput cell sorting technology combined with targeted sequencing and whole genome low-pass sequencing. Targeted sequencing of further 169 cases was performed to widen the study. miR221 and miR483-3p expression was profiled via qPCR in 112 EACs and correlation with clinical outcomes was investigated. RESULTS: 35/38 EACs carried at least one somatic mutation absent in stromal cells. TP53 was found mutated in 73.7% of cases. Selective sorting revealed tumor subclones with different mutational loads and copy number alterations, confirming the high intra-tumor heterogeneity of EAC. Mutations were in most cases at homozygous state, and we identified alterations that were missed with the whole-tumor analysis. Mutations in HNF1A gene, not previously associated with EAC, were identified in both cohorts. Higher expression of miR483-3p and miR221 was associated with poorer cancer specific survival (P=0.0293 and P=0.0059), and recurrence in the Lauren intestinal subtype (P=0.0459 and P=0.0002). Median expression levels of miRNAs were higher in patients with advanced tumor stages. The loss of SMAD4 immunoreactivity was significantly associated with poorer cancer specific survival and recurrence (P=0.0452; P=0.022 respectively). CONCLUSION: Combining selective sorting technology and next generation sequencing allowed to better define EAC inter- and intra-tumor heterogeneity. We identified HNF1A as a new mutated gene associated to EAC that could be involved in tumor progression and promising biomarkers such as SMAD4, miR221 and miR483-3p to identify patients at higher risk for more aggressive tumors.