Complex I (NADH:ubiquinone oxidoreductase) couples electron transfer to proton translocation across the inner mitochondrial membrane, to drive the synthesis of ATP. Its distinctive L-shaped structure comprises 45 subunits, encoded by both the mitochondrial and nuclear genomes, which are assembled by a complicated modular pathway. Complex I genetic defects are the most common cause of mitochondrial disorders and often present in early childhood, with high mortality rates. Recent high-resolution electron cryo-microscopy structures of mammalian complex I provide a foundation for both interpreting biochemical and biomedical data and understanding the catalytic mechanism. First, this thesis explores how the flavin cofactor is inserted into the NADH-binding (N-) domain of complex I. Genetic manipulation of cultured human cells, to starve them of flavin, revealed a hierarchal impact on the mitochondrial flavoproteome. High riboflavin content in the growth media ameliorated observed phenotypes, requiring cell conditioning in low riboflavin conditions. CRISPR knockout of the putative mitochondrial flavin transporter SLC25A32 demonstrated the severe impact of decreased flavin on complexes I and II, and mass spectrometry 'complexome' analyses suggest that the N-domain is still assembled onto complex I in the absence of the flavin. Second, the model organism Yarrowia lipolytica was used to assess the importance of residues in the quinone-binding site of complex I. Three residues with proposed roles in binding the quinone head-group were targeted. One variant was catalytically inactive, while two retained some activity. They showed decreased ability to reduce physiologically-relevant, long chain quinones, but their ability to reduce short-chain analogues was affected less severely. The results suggest a complicated picture in which interactions between the protein and both the hydrophilic quinone head-group and hydrophobic isoprenoid chain contribute to quinone-binding affinity and catalysis. Finally, a model for human complex I, generated from a recent high-resolution structure of mouse complex I, was used to investigate whether the pathogenicity of human variants could be predicted. Structural information on variant residues, including their secondary structure, proximity to key features and surface exposure, was collated and the power of each property to predict pathogenicity investigated. The analysis was then extended to the whole structure, to identify potential pathogenic hotspots in the enzyme, inform future studies of functionally important regions in complex I, and aid the diagnosis of clinically relevant pathogenic variants.

As the volume of transcriptomic data continues to increase, so too does its potential to deepen our understanding of disease; for example, by revealing gene expression patterns shared between diseases. However, key questions remain around the strength of the transcriptomic signal of disease and the identification of meaningful commonalities between datasets, which are addressed in this thesis as follows. The first chapter, Concordance of Microarray Studies of Parkinson's Disease, examines the agreement between differential expression signatures across 33 studies of Parkinson's disease. Comparison of these studies, which cover a range of microarray platforms, tissues, and disease models, reveals a characteristic pattern of differential expression in the most highly-affected tissues in human patients. Using correlation and clustering analyses to measure the representativeness of different study designs to human disease, the work described acts as a guideline for the comparison of microarray studies in the following chapters. In the next chapter, Using Dysregulated Signalling Paths to Understand Disease, gene expression changes are linked on the human signalling network, enabling identification of network regions dysregulated in disease. Applying this method across a large dataset of 141 common and rare diseases identifies dysregulated processes shared between diverse conditions, which relate to known disease- and drug-sharing-relationships. The final chapter, Understanding and Predicting Disease Relationships Through Similarity Fusion, explores the integration of gene expression with other data types - in this case, ontological, phenotypic, literature co-occurrence, genetic, and drug data - to understand relationships between diseases. A similarity fusion approach is proposed to overcome the differences in data type properties between each space, resulting in the identification of novel disease relationships spanning multiple bioinformatic levels. The similarity of disease relationships between each data type is considered, revealing that relationships in differential expression space are distinct from those in other molecular and clinical spaces. In summary, the work described in this thesis sets out a framework for the comparative analysis of transcriptomic data in disease, including the integration of biological networks and other bioinformatic data types, in order to further our knowledge of diseases and the relationships between them.

Mycoplasma hyopneumoniae, Mycoplasma flocculare e Mycoplasma hyorhinis são capazes de aderir e colonizar o trato respiratório de suínos. Enquanto a presença de M. flocculare é considerada assintomática, M. hyopneumoniae e M. hyorhynis são relacionados ao desenvolvimento de patologias. M. hyopneumoniae é o agente etiológico da pneumonia enzoótica suína e M. hyorhynis além dos pulmões pode atingir outros sítios e hospedeiros, estando relacionado a artrites, poliserosites e desenvolvimento de vários tipos de câncer em humanos. Apesar dos avanços tecnológicos na área de genômica, raros são os dados quanto ao papel de M. flocculare no trato respiratório suíno. Além do mais, informações relativas à transcrição gênica nessas espécies são escassas, apesar da importância desses microrganismos. Neste estudo são apresentados os dados da sequência do genoma de uma linhagem de M. flocculare, bem como do genoma de um novo isolado de M. hyopneumoniae. Com estas novas sequências foram realizadas análises de genômica comparativa visando a identificação de características que pudessem explicar os diferentes comportamentos quanto à patogenicidade dessas espécies. Além disso, a análise global dos transcritomas de cada uma das espécies foi realizada e o perfil transcricional entre M. hyopneumoniae, M. flocculare e M. hyorhynis foi analisado comparativamente objetivando identificar características peculiares para cada um dos mapas transcricionais, além de compreender a coordenação do modo de transcrição gênica em Mycoplasma. De um modo geral, as três espécies de Mycoplasma que habitam o trato respiratório suíno possuem grandes semelhanças na composição gênica, assim como na abundância de transcritos. A análise do repertório transcricional, mostra que os genomas são transcritos quase que em sua totalidade, incluindo as regiões intergênicas, nas três espécies. M. hyopneumoniae e M. flocculare apresentam conteúdo gênico e perfil transcricional muito semelhantes. Uma importante diferença encontrada entre estas duas espécies refere-se à presença exclusiva de genes e transcritos de adesinas específicas. M. hyorhynis possui genes e transcritos exclusivos, os quais sabidamente estão relacionados à sua capacidade mutacional, de invasividade e infecção de diferentes sítios. Por fim, a análise comparativa dos genomas, e a obtenção dos mapas transcricionais para M. hyopneumoniae, M. flocculare e M. hyorhynis, foram abordagens que resultaram em um grande número de informações, as quais são importantes para embasamento de futuros estudos de caracterização dos mecanismos moleculares, como os eventos de regulação da transcrição gênica, no gênero Mycoplasma. / Mycoplasma hyopneumoniae, Mycoplasma hyorhinis and Mycoplasma flocculare are able to adhere and to colonize the swine respiratory tract. While M. flocculare presence is virtually assymptomatic, M. hyopneumoniae and M. hyorhynis infections may cause respiratory disease. M. hyopneumoniae is the causative agent of swine enzootic pneumonia and M. hyorhynis may affect the lungs and other sites in a diversity of hosts and has been related to arthritis, poliserosites and to the development of several types of human cancer. Despite genomics technological advances, there are very few data about the possible role of M. flocculare in the swine respiratory tract. Moreover, little information about gene transcription is available in these species, despite the importance of these microorganisms. In this work the genome sequences of M. flocculare and a new isolate of M. hyopneumoniae are presented. A comparative genomic analyzes was performed to identify possible characteristics that may help to explain the different behaviors of these species in the swine respiratory tracts. Furthermore, a transcriptome map of each species was performed and a comparative transcriptional profile analysis between M. hyopneumoniae, M. flocculare and M. hyorhynis was undertaken to identify the exclusive features for each of the transcriptional maps, in addition to understanding the coordination mode of gene transcription in Mycoplasma. In general, the three Mycoplasma species that inhabit the swine respiratory tract have a similar gene composition as well as the abundance of transcripts. The transcriptome maps showed that most of the predicted genes are transcribed from these Mycoplasma genomes, as well as some intergenic regions. M. hyopneumoniae and M. flocculare present very similar gene content and transcriptional profile. However, an important difference between these two species is related to the exclusive presence of genes and transcripts of some specific adhesins. M. hyorhynis presents exclusive genes and transcripts that have been related to its invasiveness, mutation rate and infection of different sites. Finally, the comparative analysis of the genomes and transcriptional maps between M. hyopneumoniae, M. flocculare and M. hyorhynis have resulted in a large amount of information, which are important for future studies of the molecular characterization, as transcriptional regulation in the Mycoplasma spp.

Listeria monocytogenes é um bastonete Gram-positivo, anaeróbio facultativo, psicrotrófico, patogênico a humanos e transmitido por alimentos. É causador da listeriose, doença severa que acomete grupos de risco específicos, tais como idosos, imunocomprometidos, gestantes, crianças e recém-nascidos. Neste trabalho foi investigada a expressão diferencial de L. monocytogenes na presença de óleo essencial extraído de Baccharis psiadioides, planta da família Asteraceae popularmente chamada de “alecrim-do-campo”, “vassoura” ou “erva formiga”, utilizada pela população como planta medicinal. Além disso, os genomas de dois diferentes sorotipos de L. monocytogenes, frequentemente associados a surtos de listeriose, foram sequenciados através de plataforma MiSeq Illumina, sequências estas depositadas no GenBank, e comparados com genomas de referência. Anteriormente à execução das análises genômica e transcriptômica, foi determinada a composição do óleo essencial extraído de B. psiadioides utilizado nos experimentos, através de cromatografia gasosa com espectrômetro de massa (GC – MS), a qual demonstrou uma maior quantidade de β-pineno na fração composta majoritariamente por monoterpenos, composto este frequentemente encontrado em plantas medicinais aromáticas e apontado como um dos responsáveis pelo potencial antimicrobiano das mesmas. Os demais resultados obtidos no presente trabalho indicam que o óleo essencial testado apresenta potencial ação bacteriostática na concentração estudada, sendo que genes relacionados à virulência do microrganismo foram menos transcritos na sua presença, ao contrário do que foi observado para genes de resposta ao estresse, que apresentaram maiores níveis de transcrição nesta condição. A comparação genômica entre os genomas bacterianos sequenciados neste trabalho e as cepas referência sugere um maior número de proteínas expressas em L. monocytogenes do sorotipo 4b relacionadas à defesa e metabolismo do microrganismo, indicando mecanismos que podem estar envolvidos com a capacidade deste sorotipo estar mais envolvido nos casos humanos de listeriose. / LLiisstteerriiaa mmoonnooccyyttooggeenneess is a Gram-positive rod-shaped microorganism, facultative anaerobic, psychrotrophic, pathogenic to humans and transmitted by food. It causes listeriosis, a severe disease that affects specific risk groups such as elderly, immunocompromised, pregnant women, children and newborns. In this study, differential expression of LL.. mmoonnooccyyttooggeenneess in the presence of essential oil extracted from BBaacccchhaarriiss ppssiiaaddiiooiiddeess, a plant from AAsstteerraacceeaaee family popularly named as "alecrim-do-campo", "vassoura" or "erva formiga" used by population as a medicinal plant, was investigated. In addition, the genomes of two different LL.. mmoonnooccyyttooggeenneess serotypes, often associated with listeriosis outbreaks, were sequenced through the MiSeq Illumina platform. These sequences were deposited in GenBank and compared with reference genomes. Prior to the execution of genomic and transcriptomic analyzes, composition of the essential oil extracted from BB.. ppssiiaaddiiooiiddeess used in the experiments was determined by gas chromatography with mass spectrometer (GC-MS), which demonstrated a higher amount of β-pinene in the fraction composed mainly by monoterpenes. This compound is often found in aromatic medicinal plants and also pointed as one of those responsible for their antimicrobial potential. The other results obtained in the present study indicate that the essential oil tested has a potential bacteriostatic activity at the concentration studied, and genes related to the virulence of the microorganism were less transcribed in its presence, contrary to what was observed for stress response genes, which presented higher transcription levels on that condition. Comparative genomics between the bacterial genomes sequenced in this work and the reference strains suggests a higher number of proteins expressed in LL.. mmoonnooccyyttooggeenneess serotype 4b related to the defense and metabolism of the microorganism, indicating mechanisms that may be involved with the greater ability of this serotype to cause human listeriosis.

Type 2 Diabetes (T2D) affects over 415 million people globally and is characterised by cellular stresses including: poor glucose homeostasis, dyslipidaemia, inflammation, hypoxia and ER stress. Studies in mice have shown that exposure to these stresses influences beta cell differentiation status as well as cell survival and may explain the extent of beta cell mass loss that is seen in the disease. To date, studies of altered beta cell differentiation have largely been confined to murine models. I used the EndoC-bH1 human beta cell line, along with human pancreatic tissue sections, to better characterise this mechanism in human disease. To elucidate these mechanisms, I firstly established a humanised version of cell culture techniques for the EndoC βH1 cell model and assessed the influence on cell function. Secondly, I evaluated the effects of the diabetic microenvironment on beta cell differentiation and gene expression patterns. Finally, I investigated whether a diabetomimetic microenvironment induced differences in microRNA regulation in the cells. I found that the humanised EndoC-βH1 culture techniques improved glucose sensitive insulin release in the cell model. EndoC-βH1 cells exposed to a Diabetic microenvironment showed some degree of transdifferentiation and this may be due to dysregulation of splicing factor expression. These effects may be compounded by altered microRNA regulation in response to these cell stresses. These data suggest that altered gene regulation caused by a diabetic microenvironment may alter gene regulation to produce a reversible delta-like phenotype in human beta cells.

Le but de mon projet est de créer et d’appliquer des outils pour l’analyse de la biologie des systèmes musculaires en utilisant différentes données OMICS. Ce projet s’intéresse plus particulièrement à la dysferlinopathie due la déficience d’une protéine appelée dysferline qui est exprimée principalement dans les muscles squelettiques et cardiaque. La perte du dysferline due à la mutation (autosomique-récessive) du gène DYSF entraîne une dystrophie musculaire progressive (LGMD2B, MM, DMAT). Nous avons déjà développé des outils bio-informatiques qui peuvent être utilisés pour l’analyse fonctionnelle de données OMICS, relative à la dyspherlinopathie. Ces derniers incluent le test dit «gene set enrichment analysis», test comparant les profils OMICS d’intérêts aux données OMICS musculaires préalablement publiées ; et l’analyse des réseaux impliquant les diffèrent(e)s protéines et transcrits entre eux/elles. Ainsi, nous avons analysé des centaines de données omiques publiées provenant d’archives publiques. Les outils informatiques que nous avons développés sont CellWhere et MyoMiner. CellWhere est un outil facile à utiliser, permettant de visualiser sur un graphe interactif à la fois les interactions protéine-protéine et la localisation subcellulaire des protéines. Myominer est une base de données spécialisée dans le tissu et les cellules musculaires, et qui fournit une analyse de co-expression, aussi bien dans les tissus sains que pathologiques. Ces outils seront utilisés dans l'analyse et l'interprétation de données transcriptomiques pour les dyspherlinopathies mais également les autres pathologies neuromusculaires. / The aim of this project was to build and apply tools for the analysis of muscle omics data, with a focus on Dysferlin deficiency. This protein is expressed mainly in skeletal and cardiac muscles, and its loss due to mutation (autosomal-recessive) of the DYSF gene, results in a progressive muscular dystrophy (Limb Girdle Muscular Dystrophy type 2B (LGMD2B), Miyoshi myopathy and distal myopathy with tibialis anterior onset (DMAT)). We have developed various tools and pipelines that can be applied towards a bioinformatics functional analysis of omics data in muscular dystrophies and neuromuscular disorders. These include: tests for enrichment of gene sets derived from previously published muscle microarray data and networking analysis of functional associations between altered transcripts/proteins. To accomplish this, we analyzed hundreds of published omics data from public repositories. The tools we developed are called CellWhere and MyoMiner. CellWhere is a user-friendly tool that combines protein-protein interactions and protein subcellular localizations on an interactive graphical display (https://cellwhere-myo.rhcloud.com). MyoMiner is a muscle cell- and tissue-specific database that provides co-expression analyses in both normal and pathological tissues. Many gene co-expression databases already exist and are used broadly by researchers, but MyoMiner is the first muscle-specific tool of its kind (https://myominer-myo.rhcloud.com). These tools will be used in the analysis and interpretation of transcriptomics data from dysferlinopathic muscle and other neuromuscular conditions and will be important to understand the molecular mechanisms underlying these pathologies.

La vigne (Vitis vinifera) possède un métabolisme secondaire particulièrement riche donnant naissance à une large palette de molécules dont certaines sont impliquées dans les défenses contre les pathogènes et d'autres dans la grande diversité d’arômes qui fait la renommée des vins. L’analyse de la séquence de référence du génome de la vigne a permis de mettre en évidence une remarquable expansion de certaines familles de gènes liés au métabolisme secondaire par rapport aux autres plantes. Dans ce travail, j'ai étudié les familles gènes codant pour les cytochromes P450, dont certains sont impliqués dans la production d’arômes, les gènes codant pour les stilbènes synthases (STS), les endo-β-1,3-glucanases et les gènes de résistance de type NBS impliqués dans les défenses de la vigne. Ma thèse vise à proposer des hypothèses expliquant l’organisation structurale de ces familles de gènes et ainsi à mieux comprendre pourquoi certaines familles présentent une amplification dans le génome de la vigne. Des approches bioinformatiques ont été utilisées afin d’étudier ces différentes familles de gènes. Les gènes cytochromes P450 et gènes R de type NBS ont tout d'abord été annotés de manière manuelle dans le génome de référence de la vigne. L’expression des gènes endo-β-1,3- glucanases, STS et cytochromes P450 a été analysée en utilisant une approche transcriptomique à grande échelle. Pour ce faire, un outil a été développé durant cette thèse pour estimer le niveau d’expression des gènes à partir de données RNA-Seq disponibles dans les banques de données publiques. Parallèlement, des données de reséquençage d’ADN de 56 cépages et espèces de vigne ont été analysées, afin de déterminer les variations structurales de type CNV au sein des familles de gènes à domaine NBS et de gènes STS. Ces différents travaux ont permis de montrer que l’amplification des familles de gènes étudiées n’est pas spécifique du génome de référence mais est retrouvée dans l'ensemble du genre Vitis, mais également de mettre en évidence des variations structurales au sein des différents génomes étudiés. L'analyse de la famille STS a montré que ces gènes sont organisés en blocs de duplication, et que les gènes plus conservés sont aussi les plus exprimés. Nous avons également montré que les gènes à domaine NBS sont organisés en cluster, dont certains sont particulièrement soumis à variation. Ces travaux contribuent à une meilleure connaissance de facteurs de défense efficaces et durables ainsi que des gènes impliqués dans la synthèse d’arômes dans la vigne. Ces connaissances pourront bénéficier aux programmes de création variétale mis en œuvre à l’INRA de Colmar. / Grapevine (Vitis vinifera) has a particularly rich secondary metabolism, giving rise to a wide range of molecules, some of which are involved in defences against pathogens and others in the great diversity of aromas that make wines famous. Analysis of grapevine reference genome has shown a remarkable expansion of certain families of genes linked to secondary metabolism in comparison with the other plants. In this work, I have analysed gene families coding for cytochromes P450, some of them being involved in the production of aromas, genes coding for stilbene synthases (STS), endo-β-1,3-glucanases and NBS type resistance genes involved in grapevine defences. My thesis intends to propose hypothesis to explain the structural organisation of these families and therefore better understand why some of these families are amplified in the grapevine genome. Bioinformatic approaches have been used to study these different genes families. The cytochromes P450 and R genes of NBS type were manually annotated to improve the knowledge of these families of genes. The expression of endo-β-1,3-glucanases, STS and cytochromes P450 genes has been quantified using a large-scale transcriptomic approach. To this purpose, a tool has been developed during this thesis to estimate the level of genes expression from RNA- Seq data available in public databases. In the meantime, DNA resequencing data from 56 cultivars and grapevine species have been analysed to identify structural variations of CNV types within the genes with a NBS domain and the STS genes. These works showed that the amplification of the gene families of interest was not specific to the reference genome but occurred at the scale of the Vitis genus, but also to highlighted structural variations in different genomes. Regarding the STS genes, blocks of duplication and more conserved and expressed genes were identified. For the genes with NBS domain, a clustered organisation has been highlighted with some clusters varying more than others in the studied genotypes. These works contribute to a better knowledge of gene families for efficient and durable defence against pathogens and optimal aromas synthesis in grapevine. This knowledge will benefit to breeding programs currently in progress at INRA Colmar.

Mouse embryonic stem (ES) cells have played a crucial role in studying developmental processes and gene function in vivo. They are extremely useful in the generation of transgenic animals as they can be genetically manipulated and subsequently microinjected into blastocyst stage embryos, where they combine with the inner cell mass and contribute to the developing embryo. Some of the resulting pups are chimaeric, consisting of a mixture of cells derived from the host blastocyst and the injected ES cells. We have identified several ES cell clones arising from gene targeting experiments with an impaired capacity to generate viable chimaeras. When injected into blastocysts, these clones cause embryonic death during mid to late gestation, suggesting that the cells are able to contribute to the embryo but interfere with normal embryonic development. The aim of this work was to identify the underlying changes in the transcriptome, epigenome or cell surface markers that have occurred in these compromised ES cells and to further define the developmental phenotype of the chimaeric embryos. Different stages during development were analysed and whereas there was little difference in embryonic death at gestational day e13.5, there was a significant decrease in embryos surviving to gestational day e17.5. Additionally, severe haemorrhaging was observed in all the dead embryos and small foci of haemorrhaging could also be seen in a number of embryos that were still alive. This was also observed at e13.5, albeit to a less severe extent. Using RNA sequencing to discover differences in the transcriptome between control ES cells and the compromised ES cells, five genes were identified that were downregulated in the compromised cells. Four of these, Gtl2, Rtl1as, Rian and Mirg are all located in the imprinted Dlk1-Dio3 region on chromosome 12 and are normally expressed from the maternal genome. This pattern was also validated in tissues from e17.5 chimaeric embryos. The expression of this locus is to a large extent regulated by a differentially methylated region located approximately 13kb upstream of the Gtl2 promoter, the IG-DMR. Whereas this is usually only methylated on the paternal copy, in the compromised ES cells both the paternal and the maternal copy were fully methylated, likely causing the silencing of Gtl2, Rtl1as, Rian and Mirg. Using the DNA methyltransferase inhibitor 5-azacytidine, expression of Gtl2 could be rescued. Injection of those 5-azacytidine treated cells into blastocysts did partially rescue the embryonic lethal phenotype. Additionally, cell surface markers were analysed in a phenotypic screen using phage display. NGS analysis of the phage outputs indicates that there may be additional differences in cell surface markers between the control and compromised ES cell clones, but their specific details remain to be identified. Overall, we have identified the maternally expressed genes of the Dlk1-Dio3 region as markers that can distinguish between ES cells with normal or compromised developmental potency and propose to include these genes in the pre-blastocyst injection screening routine for experiments involving the production of chimaeras or genetically modified mouse strains.

RNA sequencing (RNA-seq) technologies facilitate the characterisation of genes and transcripts in different cell types as well as their expression analysis across various conditions. Due to its ability to provide in-depth insights into transcription and post-transcription mechanisms, RNA-seq has been extensively used in functional genetics and transcriptomics, system biology and developmental biology in animals, plants, diseases, etc. The aim of this project is to use mathematical and computational models to integrate big genomic and transcriptomic data from high-throughput technologies in plant biology and develop new methods to identify which genes or transcripts have significant expression variation across experimental conditions of interest, then to interpret the regulatory causalities of these expression changes by distinguishing the effects from the transcription and alternative splicing. We performed a high resolution ultra-deep RNA-seq time-course experiment to study Arabidopsis in response to cold treatment where plants were grown at 20^oC and then the temperature was reduced to 4^oC. We have developed a high quality <i>Arabidopsis thaliana</i> Reference Transcript Dataset (AtRTD2) transcriptome for accurate transcript and gene quantification. This high quality time-series dataset was used as the benchmark for novel method development and downstream expression analysis. The main outcomes of this project include three parts. i) A pipeline for differential expression (DE) and differential alternative splicing (DAS) analysis at both gene and transcript levels. Firstly, we implemented data pre-processing to reduce the noise/low expression, batch effects and technical biases of read counts. Then we used the limma-voom pipeline to compare the expression at corresponding time-points of 4^oC to the time-points of 20^oC. We identified 8,949 genes with altered expression of which 2,442 showed significant DAS and 1,647 were only regulated by AS. Compared with current publications, 3,039 of these genes were novel cold-responsive genes. In addition, we identified 4,008 differential transcript usage (DTU) transcripts of which the expression changes were significantly different to their cognate DAS genes. ii) A TSIS R package for time-series transcript isoform switch (IS) analysis was developed. IS refers to the time-points when a pair of transcript isoforms from the same gene reverse their relative expression abundances. By using a five metric scheme to evaluate robustly the qualities of each switch point, we identified 892 significant ISs between the high abundance transcripts in the DAS genes and about 57% of these switches occurred very rapidly between 0-6h following transfer to 4^oC. iii) A RLowPC R package for co-expression network construction was generated. The RLowPC method uses a two-step approach to select the high-confidence edges first by reducing the search space by only picking the top ranked genes from an initial partial correlation analysis, and then computes the partial correlations in the confined search space by only removing the linear dependencies from the shared neighbours, largely ignoring the genes showing lower association. In future work, we will construct dynamic transcriptional and AS regulatory networks to interpret the causalities of DE and DAS. We will study the coupling and de-coupling of expression rhythmicity to the Arabidopsis circadian clock in response to cold. We will develop new methods to improve the statistical power of expression comparative analysis, such as by taking into account the missing values of expression and by distinguishing the technical and biological variabilities.

<p>Fungi display diverse growth forms. Some grow as unicellular yeasts, some grow as multicellular hyphae, while others switch between these two growth forms, i.e., the dimorphic fungi. Dimorphism is found in many pathogenic fungi, and it is thought to be a strategy to maximize their fitness during different stages of life cycles. The corn smut fungus <i>Ustilago maydis</i> serves as a renowned model organism for studying fungal dimorphism and its role in pathogenesis. However, knowledge only from the model species may not be expanded to other species unless multispecies studies have been demonstrated. In this dissertation, I performed comparative analyses to examine if knowledge from <i>U. maydis</i> is translational to other dimorphic fungi. First, a physiological study was conducted to find what can serve as a common signal for dimorphic transition of several Ustilaginomycotina species. I found that the lipid serves as a potential common cue for yeast-to-hyphal transition in most dimorphic species, while alternate types of energy-source carbohydrate do not affect fungal dimorphism. In addition, pectin and high temperature can also trigger filamentous growth in some Ustilaginomycotina species. Second, I performed comparative transcriptomics to determine if a mechanism for yeast-to-hyphal dimorphic transition is conserved across multiple dimorphic species. Three species of Ustilaginomycotina (<i>U. maydis</i>, <i>Tilletiopsis washingtonensis </i>and <i>Meira miltonrushii</i>) plus one species from Ascomycota (<i>Ophiostoma novo-ulmi</i>) were included in the analyses. I found that the similarity of transcriptomic alteration is not dependent on phylogenetic relatedness. Genes in amino acid transport and metabolism, energy production and conversion and cytoskeleton are commonly altered during the dimorphic transition of all studied species. Moreover, I discovered several core genes which can play a conserved role in transducing signals for the dimorphic transition. Finally, I performed comparative analyses of 190 fungal genomes to determine genomic properties that are associated with types of fungal growth form. I found that small genome size is a characteristic for yeast-like fungi. Few indicator genes, such as genes encoding proteins in the NADPH oxidase complex and cytoskeletons, which are predominantly lost in yeast-like fungi in both Ascomycota and Basidiomycota. However, many other genes are associated with types of growth form in a lineage-specific manner. Findings from this dissertation will serve as fundamentals for future research in fungal cell biology, especially in fungal dimorphism. Additionally, results from this study suggest cautions when extrapolating results from model species onto non-model species.</p>