Next Article in Journal
Auxin Controlled by Ethylene Steers Root Development
Next Article in Special Issue
Quantifying BRCA1 and BRCA2 mRNA Isoform Expression Levels in Single Cells
Previous Article in Journal
Rex Retroelements and Teleost Genomes: An Overview
Previous Article in Special Issue
Mutations of Pre-mRNA Splicing Regulatory Elements: Are Predictions Moving Forward to Clinical Diagnostics?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

STE20/PAKA Protein Kinase Gene Releases an Autoinhibitory Domain through Pre-mRNA Alternative Splicing in the Dermatophyte Trichophyton rubrum

by
Eriston V. Gomes
1,2,
Julio C. Bortolossi
1,
Pablo R. Sanches
1,
Niege S. Mendes
1,
Nilce M. Martinez-Rossi
1,* and
Antonio Rossi
1
1
Department of Genetics, Ribeirão Preto Medical School, São Paulo University, Ribeirão Preto, São Paulo 14049-900, Brazil
2
Department of Biofunctional, Center of Higher Education Morgana Potrich Eireli, Morgana Potrich College, Mineiros, Goiás 75830-000, Brazil
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(11), 3654; https://doi.org/10.3390/ijms19113654
Submission received: 21 September 2018 / Revised: 30 October 2018 / Accepted: 16 November 2018 / Published: 20 November 2018
(This article belongs to the Special Issue Pre-mRNA Splicing 2018)

Abstract

:
Signaling pathways are highly diverse in filamentous fungi, allowing the cells to receive and process ambient information. Interaction of components from different pathways results in signaling networks. The mitogen-activated protein kinase (MAPK) pathway is dependent on phosphorylation that is accomplished by kinase proteins. Thus, the STE/PAK protein kinase family plays essential roles in MAPK signal transduction, regulating several cellular functions. The STE/PAK protein displays an autoinhibitory (Cdc42/Rac interactive binding—CRIB) domain on its N-terminal portion, which interacts with the C-terminal catalytic kinase domain. Based on current knowledge, for the STE/PAK kinase to be activated, molecular signals (e.g., interaction with the activated form of Rac1 and Cdc42 proteins) or proteolytic cleavage by caspase 3 is necessary. Both mechanisms release the kinase domain from the CRIB interaction. Here, we hypothesize a novel molecular mechanism for the activation of STE20/PAKA kinase in Trichophyton rubrum based on an alternative pre-mRNA splicing process. Our data suggest that, because of the retention of intron 1 of this gene, it is theoretically possible that the translation of STE20/PAKA kinase will be free of its autoinhibitory CRIB domain. These findings indicate a rapid response system to environmental changes. Furthermore, STE20/PAKA may be a potential T. rubrum virulence factor and an interesting target for new drugs against dermatophytes.

Graphical Abstract

1. Introduction

Regardless of lifestyle, living organisms including filamentous fungi monitor various environmental signals. It is imperative to sense a broad range of stimuli, including the nutritional availability status and the presence of hosts or predators. Thus, when facing different kinds of signals (e.g., light, temperature, pheromones, and chemical compounds), organisms respond adequately to environmental changes [1,2]. Kinase proteins play a fundamental role in the transduction of cellular signals, from the receptors to the active response, which may include transcription and translation of effector proteins [3]. Mitogen-activated protein kinase (MAPK) pathways are directly or indirectly involved in the molecular responses to these stimuli. These processes display an evolutionarily conserved mechanism from yeast to humans, acting as a switch, turning on or off the activity of its substrate proteins [4]. Pathogenic fungi carrying defective mutations in genes that integrate the MAPK pathway show reduced virulence in animal models of infection [5]. The p21 activated kinase (PAK), a member of the Ste20-related kinase family, was initially isolated from rat brain. The PAK protein displays a conserved N-terminal domain, which serves as a binding site for the activated form of Rac1 and Cdc42 proteins (Cdc42/Rac interactive binding—CRIB), acting as an autoinhibitory mechanism for the C-terminal kinase domain [6]. In the MAPK pathway, Ste20p usually phosphorylates Ste11p (a MAPKKK protein) which follows the MAPK cascade and regulates several cellular responses, such as mating, cell wall maintenance, osmotic regulation, filamentous growth, and virulence [7,8]. Under apoptotic stress, there is an interesting transient and constitutive activation of mammalian PAK2 by caspase 3 (CPP32) cleavage, leading to the transcriptional inhibition of several growth-related genes [9,10,11].
Intron retention (IR) in filamentous fungi is the most common alternative pre-mRNA splicing mechanism (AS) and may be modulated in response to nutrient signaling, as described for the MAPK protein coded by mak-2 gene in Neurospora crassa [12]. As a dermatophyte, Trichophyton rubrum needs to probe the surrounding environment, obtain its nutrients, and consequently develop hyphae with adequate structures to penetrate the host tissue. Moreover, several stress factors must be considered, such as temperature, lack of moisture, ultraviolet light, and the host defense system. Therefore, the MAPK signaling pathway is of fundamental importance for the development and successful interaction and tissue invasion by the fungus [13]. Here, we describe the retention of intron 1 during transcription of the pakA gene and propose a molecular mechanism for the activation of the T. rubrum PAKA kinase as an alternative to the mammalian caspase 3 cleavage through AS. We suggest that retention of intron 1 generates, based on the putative fungal STE20/PAK protein activation and RNA-Seq transcriptome analysis, both a premature stop codon in its respective mRNA and a new start codon soon after that, nonetheless preserving most of the original protein amino acid sequence. Subsequently, at least two independent polypeptides are theoretically produced. The first one contains the regulatory CRIB domain, and the second one—displaying the entire kinase catalytic domain—is free of the regulatory domain. This observation suggests that, depending on the external stimulus, the MAPK metabolic pathway may be activated with variable intensities and regulated by a different set of molecules, such as the production of a PAKA kinase independent of the Cdc42/Rac system regulation. Intron 1 retention is also modulated by the presence of the antifungal agent undecanoic acid (UDA), probably for a rapid response to environmental changes or stress conditions, which also add support to our hypothesis. These findings may help to clarify the mechanism of adaptation and host-fungi interactions, indicating specific novel targets for the development of more efficient drugs in the fight against dermatophytoses.

2. Results

2.1. T. rubrum pakA/Ste20 Kinase Gene

T. rubrum pakA/Ste20 (TERG_03042) is a 3239-bp gene. The complete splicing process of pakA transcription generates a mRNA with 2912 bases, which encodes a 970-amino acid protein. According to ScanProsite, a web-based tool for detecting signature matches in protein sequences [14,15], PAKA exhibits the CRIB regulatory domain on its N-terminal portion, between amino acids 302 and 315. Furthermore, both the catalytic site and the Gβ binding (GBB) domains are located at the C-terminal region (Figure 1). PAKA also presents other intriguing areas, such as six regions rich in specific amino acids, which may be necessary for interactions with other proteins, metabolites, small molecules, and ions (Figure S1A). Moreover, PAKA displays several putative post-translational modification sites. There are 39 phosphorylation sites: 16 sites for protein kinase C (PKC), five for cAMP/cGMP-dependent protein kinase (CAMP), and 18 for casein kinase II (CK2) (Figure S1B). There are also 27 sites for N-myristoylation (Myr), six sites for N-glycosylation (N-gly; only in N-terminal region), and one site for amidation (AMD) inside the kinase domain (Figure S1C,D). We also used EMBOSS epestfind (http://emboss.bioinformatics.nl/cgi-bin/emboss/epestfind) (accessed on 19 June 2018) another web-based tool for detecting PEST motifs as potential proteolytic cleavage sites. The results showed two potential PEST motifs: one with 15 amino acids between position 157 and 173 (PEST score: 15.52), and another with 44 amino acids between position 467 and 512 (PEST score: 8.06) (Figure S1E).

2.2. Transcription Retention of Introns of the pakA Kinase Gene of T. rubrum

Using a high-throughput RNA-sequencing (RNA-Seq) approach, we identified IR events in the transcriptional processing of the T. rubrum pakA gene after treatment with sub-inhibitory doses of UDA. These data indicated that 56 reads exhibited retention of intron 1 in the pakA transcriptional products after 12 h of UDA treatment (RNA-Seq data). To validate these data, a set of primers flanking the intron 1 region of this gene was designed and used for reverse transcription polymerase chain reaction (qRT-PCR) (Figure 2A), followed by 1.5% agarose gel electrophoresis. The amplicon generated from the genomic DNA presented a fragment of 336 bp as well as the predicted fragment that exhibits IR. The expected amplicon for the processed transcript was a 198-bp DNA fragment. The results showed that IR was present in both the control conditions (96 h pre-culture [control − 0 h] and 12 h of culture in the absence of UDA [control − 12 h]) as well in 12-h cultures with 70% and 100% of the UDA minimum inhibitory concentration (MIC) (Figure 2B). These results indicate that retention of intron 1 from pakA gene is a process that occurs naturally in the assayed conditions. These results were confirmed in biological triplicates (Figure S2).
In silico analysis indicated that the retention of intron 1, present in the alternative pakA gene transcription, generates a premature stop codon (UAA) 1098 bases downstream of the start codon. This first mRNA sequence encodes a polypeptide with 336 amino acid residues, which differs from the N-terminal portion of the native PAKA protein only in the last 42 amino acids. Interestingly, 33 bases downstream from the first stop codon defines a new start codon followed by a sequence which encodes a protein with 638 amino acid residues. The alignment of this amino acid sequence indicated that it is identical to the C-terminal portion of the native PAKA protein, with both the catalytic site and the GBB domains (Figure 3) as well as the translational modification sites but without the regulatory CRIB domain.

2.3. UDA Affects the Transcription Process of the pakA/Ste20 Gene in T. rubrum

A set of quantitative RT-PCR (qRT-PCR) primers was designed to target within the sequence of intron 1 of the pakA gene to determine the influence of UDA in its retention (Figure 4A). The results showed no significant variation in retention of intron 1 during the first 12 h in the absence of UDA, considering the control − 0 h (just after the 96-h pre-cultivation, without UDA) as the reference condition. However, an increase in retention was observed after 24 h of cultivation in the same situation (Figure 4B). Conversely, there was a strong retention peak after 3 h of culture in the presence of the drug regardless of its concentration (Figure 4B). These results indicated that UDA modulates intron 1 retention in the pakA kinase gene of T. rubrum. Also, after 12 h of culturing in the presence of UDA, the retention levels remained statistically representative, maintaining the IR higher than in the control conditions, especially with 100% MIC (Figure 4B). After 24 h of cultivation, there was a relative stabilization in IR levels considering both culture situations (with and without UDA), but with higher IR levels in culture with UDA when compared with that of the control condition (control − 0 h) (Figure 4B).
We also analyzed the relative expression levels of the pakA/Ste20 kinase gene in T. rubrum under different culture conditions, using a set of qRT-PCR primers designed within the exon 2 sequence of this gene (Figure 5A). Interestingly, there was a decrease in pakA expression after culturing the fungus for 12 h, being more pronounced in the presence of UDA, suggesting that the drug forced a stronger down-regulation. This same UDA effect was observed after 24 h (Figure 5B). These results indicated that UDA affects the MAPK pathway somewhere up-stream of the STE20/PAKA signaling cascade or perhaps the transcription of some transcription factor involved in this pathway.

3. Discussion

MAPK is a component of an evolutionarily well-conserved signaling pathway, generally involved in a range of cellular process [16]. Different stimuli may activate its signaling cascade, leading to several cellular responses, such as cytoskeletal architecture, cytokinesis, motility, cell-cycle progression, hyphal growth, differentiation of asexual development structures, stress response, sterol homeostasis, and virulence [8,17,18]. In this context, STE/PAK kinase has a fundamental role. In yeast, genetic evidence suggests that Ste20p functions as a G-protein (βγ subunit) effector of downstream signaling components, such as Ste11p [19,20]. Thus, it is possible that this signaling cascade model can be extrapolated to the T. rubrum MAPK pathway because in silico data indicated high identity of STE20/PAKA kinase with its Saccharomyces cerevisiae orthologue (72%). We also observed the presence of the GBB domain in the T. rubrum STE20/PAKA C-terminal region (Figure 1C). This specific domain was described as a linker of Ste20p to the Gβ subunit of the S. cerevisiae heterotrimeric G-protein [21]. It seems that STE20/PAKA is quite a promiscuous molecule since it displays regions rich in specific amino acids (Figure S1A) which enable its interaction with a range of other regulatory proteins and with molecules with different functions. This type of interaction has been described for the kinase signaling pathway activation. For instance, the serine-rich domain of the interleukin-2 receptor β chain is crucial for the activation of both tyrosine kinase and phosphatidylinositol-3-kinase in the transduction of the human IL-2 cellular proliferation signal [22,23,24]. Moreover, T. rubrum STE20/PAKA showed several sites for post-translational modifications (PTMs), such as phosphorylation, Myr, N-gly, and AMD (Figure S1B–D). These PTMs may allow a wide range of possible three-dimensional configurations and an excellent mechanism for the regulation of catalytic activity of the protein, being able to perform several biological functions in response to specific physiological requirements [25,26,27]. For instance, proper protein glycosylation was necessary to promote MAPK signal fidelity in S. cerevisiae [16]. Myr has been associated with protein kinase activation/deactivation, membrane interaction, or localization of the enzyme near membrane-bound substrates [28,29,30,31]. AMD may modify the physicochemical properties of the bioactive peptide, being essential for receptor recognition, signal transduction, and other biological functions [32,33,34]. Among the PTMs, phosphorylation is the most intriguing, considering that the activation of most kinase proteins depends on the autophosphorylation of their activation loop. This modification induces the conversion from an inactive to active conformation, and it still occurs in the protein’s inactive state [35]. As reported in S. cerevisiae, the binding of Cdc42p to the Ste20p CRIB site releases the kinase domain, activating its autophosphorylation process and consequently its catalytic activation [8]. Conversely, mammalian PAK2 activation occurs by the caspase 3 cleavage mechanism. This process generates two PAK2 (p27 and p34) fragments, with the entire kinase catalytic domain in the p34 fragment, free from the inhibitory CRIB domain. Consequently, this enzyme becomes active after its autophosphorylation [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36].
IR is the leading AS mechanism in fungi, with critical roles in different fungal aspects, such as development, reproduction, mitochondrial functions, and signal transduction. However, its contribution to fungal biology still requires further studies [37,38,39,40]. Recently, using the same RNA-Seq data of the present work, our group published AS events in a T. rubrum strain. An IR event was presented in the gene encoding the enzyme inosine monophosphate dehydrogenase (IMPDH; TERG_06846), which is involved in purine biosynthesis, as confirmed by qRT-PCR [41]. However, in this case, and in most circumstances described in fungi, the IR process generated an altered reading frame code and, consequently, smaller proteins with disrupted functional domains. Thus, a premature stop codon may activate the nonsense-mediated decay (NMD) pathway, thereby resulting in mRNA degradation before its translation or the formation of truncated proteins [12].
Conversely, functional alternative splicing has been described in fungi. For instance, two functional serine/threonine phosphatase (PTC7) isoforms can be produced by IR in S. cerevisiae strains. In this model, both proteins showed distinct functions, working in different cellular compartments, specifically in mitochondria and the nuclear envelope for the spliced and unspliced mRNA, respectively [42].
Here, we hypothesize a novel mechanism for the activation of T. rubrum STE20/PAKA through alternative pre-mRNA splicing. We considered that, because of alternative splicing, the translation of STE20/PAKA free from its autoinhibitory CRIB domain is possible (Figure 3). This hypothesis is supported by the RNA-Seq data that presented isoform transcripts with intron 1 retained, validated through qRT-PCR (Figure 2). Furthermore, although the T. rubrum Ste20/pakA kinase alternative transcript displays a premature stop codon, the subsequent reading frame is conserved. Thus, it is plausible that this mRNA may be preserved from degradation and subsequently translated, producing two independent proteins/polypeptides, because bicistronic and polycistronic transcripts have already been described in fungi [43,44]. Our data also suggested that IR by T. rubrum Ste20/pakA is a natural process, considering the transcription of both isoforms in all conditions assayed (Figure 2B). Otherwise, this process can be altered or modulated by the antifungal agent UDA, since its presence affected both the splicing of T. rubrum Ste20/pakA mRNA and the relative expression of this gene (Figure 4 and Figure 5). According to the in silico analysis, both STE20/PAKA isoforms appear to be transient molecules, likely with a rapid protein turnover. This hypothesis is based on the presence of the putative PEST domain in both CRIB and catalytic kinase subunits (Figure S1E), since the PEST sequence was postulated to be involved in protein degradation [45]. Moreover, this evidence could support a quick change in cellular response and an even faster change in the signaling pattern. Thus, the data presented here indicate that the retention of intron-1 in T. rubrum is modulated by external stressors, such as the presence of antifungal drugs and nutrient availability. The data suggest that in the presence of UDA, fungal cells strongly require the PAKA enzyme without the CRIB domain, probably for a faster response to environmental changes. After 12 h, several drug tolerance genes probably have already been expressed, and the intron-1 retention levels decrease, maintaining a basic level since other stressors, such as the availability of nutrients, also come into effect (Figure 4).
Besides describing a novel biological mechanism of enzymatic activation by AS, these findings suggest that Ste20/PAKA kinase is a crucial molecule in the processes of adaptation and interaction with the host tissue. Moreover, the results suggest that PAKA is a critical virulence factor for T. rubrum, as described for Candida albicans, which, in its pathogenic form, depends on its STE20 homolog (CST20) for hyphal growth, colonization of the host tissue, and consequent virulence [46,47]. All these pieces of evidence suggest that gene expression and its modifications likely occur for more efficient signaling transduction and, consequently, an appropriate cellular response in a challenging and constantly changing environment. Thus, STE20/PAKA may be a relevant target for the development of novel and specific drugs against dermatophytes. However, a complete understanding of the alternative splicing, as a cellular activation mechanism, requires further experimental substantiation. Heterologous expression and biochemical characterization of both STE20/PAKA kinase isoforms may be fundamental steps to achieve this goal.

4. Materials and Methods

4.1. Strains and Culture Conditions

T. rubrum strain CBS 118892, obtained from the Centraalbureau voor Schimmelcultures in the Netherlands, was cultured as previously described [41]. All assays were carried out in triplicate.

4.2. RNA Extraction, cDNA Library Construction, and High-Throughput Sequencing

RNA extraction, cDNA library construction, and high-throughput sequencing were conducted as previously described [41].

4.3. Data Analysis

The reads were quality filtered using the software FastQC, aligned with the T. rubrum genome available in the Broad Institute’s Dermatophyte Comparative Database (ftp://ftp.broadinstitute.org/pub/annotation/fungi) (accessed on 5 August 2017), using the Bowtie2 algorithm [48], and inspected with the software Integrative Genomics Viewer (IGV) [49]. Furthermore, the independent filter, which precedes the statistical filter, was applied to identify gene-by-gene variances among biological replicates, followed by FDR [50]. These analyses were performed using the DESeq package and manipulated in the R statistical environment [51].
A transcriptome-wide survey of alternative splicing in differentially expressed genes after exposure of T. rubrum to UDA was conducted using analysis of IR. This analysis was performed through ad hoc Perl scripts, using only the single-end libraries. Introns were identified in the reference genome and aligned against the sequenced libraries, and read counts were analyzed for each retained intron [41]. qRT-PCR was used to confirm the properties of genes of interest (primers are listed in Table S1). RNA-Seq data were deposited in the Gene Expression Omnibus (GEO) [52] database under accession number GSE102872 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc= GSE102872) (accessed on 6 February 2018).

4.4. RT-PCR and qRT-PCR Analysis

For gene analysis assays, about 1 × 106 conidia were inoculated into 100 mL of Sabouraud medium and incubated at 28 °C under shaking for 96 h. The resulting mycelia were transferred to Sabouraud medium without drugs (control) or in the presence of UDA (11.25 µg/mL or 25 μg/mL). Then, the mycelia were incubated for 3 h, 12 h, and 24 h at 28 °C, under constant agitation. Gene expression was quantified by qPCR with the StepOnePlus Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). Independent samples from each time point were used, with biological and technical replication performed in triplicate. Each qPCR reaction was performed in a final volume of 12.5 μL, containing 1 μL primer, 6.25 μL SYBR Green PCR Master Mix (Applied Biosystems, Warrington Cheshire, UK), and 50 ng cDNA. The cycling conditions included an initial PCR step of 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Expression was assessed based on the relative quantification of responsive genes using the 2−∆∆Ct method [53]. Data were normalized with two endogenous controls, DNA-dependent RNA polymerase II (rpb II) and DNA topoisomerase II [54]. The reference sample was taken at 0 h. Statistical analysis was performed by one-way analysis of variance (ANOVA) with Tukey’s post hoc test using Graph Pad Prism v. 5.1 software. Significance was indicated by p < 0.05.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/11/3654/s1. Figure S1: Graphical representation of functional STE20/PAKA protein. Figure S2: Agarose gel electrophoresis (1.5%) run with PCR products from three biological replicates for validation of the pakA/Ste20 intron-1 retention. Table S1: Primers used for qRT-PCR.

Author Contributions

Conceptualization and methodology, N.S.M., E.V.G., N.M.M.-R. and A.R.; Software and bioinformatics analysis, P.R.S., Validation, E.V.G. and J.C.B., Formal analysis, E.V.G., J.C.B. and A.R., writing, E.V.G., N.M.M.-R. and A.R., supervision, N.M.M.-R. and A.R., project administration and funding acquisition, N.M.M.-R. and A.R.

Funding

This work was supported by grants from the following Brazilian Funding Agencies: São Paulo Research Foundation (FAPESP), grant number 2014/03847-7, and Postdoctoral Fellowship to EVG No. 2016/04274-6; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), grant numbers 305797/2017-4 and 304989/2017-7; Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)—Finance Code 001; and Fundação de Apoio ao Ensino, Pesquisa e Assistência (FAEPA).

Acknowledgments

We thank V.M. Oliveira, M. Mazucato, and M.D. Martins for technical support.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AMDAmidation
CRIBCdc42/Rac interactive binding
FDRFalse discovery rate
GBBGβ binding
IRIntron retention
MAPKMitogen-activated protein kinase
MICMinimum inhibitory concentration
MYRN-myristoylation
N-glyN-glycosylation
PAKp21 activated kinase
PTMsPost-translational modifications
RNA-SeqRNA-Sequencing
RT-PCRReverse transcription polymerase chain reaction
qRT-PCRQuantitative reverse transcription-PCR
UDAUndecanoic acid

References

  1. Bahn, Y.S.; Xue, C.; Idnurm, A.; Rutherford, J.C.; Heitman, J.; Cardenas, M.E. Sensing the environment: Lessons from fungi. Nat. Rev. Microbiol. 2007, 5, 57–69. [Google Scholar] [CrossRef] [PubMed]
  2. Zeilinger, S.; Gruber, S.; Bansal, R.; Mukherjee, P.K. Secondary metabolism in Trichoderma–Chemistry meets genomics. Fungal Biol. Rev. 2016, 30, 74–90. [Google Scholar] [CrossRef]
  3. Kosti, I.; Mandel-Gutfreund, Y.; Glaser, F.; Horwitz, B.A. Comparative analysis of fungal protein kinases and associated domains. BMC Genomics 2010, 11, 133. [Google Scholar] [CrossRef] [PubMed]
  4. Widmann, C.; Gibson, S.; Jarpe, M.B.; Johnson, G.L. Mitogen-activated protein kinase: Conservation of a three-kinase module from yeast to human. Physiol. Rev. 1999, 79, 143–180. [Google Scholar] [CrossRef] [PubMed]
  5. Román, E.; Arana, D.M.; Nombela, C.; Alonso-Monge, R.; Pla, J. MAP kinase pathways as regulators of fungal virulence. Trends Microbiol. 2007, 15, 181–190. [Google Scholar] [CrossRef] [PubMed]
  6. Manser, E.; Leung, T.; Salihuddin, H.; Zhao, Z.S.; Lim, L. A brain serine/threonine protein kinase activated by Cdc42 and Rac1. Nature 1994, 367, 40–46. [Google Scholar] [CrossRef] [PubMed]
  7. Johnson, G.L.; Lapadat, R. Mitogen-activated protein kinase pathways mediated by ERK, JNK, and p38 protein kinases. Science 2002, 298, 1911–1912. [Google Scholar] [CrossRef] [PubMed]
  8. Boyce, K.J.; Andrianopoulos, A. Ste20-related kinases: Effectors of signaling and morphogenesis in fungi. Trends Microbiol. 2011, 19, 400–410. [Google Scholar] [CrossRef] [PubMed]
  9. Rudel, T.; Bokoch, G.M. Membrane and morphological changes in apoptotic cells regulated by caspase-mediated activation of PAK2. Science 1997, 276, 1571–1574. [Google Scholar] [CrossRef] [PubMed]
  10. Zhao, R.; Davey, M.; Hsu, Y.C.; Kaplanek, P.; Tong, A.; Parsons, A.B.; Krogan, N.; Cagney, G.; Mai, D.; Greenblatt, J.; et al. Navigating the chaperone network: An integrative map of physical and genetic interactions mediated by the Hsp90 chaperone. Cell 2005, 120, 715–727. [Google Scholar] [CrossRef] [PubMed]
  11. Hsu, Y.H.; Johnson, D.A.; Traugh, J.A. Analysis of conformational changes during activation of protein kinase Pak2 by amide hydrogen/deuterium exchange. J. Biol. Chem. 2008, 283, 36397–363405. [Google Scholar] [CrossRef] [PubMed]
  12. Gonzalez, C.I.; Bhattacharya, A.; Wang, W.; Peltz, S.W. Nonsense-mediated mRNA decay in Saccharomyces cerevisiae. Gene 2001, 274, 15–25. [Google Scholar] [CrossRef]
  13. Vermout, S.; Tabart, J.; Baldo, A.; Mathy, A.; Losson, B.; Mignon, B. Pathogenesis of dermatophytosis. Mycopathologia 2008, 166, 267–275. [Google Scholar] [CrossRef] [PubMed]
  14. De Castro, E.; Sigrist, C.J.; Gattiker, A.; Bulliard, V.; Langendijk-Genevaux, P.S.; Gasteiger, E.; Bairoch, A.; Hulo, N. ScanProsite: Detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006, 34, W362–W365. [Google Scholar] [CrossRef] [PubMed]
  15. Castro, O.; Movsichoff, F.; Parodi, A.J. Preferential transfer of the complete glycan is determined by the oligosaccharyltransferase complex and not by the catalytic subunit. Proc. Natl. Acad. Sci. USA 2006, 103, 14756–14760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Lien, E.C.; Nagiec, M.J.; Dohlman, H.G. Proper protein glycosylation promotes mitogen-activated protein kinase signal fidelity. Biochemistry 2013, 52, 115–124. [Google Scholar] [CrossRef] [PubMed]
  17. Sells, M.A.; Chernoff, J. Emerging from the Pak: The p21-activated protein kinase family. Trends Cell Biol. 1997, 7, 162–167. [Google Scholar] [CrossRef]
  18. Lin, M.; Grillitsch, K.; Daum, G.; Just, U.; Hofken, T. Modulation of sterol homeostasis by the Cdc42p effectors Cla4p and Ste20p in the yeast Saccharomyces cerevisiae. FEBS J. 2009, 276, 7253–7264. [Google Scholar] [CrossRef] [PubMed]
  19. Leberer, E.; Dignard, D.; Harcus, D.; Thomas, D.Y.; Whiteway, M. The Protein-Kinase Homolog Ste20p Is Required to Link the Yeast Pheromone Response G-Protein Beta-Gamma Subunits to Downstream Signaling Components. EMBO J. 1992, 11, 4815–4824. [Google Scholar] [CrossRef] [PubMed]
  20. Van Drogen, F.; O’Rourke, S.M.; Stucke, V.M.; Jaquenoud, M.; Neiman, A.M.; Peter, M. Phosphorylation of the MEKK Ste11p by the PAK-like kinase Ste20p is required for MAP kinase signaling in vivo. Curr. Biol. 2000, 10, 630–639. [Google Scholar] [CrossRef]
  21. Song, J.; Chen, Z.; Xu, P.; Gingras, R.; Ng, A.; Leberer, E.; Thomas, D.Y.; Ni, F. Molecular interactions of the Gbeta binding domain of the Ste20p/PAK family of protein kinases. An isolated but fully functional Gbeta binding domain from Ste20p is only partially folded as shown by heteronuclear NMR spectroscopy. J. Biol. Chem. 2001, 276, 41205–41212. [Google Scholar] [CrossRef] [PubMed]
  22. Merida, I.; Williamson, P.; Kuziel, W.A.; Greene, W.C.; Gaulton, G.N. The serine-rich cytoplasmic domain of the interleukin-2 receptor beta chain is essential for interleukin-2-dependent tyrosine protein kinase and phosphatidylinositol-3-kinase activation. J. Biol. Chem. 1993, 268, 6765–6770. [Google Scholar] [PubMed]
  23. Williamson, M.P. The structure and function of proline-rich regions in proteins. Biochem. J. 1994, 297, 249–260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Belvitch, P.; Adyshev, D.; Elangovan, V.R.; Brown, M.E.; Naureckas, C.; Rizzo, A.N.; Siegler, J.H.; Garcia, J.G.; Dudek, S.M. Proline-rich region of non-muscle myosin light chain kinase modulates kinase activity and endothelial cytoskeletal dynamics. Microvasc. Res. 2014, 95, 94–102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Walsh, L.; Hastwell, P.W.; Keenan, P.O.; Knight, A.W.; Billinton, N.; Walmsley, R.M. Genetic modification and variations in solvent increase the sensitivity of the yeast RAD54-GFP genotoxicity assay. Mutagenesis 2005, 20, 317–327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Audagnotto, M.; Dal Peraro, M. Protein post-translational modifications: In silico prediction tools and molecular modeling. Comput. Struct. Biotechnol. J. 2017, 15, 307–319. [Google Scholar] [CrossRef] [PubMed]
  27. Walsh, C.T.; Garneau-Tsodikova, S.; Gatto, G.J., Jr. Protein posttranslational modifications: The chemistry of proteome diversifications. Angew. Chem. Int. Ed. Engl. 2005, 44, 7342–7372. [Google Scholar] [CrossRef] [PubMed]
  28. Hantschel, O.; Nagar, B.; Guettler, S.; Kretzschmar, J.; Dorey, K.; Kuriyan, J.; Superti-Furga, G. A myristoyl/phosphotyrosine switch regulates c-Abl. Cell 2003, 112, 845–857. [Google Scholar] [CrossRef]
  29. Gaffarogullari, E.C.; Masterson, L.R.; Metcalfe, E.E.; Traaseth, N.J.; Balatri, E.; Musa, M.M.; Mullen, D.; Distefano, M.D.; Veglia, G. A myristoyl/phosphoserine switch controls cAMP-dependent protein kinase association to membranes. J. Mol. Biol. 2011, 411, 823–836. [Google Scholar] [CrossRef] [PubMed]
  30. Boutin, J.A. Myristoylation. Cell Signal 1997, 9, 15–35. [Google Scholar] [CrossRef]
  31. Patwardhan, P.; Resh, M.D. Myristoylation and membrane binding regulate c-Src stability and kinase activity. Mol. Cell Biol. 2010, 30, 4094–4107. [Google Scholar] [CrossRef] [PubMed]
  32. Eipper, B.A.; Stoffers, D.A.; Mains, R.E. The biosynthesis of neuropeptides: Peptide alpha-amidation. Annu. Rev. Neurosci. 1992, 15, 57–85. [Google Scholar] [CrossRef] [PubMed]
  33. Chufán, E.E.; De, M.; Eipper, B.A.; Mains, R.E.; Amzel, L.M. Amidation of bioactive peptides: The structure of the lyase domain of the amidating enzyme. Structure 2009, 17, 965–973. [Google Scholar] [CrossRef] [PubMed]
  34. Cui, W.; Niu, S.; Zheng, L.; Hu, L.; Huang, T.; Gu, L.; Feng, K.; Zhang, N.; Cai, Y.; Li, Y. Prediction of protein amidation sites by feature selection and analysis. Mol. Genet. Genomics 2013, 288, 391–400. [Google Scholar] [CrossRef] [PubMed]
  35. Beenstock, J.; Mooshayef, N.; Engelberg, D. How Do Protein Kinases Take a Selfie (Autophosphorylate)? Trends Biochem. Sci. 2016, 41, 938–953. [Google Scholar] [CrossRef] [PubMed]
  36. Walter, B.N.; Huang, Z.; Jakobi, R.; Tuazon, P.T.; Alnemri, E.S.; Litwack, G.; Traugh, J.A. Cleavage and activation of p21-activated protein kinase gamma-PAK by CPP32 (caspase 3). Effects of autophosphorylation on activity. J. Biol. Chem. 1998, 273, 28733–28739. [Google Scholar] [CrossRef] [PubMed]
  37. Leal, J.; Squina, F.M.; Freitas, J.S.; Silva, E.M.; Ono, C.J.; Martinez-Rossi, N.M.; Rossi, A. A splice variant of the Neurospora crassa hex-1 transcript, which encodes the major protein of the Woronin body, is modulated by extracellular phosphate and pH changes. FEBS Lett. 2009, 583, 180–184. [Google Scholar] [CrossRef] [PubMed]
  38. Trevisan, G.L.; Oliveira, E.H.; Peres, N.T.A.; Cruz, A.H.; Martinez-Rossi, N.M.; Rossi, A. Transcription of Aspergillus nidulans pacC is modulated by alternative RNA splicing of palB. FEBS Lett. 2011, 585, 3442–3445. [Google Scholar] [CrossRef] [PubMed]
  39. Mendes, N.S.; Silva, P.M.; Silva-Rocha, R.; Martinez-Rossi, N.M.; Rossi, A. Pre-mRNA splicing is modulated by antifungal drugs in the filamentous fungus Neurospora crassa. FEBS Open Bio 2016, 6, 358–368. [Google Scholar] [CrossRef] [PubMed]
  40. Jin, L.; Li, G.; Yu, D.; Huang, W.; Cheng, C.; Liao, S.; Wu, Q.; Zhang, Y. Transcriptome analysis reveals the complexity of alternative splicing regulation in the fungus Verticillium dahliae. BMC Genomics 2017, 18, 130. [Google Scholar] [CrossRef] [PubMed]
  41. Mendes, N.S.; Bitencourt, T.A.; Sanches, P.R.; Silva-Rocha, R.; Martinez-Rossi, N.M.; Rossi, A. Transcriptome-wide survey of gene expression changes and alternative splicing in Trichophyton rubrum in response to undecanoic acid. Sci. Rep. 2018, 8, 38–46. [Google Scholar] [CrossRef] [PubMed]
  42. Juneau, K.; Nislow, C.; Davis, R.W. Alternative splicing of PTC7 in Saccharomyces cerevisiae determines protein localization. Genetics 2009, 183, 185–194. [Google Scholar] [CrossRef] [PubMed]
  43. Pelechano, V.; Wei, W.; Steinmetz, L.M. Extensive transcriptional heterogeneity revealed by isoform profiling. Nature 2013, 497, 127–131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Gordon, S.P.; Tseng, E.; Salamov, A.; Zhang, J.; Meng, X.; Zhao, Z.; Kang, D.; Underwood, J.; Grigoriev, I.V.; Figueroa, M.; et al. Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule mRNA Sequencing. PLoS ONE 2015, 10, e0132628. [Google Scholar] [CrossRef] [PubMed]
  45. Rogers, S.; Wells, R.; Rechsteiner, M. Amino acid sequences common to rapidly degraded proteins: The PEST hypothesis. Science 1986, 234, 364–368. [Google Scholar] [CrossRef] [PubMed]
  46. Leberer, E.; Harcus, D.; Broadbent, I.D.; Clark, K.L.; Dignard, D.; Ziegelbauer, K.; Schmidt, A.; Gow, N.A.R.; Brown, A.J.P.; Thomas, D.Y. Signal transduction through homologs of the Ste20p and Ste7p protein kinases can trigger hyphal formation in the pathogenic fungus Candida albicans. Proc. Natl. Acad. Sci. USA 1996, 93, 13217–13222. [Google Scholar] [CrossRef] [PubMed]
  47. Leberer, E.; Ziegelbauer, K.; Schmidt, A.; Harcus, D.; Dignard, D.; Ash, J.; Johnson, L.; Thomas, D.Y. Virulence and hyphal formation of Candida albicans require the Ste20p-like protein kinase CaCla4p. Curr. Biol. 1997, 7, 539–546. [Google Scholar] [CrossRef]
  48. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Thorvaldsdottir, H.; Robinson, J.T.; Mesirov, J.P. Integrative Genomics Viewer (IGV): High-performance genomics data visualization and exploration. Brief Bioinform. 2012, 14, 178–192. [Google Scholar] [CrossRef] [PubMed]
  50. Bourgon, R.; Gentleman, R.; Huber, W. Independent filtering increases detection power for high-throughput experiments. Proc. Natl. Acad. Sci. USA 2010, 107, 9546–9551. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed]
  52. Edgar, R.; Domrachev, M.; Lash, A.E. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30, 207–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Schefe, J.H.; Lehmann, K.E.; Buschmann, I.R.; Unger, T.; Funke-Kaiser, H. Quantitative real-time RT-PCR data analysis: Current concepts and the novel “gene expression’s CT difference” formula. J. Mol. Med. 2006, 84, 901–910. [Google Scholar] [CrossRef] [PubMed]
  54. Jacob, T.R.; Peres, N.T.; Persinoti, G.F.; Silva, L.G.; Mazucato, M.; Rossi, A.; Martinez-Rossi, N.M. rpb2 is a reliable reference gene for quantitative gene expression analysis in the dermatophyte Trichophyton rubrum. Med. Mycol. 2012, 50, 368–377. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Graphical representation of complete pakA/Ste20 gene transcription/splicing/translation in T. rubrum. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow) and four introns (black); (B) Graphical representation of mRNA transcribed/spliced from the pakA/Ste20 kinase gene; (C) Graphical representation of functional STE20/PAKA protein, with its respective domains, translated from processed mRNA.
Figure 1. Graphical representation of complete pakA/Ste20 gene transcription/splicing/translation in T. rubrum. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow) and four introns (black); (B) Graphical representation of mRNA transcribed/spliced from the pakA/Ste20 kinase gene; (C) Graphical representation of functional STE20/PAKA protein, with its respective domains, translated from processed mRNA.
Ijms 19 03654 g001
Figure 2. Intron retention analysis of the T. rubrum pakA/Ste20 kinase gene. (A) Graphical representation of the pakA/Ste20 kinase gene indicating its size (base pairs), five exons (dark yellow), four introns (black), and the primer-binding sites used for amplification of the region flanking intron 1 (blue and orange arrows); (B) Agarose gel electrophoresis (1.5%) run with qRT-PCR products from the pakA/Ste20 intron 1 region. C: culture without the presence of UDA; G: amplification from genomic DNA (gene amplification positive control); 0 h: amplification product from cDNA obtained from culture in medium lacking the drug; 70%: amplification product from cDNA obtained from culture in medium containing 70% of the UDA MIC; 100%: amplification product from cDNA obtained from culture in medium containing 100% of the UDA MIC; MM: molecular weight marker (50 bp).
Figure 2. Intron retention analysis of the T. rubrum pakA/Ste20 kinase gene. (A) Graphical representation of the pakA/Ste20 kinase gene indicating its size (base pairs), five exons (dark yellow), four introns (black), and the primer-binding sites used for amplification of the region flanking intron 1 (blue and orange arrows); (B) Agarose gel electrophoresis (1.5%) run with qRT-PCR products from the pakA/Ste20 intron 1 region. C: culture without the presence of UDA; G: amplification from genomic DNA (gene amplification positive control); 0 h: amplification product from cDNA obtained from culture in medium lacking the drug; 70%: amplification product from cDNA obtained from culture in medium containing 70% of the UDA MIC; 100%: amplification product from cDNA obtained from culture in medium containing 100% of the UDA MIC; MM: molecular weight marker (50 bp).
Ijms 19 03654 g002
Figure 3. Graphical representation of alternative pakA/Ste20 gene transcription/splicing/translation in T. rubrum. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow) and four introns (black); (B) Graphical representation of mRNA transcribed from the pakA/Ste20 kinase gene with retention of intron 1; (C) Graphical representation of both probable functional STE20/PAKA proteins, with their respective domains, translated from intron-retained mRNA.
Figure 3. Graphical representation of alternative pakA/Ste20 gene transcription/splicing/translation in T. rubrum. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow) and four introns (black); (B) Graphical representation of mRNA transcribed from the pakA/Ste20 kinase gene with retention of intron 1; (C) Graphical representation of both probable functional STE20/PAKA proteins, with their respective domains, translated from intron-retained mRNA.
Ijms 19 03654 g003
Figure 4. qRT-PCR analysis of the intron 1 retention levels of the T. rubrum pakA/ste20 gene. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow), four introns (black), and the primer-binding sites used for amplification of the intron 1 retention (blue and orange arrows). (B) Intron 1 retention levels of the T. rubrum pakA/ste20 gene, incubated in the presence of 70% and 100% UDA MIC for 3, 12, and 24 h using ANOVA followed by Tukey’s post hoc test (* p < 0.05; *** p < 0.001).
Figure 4. qRT-PCR analysis of the intron 1 retention levels of the T. rubrum pakA/ste20 gene. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow), four introns (black), and the primer-binding sites used for amplification of the intron 1 retention (blue and orange arrows). (B) Intron 1 retention levels of the T. rubrum pakA/ste20 gene, incubated in the presence of 70% and 100% UDA MIC for 3, 12, and 24 h using ANOVA followed by Tukey’s post hoc test (* p < 0.05; *** p < 0.001).
Ijms 19 03654 g004
Figure 5. qRT-PCR analysis of the T. rubrum pakA/ste20 gene expression levels. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow), four introns (black), and the primer-binding sites used for gene expression analysis (blue and orange arrows). (B) Expression levels of the T. rubrum pakA/ste20 gene, incubated in the presence of 70% and 100% UDA MIC for 3, 12, and 24 h using ANOVA followed by Tukey’s post hoc test (** p < 0.01; *** p < 0.001).
Figure 5. qRT-PCR analysis of the T. rubrum pakA/ste20 gene expression levels. (A) Graphical representation of the pakA/Ste20 kinase gene with its five exons (dark yellow), four introns (black), and the primer-binding sites used for gene expression analysis (blue and orange arrows). (B) Expression levels of the T. rubrum pakA/ste20 gene, incubated in the presence of 70% and 100% UDA MIC for 3, 12, and 24 h using ANOVA followed by Tukey’s post hoc test (** p < 0.01; *** p < 0.001).
Ijms 19 03654 g005

Share and Cite

MDPI and ACS Style

Gomes, E.V.; Bortolossi, J.C.; Sanches, P.R.; Mendes, N.S.; Martinez-Rossi, N.M.; Rossi, A. STE20/PAKA Protein Kinase Gene Releases an Autoinhibitory Domain through Pre-mRNA Alternative Splicing in the Dermatophyte Trichophyton rubrum. Int. J. Mol. Sci. 2018, 19, 3654. https://doi.org/10.3390/ijms19113654

AMA Style

Gomes EV, Bortolossi JC, Sanches PR, Mendes NS, Martinez-Rossi NM, Rossi A. STE20/PAKA Protein Kinase Gene Releases an Autoinhibitory Domain through Pre-mRNA Alternative Splicing in the Dermatophyte Trichophyton rubrum. International Journal of Molecular Sciences. 2018; 19(11):3654. https://doi.org/10.3390/ijms19113654

Chicago/Turabian Style

Gomes, Eriston V., Julio C. Bortolossi, Pablo R. Sanches, Niege S. Mendes, Nilce M. Martinez-Rossi, and Antonio Rossi. 2018. "STE20/PAKA Protein Kinase Gene Releases an Autoinhibitory Domain through Pre-mRNA Alternative Splicing in the Dermatophyte Trichophyton rubrum" International Journal of Molecular Sciences 19, no. 11: 3654. https://doi.org/10.3390/ijms19113654

APA Style

Gomes, E. V., Bortolossi, J. C., Sanches, P. R., Mendes, N. S., Martinez-Rossi, N. M., & Rossi, A. (2018). STE20/PAKA Protein Kinase Gene Releases an Autoinhibitory Domain through Pre-mRNA Alternative Splicing in the Dermatophyte Trichophyton rubrum. International Journal of Molecular Sciences, 19(11), 3654. https://doi.org/10.3390/ijms19113654

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop