Characterisation of Cladonia rangiferina transcriptome and genome, and the effects of dehydration and rehydration on its gene expression

Lichens are symbiotic organisms, which consist of the fungal partner and the photosynthetic partner, which can be either an alga or a cyanobacterium. In some lichen species the symbiosis is tripartite, where the relationship includes both an alga and a cyanobacterium alongside the primary symbiont, fungus. The lichen symbiosis is an evolutionarily old adaptation to life on land and many extant fungal species have evolved from lichenised ancestors. Lichens inhabit a wide range of habitats and are capable of living in harsh environments and on nutrient poor substrates, such as bare rocks, often enduring frequent cycles of drying and wetting. Most lichen species are desiccation tolerant, and they can survive long periods of dehydration, but can rapidly resume photosynthesis upon rehydration. The molecular mechanisms behind lichen desiccation tolerance are still largely uncharacterised and little information is available for any lichen species at the genomic or transcriptomic level. The emergence of the high-throughput next generation sequencing (NGS) technologies and the subsequent decrease in the cost of sequencing new genomes and transcriptomes has enabled non-model organism research on the whole genome level. In this doctoral work the transcriptome and genome of the grey reindeer lichen, Cladonia rangiferina, were sequenced, de novo assembled and characterised using NGS and traditional expressed sequence tag (EST) technologies. RNA extraction methods were optimised to improve the yield and quality of RNA extracted from lichen tissue. The effects of rehydration and desiccation on C. rangiferina gene expression on whole transcriptome level were studied and the most differentially expressed genes were identified. The secondary metabolites present in C. rangiferina decreased the quality – integrity, optical characteristics and utility for sensitive molecular biological applications – of the extracted RNA requiring an optimised RNA extraction method for isolating sufficient quantities of high-quality RNA from lichen tissue in a time- and cost-efficient manner. The de novo assembly of the transcriptome of C. rangiferina was used to produce a set of contiguous unigene sequences that were used to investigate the biological functions and pathways active in a hydrated lichen thallus. The de novo assembly of the genome yielded an assembly containing mostly genes derived from the fungal partner. The assembly was of sufficient quality, in size similar to other lichen-forming fungal genomes and included most of the core eukaryotic genes. Differences in gene expression were detected in all studied stages of desiccation and rehydration, but the largest changes occurred during the early stages of rehydration. The most differentially expressed genes did not have any annotations, making them potentially lichen-specific genes, but several genes known to participate in environmental stress tolerance in other organisms were also identified as differentially expressed.

De novo transcriptome assembly and annotation of the isopod Idotea balthica

The Baltic Sea is a unique environment that contains unique genetic populations. In order to study these populations on a genetic level basic molecular research is needed. The aim of this thesis was to provide a basic genetic resource for population genomic studies by de novo assembling a transcriptome for the Baltic Sea isopod Idotea balthica. RNA was extracted from a whole single adult male isopod and sequenced using Illumina (125bp PE) RNA-Seq. The reads were preprocessed using FASTQC for quality control, TRIMMOMATIC for trimming, and RCORRECTOR for error correction. The preprocessed reads were then assembled with TRINITY, a de Bruijn graph-based assembler, using different k-mer sizes. The different assemblies were combined and clustered using CD-HIT. The assemblies were evaluated using TRANSRATE for quality and filtering, BUSCO for completeness, and TRANSDECODER for annotation potential. The 25-mer assembly was annotated using PANNZER (protein annotation with z-score) and BLASTX. The 25-mer assembly represents the best first draft assembly since it contains the most information. However, this assembly shows high levels of polymorphism, which currently cannot be differentiated as paralogs or allelic variants. Furthermore, this assembly is incomplete, which could be improved by sampling additional developmental stages.

Annotation of Gyrodactylus salaris transcriptome

Relationship between organisms within an ecosystem is one of the main focuses in the study of ecology and evolution. For instance, host-parasite interactions have long been under close interest of ecology, evolutionary biology and conservation science, due to great variety of strategies and interaction outcomes. The monogenean ecto-parasites consist of a significant portion of flatworms. Gyrodactylus salaris is a monogenean freshwater ecto-parasite of Atlantic salmon (Salmo salar) whose damage can make fish to be prone to further bacterial and fungal infections. G. salaris is the only one parasite whose genome has been studied so far. The RNA-seq data analyzed in this thesis has already been annotated by using LAST. The RNA-seq data was obtained from Illumina sequencing i.e. yielded reads were assembled into 15777 transcripts. Last resulted in annotation of 46% transcripts and remaining were left unknown. This thesis work was started with whole data and annotation process was continued by the use of PANNZER, CDD and InterProScan. This annotation resulted in 56% successfully annotated sequences having parasite specific proteins identified. This thesis represents the first of Monogenean transcriptomic information which gives an important source for further research on this specie. Additionally, comparison of annotation methods interestingly revealed that description and domain based methods perform better than simple similarity search methods. Therefore it is more likely to suggest the use of these tools and databases for functional annotation. These results also emphasize the need for use of multiple methods and databases. It also highlights the need of more genomic information related to G. salaris.

Interleukin-4 induced leukocyte differentiation

Monocytes, macrophages and dendritic cells (DCs) are important mediators of innate immune system, whereas T lymphocytes are the effector cells of adaptive immune responses. DCs play a crucial role in bridging innate and adaptive immunity. Naïve CD4+ Th progenitors (Thp) differentiate to functionally distinct effector T cell subsets including Th1, Th2 and Th17 cells, which while being responsible for specific immune functions have also been implicated in pathological responses, such as autoimmunity, asthma and allergy. The main objective of this thesis is to dissect the signalling networks involved in the IL-4 induced differentiation of two important leukocyte subtypes, Th2 cells and DCs. Gene expression profiling lead to identification of over 200 genes which are differentially expressed during cytokine induced differentiation of human monocytes to DCs or macrophages and which are likely to be essential for the proper biological functions of these cell types. Transcriptome analysis demonstrated the dynamic regulation of gene expression by IL-12 and IL-4 during the initiation of Th cell differentiation, which was partly counteracted by an immunosuppressive cytokine, TGFβ, present in the culture media. Results from RNAi mediated gene knockdown experiments and global gene expression analysis elucidated that SATB1 regulates multiple genes important for Th cell polarization or function as well as may compete with GATA3 for the reciprocal regulation of IL-5 transcription. In conclusion, the results obtained have extended our system-level understanding of the immune cell differentiation processes and provide an excellent basis for the further functional studies which could lead to development of improved therapeutic approaches for a range of immunological conditions.

Transcriptional Profiling of Organ‐Specific Autoimmunity in Human

Our understanding of the pathogenesis of organ‐specific autoinflammation has been restricted by limited access to the target organs. Peripheral blood, however, as a preferred transportation route for immune cells, provides a window to assess the entire immune system throughout the body. Transcriptional profiling with RNA stabilizing blood collection tubes reflects in vivo expression profiles at the time the blood is drawn, allowing detection of the disease activity in different samples or within the same sample over time. The main objective of this Ph.D. study was to apply gene‐expression microarrays in the characterization of peripheral blood transcriptional profiles in patients with autoimmune diseases. To achieve this goal a custom cDNA microarray targeted for gene‐expression profiling of human immune system was designed and produced. Sample collection and preparation was then optimized to allow gene‐expression profiling from whole‐blood samples. To overcome challenges resulting from minute amounts of sample material, RNA amplification was successfully applied to study pregnancy related immunosuppression in patients with multiple sclerosis (MS). Furthermore, similar sample preparation was applied to characterize longitudinal genome‐wide expression profiles in children with type 1 diabetes (T1D) associated autoantibodies and eventually clinical T1D. Blood transcriptome analyses, using both the ImmunoChip cDNA microarray with targeted probe selection and genome‐wide Affymetrix U133 Plus 2.0 oligonucleotide array, enabled monitoring of autoimmune activity. Novel disease related genes and general autoimmune signatures were identified. Notably, down‐regulation of the HLA class Ib molecules in peripheral blood was associated with disease activity in both MS and T1D. Taken together, these studies demonstrate the potential of peripheral blood transcriptional profiling in biomedical research and diagnostics. Imbalances in peripheral blood transcriptional activity may reveal dynamic changes that are relevant for the disease but might be completely missed in conventional cross‐sectional studies.

Functional Genomic Approaches for Deciphering Plant-Virus Interactions

Plant-virus interactions are very complex in nature and lead to disease and symptom formation by causing various physiological, metabolic and developmental changes in the host plants. These interactions are mainly the outcomes of viral hijacking of host components to complete their infection cycles and of host defensive responses to restrict the viral infections. Viral genomes contain only a small number of genes often encoding for multifunctional proteins, and all are essential in establishing a viral infection. Thus, it is important to understand the specific roles of individual viral genes and their contribution to the viral life cycles. Among the most important viral proteins are the suppressors of RNA silencing (VSRs). These proteins function to suppress host defenses mediated by RNA silencing and can also serve in other functions, e.g. in viral movement, transactivation of host genes, virus replication and protein processing. Thus these proteins are likely to have a significant impact on host physiology and metabolism. In the present study, I have examined the plant-virus interactions and the effects of three different VSRs on host physiology and gene expression levels by microarray analysis of transgenic plants that express these VSR genes. I also studied the gene expression changes related to the expression of the whole genome of Tobacco mosaic virus (TMV) in transgenic tobacco plants. Expression of the VSR genes in the transgenic tobacco plants causes significant changes in the gene expression profiles. HC-Pro gene derived from the Potyvirus Y (PVY) causes alteration of 748 and 332 transcripts, AC2 gene derived from the African cassava mosaic virus (ACMV) causes alteration of 1118 and 251transcripts, and P25 gene derived from the Potyvirus X (PVX) causes alterations of 1355 and 64 transcripts in leaves and flowers, respectively. All three VSRs cause similar up-regulation in defense, hormonally regulated and different stress-related genes and down-regulation in the photosynthesis and starch metabolism related genes. They also induce alterations that are specific to each viral VSR. The phenotype and transcriptome alterations of the HC-Pro expressing transgenic plants are similar to those observed in some Potyvirus-infected plants. The plants show increased protein degradation, which may be due to the HC-Pro cysteine endopeptidase and thioredoxin activities. The AC2-expressing transgenic plants show a similar phenotype and gene expression pattern as HC-Pro-expressing plants, but also alter pathways related to jasmonic acid, ethylene and retrograde signaling. In the P25 expressing transgenic plants, high numbers of genes (total of 1355) were up-regulated in the leaves, compared to a very low number of down-regulated genes (total of 5). Despite of strong induction of the transcripts, only mild growth reduction and no other distinct phenotype was observed in these plants. As an example of whole virus interactions with its host, I also studied gene expression changes caused by Tobacco mosaic virus (TMV) in tobacco host in three different conditions, i.e. in transgenic plants that are first resistant to the virus, and then become susceptible to it and in wild type plants naturally infected with this virus. The microarray analysis revealed up and down-regulation of 1362 and 1422 transcripts in the TMV resistant young transgenic plants, and up and down-regulation of a total of 1150 and 1200 transcripts, respectively, in the older plants, after the resistance break. Natural TMV infections in wild type plants caused up-regulation of 550 transcripts and down-regulation of 480 transcripts. 124 up-regulated and 29 down-regulated transcripts were commonly altered between young and old TMV transgenic plants, and only 6 up-regulated and none of the down-regulated transcripts were commonly altered in all three plants. During the resistant stage, the strong down-regulation in translation-related transcripts (total of 750 genes) was observed. Additionally, transcripts related to the hormones, protein degradation and defense pathways, cell division and stress were distinctly altered. All these alterations may contribute to the TMV resistance in the young transgenic plants, and the resistance may also be related to RNA silencing, despite of the low viral abundance and lack of viral siRNAs or TMV methylation activity in the plants.

Chromatoid body mediated RNA regulation in mouse male germline

Male germ cell differentiation, spermatogenesis is an exceptional developmental process that produces a massive amount of genetically unique spermatozoa. The complexity of this process along with the technical limitations in the germline research has left many aspects of spermatogenesis poorly understood. Post-meiotic haploid round spermatids possess the most complex transcriptomes of the whole body. Correspondingly, efficient and accurate control mechanisms are necessary to deal with the huge diversity of transcribed RNAs in these cells. The high transcriptional activity in round spermatids is accompanied by the presence of an uncommonly large cytoplasmic ribonucleoprotein granule, called the chromatoid body (CB) that is conjectured to participate in the RNA post-transcriptional regulation. However, very little is known about the possible mechanisms of the CB function. The development of a procedure to isolate CBs from mouse testes was this study’s objective. Anti-MVH immunoprecipitation of cross-linked CBs from a fractionated testicular cell lysate was optimized to yield considerable quantities of pure and intact CBs from mice testes. This protocol produced reliable and reproducible data from the subsequent analysis of CB’s protein and RNA components. We found that the majority of the CB’s proteome consists of RNA-binding proteins that associate functionally with different pathways. We also demonstrated notable localization patterns of one of the CB transient components, SAM68 and showed that its ablation does not change the general composition or structure of the CB. CB-associated RNA analysis revealed a strong accumulation of PIWI-interacting RNAs (piRNAs), mRNAs and long non-coding RNAs (lncRNAs) in the CB. When the CB transcriptome and proteome analysis results were combined, the most pronounced molecular functions in the CB were related to piRNA pathway, RNA post-transcriptional processing and CB structural scaffolding. In addition, we demonstrated that the CB is a target for the main RNA flux from the nucleus throughout all steps of round spermatid development. Moreover, we provided preliminary evidence that those isolated CBs slice target RNAs in vitro in an ATPdependent manner. Altogether, these results make a strong suggestion that the CB functions involve RNA-related and RNA-mediated mechanisms. All the existing data supports the hypothesis that the CB coordinates the highly complex haploid transcriptome during the preparation of the male gametes for fertilization. Thereby, this study provides a fundamental basis for the future functional analyses of ribonucleoprotein granules and offers also important insights into the mechanisms governing male fertility.

Transcriptional and Epigenetic Regulation of Human CD4+ T Helper Lineage Specification

Activated T helper (Th) cells have ability to differentiate into functionally distinct Th1, Th2 and Th17 subsets through a series of overlapping networks that include signaling and transcriptional control and the epigenetic mechanisms to direct immune responses. However, inappropriate execution in the differentiation process and abnormal function of these Th cells can lead to the development of several immune mediated diseases. Therefore, the thesis aimed at identifying genes and gene regulatory mechanisms responsible for Th17 differentiation and to study epigenetic changes associated with early stage of Th1/Th2 cell differentiation. Genome wide transcriptional profiling during early stages of human Th17 cell differentiation demonstrated differential regulation of several novel and currently known genes associated with Th17 differentiation. Selected candidate genes were further validated at protein level and their specificity for Th17 as compared to other T helper subsets was analyzed. Moreover, combination of RNA interference-mediated downregulation of gene expression, genome-wide transcriptome profiling and chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq), combined with computational data integration lead to the identification of direct and indirect target genes of STAT3, which is a pivotal upstream transcription factor for Th17 cell polarization. Results indicated that STAT3 directly regulates the expression of several genes that are known to play a role in activation, differentiation, proliferation, and survival of Th17 cells. These results provide a basis for constructing a network regulating gene expression during early human Th17 differentiation. Th1 and Th2 lineage specific enhancers were identified from genome-wide maps of histone modifications generated from the cells differentiating towards Th1 and Th2 lineages at 72h. Further analysis of lineage-specific enhancers revealed known and novel transcription factors that potentially control lineage-specific gene expression. Finally, we found an overlap of a subset of enhancers with SNPs associated with autoimmune diseases through GWASs suggesting a potential role for enhancer elements in the disease development. In conclusion, the results obtained have extended our knowledge of Th differentiation and provided new mechanistic insights into dysregulation of Th cell differentiation in human immune mediated diseases.

Gene Regulation in The Human Immune System. Gene Expression Signatures of Th2 Cell Differentiation, Type 1 Diabetes, and Intrauterine Immune Adaptation

The human immune system is constantly interacting with the surrounding stimuli and microorganisms. However, when directed against self or harmless antigens, these vital defense mechanisms can cause great damage. In addition, the understanding the underlying mechanism of several human diseases caused by aberrant immune cell functions, for instance type 1 diabetes and allergies, remains far from being complete. In this Ph.D. study these questions were addressed using genome-wide transcriptomic analyses. Asthma and allergies are characterized by a hyperactive response of the T helper 2 (Th2) immune cells. In this study, the target genes of the STAT6 transcription factor in naïve human T cells were identified with RNAi for the first time. STAT6 was shown to act as a central activator of the genes expression upon IL-4 signaling, with both direct and indirect effects on Th2 cell transcriptome. The core transcription factor network induced by IL-4 was identified from a kinetic analysis of the transcriptome. Type 1 diabetes is an autoimmune disease influenced by both the genetic susceptibility of an individual and the disease-triggering environmental factors. To improve understanding of the autoimmune processes driving pathogenesis in the prediabetic phase in humans, a unique series of prospective whole-blood RNA samples collected from HLA-susceptible children in the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) study was studied. Changes in different timewindows of the pathogenesis process were identified, and especially the type 1 interferon response was activated early and throughout the preclinical T1D. The hygiene hypothesis states that allergic diseases, and lately also autoimmune diseases, could be prevented by infections and other microbial contacts acquired in early childhood, or even prenatally. To study the effects of the standard of hygiene on the development of neonatal immune system, cord blood samples from children born in Finland (high standard of living), Estonia (rapid economic growth) and Russian Karelia (low standard of living) were compared. Children born in Russian Karelia deviated from Finnish and Estonian children in many aspects of the neonatal immune system, which was developmentally more mature in Karelia, resembling that of older infants. The results of this thesis offer significant new information on the regulatory networks associated with immune-mediated diseases in human. The results will facilitate understanding and further research on the role of the identified target genes and mechanisms driving the allergic inflammation and type 1 diabetes, hopefully leading to a new era of drug development.