Abstract: Recent advances in Next-Generation Sequencing protocols have opened a variety of ways to generate data. However, each newly developed methodology is most suited to represent a certain phenomenon or molecule. The object of this analysis is to identify the most appropriate way to generate and process data to study the snoRNAs, or small nucleolar RNA. Recently, snoRNAs have been revealed as taking part in a variety of unexpected alternative functions such as splicing, resistance to oxidative shock and chromatin unwinding. Finding a method to generate and treat a large quantity of data containing snoRNAs and their potential interactors could highlight some of their unexplored roles within the cell. To tackle the problem, a new protocol was put forward. This new pipeline relies on a reverse transcriptase isolated from a bacterial group II intron which boasts a better representation of structured small RNAs such as tRNAs and snoRNAs. Indeed, when compared to data created by using the standard small RNA preparation protocol, the sequencing data generated through the group II intron retrotranscriptase gives a much fairer representation. These improvements are also present in the bioinformatics pipeline. The workflow was changed to facilitate the detection of ncRNAs. These modifications rescue millions of reads, further increasing the power of the analysis. Ultimately, such corrections increase the predictive power of sequencing data. / Des avancées récentes dans le domaine du séquençage de prochaine génération ont ouvert une panoplie de façons de générer des données. Toutefois, chaque nouvelle méthode dévelopée est souvent appropriée à la caractérisation d’un seul type de phénomène ou de molécules. L’objectif de cette analyse est d’identifier la manière la plus appropriée de générer et traiter les données pour étudier les petits ARNs nucléolaires, snoRNAs. Récemment, ceux-ci ont été révélés comme des acteurs dans une variété de fonctions alternatives comme l’épissage alternatif, la résistance au choc oxidatif et l’état de la chromatine. Il est donc impératif de trouver une méthode qui puisse traiter une large quantité de données contenant les snoRNAs et leurs intéracteurs pour découvrir les rôles encore inexplorés des snoRNAs. Dans cette optique, un nouveau protocole a été élaboré. Cette nouvelle suite d’analyses s’appuie sur une reverse transcriptase isolée d’un intron de groupe II bactérien qui affiche une meilleure représentation des petits ARNs structurés comme les tRNAs et les snoRNAs. En effet, quand les données générées à travers la méthode de préparation des libraries pour petits ARNs standard est comparée à celle basée sur la reverse transcriptase bactérienne, cette dernière donne une meilleure représentation du compte des espèces. Ces avancées sont aussi présentes dans la méthode d’analyse informatique. La suite d’outils a été modifiée afin de permettre une meilleure détection des petits ARN non-codants. Ces modifications permettent de récupérer des millions de lectures par ensemble de données ce qui augmente le pouvoir prédictif de l’analyse.

Cancer related mortality is almost always due to metastatic dissemination of the primary disease. While research into the biological mechanisms that drive the metastatic cascade continues to unravel its molecular underpinnings, progress in our understanding of biological phenomena such as tumour heterogeneity and its relevance to the origins of distant recurrence or the emergence of resistance to therapy has been limited.In parallel to major breakthroughs in the development of high throughput molecular techniques, researchers have begun to utilise next generation sequencing to explore the relationship between primary and matched metastatic tumours in diverse types of neoplasia. Despite small cohort sizes and often, a limited number of matched metastases for each patient, pioneering studies have uncovered hitherto unknown biological processes such as the occurrence of organ specific metastatic lineages, polyclonal seeding and homing of metastatic cells to the primary tumour bed. While yet other studies continue to highlight the potential of genomic analyses, at the time this thesis was started, an in-depth knowledge of disease progression and metastatic dissemination was currently lacking in breast cancers.Herein, we employed phylogenetic inference methods to investigate intra-tumour heterogeneity and evolution of breast cancers. A combination of whole exome sequencing, custom ultra-deep resequencing and copy number profiling were applied to primary tumours and their associated metastases from ten autopsied breast cancer patients. Two modes of metastatic progression were observed. In the majority of cases, all distant metastases clustered on a branch separate from their primary lesion. Clonal frequency analysis of somatic mutations showed that the metastases had a monoclonal origin and descended from a common metastatic precursor. Alternatively, the primary tumour was clustered alongside metastases with early branches leading to distant organs. This dichotomy coincided with the clinical history of the patients whereby multiple seeding events from the primary tumour alongside cascading metastasis-to-metastasis disseminations occurred in treatment naïve de novo metastatic patients, whereas descent from a common metastatic precursor was observed in patients who underwent primary surgery followed by systemic treatment. The data also showed that a distant metastasis can be horizontally cross-seeded and finally revealed a correlation between the extent of somatic point mutations private to the distant lesions and patient overall survival. In an unrelated dataset of relapsed breast cancer patients with matched primary and distant lesions profiled using whole genome sequencing, the landscape of somatic alterations confirmed the time dependency of copy number aberrations implying that cancer phylogenies can be dated using a molecular clock.The work presented here harnesses the strength of high throughput genomic techniques and state of the art phylogenetic tools to tell the evolutionary history of breast cancers. Our results show that the linear and parallel models of metastatic dissemination which have been held as near doctrines for many years are overstated point of views of cancer progression. Beyond the biological insights, these results suggest that surgical excision of the primary tumour in de novo metastatic breast cancers might reduce dissemination in selected cases hence providing a potential biological rationale for this practice. Similarly, there is no strong evidence of benefit in overall survival from surgical resection of oligo-metastases in breast cancer. From our analyses, metastatic lesions constitute an additional source of seeding and heterogeneity in advanced breast cancer. The data presented here is too small to derive practice-changing evidence, but supports the concept that resecting isolated metastases may be of clinical benefit in oligo-metastatic breast cancer patients. In both cases, results from larger prospective studies are warranted. / Doctorat en Sciences biomédicales et pharmaceutiques (Médecine) / info:eu-repo/semantics/nonPublished

Ce travail de thèse illustre l’intérêt d’utiliser le séquençage de nouvelle génération (Next- Generation Sequencing ou NGS) pour améliorer le diagnostic moléculaire des formes de prédispositions aux cancers. Nous avons appliqué deux stratégies différentes, l’une basée sur l’analyse simultanée d’un panel de gènes impliqués dans une pathologie, l’autre sur l’analyse de toutes les parties codante du génome ou exome.Nous avons montré que le développement d’un panel de 10 gènes impliqués dans la carcinogenèse des cancers coliques a permis d’améliorer le diagnostic moléculaire de ces formes de cancers. Cette technologie a diminué le délai de rendu des résultats et est plus sensible que le séquençage Sanger et la QMPSF (Quantitative Multiplex PCR of Short Fragment).L’application d’une stratégie d’analyse d’exomes de trio par NGS, à l’analyse de deux formes de cancers sporadiques de phénotype extrême, cancers du sein et de l’ovaire, a permis la détection de mutations de novo portant sur des gènes impliqués dans l’oncogenèse. Pour mettre en cause ces mutations dans le phénotype, nous avons réalisé des tests fonctionnels et des études de récurrence. / The working thesis Illustre interest of the USE next generation sequencing (NGS Next- Generation Sequencing e) pay Improve Molecular diagnosis of predisposition to cancer forms. Have we applied two different strategies you in June based on the simultaneous analysis of genes not involved in panel A pathology, the Other on the analysis of all the parties of the coding genome e exome.Nous Have Shown que le development of a panel of 10 genes involved in carcinogenesis of colon cancer a Permit to improve the molecular diagnosis of forms of CES dE cancers. This technology has reduced the Delai rendering results and is more reasonable That Sanger sequencing and QMPSF (Quantitative Multiplex PCR of Short Fragment) .The application OF exomes an analysis by NGS trio Strategy at analysis of two forms of sporadic cancer phenotype extreme, breast cancer and ovarian cancer, permit the detection of de novo mutations on genes involved in oncogenesis. For Questioning mutations IN THESE phenotype, We Have made Functional testing and induction studies.

Observations of floral development indicate that floral organ initiation in pentapetalous flowers more commonly results in a medially positioned abaxial petal (MAB) than in a medially positioned adaxial petal (MAD), where the medial plane is defined by the stem and the bract during early floral development. It was proposed that the dominant MAB petal initiation might impose a developmental constraint that leads to the evolution of limited patterns of floral zygomorphy in Asteridae, a family in which the floral zygomorphy develops along the medial plane and results in a central ventral (CV) petal in mature flowers. Here, I investigate whether the pattern of floral organ initiation may limit patterns of floral zygomorphy to evolve in pentapetalous angiosperms. I analyzed floral diagrams representing 405 species in 330 genera of pentapetalous angiosperms to reconstruct the evolution of floral organ initiation and the evolution of developmental processes that give rise to floral zygomorphy on a phylogenetic framework. Results indicate that MAB petal initiation is the most common; it occupies 86.2% of diversity and represents the ancestral state of floral organ initiation in pentapetalous angiosperms. The MAD petal initiation evolved 28 times independently from the ancestral MAB petal initiation. Among the 34 independent originations of floral zygomorphy, 76.5% of these clades represent MAB petal initiation, among which only 47% of the clades result a CV petal in mature flowers. The discrepancy is explained by the existence of developmental processes that result in floral zygomorphy along oblique planes of floral symmetry in addition to along the medial plane. Findings suggest that although the early floral organ initiation plays a constraining role to the evolution of patterns of floral zygomorphy, the constraint diverges along phylogenetically distantly related groups that allow the independent originations of floral zygomorphy through distinct development processes in pentapetalous angiosperms. In additional study, the butterfly-like flowers of Schizanthus are adapted to pollination by bees, hummingbirds, and moths. I investigated the genetic basis of the zygomorphic corolla, for which development is key to the explosive pollen release mechanism found in the species of Schizanthus adapted to bee pollinators. I examined differential gene expression profiles across the zygomorphic corolla of Schizanthus pinnatus, a bee-pollinated species, by analyzing RNA transcriptome sequencing (RNA- seq). Data indicated that CYC2 is not expressed in the zygomorphic corolla of Sc. pinnatus, suggesting CYC2 is not involved in the development of floral zygomorphy in Schizanthus (Solanaceae). The data also indicated that a number of genes are differentially expressed across the corolla.

Brain disorders present a significant burden on affected individuals, their families and society at large. Existing diagnostic tests suffer from a lack of genetic biomarkers, particularly for substance use disorders, such as alcohol dependence (AD). Numerous studies have demonstrated that AD has a genetic heritability of 40-60%. The existing genetics literature of AD has primarily focused on linkage analyses in small family cohorts and more recently on genome-wide association analyses (GWAS) in large case-control cohorts, fueled by rapid advances in next generation sequencing (NGS). Numerous AD-associated genomic variations are present at a common frequency in the general population, making these variants of public health significance. However, known AD-associated variants explain only a fraction of the expected heritability. In this dissertation, we demonstrate that systems biology applications that integrate evolutionary genomics, rare variants and structural variation can dissect the genetic architecture of AD and elucidate its heritability. We identified several complex human diseases, including AD and other brain disorders, as potential targets of natural selection forces in diverse world populations. Further evidence of natural selection forces affecting AD was revealed when we identified an association between eye color, a trait under strong selection, and AD. These findings provide strong support for conducting GWAS on brain disorder phenotypes. However, with the ever-increasing abundance of rare genomic variants and large cohorts of multi-ethnic samples, population stratification becomes a serious confounding factor for GWAS. To address this problem, we designed a novel approach to identify ancestry informative single nucleotide polymorphisms (SNPs) for population stratification adjustment in association analyses. Furthermore, to leverage untyped variants from genotyping arrays – particularly rare variants – for GWAS and meta-analysis through rapid imputation, we designed a tool that converts genotype definitions across various array platforms. To further elucidate the genetic heritability of brain disorders, we designed approaches aimed at identifying Copy Number Variations (CNVs) and viral insertions into the human genome. We conducted the first CNV-based whole genome meta-analysis for AD. We also designed an integrated approach to estimate the sensitivity of NGS-based methods of viral insertion detection. For the first time in the literature, we identified herpesvirus in NGS data from an Alzheimer’s disease brain sample. The work in this dissertation represents a three-faceted advance in our understanding of brain disease etiology: 1) evolutionary genomic insights, 2) novel resources and tools to leverage rare variants, and 3) the discovery of disease-associated structural genomic aberrations. Our findings have broad implications on the genetics of complex human disease and hold promise for delivering clinically useful knowledge and resources.

La diffusion du séquençage haut débit (ou NGS pour Next Generation Sequencing) représente un tel changement d’échelle par rapport aux méthodes classiques de séquençage que les indications et l’organisation du diagnostic moléculaire s’en trouvent profondément modifiées. Le NGS permet à la fois de raccourcir le temps d’analyse et de rendu de résultat et d'élargir considérablement le nombre de gènes testés. Il promet donc d’augmenter la proportion de diagnostics posés et de faciliter l'identification de nouveaux variants et de nouveaux gènes impliqués en pathologie. Cependant dans tous les cas, il génère une quantité de données importante, données qui doivent être analysées et interprétées à l’aide d’outils bioinformatiques spécifiques.Dans la première partie de ce travail, les stratégies existantes ainsi que les difficultés et les enjeux du séquençage haut débit pour le diagnostic moléculaire des maladies génétiques sont discutés. Dans la deuxième partie, la mise en place et la validation technique de cette approche diagnostique sont décrites au sein du laboratoire de Génétique Moléculaire de la Timone à Marseille et illustrées par trois exemples concrets de diagnostics moléculaires posés grâce à la technique de séquençage à haut débit. Dans le domaine spécifique des maladies rares, ces nouvelles technologies sont porteuses d’un réel espoir pour les patients atteints de maladie génétique, permettant d'améliorer globalement leur prise en charge et d'accélérer les progrès dans le domaine de la recherche. / The diffusion of Next Generation Sequencing (NGS) technologies induces an important change that modifies molecular diagnostics indications and prompts laboratories to re-think their diagnostic strategies, up-to-now based on Sanger sequencing routine. Several high throughput approaches are available from the sequencing of a gene panel, to a whole exome, or even a whole genome. In all cases, a tremendous amount of data are generated, that have to be filtered, interpreted and analyzed by the use of powerful bioinformatics tools.In part 1, existing strategies and the difficulties and challenges of high-throughput sequencing for molecular diagnosis in genetic diseases are discussed. In part 2, the set up and the technical validation of this diagnostic approach in the Molecular Genetics’ Laboratory of the Timone Hospital in Marseille is presented and illustrated by 3 examples of complex diagnostics solved thanks to NGS. NGS promises to shorten significantly the time of analysis and results reporting, and to expand the number of tested genes. It also promises to increase the proportion of positive diagnoses. Finally, the NGS can identify new variants and new genes involved in human pathology, thus will globally improve patient clinical care.

Massively parallel high throughput sequencers are transforming the scientific research by reducing the cost and time necessary to sequence entire genomes. The goal of this project is to produce preliminary genome assemblies of calliphorid flies using Life Technologies’ Ion Torrent sequencing and Illumina’s MiSeq sequencing. I located, assembled, and annotated a novel mitochondrial genome for one such fly, the little studied Chrysomya pacifica that is central to one hypothesis about blow fly evolution. With sequencing data from Chrysomya megacephala, its forensically relevant sister species, much insight can be gained by alignments, sequence and protein analysis, and many more tools within the CLC Genomics Workbench software program. I present these analyses here of these recently diverged species.

Le lymphome est un cancer des ganglions lymphatiques, lieu de prolifération et différenciation des cellules immunitaires en particulier des lymphocytes B qui sont des cellules productrices d'anticorps. Les lymphomes résultent de l’accumulation de mutations génétiques dans le génome d’une cellule B normale contribuant à la transformation en cellule B maligne. Cette cellule B dite transformée prolifère alors pour engendrer un clone de cellules B malignes qui s’accumulent au niveau des ganglions lymphatiques, formant alors une tumeur appelée ‘lymphome B’ (pour cancer lymphoïde issu de la transformation maligne des lymphocytes B). Le ganglion lymphatique normal a une structure histologique qui se décompose ainsi, du centre à la périphérie : le centre germinatif, la zone du manteau et la zone marginale. Les lymphomes B sont classés en différents sous-types histologiques en fonction de leur origine topographique au niveau du ganglion lymphatique et de leurs caractéristiques bio cliniques spécifiques. Parmi ces sous-types, une forme particulièrement agressive peut être distinguée : le lymphome à cellules du manteau. Ce sous-type de lymphome est caractérisé par des rechutes successives et une survie qui est généralement courte (médiane de survie de 4 à 5 ans) même si certains patients, avec des formes plus indolentes de lymphomes à cellules du manteau, présentent des survies prolongées (médiane de survie de 7 à 10 ans). Des biomarqueurs prédictifs de la courte survie manquent aujourd’hui, ce qui rend difficile la prise en charge optimisée des patients. Ce projet s’intéresse à cette question. Plus précisément, nous nous proposons de rechercher des mutations génétiques associées à la résistance thérapeutique. Notre approche sera basée sur le séquençage ciblé à haut débit du génome de cellules B tumorales issus de patients présentant des cas classiques du lymphomes à cellules du manteau mais aussi des cas plus particuliers comme ceux présentant des résistances thérapeutiques précoces, par exemple. Par cette approche dite de ‘cartographie’ à l’échelle du génome, nous espérons identifier des nouveaux prédicteurs moléculaires de la survie chez ces patients atteints de lymphome à cellules du manteau et également apporter de nouvelles connaissances dans l’interconnexion entre la génétique et l’épigénétique dans cette maladie. / Lymphoma is a cancer of the lymph nodes which are organs in which immune cells, particularly the antibody producing B cells, proliferate and differentiate before circulating in the blood and tissues to fight infection. B cell lymphoid cancers – ‘B cell lymphoma’ arise as a consequence of the occurrence of gene mutations in B cells. By affecting the functions of key B cell genes, these mutations drive the malignant transformation of the affected B cells which then begin to divide abnormally eventually destroying normal lymph node organization and function. The lymph node is divided into distinct micro-anatomical compartments or zones which are called (from the inner to outer most compartment – germinal centre, mantle zone, and marginal zone). B cell lymphoma classification follows this general organization and classifies tumours depending on the compartment of origin of the particular tumour B cell population. This classification thus defines lymphoma according to a ‘histological subtype’ with defined clinic-biological features. Among these subtypes, mantle cell lymphoma (MCL) is a particularly aggressive form of B lymphoid cancer. This type of lymphoma is characterised by successive relapses and short survival (median is 4 to 5 years), although some patients can show long survival. Predictive biomarkers of this clinical behavior are lacking. This project aims to address this question. More specifically we propose to perform whole ‘exome’ sequencing – i.e. sequencing of all protein coding sections of all known protein coding genes in the genome – of the tumour B cell DNA from patients who show refractory or early relapsing disease compared to patients who show relatively long survival. By doing this genome scale study we hope to identify new gene mutations that can serve as molecular predictors of survival and bring new knowledges in the understanding between genetics and epigenetics in MCL.

Hlavným cieľom tejto práce je návrh nového algoritmu k spracovaniu unikátnych molekulárnych indexov bez mapovania na referenčný genóm. O tieto náhodné oligonukleotidové sekvencie neustále vzrastá záujem, pretože uľahčujú rozpoznávať PCR chyby a skresľovanie údajov. Keďže používanie technológií sekvenovania novej generácie neustále rastie, je vynaložené veľké úsilie vyvíjať nástroje pre analýzu produkovaných dát. V súčasnosti sú nástroje na riešenie týchto chýb relatívne časovo náročné a zložité z dôvodu výpočtovo náročného zarovnania. Najdôležitejšie obmedzenie týchto nástrojov spočíva v skutočnosti, že pri spracovávaní duplikátov sú povolené multi-mapované čítania. Tieto čítania sú zvyčajne ignorované, čo môže viesť k zníženiu kvantitatívnej presnosti a spôsobiť zavádzajúcu interpretáciu výsledkov daného sekvenovania. V snahe vyriešiť tento problém je v tejto práci uvedený nový prístup, ktorý umožňuje odhad absolútneho počtu jedinečných molekúl s relatívne rýchlym a spoľahlivým spôsobom.

Hlavním cílem této diplomové práce je analýza diferenciální exprese genů na základě negativního binomického modelu. Úvodní část je věnována teoretickému základu, pojednává o sekvenování RNA, sekvenování nové generace, výhodách a možném využití, formátu fastQ aj. Následující část už se zabývá samotnou praktickou částí, zde byl vybrán vhodný set genů, které budou později analyzovány a příslušná data byla stažena. Tato data byla zarovnána k lidskému genomu verze 37 Burrowsovou-Wheelerovou transformací s využitím bowtie mapovače, byly tak vytvořeny soubory ve formátu SAM. Toto soubory dat byly později setříděny pomocí nástroje SAMtools. Následně byly v programovém prostředí Matlab (verze R2013b) vytvořeny anotované objekty genů s využitím služby Ensembl´s BioMart. Dále byla určena genová exprese a byly odhadnuty faktory velikosti knihovny. Na závěr byly odhadnuty parametry negativního binomického rozložení a byla vyhodnocena diferenciální exprese genů.