Symbiont-specific responses to environmental cues in a threesome lichen symbiosis

Photosymbiodemes are a special case of lichen symbiosis where one lichenized fungus engages in symbiosis with two different photosynthetic partners, a cyanobacterium and a green alga, to develop two distinctly looking photomorphs. We investigated differential gene expression in photosymbiodemes of the lichen Peltigera britannica at different temperatures representing mild and putatively stressful conditions and compared gene expression of thallus sectors containing cyanobacterial photobionts with thallus sectors with both green algal and cyanobacterial photobionts. Firstly, because of known ecological differences between photomorphs, we investigated symbiont-specific responses in gene expression to temperature increases. Secondly, we quantified photobiont-mediated differences in fungal gene expression. High temperatures expectedly led to an upregulation of genes involved in heat shock responses in all organisms in whole transcriptome data. As expected, the expression of genes involved i..., RNA was isolated with the innuPREP Plant RNA Kit from Analytik Jena. RNA libraries were constructed with dual indexing using the TruSeq® Stranded mRNA Library Prep kit (Illumina) which included a poly-A selection step. Libraries were sequenced on a MiSeq as well as a HiSeq 3000/4000 SR platform (both Illumina). The processed paired-end MiSeq data was used for de novo transcriptome assembly with Trinity software version 2.4.0 (Haas et al., 2013). The HiSeq data was pseudoaligned to the de novo transcriptome with the RNA-seq quantification program kallisto version 0.45.0 (Bray et al., 2016). Differential gene expression analysis was conducted with the program DESeq2 version 1.22.2 (Love et al, 2014). All genes were taxonomically assigned with MEGAN6 version 6.13.1 (Huson et al., 2007); only ascomycete, chlorophytes and cyanobacterial genes were retained for differential gene expression analysis. Functional annotation of differentially expressed genes was conducted using UniProt BLAST (The..., The 'TrinityAscomycetePhotomorphGenes.fasta', 'TrinityAscomyceteTemperatureGenes.fasta','TrinityChlorophytesGenes.fasta' and 'TrinityCyanobacteriaGenes.fasta' files contain the transcripts obtained from the de novo transcriptome assembly, comprising only the longest isoforms per 'gene'. The 'TrinityChlorophytesGenes.fasta' file contains only the first 200 chlorophytes transcripts that were identified as significantly differentially expressed with DESeq2. The same is true for the 'TrinityCyanobacteriaGenes.fasta' file which contains the 200 most significantly differentially expressed cyanobacteria genes. Differential gene expression analysis of the chlorophytes and cyanobacteria genes was based on temperature increase (25 °C vs control (4 °C)). The 'TrinityAscomycetePhotomorphGenes.fasta' file contains the top200 photomorph-mediated significantly differentially expressed ascomycete transcripts and the 'TrinityAscomyceteTemperatureGenes.fasta' file contains the top200 temperature-mediated...