Refine
Year of publication
- 2018 (4) (remove)
Document Type
- Doctoral Thesis (4)
Language
- English (4)
Has Fulltext
- yes (4)
Is part of the Bibliography
- no (4)
Institute
Heat stress transcription factors (Hsfs) have an essential role in heat stress response (HSR) and thermotolerance by controlling the expression of hundreds of genes including heat shock proteins (Hsps) with molecular chaperone functions. Hsf family in plants shows a striking multiplicity, with more than 20 members in many species. In Solanum lycopersicum HsfA1a was reported to act as the master regulator of the onset of HSR and therefore is essential for basal thermotolerance. Evidence for this was provided by the analysis of HsfA1a co-suppression (A1CS) transgenic plants, which exhibited hypersensitivity upon exposure to heat stress (HS) due to the inability of the plants to induce the expression of many HS-genes including HsfA2, HsfB1 and several Hsps. Completion of tomato genome sequencing allowed the completion of the Hsf inventory, which is consisted of 27 members, including another three HsfA1 genes, namely HsfA1b, HsfA1c and HsfA1e.
Consequently, the suppression effect of the short interference RNA in A1CS lin e was re-evaluated for all HsfA1 genes. We found that expression of all HsfA1 proteins was suppressed in A1CS protoplasts. This result suggested that the model of single master regulator needs to be re-examined.
Expression analysis revealed that HsfA1a is constitutively expressed in different tissues and in response to HS, while HsfA1c and HsfA1e are minimally expressed in general, and show an induction during fruit ripening and a weak upregulation in late HSR. Instead HsfA1b shows preferential expression in specific tissues and is strongly and rapidly induced in response to HS. At the protein level HsfA1b and HsfA1e are rapidly degraded while HsfA1a and HsfA1c show a higher stability. In addition, HsfA1a and HsfA1c show a nucleocytosolic distribution, while HsfA1b and HsfA1e a strong nuclear retention.
A major property of a master regulator in HSR is thought to be its ability to cause a strong transactivation of a wide range of genes required for the initial activation of protective mechanisms. GUS reporter assays as well as analysis of transcript levels of several endogenous transcripts in protoplasts transiently expressing HsfA1 proteins revealed that HsfA1a can stimulate the transcription of many genes, while the other Hsfs have weaker activity and only on limited set of target genes. The low activity of HsfA1c and HsfA1e can be attributed to the lower DNA capacity of the two factors as judged by a GUS reporter repressor assay.
HsfA1a has been shown to have synergistic activity with the stress induced HsfA2 and HsfB1. The formation of such complexes is considered as important for stimulation of transcription and long term stress adaptation. All HsfA1 members show synergistic activity with HsfA2, while only HsfA1a act as co-activator of HsfB1 and HsfA7. Interestingly, HsfA1b shows an exceptional synergistic activity with HsfA3, suggesting that different Hsf complexes might regulate different HS-related gene networks. Altogether these results suggest that HsfA1a has unique characteristics within HsfA1 subfamily. This result is interesting considering the very high sequencing similarity among HsfA1s, and particularly among HsfA1a and HsfA1c.
To understand the molecular basis of this discrepancy, a series of domain swapping mutants between HsfA1a and HsfA1c were generated. Oligomerization domain and C-terminal swaps did not affect the basal activity or co-activity of the proteins. Remarkably, an HsfA1a mutant harbouring the N-terminus of HsfA1c shows reduced activity and co-activity, while the reciprocal HsfA1c with the N-terminus of HsfA1a cause a gain of activity and enhanced DNA binding capacity.
Sequence analysis of the DBD of HsfA1 proteins revealed a divergence in the highly conserved C-terminus of the turn of β3-β4 sheet. As the vast majority of HsfA1 proteins, HsfA1a at this position comprises an Arg residue (R107), while HsfA1c a Leu and HsfA1e a Cys. An HsfA1a-R107L mutant has reduced DNA binding capacity and consequently activity. Therefore, the results presented here point to the essential function of this amino acid residue for DNA binding function. Interestingly, the mutation did not affect the activity of the protein on Hsp70-1, suggesting that the functionality of the DBD and consequently the transcription factor on different promoters with variable heat stress element number and architecture is dependent on structural peculiarities of the DBD.
In conclusion, the unique properties including expression pattern, transcriptional activities, stability, DBD-peculiarities are likely responsible for the dominant function of HsfA1a as a master regulator of HSR in tomato. Instead, other HsfA1-members are only participating in HSR or developmental regulations by regulating a specific set of genes. Furthermore, HsfA1b and HsfA1e are likely function as stress primers in specific tissues while HsfA1c as a co-regulator in mild HSR. Thereby, tomato subclass A1 presents another example of function diversity not only within the Hsf family but also within the Hsf-subfamily of closely related members. The diversification based on DBD peculiarities is likely to occur in potato as well. Therefore this might have eliminated the functional redundancy observed in other species such as Arabidopsis thaliana but has probably allowed the more refined regulation of Hsf networks possibly under different stress regimes, tissues and cell types.
Heat stress transcription factors (Hsfs) are required for transcriptional changes during heat stress (HS) thereby playing a crucial role in the heat stress response (HSR). The target genes of Hsfs include heat shock proteins (Hsps), other Hsfs and genes involved in protection of the cell from irreversible damages due to exposure to elevated temperatures. Among 27 Hsfs in Solanum lycopersicum, HsfA1a, HsfA2 and HsfB1 constitute a functional triad which regulates important aspects of the HSR. HsfA1a is constitutively expressed and described as the master regulator of stress response and thermotolerance. Activation of HsfA1a under elevated temperatures leads to the induction of HsfA2 and HsfB1 which further stimulate the transcription of HS-responsive genes by forming highly active complexes with HsfA1a. Despite the well-established role of these three Hsfs in tomato HSR, information about functional relevance of other Hsfs is currently missing.
The heat stress inducible HsfA7 belongs alongside with HsfA2 to a phylogenetically distinct clade. Thereby the two proteins share high homology and a functional redundancy has been assumed. However, HsfA7 function and contribution to stress responses have not been investigated into detail in any plant species.
Tomato HsfA7 protein accumulates already at moderately elevated temperatures (~35°C) while HsfA2 becomes dominant at higher temperatures (>40°C). HsfA7 pre-mRNA undergoes complex and temperature-dependent alternative splicing resulting in several transcripts that encode for three protein isoforms. HsfA7-I contains a functional nuclear export signal (NES) and shows nucleocytoplasmic shuttling while HsfA7-II and HsfA7-III have a truncated NES which leads to the strong nuclear retention of the protein. Differences in the nucleocytoplasmic equilibrium have a major impact on the stability of protein isoforms, as nuclear retention is associated with increased protein turnover. Consequently, HsfA7-I shows a higher stability and can be detected even after 24 hours of stress attenuation, while HsfA7-II is rapidly degraded. The degradation of these factors is mediated by the ubiquitin-proteasome pathway.
HsfA7 can physically interact with HsfA1a and HsfA3 and form co-activator (“superactivator”) complexes with a very high transcriptional activity as shown on different HS-inducible promoters. In order for the complex to be successfully transferred to the nucleus and confer its activity it needs a functional nuclear localization signal (NLS) of HsfA7. In contrast, the activator (AHA) motif of HsfA7 is not essential for its co-activator function. Interestingly, while interaction of HsfA7 with either HsfA3 or HsfA1a stabilizes HsfA7 isoforms, concomitantly this leads to an increased turnover of HsfA1a and HsfA3. In contrast, HsfA2 has a stabilizing effect on the master regulator HsfA1a.
Thus, HsfA7 knockout mutants generated by CRISPR/Cas9 gene editing, show increased HsfA1a levels and a stronger induction of HS-related genes at 35°C compared to wild-type plants and HsfA2 knockout mutants. Consequently, HsfA7 knockout seedlings exhibit increased thermotolerance as shown by the enhanced hypocotyl elongation under a prolonged mild stress treatment at 35°C. In summary, these results highlight the importance of HsfA7 in regulation of cellular responses at elevated temperatures. Under moderately elevated temperatures, the accumulation of HsfA7 and its subsequent interaction with HsfA1a, leads to increased turnover of the latter, thereby ensuring a milder transcriptional activation of temperature-responsive genes like Hsps. In turn, in response to further elevated temperatures, HsfA2 becomes the dominant stress-induced Hsf. HsfA2 forms co-activator complexes with HsfA1a which in contrast to HsfA7, allows the stabilization of the master regulator, leading to the stronger expression of HS-responsive genes required for survival. Thereby, this study uncovers a new regulatory mechanism, where the temperature-dependent competitive interaction of HsfA2 and HsfA7 with HsfA1a control the fate of the master regulator and consequently the activity of temperature-responsive networks.
Microsporidia are a group of parasites that infect a wide range of species, many of which play important roles in agriculture and human disease. At least 14 microsporidian species have been confirmed to cause potentially lifethreatening infectious diseases in both immunocompromised and immunocompetent humans. Approximately 1,400 species of microsporidia have been described. Depending on their host and habitat they are classified into three groups, the aquasporidia, the terresporidia and the marinosporidia.
Microsporidia were originally classified as fungi by Naegeli (1857). However, their lack of typical eukaryotic components – such as mitochondria, Golgi bodies or peroxisomes – suggested to place the microsporidia together with other amitochondriate protists within the Archezoa kingdom. This "microsporidia-early" hypothesis was further supported by molecular phylogenies inferred from individual genes. Despite this evidence, the placement of microsporidia as an early branching eukaryote remained a topic for debate. The phylogeny of microsporidia is prone to suffer from biases in their reconstruction. The high evolutionary rate of microsporidian proteins tends to place these proteins together with other fast evolving lineages, a phenomenon known as long-branch attraction. In 1996, the first molecular phylogenetic studies placed the microsporidia inside the fungi.
Subsequently, several further studies located the microsporidia at different positions inside the fungal clade. Since then, microsporidia have been considered as members of the Ascomycota, Zygomycota, Cryptomycota, or as a sister group to the Ascomycota and Basidiomycota, or even as the sister group of all fungi.
The difficulties in determining the evolutionary origin of microsporidia are not only caused by their lack of several cellular components but also by their reduced genomes and metabolism. Being obligate intracellular parasites, microsporidia successfully reduced their genome sizes, down to the range of bacteria. As the smallest eukaryotic genome described so far, the genome of Encephalitozoon intestinalis is just 2.3 Mbp, about half the size of the one of Escherichia coli. Due to their low number of protein coding genes (less than 4,000), microsporidia are thought to retain only genes essential for their survival and development. Furthermore, several key metabolic pathways are missing in the microsporidia, such as the citric acid cycle, oxidative phosphorylation, or the de novo biosynthesis of nucleotides. As a result they are in an obligatory dependence on many primary metabolites from the hosts. However, the presence of hsp70 protein suggests a more complex genome of the microsporidian ancestor. Consequently, the small microsporidian genomes and the reduced metabolism would be consequences of a secondary loss process that molded the contemporary microsporidia from a functionally more complex ancestral species. However, it remains unclear whether the last common ancestor (LCA) of the microsporidia was already reduced, or whether the genome compaction was lineage-specific and started from a more complex LCA.
We investigated the evolutionary history of the contemporary microsporidia through the reconstruction and analysis of their LCA. As a first step in our analysis, we have developed and implemented a software facilitating an intuitive data analysis of the large presence absence-patterns resulting from the tracing of microsporidian proteins in gene sets of many different species. These so called phylogenetic profiles can now be dynamically visualized and explored with PhyloProfile. The software allows the integration of other additional information layers into the phylogenetic profile, such as the similarity of feature architecture (FAS) between the protein under study and its orthologs. The FAS score can be displayed along the presence-absence pattern, which can help to identify orthologs that have likely diverged in function. PhyloProfile closes the methodological gap that existed between tools to generate large phylogenetic profiles to delineate the evolutionary history and the contemporary distribution of large – and ultimately complete – gene sets, and the more function-oriented analysis of individual protein. In the next step we tackled the problem of how to transfer functional annotation from one protein to another. We have developed HamFAS that integrates a targeted ortholog search based on the HaMStR algorithm with a weighted assessment of feature architecture similarities (FAS) between orthologs. In brief, for a seed protein we identify orthologs in reference species in which proteins have been functionally annotated based on manually curated assignments to KEGG Ortholog (KO) groups. The FAS scores between the orthologs and seed proteins are calculated. Subsequently, we compute pairwise FAS scores for all reference proteins within a KO group. A group's mean FAS score serves then as cutoff that must be exceeded to warrant transfer of its KO identifier to the seed. A benchmark using a manually curated yeast protein set showed that HamFAS yields the best precision (98.5%) when compared with two state-of-the-art annotation tools, KAAS and BlastKOALA. Furthermore, HamFAS achieves a higher sensitivity. On average HamFAS annotates almost 50% more proteins than KAAS or BlastKOALA.
With this extended bioinformatics toolbox at hand, we aimed at reconstructing the evolutionary history of the microsporidia. We generated a robust phylogeny of microsporidia using a phylogenomics approach. As a data basis, we identified a set of microsporidian proteins encoded by 80 core genes with one-to-one orthologs. A maximum likelihood analysis of this data
with 48 fungi and additionally in 13 species from more distantly related such as animals and plants combined in a supermatrix strongly supported the hypothesis that microsporidia form the sister group of the fungi. We confirmed that the data explains this microsporidia-fungi relationship significantly better than any other of the previously proposed phylogenetic hypotheses.
On the basis of this phylogeny, and of the phylogenetic profiles of microsporidian proteins, we then focused on reconstructing the dynamics microsporidian genome evolution. Between 2% of the proteins in the compact microsporidia Encephalitozoon intestinalis and up to 49% of the proteins of Edhazardia aedis are private for individual microsporidian species. A comparison of the sequence characteristics of these proteins to that of proteins with orthologs in other microsporidian species revealed individual differences. Yet, without further evidences it remains unclear whether these private genes are indeed lineage-specific innovations contributing to the adaptation of each microsporidium to its host, or whether these are artifacts introduced in the process of gene annotation. A total of 14,410 microsporidian proteins could then be grouped into 1605 orthologous groups that can be traced back to the last common ancestor of the microsporidia (LCA set). We found that 94% of the microsporidian LCA proteins could be tracked back to the last eukaryotic common ancestor. The high evolutionary age of these proteins, together with the resistance against gene loss in the microsporidia suggests that the corresponding functions are essential for eukaryotic life. Further 3% of the LCA proteins could be dated to the common ancestor microsporidia share with the fungi. Only 3% of the LCA proteins appear as microsporidia specific inventions. These proteins are potentially of importance for the evolutionary of the obligate parasitic lifestyle nowadays shared by all microsporidia.
The functional annotation and metabolic pathway analysis of the microsporidian LCA protein set gave us more insight into the adaptation of the microsporidia to their parasitic lifestyle and the origin of the microsporidian genome reduction. The presence of E1 and E3 components of the pyruvate dehydrogenase complex and the mitochondrial hsp70 protein support an ancestral presence of mitochondria in the ancestral microsporidia. In addition, several ancient proteins that complement gapped metabolic pathways were found in the microsporidian LCA. They suggested a more complex genome and metabolism in the LCA. However, our reconstruction of the metabolic network of the microsporidian LCA still lacks many main pathways. For example, the TCA cycle for effective energy production, and key enzymes that are required for in vivo synthesis of critical metabolites like purines and pyrimidines appear absent. We therefore find that the parasitic lifestyle and the genome reduction already occurred in the microsporidian LCA. This ancestral state was followed by further losses and gains during the evolution of each individual microsporidian lineage.
In times of a growing world population and the associated demand for high crop yield, the understanding and improvement of plant reproduction is of central importance. One key step of plant reproduction is the development of the male gametophyte, which is better known as pollen. In addition, the development of pollen was shown to be very sensitive to abiotic stresses, such as heat, which can cause crop damage and yield loss. To obtain new insights in the development and heat stress response of pollen, a combined transcriptome and proteome analysis was performed for three pollen developmental stages of non- and heat-stressed tomato plants.
The analysis of the transcriptomes of non-stressed pollen developmental stages enabled the determination of mRNAs accumulated in certain developmental stages. The functional analysis of these mRNAs led to the identification of protein families and functional processes that are important at different times of pollen development. A subsequent comparison of the transcriptomes of non- and heat-stressed pollen revealed a core set of 49 mRNAs, which are upregulated in all three developmental stages. The encoded proteins include among other things different heat stress transcription factors and heat shock proteins, which are known key players of the plant heat stress response.
Furthermore, 793 potential miRNAs could be identified in the transcriptome of non- and heat-stressed pollen. Interestingly, 38 out of the 793 miRNAs have already been identified in plants. For more than half of these miRNAs potential target mRNAs were identified and the interactions between miRNAs and mRNAs linked to the development and heat stress response of pollen. In total, 207 developmentally relevant interactions could be determined, out of which 34 have an effect on transcriptional-networks. In addition, 24 of the interactions contribute the heat stress response of pollen, whereby this mainly affects post-meiotic pollen.
An initial correlation of the proteome and transcriptome of the developmental stages revealed that transcriptome analyses are not sufficient to draw exact conclusions about the state of the proteome. A closer look on the relationship of the transcriptome and proteome during pollen development revealed two translational modes that are active during the development of pollen. One mode leads to a direct translation of mRNAs, while the second mode leads a delayed translation at a later point in time. Regarding the delayed translation, it could be shown that this is likely due to a short-term storage of mRNAs in so-called EPPs. The comparison of the proteome and transcriptome response to heat stress revealed that the proteome reacts much stronger and that the reaction is mainly independent from the transcriptome. Finally, the comparison of the proteome of non- and heat-stressed pollen provided first indications for changes in the ribosome composition in response to heat stress, as 57 ribosomal proteins are differentially regulated in at least one developmental stage.