Acessibilidade / Reportar erro

Genes highly overexpressed in salt-stressed Young oil palm (Elaeis guineensis) plants 1 1 Research developed at Brasília, DF, Brazil

Genes altamente superexpressos em dendezeiros jovens (Elaeis guineensis) estressados por sal

ABSTRACT

RNA-seq is a technique based on the large-scale sequencing of transcript-derived cDNAs using next-generation sequencing platforms mostly used today to characterize an organism’s transcriptome. The analysis of RNA-seq data allows for identifying genes differentially expressed in a given condition, such as salt stress. This study aimed to search and characterize genes from the African oil palm (Elaeis guineensis Jacq.) highly up-regulated during salt stress, with a long-term goal of gene promoter prospection and validation. The apical leaves from the control (electrical conductivity of ~2 dS m-1) and salt-stressed (~40 dS m-1) young oil palm plants, collected at 5 and 12 days after the beginning of the stress, were subjected to extraction of total RNA, with three plants (replicates) per treatment. The complete genome ofE. guineensis, available at the National Center for Biotechnology Information, was used as the reference genome - BioProject PRJNA192219. The differential expression analysis led to the selection for further characterization of seven genes, which had increased expressions of 37-84 times under salt stress. The strategy used in this study enabled the selection of seven salt-responsive genes highly up-regulated during salt stress, and some of them coded for proteins already reported as responsible for salinity tolerance in other plant species through over-expression or knockout.

Key words:
transcriptomics; RNAseq; differential expression analysis; salt-responsive genes

RESUMO

O RNA-seq é uma técnica baseada no sequenciamento em larga escala de cDNAs derivados de transcritos - por meio de plataformas de sequenciamento de nova geração - amplamente utilizada atualmente para caracterizar o transcriptoma de um organismo. A análise dos dados de RNA-seq permite identificar genes diferencialmente expressos em uma determinada condição, como o estresse salino. Este estudo teve como objetivo prospectar genes do dendê (Elaeis guineensis) responsivos ao estresse salino e realizar sua caracterização funcional e estrutural. A folha apical de plantas jovens de dendê controle (condutividade elétrica de ~2 dS m-1) e estressadas (~40 dS m-1), coletadas aos 5 e 12 dias após o início do estresse, foi submetida à extração de RNA total, com três plantas (repetições) por tratamento. O genoma completo de E. guineensis, disponível no “National Center for Biotechnology Information”, foi utilizado como genoma de referência - BioProject PRJNA192219. A análise da expressão diferencial levou à seleção de sete genes, cujo nível de expressão aumentou entre 37 e 84 vezes sob estresse salino, para posterior caracterização. A estratégia utilizada neste estudo permitiu a seleção de sete genes responsivos ao sal suprarregulados durante o estresse, e alguns deles codificam proteínas já relatadas como responsáveis pela tolerância à salinidade em outras espécies de plantas por superexpressão ou knock-out.

Palavras-chave:
transcritômica; RNAseq; análise de expressão diferencial; genes responsivos ao sal

HIGHLIGHTS:

The applied prospection/annotation strategy led to the selection of seven salt-responsive genes.

Three genes are candidates for the prospection/validation of salt-responsive promoter sequences.

Five genes are candidates for the generation of oil palm plants tolerant to salt stress.

Introduction

Soil salinity is a problem present in more than 100 countries spread over almost all continents; approximately 20% of the agricultural land worldwide has saline and/or sodic soils, and between 25 and 30% of the irrigated land area is affected by salt and essentially commercially unproductive (Shahid et al., 2018Shahid, S. A.; Zaman, M.; Heng, L. Soil salinity: Historical perspectives and a world overview of the problem. In: Zaman, M.; Shahid, S. A.; Heng, L. (eds.). Guideline for salinity assessment, mitigation and adaptation using nuclear and related techniques. Cham: Springer , 2018. Chapter 2, p.43-53. https://doi.org/10.1007/978-3-319-96190-3_2
https://doi.org/10.1007/978-3-319-96190-...
).

Salinity is often seen as a problem for the agricultural sector, prompting actions aimed at prevention or remediation in the affected areas; however, in the scope of biosaline agriculture, this problem is seen as an opportunity for the production of food, fibers, and bioenergy, as well as for the recovery of degraded areas and the use of marginal ones (Joshi et al., 2020Joshi, A.; Kanthaliya, B.; Arora, J. Halophytes: The nonconventional crops as source of biofuel production. In: Grigore, M. N. (eds.). Handbook of halophytes. Cham: Springer, v.1, p.1-28, 2020. https://doi.org/10.1007/978-3-030-17854-3_126-1
https://doi.org/10.1007/978-3-030-17854-...
; Tıpırdamaz et al., 2020Tıpırdamaz, R.; Karakas, S.; Dikilitas, M. Halophytes and the future of agriculture. In: Grigore, M. N. (eds.). Handbook of halophytes - From molecules to ecosystems towards biosaline agriculture. Cham: Springer , 2020. Chapter 47, p.1-15. https://doi.org/10.1007/978-3-030-17854-3_91-1
https://doi.org/10.1007/978-3-030-17854-...
; Duarte & Caçador, 2021Duarte, B.; Caçador, I. Iberian halophytes as agroecological solutions for degraded lands and biosaline agriculture. Sustainability, v.13, p.1-14, 2021. https://doi.org/10.3390/su13021005
https://doi.org/10.3390/su13021005...
). In this context, cultivation systems for saline environments are developed, using the ability of some plants to grow under saline conditions in combination with the use of saline soils and water resources, and better soil and water management (Ventura et al., 2015Ventura, Y.; Eshel, A.; Pasternak, D.; Sagi, M. The development of halophyte-based agriculture: Past and present. Annals of Botany, v.115, p.529-540, 2015. https://doi.org/10.1093/aob/mcu173
https://doi.org/10.1093/aob/mcu173...
; Duarte & Caçador, 2021Duarte, B.; Caçador, I. Iberian halophytes as agroecological solutions for degraded lands and biosaline agriculture. Sustainability, v.13, p.1-14, 2021. https://doi.org/10.3390/su13021005
https://doi.org/10.3390/su13021005...
).

The transcriptome is the complete set of RNA produced under specific circumstances or in a cell, tissue, organ, or an entire organism, in a given moment of its development. Several RNA types are present in the transcriptome, such as mRNA, tRNA, rRNA, snRNA, snoRNA, miRNA, lncRNA, and pseudogenes (Wang et al., 2017Wang, J.; Meng, X.; Dobrovolskaya, O. B.; Orlov, Y. L.; Chen, M. Non-coding RNAs and their roles in stress response in plants. Genomics, Proteomics & Bioinformatics , v.15, p.301-312, 2017. https://doi.org/10.1016/j.gpb.2017.01.007
https://doi.org/10.1016/j.gpb.2017.01.00...
). In recent years, due to the technological advances and cost reduction achieved with the RNA-seq technique, we have witnessed an explosion in the amount of transcriptome data generated and made public (Lowe et al., 2017Lowe, R.; Shirley, N.; Bleackley, M.; Dolan, S.; Shafee, T. Transcriptomics technologies. PLoS Computational Biology, v.13, p.1-23, 2017. https://doi.org/10.1371/journal.pcbi.1005457
https://doi.org/10.1371/journal.pcbi.100...
).

The current study is built upon the study by Vieira et al. (2020Vieira, L. R.; Silva, V. N. B.; Casari, R. A. das C. N.; Carmona, P. A. O.; Sousa, C. A. F. de; Souza Junior, M. T. Morphophysiological responses of young oil palm plants to salinity stress. Pesquisa Agropecuária Brasileira, v.55, p.1-15, 2020. https://doi.org/10.1590/s1678-3921.pab2020.v55.01835
https://doi.org/10.1590/s1678-3921.pab20...
), who evaluated the morphophysiological responses and ionic imbalance in the substrate, roots, and leaves of young African oil palm plants (Elaeis guineensisJacq.) under different NaCl doses following a rigorous substrate salinization protocol. The present study carried out a comprehensive large-scale transcriptome analysis aimed to prospect and annotate salt-responsive genes in the genome of the African oil palm that were highly up-regulated during salt stress, with a long-term goal of gene promoter prospection and validation.

Material and Methods

All oil palm plants used in the study were regenerated from embryogenic calluses (AM33 genotype). Young oil palm plants at the growth stage known as “bifid saplings” were subjected to different doses of NaCl in March 2018 and maintained under these conditions for 12 days in a completely randomized experimental design at a greenhouse in Brasília, DF, Brazil (15.732° S, 47.900° W, 1,030 meters of altitude) (Vieira et al., 2020Vieira, L. R.; Silva, V. N. B.; Casari, R. A. das C. N.; Carmona, P. A. O.; Sousa, C. A. F. de; Souza Junior, M. T. Morphophysiological responses of young oil palm plants to salinity stress. Pesquisa Agropecuária Brasileira, v.55, p.1-15, 2020. https://doi.org/10.1590/s1678-3921.pab2020.v55.01835
https://doi.org/10.1590/s1678-3921.pab20...
).

Based on the morphophysiological responses of the young oil palm plants to salinity stress, the apical leaves from three stressed plants (~40 dS m-1 of electrical conductivity) were collected at 5 and 12 days after the stress onset (DAT), together with the apical leaves from three control plants at 12 days (~2 dS m-1), immediately frozen in liquid nitrogen and stored at -80 °C until total RNA extraction.

Total RNA samples were subjected to RNA-Seq using an Illumina HiSeq platform at the GenOne Company (Rio de Janeiro, RJ, Brazil), using the paired-end strategy. The fastq files with the 150-nt-long sequences were deposited in the Sequence Read Archive (SRA) database of the National Center for Biotechnology Information - BioProject PRJNA573093, BioSample SAMN12799239.

All the RNAseq analyses went through the OmicsBox bioinformatics platform (OmicsBox - Bioinformatics Made Easy, BioBam Bioinformatics, March 3, 2019, https://www.biobam.com/omicsbox), in version 1.2.4., module transcriptomics, using the adjusted Reference Based Transcriptomics Workflow. Only high-quality sequences with a minimum average quality of Q30 and a minimum length of 75 nts underwent mapping against the African oil palm reference genome (Singh et al., 2013Singh, R.; Ong-Abdullah, M.; Low, E. T.; Manaf, M. A.; Rosli, R.; Nookiah, R.; Ooi, L. C.; Ooi, S. E.; Chan, K. L.; Halim, M. A.; Azizi, N.; Nagappan, J.; Bacher, B.; Lakey, N.; Smith, S. W.; He, D.; Hogan, M.; Budiman, M. A.; Lee, E. K.; DeSalle, R.; Sambanthamurthi, R. Oil palm genome sequence reveals divergence of interfertile species in old and new worlds. Nature, v.500, p.335-339, 2013. https://doi.org/10.1038/nature12309
https://doi.org/10.1038/nature12309...
) - files downloaded from the National Center for Biotechnology Information (BioProject PRJNA192219; BioSample SAMN02981535) on October 2020 - using the software STAR (Dobin et al., 2013Dobin, A.; Davis, C. A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T. R. STAR: ultrafast universal RNA-seq aligner. Bioinformatics , v.29, p.15-21, 2013. https://doi.org/10.1093/bioinformatics/bts635
https://doi.org/10.1093/bioinformatics/b...
).

The nine BAM files generated during the mapping stage were submitted to the generation of the Count Table (BAM + GFF) using the HTSeq software version 0.9.0 (Anders et al., 2015Anders, S.; Pyl, P. T.; Huber, W. HTSeq - a Python framework to work with high-throughput sequencing data. Bioinformatics, v.31, p.166-169, 2015. https://doi.org/10.1093/bioinformatics/btu638
https://doi.org/10.1093/bioinformatics/b...
), and pairwise differential expression analysis using the software package edgeR version 3.28.0 (Robinson et al., 2010Robinson, M. D.; McCarthy, D. J.; Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics , v.26, p.139-140, 2010. https://doi.org/10.1093/bioinformatics/btp616
https://doi.org/10.1093/bioinformatics/b...
).

The false discovery rate (FDR) and Log2 (Fold Change) were used as criteria to select the genes responsive to salinity stress - FDR ≤ 0.01 and Log2 (Fold Change) ≥ 5.0. FDR is a statistical approach used in multiple hypothesis testing to correct for multiple comparisons, and it is typically used in high-throughput experiments to correct for random events that falsely appear significant (Benjamini & Hochberg, 1995Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B (Methodological), v.57, p.289-300, 1995. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
https://doi.org/10.1111/j.2517-6161.1995...
). Genes differentially expressed (DEGs) at 5 and 12 days after imposing the salinity stress (DAT) underwent further characterization based on the African oil palm genome annotated and available at NCBI (Singh et al., 2013Singh, R.; Ong-Abdullah, M.; Low, E. T.; Manaf, M. A.; Rosli, R.; Nookiah, R.; Ooi, L. C.; Ooi, S. E.; Chan, K. L.; Halim, M. A.; Azizi, N.; Nagappan, J.; Bacher, B.; Lakey, N.; Smith, S. W.; He, D.; Hogan, M.; Budiman, M. A.; Lee, E. K.; DeSalle, R.; Sambanthamurthi, R. Oil palm genome sequence reveals divergence of interfertile species in old and new worlds. Nature, v.500, p.335-339, 2013. https://doi.org/10.1038/nature12309
https://doi.org/10.1038/nature12309...
).

Results and Discussion

A total of 212,241,278 pairs of high-quality sequences underwent mapping against the African oil palm genome at an average of 23,582,364 per sample. Over 99% of these sequences mapped to the African oil palm reference genome, showing no contamination by RNA from other organisms. Finally, on average, approximately 74% of the sequence pairs mapped just once to genes in the African oil palm genome (Table 1). The African oil palm reference genome has 29,567 genomic features of type ‘gene’; however, 4,213 of these features had no aligned reads detected in any of the samples generated in this study (data not shown).

Table 1
Statistics of RNAseq data from the samples of oil palm plants submitted to two salt stress treatments; three replicates per treatment. Samples from salt-stressed plants collected at 5 and 12 days after treatment (DAT); sample from control plants collected at 12 DAT. Mapping to reference genome EG5 (BioProject PRJNA192219 and BioSample SAMN02981535) available at the National Center for Biotechnology Information

When comparing the control against stressed plants on the fifth DAT, the pairwise differential expression analysis revealed that out of the 25,354 genes from theE. guineensisgenome (Singh et al., 2013Singh, R.; Ong-Abdullah, M.; Low, E. T.; Manaf, M. A.; Rosli, R.; Nookiah, R.; Ooi, L. C.; Ooi, S. E.; Chan, K. L.; Halim, M. A.; Azizi, N.; Nagappan, J.; Bacher, B.; Lakey, N.; Smith, S. W.; He, D.; Hogan, M.; Budiman, M. A.; Lee, E. K.; DeSalle, R.; Sambanthamurthi, R. Oil palm genome sequence reveals divergence of interfertile species in old and new worlds. Nature, v.500, p.335-339, 2013. https://doi.org/10.1038/nature12309
https://doi.org/10.1038/nature12309...
), 35 genes were differentially expressed using the following selection criteria: FDR ≤ 0.01 and Log2 (Fold Change) ≥ 5. On the other hand, when comparing control against stressed plants, on the 12 DAT, the pairwise differential expression analysis revealed 11 DEGs. Seven genes were present in both groups of DEGs (Figure 1).

Figure 1
Heatmap of the genes highly overexpressed in the apical leaf of young oil palm plants under salinity stress, 5 and 14 days after treatment (DAT). Venn diagram showing that 11 genes were up-regulated at 12 DAT and 35 at 5 DAT, with seven genes highly overexpressed in both scenarios

The small number of DEGs revealed in this study is a direct consequence of the high Log2 (Fold Change) used, where only DEGs with a fold change of 32 times or more were selected. Another fact that contributed to this small number of DEGs is the FDR used in this study, i.e., ≤ 0.01. FDR offers more confidence than just p.value, and, usually, results below 0.05 are considered real (Benjamini & Hochberg, 1995Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B (Methodological), v.57, p.289-300, 1995. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
https://doi.org/10.1111/j.2517-6161.1995...
).

The reason to select such values for Log2 (Fold Change) and FDR in this study was the intention to select salt-responsive genes with high increases in the expression level for future studies on gene promoter prospection and validation. The database of salt-responsive genes generated in this present study is a reservoir of data that can be used to respond to several questions regarding the response of oil palm plants to salt stress. One question, for instance, could be whether the strategy of salt stress employed by Vieira et al. (2020Vieira, L. R.; Silva, V. N. B.; Casari, R. A. das C. N.; Carmona, P. A. O.; Sousa, C. A. F. de; Souza Junior, M. T. Morphophysiological responses of young oil palm plants to salinity stress. Pesquisa Agropecuária Brasileira, v.55, p.1-15, 2020. https://doi.org/10.1590/s1678-3921.pab2020.v55.01835
https://doi.org/10.1590/s1678-3921.pab20...
) was able to render the plants at 5 DAT to the osmotic stress and at 12 DAT to the ionic stress.

These seven DEGs underwent structural and functional annotation. Just one of the seven DEGs was in a scaffold still not linked to any chromosome in the African oil palm genome. These DEGs had between one and 12 exons, and all code for proteins. Two of them were in the sense strand, and the remaining in the anti-sense strand. These seven genes coded for a thiosulfate sulfurtransferase, a non-specific lipid transfer protein, an acid endochitinase, a structural protein of the cell wall rich in glycine, a transcription factor, a probable amino acid transporter in the vacuole, and a carbonic anhydrase (Table 2).

Table 2
Gene ontology classification of the genes highly overexpressed in the apical leaf of young oil palm plants under salinity stress, after 12 days of treatment. Genes were classified into biological process, molecular function, and cellular component at the second level of gene ontology classification

LOC105055832 is a homolog of the probable vacuolar amino acid transporter YPQ1 gene in the African oil palm genome. Transport proteins from three major families - ATF (amino acid transporter family), APC (amino acid-polyamine-choline transporter family), and UMAMIT (usually multiple acids move in and out transporter family) - mediate the transport and allocation of amino acids within plants (Yang et al., 2020Yang, G.; Wei, Q.; Huang, H.; Xia, J. Amino acid transporters in plant cells: A brief review. Plants, v.9, p.1-17, 2020. https://doi.org/10.3390/plants9080967
https://doi.org/10.3390/plants9080967...
).

The ATF subfamily has different groups of proteins: AAPs (amino acid permeases), LHTs (lysine and histidine transporters), GATs (γ-aminobutyric acid transporters), ProTs (proline transporters), AUXs (indole-3-acetic acid transporters), aromatic and neutral amino acid transporters, and amino acid transporter-like proteins (Yao et al., 2020Yao, X.; Nie, J.; Bai, R.; Sui, X. Amino acid transporters in plants: Identification and function. Plants, v.9, p.1-17, 2020. https://doi.org/10.3390/plants9080972
https://doi.org/10.3390/plants9080972...
).

Several different types of amino acid transporters transport proline, an important compatible solute known to play a role in protecting plants from abiotic stresses, such as salinity, drought, and freezing, and altered expression of some of these transporters can enhance plant tolerance to salt and drought stresses (Yao et al., 2020Yao, X.; Nie, J.; Bai, R.; Sui, X. Amino acid transporters in plants: Identification and function. Plants, v.9, p.1-17, 2020. https://doi.org/10.3390/plants9080972
https://doi.org/10.3390/plants9080972...
).

LOC105055156 is a homolog of the thiosulfate sulfurtransferase 16, chloroplastic, gene in the African oil palm genome. Sulfurtransferases/rhodaneses (Str) comprise a large and complex group of enzymes that catalyze the transfer of sulfur from thiosulfate to cyanide, leading to the formation of thiocyanate and sulfite (Papenbrock et al., 2011Papenbrock, J.; Guretzki, S.; Henne, M. Latest news about the sulfurtransferase protein family of higher plants. Amino Acids, v.41, p.43-57, 2011. https://doi.org/10.1007/s00726-010-0478-6
https://doi.org/10.1007/s00726-010-0478-...
).

Hatzfeld & Saito (2000Hatzfeld, Y.; Saito, K. Evidence for the existence of rhodanese (thiosulfate:cyanide sulfurtransferase) in plants: preliminary characterization of two rhodanese cDNAs from Arabidopsis thaliana. FEBS Letters, v.470, p.147-150, 2000. https://doi.org/10.1016/S0014-5793(00)01311-9
https://doi.org/10.1016/S0014-5793(00)01...
) were the first ones to report the existence of rhodaneses in plants, and Mao et al. (2011Mao, G.; Wang, R.; Guan, Y.; Liu, Y.; Zhang, S. Sulfurtransferases 1 and 2 play essential roles in embryo and seed development in Arabidopsis thaliana. The Journal of Biological Chemistry, v.286, p.7548-7557, 2011. https://doi.org/10.1074/jbc.M110.182865
https://doi.org/10.1074/jbc.M110.182865...
) showed that STR1/STR2 sulfurtransferases play a role in plant embryo/seed development, highlighting the important biological function of this group of ubiquitous enzymes may have in plants. More recently, Moseler et al. (2019Moseler, A.; Selles, B.; Rouhier, N.; Couturier, J. Novel insights into the diversity of the sulfurtransferase family in photosynthetic organisms with emphasis on oak. New Phytologist, v.226, p.967-977, 2019. https://doi.org/10.1111/nph.15870
https://doi.org/10.1111/nph.15870...
) classified this protein family into nine clusters depending on their primary sequence and domain arrangement based on the expansion of knowledge regarding STRs in higher photosynthetic organisms.

LOC105045973 is a homolog of a non-specific lipid transfer protein 1 gene in the African oil palm genome. Plant non-specific lipid transfer proteins (LTPs) are small soluble proteins that facilitate the transfer of fatty acids, phospholipids, glycolipids, and steroids between membranes. LTPs have roles in seed storage, lipid mobilization, cuticle synthesis, somatic embryogenesis, pollen tube adhesion, and plant tolerance to abiotic - drought, cold, and salt stress - and biotic stress - bacterial and fungal pathogens (Gangadhar et al., 2016Gangadhar, B. H.; Sajeesh, K.; Venkatesh, J.; Baskar, V.; Abhinandan, K.; Yu, J. W.; Prasad, R.; Mishra, R. K. Enhanced tolerance of transgenic potato plants over-expressing non-specific lipid transfer protein-1 (StnsLTP1) against multiple abiotic stresses. Frontiers in Plant Science , v.7, p.1-12, 2016. https://doi.org/10.3389/fpls.2016.01228
https://doi.org/10.3389/fpls.2016.01228...
). Over-expression ofnsLTPgenes led to improved tolerance to salt, drought, and cold stress in Arabidopsis and potato (Guo et al., 2013Guo, L.; Yang, H.; Zhang, X.; Yang, S. Lipid transfer protein 3 as a target of MYB96 mediates freezing and drought stress in Arabidopsis. Journal of Experimental Botany, v.64, p.1755-1767, 2013. https://doi.org/10.1093/jxb/ert040
https://doi.org/10.1093/jxb/ert040...
; Zou et al., 2013Zou, H. W.; Tian, X. H.; Ma, G. H.; Li, Z. X. Isolation and functional analysis of ZmLTP3, a homologue to Arabidopsis LTP3. International Journal of Molecular Sciences, v.14, p.5025-5035, 2013. https://doi.org/10.3390/ijms14035025
https://doi.org/10.3390/ijms14035025...
; Gangadhar et al., 2016Gangadhar, B. H.; Sajeesh, K.; Venkatesh, J.; Baskar, V.; Abhinandan, K.; Yu, J. W.; Prasad, R.; Mishra, R. K. Enhanced tolerance of transgenic potato plants over-expressing non-specific lipid transfer protein-1 (StnsLTP1) against multiple abiotic stresses. Frontiers in Plant Science , v.7, p.1-12, 2016. https://doi.org/10.3389/fpls.2016.01228
https://doi.org/10.3389/fpls.2016.01228...
).

LOC105047536 is a homolog of an acidic endochitinase gene in the African oil palm genome. Chitinases are glycosyl hydrolases that catalyze the degradation of chitin, one of the most abundant biopolymers in nature, and genes encoding proteins with homology to chitinases have been identified from a wide range of organisms, including organisms that do not contain chitin, such as bacteria, viruses, animals and plants (Kwon et al., 2007Kwon, Y.; Kim, S. H.; Jung, M. S.; Kim, M. S.; Oh, J. E.; Ju, H. W.; Kim, K. I.;Vierling, E.; Lee, H.; Hong, S. W. Arabidopsis hot2 encodes an endochitinase-like protein that is essential for tolerance to heat, salt and drought stresses. The Plant Journal, v.49, p.184-193, 2007. https://doi.org/10.1111/j.1365-313X.2006.02950.x
https://doi.org/10.1111/j.1365-313X.2006...
).

Endogenous genes encoding proteins belonging to the plant chitinase families are involved in abiotic stress responses, and the overexpression of chitinase genes in plants - whether fungal - or plant-derived genes - does increase tolerance to abiotic stresses, such as salinity and heavy metals (Kwon et al., 2007Kwon, Y.; Kim, S. H.; Jung, M. S.; Kim, M. S.; Oh, J. E.; Ju, H. W.; Kim, K. I.;Vierling, E.; Lee, H.; Hong, S. W. Arabidopsis hot2 encodes an endochitinase-like protein that is essential for tolerance to heat, salt and drought stresses. The Plant Journal, v.49, p.184-193, 2007. https://doi.org/10.1111/j.1365-313X.2006.02950.x
https://doi.org/10.1111/j.1365-313X.2006...
; Brotman et al., 2012Brotman, Y.; Landau, U.; Pnini, S.; Lisec, J.; Balazadeh, S.; Mueller-Roeber, B.; Zilberstein, A.; Willmitzer, L.; Chet, I.; Viterbo, A. The LysM Receptor-Like Kinase LysM RLK1is required to activate defense and abiotic-stress responses induced by overexpression of fungal chitinases in Arabidopsis plants. Molecular Plant, v.5, p.1113-24, 2012. https://doi.org/10.1093/mp/sss021
https://doi.org/10.1093/mp/sss021...
).

LOC105058460 is a homolog of the transcription factor bHLH94 gene in the African oil palm genome. The basic helix-loop-helix (bHLH) transcription factors are a large gene family in the plant genome involved in regulating plant responses to abiotic stresses, including salt stress. Recent studies have shown that the overexpression of this type of genes increased the tolerance to drought and salinity stress in tomato and Arabidopsis plants (Waseem et al., 2019Waseem, M.; Rong, X.; Li, Z. Dissecting the role of a basic helix-loop-helix transcription factor, SlbHLH22, under salt and drought stresses in transgenic Solanum lycopersicum L. Frontiers in Plant Science , v.10, p.734, 2019. https://doi.org/10.3389/fpls.2019.00734
https://doi.org/10.3389/fpls.2019.00734...
; Qiu et al., 2020Qiu, J. R.; Huang, Z.; Xiang, X. Y.; Xu, W. X.; Wang, J. T.; Chen, J.; Song, L.; Xiao, Y.; Li, X.; Ma, J.; Cai, S. Z.; Sun, L. X.; Jiang, C. Z. MfbHLH38, a Myrothamnus flabellifolia bHLH transcription factor, confers tolerance to drought and salinity stresses in Arabidopsis. BMC Plant Biology, v.20, p.542, 2020. https://doi.org/10.1186/s12870-020-02732-6
https://doi.org/10.1186/s12870-020-02732...
).

LOC105039727 is a homolog of a glycine-rich cell wall structural protein 2 gene in the African oil palm genome. In plants, glycine-rich proteins (GRPs) are characterized by the presence of semi-repetitive glycine-rich motifs and classified into five distinct groups - Class I, II, III, IV, and V (Mangeon et al., 2010Mangeon, A.; Junqueira, R. M.; Sachetto-Martins, G. Functional diversity of the plant glycine-rich proteins superfamily. Plant Signaling & Behavior, v.5, p.99-104, 2010. https://doi.org/10.4161/psb.5.2.10336
https://doi.org/10.4161/psb.5.2.10336...
). Czolpinska & Rurek (2018Czolpinska, M.; Rurek, M. Plant glycine-rich proteins in stress response: An emerging, still prospective story. Frontiers in Plant Science, v.9, p.1-13, 2018. https://doi.org/10.3389/fpls.2018.00302
https://doi.org/10.3389/fpls.2018.00302...
) postulate that GRPs could play a promising role in agriculture through plant genetic engineering in the coming years because some studies have shown that their overexpression can confer tolerance to cold, drought, and salinity stresses.

LOC105031978 is a homolog of an alpha carbonic anhydrase 7 gene in the African oil palm genome. Plants contain three evolutionarily distinct carbonic anhydrases (CAs) families - αCAs, βCAs, and γCAs, and multiple genes encoding all three types of CAs. CAs are zinc metalloenzymes found in numerous plant tissues and different cellular locations that catalyze the interconversion of CO2 and HCO3 - (DiMario et al., 2017DiMario, R. J.; Clayton, H.; Mukherjee, A.; Ludwig, M.; Moroney, J. V. Plant carbonic anhydrases: Structures, locations, evolution, and physiological roles. Molecular Plant, v.10, p.30-46, 2017. https://doi.org/10.1016/j.molp.2016.09.001
https://doi.org/10.1016/j.molp.2016.09.0...
).

Transgenic Arabidopsis plants overexpressing aβCAgene from rice showed improved salt tolerance (Yu et al., 2007Yu, S.; Zhang, X.; Guan, Q.; Takano, T.; Liu, S. Expression of a carbonic anhydrase gene is induced by environmental stresses in rice (Oryza sativa L.). Biotechnology Letters, v.29, p.89-94, 2007. https://doi.org/10.1007/s10529-006-9199-z
https://doi.org/10.1007/s10529-006-9199-...
). αCAs are the largest CA gene family in most plants, but they are also the least studied most likely because they are not highly abundant in leaves and roots.

The expression values of the seven genes - highly up-regulated either at 5 and 12 DAT - showed a fold change between 37 and 84 times in the salt-stressed plants at 12 DAT in comparison to the control plants. The levels of expression at 5 DAT from all seven genes were higher than at 12 DAT (Figure 2).

Figure 2
Expression profiles in fold change (x-axis) of the genes highly overexpressed in the apical leaf of young oil palm plants under salinity stress, after 5 and 14 days of treatment (DAT), and selected for annotation. Gene expression in fold change in the apical leaf of control plants (blue), plants salt-stressed for 5 days (green), and 12 days (red)

If the goal of a transcriptomics study, such as this present study, is also to select candidate salt-responsive genes for future gene promoter prospection and validation, fold change alone is not a secure criterion, as it has to take into consideration the expression value of a gene in the control plant. All seven genes had expression values close to zero (data not shown).

Based on the analysis of the expression values at 12 DAT, these seven genes were divided into two groups: those with an expression value lower than 100, and those higher than 100 (Table 3). Taken together, these results put the three genes in the latter group as candidate genes for future studies aiming at the prospection and validation of salt-responsive promoter sequences.

As pointed out above, five out of these seven salt-responsive genes code for proteins with studies showing that their overexpression increased the tolerance to salinity stress in several plant species (Kwon et al., 2007Kwon, Y.; Kim, S. H.; Jung, M. S.; Kim, M. S.; Oh, J. E.; Ju, H. W.; Kim, K. I.;Vierling, E.; Lee, H.; Hong, S. W. Arabidopsis hot2 encodes an endochitinase-like protein that is essential for tolerance to heat, salt and drought stresses. The Plant Journal, v.49, p.184-193, 2007. https://doi.org/10.1111/j.1365-313X.2006.02950.x
https://doi.org/10.1111/j.1365-313X.2006...
; Yu et al., 2007Yu, S.; Zhang, X.; Guan, Q.; Takano, T.; Liu, S. Expression of a carbonic anhydrase gene is induced by environmental stresses in rice (Oryza sativa L.). Biotechnology Letters, v.29, p.89-94, 2007. https://doi.org/10.1007/s10529-006-9199-z
https://doi.org/10.1007/s10529-006-9199-...
; Brotman et al., 2012Brotman, Y.; Landau, U.; Pnini, S.; Lisec, J.; Balazadeh, S.; Mueller-Roeber, B.; Zilberstein, A.; Willmitzer, L.; Chet, I.; Viterbo, A. The LysM Receptor-Like Kinase LysM RLK1is required to activate defense and abiotic-stress responses induced by overexpression of fungal chitinases in Arabidopsis plants. Molecular Plant, v.5, p.1113-24, 2012. https://doi.org/10.1093/mp/sss021
https://doi.org/10.1093/mp/sss021...
; Guo et al., 2013Guo, L.; Yang, H.; Zhang, X.; Yang, S. Lipid transfer protein 3 as a target of MYB96 mediates freezing and drought stress in Arabidopsis. Journal of Experimental Botany, v.64, p.1755-1767, 2013. https://doi.org/10.1093/jxb/ert040
https://doi.org/10.1093/jxb/ert040...
; Zou et al., 2013Zou, H. W.; Tian, X. H.; Ma, G. H.; Li, Z. X. Isolation and functional analysis of ZmLTP3, a homologue to Arabidopsis LTP3. International Journal of Molecular Sciences, v.14, p.5025-5035, 2013. https://doi.org/10.3390/ijms14035025
https://doi.org/10.3390/ijms14035025...
; Gangadhar et al., 2016Gangadhar, B. H.; Sajeesh, K.; Venkatesh, J.; Baskar, V.; Abhinandan, K.; Yu, J. W.; Prasad, R.; Mishra, R. K. Enhanced tolerance of transgenic potato plants over-expressing non-specific lipid transfer protein-1 (StnsLTP1) against multiple abiotic stresses. Frontiers in Plant Science , v.7, p.1-12, 2016. https://doi.org/10.3389/fpls.2016.01228
https://doi.org/10.3389/fpls.2016.01228...
; Czolpinska & Rurek, 2018Czolpinska, M.; Rurek, M. Plant glycine-rich proteins in stress response: An emerging, still prospective story. Frontiers in Plant Science, v.9, p.1-13, 2018. https://doi.org/10.3389/fpls.2018.00302
https://doi.org/10.3389/fpls.2018.00302...
; Waseem et al., 2019Waseem, M.; Rong, X.; Li, Z. Dissecting the role of a basic helix-loop-helix transcription factor, SlbHLH22, under salt and drought stresses in transgenic Solanum lycopersicum L. Frontiers in Plant Science , v.10, p.734, 2019. https://doi.org/10.3389/fpls.2019.00734
https://doi.org/10.3389/fpls.2019.00734...
; Qiu et al., 2020Qiu, J. R.; Huang, Z.; Xiang, X. Y.; Xu, W. X.; Wang, J. T.; Chen, J.; Song, L.; Xiao, Y.; Li, X.; Ma, J.; Cai, S. Z.; Sun, L. X.; Jiang, C. Z. MfbHLH38, a Myrothamnus flabellifolia bHLH transcription factor, confers tolerance to drought and salinity stresses in Arabidopsis. BMC Plant Biology, v.20, p.542, 2020. https://doi.org/10.1186/s12870-020-02732-6
https://doi.org/10.1186/s12870-020-02732...
). Three of these five genes showed an expression value lower than 100, and two higher than 100 at 12 DAT (Table 3). Taken together, these results put these five genes as candidate genes for future studies aiming at the generation of oil palm plants tolerant to salt stress, via a constitutive overexpression strategy.

Table 3
Statistics from the differential expression analysis of the genes highly overexpressed in the apical leaf of young oil palm plants under salinity stress, after 5 and 12 days of treatment, in comparison to the control plants

Conclusions

  1. The prospection and annotation strategy applied allowed the selection of seven salt-responsive genes from oil palm plants showing increased expression levels as a result of stress with approximately 40 dS m-1.

  2. Three out of those seven salt-responsive genes are candidate genes for future studies aiming at the prospection and validation of salt-responsive promoter sequences.

  3. Five out of those seven salt-responsive genes are candidate genes for future studies aiming at the generation of oil palm plants tolerant to salt stress, via a strategy of heterologous constitutive overexpression.

Acknowledgments

The authors acknowledge a) funding to T.M.M.F. by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), a Foundation within the Ministry of Education in Brazil, via the Graduate Program in Plant Biotechnology, Universidade Federal de Lavras (UFLA); and b) The grant (01.13.0315.00 DendêPalm Project) for this study was awarded by the Ministério da Ciência e Tecnologia do Brasil, e Inovação (MCTI) via the Agência Brasileira de Inovação (FINEP). The authors confirm that the funder did not influence the study design, the content of the article, or the selection of this journal.

Literature Cited

  • Anders, S.; Pyl, P. T.; Huber, W. HTSeq - a Python framework to work with high-throughput sequencing data. Bioinformatics, v.31, p.166-169, 2015. https://doi.org/10.1093/bioinformatics/btu638
    » https://doi.org/10.1093/bioinformatics/btu638
  • Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B (Methodological), v.57, p.289-300, 1995. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
    » https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
  • Brotman, Y.; Landau, U.; Pnini, S.; Lisec, J.; Balazadeh, S.; Mueller-Roeber, B.; Zilberstein, A.; Willmitzer, L.; Chet, I.; Viterbo, A. The LysM Receptor-Like Kinase LysM RLK1is required to activate defense and abiotic-stress responses induced by overexpression of fungal chitinases in Arabidopsis plants. Molecular Plant, v.5, p.1113-24, 2012. https://doi.org/10.1093/mp/sss021
    » https://doi.org/10.1093/mp/sss021
  • Czolpinska, M.; Rurek, M. Plant glycine-rich proteins in stress response: An emerging, still prospective story. Frontiers in Plant Science, v.9, p.1-13, 2018. https://doi.org/10.3389/fpls.2018.00302
    » https://doi.org/10.3389/fpls.2018.00302
  • DiMario, R. J.; Clayton, H.; Mukherjee, A.; Ludwig, M.; Moroney, J. V. Plant carbonic anhydrases: Structures, locations, evolution, and physiological roles. Molecular Plant, v.10, p.30-46, 2017. https://doi.org/10.1016/j.molp.2016.09.001
    » https://doi.org/10.1016/j.molp.2016.09.001
  • Dobin, A.; Davis, C. A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T. R. STAR: ultrafast universal RNA-seq aligner. Bioinformatics , v.29, p.15-21, 2013. https://doi.org/10.1093/bioinformatics/bts635
    » https://doi.org/10.1093/bioinformatics/bts635
  • Duarte, B.; Caçador, I. Iberian halophytes as agroecological solutions for degraded lands and biosaline agriculture. Sustainability, v.13, p.1-14, 2021. https://doi.org/10.3390/su13021005
    » https://doi.org/10.3390/su13021005
  • Gangadhar, B. H.; Sajeesh, K.; Venkatesh, J.; Baskar, V.; Abhinandan, K.; Yu, J. W.; Prasad, R.; Mishra, R. K. Enhanced tolerance of transgenic potato plants over-expressing non-specific lipid transfer protein-1 (StnsLTP1) against multiple abiotic stresses. Frontiers in Plant Science , v.7, p.1-12, 2016. https://doi.org/10.3389/fpls.2016.01228
    » https://doi.org/10.3389/fpls.2016.01228
  • Guo, L.; Yang, H.; Zhang, X.; Yang, S. Lipid transfer protein 3 as a target of MYB96 mediates freezing and drought stress in Arabidopsis. Journal of Experimental Botany, v.64, p.1755-1767, 2013. https://doi.org/10.1093/jxb/ert040
    » https://doi.org/10.1093/jxb/ert040
  • Hatzfeld, Y.; Saito, K. Evidence for the existence of rhodanese (thiosulfate:cyanide sulfurtransferase) in plants: preliminary characterization of two rhodanese cDNAs from Arabidopsis thaliana FEBS Letters, v.470, p.147-150, 2000. https://doi.org/10.1016/S0014-5793(00)01311-9
    » https://doi.org/10.1016/S0014-5793(00)01311-9
  • Joshi, A.; Kanthaliya, B.; Arora, J. Halophytes: The nonconventional crops as source of biofuel production. In: Grigore, M. N. (eds.). Handbook of halophytes. Cham: Springer, v.1, p.1-28, 2020. https://doi.org/10.1007/978-3-030-17854-3_126-1
    » https://doi.org/10.1007/978-3-030-17854-3_126-1
  • Kwon, Y.; Kim, S. H.; Jung, M. S.; Kim, M. S.; Oh, J. E.; Ju, H. W.; Kim, K. I.;Vierling, E.; Lee, H.; Hong, S. W. Arabidopsis hot2 encodes an endochitinase-like protein that is essential for tolerance to heat, salt and drought stresses. The Plant Journal, v.49, p.184-193, 2007. https://doi.org/10.1111/j.1365-313X.2006.02950.x
    » https://doi.org/10.1111/j.1365-313X.2006.02950.x
  • Lowe, R.; Shirley, N.; Bleackley, M.; Dolan, S.; Shafee, T. Transcriptomics technologies. PLoS Computational Biology, v.13, p.1-23, 2017. https://doi.org/10.1371/journal.pcbi.1005457
    » https://doi.org/10.1371/journal.pcbi.1005457
  • Mangeon, A.; Junqueira, R. M.; Sachetto-Martins, G. Functional diversity of the plant glycine-rich proteins superfamily. Plant Signaling & Behavior, v.5, p.99-104, 2010. https://doi.org/10.4161/psb.5.2.10336
    » https://doi.org/10.4161/psb.5.2.10336
  • Mao, G.; Wang, R.; Guan, Y.; Liu, Y.; Zhang, S. Sulfurtransferases 1 and 2 play essential roles in embryo and seed development in Arabidopsis thaliana The Journal of Biological Chemistry, v.286, p.7548-7557, 2011. https://doi.org/10.1074/jbc.M110.182865
    » https://doi.org/10.1074/jbc.M110.182865
  • Moseler, A.; Selles, B.; Rouhier, N.; Couturier, J. Novel insights into the diversity of the sulfurtransferase family in photosynthetic organisms with emphasis on oak. New Phytologist, v.226, p.967-977, 2019. https://doi.org/10.1111/nph.15870
    » https://doi.org/10.1111/nph.15870
  • Papenbrock, J.; Guretzki, S.; Henne, M. Latest news about the sulfurtransferase protein family of higher plants. Amino Acids, v.41, p.43-57, 2011. https://doi.org/10.1007/s00726-010-0478-6
    » https://doi.org/10.1007/s00726-010-0478-6
  • Qiu, J. R.; Huang, Z.; Xiang, X. Y.; Xu, W. X.; Wang, J. T.; Chen, J.; Song, L.; Xiao, Y.; Li, X.; Ma, J.; Cai, S. Z.; Sun, L. X.; Jiang, C. Z. MfbHLH38, a Myrothamnus flabellifolia bHLH transcription factor, confers tolerance to drought and salinity stresses in Arabidopsis. BMC Plant Biology, v.20, p.542, 2020. https://doi.org/10.1186/s12870-020-02732-6
    » https://doi.org/10.1186/s12870-020-02732-6
  • Robinson, M. D.; McCarthy, D. J.; Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics , v.26, p.139-140, 2010. https://doi.org/10.1093/bioinformatics/btp616
    » https://doi.org/10.1093/bioinformatics/btp616
  • Shahid, S. A.; Zaman, M.; Heng, L. Soil salinity: Historical perspectives and a world overview of the problem. In: Zaman, M.; Shahid, S. A.; Heng, L. (eds.). Guideline for salinity assessment, mitigation and adaptation using nuclear and related techniques. Cham: Springer , 2018. Chapter 2, p.43-53. https://doi.org/10.1007/978-3-319-96190-3_2
    » https://doi.org/10.1007/978-3-319-96190-3_2
  • Singh, R.; Ong-Abdullah, M.; Low, E. T.; Manaf, M. A.; Rosli, R.; Nookiah, R.; Ooi, L. C.; Ooi, S. E.; Chan, K. L.; Halim, M. A.; Azizi, N.; Nagappan, J.; Bacher, B.; Lakey, N.; Smith, S. W.; He, D.; Hogan, M.; Budiman, M. A.; Lee, E. K.; DeSalle, R.; Sambanthamurthi, R. Oil palm genome sequence reveals divergence of interfertile species in old and new worlds. Nature, v.500, p.335-339, 2013. https://doi.org/10.1038/nature12309
    » https://doi.org/10.1038/nature12309
  • Tıpırdamaz, R.; Karakas, S.; Dikilitas, M. Halophytes and the future of agriculture. In: Grigore, M. N. (eds.). Handbook of halophytes - From molecules to ecosystems towards biosaline agriculture. Cham: Springer , 2020. Chapter 47, p.1-15. https://doi.org/10.1007/978-3-030-17854-3_91-1
    » https://doi.org/10.1007/978-3-030-17854-3_91-1
  • Ventura, Y.; Eshel, A.; Pasternak, D.; Sagi, M. The development of halophyte-based agriculture: Past and present. Annals of Botany, v.115, p.529-540, 2015. https://doi.org/10.1093/aob/mcu173
    » https://doi.org/10.1093/aob/mcu173
  • Vieira, L. R.; Silva, V. N. B.; Casari, R. A. das C. N.; Carmona, P. A. O.; Sousa, C. A. F. de; Souza Junior, M. T. Morphophysiological responses of young oil palm plants to salinity stress. Pesquisa Agropecuária Brasileira, v.55, p.1-15, 2020. https://doi.org/10.1590/s1678-3921.pab2020.v55.01835
    » https://doi.org/10.1590/s1678-3921.pab2020.v55.01835
  • Wang, J.; Meng, X.; Dobrovolskaya, O. B.; Orlov, Y. L.; Chen, M. Non-coding RNAs and their roles in stress response in plants. Genomics, Proteomics & Bioinformatics , v.15, p.301-312, 2017. https://doi.org/10.1016/j.gpb.2017.01.007
    » https://doi.org/10.1016/j.gpb.2017.01.007
  • Waseem, M.; Rong, X.; Li, Z. Dissecting the role of a basic helix-loop-helix transcription factor, SlbHLH22, under salt and drought stresses in transgenic Solanum lycopersicum L. Frontiers in Plant Science , v.10, p.734, 2019. https://doi.org/10.3389/fpls.2019.00734
    » https://doi.org/10.3389/fpls.2019.00734
  • Yang, G.; Wei, Q.; Huang, H.; Xia, J. Amino acid transporters in plant cells: A brief review. Plants, v.9, p.1-17, 2020. https://doi.org/10.3390/plants9080967
    » https://doi.org/10.3390/plants9080967
  • Yao, X.; Nie, J.; Bai, R.; Sui, X. Amino acid transporters in plants: Identification and function. Plants, v.9, p.1-17, 2020. https://doi.org/10.3390/plants9080972
    » https://doi.org/10.3390/plants9080972
  • Yu, S.; Zhang, X.; Guan, Q.; Takano, T.; Liu, S. Expression of a carbonic anhydrase gene is induced by environmental stresses in rice (Oryza sativa L.). Biotechnology Letters, v.29, p.89-94, 2007. https://doi.org/10.1007/s10529-006-9199-z
    » https://doi.org/10.1007/s10529-006-9199-z
  • Zou, H. W.; Tian, X. H.; Ma, G. H.; Li, Z. X. Isolation and functional analysis of ZmLTP3, a homologue to Arabidopsis LTP3. International Journal of Molecular Sciences, v.14, p.5025-5035, 2013. https://doi.org/10.3390/ijms14035025
    » https://doi.org/10.3390/ijms14035025
  • 1 Research developed at Brasília, DF, Brazil

Edited by

Edited by: Hans Raj Gheyi

Publication Dates

  • Publication in this collection
    08 Sept 2021
  • Date of issue
    Dec 2021

History

  • Received
    19 Dec 2020
  • Accepted
    31 May 2021
  • Published
    01 July 2021
Unidade Acadêmica de Engenharia Agrícola Unidade Acadêmica de Engenharia Agrícola, UFCG, Av. Aprígio Veloso 882, Bodocongó, Bloco CM, 1º andar, CEP 58429-140, Campina Grande, PB, Brasil, Tel. +55 83 2101 1056 - Campina Grande - PB - Brazil
E-mail: revistagriambi@gmail.com