Next Article in Journal
MIFNN: Molecular Information Feature Extraction and Fusion Deep Neural Network for Screening Potential Drugs
Next Article in Special Issue
Genome-Wide Identification of Brassica napus PEN1-LIKE Genes and Their Expression Profiling in Insect-Susceptible and Resistant Cultivars
Previous Article in Journal
Design, Synthesis, and Bioactivity of Novel Quinazolinone Scaffolds Containing Pyrazole Carbamide Derivatives as Antifungal Agents
Previous Article in Special Issue
De Novo Transcriptome Assembly, Gene Annotations, and Characterization of Functional Profiling Reveal Key Genes for Lead Alleviation in the Pb Hyperaccumulator Greek Mustard (Hirschfeldia incana L.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 44
Issue 11
10.3390/cimb44110381
cimb-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Kandelia obovata, a Typical Mangrove Plant

by
Man-Man Sun
1,†,
Xiu Liu
2,†,
Xiao-Juan Huang
1,
Jing-Jun Yang
2,
Pei-Ting Qin
1,
Hao Zhou
1,
Ming-Guo Jiang
1 and
Hong-Ze Liao
1,*
1
Guangxi Key Laboratory of Polysaccharide Materials and Modification, School of Marine Sciences and Biotechnology, Guangxi Minzu University, 158 West Daxue Road, Nanning 530008, China
2
Guangxi Key Laboratory of Special Non-Wood Forest Cultivation and Utilization, Guangxi Forestry Research Institute, 23 Yongwu Road, Nanning 530002, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2022, 44(11), 5622-5637; https://doi.org/10.3390/cimb44110381
Submission received: 6 October 2022 / Revised: 10 November 2022 / Accepted: 11 November 2022 / Published: 13 November 2022
(This article belongs to the Special Issue Functional Genomics and Comparative Genomics Analysis in Plants)
Browse Figures
Versions Notes

Abstract

:
The NAC (NAM, ATAF1/2, and CUC2) gene family, one of the largest transcription factor families in plants, acts as positive or negative regulators in plant response and adaption to various environmental stresses, including cold stress. Multiple reports on the functional characterization of NAC genes in Arabidopsis thaliana and other plants are available. However, the function of the NAC genes in the typical woody mangrove (Kandelia obovata) remains poorly understood. Here, a comprehensive analysis of NAC genes in K. obovata was performed with a pluri-disciplinary approach including bioinformatic and molecular analyses. We retrieved a contracted NAC family with 68 genes from the K. obovata genome, which were unevenly distributed in the chromosomes and classified into ten classes. These KoNAC genes were differentially and preferentially expressed in different organs, among which, twelve up-regulated and one down-regulated KoNAC genes were identified. Several stress-related cis-regulatory elements, such as LTR (low-temperature response), STRE (stress response element), ABRE (abscisic acid response element), and WUN (wound-responsive element), were identified in the promoter regions of these 13 KoNAC genes. The expression patterns of five selected KoNAC genes (KoNAC6, KoNAC15, KoNAC20, KoNAC38, and KoNAC51) were confirmed by qRT-PCR under cold treatment. These results strongly implied the putative important roles of KoNAC genes in response to chilling and other stresses. Collectively, our findings provide valuable information for further investigations on the function of KoNAC genes.

1. Introduction

Transcription factors (TFs) are of immense importance due to their crucial impact on controlling the transcription rate by binding to the cis-regulatory elements, resulting in activation or inhibition of the transcription level of target genes [1]. There are numerous types of TF families in plants, among which the NAC (NAM, ATAF1/2, and CUC2) family serves as one of the largest plant-specific TF families and is named after the Petunia hybrida E. Vilm. NO APICAL MERISTEM (NAM) [2] and Arabidopsis thaliana (L.) Heynh. genes ATAF1/2 and CUP-SHAPED COTYLEDON 2 (CUC2) [3]. A typical NAC protein contains an N-terminal conserved NAC domain for DNA binding and nuclear localization and a variable C-terminal region with transcriptional regulatory activity [4].
As a complex plant-specific family, the NAC genes with considerable quantities are present in a wide range of species. A large number of NAC TFs have been identified in various plants, including A. thaliana [5], Actinidia eriantha Benth. [6], Asparagus officinalis L. [7], Betula pendula Roth [8], Hylocereus undatus (Haw.) Britton & Rose [9], Juglans mandshurica Maxim. [10], Medicago sativa L. [11], Miscanthus sinensis Andersson [12], banana (Musa acuminata Colla) [13], Oryza sativa L. [5], Populus trichocarpa Torr. & A. Gray ex Hook. [14], Salix psammophila C. Wang & Chang Y. Yang [15], Solanum lycopersicum L. [16], Zanthoxylum bungeanum Maxim. [17], and Zea mays L. [18]. Multiple lines of evidence illustrate that NAC genes act as positive or negative regulators involved in diverse biological processes, including plant response and adaptation to cold and other abiotic stresses [19,20]. MaNAC1, one banana NAC TF, acts as a downstream target of MaICE1 and interacts with the C-repeat binding factor MaCBF1, conferring fruit cold tolerance [21]. Two overexpressed NAC genes from H. undatus, HuNAC20 and HuNAC25, confer enhanced cold tolerance of transgenic A. thaliana plants [9]. Overexpression of MbNAC25 from Malus baccata (L.) Borkh. improves the resistance against chilling stress through enhanced scavenging capability of reactive oxygen species (ROS) in transgenic A. thaliana plants [22]. The tomato NAC gene NAM3 and its upstream regulator miR164a positively modulates cold tolerance by inducing ethylene synthesis in tomato plants [23]. CaNAC035, a novel NAC gene from Capsicum annuum, was shown to positively regulate cold stress in company with its upstream TF gene CabHLH79 [24]. CaNAC064, another NAC gene from C. annuum, is strongly induced by chilling stress and positively modulates cold stress tolerance via interacting with low temperature-induced haplo-proteinase proteins [25]. Additionally, NAC TFs also function as negative regulators in response to low temperature. Overexpression of MdNAC029, an apple NAC gene, reduces cold tolerance in apple and A. thaliana via a CBF-dependent pathway [26]. MaNAC25 and MaNAC28, two NAC genes from banana, negatively regulate cold tolerance in fruits by upregulating the expression levels of phospholipid degradation genes [27].
Mangroves are a dominant halophytic vegetation with significant ecological value in various tropical and subtropical coastal wetlands and are well-adapted to these highly stressful intertidal regions [28,29]. Among them, Kandelia obovata Sheue C.R., H.Y. Liu & J.W.H. Yong is regarded as a typical true mangrove due to its highest natural distribution latitude, indicating that K. obovata possesses stronger resistance against low temperature in contrast to other mangroves [30,31]. Various physiological evidences have shown that K. obovata displays better performance when exposed to chilling stress than other mangrove plants [32,33,34]. However, the underlying molecular mechanisms of cold response and adaptation in K. obovata are largely unknown. Here, we describe the genome-wide identification and expression analysis of K. obovata NAC (KoNAC) genes in response to low temperature based on its available chromosome-level reference genome [35] with a pluri-disciplinary approach including bioinformatic and molecular analyses, hopefully providing valuable insights into the function of NAC genes in cold response and breeding for cold resistance.

2. Materials and Methods

2.1. Identification and Chromosomal Distribution of NAC TFs in K. obovata

The K. obovata chromosome-scale genome (2n = 2x= 36) was obtained from Genome Warehouse (https://bigd.big.ac.cn/gwh) (accessed on 8 March 2022) under accession number GWHACBH00000000 [35]. The Hidden Markov Model (HMM) file for NAM domain (PF02365) was downloaded from Pfam database (https://pfam.xfam.org/) [36] (accessed on 8 March 2022), and was used to retrieve the NAC proteins with a cut-off value of 0.001 by HMMER 3.3.2 (http://hmmer.org/download.html) [37] (accessed on 8 March 2022). BLASTP (basic local alignment search tool for proteins) against K. obovata genome data with A. thaliana NAC protein sequences (Table S1) retrieved from The Arabidopsis Information Resource (TAIR, https://www.arabidopsis.org/) [38] (accessed on 10 March 2022) was implemented (e-value = 0.001). Taking these two results together, the final members of the KoNAC genes were acquired and verified by PfamScan (e-value = 0.001, https://www.ebi.ac.uk/Tools/pfa/pfamscan/) (accessed on 12 March 2022) [39] and NCBI’s conserved domain database (NCBI-CDD) (e-value = 0.001, https://www.ncbi.nlm.nih.gov/cdd/) (accessed on 12 March 2022) [40]. The basic information for KoNAC gene, including chromosome localization, intron number, average intron length, protein length, and isoelectric point (pI) values was determined based on the genome database. The chromosomal distribution map of KoNAC genes was drawn using MapChart 2.32 (https://www.wur.nl/en/show/Mapchart.htm) (accessed on 15 March 2022) [41].

2.2. Phylogenetic Analysis of NAC Proteins

The amino acid sequences of the NAC members of K. obovata and A. thaliana were aligned using Clustal X, and a neighbor-joining unrooted phylogenetic tree with 1000 bootstrap replications was constructed by MEGA 7.0 (www.megasoftware.net) (accessed on 15 March 2022) [42]. Finally, the tree was further modified by iTOL v6.5.8 (https://itol.embl.de/) (accessed on 15 March 2022) [43].

2.3. Gene Structure, Motif Identification, and Collinearity Analysis

The intron/exon structure of KoNAC genes was determined with the online gene structure display server (http://gsds.gao-lab.org/) (accessed on 18 March 2022) [44]. The conserved motifs in KoNAC proteins were identified by MEME suite v5.4.1 (http://meme-suite.org/) (accessed on 18 March 2022) [45]. The collinearity relationship of the K. obovata NAC genes between A. thaliana [5] and P. trichocarpa [14] were analyzed by MCScanX (http://chibba.pgml.uga.edu/mcscan2/) (accessed on 18 March 2022) [46]. These results were presented and visualized using TBtools (https://github.com/CJ-Chen/Tbtools) (accessed on 18 March 2022) [47].

2.4. Expression Analysis of KoNAC Genes Based on Public RNA-Seq Data

Two previously released RNA-seq data sets of K. obovata were introduced here to analyze the expression profiles of KoNAC genes. The expression patterns of KoNAC genes in eight organs (root, stem, leaf, flower, pistil, stamen, sepal, and fruit) were obtained according to the previously published transcriptomic data under the NCBI BioProject accession number PRJNA416402 (https://www.ncbi.nlm.nih.gov/bioproject) (accessed on 31 March 2022) [31]. The expression levels of KoNAC genes in response to cold stress were determined based on the publicly released data from the NCBI BioProject under accession number PRJNA678025. These two RNA-seq data were remapped back to the K. obovata genome used here [35]. All expression data were normalized as fragments per kilobase of transcript per million fragments mapped (FPKM) values [48]. The differentially expressed genes (DEGs) related to chilling stress were defined under the criteria of fold change (FC) ≥ 1.5. The expression profiles of KoNAC genes were visualized as heatmaps using TBtools [47].

2.5. Plant Materials and Treatment

The healthy mature propagules of the typical viviparous mangrove plant K. obovata were sampled from Guangxi Maoweihai Mangrove Nature Reserve, Qinzhou, China (21°37′23″ N, 108°44′13″ E) and cultured in the Mangrove Germplasm Resources Center (MGRC) of Guangxi Forestry Research Institute (GFRI) (Figure S1). The seedlings were grown in plastic pots containing sand and cultivated in a growth chamber at 28 °C and 75% humidity with a photoperiod of 14 h light/10 h darkness, and watered weekly with half-strength Hoagland’s nutrient solution [49]. At the eight-leaf stage, the seedlings were treated under low temperature (4 °C) for 0 h, 6 h, 12 h, and 24 h, respectively. All treatments were performed with three replicates. The leaves were harvested, immediately frozen in liquid nitrogen, and stored at −80 °C for RNA extraction.

2.6. Cis-Regulatory Element Analysis of the KoNAC Genes

The upstream 1500 bp promoter sequences from the ATG start codon of the KoNAC genes were retrieved from the K. obovata genome, and the cis-regulatory elements in the promoter regions were predicted using Plant CARE (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/) (accessed on 31 March 2022) [50] and displayed by TBtools [47].

2.7. Quantitative Real-Time PCR Assays

Total RNA was extracted from the sampled leaves mentioned above using TRIzol (Invitrogen, http://www.invitrogen.com) (accessed on 20 April 2022). Quantitative real-time PCR (qRT-PCR) assays were conducted, as described previously [51], using an ABI PRISM 7500 Real-time PCR System (Applied Biosystem) with 2−ΔΔCT method [52]. The specific primers of KoNAC genes used here are listed in Table S2. The actin gene (GWHTACBH010383.1) was used as an internal control. Student’s t-test in statistical analysis was performed using Graphpad Prism 9.0.0 (https://www.graphpad-prism.cn/) (accessed on 28 April 2022).

3. Results

3.1. Genome-Wide Identification of the K. obovata NAC Genes

Two independent strategies for retrieval of KoNAC genes from the K. obovata genome, HMM search and BLASTP, were used here. Taken together, 68 putative KoNAC genes were identified and confirmed by PfamScan and NCBI-CDD. Based on their chromosome location, these KoNAC genes were named KoNAC1 to KoNAC68 and unevenly distributed on 17 chromosomes (Chrs), with no KoNAC gene present on Chr18 (Table 1, Figure 1). Detail-wise, nine KoNAC genes were located on Chr08, six genes were located on both Chr03 and Chr12, and five KoNAC genes each were located on Chr02, Chr04, Chr06, Chr09, and Chr17, while only one KoNAC gene each was found on Chr07, Chr14, and Chr16.
Moreover, every KoNAC gene contained one or more introns with an average length of 371 bp, while the proteins encoded by KoNAC genes ranged from 162 amino acid (aa) residues (KoNAC61) to 638 aa (KoNAC60) in length, with an average length of 344 aa. The pI values varied from 4.20 (KoNAC45) to 10.07 (KoNAC11), over half of the members (39/68) exhibiting pI > 7 (Table 1).

3.2. Phylogenetic Analysis and Classification of KoNAC Proteins

To illustrate the phylogenetic relationship among K. obovata and Arbidopsis thaliana NAC proteins, a neighbor-joining phylogenetic tree was constructed with 68 KoNAC proteins and 105 AtNAC proteins (Table S1). The result showed that the 173 NAC proteins could be classified into ten classes, namely, Class I to Class X (Figure 2). Obviously, Class VII, with 21 KoNACs and 24 ANACs was the largest class, followed by Class X with 12 KoNACs and 14 ANACs. Other classes contained no more than 10 KoNACs each. Specially, no KoNAC belonged to Class I (Figure 3).
To better understand the phylogenetic relationship and classification of KoNAC genes, the gene structure and motif organization of the 68 KoNAC genes were analyzed. Each KoNAC gene had one or more introns and contained no more than six exons, while over half of the KoNAC genes (41/68) contained three exons (Figure 4c). Additionally, a total of 10 conserved motifs were queried within all K. obovata NAC proteins. Most motifs were located within the N-terminal region (Figure 4b), and motif 1, motif 2, motif 4, and motif 5 were the common elements in KoNAC genes. Clearly, these results showed that the KoNAC genes in the same phylogenetic cluster harbored similar gene structures and motif compositions (Figure 4), which further supported the evolutionary relationship of KoNAC genes demonstrated above.

3.3. Collinearity Analysis of KoNAC Genes

It is well-known that P. trichocarpa is a typical model plant for functional genomics and molecular studies in woody species. Moreover, K. obovata (Rhizophoraceae) and P. trichocarpa (Salicaceae) belong to the same order, Malpighiales (https://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi) (accessed on 18 March 2022). Therefore, to better investigate the evolutionary relationship of NAC genes, the collinearity analysis was performed based on the genomes of A. thaliana and P. trichocarpa (Figure 5). There were 16,355 collinear gene pairs between K. obovata and A. thaliana identified, among which 52 orthologous gene pairs between KoNACs and AtNACs were obtained (Figure 5a, blue lines). Meanwhile, a total of 26,594 collinear gene pairs between K. obovata and P. trichocarpa were available, among which 54 orthologous gene pairs between KoNACs and PtNACs were determined (Figure 5b, purple lines). Taken together, 49 common KoNAC genes shared homologous relationships with both A. thaliana and P. trichocarpa NAC genes (Table S3), implying these genes might function in a similar manner. In the meantime, there were 11 non-orthologous KoNAC genes (KoNAC13, KoNAC16, KoNAC17, KoNAC21, KoNAC40, KoNAC50, KoNAC53, KoNAC63, KoNAC65, KoNAC67, and KoNAC68) compared to NAC genes of A. thaliana and P. trichocarpa. These genes displayed different structure (Figure 4), among which four genes (KoNAC16, KoNAC17, KoNAC21, KoNAC40) clustered as the subgroup of Class II (Figure 2).

3.4. Expression Patterns of KoNAC Genes in Different Organs

To gain an insight into the function of NAC genes in K. obovata, the expression levels of all KoNAC genes in various organs, including root, stem, leaf, flower, pistil, stamen, sepal, and fruit were determined based on previously published RNA-seq data of K. obovata. Noticeably, the expression patterns of KoNAC genes were not in a constitutive mode, whereas they were differentially and preferentially expressed in different organs (Figure 6, Table S4). For example, 28 out of 68 KoNAC genes were highly expressed in roots, while 18 KoNAC genes were preferentially expressed in leaves, and 13 KoNAC genes were mainly expressed in fruits.

3.5. Expression Analysis of KoNAC Genes under Cold Treatment

To gain more insight into the function of KoNAC genes, the expression profiles of these genes under cold treatment were detected based on the public transcriptomic data of K. obovata. There were 13 KoNAC genes differentially expressed in response to chilling stress, among which one down-regulated KoNAC gene (KoNAC51) and 12 up-regulated KoNACs (KoNAC6, KoNAC11, KoNAC15, KoNAC20, KoNAC24, KoNAC26, KoNAC32, KoNAC35, KoNAC38, KoNAC41, KoNAC62, and KoNAC68) were available (Figure 7a, Table S5). Specifically, four different genes, KoNAC6, KoNAC15, KoNAC20, and KoNAC38, were largely upregulated with higher and more significant values after treatment. The expression levels of these four up-regulated and one down-regulated KoNACs were confirmed by qRT-PCR (Figure 7c,d), implying that these KoNAC genes might act as positive or negative regulators in response to chilling stress.

3.6. Stress-Related Cis-Regulatory Elements Identified in KoNAC Genes

To obtain more evidence for the differentially expressed KoNAC genes on stress responses, the cis-regulatory elements in the promoter regions of these 13 KoNAC genes were predicated. Consequently, 8 well-known stress-related elements were available (Figure 7b, Table S6). LTR (low-temperature response; CCGAAA), a core cis-acting element involved in cold stress response, was present in the majority of the detected KoNAC genes. STRE (stress response element; AGGGG) and ARE (antioxidant response element; AAACCA) were two types of regulatory elements in rapid response to anaerobic stress and environmental stimuli. ABRE (abscisic acid response element; ACGTG), ERE (ethylene response element; ATTTTAAA), and TGACG-motifs were responsible for stress induction by three major stress-related hormones, ABA, ethylene, and methyl jasmonate (MeJA), respectively. Additionally, two biotic stress-responsive elements, WRE3 (wound-response element 3; CCACCT) and WUN-motif (wound-responsive element; AAATTACT) were found in several detected promoters, as well.

4. Discussion

The NAC gene family, one of the largest TF families in plants, was used as positive or negative regulators in response to environmental stimuli including cold stress [19]. Multiple investigations on the functional characterization of NAC genes were reported for A. thaliana [5,53,54], P. trichocarpa [14], and other plants [19]. However, the function of the NAC genes in the typical woody mangrove K. obovata responding to abiotic stresses remains largely unknown. Here, we identified a contracted NAC gene family with 68 members from the K. obovata genome. These KoNAC genes were differentially and preferentially expressed in various organs, and 13 KoNAC genes were differentially expressed under cold treatment based on the publicly available RNA-seq data.
KoNAC proteins, in company with A. thaliana NAC proteins were categorized into 10 classes according to phylogenetic analysis. Obviously, the K. obovata NAC family exhibited a significant contraction in number compared to the NAC families in A. thaliana [5], P. trichocarpa [14], and other plants [15,16], and the decreased KoNAC genes in class V, class VIII, and class IX mainly contributed to the contraction (Figure 3). These results are consistent with the previous findings [35], and the contraction might relate to the evolutionary adaption to the intertidal zones. Moreover, 49 genes from the contracted KoNAC family shared orthologous relationships with the NAC genes of A. thaliana and P. trichocarpa, implying these genes might have similar functions [55]. Additionally, compared to AtNAC and KoNAC genes, there existed 11 non-orthologous KoNAC genes, among which, four genes (KoNAC16, KoNAC17, KoNAC21, and KoNAC40) clustered in class II (Figure 2). Another non-orthologous gene, KoNAC68, was induced under cold treatment, implying it might potentially function in response to cold stress in K. obovata (Figure 7a). More attention should be paid to these non-orthologous genes, and functional investigations of these genes will provide valuable knowledge about mangrove species.
To explore the function of NAC genes in K. obovata, the expression patterns of all KoNACs in various organs were determined. In contrast to the constitutive expression patterns of other gene families [56,57,58], KoNAC genes were differentially and preferentially expressed in different organs. For instance, there were 28 KoNAC genes expressed highly in roots, 18 KoNAC genes expressed preferentially in leaves, and 13 KoNAC genes expressed mainly in fruits. Referentially, the root-expressed gene OsNAC2 modulated root development in rice by involving the crosstalk of auxin and cytokinin pathways [59]. The A. thaliana rosette-expressed gene, ANAC087, positively regulated rosette development and leaf senescence [60]. FaRIF, a strawberry NAC gene, was reported as one key regulator controlling fruit ripening [61]. This evidence implied that these organ-specific KoNAC genes might function as key regulators in organ development. Moreover, the organ expression patterns of NAC genes in K. obovata were not similar to that in A. thaliana and other plants. For instance, KoNAC46 and KoNAC54 were primarily expressed in roots (Figure 6), however, ANAC048 and ANAC074, the closest orthologues of these two KoNAC genes (Figure 2), respectively, were expressed in different organs. ANAC048 was involved in vascular development [62], and ANAC074 positively regulated programmed cell death of stigmatic tissue in A. thaliana [63]. Therefore, functional characterization of the KoNAC genes primarily expressed in roots distinct from other plants should be deeply covered in the future.
To better understand the roles of KoNAC genes, the expression analysis under chilling stress was performed based on the public transcriptomic data. In total, 13 out of 68 KoNAC genes were differentially expressed under cold treatment. Among them, KoNAC51 was the only down-regulated gene, whereas its closest homologue KoNAC35 was up-regulated after treatment, implying these two class V genes might function oppositely in response to cold stress. Half of the up-regulated genes (KoNAC6, KoNAC15, KoNAC20, KoNAC26, KoNAC32, and KoNAC38) belonged to the class VII subgroup (Figure 2 and Figure 7a). Particularly, KoNAC6, KoNAC15, and KoNAC26 clustered together and shared high sequence similarity to their closest orthologs ANAC002, ANAC081, and ANAC102 in A. thaliana. ANAC002 (ATAF1) was reported to serve as dual regulators responsive to abiotic and biotic stresses [64,65,66]. ANAC081 (ATAF2) was rapidly induced by pathogen attack and involved in plant defense [67,68], while ANAC102 was responsive to low-oxygen and high-light stresses [69,70]. Overexpression of MlNAC5, another closest ortholog of KoNAC26 and ANAC002 from Miscanthus lutarioriparius L.Liu ex S.L.Chen & Renvoize, led to enhanced tolerance to cold and drought stresses in A. thaliana [71]. Additionally, three closest orthologs of KoNAC32, ANAC019, ANAC055 and ANAC072, were required for drought tolerance in A. thaliana [72], among which ANAC019 and ANAC055 displayed a dual function in regulating ABA response and jasmonate response [73,74]. A. thaliana ANAC042, the closest ortholog of KoNAC38, conferred stress tolerance through regulating phytohormone metabolism and signaling [75,76,77]. Moreover, various stress-related cis-regulatory elements were identified from the promoters of these KoNAC genes (Figure 7b). The LTR element is an indispensable cis-acting element in plant response to low temperature [78,79]. Deletion of the LTR element will result in complete loss of promoter activity under cold stress [79]. STRE is a common cis-regulatory element in eukaryotes, and involved in response to multiple environmental stimuli [80]. ARE is an antioxidant response element in rapid response to anaerobic stress [81]. Meanwhile, ABRE, ERE, and TGACG motifs are three major types of elements related to plant hormones (ABA, ethylene, and MeJA) [82,83,84]. Among them, ABRE and TGACG motifs are enriched in the majority of the detected KoNAC genes, implying these KoNAC genes might respond to stresses via hormone-mediated pathways [85]. Additionally, both WRE3 and WUN motifs are biotic stress-responsive elements and present in several detected promoters as well, implying the KoNAC genes might function in response to biotic stresses [86,87]. Taken together, these findings suggest that these KoNAC genes may be involved responses to other abiotic or biotic stresses in addition to the cold response, providing auxiliary evidence for these KoNAC genes in response to abiotic and biotic stresses. To know more about the function of KoNAC genes, further investigation and more proof are required.

5. Conclusions

In the present study, a pluri-disciplinary work concerning comprehensive analysis of the KoNAC gene family was performed. We identified a contracted NAC TF family containing 68 genes from the genome of the typical mangrove plant K. obovata based on bioinformatic analysis. These KoNAC genes were unevenly distributed in 17 chromosomes of K. obovata. The NAC genes of K. obovata and A. thaliana were classified into ten classes, while no KoNAC gene belonged to class I. Obviously, the decreased members of class V, class VIII, and class IX mainly contributed to the contraction of the KoNAC family. KoNAC genes were differentially and preferentially expressed in different organs. Among them, 13 KoNAC genes were rapidly induced by chilling stress. The expression patterns of five selected KoNAC genes (KoNAC6, KoNAC15, KoNAC20, KoNAC38, and KoNAC51) were confirmed by qRT-PCR. Additionally, several stress-related cis-acting elements were detected in the promoter regions of these KoNAC genes, implying KoNAC genes might participate in multiple stress responses. Summarily, our findings will provide positive references for further investigations on functional characterization of KoNAC genes in stress responses.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb44110381/s1, Figure S1: Morphological features of K. obovata; Table S1: The basic information of A. thaliana NAC genes; Table S2: The primers of KoNAC genes in this study; Table S3: Common orthologous gene pairs of K. obovata between A. thaliana and P. trichocarpa; Table S4: The FPKM values of KoNAC genes expressed in different organs; Table S5: The FPKM values of differentially expressed KoNAC genes under cold treatment; Table S6: Number of stress-related cis-regulatory elements in the promoter regions of the differentially expressed KoNAC genes.

Author Contributions

Conceptualization, M.-M.S. and H.-Z.L.; formal analysis, M.-M.S., X.L. and H.-Z.L.; funding acquisition, X.L. and H.-Z.L.; investigation, M.-M.S., X.L., X.-J.H., J.-J.Y., P.-T.Q., H.Z. and M.-G.J.; methodology, M.-M.S. and H.-Z.L.; resources, X.L. and H.-Z.L.; supervision, H.-Z.L.; writing—original draft, M.-M.S.; writing—review and editing, H.-Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the grants from Guangxi Minzu University Research Fund (No. 2021KJQD18; funder: Guangxi Minzu University), Guangxi Forestry Science and Technology Promotion Demonstration Project (No. 2021-01; funder: Department of Forestry of Guangxi Zhuang Autonmous Region) and Guangxi Key Laboratory of Special Non-wood Forest Cultivation and Utilization (No. JB-22-02-02; funder: Guangxi Forestry Research Institute).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The chromosome-scale genome of Kandelia obovata was obtained from Genome Warehouse under accession number GWHACBH00000000 (https://bigd.big.ac.cn/gwh) (accessed on 8 March 2022). The Arabidopsis thaliana genome was obtained from TAIR database (http://www.arabidopsis.org/) (accessed on 10 March 2022). The Populus trichocarpa genome version 2.0 was downloaded from Phytozome (http://www.phytozome.net/poplar) (accessed on 18 March 2022). The RNA-seq data in eight different organs of K. obovata was available from NCBI BioProject under accession number PRJNA416402 (https://www.ncbi.nlm.nih.gov/bioproject) (accessed on 31 March 2022). The RNA-seq data related to cold treatment in K. obovata was obtained from NCBI BioProject under accession number PRJNA678025.

Acknowledgments

We sincerely thank Hong-Hui Cui for positive suggestions and comments on this work. Yan Yan and Chang-Yuan Liao are acknowledged for their warm support.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Singh, K.; Foley, R.C.; Onate-Sanchez, L. Transcription factors in plant defense and stress responses. Curr. Opin. Plant Biol. 2002, 5, 430–436. [Google Scholar] [CrossRef]
  2. Souer, E.; van Houwelingen, A.; Kloos, D.; Mol, J.; Koes, R. The No Apical Meristem gene of Petunia is required for pattern formation in embryos and flowers and is expressed at meristem and primordia boundaries. Cell 1996, 85, 159–170. [Google Scholar] [CrossRef] [Green Version]
  3. Aida, M.; Ishida, T.; Fukaki, H.; Fujisawa, H.; Tasaka, M. Genes involved in organ separation in Arabidopsis: An analysis of the Cup-Shaped cotyledon mutant. Plant Cell 1997, 9, 841–857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Olsen, A.N.; Ernst, H.A.; Leggio, L.L.; Skriver, K. NAC transcription factors: Structurally distinct, functionally diverse. Trends Plant Sci. 2005, 10, 79–87. [Google Scholar] [CrossRef]
  5. Ooka, H.; Satoh, K.; Doi, K.; Nagata, T.; Otomo, Y.; Murakami, K.; Matsubara, K.; Osato, N.; Kawai, J.; Carninci, P.; et al. Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana. DNA Res. 2003, 10, 239–247. [Google Scholar] [CrossRef] [PubMed]
  6. Jia, D.; Jiang, Z.; Fu, H.; Chen, L.; Liao, G.; He, Y.; Huang, C.; Xu, X. Genome-wide identification and comprehensive analysis of NAC family genes involved in fruit development in kiwifruit (Actinidia). BMC Plant Biol. 2021, 21, 44. [Google Scholar] [CrossRef] [PubMed]
  7. Li, C.; Zhang, J.; Zhang, Q.; Dong, A.; Wu, Q.; Zhu, X.; Zhu, X. Genome-wide identification and analysis of the NAC transcription factor gene family in garden Asparagus (Asparagus officinalis). Genes 2022, 13, 976. [Google Scholar] [CrossRef]
  8. Chen, S.; Lin, X.; Zhang, D.; Li, Q.; Zhao, X.; Chen, S. Genome-wide analysis of nac gene family in Betula pendula. Forests 2019, 10, 741. [Google Scholar] [CrossRef] [Green Version]
  9. Hu, X.; Xie, F.; Liang, W.; Liang, Y.; Zhang, Z.; Zhao, J.; Hu, G.; Qin, Y. HuNAC20 and HuNAC25, two novel NAC genes from pitaya, confer cold tolerance in transgenic Arabidopsis. Int. J. Mol. Sci. 2022, 23, 2189. [Google Scholar] [CrossRef]
  10. Li, X.; Cai, K.; Pei, X.; Li, Y.; Hu, Y.; Meng, F.; Song, X.; Tigabu, M.; Ding, C.; Zhao, X. Genome-wide identification of NAC transcription factor family in Juglans mandshurica and their expression analysis during the fruit development and ripening. Int. J. Mol. Sci. 2021, 22, 12414. [Google Scholar] [CrossRef]
  11. He, F.; Zhang, L.; Zhao, G.; Kang, J.; Long, R.; Li, M.; Yang, Q.; Chen, L. Genome-wide identification and expression analysis of the NAC gene family in Alfalfa revealed its potential roles in response to multiple abiotic stresses. Int. J. Mol. Sci. 2022, 23, 10015. [Google Scholar] [CrossRef] [PubMed]
  12. Nie, G.; Yang, Z.; He, J.; Liu, A.; Chen, J.; Wang, S.; Wang, X.; Feng, G.; Li, D.; Peng, Y.; et al. Genome-wide investigation of the NAC transcription factor family in Miscanthus sinensis and expression analysis under various abiotic stress. Front. Plant Sci. 2021, 12, 766550. [Google Scholar] [CrossRef] [PubMed]
  13. Li, B.; Fan, R.; Yang, Q.; Hu, C.; Sheng, O.; Deng, G.; Dong, T.; Li, C.; Peng, X.; Bi, F.; et al. Genome-wide identification and characterization of the NAC transcription factor family in Musa Acuminata and expression analysis during fruit ripening. Int. J. Mol. Sci. 2020, 21, 634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hu, R.; Qi, G.; Kong, Y.; Kong, D.; Gao, Q.; Zhou, G. Comprehensive analysis of NAC domain transcription factor gene family in Populus trichocarpa. BMC Plant Biol. 2010, 10, 145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Yang, H.; Fan, L.; Yu, X.; Zhang, X.; Hao, P.; Wei, D.; Zhang, G. Analysis of the NAC gene family in Salix and the identification of SpsNAC005 gene contributing to salt and drought tolerance. Forests 2022, 13, 971. [Google Scholar] [CrossRef]
  16. Jin, J.F.; Wang, Z.Q.; He, Q.Y.; Wang, J.Y.; Li, P.F.; Xu, J.M.; Zheng, S.J.; Fan, W.; Yang, J.L. Genome-wide identification and expression analysis of the NAC transcription factor family in tomato (Solanum lycopersicum) during aluminum stress. BMC Genom. 2020, 21, 288. [Google Scholar] [CrossRef] [Green Version]
  17. Hu, H.; Ma, L.; Chen, X.; Fei, X.; He, B.; Luo, Y.; Liu, Y.; Wei, A. Genome-wide identification of the NAC gene family in Zanthoxylum bungeanum and their transcriptional responses to drought stress. Int. J. Mol. Sci. 2022, 23, 4769. [Google Scholar] [CrossRef]
  18. Peng, X.; Zhao, Y.; Li, X.; Wu, M.; Chai, W.; Sheng, L.; Wang, Y.; Dong, Q.; Jiang, H.; Cheng, B. Genome-wide identification, classification and analysis of NAC type gene family in maize. J. Genet. 2015, 94, 377–390. [Google Scholar] [CrossRef]
  19. Diao, P.; Chen, C.; Zhang, Y.; Meng, Q.; Lv, W.; Ma, N. The role of NAC transcription factor in plant cold response. Plant Signal. Behav. 2020, 15, 1785668. [Google Scholar] [CrossRef]
  20. Puranik, S.; Sahu, P.P.; Srivastava, P.S.; Prasad, M. NAC proteins: Regulation and role in stress tolerance. Trends Plant Sci. 2012, 17, 369–381. [Google Scholar] [CrossRef]
  21. Shan, W.; Kuang, J.F.; Lu, W.J.; Chen, J.Y. Banana fruit NAC transcription factor MaNAC1 is a direct target of MaICE1 and involved in cold stress through interacting with MaCBF1. Plant Cell Environ. 2014, 37, 2116–2127. [Google Scholar] [CrossRef] [PubMed]
  22. Han, D.; Du, M.; Zhou, Z.; Wang, S.; Li, T.; Han, J.; Xu, T.; Yang, G. Overexpression of a Malus baccata NAC transcription factor gene MbNAC25 increases cold and salinity tolerance in Arabidopsis. Int. J. Mol. Sci. 2020, 21, 1198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Dong, Y.; Tang, M.; Huang, Z.; Song, J.; Xu, J.; Ahammed, G.J.; Yu, J.; Zhou, Y. The miR164a-NAM3 module confers cold tolerance by inducing ethylene production in tomato. Plant J. 2022, 111, 440–456. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, Z.; Zhang, Y.; Hu, H.; Chen, L.; Zhang, H.; Chen, R. CabHLH79 acts upstream of CaNAC035 to regulate cold stress in pepper. Int. J. Mol. Sci. 2022, 23, 2537. [Google Scholar] [CrossRef]
  25. Hou, X.M.; Zhang, H.F.; Liu, S.Y.; Wang, X.K.; Zhang, Y.M.; Meng, Y.C.; Luo, D.; Chen, R.G. The NAC transcription factor CaNAC064 is a regulator of cold stress tolerance in peppers. Plant Sci. 2020, 291, 110346. [Google Scholar] [CrossRef]
  26. An, J.P.; Li, R.; Qu, F.J.; You, C.X.; Wang, X.F.; Hao, Y.J. An apple NAC transcription factor negatively regulates cold tolerance via CBF-dependent pathway. J. Plant Physiol. 2018, 221, 74–80. [Google Scholar] [CrossRef]
  27. Song, C.; Wu, M.; Zhou, Y.; Gong, Z.; Yu, W.; Zhang, Y.; Yang, Z. NAC-mediated membrane lipid remodeling negatively regulates fruit cold tolerance. Hortic. Res. 2022, 9, uhac039. [Google Scholar] [CrossRef]
  28. Krauss, K.W.; McKee, K.L.; Lovelock, C.E.; Cahoon, D.R.; Saintilan, N.; Reef, R.; Chen, L. How mangrove forests adjust to rising sea level. New Phytol. 2014, 202, 19–34. [Google Scholar] [CrossRef] [Green Version]
  29. Nizam, A.; Meera, S.P.; Kumar, A. Genetic and molecular mechanisms underlying mangrove adaptations to intertidal environments. iScience 2022, 25, 103547. [Google Scholar] [CrossRef]
  30. Guo, Z.; Ma, D.; Li, J.; Wei, M.; Zhang, L.; Zhou, L.; Zhou, X.; He, S.; Wang, L.; Shen, Y.; et al. Genome-wide identification and characterization of aquaporins in mangrove plant Kandelia obovata and its role in response to the intertidal environment. Plant Cell Environ. 2022, 45, 1698–1718. [Google Scholar] [CrossRef]
  31. Su, W.; Ye, C.; Zhang, Y.; Hao, S.; Li, Q.Q. Identification of putative key genes for coastal environments and cold adaptation in mangrove Kandelia obovata through transcriptome analysis. Sci. Total Environ. 2019, 681, 191–201. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, L.; Wang, W.; Li, Q.Q.; Zhang, Y.; Yang, S.; Osland, M.J.; Huang, J.; Peng, C. Mangrove species’ responses to winter air temperature extremes in China. Ecosphere 2017, 8, e01865. [Google Scholar] [CrossRef]
  33. Peng, Y.-L.; Wang, Y.-S.; Fei, J.; Sun, C.-C.; Cheng, H. Ecophysiological differences between three mangrove seedlings (Kandelia obovata, Aegiceras corniculatum, and Avicennia marina) exposed to chilling stress. Ecotoxicology 2015, 24, 1722–1732. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, S.-M.; Wang, Y.-S.; Su, B.-Y.; Zhou, Y.-Y.; Chang, L.-F.; Ma, X.-Y.; Li, X.-M. Ecophysiological responses of five mangrove species (Bruguiera gymnorrhiza, Rhizophora stylosa, Aegiceras corniculatum, Avicennia marina, and Kandelia obovata) to chilling stress. Front. Mar. Sci. 2022, 9, 846566. [Google Scholar] [CrossRef]
  35. Hu, M.J.; Sun, W.H.; Tsai, W.C.; Xiang, S.; Lai, X.K.; Chen, D.Q.; Liu, X.D.; Wang, Y.F.; Le, Y.X.; Chen, S.M.; et al. Chromosome-scale assembly of the Kandelia obovata genome. Hortic. Res. 2020, 7, 75. [Google Scholar] [CrossRef] [PubMed]
  36. Mistry, J.; Chuguransky, S.; Williams, L.; Qureshi, M.; Salazar, G.A.; Sonnhammer, E.L.L.; Tosatto, S.C.E.; Paladin, L.; Raj, S.; Richardson, L.J.; et al. Pfam: The protein families database in 2021. Nucleic Acids Res. 2021, 49, D412–D419. [Google Scholar] [CrossRef]
  37. Waack, S.; Keller, O.; Asper, R.; Brodag, T.; Damm, C.; Fricke, W.F.; Surovcik, K.; Meinicke, P.; Merkl, R. Score-based prediction of genomic islands in prokaryotic genomes using hidden Markov models. BMC Bioinform. 2006, 7, 142. [Google Scholar] [CrossRef] [Green Version]
  38. Lamesch, P.; Berardini, T.Z.; Li, D.; Swarbreck, D.; Wilks, C.; Sasidharan, R.; Muller, R.; Dreher, K.; Alexander, D.L.; Garcia-Hernandez, M.; et al. The Arabidopsis Information Resource (TAIR): Improved gene annotation and new tools. Nucleic Acids Res. 2012, 40, D1202–D1210. [Google Scholar] [CrossRef]
  39. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef] [Green Version]
  40. 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.; et al. CDD/SPARCLE: The conserved domain database in 2020. Nucleic Acids Res. 2020, 48, D265–D268. [Google Scholar] [CrossRef]
  41. Voorrips, R.E. MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef] [PubMed]
  44. Hu, B.; Jin, J.; Guo, A.Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  46. Wang, Y.; Tang, H.; Debarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant. 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  48. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [Green Version]
  49. Du, Z.; You, S.; Zhao, X.; Xiong, L.; Li, J. Genome-wide identification of WRKY genes and their responses to chilling stress in Kandelia obovata. Front. Genet. 2022, 13, 875316. [Google Scholar] [CrossRef]
  50. Lescot, M.; Dehais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouze, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
  51. Liao, H.Z.; Zhu, M.M.; Cui, H.H.; Du, X.Y.; Tang, Y.; Chen, L.Q.; Ye, D.; Zhang, X.Q. MARIS plays important roles in Arabidopsis pollen tube and root hair growth. J. Integr. Plant Biol. 2016, 58, 927–940. [Google Scholar] [CrossRef] [PubMed]
  52. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  53. Guan, Q.; Yue, X.; Zeng, H.; Zhu, J. The protein phosphatase RCF2 and its interacting partner NAC019 are critical for heat stress-responsive gene regulation and thermotolerance in Arabidopsis. Plant Cell 2014, 26, 438–453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Yoo, S.Y.; Kim, Y.; Kim, S.Y.; Lee, J.S.; Ahn, J.H. Control of flowering time and cold response by a NAC-domain protein in Arabidopsis. PLoS ONE 2007, 2, e642. [Google Scholar] [CrossRef] [PubMed]
  55. Gabaldon, T.; Koonin, E.V. Functional and evolutionary implications of gene orthology. Nat. Rev. Genet. 2013, 14, 360–366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Han, N.; Tang, R.; Chen, X.; Xu, Z.; Ren, Z.; Wang, L. Genome-wide identification and characterization of WOX genes in Cucumis sativus. Genome 2021, 64, 761–776. [Google Scholar] [CrossRef] [PubMed]
  57. Rahman, H.; Xu, Y.P.; Zhang, X.R.; Cai, X.Z. Brassica napus genome possesses extraordinary high number of CAMTA genes and CAMTA3 contributes to PAMP triggered immunity and resistance to Sclerotinia sclerotiorum. Front. Plant Sci. 2016, 7, 581. [Google Scholar] [CrossRef] [Green Version]
  58. Yi, H.; Sardesai, N.; Fujinuma, T.; Chan, C.W.; Veena; Gelvin, S. B. Constitutive expression exposes functional redundancy between the Arabidopsis histone H2A gene HTA1 and other H2A gene family members. Plant Cell 2006, 18, 1575–1589. [Google Scholar] [CrossRef] [Green Version]
  59. Mao, C.; He, J.; Liu, L.; Deng, Q.; Yao, X.; Liu, C.; Qiao, Y.; Li, P.; Ming, F. OsNAC2 integrates auxin and cytokinin pathways to modulate rice root development. Plant Biotechnol. J. 2020, 18, 429–442. [Google Scholar] [CrossRef] [Green Version]
  60. Vargas-Hernandez, B.Y.; Nunez-Munoz, L.; Calderon-Perez, B.; Xoconostle-Cazares, B.; Ruiz-Medrano, R. The NAC transcription factor ANAC087 induces aerial rosette development and leaf senescence in Arabidopsis. Front. Plant Sci. 2022, 13, 818107. [Google Scholar] [CrossRef]
  61. Martín-Pizarro, C.; Vallarino, J.G.; Osorio, S.; Meco, V.; Urrutia, M.; Pillet, J.; Casañal, A.; Merchante, C.; Amaya, I.; Willmitzer, L.; et al. The NAC transcription factor FaRIF controls fruit ripening in strawberry. Plant Cell 2021, 33, 1574–1593. [Google Scholar] [CrossRef] [PubMed]
  62. Yang, J.H.; Lee, K.-H.; Du, Q.; Yang, S.; Yuan, B.; Qi, L.; Wang, H. A membrane-associated NAC domain transcription factor XVP interacts with TDIF co-receptor and regulates vascular meristem activity. New Phytol. 2020, 226, 59–74. [Google Scholar] [CrossRef] [PubMed]
  63. Gao, Z.; Daneva, A.; Salanenka, Y.; Van Durme, M.; Huysmans, M.; Lin, Z.; De Winter, F.; Vanneste, S.; Karimi, M.; Van de Velde, J.; et al. KIRA1 and ORESARA1 terminate flower receptivity by promoting cell death in the stigma of Arabidopsis. Nat. Plants 2018, 4, 365–375. [Google Scholar] [CrossRef] [PubMed]
  64. Alshareef, N.O.; Otterbach, S.L.; Allu, A.D.; Woo, Y.H.; de Werk, T.; Kamranfar, I.; Mueller-Roeber, B.; Tester, M.; Balazadeh, S.; Schmockel, S.M. NAC transcription factors ATAF1 and ANAC055 affect the heat stress response in Arabidopsis. Sci. Rep. 2022, 12, 11264. [Google Scholar] [CrossRef]
  65. Liu, Y.; Sun, J.; Wu, Y. Arabidopsis ATAF1 enhances the tolerance to salt stress and ABA in transgenic rice. J. Plant Res. 2016, 129, 955–962. [Google Scholar] [CrossRef] [PubMed]
  66. Wu, Y.; Deng, Z.; Lai, J.; Zhang, Y.; Yang, C.; Yin, B.; Zhao, Q.; Zhang, L.; Li, Y.; Yang, C.; et al. Dual function of Arabidopsis ATAF1 in abiotic and biotic stress responses. Cell Res. 2009, 19, 1279–1290. [Google Scholar] [CrossRef] [Green Version]
  67. Delessert, C.; Kazan, K.; Wilson, I.W.; Van Der Straeten, D.; Manners, J.; Dennis, E.S.; Dolferus, R. The transcription factor ATAF2 represses the expression of pathogenesis-related genes in Arabidopsis. Plant J. 2005, 43, 745–757. [Google Scholar] [CrossRef]
  68. Wang, X.; Culver, J.N. DNA binding specificity of ATAF2, a NAC domain transcription factor targeted for degradation by Tobacco mosaic virus. BMC Plant Biol. 2012, 12, 157. [Google Scholar] [CrossRef] [Green Version]
  69. Christianson, J.A.; Wilson, I.W.; Llewellyn, D.J.; Dennis, E.S. The low-oxygen-induced NAC domain transcription factor ANAC102 affects viability of Arabidopsis seeds following low-oxygen treatment. Plant Physiol. 2009, 149, 1724–1738. [Google Scholar] [CrossRef] [Green Version]
  70. D’Alessandro, S.; Ksas, B.; Havaux, M. Decoding β-cyclocitral-mediated retrograde signaling reveals the role of a detoxification response in plant tolerance to photooxidative stress. Plant Cell 2018, 30, 2495–2511. [Google Scholar] [CrossRef]
  71. Yang, X.; Wang, X.; Ji, L.; Yi, Z.; Fu, C.; Ran, J.; Hu, R.; Zhou, G. Overexpression of a Miscanthus lutarioriparius NAC gene MlNAC5 confers enhanced drought and cold tolerance in Arabidopsis. Plant Cell Rep. 2015, 34, 943–958. [Google Scholar] [CrossRef] [PubMed]
  72. Tran, L.S.P.; Nakashima, K.; Sakuma, Y.; Simpson, S.D.; Fujita, Y.; Maruyama, K.; Fujita, M.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Isolation and functional analysis of Arabidopsis stress-inducible NAC transcription factors that bind to a drought-responsive cis-element in the early responsive to dehydration stress 1 promoter. Plant Cell 2004, 16, 2481–2498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Bu, Q.; Jiang, H.; Li, C.B.; Zhai, Q.; Zhang, J.; Wu, X.; Sun, J.; Xie, Q.; Li, C. Role of the Arabidopsis thaliana NAC transcription factors ANAC019 and ANAC055 in regulating jasmonic acid-signaled defense responses. Cell Res. 2008, 18, 756–767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Jiang, H.; Li, H.; Bu, Q.; Li, C. The RHA2a-interacting proteins ANAC019 and ANAC055 may play a dual role in regulating ABA response and jasmonate response. Plant Signal. Behav. 2009, 4, 464–466. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Ebrahimian-Motlagh, S.; Ribone, P.A.; Thirumalaikumar, V.P.; Allu, A.D.; Chan, R.L.; Mueller-Roeber, B.; Balazadeh, S. JUNGBRUNNEN1 confers drought tolerance downstream of the HD-Zip I transcription factor AtHB13. Front. Plant Sci. 2017, 8, 2118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Sakuraba, Y.; Bulbul, S.; Piao, W.; Choi, G.; Paek, N.C. Arabidopsis EARLY FLOWERING3 increases salt tolerance by suppressing salt stress response pathways. Plant J. 2017, 92, 1106–1120. [Google Scholar] [CrossRef] [Green Version]
  77. Shahnejat-Bushehri, S.; Tarkowska, D.; Sakuraba, Y.; Balazadeh, S. Arabidopsis NAC transcription factor JUB1 regulates GA/BR metabolism and signalling. Nat. Plants 2016, 2, 16013. [Google Scholar] [CrossRef]
  78. Qian, W.; Xiao, B.; Wang, L.; Hao, X.; Yue, C.; Cao, H.; Wang, Y.; Li, N.; Yu, Y.; Zeng, J.; et al. CsINV5, a tea vacuolar invertase gene enhances cold tolerance in transgenic Arabidopsis. BMC Plant Biol. 2018, 18, 228. [Google Scholar] [CrossRef]
  79. Wu, C.; Zheng, C.; Ji, G.; Jiang, P. Synergistic effects of HSE and LTR elements from hsp70 gene promoter of Ulva prolifera (Ulvophyceae, Chlorophyta) upon temperature induction. J. Phycol. 2019, 55, 738–743. [Google Scholar] [CrossRef]
  80. Ebeed, H.T. Genome-wide analysis of polyamine biosynthesis genes in wheat reveals gene expression specificity and involvement of STRE and MYB-elements in regulating polyamines under drought. BMC Genom. 2022, 23, 734. [Google Scholar] [CrossRef]
  81. Fan, S.; Liu, A.; Zhang, Z.; Zou, X.; Jiang, X.; Huang, J.; Fan, L.; Zhang, Z.; Deng, X.; Ge, Q.; et al. Genome-wide identification and expression analysis of the metacaspase gene family in Gossypium species. Genes 2019, 10, 527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Watanabe, K.A.; Homayouni, A.; Gu, L.; Huang, K.Y.; Ho, T.D.; Shen, Q.J. Transcriptomic analysis of rice aleurone cells identified a novel abscisic acid response element. Plant Cell Environ. 2017, 40, 2004–2016. [Google Scholar] [CrossRef] [PubMed]
  83. Itzhaki, H.; Maxson, J.M.; Woodson, W.R. An ethylene-responsive enhancer element is involved in the senescence-related expression of the carnation glutathione-S-transferase (GST1) gene. Proc. Natl. Acad. Sci. USA 1994, 91, 8925–8929. [Google Scholar] [CrossRef] [Green Version]
  84. Huo, Y.; Zhang, B.; Chen, L.; Zhang, J.; Zhang, X.; Zhu, C. Isolation and functional characterization of the promoters of miltiradiene synthase genes, TwTPS27a and TwTPS27b, and interaction analysis with the transcription factor TwTGA1 from Tripterygium wilfordii. Plants 2021, 10, 418. [Google Scholar] [CrossRef]
  85. Verma, V.; Ravindran, P.; Kumar, P.P. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016, 16, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. He, C.; Liu, X.; Teixeira da Silva, J.A.; Wang, H.; Peng, T.; Zhang, M.; Si, C.; Yu, Z.; Tan, J.; Zhang, J.; et al. Characterization of LEA genes in Dendrobium officinale and one gene in induction of callus. J. Plant Physiol. 2021, 258–259, 153356. [Google Scholar] [CrossRef] [PubMed]
  87. Tanin, M.J.; Saini, D.K.; Sandhu, K.S.; Pal, N.; Gudi, S.; Chaudhary, J.; Sharma, A. Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding. Sci. Rep. 2022, 12, 13680. [Google Scholar] [CrossRef]
Cimb 44 00381 g001 550
Figure 1. Distribution of KoNAC genes in K. obovata genome. The 68 KoNAC genes were unevenly distributed on 17 chromosomes (Chr01–Chr17) denoted in different colors, while no KoNAC gene was found on Chr18. Values on the y-axis indicate the chromosome length and gene position.
Figure 1. Distribution of KoNAC genes in K. obovata genome. The 68 KoNAC genes were unevenly distributed on 17 chromosomes (Chr01–Chr17) denoted in different colors, while no KoNAC gene was found on Chr18. Values on the y-axis indicate the chromosome length and gene position.
Cimb 44 00381 g001
Cimb 44 00381 g002 550
Figure 2. Phylogenetic analysis of the NAC proteins from K. obovata and A. thaliana. The unrooted tree was constructed by MEGA 7.0 based on the neighbor-joining (NJ) method with 1000 bootstrap replications. The NAC proteins of K. obovata (KoNAC) in blue squares, and A. thaliana (AtNAC) in red circles were classified into ten classes represented with different colors.
Figure 2. Phylogenetic analysis of the NAC proteins from K. obovata and A. thaliana. The unrooted tree was constructed by MEGA 7.0 based on the neighbor-joining (NJ) method with 1000 bootstrap replications. The NAC proteins of K. obovata (KoNAC) in blue squares, and A. thaliana (AtNAC) in red circles were classified into ten classes represented with different colors.
Cimb 44 00381 g002
Cimb 44 00381 g003 550
Figure 3. Classification comparison between K. obovata (Ko) and A. thaliana (At) NAC proteins.
Figure 3. Classification comparison between K. obovata (Ko) and A. thaliana (At) NAC proteins.
Cimb 44 00381 g003
Cimb 44 00381 g004 550
Figure 4. Gene structure and motif organization of the KoNAC genes. (a) The unrooted NJ tree of KoNAC proteins was constructed using MEGA 7.0, while different classes were represented in different colors. Bootstrap values from 1000 replicates are shown on the nodes. (b) Ten different conserved motifs of KoNAC proteins were identified using MEME software. Different colored boxes indicate different motifs. (c) The structures of KoNAC genes are shown, including UTR regions (blue box), exons (grey boxes), and introns (black lines). UTR, untranslated region; CDS, coding sequence of a gene.
Figure 4. Gene structure and motif organization of the KoNAC genes. (a) The unrooted NJ tree of KoNAC proteins was constructed using MEGA 7.0, while different classes were represented in different colors. Bootstrap values from 1000 replicates are shown on the nodes. (b) Ten different conserved motifs of KoNAC proteins were identified using MEME software. Different colored boxes indicate different motifs. (c) The structures of KoNAC genes are shown, including UTR regions (blue box), exons (grey boxes), and introns (black lines). UTR, untranslated region; CDS, coding sequence of a gene.
Cimb 44 00381 g004
Cimb 44 00381 g005 550
Figure 5. Collinearity analysis of KoNAC genes. (a) Genes in 5 chromosomes of A. thaliana (AtChrs), in orange, and 18 chromosomes of K. obovata (KoChrs), in cyan, are introduced here. The orthologous pairs between KoNACs and ANACs are highlighted in blue. The gene pairs among ANACs are colored in yellow, while the gene pairs among KoNACs are colored in red. (b) Genes in 19 chromosomes of P. trichocarpa (PtChrs), in pink, and 18 chromosomes of K. obovata (KoChrs), in cyan, are introduced here. The orthologous pairs between KoNACs and PtNACs are highlighted in purple. The gene pairs among PtNACs are colored in green, while the gene pairs among KoNACs are colored in red.
Figure 5. Collinearity analysis of KoNAC genes. (a) Genes in 5 chromosomes of A. thaliana (AtChrs), in orange, and 18 chromosomes of K. obovata (KoChrs), in cyan, are introduced here. The orthologous pairs between KoNACs and ANACs are highlighted in blue. The gene pairs among ANACs are colored in yellow, while the gene pairs among KoNACs are colored in red. (b) Genes in 19 chromosomes of P. trichocarpa (PtChrs), in pink, and 18 chromosomes of K. obovata (KoChrs), in cyan, are introduced here. The orthologous pairs between KoNACs and PtNACs are highlighted in purple. The gene pairs among PtNACs are colored in green, while the gene pairs among KoNACs are colored in red.
Cimb 44 00381 g005
Cimb 44 00381 g006 550
Figure 6. Expression patterns of KoNAC genes in different organs. The transcript levels of the KoNAC genes in eight organs of K. obovata were determined based on published transcriptomic data (NCBI BioProject: PRJNA416402). The color scale indicates increasing expression levels from blue to red. Deeper red colors represent higher expression levels, while darker blue colors indicate lower values.
Figure 6. Expression patterns of KoNAC genes in different organs. The transcript levels of the KoNAC genes in eight organs of K. obovata were determined based on published transcriptomic data (NCBI BioProject: PRJNA416402). The color scale indicates increasing expression levels from blue to red. Deeper red colors represent higher expression levels, while darker blue colors indicate lower values.
Cimb 44 00381 g006
Cimb 44 00381 g007 550
Figure 7. Expression analysis of KoNAC genes under cold treatment. (a) The transcript levels of the KoNAC genes in response to cold were determined based on publicly available RNA-seq data (NCBI BioProject: PRJNA678025). Deeper red colors represent higher expression levels of up-regulated KoNAC genes, while darker blue colors indicate higher values of down-regulated KoNAC genes. Cold 1, first-time cold treatment; Cold 2, second-time cold treatment; Cold 4, fourth-time cold treatment. (b) The cis-regulatory elements in the promoters of the 13 KoNAC genes were predicated by PlantCARE. Eight well-known stress-related elements were identified. The size of the blue ball indicates the number of the elements in the KoNAC promoters. Expression levels of four upregulated KoNAC genes (c) and one down-regulated KoNAC (d) under cold treatment were confirmed by qRT-PCR. Three independent experiments were performed. The actin gene in K. obovata acted as the internal control. Asterisks indicate significant differences compared with CK by Student’s t-test. *, p < 0.05.
Figure 7. Expression analysis of KoNAC genes under cold treatment. (a) The transcript levels of the KoNAC genes in response to cold were determined based on publicly available RNA-seq data (NCBI BioProject: PRJNA678025). Deeper red colors represent higher expression levels of up-regulated KoNAC genes, while darker blue colors indicate higher values of down-regulated KoNAC genes. Cold 1, first-time cold treatment; Cold 2, second-time cold treatment; Cold 4, fourth-time cold treatment. (b) The cis-regulatory elements in the promoters of the 13 KoNAC genes were predicated by PlantCARE. Eight well-known stress-related elements were identified. The size of the blue ball indicates the number of the elements in the KoNAC promoters. Expression levels of four upregulated KoNAC genes (c) and one down-regulated KoNAC (d) under cold treatment were confirmed by qRT-PCR. Three independent experiments were performed. The actin gene in K. obovata acted as the internal control. Asterisks indicate significant differences compared with CK by Student’s t-test. *, p < 0.05.
Cimb 44 00381 g007
Table
Table 1. Basic information of K. obovata NAC genes.
Table 1. Basic information of K. obovata NAC genes.
NameGene IDClassChromosome PositionIntron NumberAverage Intron Length (bp)Protein Length (aa)pI
KoNAC1GWHPACBH000260.1VIIChr01: 1952420-195454922103758.33
KoNAC2GWHPACBH000261.1XChr01: 1966355-196829322133816.38
KoNAC3GWHPACBH001011.1IVChr01: 9324213-932593621792689.76
KoNAC4GWHPACBH001737.1IXChr02: 545874-54792825403144.72
KoNAC5GWHPACBH002133.1VIChr02: 5018997-502239226274144.56
KoNAC6GWHPACBH002150.1VIIChr02: 5630019-563284525572916.26
KoNAC7GWHPACBH002927.1XChr02: 12147292-1214979126953089.72
KoNAC8GWHPACBH002942.1VIIChr02: 12271496-1227388622653976.91
KoNAC9GWHPACBH003351.1VIIChr03: 1575855-1576587115217010.01
KoNAC10GWHPACBH003542.1XChr03: 2845984-28479713823588.29
KoNAC11GWHPACBH003714.1IVChr03: 4257070-4258248218821510.07
KoNAC12GWHPACBH004035.1VIChr03: 8571932-857584172675734.35
KoNAC13GWHPACBH004037.1VIChr03: 8580877-858558439224325.87
KoNAC14GWHPACBH004257.1IIIChr03: 10268576-1027332564843408.05
KoNAC15GWHPACBH005193.1VIIChr04: 4833704-483571922373016.63
KoNAC16GWHPACBH005487.1IIChr04: 8899908-890172125182597.96
KoNAC17GWHPACBH005488.1IIChr04: 8902766-890402215902208.48
KoNAC18GWHPACBH005795.1VIIChr04: 10981103-1098328422293179.64
KoNAC19GWHPACBH005903.1IIIChr04: 11744358-1174711441483045.55
KoNAC20GWHPACBH005980.1VIIChr05: 171989-17357421822858.57
KoNAC21GWHPACBH006496.1IIChr05: 3585985-358807021795765.03
KoNAC22GWHPACBH006945.1VIIChr05: 8286321-828955026264257.89
KoNAC23GWHPACBH007161.1VIIChr05: 10223818-1022599621713728.40
KoNAC24GWHPACBH007671.1IXChr06: 2167096-216867723003276.13
KoNAC25GWHPACBH007697.1VIChr06: 2347486-235131943855914.47
KoNAC26GWHPACBH007806.1VIIChr06: 3191353-319341722953036.78
KoNAC27GWHPACBH007934.1IIChr06: 4507683-451220755513655.29
KoNAC28GWHPACBH008265.1XChr06: 8671993-867429723744416.35
KoNAC29GWHPACBH009271.1VIChr07: 7943036-794675252725924.37
KoNAC30GWHPACBH009626.1XChr08: 595220-599512216243007.03
KoNAC31GWHPACBH009686.1VIIChr08: 1067480-106897721143439.60
KoNAC32GWHPACBH009687.1VIIChr08: 1074862-107694321713478.70
KoNAC33GWHPACBH009711.1XChr08: 1347371-135031626363577.67
KoNAC34GWHPACBH009970.1XChr08: 3208565-3211847210962867.10
KoNAC35GWHPACBH010230.1VChr08: 6804416-680740142743455.82
KoNAC36GWHPACBH010231.1VChr08: 6808948-681009931232614.47
KoNAC37GWHPACBH010248.1IXChr08: 6936960-693864721913667.37
KoNAC38GWHPACBH010352.1VIIChr08: 7704131-770596824353056.78
KoNAC39GWHPACBH010982.1VIIChr09: 1762983-176373111601769.98
KoNAC40GWHPACBH011181.1IIChr09: 3089069-309002511122366.50
KoNAC41GWHPACBH011224.1XChr09: 3378802-338043421153567.92
KoNAC42GWHPACBH011287.1VIIIChr09: 3895422-3899009212612434.62
KoNAC43GWHPACBH011356.1XChr09: 4526977-452903725533187.21
KoNAC44GWHPACBH012132.1IIIChr10: 2632918-263787155843737.08
KoNAC45GWHPACBH012473.1VIChr10: 5929512-593229832555744.20
KoNAC46GWHPACBH012540.1XChr10: 7425092-742831835793339.19
KoNAC47GWHPACBH012930.1VIIIChr11: 4340404-434363843423265.37
KoNAC48GWHPACBH013040.1IXChr11: 5432389-543408222853546.19
KoNAC49GWHPACBH013264.1IXChr11: 7036403-703805631653316.75
KoNAC50GWHPACBH013622.1VIIChr12: 1361271-136311722152949.80
KoNAC51GWHPACBH013650.1VChr12: 1584267-158751642195344.43
KoNAC52GWHPACBH013651.1VChr12: 1588784-159108641653366.58
KoNAC53GWHPACBH013779.1XChr12: 3635952-363894525913548.79
KoNAC54GWHPACBH013961.1XChr12: 5551002-555348235712569.26
KoNAC55GWHPACBH014293.1IIIChr12: 7828847-783088824942968.61
KoNAC56GWHPACBH014606.1VIIChr13: 304056-30620222302559.29
KoNAC57GWHPACBH014849.1VIIChr13: 2243283-224483921092657.90
KoNAC58GWHPACBH015137.1IIIChr13: 8039942-804424254664614.62
KoNAC59GWHPACBH015463.1VIIChr14: 3448051-345006222222559.64
KoNAC60GWHPACBH015782.1VChr15: 43374-4736753246384.57
KoNAC61GWHPACBH015958.1IXChr15: 1512975-151370311251629.55
KoNAC62GWHPACBH016068.1VIIIChr15: 2698343-270047142613069.33
KoNAC63GWHPACBH016810.1IIIChr16: 4714504-471646425163049.25
KoNAC64GWHPACBH017115.1VIIChr17: 1665396-166661621192458.31
KoNAC65GWHPACBH017251.1IVChr17: 2901510-290256711033147.77
KoNAC66GWHPACBH017590.1VIIChr17: 5655463-565707321473588.88
KoNAC67GWHPACBH017594.1VIIChr17: 5688319-568992821463588.88
KoNAC68GWHPACBH017666.1IIIChr17: 6866174-686866527273257.67
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, M.-M.; Liu, X.; Huang, X.-J.; Yang, J.-J.; Qin, P.-T.; Zhou, H.; Jiang, M.-G.; Liao, H.-Z. Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Kandelia obovata, a Typical Mangrove Plant. Curr. Issues Mol. Biol. 2022, 44, 5622-5637. https://doi.org/10.3390/cimb44110381

AMA Style

Sun M-M, Liu X, Huang X-J, Yang J-J, Qin P-T, Zhou H, Jiang M-G, Liao H-Z. Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Kandelia obovata, a Typical Mangrove Plant. Current Issues in Molecular Biology. 2022; 44(11):5622-5637. https://doi.org/10.3390/cimb44110381

Chicago/Turabian Style

Sun, Man-Man, Xiu Liu, Xiao-Juan Huang, Jing-Jun Yang, Pei-Ting Qin, Hao Zhou, Ming-Guo Jiang, and Hong-Ze Liao. 2022. "Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Kandelia obovata, a Typical Mangrove Plant" Current Issues in Molecular Biology 44, no. 11: 5622-5637. https://doi.org/10.3390/cimb44110381

APA Style

Sun, M. -M., Liu, X., Huang, X. -J., Yang, J. -J., Qin, P. -T., Zhou, H., Jiang, M. -G., & Liao, H. -Z. (2022). Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Kandelia obovata, a Typical Mangrove Plant. Current Issues in Molecular Biology, 44(11), 5622-5637. https://doi.org/10.3390/cimb44110381

Article Metrics

Back to TopTop