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.

Quantification of Toxin Biosynthesis Genes In Cyanobacteria and Dinoflagellates - Genetic Factors as Predictors of Toxin Production in the Environment

Harmful algal blooms (HABs) are events caused by the massive proliferation of microscopic, often photosynthetic organisms that inhabit both fresh and marine waters. Although HABs are essentially a natural phenomenon, they now cause worldwide concern. Recent anthropogenic effects, such as climate change and eutrophication via nutrient runoff, can be seen in their increased prevalence and severity. Cyanobacteria and dinoflagellates are often the causative organisms of HABs. In addition to adverse effects caused by the sheer biomass, certain species produce highly potent toxic compounds: hepatotoxic microcystins are produced exclusively by cyanobacteria and neurotoxic saxitoxins, also known as paralytic shellfish toxins (PSTs), by both cyanobacteria and dinoflagellates. Specific biosynthetic genes in the cyanobacterial genomes direct the production of microcystin and paralytic shellfish toxins. Recently also the first paralytic shellfish toxin gene sequences from dinoflagellate genomes have been elucidated. The public health risks presented by HABs are evident, but the monitoring and prediction of toxic events is challenging. Characterization of the genetic background of toxin biosynthesis, including that of microcystins and paralytic shellfish toxins, has made it possible to develop highly sensitive molecular tools which have shown promise in the monitoring and study of potentially toxic microalgae. In this doctoral work, toxin-specific genes were targeted in the developed PCR and qPCR assays for the detection and quantification of potentially toxic cyanobacteria and dinoflagellates in the environment. The correlation between the copy numbers of the toxin biosynthesis genes and toxin production were investigated to assess whether the developed methods could be used to predict toxin concentrations. The nature of the correlation between gene copy numbers and amount of toxin produced varied depending on the targeted gene and the producing organism. The combined mcyB copy numbers of three potentially microcystin-producing cyanobacterial genera showed significant positive correlation to the observed total toxin production. However, the presence of PST-specific sxtA, sxtG, and sxtB genes of cyanobacterial origin was found to be a poor predictor of toxin production in the studied area. Conversely, the dinoflagellate sxtA4 was a good qualitative indicator of a neurotoxic bloom both in the laboratory and in the field, and population densities reflected well the observed toxin concentrations. In conclusion, although the specificity of each potential targeted toxin biosynthesis gene must be assessed individually during method development, the results obtained in this doctoral study support the use of quantitative PCR -based approaches in the monitoring of toxic cyanobacteria and dinoflagellates.

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.

Alternative electron transfer routes involved in photoprotection of cyanobacteria

In oxygenic photosynthesis, the highly oxidizing reactions of water splitting produce reactive oxygen species (ROS) and other radicals that could damage the photosynthetic apparatus and affect cell viability. Under particular environmental conditions, more electrons are produced in water oxidation than can be harmlessly used by photochemical processes for the reduction of metabolic electron sinks. In these circumstances, the excess of electrons can be delivered, for instance, to O2, resulting in the production of ROS. To prevent detrimental reactions, a diversified assortment of photoprotection mechanisms has evolved in oxygenic photosynthetic organisms. In this thesis, I focus on the role of alternative electron transfer routes in photoprotection of the cyanobacterium Synechocystis sp. PCC 6803. Firstly, I discovered a novel subunit of the NDH-1 complex, NdhS, which is necessary for cyclic electron transfer around Photosystem I, and provides tolerance to high light intensities. Cyclic electron transfer is important in modulating the ATP/NADPH ratio under stressful environmental conditions. The NdhS subunit is conserved in many oxygenic phototrophs, such as cyanobacteria and higher plants. NdhS has been shown to link linear electron transfer to cyclic electron transfer by forming a bridge for electrons accumulating in the Ferredoxin pool to reach the NDH-1 complexes. Secondly, I thoroughly investigated the role of the entire flv4-2 operon in the photoprotection of Photosystem II under air level CO2 conditions and varying light intensities. The operon encodes three proteins: two flavodiiron proteins Flv2 and Flv4 and a small Sll0218 protein. Flv2 and Flv4 are involved in a novel electron transport pathway diverting electrons from the QB pocket of Photosystem II to electron acceptors, which still remain unknown. In my work, it is shown that the flv4-2 operon-encoded proteins safeguard Photosystem II activity by sequestering electrons and maintaining the oxidized state of the PQ pool. Further, Flv2/Flv4 was shown to boost Photosystem II activity by accelerating forward electron flow, triggered by an increased redox potential of QB. The Sll0218 protein was shown to be differentially regulated as compared to Flv2 and Flv4. Sll0218 appeared to be essential for Photosystem II accumulation and was assigned a stabilizing role for Photosystem II assembly/repair. It was also shown to be responsible for optimized light-harvesting. Thus, Sll0218 and Flv2/Flv4 cooperate to protect and enhance Photosystem II activity. Sll0218 ensures an increased number of active Photosystem II centers that efficiently capture light energy from antennae, whilst the Flv2/Flv4 heterodimer provides a higher electron sink availability, in turn, promoting a safer and enhanced activity of Photosystem II. This intertwined function was shown to result in lowered singlet oxygen production. The flv4-2 operon-encoded photoprotective mechanism disperses excess excitation pressure in a complimentary manner with the Orange Carotenoid Protein-mediated non-photochemical quenching. Bioinformatics analyses provided evidence for the loss of the flv4-2 operon in the genomes of cyanobacteria that have developed a stress inducible D1 form. However, the occurrence of various mechanisms, which dissipate excitation pressure at the acceptor side of Photosystem II was revealed in evolutionarily distant clades of organisms, i.e. cyanobacteria, algae and plants.