My dissertation comprises three related topics of evolutionary and computational biology, which correspond to the three Chapters. Chapter 1 focuses on tempo and mode of evolution in visual genes, namely opsins, via duplication events and subsequent molecular adaptation in Odonata (dragonflies and damselflies). Gene duplication plays a central role in adaptation to novel environments by providing new genetic material for functional divergence and evolution of biological complexity. Odonata have the largest opsin repertoire of any insect currently known. In particular our results suggest that both the blue sensitive (BS) and long-wave sensitive (LWS) opsin classes were subjected to strong positive selection that greatly weakens after multiple duplication events, a pattern that is consistent with the permanent heterozygote model. Due to the immense interspecific variation and duplicability potential of opsin genes among odonates, they represent a unique model system to test hypotheses regarding opsin gene duplication and diversification at the molecular level. Chapter 2 primarily focuses on reconstruction of the phylogenetic backbone of Odonata using RNA-seq data. In order to reconstruct the evolutionary history of Odonata, we performed comprehensive phylotranscriptomic analyses of 83 species covering 75% of all extant odonate families. Using maximum likelihood, Bayesian, coalescent-based and alignment free tree inference frameworks we were able to test, refine and resolve previously controversial relationships within the order. In particular, we confirmed the monophyly of Zygoptera, recovered Gomphidae and Petaluridae as sister groups with high confidence and identified Calopterygoidea as monophyletic. Fossil calibration coupled with diversification analyses provided insight into key events that influenced the evolution of Odonata. Specifically, we determined that there was a possible mass extinction of ancient odonate diversity during the P-Tr crisis and a single odonate lineage persisted following this extinction event. Lastly, Chapter 3 focuses on identification of erroneously assigned sequence homology using the intelligent agents of machine learning techniques. Accurate detection of homologous relationships of biological sequences (DNA or amino acid) amongst organisms is an important and often difficult task that is essential to various evolutionary studies, ranging from building phylogenies to predicting functional gene annotations. We developed biologically informative features that can be extracted from multiple sequence alignments of putative homologous genes (orthologs and paralogs) and further utilized in context of guided experimentation to verify false positive outcomes.

Malgré l’usage de médicaments antipsychotiques depuis plusieurs décennies, leur mécanisme d’action précis, autre que leur interaction avec les récepteurs dopaminergiques et sérotoninergiques, demeure peu connu. Nur77 (Nr4a1 ou NGFI-B) est un facteur de transcription de la famille des récepteurs nucléaires associé aux effets des antipsychotiques. Ceci étant dit, le mécanisme d’action de Nur77 est également peu connu. Afin de mieux comprendre les éléments impliqués avec les antipsychotiques et l’activité de Nur77, nous avons comparé les niveaux de transcrits dans le striatum suite à un traitement avec l’halopéridol chez des rats sauvages et déficients en Nur77 à l’aide de la technique de séquençage à haut débit (RNAseq) et d’une analyse bio-informatique. L’halopéridol et Nur77 ont modulé d’importants groupes de gènes associés avec la signalisation des récepteurs dopaminergiques et la synapse glutamatergique. L’analyse a révélé des modulations de gènes clés reliés à la signalisation des protéines G. Parmi les transcrits modulés significativement chez les rats traités avec halopéridol et ceux déficients en Nur77, la dual specificity phosphatase 5 (Dusp5) représente un nouveau candidat d’intérêt. En effet, nous avons confirmé que les niveaux d’ARNm et protéiques de Dusp5 dans le striatum sont associés aux mouvements involontaires anormaux (dyskinésie) dans un modèle de primates non-humains traités chroniquement avec halopéridol. Cette analyse transcriptomique a démontré des altérations rapides et importantes d’éléments impliqués dans la signalisation des protéines G par l’halopéridol, et a permis d’identifier, pour la première fois, une expression de Dusp5 dépendante de Nur77 en tant que nouvelle composante reliée avec la dyskinésie tardive. / Despite antipsychotic drugs being used for several decades, their precise mechanism of action remains elusive. Nur77 (Nr4a1 or NGFI-B) is a transcription factor of the nuclear receptor family associated with antipsychotic drug effects. However, the mechanism of action of Nur77 is also not well understood. To better understand the signaling components implicated with antipsychotic drug use and Nur77 activity, we compared striatal gene transcripts following haloperidol in wild-type and Nur77-deficient rats using Next Generation RNA Sequencing (RNAseq) and a bioinformatics analysis. Haloperidol and Nur77 modulated important subsets of striatal genes associated with dopamine receptor signaling and glutamate synapses. The analysis revealed modulations of key components of G protein signaling that are consistent with a rapid adaptation of striatal cells that may partially explain long-term haloperidol-induced dopamine D2 receptor upregulation. Amongst significantly modulated transcripts in rats treated with haloperidol and rats deficient in Nur77, dual specificity phosphatase 5 (Dusp5) represents a new and very interesting candidate. Indeed, we confirmed that striatal Dusp5 mRNA and protein levels were associated with abnormal involuntary movements (dyskinesia) in non-human primates chronically exposed to haloperidol. This transcriptomic analysis showed important haloperidol-induced G protein-coupled receptor signaling alterations that may support a regulatory role of Nur77 in dopamine D2 receptor signaling pathways and identified, for the first time, a putative Nur77-dependent expression of Dusp5 as a new signaling component for antipsychotic drug-induced tardive dyskinesia.

No description available.

Affinity ligands are important molecules used in affinity chromatography for purification of significant substances from complex mixtures. To find affinity ligands specific to important target molecules could be a challenging process. Cytiva uses the powerful phage display technique to find new promising affinity ligands. The phage display technique is a method run in several enrichment cycles. When developing new affinity ligands, a protein scaffold library with a diversity of up to 1010-1011 different protein scaffold variants is run through the enrichment cycles.  The result from the phage display rounds is screened for target molecule binding followed by sequencing, usually with one of the conventional screening methods ELISA or Biacore followed by Sanger sequencing. However, the throughput of these analyses are unfortunately very low, often with only a few hundred screened clones. Therefore, Next Generation Sequencing or NGS, has become an increasingly popular screening method for phage display libraries which generates millions of sequences from each phage display round. This creates a need for a robust data analysis pipeline to be able to interpret the large amounts of data.  In this project, a pipeline for analysis of NGS data of phage displayed libraries has been developed at Cytiva. Cytiva uses NGS as one of their screening methods of phage displayed protein libraries because of the high throughput compared to the conventional screening methods. The purpose is to find new affinity ligands for purification of essential substances used in drugs.  The pipeline has been created using the object-oriented programming language R and consists of several analyses covering the most important steps to be able to find promising results from the NGS data. With the developed pipeline the user can analyze the data on both DNA and protein sequence level and per position residue breakdown, as well as filter the data based on specific amino acids and positions. This gives a robust and thorough analysis which can lead to promising results that can be used in the development of novel affinity ligands for future purification products.

Droplet-based microfluidics is a powerful tool permitting massive-scale single-cell analysis in pico-/nano-liter water-in-oil droplets. It has been integrated into various library preparation techniques to accomplish high-throughput scRNA-seq, scDNA-seq, scATAC-seq, scChIP-seq, as well as scMulti-omics-seq. These advanced technologies have been providing unique and novel insights into both normal differentiation and disease development at single-cell level. In this thesis, we develop four new droplet-based tools for single-cell omics profiling. First, the developed Drop-BS is the first droplet-based platform to construct single-cell bisulfite sequencing libraries for DNA methylome profiling and allows production of BS library of 2,000-10,000 single cells within 2 d. We applied the technology to separately profile mixed cell lines, mouse brain tissues, and human brain tissues to reveal cell type heterogeneity. Second, the new Drop-ChIP platform only requires two steps of droplet generation to achieve multiple steps of reactions in droplets such as single-cell lysis, chromatin fragmentation, ChIP, and barcoding. Third, we aim to establish a droplet-based platform to accomplish high-throughput full-length RNA-seq (Drop-full-seq), which both current tube-based and droplet-based methods cannot realize. Last, we constructed an in-house droplet-based tool to assist single-cell ATAC-seq library preparation (Drop-ATAC), which provided a low-cost and facile protocol to conduct scATAC-seq in laboratories without the expensive instrument. / Doctor of Philosophy / Microfluidics is a collection of techniques to manipulate fluids in the micrometer scale. One of microfluidic techniques is called "droplet-based microfluidics". It can manipulate (i.e., generate, merge, sort, split, etc) pico-/nano-liter of water-in-oil droplets. First, since the water phase is separated by the continuous oil phase, these droplets are discrete and individual reactors. Second, droplet-based microfluidics can achieve highly parallel manipulation of thousands to millions of droplets. These two advantages make droplet-based microfluidics an ideal tool to perform single-cell assays. Over the past 10 years, various droplet-based platforms have been developed to study single-cell transcriptome, genome, epigenome, as well as multi-ome. To expand droplet-based tools for single-cell analysis, we aim to develop four novel platforms in this thesis. First, Drop-BS, by integrating droplet generation and droplet fusion techniques, can achieve high-throughput single-cell bisulfite sequencing library preparation. It can generate 10,000 single-cell BS libraries within 2 days which is difficult to achieve for conventional library preparation in tubes/microwells. Second, we developed a novel and facile Drop-ChIP platform to prepare single-cell ChIP-seq library. It is easy to operate since it only requires two steps of droplet generation. It also generates higher quality of data compared to previous work. In addition, we are working on the development and characterization of the other two droplet-based tools to achieve full-length single-cell RNA-seq and single-cell ATAC-seq.

Psychiatric disorders are believed to result from the combination of genetic predisposition and many environmental triggers. While the large number of disease-associated genetic variations have been recognized by previous genome-wide association studies (GWAS), the role of epigenetic mechanisms that mediate the effects of environmental factors on CNS gene activity in the etiology of most mental illnesses is still largely unclear. A growing body of evidence suggested that the abnormalities (changes in gene expression, formation of neural circuits, and behavior) involved in most psychiatric syndromes are preserved by epigenetic modifications identified in several specific brain regions. In this thesis, we developed the second generation of one of our microfluidic technologies (MOWChIP-seq) and used it to profile genome-wide histone modifications in three mental illness-related biological studies: the effect of psychedelics in mice, schizophrenia, and the effect of maternal immune activation in mice offspring. The second generation of MOWChIP-seq was designed to generate histone modification profiles from as few as 100 cells per assay with a throughput as high as eight assays in each run. Then, we applied the new MOWChIP-seq and SMART-seq2 to profile the histone modification H3K27ac and transcriptome, respectively, using NeuN+ neuronal nuclei from the mouse frontal cortex after a single dose of psychedelic administration. The epigenomic and transcriptomic changes induced by 2,5-Dimethoxy-4-iodoamphetamine (DOI), a subtype of psychedelics, in mouse neuronal nuclei at various time points suggest that the long-lasting effects of the psychedelic are more closely related to epigenomic alterations than the changes in transcriptomic patterns. Next, we comprehensively characterized epigenomic and transcriptomic features from the frontal cortex of 29 individuals with schizophrenia and 29 individually matched controls (gender and age). We found that schizophrenia subjects exhibited thousands of neuronal vs. glial epigenetic differences at regions that included several susceptibility genetic loci, such as NRXN1, RGS4 and GRIN3A. Finally, we investigated the epigenetic and transcriptomic alterations induced by the maternal immune activation (MIA) in mice offspring's frontal cortex. Pregnant mice were injected with influenza virus at GD 9.5 and the frontal cortex from mice pups (10 weeks old) were examined later. The results offered us some insights into the contribution of MIA to the etiology of some mental disorders, like schizophrenia and autism. / Doctor of Philosophy / While this field is still in its early stage, the epigenetic studies of mental disorders present promise to expand our understanding about how environmental stimulates, interacting with genetic factors, contribute to the etiology of various psychiatric syndromes, like major depression and schizophrenia. Previous clinical trials suggested that psychedelics may represent a promising long-lasting treatment for patients with depression and other psychiatric conditions. These research presented the therapeutic potential of psychedelic compounds for treating major depression and demonstrated the capability of psychedelics in increasing dendritic density and stimulating synapse formation. However, the molecular mechanism mediating the clinical effectiveness of psychedelics remain largely unexplored. Our study revealed that epigenomic-driven changes in synaptic plasticity sustain psychedelics' long-lasting antidepressant action. Another serious mental illness is schizophrenia, which could affect how an individual feels, thinks, and behaves. Like most other mental disorders, schizophrenia results from a combination of genetic and environmental causes. Epigenetic marks allow a dynamic impact of environmental factors, including antipsychotic medications, on the access to genes and regulatory elements. Despite this, no study so far has profiled cell-type-specific genome-wide histone modifications in postmortem brain samples from schizophrenia subjects or the effect of antipsychotic treatment on such epigenetic marks. Here we show the first comprehensive epigenomic characterization of the frontal cortex of 29 individuals with schizophrenia and 29 matched controls. The process of brain development is surprisingly sensitive to a lot of environmental insults. Epidemiological studies have recognized maternal immune activation as a risk factor that may change the normal developmental trajectory of the fetal brain and increase the odds of developing a range of psychiatric disorders, including schizophrenia and autism, in its lifetime. Given the prevalence of the coronavirus, uncovering the molecular mechanism underlie the phenotypic alterations has become more urgent than before, for both prevention and treatment.

Epigenetic modifications, such as DNA methylation and histone modifications, play important roles in gene expression and regulation, and are highly involved in cellular processes such as stem cell pluripotency/differentiation and tumorigenesis. Chromatin immunoprecipitation (ChIP) is the technique of choice for examining in vivo DNA-protein interactions and has been a great tool for studying epigenetic mechanisms. However, conventional ChIP assays require millions of cells for tests and are not practical for examination of samples from lab animals and patients. Automated microfluidic chips offer the advantage to handle small sample sizes and facilitate rapid reaction. They also eliminate cumbersome manual handling. In this report, I will talk about three different projects that utilized microfluidic immunoprecipitation followed by next genereation sequencing technologies to enable low input and high through epigenomics profiling. First, I examined RNA polymerase II transcriptional regulation with microfluidic chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq) assays. Second, I probed the temporal dynamics in the DNA methylome during cancer development using a transgenic mouse model with microfluidic methylated DNA immunoprecipitation followed by next generation sequencing (MeDIP-seq) assays. Third, I explored negative enrichment of circulating tumor cells (CTCs) followed by microfluidic ChIP-seq technology for studying temporal dynamic histone modification (H3K4me3) of patient-derived tumor xenograft on an immunodeficient mouse model during the course of cancer metastasis. In the first study, I adapted microfluidic ChIP-seq devices to achieve ultrahigh sensitivity to study Pol2 transcriptional regulation from scarce cell samples. I dramatically increased the assay sensitivity to an unprecedented level (~50 K cells for pol2 ChIP-seq). Importantly, this is three orders of magnitude more sensitive than the prevailing pol2 ChIP-seq assays. I showed that MNase digestion provided better ChIP-seq signal than sonication, and two-steps fixation with MNase digestion provided the best ChIP-seq quality followed by one-step fixation with MNase digestion, and lastly, no fixation with MNase digestion. In the second study, I probed dynamic epigenomic changes during tumorigenesis using mice often require profiling epigenomes using a tiny quantity of tissue samples. Conventional epigenomic tests do not support such analysis due to the large amount of materials required by these assays. In this study, I developed an ultrasensitive microfluidics-based methylated DNA immunoprecipitation followed by next-generation sequencing (MeDIP-seq) technology for profiling methylomes using as little as 0.5 ng DNA (or ~100 cells) with 1.5 h on-chip process for immunoprecipitation. This technology enabled me to examine genome-wide DNA methylation in a C3(1)/SV40 T-antigen transgenic mouse model during different stages of mammary cancer development. Using this data, I identified differentially methylated regions and their associated genes in different periods of cancer development. Interestingly, the results showed that methylomic features are dynamic and change with tumor developmental stage. In the last study, I developed a negative enrichment of CTCs followed by ultrasensitive microfluidic ChIP-seq technology for profiling histone modification (H3K4Me3) of CTCs to resolve the technical challenges associated with CTC isolation and difficulties related with tools for profiling whole genome histone modification on tiny cell samples. / Ph. D.

RNA-Seq is an increasingly popular approach to transcriptome profiling that uses the capabilities of next generation sequencing technologies and provides better measurement of levels of transcripts and their isoforms. In this thesis, we apply RNA-Seq protocol and transcriptome quantification to estimate gene expression and pathway activity levels. We present a novel method, called IsoDE, for differential gene expression analysis based on bootstrapping. In the first version of IsoDE, we compared the tool against four existing methods: Fisher's exact test, GFOLD, edgeR and Cuffdiff on RNA-Seq datasets generated using three different sequencing technologies, both with and without replicates. We also introduce the second version of IsoDE which runs 10 times faster than the first implementation due to some in-memory processing applied to the underlying gene expression frequencies estimation tool and we also perform more optimization on the analysis. The second part of this thesis presents a set of tools to differentially analyze metabolic pathways from RNA-Seq data. Metabolic pathways are series of chemical reactions occurring within a cell. We focus on two main problems in metabolic pathways differential analysis, namely, differential analysis of their inferred activity level and of their estimated abundance. We validate our approaches through differential expression analysis at the transcripts and genes levels and also through real-time quantitative PCR experiments. In part Four, we present the different packages created or updated in the course of this study. We conclude with our future work plans for further improving IsoDE 2.0.

The aim of the current thesis was the analysis of the genetics of sex determination of farmed fish with sexual dimorphism, using Next Generation Sequencing. Three different species of farmed fish with sex-determining systems of varying complexity were studied. Both full-sibs and more distantly related specimens of Atlantic halibut (Hippoglossus hippoglossus), Nile tilapia (Oreochromis niloticus) and European sea bass (Dicentrarchus labrax) were used for this study. Application of Restriction-site Associated DNA sequencing (RAD-seq) and double digest Restriction-site Associated DNA sequencing (ddRAD-seq), two related techniques based on next generation sequencing, allowed the identification of thousands of Single Nucleotide Polymorphisms (SNPs; > 3,000) for each of the above species. The first SNP-based genetic maps for the above species were constructed during the current study. The first evidence concerning the location of the sex-determining region of Atlantic halibut is provided in this study. In the case of Nile tilapia both novel sex-determining regions and fine mapping of the major sex-determining region are presented. In the study of European sea bass evidence concerning the absence of a major sex-determining gene was provided. Indications of putative sex-determining regions in this species are also provided. The results of the current thesis help to broaden current knowledge concerning sex determination in three important farmed fish. In addition the results of the current thesis have practical applications as well, towards the production of mono-sex stocks of those species for the aquaculture industry.

Récemment, une nouvelle classe d'ADN circulaire extrachromosomique (eccDNA) appelée microADN a été identifiée dans des tissus humains et murins. Ces microADNs ont une longueur de 100 à 400 pb, sont dérivés de régions génomiques non répétitives uniques et présentent un enrichissement au niveau des régions géniques et riches en GC. Bien qu'il ait été proposé qu'ils puissent provenir du métabolisme de l'ARN ou des défauts de réplication, leurs mécanismes de production et leur éventuelle fonctionnalité restent à déterminer. Grâce à l'analyse des microADNs extraits d'une série de 10 lignées cellulaires lymphoblastoïdes humaines (LCL), nous avons confirmé la distribution nonaléatoire des microADNs vers les régions actives du génome. Les microADNs identifiés présentaient des loci d'origine redondants et une périodicité de taille de 190 pb pouvant correspondre à la fragmentation de l'ADN lors de l'apoptose caspase-dépendante. L'apoptose induite de ces LCLs par des drogues chimiothérapeutiques (méthotrexate ou L-asparaginase) a entrainé la modulation de la diversité et de la taille des microADNs, suggérant qu'une partie de ces entités pourrait être des produits résiduels de la mort cellulaire apoptotique. Ainsi, bien que compatible avec l'observation initiale suggérant que les microADNs proviennent d'un processus physiologique normal, ces résultats impliquent une source de production alternative ou complémentaire. / Recently, a new class of extrachromosomal circular DNA (eccDNA) called microDNA was identified in mouse and human tissues. These microDNAs are 100 to 400 bp long, derive from unique nonrepetitive genomic regions and show an enrichment in GC rich and genic sequences. While it has been proposed that they could arise from RNA metabolism or replication defects, their production mechanisms and eventual functionality remain unclear. Through the analysis of microDNAs extracted from a series of 10 human lymphoblastoid cell lines (LCLs), we confirmed the non-random distribution of microDNA towards active regions of the genome. Identified microDNAs showed redundant loci of origin and a size periodicity of 190 bp that matched caspase-dependant DNA fragmentation of apoptotic cells. Strikingly, the chemotherapeutic drug-induced apoptosis (using methotrexate or Lasparaginase) of these LCLs modulated both diversity and size of microDNAs further suggesting that a part of microDNAs could represent circularized by-products of the programmed cell death. Thus, while compatible with the original observation that microDNAs originated from a normal physiological process, these results imply an alternative or complementary source of production.