Revista UNICIENCIA
Uniciencia Vol. 38(1), January-December, 2024
E-ISSN: 2215-3470
DOI: https://dx.doi.org/10.15359/ru.38-1.16
Cloning and Characterization of the Gene Encoding the Heat Shock Protein HSP83 from Trypanosoma cruzi
Clonación y caracterización del gen codificante de la proteína de choque térmico HSP83 de Trypanosoma cruzi
Clonagem e caracterização do gene que codifica a proteína de choque térmico HSP83 do Trypanosoma cruzi
Ana Fernández1, 2, Ana Rita De Lima1, Elizabeth Ferrer2*
Received: Aug/15/2023 • Accepted: Mar/7/2024 • Published: Jul/31/2024
Abstract [Objective] Trypanosoma cruzi, a causal agent of Chagas’ disease, is a parasite whose life cycle alternates between an invertebrate (triatomine) and a vertebrate (mammal) host. Various studies have shown that in T. cruzi, HSP83 (a homolog of HSP90) is essential for cell division and control of the response to thermal stress. This investigation focused on studying the cloning, bioinformatic characterization, and expression of the heat shock protein HSP83 T. cruzi gene for further cell signaling studies. [Methodology] RNA was extracted from T. cruzi epimastigotes (EPm6 clone, MHOM/VE/2007/ 6c) using a commercial kit. The cDNA encoding HSP83 was determined using RT-PCR on extracted mRNA, for which the primers were designed based on the HSP83 sequence of T. cruzi strain CL Brener. Cloning was performed using pGEM®T-Easy and subcloned into the expression vector pQE30. Sequence and bioinformatic characterization were performed. The gene was expressed, and the recombinant protein was purified using affinity chromatography and identified through immunoblotting. [Results] Sequence analysis showed similarity to the gene encoding HSP83 from Trypanosoma cruzi, and HSP domains and B epitopes in the sequence were also observed. After 3 hours of induction with IPTG, a recombinant protein with an approximate weight of 83 kDa was obtained. The immunoblotting reaction with hyperimmune anti-T. cruzi epimastigote serum helped detect a single band with a molecular weight of nearly 83 kDa. [Conclusions] All results indicate that the cloning and characterization of HSP83 from Trypanosoma cruzi was achieved. Keywords: Trypanosoma cruzi; heat shock protein; HSP83; cloning; expression. Resumen [Objetivo] Tryapanosoma cruzi, un agente causal de la enfermedad de Chagas, es un parásito que presenta una alternancia en su ciclo de vida entre un hospedador invertebrado (triatomino) y uno vertebrado (mamífero). Varios estudios han demostrado que en T. cruzi, HSP83 (homólogo de HSP90) es esencial para la división celular y el control de la respuesta al estrés térmico. Esta investigación se concentró en el estudio de la clonación, la caracterización bioinformática y la expresión del gen de la proteína de choque térmico HSP83 de T. cruzi para estudios posteriores sobre señalización. [Metodología] El ARN fue extraído de epimagostigotes T. cruci (clon EPm6/MOHM/VE/2007/ 6c), utilizando un kit comercial. El ADNc que codifica HSP83 se obtuvo utilizando RT-PCR, a partir del ARNm extraído, para lo cual se diseñaron los cebadores con base en la secuencia HSP83 de la cepa T. cruzi CL Brener. La clonación se realizó usando pGEM®T-Easy y se subclonó en el vector de expresión pQE30. Se efectuó la caracterización bioinformática y de secuencias. El gen se expresó y la proteína recombinante se purificó por cromatografía de afinidad, igual que se identificó por inmunotransferencia. [Resultados] El análisis de secuencia mostró similitud con el gen HSP83 de Trypanosoma cruzi y se observaron dominios HSP, así como epítopos B en la secuencia. Después de 3 horas de inducción con IPTG, se obtuvo una proteína recombinante con un peso aproximado de 83 kDa. La reacción de inmunotransferencia con suero hiperinmune anti-epimastigote de T. cruzi permitió la detección de una sola banda con un peso molecular de aproximadamente 83 kDa. [Conclusiones] Todos los resultados indican que se logró la clonación y caracterización de HSP83 de Trypanosoma cruzi. Keywords: Trypanosoma cruzi; proteína de choque térmico; HSP83; clonación; expresión. Resumo [Objetivo] O Tryapanosoma cruzi, agente causador da doença de Chagas, é um parasita que apresenta uma alternância em seu ciclo de vida entre um hospedeiro invertebrado (triatomíneo) e um vertebrado (mamífero). Vários estudos demonstraram que, no T. cruzi, a HSP83 (homóloga da HSP90) é essencial para a divisão celular e o controle da resposta ao estresse térmico. Esta pesquisa concentrou-se no estudo da clonagem, caracterização bioinformática e expressão do gene da proteína de choque térmico HSP83 do T. cruzi para estudos adicionais de sinalização. [Metodologia] O RNA foi extraído de epimastigotas de T. cruzi (clone EPm6/MOHM/VE/2007/ 6c) usando um kit comercial. O cDNA que codifica a HSP83 foi obtido por meio de RT-PCR a partir do mRNA extraído, para o qual os primers foram projetados com base na sequência da HSP83 da cepa T. cruzi CL Brener. A clonagem foi realizada com o pGEM®T-Easy e subclonada no vetor de expressão pQE30. Foi realizada a caracterização bioinformática e de sequência. O gene foi expresso e a proteína recombinante foi purificada por cromatografia de afinidade e identificada por imunotransferência. [Resultados] A análise da sequência mostrou similaridade com o gene HSP83 do Trypanosoma cruzi e os domínios HSP, bem como os epítopos B, foram observados na sequência. Após 3 horas de indução com IPTG, foi obtida uma proteína recombinante com peso aproximado de 83 kDa. A imunotransferência com soro hiperimune antiepimastigota de T. cruzi permitiu a detecção de uma única banda com peso molecular de aproximadamente 83 kDa. [Conclusões] Todos os resultados indicam que a clonagem e a caracterização da HSP83 do Trypanosoma cruzi foram alcançadas. Palavras-chave: Trypanosoma cruzi; proteína de choque térmico; HSP83; clonagem; expressão. |
Trypanosoma cruzi, a causal agent of Chagas’ disease, is a parasite that presents alternation in its life cycle between an invertebrate (triatomine) and vertebrate (mammal) host. Individuals with Chagas’ disease may present cardiological, digestive, and neurological symptoms and might die. Currently, 6 to 7 million people are affected by the disease, and 25 million are at risk of infection (WHO, 2023). In Venezuela, increases in prevalence and oral transmission outbreaks have been described in several states (Alarcón de Noya et al., 2015; Añez et al., 2020).
The parasite response is essential for viability in the heat and osmotic stress during the transition from the insect to the mammalian host, as well as when invading mammalian cells. The contribution of heat shock proteins in morphological and functional changes during stress situations has been implicated. This class of proteins is divided into several conserved families according to molecular weight and broad function, among which the five major families are HSP104, HSP90, HSP70, chaperonins, and small HSPs (Folgueira & Requena, 2007; Urményi et al., 2014). They participate in protein folding and cell signaling processes and play an important role in the host’s immune response. They have been identified as dominant antigens in diseases caused by protozoa (De Andrade et al., 1992; Requena et al., 1993; Angel et al., 1996; Menezes-Souza et al., 2014) and helminths (Rothstein et al., 1989; Moser et al., 1990; Colebrook & Lightowlers, 1997; Ferrer et al., 2005).
There is evidence that indicates that in T. cruzi Hsp are constitutively expressed and induced by different types of stress. In particular, 90 kDa heat shock proteins (HSP90) are involved in numerous intracellular processes in eukaryotic cells. Studies have shown that in T. cruzi HSP83 (homolog of HSP90) is essential for cell division and control of the response to thermal stress (Jhonson, 2012). The overexpression and the participation of HSP83 in the differentiation processes of T. cruzi have been suggested, and the inhibition of HSP83 of T. cruzi produces arrest in proliferation, demonstrating the importance in the control of the cell cycle (Graefe et al., 2002).
To date, there are merely two old reports about cloning this molecule (Dragon et al., 1987; Nadeau et al., 1992), and only its genomic structure and ATPase activity have been studied. This study aimed to perform the cloning, bioinformatics characterization, and expression of the heat shock protein HSP83 of T. cruzi (EPm6 clon, MHOM/VE/2007/Elpidio Padrón 6c) coding gene. The study would enable subsequent studies on the adaptation processes of the parasite to different microenvironments, the identification of new therapeutic or diagnostic targets, and the study of host-parasite interaction in Chagas disease.
Parasite Material
T. cruzi epimastigotes obtained from a human Chagas disease case (EPm6 clon, MHOM/VE/2007/Elpidio Padrón 6c) were cultivated at 27 ºC in liver infusion tryptose broth (LITB) (Contreras et al., 1994). Parasites were counted employing a Neubauer chamber and phase contrast microscope (Nikkon, Optiphot). They were adjusted to 1.5 x 109 epimastigotes/mL. Subsequently, the parasites were pelleted by centrifugation at 8,000 rpm for 15 minutes at four °C (IEC Centrifuge, B22M) and washed twice in 0.15 M phosphate saline buffer (PBS) pH 7.2 (57 mM Na2HP04, 18 mM KH2P04, 76.9 mMNaCl) to eliminate the remains of the culture medium. Finally, they were stored at -80 ºC until their use.
RNA Extraction and the Complementary DNA (cDNA) Obtaining
RNA extraction was carried out through the FavorPrep™ Viral Nucleic Acid Extraction Kit (Favorgene), where the parasite masses were previously treated with lysis buffer supplemented with Proteinase K (0.45 mg/mL) and RNase (80 U/µL). Then, the extraction continued with the instructions described in the commercial kit. The complementary DNA (cDNA) was obtained by implementing the Reverse Transcription-Polymerase Chain Reaction technique and using the Access RT-PCR System (Promega) commercial kit. Specific primers were designed manually using BioEdit® Software v7.1.11 (Hall, 2013) and, as a template sequence, the hsp83 of a reference strain of T. cruzi of 2115 bp (Trypanosoma cruzi strain CL Brener hsp83) obtained from GenBank (Accession No. XM_809799) strain used in the T. cruzi genome project (El-Sayed et al., 2005). The primers used were TcHSP83-D (5’- ATGACCGAGACATTCGCA -3’) and TcHSP83-R (5’- CTAGTCAACCTGCTCCAT -3’). The reaction mixture consisted of AMV/Tfl reaction buffer, deoxynucleotide triphosphates (0.2 mM), forward and reverse primers (1 µM) each, MgSO4 (1 mM), AMV Reverse Transcriptase (0.1 U /µL), Thermus flavus (Tfl) DNA polymerase (0.1 U / µL), 10 µL of the extracted RNA sample, resulting in a final reaction volume of 50 µL. The amplification reaction was carried out in an automatic C1000™ Thermal Cycler (BioRad) using the amplification program. The reverse transcription was conducted at 45 °C x 45 min; the initial denaturation, at 94 °C x 2 minutes, 30 cycles (denaturation at 94 °C x 1 min, hybridization at 60 °C x 1 min, extension at 68 °C x 2 min); and the final extension, at 68 °C x 7 min.
Electrophoresis
The amplification products of RT-PCR reactions were detected by electrophoresis (2% agarose gel, stained with 0.5 μg/mL ethidium bromide). The DNA migration was conducted in a horizontal electrophoresis chamber (Minicell * EC 370M), with a constant 60-100 V voltage and TAE buffer (40 mM Tris-acetic acid, 0.5 M EDTA, pH 8). DNA bands were visualized on a photodocumentation system, Gel Doc 1000 (Bio-Rad®), using the Multi-Analyst Program (Bio-Rad®); the size of the amplification bands was compared with a molecular marker of 1 kb DNA Ladder (Bioneer).
cDNA Purification and Cloning
The cDNA samples to be purified were fractionated on agarose gels and prepared in the TAE buffer, and the bands were visualized on the ultraviolet light transilluminator. The fraction of interest was cut with a sterile scalpel, and the DNA was eluted with the commercial Wizard®SV Gel and PCR Clean-Up System kit (Promega), according to the manufacturer’s instructions. For the ligation reaction of the purified PCR products and the commercial vector pGEM®-T Easy, 1 U of the T4 DNA ligase enzyme was used. A plasmid-insert molar ratio 1: 5 was used in [30 mM Tris-HCl (pH 7.8), ten mM MgCl2, ten mM DTT, one mM ATP] buffer and incubated at four °C for 12 h. After confirming the hsp83 into pGEM®-T Easy, this recombinant molecule was digested, and the purified cDNA was subcloned in the expression plasmid pQE30 (Qiagen, Valencia, CA).
cDNA Sequencing and Bioinformatic Analysis
Sequencing of cloned cDNA was achieved at the Sequencing Service of the National Center for Microbiology, Instituto de Salud Carlos III, Madrid, Spain, using 373A system, Model 377, Applied Biosystem, following the protocol BigDye Terminator Cycle Sequencing Ready Reaction Kit (ABI-PRISM, PE Biosystems). The cDNA sequence obtained was compared with the sequences included in the databases (nucleic acid and proteins) (GenBank, EMBL) and analyzed by bioinformatics programs. Multiple alignments were done using the Clustal W. The alignment was revised using BioEdit Sequence Alignment Editor, version 7.0.5.3 (Hall, 2013). Nucleotide and amino acid sequences were analyzed using the EditSeq program of DNAstar (Lasergene®, Madison, USA). Identity and similarities were looked for in the databases (protein and nucleic acid) (EMBL, GenBank) by BLAST (Boratyn et al., 2019). Other analyses were performed using CDD-Search of the NCBI (National Center for Biotechnology Information) (Lu et al., 2020), Interpro of the EBI (European Bioinformatics Institute) (Mitchell et al., 2019), Motif scan, and ExPASy (Expert Protein Analysis System) of the SIB (Swiss Institute of Bioinformatics) (Proteomics Server) (Artimo et al., 2012). Epitopes B prediction was performed using the Protean program from DNAstar (Lasergene®, Madison, USA). In addition, the program BcePred (Prediction of continuous B-cell epitope in antigenic sequences using physico-chemical properties) (Saha et al., 2005) was employed.
cDNA Expression
Recombinant His-fusion proteins were expressed in Escherichia coli M15 (Qiagen, Valencia, CA) by induction with isopropyl-β-D-thiogalactopyranoside (IPTG, one mM) (Sigma, Madrid, Spain). The cultures were lysed using sonication and examined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and the recombinant proteins were purified by affinity chromatography with nickel agarose as described by the manufacturer (Qiagen, Valencia, CA). The protein concentration was determined by the Bradford technique (Bradford, 1976). The purity of the proteins was analyzed by SDS polyacrylamide gel electrophoresis (SDS-PAGE), done under denaturing conditions according to Laemmli (1970), on 10% separation gel and 4.5% packing gel at 15 mA. After separating the sample, the gel was stained with Comassie Blue-Silver (De Moreno et al., 1985) for 45 min, then placed in the bleaching solution until the sharp bands were visible. The molecular weight (MW) of the expressed proteins was compared with the Broad Range Protein Molecular Weigh Markers (Promega).
Western Blot (Immunoblotting)
After the proteins were separated by electrophoresis, they were transferred to a nitrocellulose membrane, according to Towbin et al. (1979), using a wet system and the Mini Trans-Blot® Cell (BioRad) equipment. For this, two sponges and two filter papers the size of the gel, impregnated in transfer buffer and the nitrocellulose membrane, were placed, and the gel was placed on the membrane. Then, two filter papers and a sponge equally impregnated with the buffer were placed in the same way and transferred at 100 mA for one hour at four °C. After the transfer had elapsed, the nitrocellulose membrane was incubated in blocking solution (0.15M PBS pH 7.2, 7.5% skim milk) overnight at four °C. Subsequently, it was incubated with anti-Epi rabbit hyperimmune serum (against T. cruzi epimastigotes) in a 1/1500 dilution in reaction solution (0.15 M PBS pH 7.2, 150 mM NaCl, ten mM Tris-HCl pH 7.4) for 90 min at 37 °C. Then, eight washes of 5 min each were performed with the wash solution (0.15 M PBS pH 7.2, 0.1% Tween 20) and subsequently incubated with the conjugate, anti-rabbit IgG coupled to peroxidase at a dilution of 1/3000 in reaction solution (PBS 0.15 M pH 7.2, NaCl 150 mM, Tris-HCl 10 mM pH 7.4) and under the same conditions described above; then, the washes were repeated. The visualization of the immune complexes was carried out by chemiluminescence incubating the membrane with the mixture provided by the commercial Super Signal® West Pico kit (Pierce) for one minute in the dark and developing using a photographic plate according to the manufacturer’s instructions with an exposure time of three minutes.
By RT-PCR, the cDNA was obtained with an approximate size of 2000 bp, which corresponds approximately to the expected size of the hsp83 molecule of T. cruzi. Figure 1A shows the bands obtained by RT-PCR from the RNA samples extracted with the different commercial kits. For gene expression, M15 transformed cells were induced with IPTG, control aliquots were taken from uninduced cells and post-1, 2, and 3 hours of induction, and analyzed on Coomassie-Silver stained acrylamide gel. The expected size of the recombinant protein, according to the deduced amino acid sequence, is 80.7 kDa + 1 kDa of the histidine tail; so, a molecular weight of almost 83 kDa would be expected. Figure 1B shows the discrete expression of a protein of nearly 83 kDa in sample two at the second hour of incubation and slightly higher at the third hour. Additionally, bands that are not present in the uninduced control are observed and appear at the second and third hour of post-induction (Figure 1B); perhaps they are truncated forms of the protein, which the bacteria tend to eliminate.
The recombinant protein was purified by affinity chromatography from the cultures induced at 3 hours with IPTG. Each fraction obtained by chromatography and the sediment and the supernatant before and after passing through the Ni-NTA nickel column were analyzed using a polyacrylamide gel. For the identification of the recombinant protein, the fractions obtained from the purification of the recombinant protein HSP83 from T. cruzi were electrophoresis (SDS-PAGE 10%) and, subsequently, electroblotting on a nitrocellulose membrane. The nitrocellulose membrane was incubated with hyperimmune rabbit serum anti-Epimastigotes of T. cruzi; after development with luminol, a specific reaction with the recombinant protein is evidenced (Figure 1C).
Figure 1. Cloning and expression of HSP83 of T. cruzi. A) Bands obtained by RT-PCR from the RNA samples. B) Transformed M15 cells induced with IPTG (1, 5), control aliquots taken from uninduced cells and post-1, 2, and 3 hours of induction (2, 3, 4 and 6, 7, 8, respectively), and analyzed on Coomassie-Silver stained acrylamide gel. C) Immunoblotting using hyperimmune rabbit serum anti-Epimastigotes of T. cruzi to detect purified recombinant protein.
The complete sequence of hsp83 cDNA was a fragment of 2115 bp corresponding to the open reading frame (ORF). This frame coded for a peptide of 704 amino acids, with a molecular mass of nearly 83 kDa and an isoelectric point of 5.0. The deduced amino acid sequence showed an amidation site (10-13), four possible N-glycosylation sites (48-51, 116-119, 151-154, 452-455), three cAMP and cGMP-dependent protein kinase phosphorylation sites (246-249, 273-276, 329-332), 11 Casein kinase II phosphorylation sites (104-107, 155-158, 164-167, 291-294, 296-299, 306-309, 360-363, 525-528, 531-534, 609-612, 621-624), four N-myristoylation sites (160-164, 190-195, 232-237,442-447), nine Protein kinase C phosphorylation sites (104-106, 291-293, 332-334, 380-382, 410-412, 531-533, 553-555, 644-646, 685-687), and four Tyrosine kinase phosphorylation sites (57-65, 218-225, 273-281, 616-623). A motif and several domains were found, such as a Chaperone protein htpG [htpG] motif (80-704), a Heat shock hsp90 proteins family signature (103-112), a TonB-dependent receptor proteins signature 1 (1-40), two Bipartite nuclear localization signals (400-414, 554-571), a Glutamic acid-rich region profile (257-327), a HATPase_cHistidine kinase-, DNA gyrase B-, a HSP90-like ATPase domain (105-258), and a HSP90 protein domain (261-704) (InterPro IPR001404 pFam PF00183). With respect to the possible immunogenicity of the complete sequence of hsp83 cDNA, five B epitopes were predicted in the molecule (226TKNEE230, 244EEGEKKKKT252, 327FEPSKK332, 470TGDSKK475, and 527EETEEEKKQR536) (Figure 2).
Figure 2. Schematic representation of the cloned cDNAs analyzed. The bars below deduced proteins represent the domains predicted in the sequence. The short grey lines below numbers represent the positions of the B epitopes; bp stands for base par, and aa means amino acid.
This sequence showed high identity with heat shock protein 85 [Trypanosoma cruzi Dm28c, Venezuela] (ESS64949.1), with Heat shock-like 85 kDa protein [Trypanosoma cruzi Dracon, Perú] (P06660.1) and lower identity with heat shock protein 85, putative [Trypanosoma cruzi CL Brener El Sayed, Brazil] (XP_811791.1) (Figure 3).
Figure 3. Clustal Alignment with heat shock protein 85 [Trypanosoma cruzi Dm28c, Venezuela] (ESS64949.1), Heat shock-like 85 kDa protein [Trypanosoma cruzi Dracon, Perú] (P06660.1), and heat shock protein 85, putative [Trypanosoma cruzi CL Brener El Sayed, Brazil] (XP_811791.1).
The cloning of genes and their expression as recombinant proteins facilitates the study of both the genes and the proteins they encode, as well as their use as tools in tests that help to elucidate the cellular mechanisms in which they are involved. In the present study, it was possible to clone and express the gene encoding the heat shock protein HSP83 from T. cruzi. To date, there are merely two old reports about cloning this molecule (Dragon et al., 1987; Nadeau et al., 1992), and only its genomic structure and ATPase activity have been studied. In addition, cloning and expression of the gene encoding HSP83 has been obtained in other tripanosomatids, such as Trypanosoma brucei (Mottram et al., 1989), Hsp83 from Leishmania mexicana amazonensis (Shapira & Pinelli, 1989), Hsp83 of Leishmania infantum (Angel et al., 1996), and Hsp83 from Leishmania braziliensis (Silva et al., 2013; Menezes-Souza et al., 2014).
Since HSP83 has been reported to be a protein highly expressed by T. cruzi, even under basal conditions (Dragon et al., 1987), and given its importance in the mechanisms of cellular metabolism of the parasite and the need to expand studies in this area, it was decided to address the cloning strategy using the mRNA as a starting sample and the amplification of the coding region by means of the RT-PCR technique. Thanks to the information from the T. cruzi Genome project, this sequence was made available to design the cloning strategy.
In all cases, a band with a size of 2139 bp was amplified, the expected size according to the GenBank sequence. The verification of the cloning by colony PCR, using the primers for the coding region of HSP83, helped amplification of a band of 2139 bp in eight of the 19 colonies analyzed. The lack of expression in a high percentage of colonies analyzed could be due to the conditions of the cells used in cloning (which were not as competent at accepting the plasmid with the gene).
Likewise, the cDNA was subcloned into different expression vectors to facilitate its expression. However, despite successful verification of the insert’s presence in the recent subcloning, the HSP83 sequence or the entire plasmid disappeared in some cases when trying to express it. This phenomenon suggests that the HSP83-plasmid construct may be unstable within these complex plasmids, and bacteria may degrade the plasmid (with the loss of the HSP83 sequence) or even eliminate it (Maizels et al., 1991).
Subsequently, it was decided to carry out the subcloning process in the vector pQE30, which was simpler and used for the cloning and expression of HSP83 from Leishmania infantum (Angel et al., 1996). However, the same result occurred: subcloning was achieved but followed by loss. This reinforces the recurrence of the hypothesis that the recombinant plasmid is expelled from the bacterium, perhaps because the HSP83 sequence, being evolutionarily conserved, can interfere with the cell’s metabolism, and the cell tends to eliminate it (Huynh et al., 1985). On the other hand, sometimes fragments of genetic material do not always integrate precisely; this may lead to artifacts that the sera would not recognize.
Considering the high similarity reported between Hsp83 sequences from T. cruzi and those from T. brucei and Hsp83 from L. m. amazonensis—being 95% and 85%, respectively (Rondinelli, 1994)—, along with the reported expression of L. infantum Hsp83 in cells transformed with pQE30 (Angel et al., 1996), and thinking that the problem may be the residence time in the cell, a decision was made to clone into the latter vector directly, without waiting for the processes to verify the presence of the insert, that is, a “blind” strategy. Thus, the residence time in the bacterium would be reduced, which could prevent expulsion before the respective expression. In this case, it was possible to obtain positive results, with the appearance of a nearly 83 kDa protein in chromatography-eluted fractions from lysates of competent cells transformed E. coli M15.
Reviewing published articles on cloning genes encoding Hsp90 heat shock proteins in different parasitic models, either from genomic DNA or mRNA, reveals that the expression of proteins from the Hsp90 family of Trypanosomatids was possible through cloning in expression vectors. Examples of these vectors are the Hsp83 from L. braziliensis with the vector pET28a-TEV (Menezes-Souza et al., 2014), Hsp90 from T. evansi with the vector pRSETA (Pallavi et al., 2010), Hsp83 from Leishmania braziliensis from pET28a (Silva et al., 2013), and Hsp83 from Leishmania infantum with pQE30 (Angel et al., 1996). Among the expression vectors used in these studies, only pRSETA was more complex and employed in the current investigation. However, the results obtained were unsatisfactory, perhaps because the time lapses between subcloning in the expression vectors and expression were several weeks, and bacteria might have shed the recombinant plasmid due to its instability during that time. Using a direct “blind” subcloning strategy may have yielded better results in these vectors. Alternatively, this could be a particular characteristic of the Hsp83 of T. cruzi since it is strange that such an important molecule from a parasite causing serious public health problems has not had more cloning and expression studies, possibly due to the difficulty of its expression in prokaryotic systems.
In this sense, it is striking that during the review of the cloning and expression procedure of Plasmodium falciparum Hsp90, it was verified that the researchers directly cloned into the expression vector pET-23a and obtained the recombinant protein by using competent cells E. coli BL21 (Kumar et al., 2003). Although the authors do not mention this aspect, they might have encountered circumstances similar to those discussed in this paper.
Regarding the expression of the recombinant protein, a low expression level and smaller bands were observed, which could be truncated forms of the protein, possibly due to bacterial attempts to degrade the recombinant protein. This degradation might occur due to sequence similarities, interfering with its metabolism (Huynh et al., 1985). This could explain the limited observation of the recombinant protein.
The expected size of the recombinant protein according to the deduced amino acid sequence is 80.7 kDa + 1 kDa for the histidine tail; so, a molecular weight of nearly 82 kDa would be expected, which was observed in the results. It is important to note that HSP83 or HSP85 was so named since the native protein has a weight of 83-85 kDa. The small difference in size of nearly three kDa may be due to glycosylation of the native protein, as the prediction analysis of possible post-translational modifications of the sequence by Motif Scan of ISB-ISREC (Pagni et al., 2001) evidenced four potential N-glycosylation sites. One of these sites can be considered functional since it is estimated that an N-glycosylation contributes an additional three kDa to the protein’s molecular weight (Obregon-Henao et al., 2001).
As expected, problems were also encountered during the protein purification process because it was a large molecule of nearly 83 kDa, and its same weight was difficult to bind to the resin of the column only by a six histidine tail. It may be convenient to use another purification protocol or subclone the cDNA into a vector capable of binding to a larger fusion molecule that serves as a support and can then be eliminated, such as GST or MBP.
Considering the poor expression of hsp 83 and aiming to confirm that the band of nearly 83 kDa, whose expression increased with induction time, corresponded to the HSP83 protein of the parasite, a rabbit’s hyperimmune serum against T. cruzi epimastigotes (from the same isolate used in the present work) was used. This serum has been employed in several experimental protocols, with recognition of antigens with approximate molecular weights of 85, 80, 65, 62, 52, 38, 35, 34, 30, and 27 kDa in epimastigotes of the EP strain (Contreras et al., 1998; Graterol et al., 2013).
Immunoblotting analysis of the fractions obtained by nickel Ni-NTA column chromatography using hyperimmune anti-epimastigote serum facilitated the identification of a single band with an approximate molecular weight of 83 kDa, indicating, in this case, a reaction with the recombinant protein obtained in this study. Due to the low concentration of the recombinant protein used in the immunoblotting technique, the researchers decided to use a highly sensitive development system such as luminography, in particular with the Super Signal® West pico system, enabling the detection of immune complexes in pictograms. This explains the strong signal of the band observed on the X-ray film used to document the reaction.
Heat shock proteins have been shown to be immunogenic in various parasitic models (Polla, 1991), and most are important antigens recognized by patient sera. In addition, it has been reported that the immune response is directed to the less conserved areas of the protein, ensuring the specificity of the response (Skeiky et al., 1995; Angel et al., 1996). For example, specific anti-HSP90 antibodies have been documented in schistosomiasis (Johnson et al., 1989), leishmaniasis (de Andrade et al., 1992), and malaria (Zhang et al., 2001). Immunodominant HSP70 antigens have also been reported in patients with hydatidosis (Colebrook & Lightowlers, 1997), schistosomiasis (Moser et al., 1990), onchocerciasis (Rothstein et al., 1989), leishmaniasis (de Andrade et al., 1992; Wallace et al., 1992; Skeiky et al., 1995; Quijada et al., 1996), Chagas disease (Engman et al., 1989; Requena et al., 1993) and malaria (Kumar et al., 2003), and HSP60 in the case of leishmaniasis (Rey-Ladino et al., 1997). Regarding low molecular weight HSPs, their diagnostic value in schistosomiasis and cysticercosis has been studied and confirmed (Nene et al., 1986; Ferrer et al., 2005).
It is important to note that no reactivity of the T. cruzi anti-epimastigote hyperimmune serum was observed with any other band of any size, only the 83 kDa, which would indicate that the possible truncated molecules lack the epitopes recognized by these antibodies. Using the complete protein to immunize rabbits and obtain hyperimmune sera that allow epitope mapping assays would be interesting.
Even at low concentrations, this recombinant protein HSP83 from T. cruzi can be used to elaborate a monospecific polyclonal antibody, which constitutes a very valuable tool to study the participation of this stress protein in the parasite’s differentiation processes, such as epimastigogenesis, metacyclogenesis, amastigogenesis, and trypomastigogenesis.
Despite the difficulties in the expression and purification of the recombinant protein, it could be expressed and identified with the hyperimmune serum. Therefore, this cloning and expression strategy is recommended for molecules with difficult cloning and expression.
This study has been supported by Universidad de Carabobo Projects: LOCTI-001.
The authors declare no competing interests.
All the authors declare that the final version of this paper was read and approved.
The total contribution percentage for the conceptualization, preparation, and correction of this paper was as follows: A.F. 33.3 %., A.R.D.L. 33.3 % and E.F. 33.3 %.
Data sharing is not applicable since no new data was created or analyzed in this study.
A Preprint version of this paper was deposited in: https://doi.org/10.5281/zenodo.8252927
Alarcón de Noya, B.A., Díaz-Bello, Z., Colmenares, C., Ruiz-Guevara, R., Murillo, L., Muñoz-Calderón, A., Noya, O. (2015). Update on oral Chagas disease outbreaks in Venezuela: epidemiological, clinical and diagnostic approaches. Memórias do Instituto Oswaldo Cruz 110, (3) 77-386. https://doi: 10.1590/0074-02760140285
Angel, S.O., Requena, J.M., Soto, M., Criado, D., Alonso, C. (1996). During canine leishmaniasis a protein belonging to the 83-kDa heat-shock protein family elicits a strong humoral response. Acta Tropica, 62(1), 45-56. https://doi.org/10.1016/S0001-706X(96)00020-4
Añez, N., Crisante, G., Rojas, A., Segninib, S., Espinoza-Álvarez, O., Teixeira, M.M.G. (2020). Update on Chagas disease in Venezuela during the period 2003–2018. Acta Tropica 203, 105310. https://doi.org/10.1016/j.actatropica.2019.105310.
Artimo, P., Jonnalagedda, M., Arnold, K., Baratin, D., Csardi, G., de Castro, E., Duvaud, S., Flegel, V., Fortier, A., Gasteiger, E., Grosdidier, A., Hernandez, C., Ioannidis, V., Kuznetsov, D., Liechti, R., Moretti, S., Mostaguir, K., Redaschi, N., Rossier, G. Xenarios, I. Stockinger, H. (2012). ExPASy: SIB bioinformatics resource portal, Nucleic Acids Research. 40 (1), W597-W603, https://dx.doi: 10.1093/nar/gks400.
Boratyn, G.M., Thierry-Mieg, J., Thierry-Mieg, D., Busby, B., Madden, T.L. (2019). Magic-BLAST, an accurate RNA-seq aligner for long and short reads, BMC Bioinformatics. 20, 405. https://dx.doi: 10.1186/s12859-019-2996-x.
Bradford, M.M. (1976). A rapid and sensivity method for quantitation of micrograms quantities of protein utilizing the principle of protein dye binding. Analytical of Biochemistry, 72 (1-2), 248-254. https://doi.org/10.1016/0003-2697(76)90527-3
Colebrook, A. Lightowlers, M., (1997). Serological reactivity to heat shock protein 70 in patients with hydatid disease. Parasite Immunology, 19, 41-46. https://doi.org/10.1046/j.1365-3024.1997.d01-141.x
Contreras, V.T., Araque, W. Delgado, V. (1994). Trypanosoma cruzi: Metacyclogenesis in vitro. I. Changes in the properties of Metacyclic Tripomastigotes maintained in the laboratory by different methods. Memórias do Instituto Oswaldo Cruz, 89 (2), 253-359. https://doi.org/10.1590/S0074-02761994000200026
Contreras, V.T., De Lima, A.R, Zorrilla, G. (1998). Trypanosoma cruzi: Maintenance in Culture Modify Gene and Antigenic Expression of Metacyclic Trypomastigotes. Memórias do Instituto Oswaldo Cruz. 93 (6),753-60. https://doi.org/10.1590/S0074-02761998000600011
De Andrade, C.R., Kirchhoff, L.V., Donelson, J.E., Otsu, K. (1992). Recombinant Leishmania Hsp90 and Hsp70 are recognized by sera from visceral leishmaniasis patients but not Chagas’ disease patients. Journal of Clinical Microbiology, 30(2), 330-335. https://doi.org/10.1128/jcm.30.2.330-335.1992
De Moreno, M.R., Smith, J.F., Smith, R.V. (1985). Silver staining of proteins in polyacrilamide gels: increased sensitivy through a combined blue-silver stain procedure. Analytical of Biochemistry, 151 (2), 466-470. https://doi.org/10.1016/0003-2697(85)90206-4
Dragon, E.A., Sias, S.A., Kato, E.A., Gabe, J.D. (1987). The genome of Trypanosoma cruzi contains a constitutively expressed arranged multicopy gene homologous to major heart shock protein. Molecular Cell Biology, 7, 1271-1275.
El-Sayed, N.M., Myler, P.J., Bartholomeu, D.C., Nilsson, D., Aggarwal, G., Tran, A.N., Ghedin, E., Worthey, E.A., Delcher, A.L., Blandin, G., Westenberger, S.J., Caler, E., Cerqueira, G.C., Branche, C., Haas, B., Anupama, A., Arner, E., Aslund, L., Attipoe, P., Bontempi, E., Bringaud, F., Burton, P., Cadag, E., Campbell, D.A., Carrington, M., Crabtree, J., Darban, H., da Silveira, J.F., de Jong, P., Edwards, K., Englund, P.T., Fazelina, G., Feldblyum, T., Ferella, M., Frasch, A.C., Gull, K., Horn, D., Hou, L., Huang, Y., Kindlund, E., Klingbeil, M., Kluge, S., Koo, H., Lacerda, D., Levin, M.J., Lorenzi, H., Louie, T., Machado, C.R., McCulloch, R., McKenna, A., Mizuno, Y., Mottram, J.C., Nelson, S., Ochaya, S., Osoegawa, K., Pai, G., Parsons, M., Pentony, M., Pettersson, U., Pop, M., Ramirez, J.L., Rinta, J., Robertson, L., Salzberg, S.L., Sanchez, D.O., Seyler, A., Sharma, R., Shetty, J., Simpson, A.J., Sisk, E., Tammi, M.T., Tarleton, R., Teixeira, S., Van Aken, S., Vogt, C., Ward, P.N., Wickstead, B., Wortman, J., White, O., Fraser, C.M., Stuart, K.D. y Andersson, B. (2005). The genome sequence of Trypanosoma cruzi, etiologic agent of Chagas disease. Science; 309(5733), 409-415. https://doi.org/10.1126/science.1112631
Engman, D.M., Kirchhoff, L.V. Donelson, J.E. (1989). Molecular cloning of Hsp70 mitochondrial member of the Hsp70 family. Molecular Cell Biology, 9 (11), 5163-5168. https://doi.org/10.1128/mcb.9.11.5163-5168.1989
Ferrer, E., González, L., Foster-Cuevas, M., Cortéz, M., Dávila, I., Rodríguez, M., Sciutto, E., Harrison, L.J.S., Parkhouse, R.M.E., y Gárate, T. (2005). Taenia solium: characterization of a small heat shock protein (Tsol-sHSP35.6) and its possible relevance to the diagnosis and pathogenesis of neurocysticercosis. Experimental Parasitology, 110 (1), 1-11. https://doi.org/10.1016/j.exppara.2004.11.014
Folgueira, C. y Requena, J.M. (2007). A post genomic view of the heat shock proteins in kinetoplastids. FEMS. Microbiology Reviews, 31 (4), 359-377. https://doi.org/10.1111/j.1574-6976.2007.00069.x
Graefe, S.E.B., Wiesgigl, M., Gaworski, I., Macdonald, A. Clos J. (2002). Inhibition of Hsp90 in Trypanosoma cruzi induces a strees response but no stage differentiation. Eucaryotic Cell, 11, 936-943. https://doi.org/10.1128/EC.1.6.936-943.2002
Graterol, D., Ramírez, D., Ramos, A., Arteaga, R., Mundaray, O., Pineda, W., Navarro, M., Domínguez, M., De Lima, A. y Contreras, V. (2013). Epimastigogénesis de Trypanosoma cruzi Dm28c en medio ML15HA: análisis proteico, glicoproteico y antigénico. Salus. 17, 12-18.
Hall, T. (2013). BioEdit (v7.1.11). [http://www.mbio.ncsu.edu/BioEdit/bioedit.html]
Huynh, T.V., Young, R., Davis, R.W. (1985). In: DNA cloning: a practical approach (Glover, D.M., de.). IRL Press. Oxford, 1, 49-78.
Jhonson, J.L. (2012). Evolution and function of diverse Hsp90 homologs and cochaperone proteins. Biochimica et Biophysica Acta, 1823, 607–613. https://doi.org/10.1016/j.bbamcr.2011.09.020
Johnson, K.S., Wells, K., Bock, J.V., Nene, V., Taylor, D.W., Cordingley, J.S. (1989). The 86-kilodalton antigen from Schistosoma mansoni is a heat-shock protein homologous to yeast HSP-90. Molecular and Biochemical Parasitology, 36(1), 19-28. https://doi.org/10.1016/0166-6851(89)90196-5
Kumar, R., Musiyenko, A., Barik, S. (2003). The heat shoch protein Hsp90 of Plasmodium falciparum and antimalarial activity of its inhibitor, geldanamycin. Malaria Journal, 2, 30-40. https://doi.org/10.1186/1475-2875-2-30
Laemmli, U.K. (1970). Cleavage of structurals Proteins during the assembly of the head of bacteriophago T4. Nature, 227, 680-685. https://doi.org/10.1038/227680a0
Lu, S., Wang, J., Chitsaz, F., Derbyshire, M.K., Geer, R.C., Gonzales, N.R., Gwadz, M., Hurwitz, D.I., Marchler, G.H., Song, J.S., Thanki, N., Yamashita, R.A., Yang, M., Zheng, D.Z.C., Lanczycki, C.J., Marchler-Bauer, A., (2020). CDD/SPARCLE: The Conserved Domain Database in 2020, Nucleic Acids Research. 48, D265-D268. https://doi.org/10.1093/nar/gkz991
Maizels, R.M., Blaxter, M.L., Robertson, B.D., Selkirk, M.E. (1991). Parasite antigens. Parasite genes: a laboratory manual for Molecular parasitology. Cambridge, United States of Amerca: University Press.
Menezes-Souza, D., de Oliveira, T., de Souza, M., Reis-Cunha, J., Pinto, R., Martins, C., Ferraz, E., da Cunha, L., Toshio, R. Castanheira, D. (2014). Epitope mapping of the Hsp83.1 protein of Leishmania braziliensis discloses novel targets for immunodiagnosis of tegumentary and visceral clinical forms of leishmaniasis. Clinical and Vaccine Immunology. 21(7), 949-959. https://doi.org/10.1128/CVI.00151-14
Mitchell, A.L., Attwood, T.K., Babbitt, P.C., Blum, M., Bork, P., Bridge, A., Brown S.D., Chang, H.Y., El-Gebali, S., Fraser, M.I., Gough, J., Haft, D.R., Huang, H., Letunic, I., Lopez, R., Luciani, A., Madeira, F., Marchler-Bauer, A., Mi, H., Natale, D.A., Necci, M., Nuka, G., Orengo, C., Pandurangan, A.P., Paysan-Lafosse, T., Pesseat, S., Potter, S.C., Qureshi, M.A., Rawlings, N.D., Redaschi, N., Richardson, L.J., Rivoire, C., Salazar, G.A., Sangrador-Vegas, A., Sigrist, C.J.A., Sillitoe, I., Sutton, G.G., Thanki, N., Thomas, P.D., Tosatto, S.C.E., Yong, S.Y., Finn, R.D. (2019). InterPro in 2019: improving coverage, classification and access to protein sequence annotations, Nucleic Acids Research. 47, D351–D360. https://dx.doi.org/10.1093/nar/gky1100
Moser, D., Doumbo, O., Klinkrt, M., (1990). The humoral response to heat shock protein 70 in human and murine Schistosomiasis mansoni. Parasite Immunology, 12, 341-352. https://doi.org/10.1111/j.1365-3024.1990.tb00973.x
Mottram, J.C., Murphy, W.J., Agabian, N. (1989). A transcriptional analysis of the Trypanosoma brucei Hsp83 gene cluster. Molecular and Biochemical Parasitology, 37, 115-128. https://doi.org/10.1016/0166-6851(89)90108-4
Nadeau, K., Sullivan, M.A., Bradley, M., Engman, D.M., Walsh, C.T. (1992). 83-kilodalton heat shock proteins of trypanosomes are potent peptide-stimulated ATPases. Protein Science, 1(8), 970-979. https://doi.org/10.1002/pro.5560010802
Nene, V., Dunne, D.W., Johnson, K.S., Taylor, D.W. y Cordingley, J.S. (1986). Sequence and expression of a major egg antigen from Schistosoma mansoni. Homologies to heat shock proteins and alpha-crystallins. Molecular and Biochemical Parasitology, 21(2), 179-188. https://doi.org/10.1016/0166-6851(86)90021-6
Obregon-Henao, A., Gil, D.L., Gómez, D.I., Sanzón, F., Teale, J.M., Restrepo, B. (2001). The role of N-linked carbohydrates in the antigenicity of Taenia solium metacestode glycoproteins of 12, 16 and 18 kD. Molecular and Biochemical Parasitology 114 (2), 209-215. https://doi.org/10.1016/S0166-6851(01)00256-0
Pagni, M., Iseli, C., Junier, T., Falquet, L., Jongeneel, V., Bucher, P. (2001). TrEST, trGEN and Hits: Access to databases of predicted protein sequences. Nucleic Acids Research 29 (1), 148-151. https://doi.org/10.1093/nar/29.1.148
Pallavi, R., Roy, N., Nageshan, R.K., Talukdar, P., Pavithra, S.R., Reddy, R., Venketesh, S., Kumar, R., Gupta, A.K., Singh, R.K., Yadav, S.C., Tatu, U. (2010). Heat shock protein 90 as a drug target against protozoan infections: biochemical characterization of HSP90 from Plasmodium falciparum and Trypanosoma evansi and evaluation of its inhibitor as a candidate drug. Journal of Biological Chemistry, 285(49), 37964-37975. https://doi.org/10.1074/jbc.M110.155317
Polla, B. (1991). Heat shock proteins in host-parasite interactions. Immunology Today, 12(3), A38-A41. https://doi.org/10.1016/S0167-5699(05)80011-8
Quijada, L., Requena, J.M., Soto, M., Gomez, L.C., Guzmán, F., Patarroyo, M.E., Alonso, C. (1996). Mapping of the linear antigenic determinants of the Leishmania infantum hsp70 recognized by leishmaniasis sera. Immunology Letters, 52(2-3), 73-79. https://doi.org/10.1016/0165-2478(96)02585-0
Requena, J.M., Soto, M., Guzmán, F., Maekelt, A., Noya, O., Patarroyo, M.E., Alonso, C. (1993). Mapping of antigenic determinants of the T. cruzi hsp70 in chagasic and healthy individuals. Molecular Immunology, 30(12), 1115-1121. https://doi.org/10.1016/0161-5890(93)90158-8
Rey-Ladino, J.A., Joshi, P.B., Singh, B., Gupta, R., Reiner, N.E. (1997). Leishmania major: molecular cloning, sequencing, and expression of the heat shock protein 60 gene reveals unique carboxy terminal peptide sequences. Experimental Parasitology, 85(3), 249-263. https://doi.org/10.1006/expr.1996.4137
Rondinelli, E. (1994). Conservation of heat-shock proteins in Trypanosoma cruzi. Parasitology Today, 10(5), 172-176. https://doi.org/10.1016/0169-4758(94)90020-5
Rothstein, N., Yates, J., Higashi, G. Rajan, T. (1989). Onchocerca volvulus heat shock protein 70 is a major immunogen in amicrofilaremic individuals from a filariasis-endemic area. Molecular and Biochemical Parasitology, 33(3), 229-236. https://doi.org/10.1016/0166-6851(89)90084-4
Saha, S., Bhasin, M., Raghava, G.P.S. (2005). Bcipep: A database of B-cell epitopes, BMC Genomics, 6, 79. https://doi.org/10.1186/1471-2164-6-79
Shapira, M., Pinelli, E. (1989). Heat-shock protein 83 of Leishmania mexicana amazonensis is an abundant cytoplasmic protein with a tandemly repeated genomic arrangement. European Journal of Biochemestry, 185 (2), 231-236. https://doi.org/10.1111/j.1432-1033.1989.tb15107.x
Skeiky, Y.A., Benson, D., Guderian, J.A., Whittle, J.A., Bacelar, O., Carvalho, E.M., y Reed, S.G. (1995). Inmune responses of leishmaniasis patients to heat shock proteins of Leishmania species and humans. Infection and Immunity, 63 (10), 4105-4114. https://doi.org/10.1128/iai.63.10.4105-4114.1995
Silva, K.P., Seraphim, T.V. Borges, J.C. (2013). Structural and functional studies of Leishmania braziliensis Hsp90. Biochimica et Biophysica Acta, 1834(1), 351-361. https://doi.org/10.1016/j.bbapap.2012.08.004
Towbin, H., Staehelin, T., Gordon, J. (1979). Electrophoretic transfer of proteinsfrom polyacrylamide gels tonitrocelulose sheets: procedure and some applications. Proceedings of the National Academy of Sciences, 76(9), 4350-4354. https://doi.org/10.1073/pnas.76.9.4350
Urményi TP, Silva R, Rondinelli E. (2014). The heat shock proteins of Trypanosoma cruzi. Subcellular Biochemistry, 74, 119-135. https://doi.org/10.1007/978-94-007-7305-9_5
Wallace, G.R., Ball, A.E., MacFarlane, J., el Safi, S.H., Miles, M.A., Kelly, J.M. (1992). Mapping of a visceral leishmaniasis-specific immunodominant B-cell epitope of Leishmania donovani Hsp70. Infection and Immunity, 60(7), 2688-2693. https://doi.org/10.1128/iai.60.7.2688-2693.1992
WHO. (2023). Chagas Disease (American trypanosomiasis) Fact sheet. World Health Organization, Geneva. http://www.who.int/news-room/fact-sheets/detail/chagas-disease-(american-trypanosomiasis).
Zhang, M., Hisaeda, H., Kano, S., Matsumoto, Y., Hao, Y.P., Looaresuwan, S, Aikawa, M., Himeno, K. (2001). Antibodies specific for heat shock proteins in human and murine malaria. Microbes and Infection, 3(5), 363-367. https://doi.org/10.1016/S1286-4579(01)01391-0
1* Corresponding author
Ana Fernández, anakfdz@gmail.com, https://orcid.org/0009-0006-9466-2401
Ana Rita De Lima, delima.anarita@gmail.com, https://orcid.org/0000-0002-6817-9495
Elizabeth Ferrer, elizabeth.ferrer@gmail.com, https://orcid.org/0000-0002-4173-6642
Laboratorio de Protozoología, Instituto BioMolP, Facultad de Ciencias de la Salud, Universidad de Carabobo, Valencia, Venezuela.
2 Instituto de Investigaciones Biomédicas “Dr. Francisco J. Triana Alonso” (BIOMED). Facultad de Ciencias de la Salud Sede Aragua, Universidad de Carabobo, Maracay, estado Aragua, Venezuela.
Cloning and Characterization of the Gene Encoding the Heat Shock Protein HSP83 from Trypanosoma cruzi (Ana Fernández • Ana Rita De Lima • Elizabeth Ferrer) Uniciencia is protected by Attribution-NonCommercial-NoDerivs 3.0 Unported (CC BY-NC-ND 3.0)
URL: www.revistas.una.ac.cr/uniciencia
Correo electrónico: revistauniciencia@una.cr