Next Article in Journal
Immunophenotyping of a Stromal Vascular Fraction from Microfragmented Lipoaspirate Used in Osteoarthritis Cartilage Treatment and Its Lipoaspirate Counterpart
Previous Article in Journal
Genome-Wide Identification and Expression Analyses of the Chitinases under Cold and Osmotic Stress in Ammopiptanthus nanus
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 6
10.3390/genes10060473
genes-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 Analysis of Class III Peroxidases in Allotetraploid Cotton (Gossypium hirsutum L.) and their Responses to PK Deficiency

by
Pengfei Duan
1,2,†,
Guo Wang
1,†,
Maoni Chao
1,
Zhiyong Zhang
1,2,* and
Baohong Zhang
1,3
1
Henan Collaborative Innovation Center of Modern Biological Breeding and Henan Key Laboratory for Molecular Ecology and Germplasm Innovation of Cotton and Wheat, Henan Institute of Science and Technology, Xinxiang 453003, Henan, China
2
Henan Collaborative Innovation Center of Water Security for Water Source Region of Mid-line of South-to-North Diversion Project and Henan Key Laboratory of Ecological Security for Water Source Region of Mid-line of South-to-North Diversion Project, Nanyang Normal University, Nanyang 473061, Henan, China
3
Department of Biology, East Carolina University, Greenville, NC 27858, USA
*
Author to whom correspondence should be addressed.
†
Co-first authors.
Genes 2019, 10(6), 473; https://0-doi-org.brum.beds.ac.uk/10.3390/genes10060473
Submission received: 14 May 2019 / Revised: 15 June 2019 / Accepted: 17 June 2019 / Published: 21 June 2019
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Class III peroxidases (PODs), commonly known as secretable class III plant peroxidases, are plant-specific enzymes that play critical roles in not only plant growth and development but also the responses to biotic and abiotic stress. In this study, we identified 198 nonredundant POD genes, designated GhPODs, with 180 PODs being predicted to secrete into apoplast. These POD genes were divided into 10 sub-groups based on their phylogenetic relationships. We performed systematic bioinformatic analysis of the POD genes, including analysis of gene structures, phylogenetic relationships, and gene expression profiles. The GhPODs are unevenly distributed on both upland cotton sub-genome A and D chromosomes. Additionally, these genes have undergone 15 segmental and 12 tandem duplication events, indicating that both segmental and tandem duplication contributed to the expansion of the POD gene family in upland cotton. Ka/Ks analysis suggested that most duplicated GhPODs experienced negative selection, with limited functional divergence during the duplication events. High-throughput RNA-seq data indicated that most highly expressed genes might play significant roles in root, stem, leaf, and fiber development. Under K or P deficiency conditions, PODs showed different expression patterns in cotton root and leaf. This study provides useful information for further functional analysis of the POD gene family in upland cotton.

1. Introduction

Peroxidases (EC 1.11.1.x) are encoded by multigenic families and are involved in several important physiological and developmental processes. Among them, class III peroxidases (EC 1.11.1.7), belonging to the haem peroxidase subfamily, exist only in plants and have an extremely widespread presence in the plant kingdom [1]. They are members of a large multigenic family with more than 200 members in switchgrass [2], 93 in Populus [3], 138 in rice [4], and 73 in Arabidopsis [5]. The nomenclature of Class III plant peroxidase is not unified and various abbreviations were reported, such as POX [6,7], GPX [8], Prx [5], ClassIIIPRX [2], and POD [9,10]. Hereafter, class III peroxidases were abbreviated as PODs.
PODs are involved in a broad range of physiological processes such as auxin metabolism, lignin and suberin formation, cross-linking of cell wall components, phytoalexin synthesis, defense against biotic or abiotic stress, cell elongation, and the metabolism of reactive nitrogen species and reactive oxygen species (ROS), throughout the plant life cycle from the early stage of germination to the final step of senescence [11,12,13]. It is probably due to their high number of enzymatic isoforms and to the versatility of their enzyme-catalyzed reactions.
PODs are mainly considered as secreted/apoplastic/cell wall proteins, but vacuolar isoforms also exist [14]. Apoplastic POD can be further classified into three major categories based on their chemical and physical association with cell walls and available extraction methods: water soluble and loosely ionically bound, tightly ionically bound, and covalently bound [15]. The specific function of each member of the family is still elusive. Furthermore, they usually show dual enzymatic activities. For example, some consume ROS and others produce ROS; some loosen the cell wall and others stiffen the cell wall [14]. Therefore, they play a pivotal role in cellular growth and response to biotic and abiotic stresses. The comprehensive researches are necessary to explore the role of POD in plant growth and defense.
Widely cultivated in more than 100 countries, cotton is considered one of the most important fiber-producing and economic crops around the world. Suboptimal phosphorus (P) and potassium (K) availability, widely present in agriculture, negatively influences cotton growth and development and reduces cotton fiber yield and quality [16,17]. PODs participate in cotton growth and development, cotton defense against biotic and abiotic stresses, and fiber development. For example, GhPOX1 play an important role during fiber cell elongation possibly mediated by ROS production [6]. The cotton flower-specifically expressed pod, predominantly in pollen, suggested that peroxidase is involved in the male reproductive processes of angiosperms [18]. Two PODs from cotton play a role in the oxidative burst response of cotton to bacterial blight [19]. Up to now, however, no genome-wide characterization of the pod family and their responses to PK deficiency has been performed in cotton. The recently published genome sequence of Gossypium hirsutum L. acc. TM-1, a tetraploid cotton species [20], provides us with a great opportunity to identify and characterize pods in the cotton genome and to explore the expression profiles of pods under PK deficiency conditions.
In the present study, we performed for the first time the comprehensive analysis and responses to PK deficiency of the pod family in G. hirsutum. A total of 198 non-redundant POD encoding genes were identified in the genome of G. hirsutum. and were subsequently subjected to a systematic genomic analysis, including studies on phylogenetic relationships, on chromosome location, on gene duplication status, on substitution rates, on gene structures, on expression profiling and secretion traits, and on responses to PK deficiency in cotton leaf and root. The differentiation of functions of GhPODs were predicted on the basis of the expression profiles of pod members and the phylogenetic analysis among the POD proteins in G. hirsutum. Additionally, we analyzed whether the expansion of the pod family in G. hirsutum was caused by segmental duplication and/or tandem duplication. In summary, our genome-wide analysis of the POD gene family will contribute to future studies on the functional differentiation of POD proteins in different physiological processes of G. hirsutum; the differential responses to PK deficiency will benefit the elucidation of the relationship of physiological processes such as root elongation and branching, leaf senescence, and ROS modulation and with specific POD isoenzymes under PK deficiency.

2. Materials and Methods

2.1. Sequence Retrieval for POD Proteins in Cotton

The local BLAST database was established with protein sequences of upland cotton (G. hirsutum L. acc. TM-1) whole genome (download from http://mascotton.njau.edu.cn). The protein sequences of POD family members in the genome of Arabidopsis were retrieved from the TAIR database (http://www.arabidopsis. org/). The candidate sequences of POD in cotton were acquired by BLASTP with each of the 73 different amino acid sequences of Arabidopsis POD gene family as query sequences (screening threshold value/E-Value: 1e−10). To verify the reliability of the initial results, the acquired candidate sequences were further submitted against PFAM (http://pfam.xfam.org/) to verify the domains for identifying the POD gene family members in cotton. The theoretical molecular weights (MWs) and isoelectric points (pIs) of the proteins were collected through an online program (http://www.ebi.ac.uk/Tools/seqstats/emboss_pepstats/).

2.2. Phylogenetic Analysis

Multiple sequence alignments were conducted on the amino acid sequences of POD proteins in G. hirsutum genomes using Cluster W of MEGA 5.0 software with the default settings [21]. Subsequently, the software was employed to construct an unrooted phylogenetic tree based on alignments using the Neighbor-Joining (NJ) method with the following parameters: model (p-distance), bootstrap (1000 replicates), and gap/missing data (pairwise deletion).

2.3. Gene Structure Analysis

The genomic and CDS sequences of cotton PODs, extracted from G. hirsutum genome databases, were compared by using the Gene Structure Display Server program (http://gsds.cbi.pku.edu.cn/) to infer the exon/intron organization of POD genes.

2.4. Analysis of Chromosomal Location and Gene Duplication

Information about the physical locations of all POD genes on chromosomes was obtained through BLASTn searches against the G. hirsutum genome database. All GhPOD genes were then mapped on the chromosomes using the software MapInspect (http://mapinspect.software.informer.com). The detection of POD gene duplication events was also carried out and paralogous POD gene pairs were identified based on the alignment results. The criteria were as follows: the shorter sequence covers over 80% of the longer sequence after alignment and the minimum identity of aligned regions is equal to or above 80%. In addition, to explore the selection pressures among POD duplicated genes, we calculated the nonsynonymous mutation rate (Ka), synonymous mutation rate (Ks), and Ka/Ks values for the duplicated gene pairs with Mega 5.0.

2.5. Cotton Culture and Expression Analysis of POD Genes under PK Deficiency

Cotton (G. hirsutum L. TM-1) was planted in a growth chamber (day/night of 14/10 h with temperature 30/25 °C and photo intensity 450 µmol/m2·s) under liquid culture. The solution composition was as follows (mmol/L): 2.5 Ca(NO3)2, 1 MgSO4, 0.5 NH4H2PO4, 2.5 KCl, 2 NaCl, 2 × 10−4 CuSO4, 1 × 10−3 ZnSO4, 0.1 EDTA-FeNa, 2 × 10−2 H3BO3, 5 × 10−6 (NH4)6Mo7O24 and 1 × 10−3 MnSO4. The seedlings with four expanded leaves were treated separately with the original solution (control), low K solution (0.05 KCl with NaCl to balance the Cl ion and others the same as in the original solution), and low P solution (0.005 NH4H2PO4 with NH4Cl to balance the NH4+ and others the same as in the original solution). On the 7th day of treatment, the third leaf from the uppermost was counted and the young roots of all treatments were sampled and stored in −80 °C for RNA extraction and gene expression analysis.
Expression profiles of POD genes response to PK deficiency were analyzed by using the Illumina Hiseq2000 (Illumina, San Diego, CA, USA) to perform high-throughput RNA-seq of the root and leaf of control, P deficiency, and K-deficiency. In total, 26.95 Gb of raw RNA-seq data were generated (BGI-Tech., Shenzhen, China). RNA-seq reads were mapped to the cotton genotype TM-1 genome using Tophat (v2.0.8; Top Hat, Toronto, Canada). To measure the gene expression level in sampled tissues, we calculated the expression of each gene using FPKM (Fragments per Kilobase of exon model per Million mapped reads) with Cufflinks (v2.1.1; http://cole-trapnell-lab.github.io/cufflinks/). We analyzed the POD gene expression changes in root and leaf under control, P-deficiency, and K-deficiency by using software MultiExperiment Viewer (MeV; http://mev.tm4.org).

2.6. Localization of POD Proteins

Secretion of POD proteins to the apoplast or to the vacuole were predicted by combinations of using SignalP (www.cbs.dtu.dk/services/SignalP/) with a signal peptide, SecretomeP (www.cbs.dtu.dk/services/SecretomeP) without a signal peptide and TargetP (www.cbs.dtu.dk/services/TargetP). The secreted POD proteins were further investigated in the xylem saps, separately, from field cotton [8] and chamber cotton.

3. Results

3.1. Identification of POD Genes

We used the 73 Arabidopsis POD genes to acquire 264 cotton POD genes by BLASTP and further verified their domains with PFAM. A total of 198 non-redundant POD genes with conserved POD domains were identified in cotton. This number is greater than that in Arabidopsis (73) [5], Populus (93) [3], Chinese Pear (94) [22], maize (119) [23], and rice (138) ([4]; but it was similar with that in switchgrass (200) [2]. For convenience, we assigned names to these POD genes (GhPOD01-198) according to their chromosomal positions. The length of the 198 newly identified POD proteins varies from 160 to 1098 amino acid (aa) with an average of 332 aa. There is only one POD containing more than 672 aa. The isoelectric point (PI) varied from 4.12 to 10.50 with a mean of 7.73 and >7.0 of 67.2% POD proteins. Other information of chromosomal location, molecular weight (MW) gene size, coding sequence (CDs) size of each GhPOD gene/protein is shown in Table 1.

3.2. Phylogenetic Analysis

The phylogenetic tree was constructed with the NJ method based on multiple sequence alignment of entire amino acid sequence of 198 upland cotton POD protein sequences in order to acquire a better understanding of evolutionary history and the phylogenetic relationship of POD in upland cotton. Based on the phylogenetic tree (Figure 1), we identified 10 major clusters with high bootstrap probabilities (BPs) ranging from 59 to 100%, among them six clusters had 100% BPs, one had 95% BPs, and two had 71% or nearly 71% BPs (Figure 1). The POD genes were not evenly distributed in some groups in upland cotton, Cluster I had the most members (68), which could be divided into seven subgroups, but Cluster VII only had two members (Figure 1 and Figure 2).

3.3. Gene Structures

POD genes with classical conserved intron/exon gene structure were observed. The coding sequence of 100 of the 198 peroxidase genes are disrupted by three introns at conserved positions (Figure 2). However, variations in this basic gene structure were observed for another 98 of the family, implicating loss of one or more introns (64) or gain of one or more introns (34). Forty genes lost one of the three putative ancestral introns, while fifteen genes lost two introns. Additionally, nine genes (among them, eight genes were closely related, belonged to VI subgroup) were devoid of any introns. In comparison with the classical three introns, the number of POD genes gaining one to nine more introns (except eight), were 5, 3, 2, 4, 13, 2, 4, and 1, respectively. These differences may be derived from a single intron loss or gain events during the long evolutionary period. In addition, 23 of 24 genes with more than seven introns constitute group X, which contains the largest numbers of introns. Sub-clusters with conserved intron/exon gene structure were also observed.

3.4. Chromosomal Location and Gene Duplication

To investigate the genome organization and distribution of GhPODs on different subgenomes and chromosomes (CH) in upland cotton, a chromosome map was constructed. Among 198 GhPOD genes, 94 and 101 were located in subgenomes A and D, respectively. The other three POD genes were located on scaffolds. GhPOD genes unequally distributes in both subgenome A and D. For subgenome A, CHA5 has the most POD members (15), and ChA1 and CHA4 have the least but same members (3); for subgenome D, CHD12 has the most POD members (12), and CHD6 has the least members (4). In addition, the majority of chromosomes, especially for subgenome A, exhibit a relatively high density of GhPOD genes, which tends to assemble at the bottoms or the tops, such as CHA3, 5, 7, 8, 12 and CHD3, 5, 8, 9.
Gene duplications, including segmental and tandem duplication, are one of the primary driving forces in genome evolution [24]. In this study, 27 duplicated gene pairs were identified (Figure 3, Table 2); among them, 12 and 15 GhPOD gene pairs were in subgenome A and subgenome D, respectively. A total of 15 gene pairs (29 genes) were localized to segmentally duplicated regions, while 9 gene clusters (12 gene pairs or duplication events; 21 genes) are arranged in tandem repeats [Segmental duplication (Gh_D07G1630/Gh_D09G2418) (I); (Gh_D07G1630/Gh_D09G2046) (I); Gh_A08G2028/Gh_A11G3132) (I); Gh_D08G2420/Gh_D11G2183 (I); Gh_D04G0130/Gh_D05G2666 (I); Gh_A05G1577/Gh_A06G0019 (IV); Gh_D03G0246/Gh_D12G0699 (VI); Gh_D04G1101/Gh_D09G1420 (VI); Gh_A04G0639/Gh_A09G1415 (VI); Gh_A05G2945/Gh_A07G2109 (VIII); Gh_A05G3141/Gh_A07G1997 (VIII); Gh_D01G1632/Gh_D02G2245 (X); Gh_D05G3875/Gh_D08G2093 (X); Gh_A01G1388/Gh_A03G1812 (X); Gh_A05G0863/Gh_A08G1744 (X); Tandem duplication (Gh_D09G2046/Gh_D09G2047 (I); Gh_D11G0463/Gh_D11G0607 (I); Gh_D13G0906/Gh_D13G0907/Gh_D13G0909 (IV); Gh_A10G2288/Gh_A10G2290 (IV); Gh_A12G0055/Gh_A12G0056 (IV); Gh_A13G0772/Gh_A13G0773 (IV); Gh_D01G1310/Gh_D01G1317 (VI); Gh_A08G1744/Gh_A08G1745/Gh_A08G1746 (X); Gh_D08G2093/Gh_D08G2094/Gh_D082095 (X)].
To explore the selection pressures among duplicated POD genes, we calculated the Ka, Ks, and Ka/Ks values for 27 identified gene pairs (Table 2). In general, Ka/Ks > 1 indicates positive selection, Ka/Ks = 1 indicates neutral selection, and Ka/Ks < 1 indicates negative selection. The Ka/Ks ratios of most GhPOD gene pairs were <1 except for two pairs (GhPOD55 and GhPOD56; GhPOD56 and GhPOD57) with Ka/Ks = 1, suggesting that these gene pairs were evolved under negative selection in upland cotton.

3.5. Responses of POD Genes to PK Deficiency in Cotton Roots and Leaves

Among 198 GhPOD genes, the expression of 30 genes was not detected in all samples including roots and leaves under control conditions, and under K and P deficient conditions. In roots, K and P deficiency, respectively, induced 10 and 11 POD gene expression from zero to slight level (FPKM value < 1). In leaves, K and P deficiency, respectively, induced 28 and 11 POD gene expression from zero to slight level (FPKM value < 1 for most of them). For roots, compared with controls, 12 and 42 POD gene expression, respectively, was above 2-fold and below 0.5-fold under K deficiency; the expression of 18 POD genes, respectively, was above 2-fold and below 0.5-fold under P deficiency. In leaves, compared with controls, the expression of 25 and 16 POD genes, respectively, was above 2-fold and below 0.5-fold being subjected to K deficiency; the expression of 12 and 33 POD genes, respectively, was above 2-fold and below 0.5-fold being subjected to P deficiency. The same POD gene expressed itself with obviously different patterns in leaves or roots under K or P deficiency. For example, under K deficiency, the expression of Gh_A08G1806 was 0.2 times as much as the controls; however, under P deficiency, it was 6.1 times as much as controls in roots. K deficiency induced the expression of Gh_A08G1806 from zero to 1.3 FPKM, but P deficiency did not change its transcription level.

3.6. Secret Traits of Cotton POD

Among 198 GhPOD genes, 142 POD enzymes were predicted with signal peptide by using SignalP, 147 with signal peptide by using TargetP, another 38 POD enzymes were predicted being secreted into apoplast with SecretomeP among no signal peptides by SignalP prediction (Table 1). In xylem sap, 61 POD enzyme isoforms were identified. Among them, 31 isoforms were found not only in field conditions but also in greenhouse conditions (Figure 4).

4. Discussion

4.1. Identification of Cotton POD Genes and Their Expansion

Members of POD gene family are involved in the regulation of a variety of biological processes. POD proteins are classified into apoplast type and vacuole type [25]. Apoplast type PODs participate in plant cell wall lignification, defense to abiotic and biotic stresses, plant growth and development, etc. The majority of PODs (90%) was predicted to be secreted to apoplast by using SignalP plus SecretomeP, and 61 PODs in total were detected in cotton xylem sap from adult cotton plants in the field and cotton seedlings in the greenhouse, indicating their different roles in cotton growth and development. The predictive tools for localization show very different results, indicating the fact that plant localization signals are very variable.
Previously, Delannoy et al. (2003) characterized nine POD genes, found them showing differential expressions in response to the pathogen and suggested that they may have various functions in cotton defense to bacterial blight disease [9]. Furthermore, Delannoy et al. (2006) analyzed 12 POD genes from cotton and found two of them played a role in the oxidative burst response of cotton to bacterial blight [19]. Mei et al. (2009) also investigated 10 POD genes in cotton fiber development [6].
Systematic and comprehensive analyses of POD gene families have been published for Populus trichocarpa [3], Zea mays [23], Arabidopsis thaliana [5], and Oryza sativa [4]. The genome data of allotetraploid upland cotton [20] provides a useful tool for analysis of the upland cotton POD gene family. In our study, 198 POD genes were identified and characterized in upland cotton. The number in the upland cotton is higher than that in Arabidopsis (73) [5], Poplar (93) [3], Chinese Pear (94) [22], maize (119) [23], and rice (138) [4] and similar to that in switchgrass (more than 200) [2]. This is probably due to the fact that upland cotton and switchgrass are tetraploid with larger genomes containing two-type sub-genomes (respectively, 26 chromosomes, A and D; 18, A and B), and that Arabidopsis, maize and rice were diploid with smaller genomes (respectively, 5 chromosomes; 10; 12). However, this cannot explain the fact that the tetraploid Populus with 19 chromosomes in the genome has only 93 POD isoforms.
Gene duplications are one of the primary driving forces in the evolution of genomes and genetic systems [24]. Certain studies have shown that segmental duplication was largely responsible for the expansion of cotton gene families such as the TCP transcription factors in G. raimondii, YABBY and GhHsp20 in G. hirsutum [26,27,28]. By contrast, tandem duplication has contributed significantly to the expansion of this gene family in poplar [3]. However, for nsLTPs, both tandem and segmental duplication contributed to its expansion in G. arboreum and G. hirsutum, while tandem duplication was the dominant pattern in G. raimondii [29]. Interestingly, in this study, we determined that the number of GhPOD genes involved in segmental duplication and tandem duplication is similar, suggesting that both segmental and tandem duplications were equal contributors to the expansion of the POD gene family in upland cotton. It showed a similar gene duplication with POD genes in maize [23]. The Ka/Ks < 1 of the most GhPOD duplicated pairs showed that negative selection may be largely responsible for maintaining the functions of upland cotton POD enzymes.
Phylogenetic analysis of the GHPOD gene family revealed that the exon/intron structures of these genes are relatively conserved due to one half of the GhPOD genes with the 4 exons/3 introns structures. Similarly, 48 of the 73 (65.8%) peroxidase encoding genes in Arabidopsis consist of 4 exons/3 introns [5], and 38 of 138 (27.5%) in rice constitute this structure. It suggests a common ancestral gene with a classical pattern of 4 exons/3 introns [4]. Many studies have shown that introns were specifically inserted into plants and were retained in the genome during the course of evolution [30]. Another half of GhPOD genes gained or lost one or more introns from the POD coding region in a subfamily specific manner, and what is most serious is that some genes contain no introns. This extreme case for POD genes also exist in Arabidopsis [5], rice [4], maize [23], which might be explained either by the loss of all introns or by the occurrence of a reverse transcription event followed by the integration of the cDNA copy back in the genome, as described in mammals, yeast, and maize [31,32]. It is well known that the structural diversity of genes drives the evolution of multigene families. Also, the differences in these characteristics detected between different subfamilies suggest that upland cotton POD members are functionally diversified.

4.2. Expression Profiles of GhPOD Genes

Gene expression patterns can provide important clues about gene function. We used publicly available [20] (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/?term=PRJNA248163) and our own genome-wide transcripts profiling data from upland cotton tissues as a resource to investigate the expression patterns of GhPODs. Based on the public data, 28 of the 198 identified GhPOD genes were not expressed in leaves, roots, stems, petals or fibers (Figure 4 and Figure 5), indicating their functional loss among all organs. In comparison between our results and the results of public data by using data from normal culture conditions, 56, 35, and 19 of the 198 GhPOD genes, respectively, exhibited no expression in leaves or/and roots (Figure 4 and Figure 5). It indicated that a part of GhPODs are expressed coincidently under different conditions or at different developmental stages. Few POD genes demonstrate tissue or organ specificity. In the Arabidopsis genome, 73 POD genes have been annotated, 65 of which were expressed in various tissues, and only three (AtPrx12, AtPrx62, AtPrx65) identified as specific to roots [33]. In the upland cotton genome, 17 of 198 POD genes were identified as specific to roots. However, only 4 and 1 were expressed in leaves and fibers, respectively (Figure 4 and Figure 5).
PODs are expressed in different patterns when facing different biotic and abiotic stresses [19,23,34,35]. This was also confirmed in roots and leaves when plants were subjected to K or P deficiency. For example, expression level of gene Gh_D10G1158 decreased obviously in roots under K deficiency and was about 30% of the controls; under P deficiency, its expression was increased obviously and was 202% as much as controls. Conversely, its expression level was increased by 20.9-fold in the leaf under K deficiency, but very few of them show expression changes under P deficiency. Additionally, different subfamily genes showed obviously different responses to K or P deficiency. In comparison with the controls, 16.7% and 45.2% of GhPOD genes had, respectively, more than 2 and less than 0.5 times the expression level in the root under K deficit in the I subfamily; 48.0% and 33.3% of GhPOD genes had, respectively, more than 2 and less than 0.5 times the expression level in the root under P deficit in the I subfamily. Interestingly, few GhPOD genes with higher expression level with FPKM > 30 showed changes of more than 2 or less than 0.5 times in leaf or root under K or P deficit, indicating that these genes play important roles in the maintenance of basic plant growth.

Author Contributions

P.D. and Z.Z. designed the experiments and wrote the manuscript, B.Z., G.W. and M.C. contribute substantially in data analysis and interpretation of the data. All authors reviewed and approved the final manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (31571600 and 31271648).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Passardi, F.; Theiler, G.; Zamocky, M.; Cosio, C.; Rouhier, N.; Teixera, F.; Margis-Pinheiro, M.; Ioannidis, V.; Penel, C.; Falquet, L.; et al. PeroxiBase: The peroxidase database. Phytochemistry 2007, 68, 1605–1611. [Google Scholar] [CrossRef] [PubMed]
  2. Moural, T.W.; Lewis, K.M.; Barnaba, C.; Zhu, F.; Palmer, N.A.; Sarath, G.; Scully, E.D.; Jones, J.P.; Sattler, S.E.; Kang, C. Characterization of class III peroxidases from switchgrass. Plant Physiol. 2017, 173, 417–433. [Google Scholar] [CrossRef] [PubMed]
  3. Ren, L.L.; Liu, Y.J.; Liu, H.J.; Qian, T.T.; Qi, L.W.; Wang, X.R.; Zeng, Q.Y. Subcellular relocalization and positive selection play key roles in the retention of duplicate genes of Populus class III peroxidase family. Plant Cell 2014, 26, 2404–2419. [Google Scholar] [CrossRef] [PubMed]
  4. Passardi, F.; Longet, D.; Penel, C.; Dunand, C. The class III peroxidase multigenic family in rice and its evolution in land plants. Phytochemistry 2004, 65, 1879–1893. [Google Scholar] [CrossRef] [PubMed]
  5. Tognolli, M.; Penel, C.; Greppin, H.; Simon, P. Analysis and expression of the class III peroxidase large gene family in Arabidopsis thaliana. Gene 2002, 288, 129–138. [Google Scholar] [CrossRef]
  6. Mei, W.; Qin, Y.; Song, W.; Li, J.; Zhu, Y. Cotton GhPOX1 encoding plant class III peroxidase may be responsible for the high level of reactive oxygen species production that is related to cotton fiber elongation. J. Genet. Genom. 2009, 36, 141–150. [Google Scholar] [CrossRef]
  7. Martinez-Rubio, R.; Acebes, J.L.; Encina, A.; Karkonen, A. Class III peroxidases in cellulose deficient cultured maize cells during cell wall remodeling. Physiol. Plant. 2018, 164, 45–55. [Google Scholar] [CrossRef] [PubMed]
  8. Zhang, Z.; Xin, W.; Wang, S.; Zhang, X.; Dai, H.; Sun, R.; Frazier, T.; Zhang, B.; Wang, Q. Xylem sap in cotton contains proteins that contribute to environmental stress response and cell wall development. Funct. Integr. Genom. 2015, 15, 17–26. [Google Scholar] [CrossRef] [PubMed]
  9. Delannoy, E.; Jalloul l, A.; Assigbetse, K.; Marmey, P.; Geiger, J.P.; Lherminier, J.; Daniel, J.F.; Martinez, C.; Nicole, M. Activity of class III peroxidases in the defense of cotton to bacterial blight. Mol. Plant Microbe Interact. 2003, 16, 1030–1038. [Google Scholar] [CrossRef]
  10. Racz, A.; Hideg, E.; Czegeny, G. Selective responses of class III plant peroxidase isoforms to environmentally relevant UV-B doses. J. Plant Physiol. 2018, 221, 101–106. [Google Scholar] [CrossRef]
  11. Passardi, F.; Cosio, C.; Penel, C.; Dunand, C. Peroxidases have more functions than a Swiss army knife. Plant Cell Rep. 2005, 24, 255–265. [Google Scholar] [CrossRef] [PubMed]
  12. Almagro, L.; Gomez Ros, L.V.; Belchi-Navarro, S.; Bru, R.; Ros Barcelo, A.; Pedreno, M.A. Class III peroxidases in plant defence reactions. J. Exp. Bot. 2009, 60, 377–390. [Google Scholar] [CrossRef] [PubMed]
  13. Marjamaa, K.; Hilden, K.; Kukkola, E.; Lehtonen, M.; Holkeri, H.; Haapaniemi, P.; Koutaniemi, S.; Teeri, T.H.; Fagerstedt, K.; Lundell, T. Cloning, characterization and localization of three novel class III peroxidases in lignifying xylem of Norway spruce (Picea abies). Plant Mol. Biol. 2006, 61, 719–732. [Google Scholar] [CrossRef] [PubMed]
  14. Francoz, E.; Ranocha, P.; Nguyen-Kim, H.; Jamet, E.; Burlat, V.; Dunand, C. Roles of cell wall peroxidases in plant development. Phytochemistry 2015, 112, 15–21. [Google Scholar] [CrossRef] [PubMed]
  15. Zhu, J.; Chen, S.; Alvarez, S.; Asirvatham, V.S.; Schachtman, D.P.; Wu, Y.; Sharp, R.E. Cell wall proteome in the maize primary root elongation zone. I. Extraction and identification of water-soluble and lightly ionically bound proteins. Plant Physiol. 2006, 140, 311–325. [Google Scholar] [CrossRef]
  16. Girma, K.; Teal, R.K.; Freeman, K.W.; Boman, R.K.; Raun, W.R. Cotton lint yield and quality as affected by applications of N, P, and K fertilizers. J. Cotton Sci. 2007, 11, 12–19. [Google Scholar]
  17. Dong, H.; Kong, X.; Li, W.; Tang, W.; Zhang, D. Effects of plant density and nitrogen and potassium fertilization on cotton yield and uptake of major nutrients in two fields with varying fertility. Field Crops Res. 2010, 119, 106–113. [Google Scholar] [CrossRef]
  18. Chen, D.; Ding, Y.; Guo, W.; Zhang, T. Molecular cloning and characterization of a flower-specific class III peroxidase gene in G. hirsutum. Mol. Biol. Rep. 2009, 36, 461–469. [Google Scholar] [CrossRef]
  19. Delannoy, E.; Marmey, P.; Jalloul, A.; Etienne, H.; Nicole, M. Molecular analysis of class III peroxidases from cotton. J. Cotton Sci. 2006, 10, 53–60. [Google Scholar]
  20. Zhang, T.; Hu, Y.; Jiang, W.; Fang, L.; Guan, X.; Chen, J.; Zhang, J.; Saski, C.A.; Scheffler, B.E.; Stelly, D.M.; et al. Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat. Biotechnol. 2015, 33, 531–537. [Google Scholar] [CrossRef]
  21. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed]
  22. Cao, Y.; Han, Y.; Meng, D.; Li, D.; Jin, Q.; Lin, Y.; Cai, Y. Structural, evolutionary, and functional analysis of the class III peroxidase gene family in Chinese Pear (Pyrus bretschneideri). Front. Plant Sci. 2016, 7, 1874. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, Y.; Wang, Q.; Zhao, Y.; Han, G.; Zhu, S. Systematic analysis of maize class III peroxidase gene family reveals a conserved subfamily involved in abiotic stress response. Gene 2015, 566, 95–108. [Google Scholar] [CrossRef]
  24. Moore, R.C.; Purugganan, M.D. The early stages of duplicate gene evolution. Proc. Natl. Acad. Sci. USA 2003, 100, 15682–15687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Matsui, T.; Tabayashi, A.; Iwano, M.; Shinmyo, A.; Kato, K.; Nakayama, H. Activity of the C-terminal-dependent vacuolar sorting signal of horseradish peroxidase C1a is enhanced by its secondary structure. Plant Cell Physiol. 2011, 52, 413–420. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, Z.; Gong, Q.; Wang, L.; Jin, Y.; Xi, J.; Li, Z.; Qin, W.; Yang, Z.; Lu, L.; Chen, Q.; et al. Genome-wide study of YABBY genes in upland cotton and their expression patterns under different stresses. Front. Genet. 2018, 9, 33. [Google Scholar] [CrossRef]
  27. Ma, J.; Wang, Q.; Sun, R.; Xie, F.; Jones, D.C.; Zhang, B. Genome-wide identification and expression analysis of TCP transcription factors in Gossypium raimondii. Sci. Rep. 2014, 4, 6645. [Google Scholar] [CrossRef] [PubMed]
  28. Ma, W.; Zhao, T.; Li, J.; Liu, B.; Fang, L.; Hu, Y.; Zhang, T. Identification and characterization of the GhHsp20 gene family in Gossypium hirsutum. Sci. Rep. 2016, 6, 32517. [Google Scholar] [CrossRef]
  29. Li, F.; Fan, K.; Ma, F.; Yue, E.; Bibi, N.; Wang, M.; Shen, H.; Hasan, M.M.; Wang, X. Genomic identification and comparative expansion analysis of the non-specific lipid transfer protein gene family in Gossypium. Sci. Rep. 2016, 6, 38948. [Google Scholar] [CrossRef]
  30. Rogozin, I.B.; Wolf, Y.I.; Sorokin, A.V.; Mirkin, B.G.; Koonin, E.V. Remarkable interkingdom conservation of intron positions and massive, lineage-specific intron loss and gain in eukaryotic evolution. Curr. Biol. 2003, 13, 1512–1517. [Google Scholar] [CrossRef]
  31. Robertson, H.M. The large srh family of chemoreceptor genes in Caenorhabditis nematodes reveals processes of genome evolution involving large duplications and deletions and intron gains and losses. Genome Res. 2000, 10, 192–203. [Google Scholar] [CrossRef] [PubMed]
  32. Trotman, C.N. Introns-early: Slipping lately? Trends Genet. 1998, 14, 132–134. [Google Scholar] [CrossRef]
  33. Valerio, L.; De Meyer, M.; Penel, C.; Dunand, C. Expression analysis of the Arabidopsis peroxidase multigenic family. Phytochemistry 2004, 65, 1331–1342. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, Z.; Chao, M.; Wang, S.; Bu, J.; Tang, J.; Li, F.; Wang, Q.; Zhang, B. Proteome quantification of cotton xylem sap suggests the mechanisms of potassium-deficiency-induced changes in plant resistance to environmental stresses. Sci. Rep. 2016, 6, 21060. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Shigeto, J.; Tsutsumi, Y. Diverse functions and reactions of class III peroxidases. New Phytol. 2016, 209, 1395–1402. [Google Scholar] [CrossRef]
Genes 10 00473 g001 550
Figure 1. Phylogenetic relationship of the 198 identified upland cotton POD genes. Unrooted tree constructed using MEGA5.02 by the Neighbor-Joining (NJ) method. Bootstrap values (above 50%) from 1000 replicates are indicated at each node. The tree shows 10 major phylogenetic subfamilies (subfamilies I to X).
Figure 1. Phylogenetic relationship of the 198 identified upland cotton POD genes. Unrooted tree constructed using MEGA5.02 by the Neighbor-Joining (NJ) method. Bootstrap values (above 50%) from 1000 replicates are indicated at each node. The tree shows 10 major phylogenetic subfamilies (subfamilies I to X).
Genes 10 00473 g001
Genes 10 00473 g002 550
Figure 2. Exon–intron structures of the 198 identified upland cotton POD genes. Exons and introns are indicated by yellow cylinder bars and black lines, respectively.
Figure 2. Exon–intron structures of the 198 identified upland cotton POD genes. Exons and introns are indicated by yellow cylinder bars and black lines, respectively.
Genes 10 00473 g002
Genes 10 00473 g003 550
Figure 3. Chromosomal location and gene duplication events of 198 POD genes on 26 upland cotton chromosomes. Chromosome types and numbers are indicated at the top of each bar. The scale on the left is in mega-bases. The gene ID on the left side of each chromosome correspond to the approximate locations of each POD gene. The segmentally duplicated genes are connected by dashed blue lines, and the tandemly duplicated gene clusters are marked by red square bracket with dashed line.
Figure 3. Chromosomal location and gene duplication events of 198 POD genes on 26 upland cotton chromosomes. Chromosome types and numbers are indicated at the top of each bar. The scale on the left is in mega-bases. The gene ID on the left side of each chromosome correspond to the approximate locations of each POD gene. The segmentally duplicated genes are connected by dashed blue lines, and the tandemly duplicated gene clusters are marked by red square bracket with dashed line.
Genes 10 00473 g003
Genes 10 00473 g004 550
Figure 4. Responses of gene expression in cotton root and leaf to PK deficiency. Values are log10 FPKM. * the POD protein was detected in the field cotton xylem sap; # the POD protein was detected in the greenhouse cotton.
Figure 4. Responses of gene expression in cotton root and leaf to PK deficiency. Values are log10 FPKM. * the POD protein was detected in the field cotton xylem sap; # the POD protein was detected in the greenhouse cotton.
Genes 10 00473 g004
Genes 10 00473 g005 550
Figure 5. Expression profiles of POD genes across different cotton tissues. Values are log10 FPKM. * the POD protein was detected in the field cotton xylem sap; # the POD protein was detected in the greenhouse cotton. Sources of the samples are as follows: leaf-1 (true leaves, Accession: SRX849561); leaf-2 (leaves at 2-week-old plants, Accession: SRX797901); root (roots at 2-week-old plants, Accession: SRX797899); stem (stems at 2-week-old plants, Accession: SRX797900); petals (petals of mature flowers, Accession: SRX797903); fibers (fibers of 25 days post-anthesis).
Figure 5. Expression profiles of POD genes across different cotton tissues. Values are log10 FPKM. * the POD protein was detected in the field cotton xylem sap; # the POD protein was detected in the greenhouse cotton. Sources of the samples are as follows: leaf-1 (true leaves, Accession: SRX849561); leaf-2 (leaves at 2-week-old plants, Accession: SRX797901); root (roots at 2-week-old plants, Accession: SRX797899); stem (stems at 2-week-old plants, Accession: SRX797900); petals (petals of mature flowers, Accession: SRX797903); fibers (fibers of 25 days post-anthesis).
Genes 10 00473 g005
Table
Table 1. The 198 POD genes identified in cotton and their sequence characteristics and location.
Table 1. The 198 POD genes identified in cotton and their sequence characteristics and location.
Protein Name/IDChr LocationGene/CDS Size (bp)PL (aa)/MW (Kda)/PISignalP/SeretomeP/TargetP
GhPOD01/Gh_A01G1388chrA01:87059677-870668267150/867288/31.86/5.64N/S/-
GhPOD02/Gh_A01G1487chrA01:90340312-903419271616/972323/34.9/5.69S/N/S
GhPOD03/Gh_A01G2012chrA01:175229-1770951867/975324/35.23/9.57S/N/-
GhPOD04/Gh_A02G0542chrA02:8131024-81322671244/1002333/36.13/4.12S/N/S
GhPOD05/Gh_A02G0651chrA02:10404884-104073702487/1029342/37.15/5.16S/N/S
GhPOD06/Gh_A02G0927chrA02:35288344-352903452002/1020339/37.96/8.3S/N/S
GhPOD07/Gh_A02G1203chrA02:69325920-693271301211/984327/35.33/8.2N/N/S
GhPOD08/Gh_A02G1466chrA02:80770130-807728802751/993330/36.02/9.44S/N/S
GhPOD09/Gh_A02G1648chrA02:82884707-828875182812/909302/33.98/9.11N/S/-
GhPOD10/Gh_A02G1663chrA02:82988074-829893011228/954317/34.25/9.27S/N/S
GhPOD11/Gh_A03G0199chrA03:3093785-30959822198/969322/34.57/9.93S/N/S
GhPOD12/Gh_A03G0200chrA03:3102427-31039371511/744247/26.58/5.73N/S/C
GhPOD13/Gh_A03G0944chrA03:60522384-605235641181/1008335/36.76/5.19S/N/S
GhPOD14/Gh_A03G0960chrA03:61923097-6197765254556/1971656/71.66/8.68S/N/S
GhPOD15/Gh_A03G1517chrA03:95777005-957781701166/984327/35.68/9.58S/N/S
GhPOD16/Gh_A03G1519chrA03:95814422-958160621641/723240/25.77/6.99S/N/S
GhPOD17/Gh_A03G1812chrA03:99201877-992059794103/867288/31.94/7.2N/N/-
GhPOD18/Gh_A03G2152chrA03:16442-176061165/984327/35.9/8.94S/N/S
GhPOD19/Gh_A03G2153chrA03:28781-299581178/987328/35.15/7.68S/N/S
GhPOD20/Gh_A04G0639chrA04:44156365-441585232159/990329/37.28/7.27S/N/S
GhPOD21/Gh_A04G0963chrA04:59248781-592499391159/996331/36/8.26S/N/S
GhPOD22/Gh_A04G1453chrA04:7944-100752132/1020339/37.62/7.73N/S/S
GhPOD23/Gh_A05G0093chrA05:1157788-11590201233/963320/34.93/7.34S/N/S
GhPOD24/Gh_A05G0507chrA05:5454946-54562831338/999332/36.14/5.05N/N/S
GhPOD25/Gh_A05G0661chrA05:6938388-695033711950/1014337/36.17/6.73S/N/S
GhPOD26/Gh_A05G0863chrA05:8607073-86089561884/753250/27.59/6.35N/S/-
GhPOD27/Gh_A05G1328chrA05:13602610-136043581749/1053350/38.23/9.51S/N/S
GhPOD28/Gh_A05G1452chrA05:15020810-150223601551/999332/37.66/8.54S/N/S
GhPOD29/Gh_A05G1479chrA05:15196748-151978631116/1011336/36.82/8.25N/S/-
GhPOD30/Gh_A05G1577chrA05:16171654-161729911338/951316/34.19/9.5S/N/S
GhPOD31/Gh_A05G1635chrA05:16892245-168936171373/1050349/37.73/7.7S/N/S
GhPOD32/Gh_A05G2401chrA05:29807925-298092491325/972323/34.96/9.47S/N/S
GhPOD33/Gh_A05G2945chrA05:72277688-722792451558/972323/35.2/4.67S/N/S
GhPOD34/Gh_A05G3141chrA05:81059455-81060420966/966321/35.48/8.86S/N/S
GhPOD35/Gh_A05G3239chrA05:84740299-847418751577/972323/35.83/5.06N/S/S
GhPOD36/Gh_A05G3489chrA05:90428029-904296591631/1014337/36.06/9.57S/N/S
GhPOD37/Gh_A05G3726chrA05:122462-1259283467/1182393/42.73/9.35N/S/C
GhPOD38/Gh_A06G0019chrA06:84546-858071262/954317/34.23/9.28S/N/S
GhPOD39/Gh_A06G0383chrA06:6376778-63809064129/1167388/42.25/7.68N/S/C
GhPOD40/Gh_A06G0929chrA06:38044180-380455721393/1038345/37.15/4.4S/N/S
GhPOD41/Gh_A06G1006chrA06:49554553-495569082356/987328/35.73/8.27S/N/S
GhPOD42/Gh_A06G2046chrA06:45925-489813057/960319/34.42/8.88N/S/C
GhPOD43/Gh_A07G0012chrA07:136360-1376131254/996331/36.45/9.58S/N/S
GhPOD44/Gh_A07G0275chrA07:3399593-34008981306/996331/35.82/4.82S/N/-
GhPOD45/Gh_A07G1997chrA07:76170501-761722041704/1026341/37.65/8.37S/N/S
GhPOD46/Gh_A07G2090chrA07:77469444-774715522109/1101366/40.24/6.74N/S/-
GhPOD47/Gh_A07G2109chrA07:77657753-776588581106/966321/34.93/4.26S/N/S
GhPOD48/Gh_A07G2110chrA07:77659654-776608041151/975324/35.43/7.68S/N/S
GhPOD49/Gh_A08G0345chrA08:4189416-41905031088/1002333/35.9/7.9S/N/S
GhPOD50/Gh_A08G0347chrA08:4289654-4290385732/732243/25.94/7.89N/S/-
GhPOD51/Gh_A08G0711chrA08:17779406-177811561751/999332/35.79/5.71N/S/S
GhPOD52/Gh_A08G0712chrA08:18361761-183675215761/1056351/38.2/4.88S/N/S
GhPOD53/Gh_A08G0714chrA08:18413356-184185315176/1032343/36.85/4.51S/N/S
GhPOD54/Gh_A08G0747chrA08:24039157-240425293373/990329/36.31/10.43S/N/S
GhPOD55/Gh_A08G1744chrA08:97409160-974112192060/753250/27.54/5.6N/N/-
GhPOD56/Gh_A08G1745chrA08:97440700-974426201921/726241/26.72/7.74N/N/-
GhPOD57/Gh_A08G1746chrA08:97442829-974464113583/867288/32.07/6.91N/N/-
GhPOD58/Gh_A08G1806chrA08:98615798-986172171420/960319/34.42/5.68S/N/S
GhPOD59/Gh_A08G1950chrA08:100557586-1005586161031/951316/34.48/9.77S/N/S
GhPOD60/Gh_A08G2028chrA08:101461027-1014624411415/969322/34.96/7.75S/N/S
GhPOD61/Gh_A09G0591chrA09:48054359-480581813823/777258/28.28/8.82N/S/-
GhPOD62/Gh_A09G1202chrA09:63774095-637754481354/1020339/36.01/5.24S/N/S
GhPOD63/Gh_A09G1415chrA09:67412511-674143431833/990329/37.12/6.58S/N/S
GhPOD64/Gh_A09G2334chrA09:93316-944821167/951316/34.58/10.02S/N/S
GhPOD65/Gh_A09G2396chrA09:35011-35983973/906301/33.44/8.46S/N/S
GhPOD66/Gh_A10G0565chrA10:7031511-70327081198/1041346/38.29/6.03S/N/S
GhPOD67/Gh_A10G0810chrA10:16610264-166115181255/993330/35.91/9.88S/N/M
GhPOD68/Gh_A10G1317chrA10:69770818-697727541937/966321/33.76/8.33S/N/S
GhPOD69/Gh_A10G1318chrA10:69892772-698939201149/954317/34.64/9.07S/N/S
GhPOD70/Gh_A10G1537chrA10:84415649-844179602312/972323/35.59/8.64S/N/S
GhPOD71/Gh_A10G1626chrA10:87538433-875444125980/909302/33.3/10.36N/S/-
GhPOD72/Gh_A10G1627chrA10:87544899-875460951197/987328/35.3/9.64S/N/S
GhPOD73/Gh_A10G2288chrA10:1167442-11696472206/984327/35.78/10.03S/N/S
GhPOD74/Gh_A10G2290chrA10:1189614-11909831370/993330/36.03/8.94S/N/S
GhPOD75/Gh_A11G0400chrA11:3718938-37202821345/996331/36.51/9.37S/N/S
GhPOD76/Gh_A11G0523chrA11:4917637-49188701234/927308/34.89/9.73N/S/S
GhPOD77/Gh_A11G1669chrA11:25291980-252932321253/1005334/36.93/8.46S/N/S
GhPOD78/Gh_A11G1859chrA11:45617267-456226825416/930309/34.02/4.72N/N/-
GhPOD79/Gh_A11G2043chrA11:63260806-63261288483/483160/17.1/4.89N/S/-
GhPOD80/Gh_A11G3132chrA11:6967-82081242/966321/34.97/8S/N/S
GhPOD81/Gh_A12G0055chrA12:735887-75665920773/972323/35.31/8.26S/N/S
GhPOD82/Gh_A12G0056chrA12:760206-7614261221/972323/35.56/8.99S/N/S
GhPOD83/Gh_A12G0695chrA12:26988350-26989090741/561186/20.48/6.66N/S/-
GhPOD84/Gh_A12G0795chrA12:45772636-457762923657/1068355/39.81/5.12S/N/S
GhPOD85/Gh_A12G1441chrA12:73250331-732520011671/1023340/37.17/9.36N/N/M
GhPOD86/Gh_A12G1915chrA12:81887117-818884291313/1014337/37.91/6S/N/S
GhPOD87/Gh_A12G2221chrA12:84939580-849466257046/32911096/122.66/7.51N/N/-
GhPOD88/Gh_A12G2370chrA12:86207288-862089881701/972323/35.63/8.88S/N/S
GhPOD89/Gh_A12G2508chrA12:87451770-874527741005/1005334/37.17/7.73S/N/S
GhPOD90/Gh_A12G2622chrA12:49766-510611296/1011336/37.08/5.05S/N/S
GhPOD91/Gh_A13G0039chrA13:393916-3960052090/978325/35.42/9.01S/N/S
GhPOD92/Gh_A13G0772chrA13:33291935-332934561522/975324/34.9/6.65S/N/S
GhPOD93/Gh_A13G0773chrA13:33310901-333121631263/987328/35.42/7.36S/N/S
GhPOD94/Gh_A13G2003chrA13:79484281-794873073027/762253/28/5.19N/N/-
GhPOD95/Gh_D01G0112chrD01:855518-8573681851/975324/35.25/9.56S/N/-
GhPOD96/Gh_D01G1310chrD01:36863811-36864806996/996331/36.9/8.19S/N/S
GhPOD97/Gh_D01G1317chrD01:37193921-37194913993/993330/36.56/7.96S/N/S
GhPOD98/Gh_D01G1359chrD01:40012715-400140051291/969322/34.7/7.93S/N/S
GhPOD99/Gh_D01G1632chrD01:51431707-514370245318/867288/31.78/5.81N/S/-
GhPOD100/Gh_D01G1726chrD01:53949584-539511281545/972323/34.9/5.4S/N/S
GhPOD101/Gh_D02G0606chrD02:8255229-82564741246/1014337/36.66/4.13S/N/S
GhPOD102/Gh_D02G0695chrD02:9863175-98657612587/1116371/40.62/7.7N/S/M
GhPOD103/Gh_D02G1110chrD02:31109209-311111601952/1020339/37.91/8.45S/N/S
GhPOD104/Gh_D02G1327chrD02:43855329-438565081180/1008335/36.76/5.38S/N/S
GhPOD105/Gh_D02G1345chrD02:44875923-448773731451/990329/36.19/9.18S/N/S
GhPOD106/Gh_D02G1346chrD02:44991381-449926201240/984327/35.44/7.72S/N/S
GhPOD107/Gh_D02G1997chrD02:64059836-640653865551/987328/35.24/7.94S/N/S
GhPOD108/Gh_D02G1998chrD02:64081876-640830401165/984327/35.7/9.37S/N/S
GhPOD109/Gh_D02G2245chrD02:66344682-663489424261/867288/31.96/7.2N/N/-
GhPOD110/Gh_D03G0059chrD03:390493-3917191227/954317/34.27/9.27S/N/S
GhPOD111/Gh_D03G0074chrD03:531888-534705 2818/912303/33.96/8.22N/S/-
GhPOD112/Gh_D03G0246chrD03:2618435-26202711837/993330/35.99/9.44S/N/S
GhPOD113/Gh_D03G0641chrD03:19207034-19208004971/699232/25.15/8.87N/N/S
GhPOD114/Gh_D03G1381chrD03:42417094-424205363443/1158385/42.12/7.65N/S/-
GhPOD115/Gh_D03G1382chrD03:42426148-424283822235/969322/34.55/9.72S/N/S
GhPOD116/Gh_D04G0130chrD04:1796590-17982441655/1014337/36.04/9.22S/N/S
GhPOD117/Gh_D04G0735chrD04:15156949-151580961148/972323/35.41/8.64S/N/S
GhPOD118/Gh_D04G1101chrD04:36441991-364441392149/987328/37.19/6.73S/N/S
GhPOD119/Gh_D04G1116chrD04:36784560-367860181459/723240/26.55/6.51N/N/-
GhPOD120/Gh_D04G1506chrD04:46932159-4696997837820/2016671/72.76/7.71S/N/S
GhPOD121/Gh_D04G1593chrD04:48023771-480258562086/1020339/37.4/7.06N/N/S
GhPOD122/Gh_D05G0154chrD05:1545607-15468401234/978325/35.03/7.99S/N/S
GhPOD123/Gh_D05G0626chrD05:5026931-5028271 1341/1002333/36.15/5.05N/N/S
GhPOD124/Gh_D05G0807chrD05:6771535-67728311297/1014337/36.1/6.73S/N/S
GhPOD125/Gh_D05G1498chrD05:13463061-134648141754/1059352/38.37/9.62N/S/C
GhPOD126/Gh_D05G1626chrD05:14640100-146416471548/999332/37.65/8.54S/N/S
GhPOD127/Gh_D05G1755chrD05:15837493-158388251333/951316/34.22/9.47S/N/S
GhPOD128/Gh_D05G1817chrD05:16541131-165425041374/1050349/37.76/7.7S/N/S
GhPOD129/Gh_D05G2244chrD05:21543390-215477474358/1218405/43.82/9.15N/S/C
GhPOD130/Gh_D05G2666chrD05:27862573-278638981326/972323/34.85/8.92S/N/S
GhPOD131/Gh_D05G3256chrD05:51004089-510056271539/951316/33.69/8.5S/N/S
GhPOD132/Gh_D05G3875chrD05:76982-788481867/753250/27.57/6.04N/S/-
GhPOD133/Gh_D06G0413chrD06:5889650-5894261 4612/1353450/49.69/7.3N/S/C
GhPOD134/Gh_D06G1049chrD06:22550424-225536403217/1008335/36.1/9.95N/S/C
GhPOD135/Gh_D06G1170chrD06:28537753-285391521400/1038345/37.22/4.46S/N/S
GhPOD136/Gh_D06G1200chrD06:30337737-303408273091/1050349/37.92/8.39S/N/S
GhPOD137/Gh_D07G0019chrD07:200078-2013301253/996331/36.42/9.11S/N/S
GhPOD138/Gh_D07G0331chrD07:3515650-35169551306/996331/35.86/4.66S/N/S
GhPOD139/Gh_D07G1629chrD07:32319582-32320538957/792263/28.65/9.4N/S/-
GhPOD140/Gh_D07G1630chrD07:32393320-323944311112/939312/34.24/9.89S/N/S
GhPOD141/Gh_D07G2217chrD07:53125200-531269091710/1026341/37.7/8.37S/N/S
GhPOD142/Gh_D07G2320chrD07:54647068-546481671100/966321/35.08/4.35S/N/S
GhPOD143/Gh_D07G2321chrD07:54648867-546500181152/975324/35.39/7.24S/N/S
GhPOD144/Gh_D08G0441chrD08:4640562-4641585 1024/858285/30.9/8.72N/S/C
GhPOD145/Gh_D08G0442chrD08:4670537-46716271091/1005334/35.7/8.55S/N/S
GhPOD146/Gh_D08G0443chrD08:4687890-46889791090/1005334/35.66/8.55S/N/S
GhPOD147/Gh_D08G0829chrD08:13789103-137908201718/999332/35.8/6.5S/N/S
GhPOD148/Gh_D08G0832chrD08:14037040-140383661327/1065354/38.22/4.6S/N/S
GhPOD149/Gh_D08G2093chrD08:59783046-597849951950/741246/27.13/6.04N/S/-
GhPOD150/Gh_D08G2094chrD08:59804611-598066832073/738245/27.03/5.9N/N/-
GhPOD151/Gh_D08G2095chrD08:59806893-598100583166/879292/32.56/6.95N/N/-
GhPOD152/Gh_D08G2167chrD08:61090294-610917171424/960319/34.38/5.68S/N/S
GhPOD153/Gh_D08G2330chrD08:63032857-630338811025/945314/34.22/9.44S/N/S
GhPOD154/Gh_D08G2420chrD08:64145436-641468461411/969322/34.92/7.75S/N/S
GhPOD155/Gh_D09G0590chrD09:28473937-284775623626/990329/36.6/9.74S/N/S
GhPOD156/Gh_D09G1208chrD09:39009612-390109321321/1020339/35.98/4.95S/N/S
GhPOD157/Gh_D09G1420chrD09:41767263-417691011839/987328/36.92/6.28S/N/S
GhPOD158/Gh_D09G1611chrD09:43695508-43696494987/987328/36.13/6.76S/N/S
GhPOD159/Gh_D09G2046chrD09:47806422-478075791158/951316/34.5/9.58S/N/S
GhPOD160/Gh_D09G2047chrD09:47809072-478101961125/885294/31.91/9.35S/N/S
GhPOD161/Gh_D09G2048chrD09:47810407-4782121210806/2019672/73.64/9.25S/N/S
GhPOD162/Gh_D09G2418chrD09:1692-28151124/1026341/37.69/9.36S/N/S
GhPOD163/Gh_D10G0605chrD10:6317691-63188911201/1041346/38.32/6.12S/N/S
GhPOD164/Gh_D10G0951chrD10:12771913-127731861274/993330/36.01/9.99S/N/M
GhPOD165/Gh_D10G1157chrD10:19241093-192422361144/954317/34.55/9.06S/N/S
GhPOD166/Gh_D10G1158chrD10:19299752-193007721021/705234/24.46/8.21N/S/-
GhPOD167/Gh_D10G1643chrD10:45423530-454248911362/993330/36.1/8.98S/N/S
GhPOD168/Gh_D10G1644chrD10:45445937-454481112175/984327/35.72/10.04S/N/S
GhPOD169/Gh_D10G1784chrD10:50389073-503913422270/972323/35.64/8.8S/N/S
GhPOD170/Gh_D10G1880chrD10:52703553-52704762 1210/1041346/37.99/8.78S/N/S
GhPOD171/Gh_D10G1881chrD10:52813459-528149111453/972323/35.59/10.39S/N/S
GhPOD172/Gh_D10G1882chrD10:52815399-528166031205/987328/35.25/9.64S/N/S
GhPOD173/Gh_D11G0463chrD11:3963558-39648991342/996331/36.57/9.22S/N/S
GhPOD174/Gh_D11G0607chrD11:5302020-5303236 1217/978325/36.13/8.95N/S/S
GhPOD175/Gh_D11G1824chrD11:20830683-208319331251/969322/35.55/8.46S/N/S
GhPOD176/Gh_D11G2151chrD11:32727567-327301172551/972323/35.51/4.66S/N/S
GhPOD177/Gh_D11G2183chrD11:34071064-340723041241/966321/34.9/7.74S/N/S
GhPOD178/Gh_D11G2348chrD11:45744028-457450721045/957318/35.11/10.5S/N/M
GhPOD179/Gh_D12G0069chrD12:815867-853649 37783/1764587/64.17/7.56N/S/-
GhPOD180/Gh_D12G0070chrD12:860573-8618411269/972323/35.07/6.5S/N/S
GhPOD181/Gh_D12G0071chrD12:866590-8678671278/1029342/37.73/8.99N/S/S
GhPOD182/Gh_D12G0072chrD12:876401-8774381038/789262/28.85/8.06N/S/-
GhPOD183/Gh_D12G0699chrD12:15282699-152839711273/993330/35.78/8.49S/N/S
GhPOD184/Gh_D12G0788chrD12:22644947-226461011155/1068355/39.65/4.9S/N/S
GhPOD185/Gh_D12G0853chrD12:27599646-276010091364/1050349/38.71/6.76S/N/S
GhPOD186/Gh_D12G1559chrD12:46570546-465722261681/1023340/37.12/9.36N/S/M
GhPOD187/Gh_D12G1577chrD12:46856738-468646947957/1803600/68.11/6.77N/S/-
GhPOD188/Gh_D12G2095chrD12:53969820-539711141295/1002333/37.49/5.48S/N/S
GhPOD189/Gh_D12G2504chrD12:57927732-579294471716/972323/35.48/8.72S/N/S
GhPOD190/Gh_D12G2635chrD12:59085174-590861781005/1005334/37.2/7.99S/N/S
GhPOD191/Gh_D13G0054chrD13:451962-4540352074/978325/35.63/9.66S/N/S
GhPOD192/Gh_D13G0906chrD13:18641774-186433021529/978325/34.95/6.77S/N/S
GhPOD193/Gh_D13G0907chrD13:18659873-186611351263/987328/35.47/7.41S/N/S
GhPOD194/Gh_D13G0909chrD13:18677172-186784341263/987328/35.5/7.41S/N/S
GhPOD195/Gh_D13G2402chrD13:60076810-600797622953/750249/27.31/5.17N/N/-
GhPOD196/Gh_Sca005203G02chrSca:12082-12672591/591196/22.12/8.79N/S/-
GhPOD197/Gh_Sca005268G01chrSca:3248-46621415/960319/34.49/7.93S/N/S
GhPOD198/Gh_Sca007173G01chrSca:8026-8850825/825275/30.21/6.74S/N/S
Chr: chromosome; CDS: coding sequence; PL: peptide length; MW: molecular weight; N: non-secreted into the apoplast; S: secreted into apoplast with signal peptide in SignalP, without signal peptide but with non-classical secretion mode in SeretomeP and with signal peptide in TargetP; C: chloroplast; M: mitochondria; -: unsure location).
Table
Table 2. The Ka, Ks, and Ka/Ks values for the 27 gene pairs.
Table 2. The Ka, Ks, and Ka/Ks values for the 27 gene pairs.
Paralogous PairsKsKaKa/KsDuplicate Type
Gh_A01G1388-Gh_A03G18120.4460.0650.15Segmental
Gh_D03G0246-Gh_D12G06991.0370.0990.10Segmental
Gh_D04G0130-Gh_D05G26660.5730.0860.15Segmental
Gh_D04G1101-Gh_D09G14200.4520.0950.21Segmental
Gh_D05G3875-Gh_D08G20930.3730.030.08Segmental
Gh_D07G1630-Gh_D09G20460.2980.0630.21Segmental
Gh_D07G1630-Gh_D09G24180.2050.0460.22Segmental
Gh_D08G2420-Gh_D11G21830.4790.0660.14Segmental
Gh_A04G0639-Gh_A09G14150.4860.0950.20Segmental
Gh_A05G0863-Gh_A08G17440.1480.0740.50Segmental
Gh_A05G1577-Gh_A06G00190.5080.070.14Segmental
Gh_A05G2945-Gh_A07G21090.790.0960.12Segmental
Gh_A05G3141-Gh_A07G19970.7210.0790.11Segmental
Gh_A08G2028-Gh_A11G31320.4830.0730.15Segmental
Gh_D01G1632-Gh_D02G22450.4520.0630.14Segmental
Gh_D08G2093-Gh_D08G20940.0870.0110.13Tandem
Gh_D08G2094-Gh_D08G20950.0870.0170.20Tandem
Gh_D09G2046-Gh_D09G20470.3550.0460.13Tandem
Gh_D11G0463-Gh_D11G06070.1290.0390.30Tandem
Gh_D13G0906-Gh_D13G09070.2050.070.34Tandem
Gh_D13G0907-Gh_D13G09090.0090.0010.11Tandem
Gh_A08G1744-Gh_A08G17450.0390.0391.00Tandem
Gh_A08G1745-Gh_A08G17460.0330.0331.00Tandem
Gh_A10G2288-Gh_A10G22900.3060.0930.30Tandem
Gh_A12G0055-Gh_A12G00560.2890.0610.21Tandem
Gh_A13G0772-Gh_A13G07730.2050.0760.37Tandem
Gh_D01G1310-Gh_D01G13170.0180.0080.44Tandem

Share and Cite

MDPI and ACS Style

Duan, P.; Wang, G.; Chao, M.; Zhang, Z.; Zhang, B. Genome-Wide Identification and Analysis of Class III Peroxidases in Allotetraploid Cotton (Gossypium hirsutum L.) and their Responses to PK Deficiency. Genes 2019, 10, 473. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10060473

AMA Style

Duan P, Wang G, Chao M, Zhang Z, Zhang B. Genome-Wide Identification and Analysis of Class III Peroxidases in Allotetraploid Cotton (Gossypium hirsutum L.) and their Responses to PK Deficiency. Genes. 2019; 10(6):473. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10060473

Chicago/Turabian Style

Duan, Pengfei, Guo Wang, Maoni Chao, Zhiyong Zhang, and Baohong Zhang. 2019. "Genome-Wide Identification and Analysis of Class III Peroxidases in Allotetraploid Cotton (Gossypium hirsutum L.) and their Responses to PK Deficiency" Genes 10, no. 6: 473. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10060473

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

Article Metrics

Back to TopTop