Xu Hu1,2†
Jiangshan Hao1,3†
Ling Pan2
Tao Xu1,2
Longzhou Ren1,2
Yu Chen4
Minqiang Tang2*
Li Liao1*
Zhiyong Wang2*- 1College of Tropical Crops, Hainan University, Haikou, China
- 2Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants, Ministry of Education, College of Forestry, Hainan University, Haikou, China
- 3School of Agriculture, Jinhua Polytechnic, Jinhua, China
- 4College of Agro-Grassland Science, Nanjing Agricultural University, Nanjing, China
Seashore paspalum (Paspalum vaginatum) is a halophytic, warm-season grass which is closely related to various grain crops. Gene duplication plays an important role in plant evolution, conferring significant plant adaptation at the genomic level. Here, we identified 2,542 tandem duplicated genes (TDGs) in the P. vaginatum genome and estimated the divergence time of pairs of TDGs based on synonymous substitution rates (Ks). Expression of P. vaginatum TDGs resulted in enrichment in many GO terms and KEGG pathways when compared to four other closely-related species. The GO terms included: “ion transmembrane transporter activity,” “anion transmembrane transporter activity” and “cation transmembrane transport,” and KEGG pathways included “ABC transport.” RNA-seq analysis of TDGs showed tissue-specific expression under salt stress, and we speculated that P. vaginatum leaves became adapted to salt stress in the earlier whole-genome duplication (WGD; ~83.3 million years ago; Ma), whereas the entire P. vaginatum plant acquired a large number of TDGs related to salt stress in the second WGD (~23.3 Ma). These results can be used as a reference resource to accelerate salt-resistance research in other grasses and crops.
Introduction
High salinity is a major abiotic environmental stress that is reported to be responsible for reductions in plant growth and crop production worldwide (Roy et al., 2014; Kumari et al., 2015). The production of salt-tolerant crops is potentially a cost-effective approach to provide improved growth in saline soils (Rasheed et al., 2022). It is therefore of critical interest to unravel the salt-resistance mechanisms of halophytes and transfer these, if possible, to glycophytes (Kumari et al., 2015). Seashore paspalum (Paspalum vaginatum Sw.) is a halophytic, warm-season, perennial grass that has been utilized as turf for almost 100 years, especially in coastal and salt-affected regions across the world (Wu et al., 2018; Qi et al., 2019). P. vaginatum is one of the most saline-tolerant turfgrass species (Liu et al., 2011; Uddin et al., 2012; Spiekerman and Devos, 2020) and, as it is closely related to some of the world’s most important grain crops, including maize, sorghum and millet (Qi et al., 2019), may provide a gateway for cereal crop improvements in salt resistance.
Gene duplication is a fundamental process in genome evolution (Holland, 1999; Paterson et al., 2010), and can result from whole-genome duplication, tandem duplication, duplication mediated by transposable elements, segmental duplication, and/or retroduplication (Panchy et al., 2016; Guo et al., 2019; Qiao et al., 2019). Tandem duplication refers to the generation of tandem arrays consisting of identical sequences in close genomic proximity. This occurs due to unequal chromosomal recombination, a widespread phenomenon in plant genomes which plays significant evolutionary roles including in adaptation to changing environments (Cannon et al., 2004; Yu et al., 2015). Tandem duplication events have been implicated in various plant traits such as stress resistance and membrane function in arabidopsis and rice (Rizzon et al., 2006), disease resistance in Solanaceae and Brassicaceae species (Leister, 2004), signal transduction in legumes (Bellieny-Rabelo et al., 2013) and glucosinolate biosynthesis diversification in the mustard family (Hofberger et al., 2013).
However, little is known about tandem duplicated genes (TDGs) and their possible contributions to the genome evolution of salt-stress resistance in P. vaginatum. Here, we report a comprehensive study of genome-wide TDGs present in the genome of P. vaginatum and surmise their evolutionary contributions. The functions of the TDGs were proposed using GO and KEGG enrichment analyses. The use of RNA-seq data made it feasible to identify TDGs that respond to salt stress in various tissues. Our findings can provide further insights into P. vaginatum evolution, particularly in relation to salt resistance. In addition, this genetic resource might also be useful for salt resistance research into other grasses and crops.
Materials and methods
Data sources
Protein sequence and General Feature Format (GFF) files of Oryza sativa, Setaria italica, Sorghum bicolor, and Zea mays were downloaded from EnsemblPlants.1 We obtained a chromosome-level reference genome of diploid cultivated P. vaginatum: The reference genome size was 517.98 million bases (Mb), including 28,712 codable proteins. For genes with multiple transcripts, the longest transcript was selected for subsequent analysis.
TDG analysis
TDGs were mainly identified based on protein sequences and the GFF file. Firstly, BLASTP (Altschul et al., 1997; settings: E < 1e–10; first 10 matches) was performed with protein sequences to search for all potential homologous gene pairs within each genome. Secondly, the blast results were analyzed in MCScanX (Wang et al., 2012) using a modified MCScan algorithm. Thirdly, duplicated gene pairs were identified using the downstream analysis tool (duplicate_gene_classifier) which is incorporated into the MCScanX package. Finally, duplicated gene pairs with code 3 (representing TDGs) were extracted. The non-synonymous (Ka) and synonymous substitution (Ks) frequencies of duplicated genes were calculated using ParaAT (Zhang et al., 2012), which was further used to compute the approximate dates of duplication and divergence events using the formula T = Ks/2λ, assuming a clocklike rate (λ) of 1.5 synonymous substitutions per 10−8 years (Yang et al., 2020; Wang et al., 2021b) for P. vaginatum. Moreover, the Ka/Ks ratio was also employed to show the selection pressure for the duplicated genes.
GO term and KEGG pathway enrichment analyses
Based on the annotation information2 (Huerta-Cepas et al., 2017), TDGs from P. vaginatum were analyzed for GO term and KEGG pathway functional enrichment using the R package clusterProfiler (Wu et al., 2021). “Rich factor” is the ratio of the number of differentially-expressed genes annotated in a term (or pathway) to the number of all genes annotated in this term (or pathway).
Identification and sequence analysis of the ABC gene family
The Hidden Markov Model (HMM) profile of the ABC domain (PF00664, PF00005, PF01061, PF02470 and PF01458) from the Pfam database3 (Mistry et al., 2021) was utilized to identify the ABC members in P. vaginatum by using the HMMER software (Mistry et al., 2013)(settings: E < 1e–10), and redundant sequences were removed manually. Additionally, all obtained ABC protein sequences were analyzed online using the CDD website4 (Marchler-Bauer et al., 2017) to verify conserved ABC domains. The members of the P. vaginatum ABC family were named according to chromosomal position. The number of amino acids, theoretical molecular weight (MW), and theoretical isoelectric point (pI) of the ABC family were obtained from the ExPASy web resource5 (Gasteiger et al., 2003). In order to explore the phylogenetic relationships of P. vaginatum ABC-family genes, data from 105 ABC proteins of O. sativa were downloaded for multiple sequence alignment using MAFFT (Rozewicki et al., 2019), and multiple sequences were trimmed using TBtools (Chen et al., 2020). Phylogenetic trees (ML) were constructed using IQ-TREE (Minh et al., 2020) (settings: -m MFP -bb 1000 -alrt 1000) and generated using iTol online tools6 (Letunic and Bork, 2021).
RNA-seq and bioinformatics analysis
The P. vaginatum salt-resistant ecotypes USA17-18 were analyzed by RNA-seq. For the RNA-seq experiments, similar stolons were cultivated in plastic pots under typical conditions. Two-month-old plants were treated with 400 mM NaCl or water as a control for 8, 12, 24, or 48 h or 5 days before harvesting tissues. Total RNAs were isolated from leaves with three biological replicates at each stress stage. RNA-seq and de novo assembly Paired-end sequencing of cDNA libraries were performed using the HiSeq 2000 platform (Supplementary Table 10). Clean reads were aligned to the reference genome using HISAT2 (Kim et al., 2015). FPKM (fragments per kilobase per million mapped reads) was used to estimate the expression levels of individual genes (Supplementary Table 11). Differentially expressed genes (DEGs) were determined using the R package DESeq2 (Love et al., 2014) with a false discovery rate (FDR) of ≤ 0.05 and |log2 FC| ≥ 1 used as the threshold to determine statistically significant differences in gene expression. To determine whether these genes were tissue-specific to leaves or roots, gene expression patterns were compared in both transcriptomic data of control and treatment plants. Tissue-specifically expressed genes and co-expressed genes were defined based on the following rules: If the gene appeared up-regulated at least four times in one tissue but less than (or equal to) one time in another tissue, we considered it to be a tissue-specifically expressed gene; if the gene appeared up-regulated at least seven times in both tissues, it was considered to be a co-expressed gene.
Results
Identification and genomic distribution of TDGs
The P. vaginatum genome sequence consists of 517.98 Mb of DNA sequence and 28,712 protein coding genes. A total of 2,542 TDGs (8.85% of the gene set) were identified in the P. vaginatum genome, with a higher frequency than O. sativa (7.78%) and Z. mays (4.74%) but lower frequency than S. italica (11.55%) and S. bicolor (10.82%) (Supplementary Table 1). In addition, the distribution of TDGs was not uniform across chromosomes (Figure 1A): The most, 428, TDGs (16.84% of total TDGs), were located on chromosome 1 (chr01); the least, 134 (5.27% of the total TDGs) on chromosome 8 (chr08) (Supplementary Table 2). The synonymous substitution rates (Ks) in the TDG pairs were calculated for 1,581 gene pairs, and the distribution of Ks showed a single peak value ranging from 0.6 to 0.7 (Figure 1B; Supplementary Table 3). To infer the speciation occurrence time of P. vaginatum, the T = Ks/2λ method was used. Results indicated an estimated time of approximately 20.0–23.3 Ma, suggesting that P. vaginatum and S. italica might have shared common whole-genome duplicates. Ka/Ks values were then used to determine possible selection pressure between individual genes in a pair. The Ka/Ks values for P. vaginatum TDGs showed that there were few TDG pairs with Ka/Ks values much greater than 1, with most TDG pairs having Ka/Ks values much less than 1 (Supplementary Table 3), suggesting that most TDGs were under negative selection and that only a small number were under positive selection.
Figure 1. Overview of the Paspalum vaginatum genome and comparison of Ks. (A) Tracks from outer (a) to inner (e) rings indicate the following: (a) Chromosome size with units in Mb; (b) density of genes; (c) density of TDGs; (d) Ks of TDG pairs; (e) Ka/Ks values of TDG pairs and normalization by log2 (Ka/Ks). (B) Distribution of Ks calculated by TDG pairs between study species.
The number of TDGs in the same tandem cluster ranged from two to nine with the longest tandem clusters found with “UDP-glycosyltransferase” and “GRAS domain” (Table 1). We focused on TDGs with more than six genes in the same cluster and their functions (Supplementary Table 4). The majority of these long TDG clusters were found to be associated with salt resistance after reviewed evidence and included: “UDP-glycosyltransferase,” “GRAS domain,” “Auxin responsive protein,” “Cytochrome P450,” “FCS-type zinc-finger,” “Dirigent protein,” “BTB/POZ domain,” “Expansin, cellulose-binding-like domain,” “Glutathione S-transferase,” “RING-type zinc-finger” and “FAD-linked oxidoreductase” (Table 1).
Table 1. TDGs with more than six genes in the same cluster and their functions.
The TDGs in Paspalum vaginatum contribute to adaptability, according to enrichment analyses
To gain insights into the biological processes necessary for the adaptation to the environment, the 2,542 TDGs in P. vaginatum were analyzed for GO enrichment. The set of TDGs was involved in 177 significant biological processes. The maximum Rich Factor was 0.75 (with xyloglucan biosynthesis), followed by 0.70 (with cellular response to high light intensity; Figure 2B; Supplementary Table 5). Then, a comparative analysis of four grass species (O. sativa, S. italica, S. bicolor, and Z. mays) with P. vaginatum was performed. Based on the GO terms, the P. vaginatum TDGs were found to be enriched in 62 unique GO terms compared with other analyzed species (Figure 2A). These GO terms included “cellular response to light intensity,” “cellular response to UV” and “cellular response to heat” etc. (Figure 2C), which are related to adaptation to tropical climates. Other enriched GO terms included “ion transmembrane transporter activity,” “anion transmembrane transporter activity” and “cation transmembrane transport” etc. (Figure 2C), which might be associated with adaptation to a saline environment.
Figure 2. GO enrichment analysis of TDGs. (A) The venn diagram represents the shared and unique enriched GO terms among five species. (B) The size of enriched GO terms in each species. (C) The bubble diagram shows 30 significantly enriched GO terms. The terms in red are unique to P. vaginatum compared to the other four species.
KEGG pathway analysis provides classifications that are valuable for studying the complex biological functions of genes. Therefore, KEGG enrichment analysis was also performed on TDGs of the studied species. The results showed that TDGs in P. vaginatum were significantly enriched in 34 pathways (Figure 3B; Supplementary Table 6). The pathways enriched with the top three numbers of TDGs were “Phenylpropanoid biosynthesis,” “MAPK signaling pathway” and “Metabolism of xenobiotics by cytochrome P450.” In comparison to the other four species, P. vaginatum had six unique pathways (Figures 3A,C), five of which were related to KEGG pathway Metabolism (“Prodigiosin biosynthesis,” “Tropane, piperidine and pyridine alkaloid biosynthesis,” “Caprolactam degradation,” “Nitrogen metabolism,” “Biosynthesis of unsaturated fatty acids”), and one of which was related to Environmental Information Processing (“ABC transporters”).
Figure 3. KEGG enrichment analysis of TDGs. (A) The venn diagram represents the shared and unique enriched KEGG pathways among five species. (B) The size of enriched KEGG pathways in each species. (C) The bubble diagram shows 30 significantly enriched KEGG pathways. The pathways in red are unique to P. vaginatum compared to the other four species.
TDGs were tissue-specific under salt stress
Time-ordered, comparative, transcriptome analyses were performed, which focused on the expression of TDGs under salt stress. A total of 302 TDGs were identified as having up-regulated expression in response to salt stress (Supplementary Table 7), of which 48 were specifically expressed in leaves (Figure 4A); 116 were specifically expressed in roots (Figure 4B), and 138 were co-expressed in both roots and leaves (Figure 4C). These genes were annotated according to the NR, Swissprot and Pfam databases, and then literature searches were conducted to confirm whether their functions were related to salt resistance. It was found that most TDGs appeared to respond to salt stress. Some of these genes had been annotated as “Uncharacterized protein” (Supplementary Table 7), and it can be speculated that these genes also play an important role in the adaptation of P. vaginatum to salt stress. Even in the same gene family, some family members were expressed only in leaves but others in roots, such as “Cytochrome P450,” “UDP glycosyltransferase” and “Multicopper oxidase family” etc.
Figure 4. Tissue-specific expression of TDGs under salt stress. (A) Specifically expressed TDGs in leaves. Genes related to salt resistance with reviewed evidence were marked in red. (B) Specifically expressed TDGs in roots. (C) Expressed TDGs in both roots and leaves. (D) Distribution of Ks with different specific expression TDGs.
The frequencies of synonymous substitution (Ks) were calculated to estimate the age of duplication events. There were similar peaks (~0.7) among the tissue-specifically expressed TDGs and the TDGs in whole genome. However, there was another peak (~2.5) with only tissue-specifically-expressed TDGs in leaves (Figure 4D), which included the gene functions for “Cation transporter,” “Sulfate transporter,” and “UDP-glycosyltransferase” etc. (Table 2).
Table 2. Specifically expressed TDGs in leaves with Ks of about 2.5.
Analysis of ABC gene family in Paspalum vaginatum
A total of 131 P. vaginatum ABC protein sequences were identified. To further understand the relationship between the ABC family genes, full-length protein sequences of P. vaginatum and O. sativa ABC proteins were aligned to construct a phylogenetic tree. The result showed that the ABC genes of these two species can be divided into eight subfamilies (Supplementary Figure 1). The ABC family genes in P. vaginatum were denoted as PvABCA1-PvABCA6, PvABCB1-PvABCB33, PvABCC1-PvABCC17, PvABCD1-PvABCD5, PvABCE1-PvABCE4, PvABCF1-PvABCF9, PvABCG1-PvABCG45, and PvABCI1-PvABCI12, according to their chromosomal location (Figure 5). Among these subfamilies, PvABCG had the largest number of members with 45 genes; PvABCB was the second largest subgroup, containing 33 genes; and the smallest subgroup was PvABCE, containing only 4 genes. Molecular masses ranged from 15.81 kDa to 596.82 kDa and pIs ranged from 4.81 to 9.87 (Supplementary Table 8).
Figure 5. Phylogenetic tree of ABC transporter proteins from P. vaginatum. Eight subfamilies are highlighted in different colors. The genes marked with black dots represent TDGs. Heatmap shows transcriptome expression under salt stress.
To better understand the biological functions of the ABC family genes, the combined P. vaginatum transcriptome under salt stress was analyzed for expression patterns. As shown in Figure 5, the PvABCs exhibited different expression levels under salt stress. PvABCB14 and PvABCB15, identified as TDGs, were highly upregulated in roots, whereas PvABCC8 and PvABCC9, also identified as TDGs, were highly upregulated in leaves under salt stress. Many PvABCs were not differentially expressed.
Discussion
Long-tandem clusters in Paspalum vaginatum are strongly associated with abiotic stress
With long-tandem clusters (more than six genes in the same cluster), chromosomes 3 and 6 had the largest numbers of long clusters (four and three, respectively). A total of 14 long-tandem clusters (containing 97 TDGs) were found to be associated with abiotic stress, with 12 specifically associated with salt stress (Table 1). Two clusters in particular (with nine and eight genes in each cluster) with function “UDP-glucosyl transferase” were of particular interest. Further research showed that “UDP-glucosyl transferase” can be associated with elevated antioxidant enzyme activity and reduced production of reactive oxygen species, which could control the oxidative burst under stress situations (Wang et al., 2021a). In the present study, many “UDP-glucosyl transferase” genes were identified resulting from tandem-duplication events, which have enriched the antioxidant capacity of P. vaginatum.
The unique GO terms and KEGG pathways of TDGs in Paspalum vaginatum seem to indicate contributions to salt tolerance
GO and KEGG enrichment analyses were used to determine the possible roles of TDGs in P. vaginatum, and were also compared in four related species (O. sativa, S. italica, S. bicolor and Z. mays). It was discovered that P. vaginatum TDGs were associated with some unique GO terms (such as “ion transmembrane transporter activity,” “anion transmembrane transporter activity” and “cation transmembrane transport”). These GO terms indicate functions promoting osmoregulation, for example by accumulating compatible solutes to avoid osmotic stress caused by salinity (Kumari et al., 2015; Rahman et al., 2021). P. vaginatum TDGs were also found to be associated with some unique KEGG pathways, and one being “ABC transport.” ABC transporters belong to a large protein family that utilize the energy released by ATP hydrolysis to transport a wide range of substrates across biological membranes. ABC transporters are involved in diverse cellular processes such as biotic and abiotic stress responses through plant hormone transport, ion transport, lipid transport and redox homeostasis (Do et al., 2018; Dahuja et al., 2021). In the present study, 131 ABC transporter genes were identified in P. vaginatum, with 13 (40%) of the 33 genes identified being TDGs. PvABCB was the second largest subgroup in the ABC transporter family and contained the most TDGs. PvABCB14 and PvABCB15, identified from TDGs, were highly upregulated in roots, which might have an important role in enhancing the salt resistance of P. vaginatum.
The evolution of roots and leaves in response to salt stress
With transcriptome expression available at multiple time points under salt stress, P. vaginatum TDGs were divided into those with leaf-specific expression (48 TDGs), root-specific expression (116 TDGs), or co-expression in both leaves and roots (138 TDGs). Compared to leaves, roots had more TDGs in response to salt stress, possibly indicating that roots are more important than leaves in adapting to saline environments. Interestingly, when Ks was calculated for tissue-specifically expressed TDGs and compared with that for TDGs in whole genome, a similar peak (~0.7) was found, whereas another peak (~2.5) was found only in tissue-specifically expressed TDGs in leaves. Therefore, it was speculated that leaves became adapted to salt stress in an earlier whole genome duplication event (WGD; ~83.3 million years ago; Ma), whereas the entire P. vaginatum plant acquired a large number of TDGs related to salt stress in a second WGD (~23.3 Ma). Annotations of these expressed TDGs were found from multiple databases and it was found that most were related to salt stress. We therefore speculate that expressed TDGs annotated as “Uncharacterized protein” also contribute to salt resistance in P. vaginatum, and these genes could serve as potential and novel genes related to salt resistance.
Subfunctionalization of tandem duplicated pairs
Duplicate genes can diverge in function, representing a potential source of response mechanisms to survive stressful environments (Zou et al., 2009; Arsovski et al., 2015) and the TDGs in P. vaginatum have shown functional divergence. For example, emFS26.201 and emFS26.202 are tandem duplicated pairs (annotated as “UDP-glycosyltransferase”), but emFS26.201 has leaf-specific expression while emFS26.202 has root-specific expression under salt stress. As another example, from TDGs in the same tandem cluster, emOS168.286, emOS168.288, emOS168.289, and emOS168.290 (annotated as “Peroxidase”), only emOS168.290 responds to salt stress (Supplementary Table 9). A total of 302 TDGs were identified in response to salt stress, and these TDGs were part of tandem clusters comprising 623 TDGs (Supplementary Table 9). In other words, nearly half of these duplicate lineages retain stress responsiveness.
The role of TDGs in evolution and adaptation
Evolutionary innovation is often built on variations from redundant genetic materials generated by gene duplication, and tandem duplication represents a potential source for such innovation (Arsovski et al., 2015; Guo et al., 2019). Previous studies have shown that: TDGs have driven diversification of the protein modifier SUMO in angiosperms (Hammoudi et al., 2016); TDGs have driven the divergent evolution of caffeine and crocin biosynthetic pathways in Gardenia jasminoides (Xu et al., 2020); TDGs are involved in defense and pollinator attraction in Tectona grandis (Zhao et al., 2019); and lastly that TDGs are significantly enriched in resistance-related pathways and more abundant in retrotransposon-related genes in Cajanus cajan (Liu et al., 2021). In this study, we identified TDGs in the P. vaginatum and found that TDGs of P. vaginatum were enriched in many unique GO terms and KEGG pathways which were associated with resistance, especially salt resistance. The TDGs associated with response to salt stress were identified by transcriptome analyses at multiple time points. This study provides insights into the roles of tandem duplications in the evolution and adaptation of P. vaginatum and lays the foundation for the genomics-based breeding of other grasses.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material. RNA-Seq data obtained have been uploaded to the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/) and can be accessed with accession number PRJNA874860.
Author contributions
XH and J-SH contributed equally to this work performed data analyses, and wrote the manuscript. Z-YW and LL conceived and managed the project. MQ-T and YC designed the experiments. TX and L-ZR identified genes related to salt resistance through literature search. LP interpreted the results and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (No.32060409), Key scientific research projects of colleges and universities in Hainan Province (No. Hnky2019ZD-3), the Construction of World First Class Discipline of Hainan University (No.RZZX201905), and National Project on Sci-Tec Foundation Resources Survey (2017FY100600).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.971999/full#supplementary-material
SUPPLEMENTARY FIGURE 1 | Phylogenetic tree of ABC transporter proteins from P. vaginatum and O. sativa. Eight subfamilies are highlighted in different colors.
SUPPLEMENTARY TABLE 1 | Statistics of TDGs in P. vaginatum and other studied species.
SUPPLEMENTARY TABLE 2 | Number of TDGs in chromosome of P. vaginatum genome.
SUPPLEMENTARY TABLE 3 | Calculation of synonymous and nonsynonymous substitutions per site of TDG pairs.
SUPPLEMENTARY TABLE 4 | TDGs with more than six genes in the same cluster and their functions.
SUPPLEMENTARY TABLE 5 | GO enrichment results of TDGs in P. vaginatum.
SUPPLEMENTARY TABLE 6 | KEGG enrichment results of TDGs in P. vaginatum.
SUPPLEMENTARY TABLE 7 | TDGs which response to salt stress in P. vaginatum.
SUPPLEMENTARY TABLE 8 | Basic information about the ABC transporter gene family in P. vaginatum.
SUPPLEMENTARY TABLE 9 | Information of tandem clusters containing 302 tissue-specific TDGs.
SUPPLEMENTARY TABLE 10 | Quality control of RNA-seq.
SUPPLEMENTARY TABLE 11 | Expression (FPKM) of genes in P. vaginatum under salt stress.
Footnotes
1. ^http://plants.ensembl.org/index.html
2. ^http://eggnog-mapper.embl.de/
4. ^https://www.ncbi.nlm.nih.gov/cdd
References
Agarwal, P., and Khurana, P. (2020). TaZnF, a C3HC4 type RING zinc finger protein from Triticum aestivum is involved in dehydration and salinity stress. J. Plant Biochem. Biotechnol. 29, 395–406. doi: 10.1007/s13562-019-00546-8
Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Arsovski, A. A., Pradinuk, J., Guo, X. Q., Wang, S., and Adams, K. L. (2015). Evolution of Cis-regulatory elements and regulatory networks in duplicated genes of Arabidopsis. Plant Physiol. 169, 2982–2991. doi: 10.1104/pp.15.00717
Bellieny-Rabelo, D., Oliveira, A. E., and Venancio, T. M. (2013). Impact of whole-genome and tandem duplications in the expansion and functional diversification of the F-box family in legumes (Fabaceae). PLoS One 8:e55127. doi: 10.1371/journal.pone.0055127
Bouzroud, S., Gouiaa, S., Hu, N., Bernadac, A., Mila, I., Bendaou, N., et al. (2018). Auxin response factors (ARFs) are potential mediators of auxin action in tomato response to biotic and abiotic stress (Solanum lycopersicum). PLoS One 13:e0193517. doi: 10.1371/journal.pone.0193517
Cannon, S. B., Mitra, A., Baumgarten, A., Young, N. D., and May, G. (2004). The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 4:10. doi: 10.1186/1471-2229-4-10
Chen, C. J., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y. H., et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194–1202. doi: 10.1016/j.molp.2020.06.009
Chen, Y. H., Han, Y. Y., Kong, X. Z., Kang, H. H., Ren, Y. Q., and Wang, W. (2017). Ectopic expression of wheat expansin gene TaEXPA2 improved the salt tolerance of transgenic tobacco by regulating Na+/K+ and antioxidant competence. Physiol. Plant. 159, 161–177. doi: 10.1111/ppl.12492
Dahuja, A., Kumar, R. R., Sakhare, A., Watts, A., Singh, B., Goswami, S., et al. (2021). Role of ATP-binding cassette transporters in maintaining plant homeostasis under abiotic and biotic stresses. Physiol. Plant. 171, 785–801. doi: 10.1111/ppl.13302
Do, T. H. T., Martinoia, E., and Lee, Y. (2018). Functions of ABC transporters in plant growth and development. Curr. Opin. Plant Biol. 41, 32–38. doi: 10.1016/j.pbi.2017.08.003
Gasteiger, E., Gattiker, A., Hoogland, C., Ivanyi, I., Appel, R. D., and Bairoch, A. (2003). ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 31, 3784–3788. doi: 10.1093/nar/gkg563
Guo, H., Jiao, Y., Tan, X., Wang, X., Huang, X., Jin, H., et al. (2019). Gene duplication and genetic innovation in cereal genomes. Genome Res. 29, 261–269. doi: 10.1101/gr.237511.118
Hammoudi, V., Vlachakis, G., Schranz, M. E., and van den Burg, H. A. (2016). Whole-genome duplications followed by tandem duplications drive diversification of the protein modifier SUMO in angiosperms. New Phytol. 211, 172–185. doi: 10.1111/nph.13911
Hofberger, J. A., Lyons, E., Edger, P. P., Chris Pires, J., and Eric Schranz, M. (2013). Whole genome and tandem duplicate retention facilitated glucosinolate pathway diversification in the mustard family. Genome Biol. Evol. 5, 2155–2173. doi: 10.1093/gbe/evt162
Holland, P. W. (1999). Gene duplication: past, present and future. Semin. Cell Dev. Biol. 10, 541–547. doi: 10.1006/scdb.1999.0335
Huerta-Cepas, J., Forslund, K., Coelho, L. P., Szklarczyk, D., Jensen, L. J., von Mering, C., et al. (2017). Fast genome-wide functional annotation through orthology assignment by egg NOG-mapper. Mol. Biol. Evol. 34, 2115–2122. doi: 10.1093/molbev/msx148
Jung, Y. J., Lee, I. H., Nou, I. S., Lee, K. D., Rashotte, A. M., and Kang, K. K. (2013). BrRZFP1 a Brassica rapa C3HC4-type RING zinc finger protein involved in cold, salt and dehydration stress. Plant Biol. 15, 274–283. doi: 10.1111/j.1438-8677.2012.00631.x
Kang, C., He, S. Z., Zhai, H., Li, R. J., Zhao, N., and Liu, Q. C. (2018). A sweetpotato auxin response factor gene (IbARF5) is involved in carotenoid biosynthesis and salt and drought tolerance in transgenic Arabidopsis. Front. Plant Sci. 9:1307. doi: 10.3389/fpls.2018.01307
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Krishnamurthy, P., Vishal, B., Ho, W. J., Lok, F. C. J., Lee, F. S. M., and Kumar, P. P. (2020). Regulation of a cytochrome P450 gene CYP94B1 by WRKY33 transcription factor controls Apoplastic barrier formation in roots to confer salt tolerance. Plant Physiol. 184, 2199–2215. doi: 10.1104/pp.20.01054
Kumari, A., Das, P., Parida, A. K., and Agarwal, P. K. (2015). Proteomics, metabolomics, and ionomics perspectives of salinity tolerance in halophytes. Front. Plant Sci. 6:537. doi: 10.3389/fpls.2015.00537
Leister, D. (2004). Tandem and segmental gene duplication and recombination in the evolution of plant disease resistance gene. Trends Genet. 20, 116–122. doi: 10.1016/j.tig.2004.01.007
Letunic, I., and Bork, P. (2021). Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293–W296. doi: 10.1093/nar/gkab301
Li, P., Li, Y. J., Zhang, F. J., Zhang, G. Z., Jiang, X. Y., Yu, H. M., et al. (2017). The Arabidopsis UDP-glycosyltransferases UGT79B2 and UGT79B3, contribute to cold, salt and drought stress tolerance via modulating anthocyanin accumulation. Plant J. 89, 85–103. doi: 10.1111/tpj.13324
Liu, L. M., Du, H. M., Wang, K., Huang, B. R., and Wang, Z. L. (2011). Differential photosynthetic responses to salinity stress between two perennial grass species contrasting in salinity tolerance. HortScience 46, 311–316. doi: 10.21273/Hortsci.46.2.311
Liu, C., Wu, Y. H., Liu, Y. X., Yang, L. Y., Dong, R. S., Jiang, L. Y., et al. (2021). Genome-wide analysis of tandem duplicated genes and their contribution to stress resistance in pigeonpea (Cajanus cajan). Genomics 113, 728–735. doi: 10.1016/j.ygeno.2020.10.003
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Magwanga, R. O., Lu, P., Kirungu, J. N., Dong, Q., Cai, X. Y., Zhou, Z. L., et al. (2019). Knockdown of cytochrome P450 genes Gh_D07G1197 and Gh_A13G2057 on chromosomes D07 and A13 reveals their putative role in enhancing drought and salt stress tolerance in Gossypium hirsutum. Genes (Basel) 10:226. doi: 10.3390/genes10030226
Marchler-Bauer, A., Bo, Y., Han, L., He, J., Lanczycki, C. J., Lu, S., et al. (2017). CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 45, D200–D203. doi: 10.1093/nar/gkw1129
Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., von Haeseler, A., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015
Mistry, J., Chuguransky, S., Williams, L., Qureshi, M., Salazar, G. A., Sonnhammer, E. L. L., et al. (2021). Pfam: the protein families database in 2021. Nucleic Acids Res. 49, D412–D419. doi: 10.1093/nar/gkaa913
Mistry, J., Finn, R. D., Eddy, S. R., Bateman, A., and Punta, M. (2013). Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 41:e121. doi: 10.1093/nar/gkt263
Panchy, N., Lehti-Shiu, M., and Shiu, S. H. (2016). Evolution of gene duplication in plants. Plant Physiol. 171, 2294–2316. doi: 10.1104/pp.16.00523
Paterson, A. H., Freeling, M., Tang, H., and Wang, X. (2010). Insights from the comparison of plant genome sequences. Annu. Rev. Plant Biol. 61, 349–372. doi: 10.1146/annurev-arplant-042809-112235
Qi, P., Eudy, D., Schnable, J. C., Schmutz, J., Raymer, P. L., and Devos, K. M. (2019). High density genetic maps of seashore paspalum using genotyping-by-sequencing and their relationship to the sorghum bicolor genome. Sci. Rep. 9:12183. doi: 10.1038/s41598-019-48257-3
Qiao, X., Li, Q. H., Yin, H., Qi, K. J., Li, L. T., Wang, R. Z., et al. (2019). Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants. Genome Biol. 20:38. doi: 10.1186/s13059-019-1650-2
Qin, Y. X., Cui, S. F., Cui, P., Zhang, B., and Quan, X. Y. (2021). TaFLZ2D enhances salinity stress tolerance via superior ability for ionic stress tolerance and ROS detoxification. Plant Physiol. Biochem. 168, 516–525. doi: 10.1016/j.plaphy.2021.11.014
Rahman, M. M., Mostofa, M. G., Keya, S. S., Siddiqui, M. N., Ansary, M. M. U., Das, A. K., et al. (2021). Adaptive mechanisms of halophytes and their potential in improving salinity tolerance in plants. Int. J. Mol. Sci. 22:10733. doi: 10.3390/ijms221910733
Rasheed, A., Li, H., Nawaz, M., Mahmood, A., Hassan, M. U., Shah, A. N., et al. (2022). Molecular tools, potential frontiers for enhancing salinity tolerance in rice: A critical review and future prospective. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.966749
Rizzon, C., Ponger, L., and Gaut, B. S. (2006). Striking similarities in the genomic distribution of tandemly arrayed genes in Arabidopsis and rice. PLoS Comput. Biol. 2:e115. doi: 10.1371/journal.pcbi.0020115
Roy, S. J., Negrao, S., and Tester, M. (2014). Salt resistant crop plants. Curr. Opin. Biotechnol. 26, 115–124. doi: 10.1016/j.copbio.2013.12.004
Rozewicki, J., Li, S. L., Amada, K. M., Standley, D. M., and Katoh, K. (2019). MAFFT-DASH: integrated protein sequence and structural alignment. Nucleic Acids Res. 47, W5–W10. doi: 10.1093/nar/gkz342
Spiekerman, J. J., and Devos, K. M. (2020). The halophyte seashore paspalum uses adaxial leaf papillae for sodium sequestration. Plant Physiol. 184, 2107–2119. doi: 10.1104/pp.20.00796
Uddin, M. K., Juraimi, A. S., Ismail, M. R., Hossain, M. A., Othman, R., and Rahim, A. A. (2012). Physiological and growth responses of six Turfgrass species relative to salinity tolerance. Sci. World J 2012:905468. doi: 10.1100/2012/905468
Wan, X., Peng, L., Xiong, J., Li, X. Y., Wang, J. M., Li, X. F., et al. (2019). AtSIBP1, a novel BTB domain-containing protein, positively regulates salt signaling in Arabidopsis thaliana. Plants Basel 8:573. doi: 10.3390/plants8120573
Wang, Y. Y., Cao, Y. B., Liang, X. Y., Zhuang, J. H., Wang, X. F., Qin, F., et al. (2022). A dirigent family protein confers variation of casparian strip thickness and salt tolerance in maize. Nat. Commun. 13:2222. doi: 10.1038/s41467-022-29809-0
Wang, Y., Chen, F., Ma, Y., Zhang, T., Sun, P., Lan, M., et al. (2021b). An ancient whole-genome duplication event and its contribution to flavor compounds in the tea plant (Camellia sinensis). Hortic Res. 8:176. doi: 10.1038/s41438-021-00613-z
Wang, Z. B., Li, H., Wei, Z. P., Sun, H. R., He, Y., Gao, J., et al. (2021c). Overexpression of UDP-glycosyltransferase genes enhanced aluminum tolerance through disrupting cell wall polysaccharide components in soybean. Plant Soil 469, 135–147. doi: 10.1007/s11104-021-05157-8
Wang, T., Ma, Y. Q., Huang, X. X., Mu, T. J., Li, Y. J., Li, X. K., et al. (2021a). Overexpression of OsUGT3 enhances drought and salt tolerance through modulating ABA synthesis and scavenging ROS in rice. Environ. Exp. Bot. 192:104653. doi: 10.1016/j.envexpbot.2021.104653
Wang, Y., Tang, H., Debarry, J. D., Tan, X., Li, J., Wang, X., et al. (2012). MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40:e49. doi: 10.1093/nar/gkr1293
Wang, T. T., Yu, T. F., Fu, J. D., Su, H. G., Chen, J., Zhou, Y. B., et al. (2020). Genome-wide analysis of the GRAS gene family and functional identification of GmGRAS37 in drought and salt tolerance. Front. Plant Sci. 11:604690. doi: 10.3389/fpls.2020.604690
Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., et al. (2021). clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovations 2:100141. doi: 10.1016/j.xinn.2021.100141
Wu, J. Q., Liu, W. K., Jahan, M. S., Shu, S., Sun, J., and Guo, S. R. (2022). Characterization of polyamine oxidase genes in cucumber and roles of CsPAO3 in response to salt stress. Environ. Exp. Bot. 194:104696. doi: 10.1016/j.envexpbot.2021.104696
Wu, X. L., Shi, H. F., and Guo, Z. F. (2018). Overexpression of a NF-YC gene results in enhanced drought and salt tolerance in transgenic seashore paspalum. Front. Plant Sci. 9:1355. doi: 10.3389/fpls.2018.01355
Xu, Z., Pu, X., Gao, R., Demurtas, O. C., Fleck, S. J., Richter, M., et al. (2020). Tandem gene duplications drive divergent evolution of caffeine and crocin biosynthetic pathways in plants. BMC Biol. 18:63. doi: 10.1186/s12915-020-00795-3
Xu, J., Xing, X. J., Tian, Y. S., Peng, R. H., Xue, Y., Zhao, W., et al. (2015). Transgenic Arabidopsis plants expressing tomato glutathione S-Transferase showed enhanced resistance to salt and drought stress. PLoS One 10:e0136960. doi: 10.1371/journal.pone.0136960
Yang, Q., Liu, Y. J., and Zeng, Q. Y. (2019). Overexpression of three orthologous glutathione S-transferases from populus increased salt and drought resistance in Arabidopsis. Biochem. Syst. Ecol. 83, 57–61. doi: 10.1016/j.bse.2019.01.001
Yang, F. S., Nie, S., Liu, H., Shi, T. L., Tian, X. C., Zhou, S. S., et al. (2020). Chromosome-level genome assembly of a parent species of widely cultivated azaleas. Nat. Commun. 11:5269. doi: 10.1038/s41467-020-18771-4
Yu, J., Ke, T., Tehrim, S., Sun, F., Liao, B., and Hua, W. (2015). PTGBase: an integrated database to study tandem duplicated genes in plants. Database (Oxford) 2015:bav017. doi: 10.1093/database/bav017
Zhang, S., Li, X. W., Fan, S. D., Zhou, L. J., and Wang, Y. (2020). Overexpression of HcSCL13, a Halostachys caspica GRAS transcription factor, enhances plant growth and salt stress tolerance in transgenic Arabidopsis. Plant Physiol. Biochem. 151, 243–254. doi: 10.1016/j.plaphy.2020.03.020
Zhang, Z., Xiao, J. F., Wu, J. Y., Zhang, H. Y., Liu, G. M., Wang, X. M., et al. (2012). ParaAT: a parallel tool for constructing multiple protein-coding DNA alignments. Biochem. Biophys. Res. Commun. 419, 779–781. doi: 10.1016/j.bbrc.2012.02.101
Zhao, D., Hamilton, J. P., Bhat, W. W., Johnson, S. R., Godden, G. T., Kinser, T. J., et al. (2019). A chromosomal-scale genome assembly of Tectona grandis reveals the importance of tandem gene duplication and enables discovery of genes in natural product biosynthetic pathways. Gigascience 8:giz005. doi: 10.1093/gigascience/giz005
Zhao, Y., Li, X., Zhang, Z. X., Pan, W. J., Li, S. N., Xing, Y., et al. (2021). GmGPDH12, a mitochondrial FAD-GPDH from soybean, increases salt and osmotic stress resistance by modulating redox state and respiration. Crop J. 9, 79–94. doi: 10.1016/j.cj.2020.05.008
Keywords: abiotic stress, salt, seashore paspalum, tandem duplicated genes, transcriptome
Citation: Hu X, Hao J-S, Pan L, Xu T, Ren L-Z, Chen Y, Tang M-Q, Liao L and Wang Z-Y (2022) Genome-wide analysis of tandem duplicated genes and their expression under salt stress in seashore paspalum. Front. Plant Sci. 13:971999. doi: 10.3389/fpls.2022.971999
Edited by:
Jin-Lin Zhang, Lanzhou University, ChinaReviewed by:
Xiaojun Nie, Northwest A&F University, ChinaLiang Chen, University of Chinese Academy of Sciences, China
Copyright © 2022 Hu, Hao, Pan, Xu, Ren, Chen, Tang, Liao and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Minqiang Tang, tangminqiang@hainanu.edu.cn; Li Liao, liaoli@hainanu.edu.cn; Zhiyong Wang, wangzhiyong@hainanu.edu.cn
†These authors have contributed equally to this work