Next Article in Journal
Characterization and Dye Decolorization Potential of Two Laccases from the Marine-Derived Fungus Pestalotiopsis sp.
Next Article in Special Issue
Phylogenetic, Molecular, and Functional Characterization of PpyCBF Proteins in Asian Pears (Pyrus pyrifolia)
Previous Article in Journal
An In Vitro Study on Mitochondrial Compensatory Response Induced by Gliadin Peptides in Caco-2 Cells
Previous Article in Special Issue
HvPAA1 Encodes a P-Type ATPase, a Novel Gene for Cadmium Accumulation and Tolerance in Barley (Hordeum vulgare L.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 8
10.3390/ijms20081863
ijms-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

Deep Transcriptome Analysis Reveals Reactive Oxygen Species (ROS) Network Evolution, Response to Abiotic Stress, and Regulation of Fiber Development in Cotton

by
Yanchao Xu
1,†,
Richard Odongo Magwanga
1,2,†,
Xiaoyan Cai
1,†,
Zhongli Zhou
1,†,
Xingxing Wang
1,
Yuhong Wang
1,
Zhenmei Zhang
1,
Dingsha Jin
1,
Xinlei Guo
1,
Yangyang Wei
1,3,
Zhenqing Li
1,
Kunbo Wang
1,* and
Fang Liu
1,*
1
State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences (ICR, CAAS), Anyang 455000, China
2
Jaramogi Oginga Odinga University of Science and Technology (JOOUST), School of Biological and Physical Sciences (SPBS), P.O BOX 210-40600, Bondo 210-40600, Kenya
3
Biological and Food Engineering, Anyang Institute of Technology, Anyang 455000, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2019, 20(8), 1863; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20081863
Submission received: 28 February 2019 / Revised: 8 April 2019 / Accepted: 8 April 2019 / Published: 15 April 2019
(This article belongs to the Special Issue Mapping Abiotic Stress-Tolerance Genes in Plants)
Browse Figures
Review Reports Versions Notes

Abstract

:
Reactive oxygen species (ROS) are important molecules in the plant, which are involved in many biological processes, including fiber development and adaptation to abiotic stress in cotton. We carried out transcription analysis to determine the evolution of the ROS genes and analyzed their expression levels in various tissues of cotton plant under abiotic stress conditions. There were 515, 260, and 261 genes of ROS network that were identified in Gossypium hirsutum (AD1 genome), G. arboreum (A genome), and G. raimondii (D genome), respectively. The ROS network genes were found to be distributed in all the cotton chromosomes, but with a tendency of aggregating on either the lower or upper arms of the chromosomes. Moreover, all the cotton ROS network genes were grouped into 17 families as per the phylogenetic tress analysis. A total of 243 gene pairs were orthologous in G. arboreum and G. raimondii. There were 240 gene pairs that were orthologous in G. arboreum, G. raimondii, and G. hirsutum. The synonymous substitution value (Ks) peaks of orthologous gene pairs between the At subgenome and the A progenitor genome (G. arboreum), D subgenome and D progenitor genome (G. raimondii) were 0.004 and 0.015, respectively. The Ks peaks of ROS network orthologous gene pairs between the two progenitor genomes (A and D genomes) and two subgenomes (At and Dt subgenome) were 0.045. The majority of Ka/Ks value of orthologous gene pairs between the A, D genomes and two subgenomes of TM-1 were lower than 1.0. RNA seq. analysis and RT-qPCR validation, showed that, CSD1,2,3,5,6; FSD1,2; MSD1,2; APX3,11; FRO5.6; and RBOH6 played a major role in fiber development while CSD1, APX1, APX2, MDAR1, GPX4-6-7, FER2, RBOH6, RBOH11, and FRO5 were integral for enhancing salt stress in cotton. ROS network-mediated signal pathway enhances the mechanism of fiber development and regulation of abiotic stress in Gossypium. This study will enhance the understanding of ROS network and form the basic foundation in exploring the mechanism of ROS network-involving the fiber development and regulation of abiotic stress in cotton.

Graphical Abstract

1. Introduction

Gossypium is one of the largest and most widely distributed genus with more than 50 species [1,2]. Cotton has a long evolution and domestication history. Gossypium diverged from its relatives approximately 10–15 million years ago (MYA), and evolved into three major diploid lineages, the New World clade, the African-Asian clade and the Australian clade, about 5–10 MYA [3]. The allopolyploid cottons emerged about 1–2 MYA due to an intergenomic hybridization event between A and D genomes [3,4], and were domesticated at least 4000 to 5000 years ago [5]. Both G. arboreum and G. hirsutum have a natural long spinable fiber, although, G. raimondii also generates fiber, its fibers are short and not spinable [6]. During human domestication history of cotton, the central focus was on the fiber quality and quantity [7]. The changes in the environment has resulted in to paradigm shift not only to explore on fiber and fiber attributes, but also on the survival strategies developed by plants, whereas plants try to adapt to various environmental conditions in the course of their evolutionary history. Thus, positive selections of genes were different during evolution and human domestication. Increasing number of genome sequencing and resequencing, mRNA sequencing and phenotype assessment of cotton [8] provided an important resource for studying potential biological mechanism related to abiotic stress tolerance in cotton.
Reactive oxygen species (ROS), including singlet oxygen (1O2), superoxide anion (O2−), hydrogen peroxide (H2O2), and hydroxyl radical (•OH), produced by many of the metabolic pathways in aerobic cells, most likely appeared on Earth together with the first atmospheric oxygen molecules about 2.4–3.8 billion years ago and have been a constant companion of aerobic life [9]. Many biological processes are controlled by ROS, which include but are not limited to programmed cell death, biotic and abiotic stress responses, and plant growth and development [10]. However, ROS have been thought to be toxic and unwanted molecules for a long time, but now are widely known as indispensable molecule for plant lives [11]. ROS level is regulated by an intricate network composed of at least 152 genes in Arabidopsis [12,13,14,15], but ROS regulation, mechanisms and the genes involved is not clear in cotton. The ROS regulation network consists of genes which functions as inductors of enzymes which do scavenge on the ROS molecules thereby reducing their toxicity to the plant cell [12].
Cotton fiber development and abiotic stress responses are controlled by many molecules and complex pathways, including ROS [16]. The evolution of long spinnable fibers in cotton is accompanied by novel expression of genes assisting in the regulation of ROS levels [17]. A calcium sensor, GhCaM7, might modulate ROS production and act as a molecular link between Ca2+ and ROS signal pathways in early fiber development stages [18]. ROS-producing genes, including gene encoding NADPH oxidases/respiratory burst oxidase homologue (RBOH), regulated cell expansion through activation of Ca2+ channels, and are also involved in ROS regulation in the plant cell [19,20]. The detection of these genes found to have an integral role in the ROS regulation and fiber development, clearly shows that ROS-scavenging pathway may also be important in the regulation of fiber development. A cytosolic ascorbate peroxidase (APX) gene, GhAPX1AT/DT, regulating optimal H2O2 level, is the key mechanism regulating fiber elongation and deregulation of the increase in H2O2 results in shorter fiber by initiating secondary cell wall-related gene expression [21]. An ascorbate oxidase (AO) gene, GhAO1, controls cotton fiber elongation via the auxin-mediated signaling pathway [22]. A NAC transcription factor, GhFSN1, acts as a positive regulator in controlling second cell wall formation of cotton fiber [23]. Many of the abiotic stresses cause oxidative stress and genes of ROS network also involve in abiotic stress responses [24]. Furthermore, transcriptome analysis of an upland cotton, G. hirsutum, has provided the basis to elucidate the various mechanisms adopted by plants in dealing with salt stress [25].
In this study, genes in ROS network were identified in cotton. Evolutionary analysis proved that ROS network was a conservative and stable system during Gossypium evolution and domestication. Genes of the ROS network that are involved in fiber development and abiotic stress were detected in this research work. A simple ROS signal pathway was used to describe the potential mechanism of fiber development and regulation of abiotic stress in Gossypium. In summary, this study will enhance the understanding of ROS network and pave the path in studying the mechanism of ROS network-involving fiber development and regulation of abiotic stress in cotton.

2. Results

2.1. ROS Network Genes Identification in Gossypium Species

A total of 515, 260, and 261 ROS candidate proteins were identified in G. hirsutum, G. arboreum, and G. raimondii, respectively (Table S1). The proportions of the ROS proteins obtained in the three cotton species, supports the evolution theory of the tetraploid cotton. The allotetraploid cotton, G. hirsutum harbored twice the number individual numbers of the ROS proteins obtained for either of the two diploid cotton species. In all the ROS proteins obtained from the three cotton species, they were further subdivided into different classes or sub families, in which the highest ROS proteins were observed for peroxidase (PER) with 279 proteins encoding the ROS genes, accounting for 26.9% of all the ROS proteins in the three cotton species. These proteins included 36 superoxide dismutase (SOD), 47 ascorbate peroxidase (APX), 24 monodehydroascorbate reductase (MDAR), 20 dehydroascorbate reductase (DHAR), 12 glutathione reductase (GR), 17 catalase (CAT), 32 glutathione peroxidase (GPX), 21 ferritin (FER), 58 respiratory burst oxidase homologue (RBOH), 27 NADPH-like oxidase (FRO), 19 Ubiquinol oxidase (AOX), 37 peroxiredoxin (PrxR), 224 thioredoxin reductase (Trx), 183 Glutaredoxin (GLR) and 279 peroxidase (PER). Beside NADPH oxidases and NADPH oxidase-like (ROS producer), all of the genes encode ROS-scavenging enzymes. There are 152 genes, so far identified in Arabidopsis. Moreover, in genome-wide identification of the ROS genes in Mizuna plants, a total of 32 ROS genes were identified out of the 22,428 transcript, in which 8258 and 14,170 transcripts were up- and down-regulated, respectively [26]. The high number of the ROS genes in Gossypium, suggested a more complex ROS network in relation to fiber development and enhancing abiotic stress tolerance.

2.2. Chromosome Mapping of the ROS Network Gene Family

The ROS genes were mapped and found to be widely distributed across the entire cotton chromosomes. In tetraploid cotton, G. hirsutum, the ROS genes were mapped in all the 26 chromosomes, in which the highest gene loci density were observed in chromosome At05 and its homeolog chromosome with 46 (18.11%) and 38 (14.28%) ROS genes, respectively, while the lowest gene loci density in chromosome At04 with only 10 ROS genes accounting for 1.94% of all the ROS genes identified in G. hirsutum (Figure 1A). The gene distribution among the diploid cotton species, G. arboreum and G. raimondii were found to be present in each of the 13 chromosomes, in which the highest gene loci were detected in chromosome A10 (35, 12.36%) and D09 (40, 13.89%) in G. arboreum and G. raimondii, respectively (Figure 1B,C). Moreover, the genes mapping in the individual chromosomes showed the tendencies of the genes aggregating on either upper or the lower arms of the chromosomes, similar observation have been previously observed in the distribution of the LEA genes in cotton, with high gene loci densities being observed in chromosomes A10 and D09 [27]. The detection of the high number of the ROS genes within the chromosomes showed that the three sets of chromosomes harbor significant genes with diverse roles in enhancing abiotic stress tolerance in plants. This result suggests that chromosome D09, chromosome Dt05 and its homeolog chromosome At05 play a core role in the evolution and dynamic balance in the ROS genes balance in cotton.

2.3. Phylogenetic Tree, Gene Duplication, and Synteny Analysis of the ROS Network Gene Family

In evaluating the evolution of the ROS network genes, the entire ROS network genes were aligned and phylogenetic tree constructed in order to understand the pattern and nature of their evolution. All the genes were classified into seventeen (17) different clades and each clade was made up of members of a particular family. The 17 clades were superoxide dismutase (SOD), ascorbate peroxidase (APX), monodehydroascorbate reductase (MDAR), dehydroascorbate reductase (DHAR), glutathione reductase (GR), catalase (CAT), glutathione peroxidase (GPX), ferritin (FER), respiratory burst oxidase homologue (RBOH), NADPH-like oxidase (FRO), Ubiquinol oxidase (AOX), peroxiredoxin (PrxR), thioredoxin reductase (Trx), Glutaredoxin (GLR), and peroxidase (PER) (Figure 2A). A total of 683 ROS network genes were found to have undergone duplications, in which 631 genes through segmental duplication while only 52 genes were tandemly duplicated (Table S1). The results were in agreement to previous findings in which a number stress responsive genes evolution were majorly governed by segmental type of gene duplication, for instance the LEA genes [27]. Gene duplication is considered as the main processes which do lead to the expansion of various gene families in organisms, and it is a prevailing feature in plant genomes [28,29]. Generally, a high number of gene copies in a gene family are maintained through large-scale segmental duplication or small-scale tandem duplication during evolution.
Synteny analysis, showed a total of 161, 166, 99, and 105 orthologous genes synteny blocks between At and A, Dt and A, At and D, and Dt and D, respectively (Figure 2B–D). It is interesting to note that more ROS genes were observed between the syntenic blocks among the three cotton species, for instance, at most five genes were found within the syntenic regions between At and A, and the same number of genes were also observed between Dt and A. However, at most 20 genes have been found within the syntenic regions between Dt and D (Table 1 and Table S2). The high number of ROS genes detected within the same syntenic regions showed that this syntenic regions have been highly conserved, and this indicates their integral role they play within the plants. This result supports the principle of asymmetric selection of genes within the genomes [8].

2.4. Orthologous ROS Gene Analysis and Evolution Prediction through Synonymous (Ks) and Non-Synonymous (Ka) Rate of Substituion

Allopolyploid cotton evolved following trans-oceanic dispersal of an A genome diploid to the Americas, where the immigrant underwent hybridization, as female, with a native D genome diploid similar to modern G. raimondii [30]. The formation of ROS network might be earlier than that of Gossypium diverged from its closest relatives and Allopolyploid cottons formation. Phylogenetic analysis (Table S3) revealed that, a total of 241 paired genes, were orthologous among the three cotton species, G. raimondii, G. arboreum and G. hirsutum, found to be either directly or indirectly involved in the ROS network. In the analysis of the evolutionary history of the tetraploid cotton, 243 genes were orthologous between G. arboreum and G. raimondii. G. arboreum contained 14 genes, but not in G. raimondii, while only five specific genes were found to be orthologs between G. arboreum and At subgenome. In the D genome, G. raimondii had 17 genes, but not in G. arboreum and only five specific genes were found to be ortholog between G. raimondii and Dt subgenome (Figure 3A). This result suggested that 243 genes were formed evolved much earlier before the separation of the A and D cotton genomes, which are the G. arboreum and G. raimondii, respectively. Thirty-one (31) genes were found to have evolved before the separations of the two diploid cotton species, and even the emergence of the allotetraploid cotton. three genes were found to have originated from D, and introgressed into the Dt subgenome, of the allotetraploid cotton, while some element of gene loss was evident, this could have been occasioned by either loose of function or chromosomal rearrangement during the evolution of allotetraploid from the two diploid cotton species, A and D cotton genomes (Figure 3B). The synonymous substitution value (Ks) peaks of orthologous gene pairs between At/A and Dt/D was estimated at 0.004 and 0.015, respectively. The distribution of 241 putative ortholog pairs between the three cotton species/genomes (mode Ks = 0.0040–0.015) clearly indicated that the secondary peak for cotton, lies to the left of the ortholog peak, which represents a burst of gene duplications that occurred in cotton after the separation between the two lineages. Similar observation has been made in the analysis of the genome evolution of soybean, which showed that soybean genome also underwent whole genome duplication [31]. From our findings, the ROS genes network seemed to have occurred in tandem with the evolution of tetraploid cotton, approximately 0.8–2.8 million years ago (MYA) [2]. Both synonymous site (Ks) peaks of the ROS orthologous gene pairs between the two cotton progenitors, A and D genomes, and two subgenomes (At and Dt) was 0.045 (Figure 3C,D). Moreover, we estimated the ROS genes network divergence time between the A and D genomes to have occurred approximately 8.6 MYA, in the results was in agreement to the divergence time of the two diploid cotton, G. arboreum and G. raimondii which occurred 5–10 MYA [32]. Additionally, the synonymous (Ks) and non-synonymous (Ka) substitution values ratio (Ka/Ks) were less than one (1) that suggested that the ROS gene network evolution, was primarily governed by purifying selection, similar results were also obtained in the analysis of the 156 paralogous pairs of the LEA genes in cotton, in which over 99% of the orthologous gene pairs had the Ka/Ks ratio values of less than 1 [27].

2.5. Analysis of the Expression Levels of Homeolog and Ortholog Genes ROS Network Genes at 10 DPA and 20 DPA of Cotton Fiber Development

RNA-seq data of the public database (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/) were downloaded and reanalyzed, which showed 192 ROS genes pairs were expressed during cotton fiber development (Tables S3 and S4). Moreover, we analyzed the homeolog and ortholog gene pair’s expression levels at 10 and 20 DPA in the three cotton species, G. hirsutum, G. raimondii, and G. arboreum. It is interesting that the expression levels of genes encoding the SOD proteins, such as the CSD, FSD, and MSD proteins exhibited lower expression levels at 10 and 20 DPA in G. raimondii but were significantly higher in G. arboreum and G. hirsutum. The expression levels of most of the SOD genes were higher at 10 DPA than at 20 DPA (Figure 4A,B). The result showed that the SOD genes are integral at 10 DPA, being a critical stage of fiber initiation and elongation, moreover, at the stage of 10 DPA, ABA inhibits fiber cell initiation and elongation [33]. Moreover, previous studies have shown that cotton annexin proteins; AnxGh1:AAR13288, AnxGh2:AAB67993, AnxGhFx:FJ415173, and AnxGhF:AAC33305 are present in higher amount in fibers of 10 DPA wild-type plants as compared to the fuzzless-lintless mutant [34]. This result indicated that SOD genes were highly expressed in taxa with long fiber (G. arboreum and G. hirsutum), while were either down regulated or partially expressed in cotton species with shorter fibers, such G. raimondii of the D genome. Additionally, the expression level of SOD genes showed up regulation at fiber elongation phase (10 DPA) compared to the secondary cell wall biosynthesis/thickening phase (20 DPA). The results obtained for the expression pattern of the ROS genes, correlated positively to previous findings on the expression patterns of some of the genes known to be highly involved in fiber development in cotton, such as APX, MDAR, GR, GPX, FER, AOX, PrxR, Trx, GLR, and PER [21,22,35]. The SOD gene family and other ROS scavenging genes are believed to be involved in cotton fiber development through the change in their expression level in the course of fiber development period. On the other hand, most of the genes with functional role in the regulation of ROS, including RBOH and FRO genes were up regulated at the secondary cell wall biosynthesis phase (20 DPA) relative to fiber elongation phase (10 DPA) in G. hirsutum, G. arboreum, or G. raimondii. In summary, ROS level might be low at fiber elongation phase, but high at secondary cell wall biosynthesis phase in G. hirsutum and G. arboreum. Compared with G. hirsutum and G. arboreum, although the expression model of ROS production genes was similar in G. raimondii, the expression level of ROS removing genes was lower. Fine regulation of steady-state levels of ROS are essential for proper cell elongation [36], and long and short fiber properties were regulated by ROS removing genes, such as SOD gene family. Moreover, the Ascorbate peroxidase 1 (APX1), has been found to be integral in the process of cotton fiber development, and that APX gene profiling in upland cotton, G. hirsutum has provided valuable information on the redox homeostasis in the process of fiber development [37]. Furthermore, when the cytosolic H2O2-scavenging enzyme, APX1, the whole chloroplastic H2O2-scavenging process stops in A. thaliana and the level H2O2 increases resulting into protein oxidation [38]. Ferric Chelate Reductase-2 (FRO2), Ferric Chelate Reductase-5 (FRO5), and respiratory burst oxidase homologs (RBOH6) orthologous gene pairs, modulation ROS level in fiber development, were identified by GWAS analysis and RNA-seq analysis [8]. Furthermore, the RBOH do modulate the ROS and in turn promote root growth in Arabidopsis [39].

2.6. ROS Network Genes Response to Abiotic Stress

We downloaded data of ROS-gene expression levels from pubic RNA-Seq data of G. hirsutum (TM-1). After analyzing gene expression under cold, heat, dehydration, and salt stress treatment, a total of 135 ROS genes were founded to be upregulated by abiotic stresses (Table S5). Compared with the controls, 92 (68%), 79 (59%), 74 (55%), and 85 (62.9%) genes were differently expressed under cold, heat, dehydration, and salt stress, respectively. Twenty-nine genes were differentially expressed under all four different stress conditions, which included CSD, APX, MDAR, CAT, GPX, FER, RBOH, FRO, AOX, GLR, and PER. The majority of the differentially expressed genes (DEGs) showed more expression changes under cold stress than under other abiotic stresses, suggesting that ROS genes were more sensitive to cold stress. ROS-producing genes were up regulated after cold, heat, dehydration, and salt treatments. FRO5 gene shown higher expression level than other ROS-producing DEGs and was upregulated after 3 and 6 h of post stress exposure. APX1 and MDAR1 genes, showed higher expression level than other ROS-scavenging DEGs, and were upregulated after 3 and 6 h (Figure 5). These results proved that the ROS gene do respond positively towards abiotic stress during the early developmental stages. RBOH6 and APX1 genes were highly expressed in G. hirsutum D subgenome than At subgenome which suggested that some gene pairs of G. hirsutum D subgenome and A subgenome might be playing different roles in cotton growth and development [8].

2.7. Analysis of Morphological and Physiological Changes under Salt Stress

Four semi-wild (Latifolium 130, Latifolium 32, Latifolium 40, and Marie-galante 85) and two cultivated (CRI12 and CRI16) Gossypium hirsutum accessions were evaluated for phenotypic traits under salt treatment. Among them, Marie-galante 85 (M85) was highly tolerant to salinity stress, as evidenced by significantly higher relative stem length, relative root water content, relative stem water content, and lower rate of area of damaged leaves (Table 2; Figure 6). Compared with the controls, proline content was significantly enhanced under salt stress in M85, but was significantly lower in CRI12. These results suggested that M85 is highly tolerant compared to CRI12, to salt stress. The content of MDA was not significantly changed between normal condition and salt stress condition in both CRI12 and M85. However, the concentration level of MDA was relatively lower in the leaf tissues of M85 compared to the level in CRI12. Proline plays an important role in maintaining osmotic balance and protecting some important enzymes in plant cells, thus mainly acts as osmoprotectant when plants are exposed to water deficit condition, the detection of high proline content in the leaf tissues of M85, showed that the genotypes has higher adaptability to salinity compared to CRI12. Therefore, the same with salt (NaCl) stress, cotton is exposed to complex alkali-salt stress could result in osmotic stress. Malondialdehyde (MDA) is the final product of lipid peroxidation, and its content can reflect stress tolerance level of the plant [36]. However, in our study, the content of MDA was not significantly changed after salinity treatment. It suggested that alkali-salt stress has not resulted in lipid peroxidation of M85 and CRI12. Taken together, cotton regulation of complex alkali-salt stress was different from salt (NaCl) stress, and M85 had a higher tolerance to osmotic stress than CRI12. Reactive oxygen species (ROS) balance was also observed. The content of hydrogen peroxide (H2O2) and the activity of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), ascorbate peroxidase (APX), and glutathione reductase (GR) were measured. The content of hydrogen peroxide was significantly enhanced in root and leaf of CRI12 under alkali-salt stress. Hydrogen peroxide was also increased in leaf of M85, but significantly decreased in root of M85 under alkali-salt stress. Furthermore, we found that the activity of POD, CAT, APX, and GR was not significantly changed in roots of M85, suggesting that ROS production and removal was balanced in roots of M85. The content of hydrogen peroxide was higher in leaf than that in root. In addition, the activity of POD was significantly increased in leaf of M85, but POD, CAT, APX, GR, and SOD was not significantly affected. This result suggests that the leaf of M85 was better to remove excessive accumulation of ROS. Additionally, the activity of SOD was significantly increased under alkali-salt stress in root of M85, and the activity of POD and CAT was significantly enhanced in root of CRI12, suggesting that redox enzymes was act on different conditions, flexibly [40].

2.8. RNA-Seq Analysis and ROS Network Response to Alkali-Salt Stress

After Illumina Hiseq2500 sequencing, a total of 39.62–57.16 M clean reads were generated for each sample, in which approximately 85.65%–89.61% of reads were mapped to the reference genome, and 77.79%–82.27% of clean reads were mapped to the unique locations of the reference genome (Table S6). Unique mapped reads were assembled and annotated with Cufflinks and were used for analyzing mechanism of plant response to alkali-salt stress. Based on a cutoff q-value ≤ 0.01 and fold change ≥ 2, different expressed genes (DEGs) were identified. The number of DEGs was higher at 12 h post alkali-salt stress treatment than that of 3 h in both roots and leaves in CRI12 and M85 moreover; the ROS network genes analysis in the same region by principal component analysis showed genes were aggregated at a particular region (Figure 7A,B).
The ROS balance system was found to highly enriched, the finding which was consistent with the previous studies [41]. Moreover, the majority of the DEGs were found to be enriched in response to oxidative stress (GO0006979), oxidation-reduction process (GO0055114), oxidoreductase activity (GO0016717, GO0016491, and GO0046857), and peroxidase activities (GO0004601) (Figure S1). We also observed H2O2 accumulation in leaves and roots under alkali-saline stress. The result shows the importance of ROS networks in acclimation for alkali-saline stress in Gossypium hirsutum. DEGs were also enriched in protein serine/threonine kinase activity (GO0004674). This result supports that the serine/threonine kinases, act as one of signal transduction molecular, is an essential kinase for oxidative stress responses [42]. Genes of ROS network that control ROS levels, encoded most oxidation-reduction proteins. The correlation analysis (R2 > 0.7) between different cotton accessions through orthologous genes expression levels analysis were significantly higher (Figure 7C). Based on the analysis of the various ROS network genes, we further predicted the interlink between fiber and abiotic stress responsive ROS genes (Figure 8). We found that CSD1-6, FSD1-2, MSD 1-2, APX3, APX11, FRO5, FRO6, and RBOH6 have a functional role in enhancing fiber development, more specifically at 10 and 20 DPA, while CSD1, APX1, APX2, MDAR1, GPX4-6-7, FER2, RBOH6, RBOH11, and FRO5 are specifically involved in enhancing abiotic stress tolerance, more salt stress. Finally, ten genes were selected from three different libraries and validated through RT-qPCR, the expression levels were same as per the RNA sequencing data (Figure 7D).

3. Discussion

In plants, many different biology processes are controlled by ROS [43]. For example, fiber development and biotic and abiotic stresses have previously been found to be partially regulated by ROS in cotton [44,45]. During 5–10 MYA, Gossypium split into three major diploid lineages: the New World clade (D genome), the African-Asian (A, B, E, and F genomes) clade, and the Australian clade (C, G, and K genomes) [1]. Allotetraploid Upland cotton was formed approximately 1–2 MYA [46] and domesticated for at least 4000–5000 years ago [47]. During human domestication history, the central focus was the fiber quality and quantity, but in the recent past, researchers have found that fiber and fiber qualities are not only affected by the intrinsic nature of the cotton plant, but greatly influenced by the environmental conditions, and therefore in elucidating the mechanisms of fiber development and improvement of its qualities, multidimensional approach has so far been preferred by plant breeders. The change in environmental condition, occasioned by continued human interference, rainfall has become scarce and very erratic in nature, in addition to the emergence of other forms of abiotic stresses, such as salinity, heat, and cold, among others [48]. The sessile nature of plants, makes them to be much prawn to full effects of abiotic stresses, and thus in the course of plant evolution, various survival strategies have been adopted by plants, more so the evolution of various stress responsive plants transcription factors, in which ROS genes are among one of them. In our study, we found that ROS network was conserved in cotton. The Evolutionary history of ROS genes network in cotton was similar to Gossypium evolution. Thus, ROS network was positively selected during cotton evolution and the result of Ka/Ks analysis also supported this view. Gene order and colinearity of allotetraploid cotton are largely conserved between the A and D subgenomes and the extant D progenitor genome (G. raimondii) [47], and the genes colinearity of D09 (G. raimondii) was conserved with Dt05 and At05 (G. hirsutum). Interestingly, more potential ROS genes were mapped to A10 (35, 12.36%), D09 (40, 13.89%), At05 (46, 18.11%), and Dt05 (38, 14.28%) in G. arboreum, G. raimondii, G. hirsutum A subgenome, and G. hirsutum D subgenome, respectively. It suggested that the location of ROS genes was conserved during the domestication process and D09, Dt05, and At05 play a core role in regulated ROS dynamic balance. Above all, ROS plays an important role in cotton growth and development, because of ROS network genes were positive selection and conserved during evolution and domestication of cotton.
At least 15 gene families (SOD, APX, MDAR, DHAR, GR, CAT, GPX, FER, RBOH, FRO, AOX, PrxR, Trx, GLR, and PER) of the ROS network were identified in Gossypium. SOD, APX, and GPX gene families were identified and analyzed. Our results are consistent with the results of SOD gene analyses in previous studies [37] and SOD genes cover the chromosome at the same locations. However, two genes were annotated as copper chaperone for superoxide dismutase in G. arboreum and G. raimondii. Compared with a previous study, the member of GPX genes family was different because of different genomes (NBI_Gossypium hirsutum vs. BGI_Gossypium hirsutum) used for identifying GPX genes. In previous studies, ROS were found to regulate development, differentiation, redox levels, stress signaling, interactions with other organisms, systemic responses, and cell death in higher plants [49,50,51]. Wild cottons allocate greater resources to stress response pathways, while domestication led to reprogrammed resource allocation toward increased fiber growth, possibly through modulating stress-response network [52]. Modulating reactive oxygen species (ROS) production and scavenging could be regulating cotton fiber elongation [53], abiotic stress [41].
In previous studies, different developmental or environmental signals feed into the ROS signaling network and perturb ROS homeostasis in a compartment-specific or even cell-specific manner [12]. Then, the fluctuant ROS is perceived by receptor proteins (NPR, HSF), redox-sensitive transcription factors, or direct inhibition of phosphatases [54]. Serine/threonine protein kinase (OXI1) is involved in ROS sensing and the activation of mitogen-activated-protein kinases (MAPKs) [55] that plays a central role in ROS signal network [56,57,58]. ROS-producing or -scavenging pathway is activated to regulate ROS levels. Different developmental, metabolic and defense pathways are activated to regulate growth and development. In cotton, several families of transcription factors, such as MYB, WRKY, AP2/ERF, and NAC have been found to be involved in regulating plant response to abiotic stress [59], and some of them are also involved in ROS signal network [60].
Cotton fiber development is a complicated biological process. Many researches demonstrate that fiber development was regulated by ROS-mediated pathway. In the previous study, 119 association signals, 71 SNP associated with yield-related traits, 45 association with fiber qualities, and three associated with resistance to Verticillium wilt, were identified through genome-wide association study, and a total of 43 ROS genes were found [61]. Among the 43 genes, seven genes were within the associated SNP signals region of 200 kb (±100 kb), and one of the 43 genes, the following genes were found to have a direct involvement in the regulation of fiber development, for instance Gh_A13G1654 (GLR34) was at the downstream, and was highly associated to signals A13:75463332 (associated with fiber elongation), A13:75463913 (associated with fiber strength), and A13:75473149 (associated with fiber length); Gh_D05G2413 (FRO2) was at the downstream, at the boll weight association signal D05: 24,174,044 (associated with boll weight). This finding may partially support the speculation that ROS network associated with some different plant traits. ROS is involved in regulating plant cell expansion [62], which suggests ROS may also be involved in regulating cotton fiber development [17]. High concentration of ROS in wild-cotton fiber cells inhibits fiber elongation by regulating secondary-cell-wall synthesis [8]. Moreover, the H2O2, involved in cell-wall relaxation, is necessary molecule for cell elongation, however, higher ROS levels may halt cell elongation [63,64,65]. Compared with fiber elongation (10 DPA), ROS level was regulated by ROS network, was upregulated at secondary cell wall biosynthesis (20 DPA) [21]. It is interesting that the expression levels of the SOD genes and some other ROS-scavenging genes were higher G. hirsutum and G. arboreum, and they are known to have relatively longer fiber compared G. raimondii. Moreover, G. hirsutum and G. arboreum, exhibited lower ROS-scavenging enzymes, which implied that they had the higher ROS-producing ability than G. raimondii. The results therefore showed that perturbed ROS levels as a result of the SOD genes and some other genes could perhaps have a net effect on the regulation of various stages of fiber development in G. hirsutum, G. arboreum, and G. raimondii. In addition, fluctuating levels in ROS levels, caused by different expression level of ROS-scavenging genes and ROS-producing genes at fiber elongation and secondary cell wall biosynthesis phases, showed the integral role played by the ROS genes in fiber development.
Almost all the abiotic stress factors do lead to the development of secondary stress condition such as osmotic and oxidative stresses [66]. RNA-seq data, physiological, and phenotypic index analysis in the two upland cotton genotypes, M85 and CRI12 suggested that the ROS gene network is partially responsible to salinity stress tolerance. The majority of the ROS genes were expressed under abiotic stress, which demonstrated that ROS network was also interconnected with abiotic stress tolerance in cotton [67,68,69]. According to the expression values (FPKM) of the ROS network genes in G. hirsutum TM-1, the expression levels of the homologous genes in under various abiotic stress conditions were: 92 (68%), 79 (59%), 74 (55%), and 85 (62.9%) under cold, heat, dehydration, and salt stress conditions, respectively. A total of 29 genes were differentially expressed under four different stress conditions, in which out of the six genes, three genes were the ROS-producing genes, while the other three were the ROS-scavenging genes.

4. Materials and Methods

4.1. Identification of the ROS Genes in Cotton

By against the complete genomes of Gossypium hirsutum (NAU-NBI Assembly v1.1 & Annotation v1.1), G. arboreum (BGI Assembly v2.0 & Annotation v1.0), and G. raimondii (JGI assembly v2.0 & annotation v2.1) which were downloaded from CottonGen database (https://www.cottongen.org), a total of 152 ROS genes controlling ROS levels in Arabidopsis [12] were used to identify ROS genes in cotton genome using BLASTP using the following criteria [70]: the longest transcript in each gene loci was chosen to represent that locus and the sequences were filtered out when (1) CDS length with < 150 bp, (2) CDSs with percentages of ambiguous nucleotides (‘N’) > 50%, (3) CDSs with internal termination codons, and (4) the CDSs with hits (BLAST identities ≥80%) to RepBase sequences (http://www.girinst.org/repbase/index.html). The conservative domains of candidate genes were obtained from Pfam website (http://pfam.xfam.org/) [71] and NCBI Web CD-Search Tool (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/Structure/bwrpsb/bwrpsb.cgi) [72]. The biophysical characters of the encoded proteins were computed using the ExPASy ProtParam tool (http://us.expasy.org/tools/protparam.html).

4.2. Phylogenetic Tree, Orthologous Gene Pairs and Values of Ka, Ks and Ka/Ks

Full-length sequences of ROS-network proteins were first aligned with the ClustalX program (http://www.clustal.org/clustal2/) [73], then phylogenetic trees were constructed by using neighbor-joining (NJ) method with 1,000 bootstrap replicas [74] and the Poisson model by using MEGA 7.0 software (http://www.megasoftware.net). Meanwhile, the orthologous gene pairs of the ROS network in A, D genomes, At and Dt subgenomes were searched by InParanoid software (http://inparanoid.sbc.su.se/cgi-bin/index.cgi). Additionally, the evolutionary rates Ka, Ks, and Ka/Ks ratio were estimated by KaKs_Calculator package (https://kakscalculator.herokuapp.com/). On the basis of the synonymous substitutions per year (λ) of 2.6 × 10−9 for cotton, the divergent time of ROS orthologous gene pairs were estimated by T = Ks/2λ × 10−6 MYA [47].

4.3. Chromosomal Locations, Conserved Synteny Blocks and Gene Duplication

Chromosomal positions of the ROS genes were obtained from the annotation information of cotton genomes (G. arboreum, G. raimondii and G. hirsutum). The chromosomal distributions of genes were drafted on the cotton chromosomes according to the gene positions by Mapchart software [75]. The conserved synteny blocks between A genome and At subgenome, D genome and Dt genome in cotton were inferred using the MCSCANX program with the default parameters (http://chibba.pgml.uga.edu/mcscan2/) [76]. A local BLASTP program was run on different genomes with parameters (e−10). The syntenic relationships were illustrated with the CIRCOS program (http://circos.ca/) [77]. The duplicate gene classifier program of MCSCANX software was used to identify each ROS network gene duplication with the default parameters (OVERLAP WINDOW: 5). All the ROS network genes were classified into various types of duplications.

4.4. Plant Materials and Salt Stress Treatment

G. hirsutum var. marie-galante 85 (M85) and G. hirsutum cv. CRI12 were selected for the salinity treatment. Semi-wild cotton marie-galante was originally collected from Yucatan, Mexico and cultivated in the National Wild Cotton Nursery (Sanya, Hainan, China). G. hirsutum cv. CRI12 was widely planted in china. The seeds of M85 and CRI12 were first germinated at 28 °C in a 16 h light/8 h dark cycle. Then, seedlings were planted in the normal solutions for three weeks. The similar growing seedlings were selected for salt-alkali stress. The salt-alkali treatment was imposed by simulating the actual saline condition of Xinjiang province, a region commonly known for salt stress [78]. We examined the soil salt composition and simulated a similar environment to impose salt-alkali stress, the composition was made of 2.664 mg:g−1 of calcium chloride (CaCl2), 0.179 mg:g−1 of sodium hydrogen carbonate (NaHCO3), 9.94 mg:g−1 of sodium sulphate (Na2SO4), 0.848 mg:g−1 of potassium sulphate (K2SO4) and finally 3.587 mg:g−1 of hydrated magnesium sulphate (MgSO4·7H2O) [79].

4.5. Measurement of Morphological, Physiological and Biochemical Traits

Samples of roots and leaves that collected at 48 h after salt-alkali stress treatment were used for measuring physiological traits. And, seedlings were also used for measuring phenotype at 48 h after salt-alkali stress treatment. Malondialdehyde (MDA) and Proline (PRO), the most important protective macromolecules in response to salt stress, were measured using the kits (Nanjing Jiancheng Bioengineering Institute). Hydrogen peroxide (H2O2), superoxide dismutase (SOD), peroxidase (POD), ascorbate peroxidase (APX), catalase (CAT) and glutathione reductase (GR), important ROS and ROS-scavenging molecules, were measured using the corresponding assay kits (Nanjing Jiancheng Bioengineering Institute). The root length (RL), stem length (SL), roots fresh weight (RWC) and stem fresh weight (SWC) were measured after 48 h salt stress. The value of SPAD, one of photosynthesis indices, was measured using chlorophyll meter model SPAD-502.

4.6. RNA Extraction, cDNA Library Construction, and RNA-Seq

Roots and leaves were collected at 0, 3, 12 and 48 h time points after salt-alkali stress treatment. Then, collected samples were used for transcriptome sequencing. Total RNA was extracted from each cotton sample using TRlzol Reagent (Life technologies, CA, USA) according to the instruction manual. RNA integrity and concentration were checked using an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA). The mRNAs were isolated by NEBNext Poly (A) mRNA Magnetic Isolation Module (NEB, E7490) by New England BioLabs (Ipswich, MA, USA). The cDNA libraries were constructed by following the manufacturer’s instructions of NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, E7530) and NEBNext Multiplex Oligos for Illumina (NEB, E7500) (Illumina, San Diego, CA, USA). Briefly, the enriched mRNA was fragmented into RNAs with approximately 200 nt, which were used to synthesize the first-strand cDNA and then the second cDNA. The double-stranded cDNAs were performed end-repair/dA-tail and adaptor ligation. The suitable fragments were isolated by Agencourt AMPure XP beads (Beckman Coulter, Inc., Indianapolis, IN, USA) and enriched by PCR amplification. Finally, the constructed cDNA libraries were sequenced on a flow cell using an Illumina HiSeq™ 2500 sequencing platform.

4.7. Expression Analysis

Whole-transcriptome sequencing data for G. hirsutum were obtained from the NCBI Sequence Read Archive (SRA) (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra) to analyze the tissue/organ-specific, stage-specific and stress-induced expression patterns of cotton ROS genes [47]. The details of the SRA are shown in (Tables S7 and S8). The row data of RNA-seq of M85 and CRI12 were separately analyzed and the clean reads were obtained by removing reads containing adapter, reads containing ploy-N and lower quality reads from raw data. At the same time, Q30 [80], GC-content and sequence duplication level of the clean data were calculated. Raw sequences were transformed into clean reads after data processing. These clean reads were then mapped to the reference genome sequence. Only reads with a perfect match or one mismatch were further analyzed and annotated based on the reference genome. Tophat2 [81] software was used to map with reference genome. Base on the reference genome, using Cufflinks software [82], mapped reads had been assembled. Quantification of gene expression levels was estimated by fragments per kilobase of transcript per million fragments mapped (FPKM). Base on mapped reads, using FPKM as the index, each gene was estimated by Cuffquant and Cufform. EBSeq software was used to identify the differential expression genes by Fold change ≥ 2 and FDR ≤ 0.01 (False Discovery Rate). FDR was corrected using Benjamini-Hochberg method by p-value. Finally based on the RNA sequence expression data, ten (10) genes were selected and their expression validated through RT-qPCR in relation to salt-alkali stress. In each, three libraries were chosen and the gene specific primers were designed by using primer premier 5 (Table S9). The RT-qPCR analysis was carried out as described by Magwanga et al [83], with GhActin as the internal control gene.

4.8. Statistical Analysis

The experimental data on phenotype and physiology were statistically analyzed using the SAS version 8.2. Graphic presentations were performed using OrginPro 8.0 program and R packages.

5. Conclusions

Due to the ever-changing environmental conditions, the available land for cultivation is declining as a result of various environmental stress factors such as drought and salinity, among others [84]. The combined effect of drought and salinity has been estimated to cause huge losses in agricultural crops as over 6% of arable lands are saline [85]. Thus non edible crops such as cotton are currently being grown in harsher conditions, in soils with high levels of salt content, but due to the importance of the cotton plant, focus is aimed of developing a more robust and salt tolerant cotton cultivars. When plants are exposed to any form stress, the ROS equilibrium shifts leading to high concentrations which results in toxicity and eventually leads to total plant death [86]. In this research work, we carried out deep transcriptome analysis of the ROS network genes, evaluated their evolution pattern and their expression levels under abiotic stress and fiber development stages. A total of 515, 260, and 261 ROS candidate proteins were identified in G. hirsutum, G. arboreum, and G. raimondii, respectively. The high numbers of the ROS network genes were found to have evolved through segmental gene duplication compared to tandem, a result which was in agreement to evolution pattern of various stress responsive genes such the LEA genes which evolved mainly through segmental type of gene duplication [27]. Moreover, RNA sequencing and RT-qPCR validation revealed that the ROS genes had a putative role in enhancing abiotic stress tolerance in cotton and in turn promote fiber development in cotton. In relation to fiber development, the following genes were found to be highly upregulated, CSD1,2,3,5,6; FSD1,2; MSD1,2; APX3,11; FRO5.6; and RBOH6 while for salt stress tolerance, CSD1, APX1, APX2, MDAR1, GPX4-6-7, FER2, RBOH6, RBOH11; and FRO5 were found to have a significant role. This research lays a significant foundation for future exploration of the ROS genes in developing a more stress tolerant cotton genotypes and with superior fiber quality.

Supplementary Materials

Supplementary Materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/8/1863/s1. Table S1. Characteristics of genes of ROS network in Gossypium hirsutum, G. arboreum, and G. raimondii; Table S2. Synteny blocks of G.hirsutum and diploid cotton (G. arboreum and G. raimondii); Table S3. Orthologous gene pairs of G. hirsutum, G. arboreum, and G. raimondii; Table S4. Expression levels of ROS network genes of G. hirsutum (TM-1), G. arboreum, and G. raimondii fiber; Table S5. Expression levels of ROS network genes of G. hirsutum (TM-1) leaves under abiotic stress; Table S6. Summary of the results of RNA-seq and map to the G. hirsutum genome; Table S7. Gene accession number and samples information of transcriptome data used in our research; Table S8. Genes of ROS network within 2 Mb (±1 Mb) of association signals; Figure S1. GO enrichment analysis of M85 and CRI12 under alkali-salt stress.

Author Contributions

F.L. and K.W. conceived and designed the study. Y.X., R.O.M., Y.W. (Yangyang Wei) and Z.Z. (Zhongli Zhou) performed the experiments. Y.X., R.O.M., X.C., Z.Z. (Zhongli Zhou), X.W., Y.W. (Yuhong Wang), Z.Z. (Zhenmei Zhang), D.J. and X.G. contributed the materials/analysis tools. Y.X. and R.O.M. wrote the manuscript. F.L., K.W., X.C., R.O.M., Z.Z. (Zhongli Zhou), Z.L. and D.J. revised the manuscript. All authors reviewed and approved the final manuscript.

Funding

This research program was financially sponsored by the National Natural Science Foundation of China (31671745, 31601352, and 31530053) and the National key research and development plan (2016YFD0100306, 2016YFD0101401).

Acknowledgments

We deeply indebted to the entire research team for their valuable effort during the period of this research work.

Conflicts of Interest

The authors declare that they have no competing interests.

Abbreviations

ROSReactive oxygen species
SODSuperoxide dismutase
KaSynonymous
MDARMonodehydroascorbate reductase
KsNonsynonymous
DHARDehydroascorbate reductase
MYAMillion years ago
GRGlutathione reductase
RBOHRespiratory burst oxidase homologue
CATCatalase
NADPHNicotinamide adenine dinucleotide phosphate
GPXGlutathione peroxidase
APXAscorbate peroxidase
FERFerritin
AOAscorbate oxidase
FRONADPH-like oxidase
AOXUbiquinol oxidase
PrxRPeroxiredoxin
GLRGlutaredoxin
TrxThioredoxin reductase
PERPeroxidase
SNPSingle nucleotide polymorphism
DPADays post anthesis
GWASGenome-wide association study.

References

  1. Wendel, J.F.; Flagel, L.E.; Adams, K.L. Jeans, genes, and genomes: Cotton as a model for studying polyploidy. In Polyploidy and Genome Evolution; Springer: Berlin/Heidelberg, Germany, 2012; pp. 181–207. ISBN 9783642314421. [Google Scholar]
  2. Bradshaw, J.E. Genome Evolution and Polyploidy. In Plant Breeding: Past, Present and Future; Springer: Dordrecht, The Netherlands, 2016; pp. 233–269. [Google Scholar]
  3. Flagel, L.E.; Wendel, J.F.; Udall, J.A. Duplicate gene evolution, homoeologous recombination, and transcriptome characterization in allopolyploid cotton. BMC Genomics 2012, 13. [Google Scholar] [CrossRef]
  4. Shan, X.; Liu, Z.; Dong, Z.; Wang, Y.; Chen, Y.; Lin, X.; Long, L.; Han, F.; Dong, Y.; Liu, B. Mobilization of the active MITE transposons mPing and Pong in rice by introgression from wild rice (Zizania latifolia Griseb.). Mol. Biol. Evol. 2005, 22, 976–990. [Google Scholar] [CrossRef]
  5. Dos Soares, L.A.A.; Fernandes, P.D.; de Lima, G.S.; Brito, M.E.B.; do Nascimento, R.; Arriel, N.H.C. Physiology and production of naturally-colored cotton under irrigation strategies using salinized water. Pesqui. Agropecu. Bras. 2018, 53, 746–755. [Google Scholar] [CrossRef]
  6. Xu, Z.; Yu, J.; Kohel, R.J.; Percy, R.G.; Beavis, W.D.; Main, D.; Yu, J.Z. Distribution and evolution of cotton fiber development genes in the fibreless Gossypium raimondii genome. Genomics 2015, 106, 61–69. [Google Scholar] [CrossRef] [PubMed]
  7. Haigler, C.H.; Betancur, L.; Stiff, M.R.; Tuttle, J.R. Cotton fiber: A powerful single-cell model for cell wall and cellulose research. Front. Plant Sci. 2012, 3. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, M.; Tu, L.; Lin, M.; Lin, Z.; Wang, P.; Yang, Q.; Ye, Z.; Shen, C.; Li, J.; Zhang, L.; et al. Asymmetric subgenome selection and cis-regulatory divergence during cotton domestication. Nat. Genet. 2017, 49, 579–587. [Google Scholar] [CrossRef]
  9. Mittler, R. ROS Are Good. Trends Plant Sci. 2017, 22, 11–19. [Google Scholar] [CrossRef] [PubMed]
  10. Mullineaux, P.; Karpinski, S. Signal transduction in response to excess light: Getting out of the chloroplast. Curr. Opin. Plant Biol. 2002, 5, 43–48. [Google Scholar] [CrossRef]
  11. Ba, Q.S.; Liu, H.Z.; Wang, J.W.; Che, H.X.; Niu, N.; Wang, J.S.; Ma, S.C.; Zhang, G.S. Relationship between metabolism of reactive oxygen species and chemically induced male sterility in wheat (Triticum aestivum L.). Can. J. Plant Sci. 2013, 93, 675–681. [Google Scholar] [CrossRef]
  12. Mittler, R.; Vanderauwera, S.; Gollery, M.; Van Breusegem, F. Reactive oxygen gene network of plants. Trends Plant Sci. 2004, 9, 490–498. [Google Scholar] [CrossRef]
  13. Perata, P.; Pucciariello, C.; Novi, G.; Parlanti, S.; Banti, V. Reactive Oxygen Species-Driven Transcription in Arabidopsis under Oxygen Deprivation. Plant Physiol. 2012, 159, 184–196. [Google Scholar] [CrossRef]
  14. Leymarie, J.; Vitkauskaité, G.; Hoang, H.H.; Gendreau, E.; Chazoule, V.; Meimoun, P.; Corbineau, F.; El-Maarouf-Bouteau, H.; Bailly, C. Role of reactive oxygen species in the regulation of arabidopsis seed dormancy. Plant Cell Physiol. 2012, 53, 96–106. [Google Scholar] [CrossRef] [PubMed]
  15. Fu, Y.; Guo, C.; Wu, H.; Chen, C. Arginine decarboxylase ADC2 enhances salt tolerance through increasing ROS scavenging enzyme activity in Arabidopsis thaliana. Plant Growth Regul. 2017, 83, 253–263. [Google Scholar] [CrossRef]
  16. Waszczak, C.; Carmody, M.; Kangasjärvi, J. Reactive Oxygen Species in Plant Signaling. Annu. Rev. Plant Biol. 2018, 69. [Google Scholar] [CrossRef] [PubMed]
  17. Hovav, R.; Udall, J.A.; Chaudhary, B.; Hovav, E.; Flagel, L.; Hu, G.; Wendel, J.F. The evolution of spinnable cotton fiber entailed prolonged development and a novel metabolism. PLoS Genet. 2008, 4. [Google Scholar] [CrossRef]
  18. Tang, W.; Tu, L.; Yang, X.; Tan, J.; Deng, F.; Hao, J.; Guo, K.; Lindsey, K.; Zhang, X. The calcium sensor GhCaM7 promotes cotton fiber elongation by modulating reactive oxygen species (ROS) production. New Phytol. 2014, 202, 509–520. [Google Scholar] [CrossRef]
  19. Damiano, S.; Fusco, R.; Morano, A.; de Mizio, M.; Paternò, R.; de Rosa, A.; Spinelli, R.; Amente, S.; Frunzio, R.; Mondola, P.; et al. Reactive oxygen species regulate the levels of dual oxidase (duox1-2) in human neuroblastoma cells. PLoS ONE 2012, 7. [Google Scholar] [CrossRef]
  20. Kwak, H.J.; Liu, P.; Bajrami, B.; Xu, Y.; Park, S.Y.; Nombela-Arrieta, C.; Mondal, S.; Sun, Y.; Zhu, H.; Chai, L.; et al. Myeloid Cell-Derived Reactive Oxygen Species Externally Regulate the Proliferation of Myeloid Progenitors in Emergency Granulopoiesis. Immunity 2015, 42, 159–171. [Google Scholar] [CrossRef]
  21. Guo, K.; Du, X.; Tu, L.; Tang, W.; Wang, P.; Wang, M.; Liu, Z.; Zhang, X. Fibre elongation requires normal redox homeostasis modulated by cytosolic ascorbate peroxidase in cotton (Gossypium hirsutum). J. Exp. Bot. 2016, 67, 3289–3301. [Google Scholar] [CrossRef]
  22. Li, R.; Xin, S.; Tao, C.; Jin, X.; Li, H. Cotton ascorbate oxidase promotes cell growth in cultured tobacco bright yellow-2 cells through generation of apoplast oxidation. Int. J. Mol. Sci. 2017, 18, 1346. [Google Scholar] [CrossRef]
  23. Zhang, J.; Huang, G.Q.; Zou, D.; Yan, J.Q.; Li, Y.; Hu, S.; Li, X.B. The cotton (Gossypium hirsutum) NAC transcription factor (FSN1) as a positive regulator participates in controlling secondary cell wall biosynthesis and modification of fibers. New Phytol. 2018, 217, 625–640. [Google Scholar] [CrossRef]
  24. Watkins, J.M.; Chapman, J.M.; Muday, G.K. Abscisic Acid-Induced Reactive Oxygen Species Are Modulated by Flavonols to Control Stomata Aperture. Plant Physiol. 2017, 175, 1807–1825. [Google Scholar] [CrossRef]
  25. Zhang, B.; Chen, X.; Lu, X.; Shu, N.; Wang, X.; Yang, X.; Wang, S.; Wang, J.; Guo, L.; Wang, D.; et al. Transcriptome Analysis of Gossypium hirsutum L. Reveals Different Mechanisms among NaCl, NaOH and Na2CO3 Stress Tolerance. Sci. Rep. 2018, 8. [Google Scholar] [CrossRef]
  26. Sugimoto, M.; Oono, Y.; Gusev, O.; Matsumoto, T.; Yazawa, T.; Levinskikh, M.A.; Sychev, V.N.; Bingham, G.E.; Wheeler, R.; Hummerick, M. Genome-wide expression analysis of reactive oxygen species gene network in Mizuna plants grown in long-term spaceflight. BMC Plant Biol. 2014, 14. [Google Scholar] [CrossRef]
  27. Magwanga, R.O.; Lu, P.; Kirungu, J.N.; Lu, H.; Wang, X.; Cai, X.; Zhou, Z.; Zhang, Z.; Salih, H.; Wang, K.; et al. Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton. BMC Genet. 2018, 19. [Google Scholar] [CrossRef]
  28. Flagel, L.E.; Wendel, J.F. Gene duplication and evolutionary novelty in plants. New Phytol. 2009, 183, 557–564. [Google Scholar] [CrossRef]
  29. Soltis, P.S.; Burleigh, J.G.; Chanderbali, A.S.; Yoo, M.-J.; Soltis, D.E. Gene and Genome Duplications in Plants. In Evolution after Gene Duplication; Wiley: London, UK, 2011; pp. 269–298. [Google Scholar]
  30. Wang, Q.; Wang, S.; Zhu, X.; Fang, L.; Hu, Y.; Chen, X.; Huang, X.; Du, X.; Chen, S.; Wan, Q.; et al. Genomic insights into divergence and dual domestication of cultivated allotetraploid cottons. Genome Biol. 2017, 18. [Google Scholar] [CrossRef]
  31. Shoemaker, R.C.; Polzin, K.; Labate, J.; Specht, J.; Brummer, E.C.; Olson, T.; Young, N.; Concibido, V.; Wilcox, J.; Tamulonis, J.P.; et al. Genome duplication in soybean (Glycine subgenus soja). Genetics 1996, 144, 329–338. [Google Scholar] [CrossRef]
  32. Small, R.L.; Wendel, J.F. Phylogeny, duplication, and intraspecific variation of Adh sequences in new world diploid cottons (Gossypium L., Malvaceae). Mol. Phylogenet. Evol. 2000, 16, 73–84. [Google Scholar] [CrossRef]
  33. Miao, Q.; Deng, P.; Saha, S.; Jenkins, J.N.; Hsu, C.-Y.; Abdurakhmonov, I.Y.; Buriev, Z.T.; Pepper, A.; Ma, D.-P. Transcriptome Analysis of Ten-DPA Fiber in an Upland Cotton (Gossypium hirsutum) Line with Improved Fiber Traits from Phytochrome A1 RNAi Plants. Am. J. Plant Sci. 2017, 8, 2530–2553. [Google Scholar] [CrossRef]
  34. Huang, Y.; Wang, J.; Zhang, L.; Zuo, K. A Cotton Annexin Protein AnxGb6 Regulates Fiber Elongation through Its Interaction with Actin 1. PLoS ONE 2013, 8. [Google Scholar] [CrossRef]
  35. Li, H.B.; Qin, Y.M.; Pang, Y.; Song, W.Q.; Mei, W.Q.; Zhu, Y.X. A cotton ascorbate peroxidase is involved in hydrogen peroxide homeostasis during fibre cell development. New Phytol. 2007, 175, 462–471. [Google Scholar] [CrossRef]
  36. Cosgrove, D.J. Enzymes and Other Agents That Enhance Cell Wall Extensibility. Annu. Rev. Plant Physiol. Plant Mol. Biol. 2002, 50, 391–417. [Google Scholar] [CrossRef]
  37. Tao, C.; Jin, X.; Zhu, L.; Xie, Q.; Wang, X.; Li, H. Genome-wide investigation and expression profiling of APX gene family in Gossypium hirsutum provide new insights in redox homeostasis maintenance during different fiber development stages. Mol. Genet. Genomics 2018, 293, 685–697. [Google Scholar] [CrossRef]
  38. Davletova, S.; Rizhsky, L.; Liang, H.; Shengqiang, Z.; Oliver, D.J.; Coutu, J.; Shulaev, V.; Schlauch, K.; Mittler, R. Cytosolic Ascorbate Peroxidase 1 Is a Central Component of the Reactive Oxygen Gene Network of Arabidopsis. Plant Cell 2005, 17, 268–281. [Google Scholar] [CrossRef]
  39. Orman-Ligeza, B.; Parizot, B.; de Rycke, R.; Fernandez, A.; Himschoot, E.; Van Breusegem, F.; Bennett, M.J.; Périlleux, C.; Beeckman, T.; Draye, X. RBOH-mediated ROS production facilitates lateral root emergence in Arabidopsis. Development 2016, 143, 3328–3339. [Google Scholar] [CrossRef]
  40. Sarwar, M.B.; Sadique, S.; Hassan, S.; Riaz, S.; Rashid, B.; Bahaeldeen Babiker Mohamed, T.H. Physio-Biochemical and Molecular Responses in Transgenic Cotton under Drought Stress. J. Agric. Sci. 2017, 23, 157–166. [Google Scholar]
  41. Wei, Y.; Xu, Y.; Lu, P.; Wang, X.; Li, Z.; Cai, X.; Zhou, Z.; Wang, Y.; Zhang, Z.; Lin, Z.; et al. Salt stress responsiveness of a wild cotton species (Gossypium klotzschianum) based on transcriptomic analysis. PLoS ONE 2017, 12. [Google Scholar] [CrossRef]
  42. Enomoto, A.; Kido, N.; Ito, M.; Takamatsu, N.; Miyagawa, K. Serine-Threonine Kinase 38 is regulated by Glycogen Synthase Kinase-3 and modulates oxidative stress-induced cell death. Free Radic. Biol. Med. 2012, 52, 507–515. [Google Scholar] [CrossRef]
  43. Vaahtera, L.; Brosché, M.; Wrzaczek, M.; Kangasjärvi, J. Specificity in ROS Signaling and Transcript Signatures. Antioxid. Redox Signal. 2013, 21, 1422–1441. [Google Scholar] [CrossRef]
  44. Jiang, P.; Zhang, X.; Zhu, Y.; Zhu, W.; Xie, H.; Wang, X. Metabolism of reactive oxygen species in cotton cytoplasmic male sterility and its restoration. Plant Cell Rep. 2007, 26, 1627–1634. [Google Scholar] [CrossRef]
  45. Deng, M.H.; Wen, J.F.; Huo, J.L.; Zhu, H.S.; Dai, X.Z.; Zhang, Z.Q.; Zhou, H.; Zou, X.X. Relationship of metabolism of reactive oxygen species with cytoplasmic male sterility in pepper (Capsicum annuum L.). Sci. Hortic. 2012, 134, 232–236. [Google Scholar] [CrossRef]
  46. Senchina, D.S.; Alvarez, I.; Cronn, R.C.; Liu, B.; Rong, J.; Noyes, R.D.; Paterson, A.H.; Wing, R.A.; Wilkins, T.A.; Wendel, J.F. Rate variation among nuclear genes and the age of polyploidy in Gossypium. Mol. Biol. Evol. 2003, 20, 633–643. [Google Scholar] [CrossRef]
  47. 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]
  48. Magwanga, R.; Lu, P.; Kirungu, J.; Diouf, L.; Dong, Q.; Hu, Y.; Cai, X.; Xu, Y.; Hou, Y.; Zhou, Z.; et al. GBS Mapping and Analysis of Genes Conserved between Gossypium tomentosum and Gossypium hirsutum Cotton Cultivars that Respond to Drought Stress at the Seedling Stage of the BC2F2 Generation. Int. J. Mol. Sci. 2018, 19, 1614. [Google Scholar] [CrossRef]
  49. Simon, U.K.; Polanschütz, L.M.; Koffler, B.E.; Zechmann, B. High Resolution Imaging of Temporal and Spatial Changes of Subcellular Ascorbate, Glutathione and H2O2 Distribution during Botrytis cinerea Infection in Arabidopsis. PLoS ONE 2013, 8. [Google Scholar] [CrossRef]
  50. Heino, P.; Kariola, T.; Pennanen, V.; Broberg, M.; Survila, M.; Palva, E.T.; Sipari, N.; Davidsson, P.R. Peroxidase-Generated Apoplastic ROS Impair Cuticle Integrity and Contribute to DAMP-Elicited Defenses. Front. Plant Sci. 2016, 7. [Google Scholar] [CrossRef]
  51. Rossi, F.R.; Krapp, A.R.; Bisaro, F.; Maiale, S.J.; Pieckenstain, F.L.; Carrillo, N. Reactive oxygen species generated in chloroplasts contribute to tobacco leaf infection by the necrotrophic fungus Botrytis cinerea. Plant J. 2017, 92, 761–773. [Google Scholar] [CrossRef] [PubMed]
  52. Yoo, M.J.; Wendel, J.F. Comparative Evolutionary and Developmental Dynamics of the Cotton (Gossypium hirsutum) Fiber Transcriptome. PLoS Genet. 2014, 10. [Google Scholar] [CrossRef]
  53. Li, S.; Zhang, F.; Wu, S.; Jin, X.; Zhang, T.; Guo, W.; Wang, L.; Cheng, C. A Cotton Annexin Affects Fiber Elongation and Secondary Cell Wall Biosynthesis Associated with Ca2+ Influx, ROS Homeostasis, and Actin Filament Reorganization. Plant Physiol. 2016, 171, 1750–1770. [Google Scholar] [CrossRef]
  54. Van Montagu, M.; Van Camp, W.; Vranova, E.; Villarroel, R.; Atichartpongkul, S.; Inze, D. Comprehensive analysis of gene expression in Nicotiana tabacum leaves acclimated to oxidative stress. Proc. Natl. Acad. Sci. USA 2002, 99, 10870–10875. [Google Scholar] [CrossRef]
  55. Rentel, M.C.; Lecourieux, D.; Ouaked, F.; Usher, S.L.; Petersen, L.; Okamoto, H.; Knight, H.; Peck, S.C.; Grierson, C.S.; Hirt, H.; et al. OXI1 kinase is necessary for oxidative burst-mediated signalling in Arabidopsis. Nature 2004, 427, 858–861. [Google Scholar] [CrossRef]
  56. Jia, H.; Hao, L.; Guo, X.; Liu, S.; Yan, Y.; Guo, X. A Raf-like MAPKKK gene, GhRaf19, negatively regulates tolerance to drought and salt and positively regulates resistance to cold stress by modulating reactive oxygen species in cotton. Plant Sci. 2016, 252, 267–281. [Google Scholar] [CrossRef]
  57. Chu, X.; Wang, C.; Chen, X.; Lu, W.; Li, H.; Wang, X.; Hao, L.; Guo, X. The cotton WRKY gene GhWRKY41 positively regulates salt and drought stress tolerance in transgenic Nicotiana benthamiana. PLoS ONE 2015, 10. [Google Scholar] [CrossRef]
  58. Tang, C.; Kang, Z.; Chen, S.; Deng, L.; Chang, D.; Wang, X. TaADF3, an Actin-Depolymerizing Factor, Negatively Modulates Wheat Resistance Against Puccinia striiformis. Front. Plant Sci. 2016, 6. [Google Scholar] [CrossRef]
  59. Shang, H.; Wang, Z.; Zou, C.; Zhang, Z.; Li, W.; Li, J.; Shi, Y.; Gong, W.; Chen, T.; Liu, A.; et al. Comprehensive analysis of NAC transcription factors in diploid Gossypium: Sequence conservation and expression analysis uncover their roles during fiber development. Sci. China Life Sci. 2016, 59, 142–153. [Google Scholar] [CrossRef]
  60. Miao, Y.; Wang, P.; Du, Y.; Song, C.-P.; Zhao, X. The MPK6-ERF6-ROS-Responsive cis-Acting Element7/GCC Box Complex Modulates Oxidative Gene Transcription and the Oxidative Response in Arabidopsis. Plant Physiol. 2013, 161, 1392–1408. [Google Scholar] [CrossRef]
  61. Fang, L.; Wang, Q.; Hu, Y.; Jia, Y.; Chen, J.; Liu, B.; Zhang, Z.; Guan, X.; Chen, S.; Zhou, B.; et al. Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits. Nat. Genet. 2017, 49, 1089–1098. [Google Scholar] [CrossRef]
  62. Apel, K.; Hirt, H. REACTIVE OXYGEN SPECIES: Metabolism, Oxidative Stress, and Signal Transduction. Annu. Rev. Plant Biol. 2004, 55, 373–399. [Google Scholar] [CrossRef] [PubMed]
  63. Rodriguez, A.A. Reactive Oxygen Species in the Elongation Zone of Maize Leaves Are Necessary for Leaf Extension. Plant Physiol. 2002, 129, 1627–1632. [Google Scholar] [CrossRef] [PubMed]
  64. Rodríguez, A.A.; Ramiro Lascano, H.; Bustos, D.; Taleisnik, E. Salinity-induced decrease in NADPH oxidase activity in the maize leaf blade elongation zone. J. Plant Physiol. 2007, 164, 223–230. [Google Scholar] [CrossRef]
  65. Yamaguchi, M.; Sharp, R.E. Complexity and coordination of root growth at low water potentials: Recent advances from transcriptomic and proteomic analyses. Plant Cell Environ. 2010, 33, 590–603. [Google Scholar] [CrossRef]
  66. Shivashankara, K.S.; Pavithra, K.C.; Geetha, G.A. Antioxidant protection mechanism during abiotic stresses. In Abiotic Stress Physiology of Horticultural Crops; Springer: New Delhi, India, 2016; pp. 47–70. ISBN 9788132227250. [Google Scholar]
  67. Lee, S.; Seo, P.J.; Lee, H.J.; Park, C.M. A NAC transcription factor NTL4 promotes reactive oxygen species production during drought-induced leaf senescence in Arabidopsis. Plant J. 2012, 70, 831–844. [Google Scholar] [CrossRef]
  68. Zhou, T.; Yang, X.; Wang, L.; Xu, J.; Zhang, X. GhTZF1 regulates drought stress responses and delays leaf senescence by inhibiting reactive oxygen species accumulation in transgenic Arabidopsis. Plant Mol. Biol. 2014, 85, 163–177. [Google Scholar] [CrossRef]
  69. Mahmood, K.; El-Kereamy, A.; Kim, S.H.; Nambara, E.; Rothstein, S.J. ANAC032 positively regulates age-dependent and stress-induced senescence in arabidopsis thaliana. Plant Cell Physiol. 2016, 57, 2029–2046. [Google Scholar] [CrossRef]
  70. Wang, J.; Wang, J.; Qiu, Q.; Wang, K.; Yang, H.; Abbott, R.J.; Liu, J.; Yu, L.; Ma, H.; Lu, Y.; et al. Erratum: Genomic insights into salt adaptation in a desert poplar. Nat. Commun. 2014, 5. [Google Scholar] [CrossRef]
  71. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 2016, 44, D279–D285. [Google Scholar] [CrossRef]
  72. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017, 45, D200–D203. [Google Scholar] [CrossRef]
  73. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; Mcgettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef]
  74. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  75. Voorrips, R.E. MapChart: Software for the Graphical Presentation of Linkage Maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef]
  76. Wang, Y.; Tang, H.; Debarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40. [Google Scholar] [CrossRef] [PubMed]
  77. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed]
  78. Wang, Y.; Xiao, D.; Li, Y.; Li, X. Soil salinity evolution and its relationship with dynamics of groundwater in the oasis of inland river basins: Case study from the Fubei Region of Xinjiang Province, China. Environ. Monit. Assess. 2008, 140, 291–302. [Google Scholar] [CrossRef] [PubMed]
  79. Yanchao, X.; Yangyang, W.; Zhenqing, L.; Xiaoyan, C.; Yuhong, W.; Xingxing, W.; Zhenmei, Z.; Kunbo, W.; Fang, L.; Zhongli, Z. Integrated Evaluation and the Physiological and Biochemical Responses of Semi-Wild Cotton under Complex Salt-Alkali Stress. Cotton Sci. 2018, 231–241. [Google Scholar] [CrossRef]
  80. Ewing, B.; Hillier, L.D.; Wendl, M.C.; Green, P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8, 175–185. [Google Scholar] [CrossRef]
  81. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [PubMed]
  82. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [PubMed]
  83. Magwanga, R.O.; Lu, P.; Kirungu, J.N.; Dong, Q.; Hu, Y.; Zhou, Z.; Cai, X.; Wang, X.; Hou, Y.; Wang, K.; et al. Cotton Late Embryogenesis Abundant (LEA2) Genes Promote Root Growth and Confers Drought Stress Tolerance in Transgenic Arabidopsis thaliana. G3 Genes Genomes Genet. 2018, 8, 2781–2803. [Google Scholar] [CrossRef]
  84. Yoshino, M. Climate Change and Its Effects on Agriculture. J. Agric. Meteorol. 2011, 46, 251–254. [Google Scholar] [CrossRef]
  85. Diouf, L.; Pan, Z.; He, S.-P.; Gong, W.-F.; Jia, Y.H.; Magwanga, R.O.; Romy, K.R.E.; Rashid, H.; Kirungu, J.N.; Du, X. High-Density Linkage Map Construction and Mapping of Salt-Tolerant QTLs at Seedling Stage in Upland Cotton Using Genotyping by Sequencing (GBS). Int. J. Mol. Sci. 2017, 18, 2622. [Google Scholar] [CrossRef]
  86. Kohli, S.K.; Handa, N.; Gautam, V.; Bali, S.; Sharma, A.; Khanna, K.; Arora, S.; Thukral, A.K.; Ohri, P.; Karpets, Y.V.; et al. ROS signaling in plants under heavy metal stress. In Reactive Oxygen Species and Antioxidant Systems in Plants: Role and Regulation under Abiotic Stress; Springer: Singapore, 2017; pp. 185–214. ISBN 9789811052545. [Google Scholar]
Ijms 20 01863 g001 550
Figure 1. Distribution of reactive oxygen species (ROS) network genes in chromosomes of G. hirsutum, G. arboreum and G. raimondii. (A) Distribution of ROS network genes in chromosomes of G. arboreum. (B) Distribution of ROS network genes in chromosomes of raimondii. (C) Distribution of ROS network genes in chromosomes of G. hirsutum.
Figure 1. Distribution of reactive oxygen species (ROS) network genes in chromosomes of G. hirsutum, G. arboreum and G. raimondii. (A) Distribution of ROS network genes in chromosomes of G. arboreum. (B) Distribution of ROS network genes in chromosomes of raimondii. (C) Distribution of ROS network genes in chromosomes of G. hirsutum.
Ijms 20 01863 g001
Ijms 20 01863 g002 550
Figure 2. Phylogenetic tree and synteny analysis. (A) Phylogenetic tress analysis of all the cotton ROS network genes in G. hirsutum, G. raimondii, and G. arboreum. (B) Synteny analysis between A genome and allotetraploid AD1 genome. (C) Synteny analysis between D genome and allotetraploid AD1 genome. (D) Synteny analysis between A and D genomes. Red, green, yellow, and blue block means A genome, D genome, Allotetraploid cotton (TM-1) A and D subgenomes, respectively. The lines with different color mean genes of different chromosome (A and D genomes).
Figure 2. Phylogenetic tree and synteny analysis. (A) Phylogenetic tress analysis of all the cotton ROS network genes in G. hirsutum, G. raimondii, and G. arboreum. (B) Synteny analysis between A genome and allotetraploid AD1 genome. (C) Synteny analysis between D genome and allotetraploid AD1 genome. (D) Synteny analysis between A and D genomes. Red, green, yellow, and blue block means A genome, D genome, Allotetraploid cotton (TM-1) A and D subgenomes, respectively. The lines with different color mean genes of different chromosome (A and D genomes).
Ijms 20 01863 g002
Ijms 20 01863 g003 550
Figure 3. Orthologous gene pairs and evolution analysis of the ROS genes network in G. arboreum, G. raimondii, and G. hirsutum. (A) Orthologous gene pairs of the ROS genes network. (B) The orthologous genes pairs changes with a temporal depiction of phenomena that characterize polyploid evolution in Gossypium. (C) Distribution of Ks values for orthologous genes pairs between the A genome (G. arboreum), D genome (G. raimondii), and two subgenomes of TM-1(G. hirsutum At subgenome and Dt subgenome). At Dt;‘ t’ indicates tetraploid. Peak values for each comparison are indicated with arrows. (D) Distribution of Ka/Ks value between the A genome, D genome and two subgenomes of TM-1.
Figure 3. Orthologous gene pairs and evolution analysis of the ROS genes network in G. arboreum, G. raimondii, and G. hirsutum. (A) Orthologous gene pairs of the ROS genes network. (B) The orthologous genes pairs changes with a temporal depiction of phenomena that characterize polyploid evolution in Gossypium. (C) Distribution of Ks values for orthologous genes pairs between the A genome (G. arboreum), D genome (G. raimondii), and two subgenomes of TM-1(G. hirsutum At subgenome and Dt subgenome). At Dt;‘ t’ indicates tetraploid. Peak values for each comparison are indicated with arrows. (D) Distribution of Ka/Ks value between the A genome, D genome and two subgenomes of TM-1.
Ijms 20 01863 g003
Ijms 20 01863 g004 550
Figure 4. Expression level obtains from RNA-seq data and was show with RPKM (Reads per Kilobase of exon model per Million mapped reads). (A) The expression level of superoxide dismutase (SOD) genes family. (B) The expression of respiratory burst oxidase homologue (RBOH)and NADPH-like oxidase (FRO).
Figure 4. Expression level obtains from RNA-seq data and was show with RPKM (Reads per Kilobase of exon model per Million mapped reads). (A) The expression level of superoxide dismutase (SOD) genes family. (B) The expression of respiratory burst oxidase homologue (RBOH)and NADPH-like oxidase (FRO).
Ijms 20 01863 g004
Ijms 20 01863 g005 550
Figure 5. Expression of ROS Network genes in G. arboreum, G. hirsutum, and G. raimondii under abiotic stress. Expression level obtained from RNA-seq data and was show with RPKM (Reads per Kilobase of exon model per Million mapped reads). (A) The expression of RBOH6/RBOH11 and FRO5 orthologous gene pairs. (B) The expression level of CSD1, APX1, and MDAR1 orthologous gene pairs.
Figure 5. Expression of ROS Network genes in G. arboreum, G. hirsutum, and G. raimondii under abiotic stress. Expression level obtained from RNA-seq data and was show with RPKM (Reads per Kilobase of exon model per Million mapped reads). (A) The expression of RBOH6/RBOH11 and FRO5 orthologous gene pairs. (B) The expression level of CSD1, APX1, and MDAR1 orthologous gene pairs.
Ijms 20 01863 g005
Ijms 20 01863 g006 550
Figure 6. Complex alkali-salinity stress effects on the growth and the balance between intracellular active oxygen production and the defense system of cotton seedlings. (A) The phenotype of M85 and CRI12 indicated by blue and red arrows, respectively. The relative content of RL (root length), SL (stem length), RWC (root water content), SWC (stem water content), LWC (leaf water content), and SPAD (chlorophyll SPAD value) could visually show the variations among the plants. (B,C) The content of MDA (malondialdehyde) and PRO (proline) were measured at 48 h post to complex alkali-saline stress. (D) The content of hydrogen peroxide (H2O2) in leaves. (E–I) The enzymatic activity of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), ascorbate peroxidase (APX), and glutathione reductase (GR) were measured in leaves and roots, respectively. Error bars represent the standard deviation (SD) of three biological replicates.
Figure 6. Complex alkali-salinity stress effects on the growth and the balance between intracellular active oxygen production and the defense system of cotton seedlings. (A) The phenotype of M85 and CRI12 indicated by blue and red arrows, respectively. The relative content of RL (root length), SL (stem length), RWC (root water content), SWC (stem water content), LWC (leaf water content), and SPAD (chlorophyll SPAD value) could visually show the variations among the plants. (B,C) The content of MDA (malondialdehyde) and PRO (proline) were measured at 48 h post to complex alkali-saline stress. (D) The content of hydrogen peroxide (H2O2) in leaves. (E–I) The enzymatic activity of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), ascorbate peroxidase (APX), and glutathione reductase (GR) were measured in leaves and roots, respectively. Error bars represent the standard deviation (SD) of three biological replicates.
Ijms 20 01863 g006
Ijms 20 01863 g007 550
Figure 7. RNA-Seq analyses and RT-qPCR validation of the ROS network genes in M85 and CRI12 under alkali-salt stress. (A) The number of differentially expressed genes (DEGs) of G. hirsutum races marie-galant (M85) and G. hirsutum cultivar (CRI12) under complex alkali-saline stress at 3, 12, and 48 h. Red block is up-regulated DEGs. Green block is down-regulated DEGs. (B) PCA (principal component analysis) plot of the first two components for genes of ROS network. The dot color scheme represents different cotton accession (TM1, red; M85, yellow; CRI12, blue.). (C) Correlation coefficients of TM-1, CRI12, and M85. (D) RT-qPCR analysis of selected ROS genes, 1, 2, and 3 refers to the libraries, yellow-upregulated, blue: down-regulated, and black: no expression.
Figure 7. RNA-Seq analyses and RT-qPCR validation of the ROS network genes in M85 and CRI12 under alkali-salt stress. (A) The number of differentially expressed genes (DEGs) of G. hirsutum races marie-galant (M85) and G. hirsutum cultivar (CRI12) under complex alkali-saline stress at 3, 12, and 48 h. Red block is up-regulated DEGs. Green block is down-regulated DEGs. (B) PCA (principal component analysis) plot of the first two components for genes of ROS network. The dot color scheme represents different cotton accession (TM1, red; M85, yellow; CRI12, blue.). (C) Correlation coefficients of TM-1, CRI12, and M85. (D) RT-qPCR analysis of selected ROS genes, 1, 2, and 3 refers to the libraries, yellow-upregulated, blue: down-regulated, and black: no expression.
Ijms 20 01863 g007
Ijms 20 01863 g008 550
Figure 8. ROS network mediated signal pathway to illuminate the mechanism of fiber development and regulation of abiotic stress in Gossypium. Red box indicates modulation of ROS level by the reactive oxygen gene network of plants. ROS network regulated fiber development and abiotic stress by control ROS level, but it is not clear about how to do it. RNA-seq data of TM-1, CRI12 and M85 proved that serine/threonine kinase (OXI), and motigen-activated protein kinase (MAPK) cascades might involve in the alkali-salt stress (one of abiotic stresses) response which ROS network mediate. Under fiber development-box, heatmap of some ROS network genes indicted genes expression levels at 10 and 20 DPA (A, G. arboreum; D, G. raimondii; At, G. hirsutum A subgenome; Dt, G. hirsutum D subgenome) by R package (Pheatmap, scale = “row”). Under abiotic stress-box, heatmap indicated value of genes expression levels fold change (log2 (FC)) under cold, hot, dehydration and salt stress by Pheatmap package.
Figure 8. ROS network mediated signal pathway to illuminate the mechanism of fiber development and regulation of abiotic stress in Gossypium. Red box indicates modulation of ROS level by the reactive oxygen gene network of plants. ROS network regulated fiber development and abiotic stress by control ROS level, but it is not clear about how to do it. RNA-seq data of TM-1, CRI12 and M85 proved that serine/threonine kinase (OXI), and motigen-activated protein kinase (MAPK) cascades might involve in the alkali-salt stress (one of abiotic stresses) response which ROS network mediate. Under fiber development-box, heatmap of some ROS network genes indicted genes expression levels at 10 and 20 DPA (A, G. arboreum; D, G. raimondii; At, G. hirsutum A subgenome; Dt, G. hirsutum D subgenome) by R package (Pheatmap, scale = “row”). Under abiotic stress-box, heatmap indicated value of genes expression levels fold change (log2 (FC)) under cold, hot, dehydration and salt stress by Pheatmap package.
Ijms 20 01863 g008
Table
Table 1. The number of synteny blocks which including 1–20 gene number, respectively.
Table 1. The number of synteny blocks which including 1–20 gene number, respectively.
Number of Synteny Block
NGEBAt vs. AAAt vs. DDDt vs. AADt vs. DD
1136 *6314058
219101913
331148
42723
51714
60303
70102
80202
90002
100100
110002
150000
160001
200001
NGEB: number of genes per synteny block; * for example, there are 136 synteny blocks which include one ROS gene between At subgenome and A genome.
Table
Table 2. Phenotype of six G. hirsutum accessions at 48 h post to complex alkali-saline stress.
Table 2. Phenotype of six G. hirsutum accessions at 48 h post to complex alkali-saline stress.
AccessionRLSLArea of Damaged Blades
First Leaf *Second Leaf *Third Leaf *Fourth Leaf *
Latifolium 1300.90 a0.80 c0.97 a0.73 a0.65 a0.00 a
Latifolium 320.92 a0.97 a,b0.17 b0.03 b0.00 b0.00 a
Latifolium 400.87 a0.89 b,c0.40 b0.20 b0.23 a,b0.03 a
Marie-galante 850.92 a1.02 a0.32 b0.03 b0.00 b0.00 a
CRI 120.98 a0.68 d0.38 b0.08 b0.02 b0.00 a
CRI 161.12 a0.94 a,b0.38 b0.10 b0.05 b0.05 a
Letters following the number represents significance at 0.05 probability. Different letters a, b, c and d in a column mean statistically significant difference at the 5% level. * Represents the value of leaves damaged rate (Area of damaged leaf/Area of leaf). First leaf means that the expansion of the first true leaf.

Share and Cite

MDPI and ACS Style

Xu, Y.; Magwanga, R.O.; Cai, X.; Zhou, Z.; Wang, X.; Wang, Y.; Zhang, Z.; Jin, D.; Guo, X.; Wei, Y.; et al. Deep Transcriptome Analysis Reveals Reactive Oxygen Species (ROS) Network Evolution, Response to Abiotic Stress, and Regulation of Fiber Development in Cotton. Int. J. Mol. Sci. 2019, 20, 1863. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20081863

AMA Style

Xu Y, Magwanga RO, Cai X, Zhou Z, Wang X, Wang Y, Zhang Z, Jin D, Guo X, Wei Y, et al. Deep Transcriptome Analysis Reveals Reactive Oxygen Species (ROS) Network Evolution, Response to Abiotic Stress, and Regulation of Fiber Development in Cotton. International Journal of Molecular Sciences. 2019; 20(8):1863. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20081863

Chicago/Turabian Style

Xu, Yanchao, Richard Odongo Magwanga, Xiaoyan Cai, Zhongli Zhou, Xingxing Wang, Yuhong Wang, Zhenmei Zhang, Dingsha Jin, Xinlei Guo, Yangyang Wei, and et al. 2019. "Deep Transcriptome Analysis Reveals Reactive Oxygen Species (ROS) Network Evolution, Response to Abiotic Stress, and Regulation of Fiber Development in Cotton" International Journal of Molecular Sciences 20, no. 8: 1863. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20081863

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