NGS y Metagenómica

The Next Generation Sequencing (NGS) allows to sequence the whole genome of an organism, compared to Maxam and Gilbert and Sanger sequencing that only allow to sequence, hardly, a single gene. Removing the separation of DNA fragments by electrophoresis, and the development of techniques that let the parallelization (analysing simultaneously several DNA fragments) have been crucial for the improvements of this process. The new companies in this ambit, Roche and Illumina, bet for different protocols to achieve these goals. Illumina bets for the sequencing by synthesis (SBS), requiring the library preparation and the use of adapters. Likewise, Illumina has replaced Roche because its lower rate of misincorporation, making it ideal for studies of genetic variability, transcriptomic, epigenomic, and metagenomic, in which this study will focus. However, it is noteworthy that the last progress in sequencing is carried out by the third generation sequencing, using nanotechnology to design small sequencers that sequence the whole genome of an organism quickly and inexpensively. Moreover, they provide more reliable data than current systems because they sequence a single molecule, solving the problem of synchronisation. In this way, PacBio and Nanopore allow a great progress in diagnostic and personalized medicine. Metagenomics provide to make a qualitative and quantitative analysis of the various species present in a sample. The main advantage of this technique is the no necessary isolation and growth of the species, allowing the analysis of nonculturable species. The Illumina protocol studies the variable regions of the 16S rRNA gene, which contains variable and not variables regions providing a phylogenetic classification. Therefore, metagenomics is a topic of interest to know the biodiversity of complex ecosystems and to study the microbiome of patients given the high involvement with certain microbial profiles on the condition of certain metabolic diseases.

Degradación bacteriana de cianuro y compuestos nitrogenados tóxicos

La ruta de asimilación de cianuro en P. pseudoalcaligenes CECT5344 transcurre a través de un nitrilo formado por la reacción química del cianuro con el oxalacetato, siendo este último acumulado como consecuencia de la acción conjunta de una malato:quinona oxidoreductasa (MQO) y la oxidasa terminal resistente a cianuro (CioAB) (Luque-Almagro et al., 2011b). Los nitrilos pueden ser convertidos en amonio por la acción de una nitrilasa o un sistema nitrilo hidratasa/amidasa. Con el objetivo de elucidar la ruta de asimilación de cianuro en P. pseudoalcalígenes CECT5344, se ha analizado el proteoma de este microorganismo en condiciones cianotróficas frente a nitrato como fuente de nitrógeno como control. En este estudio se identificaron proteínas relacionadas con la ruta de asimilación de cianuro en la estirpe CECT5344, que aparecían inducidas por cianuro, como NitB y NitG, cuyos genes se encuentran localizados en la agrupación génica nit1C. Además de NitB y NitG, de función desconocida, la agrupación génica nit1C codifica un regulador transcripcional del tipo Fis dependiente de σ54 (NitA), una nitrilasa (NitC), una proteína que pertenece a la superfamilia S-adenosilmetionina (NitD), un miembro de la superfamilia N-aciltransferasa (NitE), un polipéptido de la familia AIRS/GARS (NitF) y una oxidorreductasa dependiente de NADH (NitH). Un análisis transcripcional mediante RT-PCR determinó que los genes nitBCDEFGH se cotranscriben, mientras que el gen regulador nitA se transcribe de forma divergente. Además, resultados obtenidos por RT-PCR confirman que la expresión de los genes nitBCDEFGH está inducida por cianuro y reprimida por amonio. La relación entre el cianuro y el grupo de genes nit1C queda patente por el fenotipo de los mutantes deficientes nitA, nitB y nitC, incapaces de usar complejos cianuro-metálicos o 2-hidroxinitrilos como única fuente de nitrógeno. Todos estos datos indican que la nitrilasa NitC, junto con la proteína NitB, utilizan de forma específica determinados nitrilos alifáticos como sustrato, entre los que se encuentran el formado durante la asimilación de cianuro (Estepa et al., 2012). Además, entre las proteínas inducidas por cianuro se identificaron una dihidropicolinato sintasa (DapA), una fosfoserina transaminasa (SerC) y una proteína de función desconocida (Orf1), las tres codificadas por genes del operón cio, una cianasa (CynS), la proteína S6 de la subunidad ribosomal 30S (RpsF), una superóxido dismutasa (SodB), la ferritina (Dps), una oxidorreductasa (Fpr) y un factor de elongación P (EF-P). Una vez identificadas, estas proteínas se han analizado funcionalmente y se han localizado en el genoma de P. pseudoalcaligenes CECT5344 los genes correspondientes, así como los genes adyacentes. La inducción de estas proteínas en condiciones cianotróficas sugiere que el metabolismo del cianuro incluye, además de la resistencia y asimilación de este tóxico, otros procesos biológicos relacionados con el metabolismo del cianato y de algunos aminoácidos, el estrés oxidativo y la homeostasis de hierro, entre otros. Por otra parte, el conocimiento en profundidad y la interpretación de la secuencia génica de P. pseudoalcaligenes CECT5344, así como el análisis comparativo frente a organismos no cianotrofos ha permitido entender algunos de los mecanismos implicados en la resistencia y asimilación de cianuro, lo que permitiría conducir a la posterior mejora del proceso de biodegradación de cianuro. Además, el estudio del genoma de la estirpe CECT5344 permitirá explorar la capacidad de este organismo para ser utilizado en procesos de biorremediación de residuos cianurados en los que se encuentran metales y otros tóxicos (Luque-Almagro et al., 2013; Wibberg et al., 2014). En este trabajo se muestran y discuten los resultados de la secuenciación del genoma de P. pseudoalcaligenes, así como el estudio del análisis filogenético y evolutivo de la cepa, estableciéndose de esta manera relaciones con otras especies en base a los genomas secuenciados de las mismas, entre las que destaca P. mendocina ymp relacionada con P. pseudoalcaligenes CECT5344. El estudio de las características del genoma de P. pseudoalcaligenes CECT5344 ha sido completado con un análisis comparativo frente a los genomas de otras especies de Pseudomonas, encontrándose así semejanzas y diferencias en cuanto a la distribución génica funcional. Por último, se muestra un análisis del genoma de P. pseudoalcaligenes CECT5344 en relación con los genes implicados probablemente en los procesos de asimilación de cianuro y residuos cianurados, tales como los codificantes de nitrilasas y aquellos implicados en la resistencia a cianuro como los constituyentes del operón cio que codifican la oxidasa terminal insensible a cianuro. Finalmente, se discute la presencia de genes implicados posiblemente en otros procesos con una alto potencial biotecnológico, tales como la producción de bioplásticos y la biodegradación de diversos contaminantes.

Insights into Neural Crest Evolution

Neural crest cells are unique to vertebrates and essential to the development and evolution of the craniofacial skeleton. Using a combination of DiI cell lineage tracing, transcriptomics, and analysis of key transcription factors of the Sox Family, I examined neural crest development in the sea lamprey, Petromyzon marinus, as the most basal extant vertebrate from which it is possible to get embryos. The results have uncovered distinct cranial and trunk neural crest subpopulations along the anterior-posterior axis of the lamprey embryo, with a clear separation between the two. However, no evidence of the presence of an intermediate vagal neural crest population was uncovered. Comparing cranial neural crest genes between lamprey and chick, either by examining individual candidate genes or whole genome transcriptome analysis, reveals significant changes in the cranial neural crest gene regulatory network of lamprey compared with chick. In particular, the lamprey cranial neural crest is "missing" several gnathostome cranial crest genes. We speculate that these may underlie the evolutionary divergence of craniofacial development between jawed and jawless vertebrates. Despite the absence of vagal neural crest, DiI-labeling shows that trunk neural crest-derived cells, likely homologous to mammalian Schwann cell precursors, contribute to the lamprey enteric nervous system, potentially representing the most primitive form of neural crest cells contribution to the ENS. Finally, I characterized key members of the Sox Family (Sox B-F) due to their importance in neural crest specification in other species. In comparative studies of the SoxC genes (Sox4, Sox11, and Sox12) in both lamprey and Xenopus, I found similar expression patterns and a novel key role in early neural crest specification, suggesting a conserved role of the SoxC genes amongst vertebrates. Taken together, this work represents important progress in characterizing the early evolution of the neural crest in vertebrates and its role in the transition from jawless to jawed vertebrates.

Hormone-Mediated Gene-Specific Translation Regulation

The central role of translation regulation in the control of critical cellular processes has long been recognized. Yet the systematic exploration of quantitative changes in translation at a genome-wide scale in response to specific stimuli has only recently become technically feasible. Using a genetic approach, we have identified new Arabidopsis weak-ethylene insensitive mutants that also display defects in translation, which suggested the existence of a previously unknown molecular module involved in ethylene-mediated translation regulation of components of this signaling pathway. To explore this link in detail, we implemented for Arabidopsis the ribosome-footprinting technology, which enables the study of translation at a whole-genome level at single codon resolution[1]. Using ribosome-footprinting we examined the effects of short exposure to ethylene on the Arabidopsis translatome looking for ethylene-triggered changes in translation rates that could not be explained by changes in transcript levels. The results of this research, in combination with the characterization of a subset of the aforementioned weak-ethylene insensitive mutants that are defective in the UPF genes (core-components of the nonsense-mediated mRNA decay machinery), uncovered a translation-based branch of the ethylene signaling pathway[2]. In the presence of ethylene, translation of a negative regulator of ethylene signaling EBF2 is repressed, despite induced transcription of this gene. These translational effects of ethylene require the long 3´UTR of EBF2 (3´EBF2), which is recognized by the C-terminal end of the key ethylene-signaling protein EIN2 (EIN2C) in the cytoplasm once EIN2C is released from the ER-membrane by proteolytic cleavage. EIN2C binds the 3´EBF2, recruits the UPF proteins and moves to P-bodies, where the translation of EBF2 in inhibited despite its mRNA accumulation. Once the ethylene signal is withdrawn, the translation of the stored EBF2 mRNAs is resumed, thus rapidly dampening the ethylene response. These findings represent a mechanistic paradigm of gene-specific regulation of translation in response to a key growth regulator. Translation regulatory elements can be located in both 3′ and 5′ UTRs. We are now focusing on the ead1 and ead2 mutants, another set of ethylene-signaling mutants defective in translational regulation. Ribosome-footprinting on the ead1 mutant revealed an accumulation of translating ribosomes in the 5´UTRs of uORF-containing genes and reduction in the levels of ribosomes in the main ORF. The mutant is also impaired in the translation of GFP when this reporter is fused to WT 5´UTR of potential EAD1 targets but not when GFP is fused to the uORF-less versions of the same 5´UTRs. Our hypothesis is that EAD1/2 work as a complex that is required for the efficient translation of mRNAs that have common structural (complex 5´UTR with uORFs) and functional (regulation of key cellular processes) features. We are working towards the identification of the conditions where the EAD1 regulation of translation is required. [1] Ingolia, N. et al. (2009) Genome-Wide Analysis in Vivo of Translation with Nucleotide Resolution Using Ribosome Profiling. Science, 324; 218-222 [2] Merchante, C. et al. (2015) Gene-Specific Translation Regulation Mediated by the Hormone-Signaling Molecule EIN2. Cell, 163(3): 684-697

Coexistence of mcr-1 and blaNDM-1 in Escherichia coli from Venezuela

We studied the presence of the mobile colistin resistance gene mcr-1 in human, animal, and environmental Enterobacteriaceae samples from Cumana, Venezuela, that were collected in 2015. The mcr-1 gene was detected in 2/93 Escherichia coli isolates from swine (novel ST452) and human (ST19) samples that were resistant to colistin. Whole-genome sequencing and transformation experiments identified mcr-1 on an IncI2 plasmid. One of the isolates also bore the widely spread carbapenemase NDM-1. A One Health approach is necessary to further elucidate the flux of these high-risk genes.

Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.

Global transcriptome analysis of Lactococcus garvieae strains in response to temperature

Lactococcus garvieae is an important fish and an opportunistic human pathogen. The genomic sequences of several L. garvieae strains have been recently published, opening the possibility of global studies on the biology of this pathogen. In this study, a whole genome DNA microarray of two strains of L. garvieae was designed and validated. This DNA microarray was used to investigate the effects of growth temperature (18°C and 37°C) on the transcriptome of two clinical strains of L. garvieae that were isolated from fish (Lg8831) and from a human case of septicemia (Lg21881). The transcriptome profiles evidenced a strain-specific response to temperature, which was more evident at 18°C. Among the most significant findings, Lg8831 was found to up-regulate at 18°C several genes encoding different cold-shock and cold-induced proteins involved in an efficient adaptive response of this strain to low-temperature conditions. Another relevant result was the description, for the first time, of respiratory metabolism in L. garvieae, whose gene expression regulation was temperature-dependent in Lg21881. This study provides new insights about how environmental factors such as temperature can affect L. garvieae gene expression. These data could improve our understanding of the regulatory networks and adaptive biology of this important pathogen.

Staphylococcus aureus CC398: host adaptation and emergence of methicillin resistance in livestock

UNLABELLED Since its discovery in the early 2000s, methicillin-resistant Staphylococcus aureus (MRSA) clonal complex 398 (CC398) has become a rapidly emerging cause of human infections, most often associated with livestock exposure. We applied whole-genome sequence typing to characterize a diverse collection of CC398 isolates (n = 89), including MRSA and methicillin-susceptible S. aureus (MSSA) from animals and humans spanning 19 countries and four continents. We identified 4,238 single nucleotide polymorphisms (SNPs) among the 89 core genomes. Minimal homoplasy (consistency index = 0.9591) was detected among parsimony-informative SNPs, allowing for the generation of a highly accurate phylogenetic reconstruction of the CC398 clonal lineage. Phylogenetic analyses revealed that MSSA from humans formed the most ancestral clades. The most derived lineages were composed predominantly of livestock-associated MRSA possessing three different staphylococcal cassette chromosome mec element (SCCmec) types (IV, V, and VII-like) including nine subtypes. The human-associated isolates from the basal clades carried phages encoding human innate immune modulators that were largely missing among the livestock-associated isolates. Our results strongly suggest that livestock-associated MRSA CC398 originated in humans as MSSA. The lineage appears to have undergone a rapid radiation in conjunction with the jump from humans to livestock, where it subsequently acquired tetracycline and methicillin resistance. Further analyses are required to estimate the number of independent genetic events leading to the methicillin-resistant sublineages, but the diversity of SCCmec subtypes is suggestive of strong and diverse antimicrobial selection associated with food animal production. IMPORTANCE Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production.

FROM GENES TO GENOMES: LOCAL ADAPTATION AND ADAPTIVE POTENTIAL IN TWO ARCTIC SEABIRDS

By investigating the mechanisms underlying the evolution and the maintenance of local adaptations we can help predict how species will adapt to future environmental change. In this thesis I investigate local adaptation and adaptive potential in thick-billed and common murres (Uria lomvia and U. aalge), two arctic seabirds of international conservation concern. Thanks to the recent development of new genomic methods, I address three major themes that are relevant for both the development of evolutionary theory and conservation: 1) the role of gene flow in the origin and maintenance of adaptation; 2) levels and distribution of standing genetic variation, and their contribution to adaptive potential; and 3) the genomic mechanisms maintaining an adaptive dimorphism within a single interbreeding population. First, I review the literature on genomics of local adaptation with gene flow and find that adaptation can be maintained despite gene flow, that gene flow itself can promote adaptation, and that genetic architecture is important in the origin and maintenance of local adaptations. Second, I genotype genome-wide markers and toll-like receptor genes (TLRs) to investigate local adaptation and adaptive potential in thick-billed murres. Thick-billed murres do not show signatures of local adaptation to their breeding grounds, but outlier loci group birds according to their non-breeding distributions, suggesting that selection and/or demographic connectivity in the winter may explain patterns of differentiation in this species. Genetic variation at TLRs does not decrease with increasing latitude as predicted, but tests of selection and measures of genetic diversity suggest differences in local selective regimes at most genes. Thick-billed murres show high levels of standing genetic variation and their adaptive potential will mostly depend on rate and magnitude of environmental change. Finally, I improve and annotate the assembly of the highly heterozygous genome of the thick-billed murre. Using this assembly as a reference, I perform whole genome analyses to investigate the genomic basis of an adaptive dimorphism in Atlantic common murres. I show for the first time that a 60 kb complex copy number variant in a non-coding region maintains differences in plumage and cold adaptation despite high gene flow.

Genomics to elucidate the molecular basis of calcific aortic valve disease

Le rétrécissement valvulaire aortique (RVA) est causé par une calcification et une fibrose progressive de la valve aortique. Le risque de développer la maladie augmente avec l’âge. À cause de l’augmentation de l’espérance de vie, le RVA est devenu un problème de santé publique. Le RVA est fatal en absence de traitement médical. Actuellement, la chirurgie est le seul traitement pour le stade sévère de la maladie, mais près de 50% des individus avec RVA n’y sont pas éligibles, principalement due à la présence de comorbidités. Plusieurs processus biologiques ont été associés à la maladie, mais les voies moléculaires spécifiques et les gènes impliqués dans le développement et la progression du RVA ne sont pas connus. Il est donc urgent de découvrir les gènes de susceptibilité pour le RVA afin d’identifier les personnes à risque ainsi que les biomarqueurs et les cibles thérapeutiques pouvant mener au développement de médicaments pour inverser ou limiter la progression de la maladie. L’objectif de cette thèse de doctorat était d’identifier la base moléculaire du RVA. Des approches modernes en génomique, incluant l’étude de gènes candidats et le criblage génomique par association (GWAS), ont été réalisées à l’aide de collections d’ADN provenant d’un grand nombre de patients bien caractérisés pour le RVA. Des études complémentaires en transciptomique ont comparé le profil d’expression global des gènes entre des valves calcifiées et non-calcifiées à l’aide de biopuces à ADN et de séquençage de l’ARN. Une première étude a identifié des variations dans le gène NOTCH1 et suggère pour la première fois la présence d’un polymorphisme commun dans ce gène conférant une susceptibilité au RVA. La deuxième étude a combiné par méta-analyse deux GWAS de patients provenant de la ville de Québec et Paris (France) aux données transcriptomiques. Cette étude de génomique intégrative a confirmé le rôle de RUNX2 dans le RVA et a permis l’identification d’un nouveau gène de susceptibilité, CACNA1C. Les troisième et quatrième études sur l’expression des gènes ont permis de mieux comprendre les bases moléculaires de la calcification des valves aortiques bicuspides et ainsi d’identifier de nouvelles cibles thérapeutiques pour le RVA. Les données générées par ce projet sont la base de futures découvertes importantes qui permettront d’améliorer les options de traitement et la qualité de vie des patients atteints du RVA.

Evidences for an opportunistic and endophytic lifestyle of the Bursaphelenchus xylophilus -associated bacteria Serratia marcescens PWN146 isolated from wilting Pinus pinaster

Pine wilt disease (PWD) results from the interaction of three elements: the pathogenic nematode, Bursaphelenchus xylophilus; the insect-vector, Monochamus sp.; and the host tree, mostly Pinus species. Bacteria isolated from B. xylophilus may be a fourth element in this complex disease. However, the precise role of bacteria in this interaction is unclear as both plant-beneficial and as plant-pathogenic bacteria may be associated with PWD. Using whole genome sequencing and phenotypic characterization, we were able to investigate in more detail the genetic repertoire of Serratia marcescens PWN146, a bacterium associated with B. xylophilus. We show clear evidence that S. marcescens PWN146 is able to withstand and colonize the plant environment, without having any deleterious effects towards a susceptible host (Pinus thunbergii), B. xylophilus nor to the nematode model C. elegans. This bacterium is able to tolerate growth in presence of xenobiotic/organic compounds, and use phenylacetic acid as carbon source. Furthermore, we present a detailed list of S. marcescens PWN146 potentials to interfere with plant metabolism via hormonal pathways and/or nutritional acquisition, and to be competitive against other bacteria and/or fungi in terms of resource acquisition or production of antimicrobial compounds. Further investigation is required to understand the role of bacteria in PWD. We have now reinforced the theory that B. xylophilus-associated bacteria may have a plant origin.

Phenotypic and genotypic characterisation of antimicrobial resistance in Escherichia coli indicator of animal, food and human origin.

This thesis presents AMR phenotypic evaluation and whole genome sequencing analysis of 288 Escherichia coli strains isolated from different sources (livestock, companion animal, wildlife, food and human) in Italy. Our data reflects general resistance trends in Europe, reporting tetracycline, ampicillin, sulfisoxazole and aminoglycosides resistance as the most common phenotypic AMR profile among livestock, pets, wildlife and humans. Identification of human and animal (livestock and companion animal) AMR profiles in niches with a rare (fishery, mollusc) or absent (vegetable, wild animal, wild boar) direct exposure to antimicrobials, suggests widespread environmental pollution with ARGs conferring resistance to these antimicrobials. Phenotypic resistance to highest priority critically important antimicrobials was mainly observed in food-producing animals and related food such as rabbit, poultry, beef and swine. Discrepancies between AMR phenotypic pattern and genetic profile were observed. In particular, phenotypic aminoglycoside, cephalosporin, meropenem, colistin resistance and ESBL profile did not have a genetic explanation in different cases. This data could suggest the diffusion of new genetic variants of ARGs, associated to these antimicrobial classes. Generally, our collection shows a virulence profile typical of extraintestinal pathogenic Escherichia coli (ExPEC) pathotype. Different pandemic and emerging ExPEC lineages were identified, in particular in poultry meat (ST10; ST23; ST69, ST117; ST131). Rabbit was suggested as a source of ST20-ST40 potential hybrid pathogens. Wildlife carried a high average number (10) of VAGs (mostly associated to ExPEC pathotype) and different predominant ExPEC lineages (ST23, ST117, ST648), suggesting its possible involvement in maintenance and diffusion of virulence determinants. In conclusion, our study provides important knowledge related to the phenotypic/genetic AMR and virulence profiles circulating in E. coli in Italy. The role of different niches in AMR dynamics has been discussed. In particular, food-producing animals are worthy of continued investigation as a source of potential zoonotic pathogens, meanwhile wildlife might contribute to VAGs spread.

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.

Contribution of rare damaging variants in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a heterogeneous and highly heritable neurodevelopmental disorder with a complex genetic architecture, consisting of a combination of common low-risk and more penetrant rare variants. This PhD project aimed to explore the contribution of rare variants in ASD susceptibility through NGS approaches in a cohort of 106 ASD families including 125 ASD individuals. Firstly, I explored the contribution of inherited rare variants towards the ASD phenotype in a girl with a maternally inherited pathogenic NRXN1 deletion. Whole exome sequencing of the trio family identified an increased burden of deleterious variants in the proband that could modulate the CNV penetrance and determine the disease development. In the second part of the project, I investigated the role of rare variants emerging from whole genome sequencing in ASD aetiology. To properly manage and analyse sequencing data, a robust and efficient variant filtering and prioritization pipeline was developed, and by its application a stringent set of rare recessive-acting and ultra-rare variants was obtained. As a first follow-up, I performed a preliminary analysis on de novo variants, identifying the most likely deleterious variants and highlighting candidate genes for further analyses. In the third part of the project, considering the well-established involvement of calcium signalling in the molecular bases of ASD, I investigated the role of rare variants in voltage-gated calcium channels genes, that mainly regulate intracellular calcium concentration, and whose alterations have been correlated with enhanced ASD risk. Specifically, I functionally tested the effect of rare damaging variants identified in CACNA1H, showing that CACNA1H variation may be involved in ASD development by additively combining with other high risk variants. This project highlights the challenges in the analysis and interpretation of variants from NGS analysis in ASD, and underlines the importance of a comprehensive assessment of the genomic landscape of ASD individuals.

Exploring mediterranean durum wheat wild relatives and elite diversity panel for future breeding

Durum wheat (Triticum durum) is an important crop that has been used for millennia for human consumption, and modern breeding can take advantage of the wide variability useful for the adaptation to new challenges. Novel beneficial alleles can be found in wild relatives and landraces thus enhancing crop adaptation to many biotic and abiotic stresses. This dissertation considers the source of variability from both before and after wheat domestication, that caused a loss of potentially useful alleles. Chapter 1. is the thesis introduction, which outlines the importance of wheat in the world, providing an historical overview of the domestication, the evolution mechanisms that led to the current forms of durum wheat and the use of wild relatives as a source of germplasm for future breeding programs is crucial. Moreover, the emergence of Z. tritici has been considered as the main pathogen of wheat since it contains extremely high levels of genetic variability and is thus difficult to control. Chapter 2. Considers the contribution of the phenotypic diversity of 242 accessions of Aegilops tauschii from the Open Wild Wheat Consortium, involved in wheat domestication, provided with whole-genome resequencing. The accessions were phenotyped both in the field and in controlled conditions and A k-mer-based GWAS was performed to identify genomic regions involved in useful traits. Chapter 3. Describes the genetic basis of resistance to Z. tritici in a durum wheat elite diversity panel representative of the germplasm bred in Mediterranean. Quantitative trait loci (QTL) analysis results revealed several loci involved in the STB response that were found in several chromosome regions with a high infection rate. The genomic regions associated with STB resistance identified in this study could be of interest for marker assisted selection (MAS) in durum wheat breeding programs.