Search

Establishing a sustainable industrial structure system is imperative for transforming the society to high-quality development and achieve the ''carbon peak'' and ''carbon neutrality'' goals. In this study, we explored a low-carbon pathway of regional industrial structure optimization to achieve a balance among economic growth, social security, and environmental conservation. A multi-objective optimization model of industrial structure of economic-social-environmental system was constructed by minimizing carbon emissions and maximizing economic benefits and social employment. In this study, we explored the scientific allocation path for industrial structure and set various scenarios to meet the diverse preferences and development requirements of decision-makers. Southwest China was used as a practical case to validate the model's efficacy, offering a decision-making framework for the sustainable development of the regional economy. After optimization, carbon emissions reduced by 0.9 %∼3.2 %, economic benefits increased by 2.5 %∼4.8 %, and employment levels nearly doubled. The proportion of the transportation industry in southwest China decreased by 11.58 %. The raw coal was the energy source with the highest contribution to carbon emissions, 2494.2 × 108 kgCO2eq.

Abstract The disaster caused by landslide is huge. To prevent the spread of the disaster to the maximum extent, it is particularly important to carry out landslide disaster treatment work. The selection of landslide disaster treatment alternative is a large scale group decision-making (LSGDM) problem. Because of the wide application of social media, a large number of experts and the public can participate in decision-making process, which is conducive to improving the efficiency and correctness of decision-making. A IF-TW-LSGDM method based on three-way decision (TWD) and intuitionistic fuzzy set (IFS) is proposed and applied to the selection of landslide treatment alternatives. First of all, considering that experts and the public participate in the evaluation of LSGDM events, respectively, the method of obtaining and handling the public evaluation information is given, and the information fusion approach of the public and experts evaluation information is given. Second, evaluation values represented by fuzzy numbers are converted into intuitionistic fuzzy numbers (IFNs), and the intuitionistic fuzzy evaluation decision matrix described by IFNs is obtained. Then, a new LSGDM method of alternatives classification and ranking based on IFS and TWD is proposed, the calculation steps and algorithm description are given. In this process, we first cluster the experts, then consider the identification and management of non-cooperative behavior of expert groups. This work provides an effective method based on LSGDM for the selection of landslide treatment alternatives. Finally, the sensitivity of parameters is analyzed, and the feasibility and effectiveness of this method are compared and verified.

Abstract Background The cotton industry suffers significant yield losses annually due to Verticillium wilt, which is considered the most destructive disease affecting the crop. However, the precise mechanisms behind this disease in cotton remain largely unexplored. Methods Our approach involved utilizing transcriptome data from G. australe which was exposed to Verticillium dahliae infection. From this data, we identified ethylene-responsive factors and further investigated their potential role in resistance through functional validations via Virus-induced gene silencing (VIGS) in cotton and overexpression in Arabidopsis. Results A total of 23 ethylene response factors (ERFs) were identified and their expression was analyzed at different time intervals (24 h, 48 h, and 72 h post-inoculation). Among them, GauERF105 was selected based on qRT-PCR expression analysis for further investigation. To demonstrate the significance of GauERF105, VIGS was utilized, revealing that suppressing GauERF105 leads to more severe infections in cotton plants compared to the wild-type. Additionally, the silenced plants exhibited reduced lignin deposition in the stems compared to the WT plants, indicating that the silencing of GauERF105 also impacts lignin content. The overexpression of GauERF105 in Arabidopsis confirmed its pivotal role in conferring resistance against Verticillium dahliae infection. Our results suggest that WT possesses higher levels of the oxidative stress markers MDA and H2O2 as compared to the overexpressed lines. In contrast, the activities of the antioxidant enzymes SOD and POD were higher in the overexpressed lines compared to the WT. Furthermore, DAB and trypan staining of the overexpressed lines suggested a greater impact of the disease in the wild-type compared to the transgenic lines. Conclusions Our findings provide confirmation that GauERF105 is a crucial candidate in the defense mechanism of cotton against Verticillium dahliae invasion, and plays a pivotal role in this process. These results have the potential to facilitate the development of germplasm resistance in cotton.

Abstract Background The cyclic nucleotide-gated channel (CNGC) gene family plays a significant role in the uptake of both essential and toxic cations, and has a role in enhancing tolerance to various forms of abiotic stresses as well as the modulation of the heavy metal toxicity to plant through the absorption of heavy metals. Results A complete genome-wide identification and functional characterization of the cotton CNGC genes was carried out, in which 55, 28, and 29 CNGC genes were identified in Gossypium hirsutum, G. raimondii, and G. arboreum, respectively. The protein encoded by the CNGC genes exhibited GRAVY value below zero, indicating their hydrophilic property. CNGC genes were unevenly distributed in 19 out of 26 chromosomes, in which the highest density were observed on Ah05, with 8 genes. High gene coverage was observed among the diploid cotton species, with CNGC genes mapped on all A chromosomes and on 11 out of 13 of D chromosomes. The majority of CNGC proteins were localized in the endoplasmic reticulum, nucleus, and plasma membrane. Gene expression analysis revealed the up-regulation of Gh_A01G0520 (CNGC4) and Gh_D13G1974 (CNGC5) across various forms of abiotic stresses. Moreover, down-regulation of Gh_A01G0520 (CNGC4) and Gh_D13G1974 (CNGC5) in CNGC s silenced plants caused the significantly reduced ability to tolerate drought and salt stresses. All CNGCs silenced plants were recorded to have significantly low content of antioxidants but relatively higher content of oxidant, including MDA and H2O2. Furthermore, SPAD, CMS (cell membrane stability), ELWL (excised leaf water loss), SDW (shoot dry matter weight), and RDW (root dry matter weight) were all lower in CNGCs silenced plants compared with the wild type plants. Conclusion Significant reduction in antioxidant content and negative effects of physiological and morphological characters in CNGCs silenced plants has revealed the novel role of CNGC genes in enhancing cell integrity under abiotic stress conditions. These results provide vital information that will expand our understanding of the CNGC gene family in cotton and other plants, thus promoting the integration of these genes in the development of the environmental resilient plants.

Summary: The 5′-deoxyadenosine deaminase (DADD), a member of the amidohydrolase family regulates biological purine metabolism. In this study, bioinformatic analysis, overexpression and knockdown of GhdadD gene were detected to identify its potential role in drought and salt stress tolerance. The results revealed that GhdadD was induced by ABA, Auxin, MBS and light responsive elements. In transgenic Arabidopsis, seed germination rate and root length were increased under drought or salt stress. GhdadD overexpressed seedlings resulted in higher plant height, less leaf damage and lower ion permeability. The expression of osmotic stress and ABA-responsive genes were up regulated. While in GhdadD-silenced cotton seedlings, CAT, SOD activity and soluble sugar content were reduced, MDA content was increased, and the stoma opening was depressed under drought or salt stress. Some osmics stress marker genes were also up regulated. These data indicating that GhdadD enhanced plant resistance to drought and salt stress through ABA pathways.

Abstract Background Plant pollen has diverse morphological characteristics that can be consistently passed down from generation to generation. Information on pollen morphology is thus immensely important for plant classification and identification. In the genus Gossypium, however, in-depth research on pollen morphology is lacking, with only few reports on limited cotton species. To evaluate the diversity of pollen in Gossypium, we therefore conducted a comprehensive analysis of the pollen morphology of 33 cotton species and varieties using scanning electron microscopy. Results The 33 analyzed cotton samples exhibited common pollen morphological features, including spherical shapes, radial symmetry, echination, panporation, and operculation, while the pollen size, spine shape, spine density and length showed distinctive features. Pollen size varied significantly among species, with diameters ranging from 62.43 μm in G. harknessii to 103.41 μm in G. barbadense. The exine had an echinate sculptural texture, and spines were mostly conical or sharply conical but occasionally rod-like. Spine density varied from 173 in G. incanum to 54 in G. gossypioides, while spine length ranged from 3.53 μm in G. herbaceum to 9.47 μm in G. barbadense. In addition, the 33 cotton species and varieties were grouped at a genetic distance of 3.83 into three clusters. Cluster I comprised five allotetraploid AD-genome cotton species, four D-genome species, and one K-genome species. Cluster II included 13 diploid species from A, B, D, E, and G genomes, whereas Cluster III only consisted one E-genome species G. incanum. Conclusions Although pollen characteristics alone are not enough to resolve taxonomic and systematic relationships within the genus Gossypium, our results add to knowledge on palynomorphology and contribute to phenological information on these taxa. Our findings should aid future systematic and phylogenetic studies of the Gossypium genus.

Abstract Geological hazards caused by strong earthquakes have caused continuous social and economic losses and destruction of the ecological environment in the hazard area, and are mostly manifested in the areas with frequent occurrence of geological hazards or the clustering of geological hazards. Considering the long-term nature of earthquakes and geological disasters in this region, this paper takes ten earthquake-stricken areas in Wenchuan earthquake zone as examples to collect shallow landslide data in 2010, combined with the spatial location of landslides and other factors. Kernel density estimation (KDE) method is used to analyze the spatial characteristics of shallow landslide. Taking the space of shallow landslide as the characteristic variable and fully considering the regulating factors of earthquake-induced landslide: terrain complexity, distance to river, distance to fault, distance to road, lithology, normalized vegetation difference index (NDVI) and ground peak acceleration (PGA) as independent variables, based on KDE and polynomial logistic regression (MLR), A quantitative model of shallow landslide in the earthquake area is constructed. The results show that: (1) PGA has the greatest impact on landslide in the study area. (2) Compared with the two-category logistic regression (two-category LR) model, the susceptibility map of landslide prediction results based on the KDE-MLR landslide susceptibility prediction model is more consistent with the actual situation. (3) The prediction accuracy of the model validation set is 70.7%, indicating that the landslide susceptibility prediction model based on KDE-MLR can effectively highlight the spatial characteristics of shallow landslides in 10 extreme disaster areas. The research results can provide decision-making basis for shallow landslide warning and post-disaster reconstruction in earthquake-stricken areas.

Abstract Background Current climate change scenarios are posing greater threats to the growth and development of plants. Thus, significant efforts are required that can mitigate the negative effects of drought on the cotton plant. GDSL esterase/lipases can offer an imperative role in plant development and stress tolerance. However, thesystematic and functional roles of the GDSL gene family, particularly in cotton under water deficit conditions have not yet been explored. Results In this study, 103, 103, 99, 198, 203, 239, 249, and 215 GDSL proteins were identified in eight cotton genomes i.e., Gossypium herbaceum (A1), Gossypium arboretum (A2), Gossypium raimondii (D5), Gossypium hirsutum (AD1), Gossypium barbadense (AD2), Gossypium tomentosum (AD3), Gossypium mustelinum (AD4), Gossypium darwinii (AD5), respectively. A total of 198 GDSL genes of Gossypium hirsutum were divided into eleven clades using phylogenetic analysis, and the number of GhirGDSL varied among different clades. The cis-elements analysis showed that GhirGDSL gene expression was mainly related to light, plant hormones, and variable tense environments. Combining the results of transcriptome and RT-qPCR, GhirGDSL26 (Gh_A01G1774), a highly up-regulated gene, was selected for further elucidating its tole in drought stress tolerance via estimating physiological and biochemical parameters. Heterologous expression of the GhirGDSL26 gene in Arabidopsis thaliana resulted in a higher germination and survival rates, longer root lengths, lower ion leakage and induced stress-responsive genes expression under drought stress. This further highlighted that overexpressed plants had a better drought tolerance as compared to the wildtype plants. Moreover, 3, 3’-diaminobenzidine (DAB) and Trypan staining results indicated reduced oxidative damage, less cell membrane damage, and lower ion leakage in overexpressed plants as compared to wild type. Silencing of GhirGDSL26 in cotton via VIGS resulting in a susceptible phenotype, higher MDA and H2O2 contents, lower SOD activity, and proline content. Conclusion Our results demonstrated that GhirGDSL26 plays a critical role in cotton drought stress tolerance. Current findings enrich our knowledge of GDSL genes in cotton and provide theoretical guidance and excellent gene resources for improving drought tolerance in cotton.

Salinity is a major threat to the yield and productivity of cotton seedlings. In the present study, we developed a BC2F2 population of cotton plants from Gossypium darwinii (5–7) and Gossypium hirsutum (CCRI 12–4) salt-susceptible parents to identify salt-resistant candidate genes. The Illumina HiSeq™ strategy was used with bulked segregant analysis. Salt-resistant and salt-susceptible DNA bulks were pooled by using 30 plants from a BC2F2 population. Next-generation sequencing (NGS) technology was used for the sequencing of parents and both bulks. Four significant genomic regions were identified: the first genomic region was located on chromosome 18 (1.86 Mb), the second and third genomic regions were on chromosome 25 (1.06 Mb and 1.94 Mb, respectively), and the fourth was on chromosome 8 (1.41 Mb). The reads of bulk1 and bulk2 were aligned to the G. darwinii and G. hirsutum genomes, respectively, leading to the identification of 20,664,007 single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels). After the screening, 6,573 polymorphic markers were obtained after filtration of the candidate regions. The SNP indices in resistant and susceptible bulks and Δ(SNP-index) values of resistant and susceptible bulks were measured. Based on the higher Δ(SNP-index) value, six effective polymorphic SNPs were selected in a different chromosome. Six effective SNPs were linked to five candidate genes in four genomic regions. Further validation of these five candidate genes was carried out using reverse transcription-quantitative polymerase chain reaction (RT-qPCR), resulting in an expression profile that showed two highly upregulated genes in the salt-tolerant species G. darwinii, i.e., Gohir.D05G367800 and Gohir.D12G239100; however, the opposite was shown in G. hirsutum, for which all genes, except one, showed partial expression. The results indicated that Gohir.D05G367800 and Gohir.D12G239100 may be salt-tolerant genes. We are confident that this study could be helpful for the cloning, transformation, and development of salt-resistant cotton varieties.

Abstract Background Crops face several environmental stresses (biotic and abiotic), thus resulting in severe yield losses. Around the globe abiotic stresses are the main contributors of plant damages, primarily drought and salinity. Many genes and transcription factors are involved in abiotic and biotic stress responses. NAC TF (Transcription Factors) improves tolerance to stresses by controlling the physiological and enzyme activities of crops. Results In current research, GhNAC072 a highly upregulated TF in RNA-Seq was identified as a hub gene in the co-expression network analysis (WGCNA). This gene was transformed to Arabidopsis thaliana to confirm its potential role in drought and salt stress tolerance. Significant variations were observed in the morpho-physiological traits with high relative leaf water contents, chlorophyll contents, higher germination and longer root lengths of the overexpressed lines and low excised leaf loss and ion leakage as compared to the wildtype plants. Besides, overexpressed lines have higher amounts of antioxidants and low oxidant enzyme activities than the wildtype during the period of stress exposure. Conclusions In summary, the above analysis showed that GhNAC072 might be the true candidate involved in boosting tolerance mechanisms under drought and salinity stress.

The development of resource-based cities faces with stagnation due to the exhaustion of resources. Scientific and rational characterization of the temporal and spatial characteristics of resource environmental carrying capacity (RECC) is conducive to the sustainable development of these cities. In the present study, 8 resource-based cities in Sichuan Province were used as the research objects. A RECC research framework was constructed based on natural resources, environmental resources, economic resources and social resources dimensions (NEES), and a RECC assessment system was then extracted. This paper constructed a RECC calculation method based on the EWM-CPM (EWM: entropy weight method, CPM: catastrophe progression method). Mean and standard deviation were used for statistical analysis of RECC. The temporal trends of the NEES and RECC of resource-based cities in Sichuan Province were analyzed from the perspective of time, and the dynamic characteristics were analyzed combined with the dynamic evaluation measure model. The spatial pattern of RECC in the resource-based cities in Sichuan Province was analyzed from a spatial perspective. The relationship between NEES and RECC was evaluated through a path analysis method. The research results were as follows: (1) In terms of time, the average value of RECC in the resource-based cities in Sichuan Province was 0.628, and the standard deviation was 0.200, indicating that the stability of RECC in the entire region was poor. Natural resources are the main subsystem of RECC impacts in NEES. From the calculation results of the dynamic evaluation measure model, the level of RECC in PZH (Panzhihua) was stable, and the level of RECC in GA (Guang'an) increased with time, whereas the level of RECC in NC (Nanchong) was poor. (2) In terms of space, the eight resource-based cities showed uneven development in regard to RECC from 2010 to 2019, and the urban advantageous industries were related to the types of resources in the NEES. (3) The path analysis method results indicated that the directional impact of natural resources on RECC was the highest (0.658), followed by the effect of social resources, environmental resources and economic resources. The findings of the present study provide a basis for the establishment of methods for the sustainable development of urban areas, effective information extraction and the integrated assessment of factors that affect the development of resource-based cities.

Abstract Background Cotton is a valuable economic crop and the main significant source of natural fiber for textile industries globally. The effects of drought and salt stress pose a challenge to strong fiber and large-scale production due to the ever-changing climatic conditions. However, plants have evolved a number of survival strategies, among them is the induction of various stress-responsive genes such as the ribosomal protein large (RPL) gene. The RPL gene families encode critical proteins, which alleviate the effects of drought and salt stress in plants. In this study, comprehensive and functional analysis of the cotton RPL genes was carried out under drought and salt stresses. Results Based on the genome-wide evaluation, 26, 8, and 5 proteins containing the RPL14B domain were identified in Gossypium hirsutum, G. raimondii, and G. arboreum, respectively. Furthermore, through bioinformatics analysis, key cis-regulatory elements related to RPL14B genes were discovered. The Myb binding sites (MBS), abscisic acid-responsive element (ABRE), CAAT-box, TATA box, TGACG-motif, and CGTCA-motif responsive to methyl jasmonate, as well as the TCA-motif responsive to salicylic acid, were identified. Expression analysis revealed a key gene, Gh_D01G0234 (RPL14B), with significantly higher induction levels was further evaluated through a reverse genetic approach. The knockdown of Gh_D01G0234 (RPL14B) significantly affected the performance of cotton seedlings under drought/salt stress conditions, as evidenced by a substantial reduction in various morphological and physiological traits. Moreover, the level of the antioxidant enzyme was significantly reduced in VIGS-plants, while oxidant enzyme levels increased significantly, as demonstrated by the higher malondialdehyde concentration level. Conclusion The results revealed the potential role of the RPL14B gene in promoting the induction of antioxidant enzymes, which are key in oxidizing the various oxidants. The key pathways need to be investigated and even as we exploit these genes in the developing of more stress-resilient cotton germplasms.

Lint percentage is one of the most essential yield components and an important economic index for cotton planting. Improving lint percentage is an effective way to achieve high-yield in cotton breeding worldwide, especially upland cotton (Gossypium hirsutum L.). However, the genetic basis controlling lint percentage has not yet been systematically understood. Here, we performed a genome-wide association mapping for lint percentage using a natural population consisting of 189 G. hirsutum accessions (188 accessions of G. hirsutum races and one cultivar TM-1). The results showed that 274 single-nucleotide polymorphisms (SNPs) significantly associated with lint percentage were detected, and they were distributed on 24 chromosomes. Forty-five SNPs were detected at least by two models or at least in two environments, and their 5 Mb up- and downstream regions included 584 makers related to lint percentage identified in previous studies. In total, 11 out of 45 SNPs were detected at least in two environments, and their 550 Kb up- and downstream region contained 335 genes. Through RNA sequencing, gene annotation, qRT-PCR, protein–protein interaction analysis, the cis-elements of the promotor region, and related miRNA prediction, Gh_D12G0934 and Gh_A08G0526 were selected as key candidate genes for fiber initiation and elongation, respectively. These excavated SNPs and candidate genes could supplement marker and gene information for deciphering the genetic basis of lint percentage and facilitate high-yield breeding programs of G. hirsutum ultimately.

Abstract Background Cotton is an important commercial crop for being a valuable source of natural fiber. Its production has undergone a sharp decline because of abiotic stresses, etc. Drought is one of the major abiotic stress causing significant yield losses in cotton. However, plants have evolved self-defense mechanisms to cope abiotic factors like drought, salt, cold, etc. The evolution of stress responsive transcription factors such as the trihelix, a nodule-inception-like protein (NLP), and the late embryogenesis abundant proteins have shown positive response in the resistance improvement to several abiotic stresses. Results Genome wide identification and characterization of the effects of Light-Harvesting Chloro a/b binding (LHC) genes were carried out in cotton under drought stress conditions. A hundred and nine proteins encoded by the LHC genes were found in the cotton genome, with 55, 27, and 27 genes found to be distributed in Gossypium hirsutum, G. arboreum, and G. raimondii, respectively. The proteins encoded by the genes were unevenly distributed on various chromosomes. The Ka/Ks (Non-synonymous substitution rate/Synonymous substitution rate) values were less than one, an indication of negative selection of the gene family. Differential expressions of genes showed that majority of the genes are being highly upregulated in the roots as compared with leaves and stem tissues. Most genes were found to be highly expressed in MR-85, a relative drought tolerant germplasm. Conclusion The results provide proofs of the possible role of the LHC genes in improving drought stress tolerance, and can be explored by cotton breeders in releasing a more drought tolerant cotton varieties.

Trichomes protect plants against insects, microbes, herbivores, and abiotic damages and assist seed dispersal. The function of CPR5 genes have been found to be involved in the trichome development but the research on the underlying genetic and molecular mechanisms are extremely limited. Herein, genome wide identification and characterization of CPR5 genes was performed. In total, 26 CPR5 family members were identified in Gossypium species. Phylogenetic analysis, structural characteristics, and synteny analysis of CPR5s showed the conserved evolution relationships of CPR5. The promoter analysis of CPR5 genes revealed hormone, stress, and development-related cis-elements. Gene ontology (GO) enrichment analysis showed that the CPR5 genes were largely related to biological regulation, developmental process, multicellular organismal process. Protein-protein interaction analysis predicted several trichome development related proteins (SIM, LGO, and GRL) directly interacting with CPR5 genes. Further, nine putative Gossypium-miRNAs were also identified, targeting Gossypium CPR5 genes. RNA-Seq data of G. arboreum (with trichomes) and G. herbaceum (with no trichomes) was used to perform the co-expression network analysis. GheCPR5.1 was identified as a hub gene in a co-expression network analysis. RT-qPCR of GheCPR5.1 gene in different tissues suggests that this gene has higher expressions in the petiole and might be a key candidate involved in the trichome development. Virus induced gene silencing of GheCPR5.1 (Ghe02G17590) confirms its role in trichome development and elongation. Current results provide proofs of the possible role of CPR5 genes and provide preliminary information for further studies of GheCPR5.1 functions in trichome development.

Drought seriously threats the growth and development of Gossypium hirsutum L. To dissect the genetic basis for drought tolerance in the G. hirsutum L. germplasm, a population, consisting of 188 accessions of G. hirsutum races and a cultivar (TM-1), was genotyped using the Cotton80KSNP biochip, and 51,268 high-quality single-nucleotide polymorphisms (SNPs) were obtained. Based on the phenotypic data of eight drought relative traits from four environments, we carried out association mapping with five models using GAPIT software. In total, thirty-six SNPs were detected significantly associated at least in two environments or two models. Among these SNPs, 8 and 28 (including 24 SNPs in 5 peak regions) were distributed in the A and D subgenome, respectively; eight SNPs were found to be distributed within separate genes. An SNP, TM73079, located on chromosome D10, was simultaneously associated with leaf fresh weight, leaf wilted weight, and leaf dry weight. Another nine SNPs, TM47696, TM33865, TM40383, TM10267, TM59672, TM59675, TM59677, TM72359, and TM72361, on chromosomes A13, A10, A12, A5, D6, and D9, were localized within or near previously reported quantitative trait loci for drought tolerance. Moreover, 520 genes located 200 kb up- and down-stream of 36 SNPs were obtained and analyzed based on gene annotation and transcriptome sequencing. The results showed that three candidate genes, Gh_D08G2462, Gh_A03G0043, and Gh_A12G0369, may play important roles in drought tolerance. The current GWAS represents the first investigation into mapping QTL for drought tolerance in G. hirsutum races and provides important information for improving cotton cultivars.

Drought has been identified as a major threat for global crop production worldwide. Phosphofructokinase (PFK) is vital for sugar metabolism. During phosphorylation, plants have two enzymes: ATP-dependent phosphofructokinase (PFK) and pyrophosphate-dependent fructose-6-phosphate phosphotransferase (PFP). Genome-wide identification led to the identification of 80 PFK genes, 26 genes in G. hirsutum and G. barbadense, and 14 genes in G. arboreum and G. raimondii. Phylogenetic, gene structure, and motif analyses showed that PFK genes were grouped into two main categories, namely, PFK and PFP, with 18 and 8 genes in the allotetraploid species and 10 PFK and 4 PFP genes in the diploid species, respectively. Using the RNA-seq expressions of 26 genes from GhPFK, a co-expression network analysis was performed to identify the hub genes. GhPFK04, GhPFK05, GhPFK09, GhPFK11, GhPFK13, GhPFK14, and GhPFK17 in leaves and GhPFK02, GhPFK09, GhPFK11, GhPFK15, GhPFK16, and GhPFK17 in root tissues were found as hub genes. RT-qPCR analysis validated the expressions of identified hub genes. Interestingly, GhPFK11 and GhPFK17 were identified as common hub genes, and these might be the true candidate genes involved in the drought stress tolerance. In the KEGG enrichment analysis, amino acids such as L-valine, L-histidine, L-glutamine, L-serine, L-homoserine, L-methionine, L-cysteine, and gluconic acid were significantly upregulated, whereas sugars, mainly fructose-1-phosphate, D-mannitol, D-sorbitol, dulcitol, and lactose, were significantly downregulated during drought stress. Genome-wide analysis paves the way for a deeper understanding of the PFK genes and establishes the groundwork for future research into PFK’s role in enhancing drought stress tolerance and sugar metabolism in cotton.

The most devastating abiotic factors worldwide are drought and salinity, causing severe bottlenecks in the agricultural sector. To acclimatize to these harsh ecological conditions, plants have developed complex molecular mechanisms involving diverse gene families. Among them, S-adenosyl-L-methionine synthetase (SAMS) genes initiate the physiological, morphological, and molecular changes to enable plants to adapt appropriately. We identified and characterized 16 upland cotton SAMS genes (GhSAMSs). Phylogenetic analysis classified the GhSAMSs into three major groups closely related to their homologs in soybean. Gene expression analysis under drought and salt stress conditions revealed that GhSAMS2, which has shown the highest interaction with GhCBL10 (a key salt responsive gene), was the one that was most induced. GhSAMS2 expression knockdown via virus-induced gene silencing (VGIS) enhanced transgenic plants’ susceptibility to drought and salt stress. The TRV2:GhSAMS2 plants showed defects in terms of growth and physiological performances, including antioxidative processes, chlorophyll synthesis, and membrane permeability. Our findings provide insights into SAMS genes’ structure, classification, and role in abiotic stress response in upland cotton. Moreover, they show the potential of GhSAMS2 for the targeted improvement of cotton plants’ tolerance to multiple abiotic stresses.

Abstract Background Segregation distortion (SD) is a common phenomenon among stable or segregating populations, and the principle behind it still puzzles many researchers. The F2:3 progenies developed from the wild cotton species of the D genomes were used to investigate the possible plant transcription factors within the segregation distortion regions (SDRs). A consensus map was developed between two maps from the four D genomes, map A derived from F2:3 progenies of Gossypium klotzschianum and G. davidsonii while Map B from G. thurberi and G. trilobum F2:3 generations. In each map, 188 individual plants were used. Results The consensus linkage map had 1 492 markers across the 13 linkage groups with a map size of 1 467.445 cM and an average marker distance of 1.037 0 cM. Chromosome D502 had the highest percentage of SD with 58.6%, followed by Chromosome D507 with 47.9%. Six thousand and thirty-eight genes were mined within the SDRs on chromosome D502 and D507 of the consensus map. Within chromosome D502 and D507, 2 308 and 3 730 genes were mined, respectively, and were found to belong to 1 117 gourp out of which 622 groups were common across the two chromosomes. Moreover, genes within the top 9 groups related to plant resistance genes (R genes), whereas 188 genes encoding protein kinase domain (PF00069) comprised the largest group. Further analysis of the dominant gene group revealed that 287 miRNAs were found to target various genes, such as the gra-miR398, gra-miR5207, miR164a, miR164b, miR164c among others, which have been found to target top-ranked stress-responsive transcription factors such as NAC genes. Moreover, some of the stress- responsive cis-regulatory elements were also detected. Furthermore, RNA profiling of the genes from the dominant family showed that higher numbers of genes were highly upregulated under salt and osmotic stress conditions, and also they were highly expressed at different stages of fiber development. Conclusion The results indicated the critical role of the SDRs in the evolution of the key regulatory genes in plants.

Abstract Background Wild species of cotton are excellent resistance to abiotic stress. Diploid D-genome cotton showed abundant phenotypic diversity and was the putative donor species of allotetraploid cotton which produce the largest textile natural fiber. Results A total of 41,053 genes were expressed in all samples by mapping RNA-seq Illumina reads of G. thurberi (D1), G. klotzschianum (D3-k), G. raimondii (D5) and G. trilobum (D8) to reference genome. The numbers of differently expressed genes (DEGs) were significantly higher under cold stress than salt stress. However, 34.1% DEGs under salt stress were overlapped with cold stress in four species. Notably, a potential shared network (cold and salt response, including 16 genes) was mined out by gene co-expression analysis. A total of 47,180–55,548 unique genes were identified in four diploid species by De novo assembly. Furthermore, 163, 344, 330, and 161 positively selected genes (PSGs) were detected in thurberi, G. klotzschianum, G. raimondii and G. trilobum by evolutionary analysis, respectively, and 9.5–17% PSGs of four species were DEGs in corresponding species under cold or salt stress. What’s more, most of PSGs were enriched GO term related to response to stimulation. G. klotzschianum showed the best tolerance under both cold and salt stress. Interestingly, we found that a RALF-like protein coding gene not only is PSGs of G. klotzschianum, but also belongs to the potential shared network. Conclusion Our study provided new evidence that gene expression variations of evolution by natural selection were essential drivers of the morphological variations related to environmental adaptation during evolution. Additionally, there exist shared regulated networks under cold and salt stress, such as Ca2+ signal transduction and oxidation-reduction mechanisms. Our work establishes a transcriptomic selection mechanism for altering gene expression of the four diploid D-genome cotton and provides available gene resource underlying multi-abiotic resistant cotton breeding strategy.

Abstract Background Cotton is mainly grown for its natural fiber and edible oil. The fiber obtained from cotton is the indispensable raw material for the textile industries. The ever changing climatic condition, threatens cotton production due to a lack of sufficient water for its cultivation. Effects of drought stress are estimated to affect more than 50% of the cotton growing regions. To elucidate the drought tolerance phenomenon in cotton, a backcross population was developed from G. tomentosum, a drought tolerant donor parent and G. hirsutum which is highly susceptible to drought stress. Results A genetic map of 10 888 SNP markers was developed from 200 BC2F2 populations. The map spanned 4 191.3 centi-Morgan (cM), with an average distance of 0.104 7 cM, covering 51% and 49% of At and Dt sub genomes, respectively. Thirty stable Quantitative trait loci (QTLs) were detected, in which more than a half were detected in the At subgenome. Eighty-nine candidate genes were mined within the QTL regions for three traits: cell membrane stability (CMS), saturated leaf weight (SLW) and chlorophyll content. The genes had varied physiochemical properties. A majority of the genes were interrupted by introns, and only 15 genes were intronless, accounting for 17% of the mined genes. The genes were found to be involved molecular function (MF), cellular component (CC) and biological process (BP), which are the main gene ontological (GO) functions. A number of miRNAs were detected, such as miR164, which is associated with NAC and MYB genes, with a profound role in enhancing drought tolerance in plants. Through RT-qPCR analysis, 5 genes were found to be the key genes involved in enhancing drought tolerance in cotton. Wild cotton harbors a number of favorable alleles, which can be exploited to aid in improving the narrow genetic base of the elite cotton cultivars. The detection of 30 stable QTLs and 89 candidate genes found to be contributed by the donor parent, G. tomentosum, showed the significant genes harbored by the wild progenitors which can be exploited in developing more robust cotton genotypes with diverse tolerance levels to various environmental stresses. Conclusion This was the first study involving genome wide association mapping for drought tolerance traits in semi wild cotton genotypes. It offers an opportunity for future exploration of these genes in developing highly tolerant cotton cultivars to boost cotton production.

Abstract Background Cotton grows in altering environments that are often unfavorable or stressful for its growth and development. Consequently, the plant must cope with abiotic stresses such as soil salinity, drought, and excessive temperatures. Alkali-salt stress response remains a cumbersome biological process and is regulated via a multifaceted transcriptional regulatory network in cotton. Results To discover the molecular mechanisms of alkali-salt stress response in cotton, a comprehensive transcriptome analysis was carried out after alkali-salt stress treatment in three accessions of Gossypium hirsutum with contrasting phenotype. Expression level analysis proved that alkali-salt stress response presented significant stage-specific and tissue-specific. GO enrichment analysis typically suggested that signal transduction process involved in salt-alkali stress response at SS3 and SS12 stages in leaf; carbohydrate metabolic process and oxidation-reduction process involved in SS48 stages in leaf; the oxidation-reduction process involved at all three phases in the root. The Co-expression analysis suggested a potential GhSOS3/GhCBL10-SOS2 network was involved in salt-alkali stress response. Furthermore, Salt-alkali sensitivity was increased in GhSOS3 and GhCBL10 Virus-induced Gene Silencing (VIGS) plants. Conclusion The findings may facilitate to elucidate the underlying mechanisms of alkali-salt stress response and provide an available resource to scrutinize the role of candidate genes and signaling pathway governing alkali-salt stress response.

Low temperature is a common biological abiotic stress in major cotton-growing areas. Cold stress significantly affects the growth, yield, and yield quality of cotton. Therefore, it is important to develop more robust and cold stress-resilient cotton germplasms. In response to climate change and erratic weather conditions, plants have evolved various survival mechanisms, one of which involves the induction of various stress responsive transcript factors, of which the C-repeat-binding factors (CBFs) have a positive effect in enhancing plants response to cold stress. In this study, genomewide identification and functional characterization of the cotton CBFs were carried out. A total of 29, 28, 25, 21, 30, 26, and 15 proteins encoded by the CBF genes were identified in seven Gossypium species. A phylogenetic evaluation revealed seven clades, with Clades 1 and 6 being the largest. Moreover, the majority of the proteins encoded by the genes were predicted to be located within the nucleus, while some were distributed in other parts of the cell. Based on the transcriptome and RT-qPCR analysis, Gthu17439 (GthCBF4) was highly upregulated and was further validated through forward genetics. The Gthu17439 (GthCBF4) overexpressed plants exhibited significantly higher tolerance to cold stress, as evidenced by the higher germination rate, increased root growth, and high-induction levels of stress-responsive genes. Furthermore, the overexpressed plants under cold stress had significantly reduced oxidative damage due to a reduction in hydrogen peroxide (H2O2) production. Moreover, the overexpressed plants under cold stress had minimal cell damage compared to the wild types, as evidenced by the Trypan and 3,3′-Diaminobenzidine (DAB) staining effect. The results showed that the Gthu17439 (GthCBF4) could be playing a significant role in enhancing cold stress tolerance in cotton and can be further exploited in developing cotton germplasm with improved cold-stress tolerance.

Drought and low-temperature stresses are the most prominent abiotic stresses affecting cotton. Wild cotton being exposed to harsh environments has more potential to cope with both biotic and abiotic stresses. Exploiting wild cotton material to induce resistant germplasm would be of greater interest. The candidate gene was identified in the BC2F2 population among Gossypium tomentosum and Gossypium hirsutum as wild male donor parent noted for its drought tolerance and the recurrent parent and a high yielding but drought susceptible species by genotyping by sequencing (GBS) mapping. Golden2-like (GLK) gene, which belongs to the GARP family, is a kind of plant-specific transcription factor (TF) that was silenced by virus-induced gene silencing (VIGS). Silencing of GhGLK1 in cotton results in more damage to plants under drought and cold stress as compared with wild type (WT). The overexpression of GhGLK1 in Arabidopsis thaliana showed that the overexpressing plants showed more adaptability than the WT after drought and cold treatments. The results of trypan blue and 3,3′-diaminobenzidine (DAB) staining showed that after drought and cold treatment, the leaf damage in GhGLK1 overexpressed plants was less as compared with the WT, and the ion permeability was also lower. This study suggested that the GhGLK1 gene may be involved in the regulation of drought and cold stress response in cotton. Our current research findings add significantly to the existing knowledge of cold and drought stress tolerance in cotton.

Cotton is one of the most important fiber crops globally. Despite this, various abiotic stresses, including drought, cause yield losses. We used transcriptome profiles to investigate the co-expression patterns of gene networks associated with drought stress tolerance. We identified three gene modules containing 3,567 genes highly associated with drought stress tolerance. Within these modules, we identified 13 hub genes based on intramodular significance, for further validation. The yellow module has five hub genes (Gh_A07G0563, Gh_D05G0221, Gh_A05G3716, Gh_D12G1438, and Gh_D05G0697), the brown module contains three hub genes belonging to the aldehyde dehydrogenase (ALDH) gene family (Gh_A06G1257, Gh_A06G1256, and Gh_D06G1578), and the pink module has five hub genes (Gh_A02G1616, Gh_D12G2599, Gh_D07G2232, Gh_A02G0527, and Gh_D07G0629). Based on RT-qPCR results, the Gh_A06G1257 gene has the highest expression under drought stress in different plant tissues and it might be the true candidate gene linked to drought stress tolerance in cotton. Silencing of Gh_A06G1257 in cotton leaves conferred significant sensitivity in response to drought stress treatments. Overexpression of Gh_A06G1257 in Arabidopsis also confirms its role in drought stress tolerance. L-valine, Glutaric acid, L-proline, L-Glutamic acid, and L-Tryptophan were found to be the most significant metabolites playing roles in drought stress tolerance. These findings add significantly to existing knowledge of drought stress tolerance mechanisms in cotton.

The acyl-coenzyme A oxidase 3 (ACX3) gene involved in the β-oxidation pathway plays a critical role in plant growth and development as well as stress response. Earlier on, studies focused primarily on the role of β-oxidation limited to fatty acid breakdown. However, ACX3 peroxisomal β-oxidation pathways result in a downstream cascade of events that act as a transduction of biochemical and physiological responses to stress. A role that is yet to be studied extensively. In this study, we identified 20, 18, 22, 23, 20, 11, and 9 proteins in Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum, G. darwinii, G. arboretum, and G. raimondii genomes, respectively. The tetraploid cotton genome had protein ranging between 18 and 22, while diploids had between 9 and 11. After analyzing the gene family evolution or selection pressure, we found that this gene family undergoes purely segmental duplication both in diploids and tetraploids. W-Box (WRKY-binding site), ABRE, CAAT–Box, TATA-box, MYB, MBS, LTR, TGACG, and CGTCA-motif are abiotic stress cis-regulatory elements identified in this gene family. All these are the binding sites for abiotic stress transcription factors, indicating that this gene is essential. Genes found in G. hirsutum showed a clear response to drought and salinity stress, with higher expression under drought and salt stress, particularly in the leaf and root, according to expression analysis. We selected Gh_DO1GO186, one of the highly expressed genes, for functional characterization. We functionally characterized the GhACX3 gene through overexpression and virus-induced gene silencing (VIGS). Overexpression of this gene enhanced tolerance under stress, which was exhibited by the germination assay. The overexpressed seed growth rate was faster relative to control under drought and salt stress conditions. The survival rate was also higher in overexpressed plants relative to control plants under stress. In contrast, the silencing of the GhACX3 gene in cotton plants resulted in plants showing the stress susceptibility phenotype and reduced root length compared to control. Biochemical analysis also demonstrated that GhACX3-silenced plants experienced oxidative stress while the overexpressed plants did not. This study has revealed the importance of the ACX3 family during stress tolerance and can breed stress-resilient cultivar.

Abstract Background The efficient detection and initiation of appropriate response to abiotic stresses are important to plants survival. The plant G-protein coupled receptors (GPCRs) are diverse membranous proteins that are responsible for signal transduction. Results In this research work, we identified a novel gene of the GPCR domain, transformed and carried out the functional analysis in Arabidopsis under drought and cold stresses. The transgenic lines exposed to drought and cold stress conditions showed higher germination rate, increased root length and higher fresh biomass accumulation. Besides, the levels of antioxidant enzymes, glutathione (GSH) and ascorbate peroxidase (APX) exhibited continuously increasing trends, with approximately threefold higher than the control, implying that these ROS-scavenging enzymes were responsible for the detoxification of ROS induced by drought and cold stresses. Similarly, the transgenic lines exhibited stable cell membrane stability (CMS), reduced water loss rate in the detached leaves and significant values for the saturated leaves compared to the wild types. Highly stress-responsive miRNAs were found to be targeted by the novel gene and based on GO analysis; the protein encoded by the gene was responsible for maintaining an integral component of membrane. In cotton, the virus-induced gene silencing (VIGS) plants exhibited a higher susceptibility to drought and cold stresses compared to the wild types. Conclusion The novel GPCR gene enhanced drought and cold stress tolerance in transgenic Arabidopsis plants by promoting root growth and induction of ROS scavenging enzymes. The outcome showed that the gene had a role in enhancing drought and cold stress tolerance, and can be further exploited in breeding for more stress-resilient and tolerant crops.

Abstract Background Auxins play an important role in plant growth and development; the auxins responsive gene; auxin/indole-3-acetic acid (Aux/IAA), small auxin-up RNAs (SAUR) and Gretchen Hagen3 (GH3) control their mechanisms. The GH3 genes function in homeostasis by the catalytic activities in auxin conjugation and bounding free indole-3-acetic acid (IAA) to amino acids. Results In our study, we identified the GH3 genes in three cotton species; Gossypium hirsutum, Gossypium arboreum and Gossypium raimondii, analyzed their chromosomal distribution, phylogenetic relationships, cis-regulatory element function and performed virus induced gene silencing of the novel Gh_A08G1120 (GH3.5) gene. The phylogenetic tree showed four clusters of genes with clade 1, 3 and 4 having mainly members of the GH3 of the cotton species while clade 2 was mainly members belonging to Arabidopsis. There were no paralogous genes, and few orthologous genes were observed between Gossypium and other species. All the GO terms were detected, but only 14 genes were found to have described GO terms in upland cotton, more biological functions were detected, as compared to the other functions. The GH3.17 subfamily harbored the highest number of the cis-regulatory elements, most having promoters towards dehydration-responsiveness. The RNA expression analysis revealed that 10 and 8 genes in drought and salinity stress conditions respectively were upregulated in G. hirsutum. All the genes that were upregulated in plants under salt stress conditions were also upregulated in drought stress; moreover, Gh_A08G1120 (GH3.5) exhibited a significant upregulation across the two stress factors. Functional characterization of Gh_A08G1120 (GH3.5) through virus-induced gene silencing (VIGS) revealed that the VIGS plants ability to tolerate drought and salt stresses was significantly reduced compared to the wild types. The chlorophyll content, relative leaf water content (RLWC), and superoxide dismutase (SOD) concentration level were reduced significantly while malondialdehyde concentration and ion leakage as a measure of cell membrane stability (CMS) increased in VIGS plants under drought and salt stress conditions. Conclusion This study revealed the significance of the GH3 genes in enabling the plant’s adaptation to drought and salt stress conditions as evidenced by the VIGS results and RT-qPCR analysis.

Abstract Background Cotton (Gossypium hirsutum) provides the largest natural fiber for the textile manufacturing industries, but its production is on the decline due to the effects of salinity. Soil salt-alkalization leads to damage in cotton growth and a decrease in yields. Hyperosmolality-gated calcium-permeable channels (OSCA) have been found to be involved in the detection of extracellular changes which trigger an increase in cytosolic free calcium concentration. Hyperosmolality-induced calcium ion increases have been widely speculated to be playing a role in osmosensing in plants. However, the molecular nature of the corresponding calcium ion channels remains unclearly. In this research work, we describe the OSCA genes and their putative function in osmosensing in plants by carrying out genome-wide identification, characterization and functional analysis of the significantly up-regulated OSCA gene, GhOSCA1.1 through reverse genetics. Result A total of 35, 21 and 22 OSCA genes were identified in G. hirsutum, G. arboreum, and G. raimondii genomes, respectively, and were classified into four different clades according to their gene structure and phylogenetic relationship. Gene and protein structure analysis indicated that 35 GhOSCA genes contained a conserved RSN1_7TM (PF02714) domain. Moreover, the cis-regulatory element analysis indicated that the OSCA genes were involved in response to abiotic stress. Furthermore, the knockdown of one of the highly up-regulated genes, Gh_OSCA1.1 showed that the virus-induced gene silenced (VIGS) plants were highly sensitive to dehydration and salinity stresses compared with the none VIGS plants as evident with higher concentration levels of oxidant enzymes compared with the antioxidant enzymes on the leaves of the stressed plants. Conclusion This study provides the first systematic analysis of the OSCA gene family and will be important for understanding the putative functions of the proteins encoded by the OSCA genes in cotton. These results provide a new insight of defense responses in general and lay the foundation for further investigation of the molecular role played by the OSCA genes, thereby providing suitable approaches to improve crop performance under salinity and drought stress conditions.

Abstract Background Fluorescence in situ hybridization (FISH) is an efficient cytogenetic technology to study chromosome structure. Transposable element (TE) is an important component in eukaryotic genomes and can provide insights in the structure and evolution of eukaryotic genomes. Results A FISH probe derived from bacterial artificial chromosome (BAC) clone 299N22 generated striking signals on all 26 chromosomes of the cotton diploid A genome (AA, 2x=26) but very few on the diploid D genome (DD, 2x=26). All 26 chromosomes of the A sub genome (At) of tetraploid cotton (AADD, 2n=4x=52) also gave positive signals with this FISH probe, whereas very few signals were observed on the D sub genome (Dt). Sequencing and annotation of BAC clone 299N22, revealed a novel Ty3/gypsy transposon family, which was named as ‘CICR’. This family is a significant contributor to size expansion in the A (sub) genome but not in the D (sub) genome. Further FISH analysis with the LTR of CICR as a probe revealed that CICR is lineage-specific, since massive repeats were found in A and B genomic groups, but not in C–G genomic groups within the Gossypium genus. Molecular evolutionary analysis of CICR suggested that tetraploid cottons evolved after silence of the transposon family 1–1.5 million years ago (Mya). Furthermore, A genomes are more homologous with B genomes, and the C, E, F, and G genomes likely diverged from a common ancestor prior to 3.5–4 Mya, the time when CICR appeared. The genomic variation caused by the insertion of CICR in the A (sub) genome may have played an important role in the speciation of organisms with A genomes. Conclusions The CICR family is highly repetitive in A and B genomes of Gossypium, but not amplified in the C–G genomes. The differential amount of CICR family in At and Dt will aid in partitioning sub genome sequences for chromosome assemblies during tetraploid genome sequencing and will act as a method for assessing the accuracy of tetraploid genomes by looking at the proportion of CICR elements in resulting pseudochromosome sequences. The timeline of the expansion of CICR family provides a new reference for cotton evolutionary analysis, while the impact on gene function caused by the insertion of CICR elements will be a target for further analysis of investigating phenotypic differences between A genome and D genome species.

Segregation distortion (SD) is a genetic mechanism commonly found in segregating or stable populations. The principle behind this puzzles many researchers. The F2 generation developed from wild Gossypium darwinii and G. hirsutum CCRI12 species was used to investigate the possible transcription factors within the segregation distortion regions (SDRs). The 384 out of 2763 markers were distorted in 29 SDRs on 18 chromosomes. Good collinearity was observed among genetic and physical maps of G. hirsutum and G. barbadense syntenic blocks. Total 568 genes were identified from SDRs of 18 chromosomes. Out of these genes, 128 belonged to three top-ranked salt-tolerant gene families. The DUF597 contained 8 uncharacterized genes linked to Pkinase (PF00069) gene family in the phylogenetic tree, while 15 uncharacterized genes clustered with the zinc finger gene family. Two hundred thirty four miRNAs targeted numerous genes, including ghr-miR156, ghr-miR399 and ghr-miR482, while others targeted top-ranked stress-responsive transcription factors. Moreover, these genes were involved in the regulation of numerous stress-responsive cis-regulatory elements. The RNA sequence data of fifteen upregulated genes were verified through the RT-qPCR. The expression profiles of two highly upregulated genes (Gh_D01G2015 and Gh_A01G1773) in salt-tolerant G. darwinii showed antagonistic expression in G. hirsutum. The results indicated that salt-tolerant genes have been possibly transferred from the wild G. darwinii species. A detailed functional analysis of these genes can be carried out which might be helpful in the future for gene cloning, transformation, gene editing and the development of salt-resistant cotton varieties.

As one of the pioneer crops widely planted in saline-alkaline areas, Gossypium provides daily necessities, including natural fiber, vegetable proteins, and edible oils. However, cotton fiber yield and quality are highly influenced by salt stress. Therefore, elucidating the molecular mechanisms of cotton in response to salinity stress is importance to breed new cultivars with high tolerance. In this study, we first developed a method for single-cell RNA-seq based on isolating protoplast from cotton root tips; then, we studied the impact of salinity stress on gene expression profiling and their dynamic changes using the developed high-efficiency method for protoplast dissociation suitable for single-cell RNA-seq. A total of 3391 and 2826 differentially expressed genes (DEGs) were identified in salt-treated samples before and after protoplast dissociation, respectively, which were enriched into several molecular components, including response to stimulus, response to stress, and cellular macromolecule metabolic process by gene ontology (GO) analysis. Plant hormone signal transduction, phenylpropanoid biosynthesis, and MAPK signaling pathway were found to be enriched via Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Twenty-two and nine salinity-responsive DEGs participated in plant hormone signaling and MAPK signaling in roots, before and after protoplast dissociation, respectively; six upregulated DEGs were involved in ABA signaling transduction, namely, Ga04G2111, Ga07G0142, Ga09G2061, Ga10G0262, Ga01G0063, and Ga08G1915 which indicates their potential functions on plants adapting to salt stress. Additionally, 384 and 257 transcription factors (TFs) were differentially expressed in salt-stress roots before and after protoplast dissociation, respectively, of which significantly up-regulated TFs mainly belonged to the AP2/ERF-ERF family, which implied their potential roles responding to salt stress. These results not only provide novel insights to reveal the regulatory networks in plant’s root response to salt stress, but also lay the solid foundation for further exploration on cellular heterogeneity by single-cell transcriptome sequencing.

Background: In recent studies, microRNAs have been demonstrated as stable detectable biomarkers in blood for cancer. In addition, computer-aided biomarker discovery has now become an attractive paradigm for precision diagnosis. Methods: In this study, we identified and evaluated miR-139-3p as a biomarker for screening of esophageal squamous cell carcinoma using the Cancer Genome Atlas and Gene Expression Omnibus database analyses. We identified possible miR-139-3p target genes through the predicted database and esophageal squamous cell carcinoma upregulated genes from the Cancer Genome Atlas and Gene. Bioinformatics analysis was performed to determine key miR-139-3p targets and pathways associated with esophageal carcinoma. Finally, the expression and expected significance of hub genes were evaluated via the Genotype-Tissue Expression project. Results: MiR-139-3p was significantly downregulated in patients with esophageal squamous cell carcinoma/esophageal carcinoma. In GSE 122497, the area under the curve-receiver operating characteristic value, sensitivity, and specificity for serum miR-139-3p were 0.754, 67.49%, and 80.00%, respectively. The pattern specification process, skeletal system development, and regionalization process were the most enriched interactions in esophageal carcinoma. In addition, Epstein-Barr virus infection, human T-cell leukemia virus 1 infection, and human cytomegalovirus infection were identified as crucial pathways. Six hub genes (CD1A, FCGR2A, ANPEP, CD1B, membrane metalloendopeptidase, and TWIST1) were found, and FCGR2A and membrane metalloendopeptidase were further confirmed by genotype-tissue expression. High expression of membrane metalloendopeptidase correlated with a better overall survival but not with disease-free survival of patients with esophageal carcinoma. Conclusions: MiR-139-3p was identified as a candidate biomarker for predicting esophageal squamous cell carcinoma based on network analysis. MiR-139-3p acted as a tumor suppressor by targeting membrane metalloendopeptidase in esophageal carcinoma, and low expression of membrane metalloendopeptidase may indicate a better prognosis of patients with esophageal carcinoma.

The activity of genome-specific repetitive sequences is the main cause of genome variation between Gossypium A and D genomes. Through comparative analysis of the two genomes, we retrieved a repetitive element termed ICRd motif, which appears frequently in the diploid Gossypium raimondii (D5) genome but rarely in the diploid Gossypium arboreum (A2) genome. We further explored the existence of the ICRd motif in chromosomes of G. raimondii, G. arboreum, and two tetraploid (AADD) cotton species, Gossypium hirsutum and Gossypium barbadense, by fluorescence in situ hybridization (FISH), and observed that the ICRd motif exists in the D5 and D-subgenomes but not in the A2 and A-subgenomes. The ICRd motif comprises two components, a variable tandem repeat (TR) region and a conservative sequence (CS). The two constituents each have hundreds of repeats that evenly distribute across 13 chromosomes of the D5genome. The ICRd motif (and its repeats) was revealed as the common conservative region harbored by ancient Long Terminal Repeat Retrotransposons. Identification and investigation of the ICRd motif promotes the study of A and D genome differences, facilitates research on Gossypium genome evolution, and provides assistance to subgenome identification and genome assembling.

The extrusion of toxins and substances at a cellular level is a vital life process in plants under abiotic stress. The multidrug and toxic compound extrusion (MATE) gene family plays a large role in the exportation of toxins and other substrates. We carried out a genome-wide analysis of MATE gene families in Gossypium raimondii and Gossypium arboreum and assessed their expression levels under salt, cadmium and drought stresses. We identified 70 and 68 MATE genes in G. raimondii and G. arboreum, respectively. The majority of the genes were predicted to be localized within the plasma membrane, with some distributed in other cell parts. Based on phylogenetic analysis, the genes were subdivided into three subfamilies, designated as M1, M2 and M3. Closely related members shared similar gene structures, and thus were highly conserved in nature and have mainly evolved through purifying selection. The genes were distributed in all chromosomes. Twenty-nine gene duplication events were detected, with segmental being the dominant type. GO annotation revealed a link to salt, drought and cadmium stresses. The genes exhibited differential expression, with GrMATE18, GrMATE34, GaMATE41 and GaMATE51 significantly upregulated under drought, salt and cadmium stress, and these could possibly be the candidate genes. Our results provide the first data on the genome-wide and functional characterization of MATE genes in diploid cotton, and are important for breeders of more stress-tolerant cotton genotypes.

We found 33, 17, and 20 Alba genes in Gossypium hirsutum, Gossypium arboretum, and Gossypium raimondii, respectively. The Alba protein lengths ranged from 62 to 312 aa, the molecular weight (MW) from 7.003 to 34.55 kDa, grand average hydropathy values of −1.012 to 0.609 and isoelectric (pI) values of −3 to 11. Moreover, miRNAs such as gra-miR8770 targeted four genes, gra-miR8752 and gra-miR8666 targeted three genes, and each and gra-miR8657 a, b, c, d, e targeted 10 genes each, while the rests targeted 1 to 2 genes each. Similarly, various cis-regulatory elements were detected with significant roles in enhancing abiotic stress tolerance, such as CBFHV (RYCGAC) with a role in cold stress acclimation among others. Two genes, Gh_D01G0884 and Gh_D01G0922, were found to be highly induced under water deficit and salt stress conditions. Through virus-induced gene silencing (VIGS), the VIGS cotton plants were found to be highly susceptible to both water deficit and salt stresses; the VIGS plants exhibited a significant reduction in root growth, low cell membrane stability (CMS), saturated leaf weight (SLW), chlorophyll content levels, and higher excised leaf water loss (ELWL). Furthermore, the stress-responsive genes and ROS scavenging enzymes were significantly reduced in the VIGS plants compared to either the wild type (WT) and or the positively controlled plants. The VIGS plants registered higher concentration levels of hydrogen peroxide and malondialdehyde, with significantly lower levels of the various antioxidants evaluated an indication that the VIGS plants were highly affected by salt and drought stresses. This result provides a key foundation for future exploration of the Alba proteins in relation to abiotic stress.

Abstract Background Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally. Results In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars. Conclusion The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

Cytogenetic maps of Gossypium hirsutum (Linnaeus, 1753) homoeologous chromosomes Ah01 and Dh01 were constructed by fluorescence in situ hybridization (FISH), using eleven homoeologous-chromosomes-shared bacterial artificial chromosomes (BACs) clones and one chromosome-specific BAC clone respectively. We compared the cytogenetic maps with the genetic linkage and draft genome assembly maps based on a standardized map unit, relative map position (RMP), which allowed a global view of the relationship of genetic and physical distances along each chromosome, and assembly quality of the draft genome assembly map. By integration of cytogenetic maps with sequence maps of the two chromosomes (Ah01 and Dh01), we inferred the locations of two scaffolds and speculated that some homologous sequences belonging to homoeologous chromosomes were removed as repetitiveness during the sequence assembly. The result offers molecular tools for cotton genomics research and also provides valuable information for the improvement of the draft genome assembly.

Abstract Background Chromosome microdissection is one of the most important techniques in molecular cytogenetic research. Cotton (Gossypium Linnaeus, 1753) is the main natural fiber crop in the world. The resistance gene analog (RGA) cloning after its single chromosome microdissection can greatly promote cotton genome research and breeding. Results Using the linker adaptor PCR (LA-PCR) with the primers of rice disease-resistance homologues, three nucleotide sequences PS016 (KU051681), PS054 (KU051682), and PS157 (KU051680) were obtained from the chromosome Ah01 of upland cotton (cv. TM-1). The Blast results showed that the three sequences are the nucleotide binding site-leucine rich repeat (NBS-LRR) type RGAs. Clustering results indicated that they are homologous to these published RGAs. Thus, the three RGAs can definitely be confirmed as NBS-LRR class of RGAs in upland cotton. Conclusions Using single chromosome microdissection technique, DNA libraries containing cotton RGAs were obtained. This technique can promote cotton gene cloning, marker development and even the improvement of cotton genome research and breeding.

Abiotic stresses have negative effects on plants growth and development. Plants, being sessile, have developed specific adaptive strategies that allow them to rapidly detect and respond to abiotic stress factors. The detoxification efflux carriers (DTX)/multidrug and toxic compound extrusion (MATE) transporters are of significance in the translocation of abscisic acid (ABA), a phytohormone with profound role in plants under various abiotic stress conditions. The ABA signaling cascades are the core regulators of abiotic stress responses in plants, triggering major changes in gene expression and adaptive physiological responses. We therefore carried out genome-wide analysis of the DTX/MATE gene family, transformed a DTX/MATE gene in Arabidopsis and carried out functional analysis under drought, salt, and cold stress conditions. We identified 128, 70, and 72 DTX/MATE genes in Gossypium hirsutum, Gossypium arboreum, and Gossypium raimondii, respectively. The proteins encoded by the DTX/MATE genes showed varied physiochemical properties but they all were hydrophobic. The Gh_D06G0281 (DTX/MATE) over-expressing Arabidopsis lines were highly tolerant under drought, salt, and cold stress with high production of antioxidant enzymes and significantly reduced levels of oxidants. Lipid peroxidation, as measured by the level of malondialdehyde concentrations was relatively low in transgenic lines compared to wild types, an indication of reduced oxidative stress levels in the transgenic plants. Based on physiological measurements, the transgenic plants exhibited significantly higher relative leaf water content, reduced excised leaf water loss and a significant reduction in ion leakage as a measure of the cell membrane stability compared to the wild types. Abiotic stress responsive genes, ABF4, CBL1, SOS1, and RD29B were highly expressed in the transgenic lines compared to the non-transformed wild type plants. The protein encoded by the Gh_D06G0281 (DTX/MATE) gene was predicted to be located within the plasma membrane. Since signals from extracellular stimuli are transmitted through the plasma membrane most of which are conducted by plasma membrane proteins it is possible the Gh_D06G0281 (DTX/MATE) gene product could be important for this process.

The high recombination rate of the electron-hole pair on the surface of rutile TiO2 (RT) reduces its photocatalytic performance, although it has high thermodynamic stability and few internal grain defects. Therefore, it is necessary for RT to develop effective methods to reduce electron-hole pair recombination. In this study, magnetic α-Fe2O3/Rutile TiO2 self-assembled hollow spheres were fabricated via a facile hydrothermal reaction and template-free method. Based on the experimental result, phosphate concentration was found to play a crucial role in controlling the shape of these hollow α-Fe2O3/RT nanospheres, and the optimal concentration is 0.025 mM. Due to a heterojunction between α-Fe2O3 and RT, the electron-hole pair recombination rate was reduced, the as-synthesized hollow α-Fe2O3/RT nanospheres exhibited excellent photocatalysis in rhodamine B (RhB) photodegradation compared to α-Fe2O3 and RT under visible-light irradiation, and the degradation rate was about 16% (RT), 60% (α-Fe2O3), and 93% (α-Fe2O3/RT) after 100 min. Moreover, α-Fe2O3/RT showed paramagnetism and can be recycled to avoid secondary environmental pollution.

Nitrogen is a key macronutrient needed by plants to boost their production, but the development of cotton genotypes through conventional approaches has hit a bottleneck due to the narrow genetic base of the elite cotton cultivars, due to intensive selection and inbreeding. Based on our previous research, in which the BC2F2 generations developed from two upland cotton genotypes, an abiotic stress-tolerant genotype, G. tomentosum (donor parent) and a highly-susceptible, and a highly-susceptible, but very productive, G. hirsutum (recurrent parent), were profiled under drought stress conditions. The phenotypic and the genotypic data generated through genotyping by sequencing (GBS) were integrated to map drought-tolerant quantitative trait loci (QTLs). Within the stable QTLs region for the various drought tolerance traits, a nodule-inception-like protein (NLP) gene was identified. We performed a phylogenetic analysis of the NLP proteins, mapped their chromosomal positions, intron-exon structures and conducted ds/dn analysis, which showed that most NLP genes underwent negative or purifying selection. Moreover, the functions of one of the highly upregulated genes, Gh_A05G3286 (Gh NLP5), were evaluated using the virus gene silencing (VIGS) mechanism. A total of 226 proteins encoded by the NLP genes were identified, with 105, 61, and 60 in Gossypium hirsutum, G. raimondii, and G. arboreum, respectively. Comprehensive Insilico analysis revealed that the proteins encoded by the NLP genes had varying molecular weights, protein lengths, isoelectric points (pI), and grand hydropathy values (GRAVY). The GRAVY values ranged from a negative one to zero, showing that proteins were hydrophilic. Moreover, various cis-regulatory elements that are the binding sites for stress-associated transcription factors were found in the promoters of various NLP genes. In addition, many miRNAs were predicted to target NLP genes, notably miR167a, miR167b, miR160, and miR167 that were previously shown to target five NAC genes, including NAC1 and CUC1, under N-limited conditions. The real-time quantitative polymerase chain reaction (RT-qPCR) analysis, revealed that five genes, Gh_D02G2018, Gh_A12G0439, Gh_A03G0493, Gh_A03G1178, and Gh_A05G3286 were significantly upregulated and perhaps could be the key NLP genes regulating plant response under N-limited conditions. Furthermore, the knockdown of the Gh_A05G3286 (GhNLP5) gene by virus-induced silencing (VIGS) significantly reduced the ability of these plants to the knockdown of the Gh_A05G3286 (GhNLP5) gene by virus-induced gene silencing (VIGS) significantly reduced the ability of the VIGS-plants to tolerate N-limited conditions compared to the wild types (WT). The VIGS-plants registered lower chlorophyll content, fresh shoot biomass, and fresh root biomass, addition to higher levels of malondialdehyde (MDA) and significantly reduced levels of proline, and superoxide dismutase (SOD) compared to the WT under N-limited conditions. Subsequently, the expression levels of the Nitrogen-stress responsive genes, GhTap46, GhRPL18A, and GhKLU were shown to be significantly downregulated in VIGS-plants compared to their WT under N-limited conditions. The downregulation of the nitrogen-stress responsive genes provided evidence that the silenced gene had an integral role in enhancing cotton plant tolerance to N-limited conditions.

The problems of mode mixing, mode splitting, and pseudocomponents caused by intermittence or white noise signals during empirical mode decomposition (EMD) are difficult to resolve. The partly ensemble EMD (PEEMD) method is introduced first. The PEEMD method can eliminate mode mixing via the permutation entropy (PE) of the intrinsic mode functions (IMFs). Then, bilateral permutation entropy (BPE) of the IMFs is proposed as a means to detect and eliminate mode splitting by means of the reconstructed signals in the PEEMD. Moreover, known ingredient component signals are comparatively designed to verify that the PEEMD method can effectively detect and progressively address the problem of mode splitting to some degree and generate IMFs with better performance. The microseismic signal is applied to prove, by means of spectral analysis, that this method is effective.

Memory in finance is the foundation of a well-established forecasting model, and new financial theory research shows that the stochastic memory model depends on different time windows. To accurately identify the multivariate long memory model in the financial market, this paper proposes the concept of a moving V-statistic on the basis of a modified R/S method to determine whether the time series has a long-range dependence and subsequently to apply wavelet-based multiresolution analysis to study the multifractality of the financial time series to determine the initial data windows. Finally, we check the moving V-statistic estimation in wavelet analysis in the same condition; the paper selects the volatilities of the gold foreign exchange rates to evaluate the moving V-statistic. According to the results, the method of testing memory established in this paper can identify the breakpoint of the memories effectively. Furthermore, this method can provide support for forecasting returns in the financial market.

Wild cotton species have significant agronomic traits that can be introgressed into elite cultivated varieties. The use of a genetic map is important in exploring, identification and mining genes which carry significant traits. In this study, 188 F2mapping individuals were developed from Gossypium thurberi (female) and Gossypium trilobum (male), and were genotyped by using simple sequence repeat (SSR) markers. A total of 12,560 simple sequence repeat (SSR) markers, developed by Southwest University, thus coded SWU were screened out of which only 994 were found to be polymorphic, and 849 markers were linked in all the 13 chromosomes. The map had a length of 1,012.458 cM with an average marker distance of 1.193 cM. Segregation distortion regions (SDRs) were observed on Chr01, Chr02, Chr06, Chr07 Chr09, Chr10 and Chr11 with a large proportion of the SDR regions segregating towards the heterozygous allele. There was good syntenic block formation that revealed good collinearity between the genetic and physical map of G. raimondii, compared to the Dt_sub genome of the G. hirsutum and G. barbadense. A total of 2,496 genes were mined within the SSR related regions. The proteins encoding the mined genes within the SDR had varied physiochemical properties; their molecular weights ranged from 6.586 to 252.737 kDa, charge range of -39.5 to 52, grand hydropathy value (GRAVY) of -1.177 to 0.936 and isoelectric (pI) value of 4.087 to 12.206. The low GRAVY values detected showed that the proteins encoding these genes were hydrophilic in nature, a property common among the stress responsive genes. The RNA sequence analysis revealed more of the genes were highly upregulated in various stages of fiber development for instance; Gorai.002G241300 was highly up regulated at 5, 10, 20 and 25 day post anthesis (DPA). Validation through RT-qPCR further revealed that these genes mined within the SDR regions might be playing a significant role under fiber development stages, therefore we infer that Gorai.007G347600 (TFCA), Gorai.012G141600 (FOLB1), Gorai.006G024500 (NMD3), Gorai.002G229900 (LST8) and Gorai.002G235200 (NSA2) are significantly important in fiber development and in turn the quality, and further researches needed to be done to elucidate their exact roles in the fiber development process. The construction of the genetic map between the two wild species paves away for the mapping of quantitative trait loci (QTLs) since the average distance between the markers is small, and mining of genes on the SSR regions will provide an insight in identifying key genes that can be introgressed into the cultivated cotton cultivars.

Cotton is a pioneer of saline land crop, while salt stress still causes its growth inhibition and fiber production decrease. Phenotype identification showed better salt tolerance of a wild diploid cotton species Gossypium klotzschianum. To elucidate the salt-tolerant mechanisms in G. klotzschianum, we firstly detected the changes in hormones, H2O2 and glutathione (GSSH and GSH), then investigated the gene expression pattern of roots and leaves treated with 300 mM NaCl for 0, 3, 12, 48 h, and each time control by RNA-seq on the Illumina-Solexa platform. Physiological determination proved that the significant increase in hormone ABA at 48 h, while that in H2O2 was at 12 h, likewise, the GSH content decrease at 48 h and the GSSH content increase at 48 h, under salt stress. In total, 37,278 unigenes were identified from the transcriptome data, 8,312 and 6,732 differentially expressed genes (DEGs) were discovered to be involved in salt stress tolerance in roots and leaves, respectively. Gene function annotation and expression analysis elucidated hormone biosynthesis and signal transduction, reactive oxygen species (ROS), and salt overly sensitive (SOS) signal transduction related genes revealed the important roles of them in signal transmission, oxidation balance and ion homeostasis in response to salinity stress. This is a report which focuses on primary response to highly salty stress (upto 300 mM NaCl) in cotton using a wild diploid Gossypium species, broadening our understanding of the salt tolerance mechanism in cotton and laying a solid foundation of salt resistant for the genetic improvement of upland cotton with the resistance to salt stress.

Abiotic stress is an important limiting factor in crop growth and yield around the world. Owing to the continued genetic erosion of the upland cotton germplasm due to intense selection and inbreeding, attention has shifted towards wild cotton progenitors which offer unique traits that can be introgressed into the cultivated cotton to improve their genetic performance. The purpose of this study was to characterize the Pkinase gene family in a previously developed genetic map of the F2 population derived from a cross between two cotton species: Gossypium hirsutum (CCRI 12-4) and Gossypium darwinii (5-7). Based on phylogenetic analysis, Pkinase (PF00069) was found to be the dominant domain with 151 genes in three cotton species, categorized into 13 subfamilies. Structure analysis of G. hirsutum genes showed that a greater percentage of genes and their exons were highly conserved within the group. Syntenic analysis of gene blocks revealed 99 duplicated genes among G. hirsutum, Gossypium arboreum and Gossypium raimondii. Most of the genes were duplicated in segmental pattern. Expression pattern analysis showed that the Pkinase gene family possessed species-level variation in induction to salinity and G. darwinii had higher expression levels as compared to G. hirsutum. Based on RNA sequence analysis and preliminary RT-qPCR verification, we hypothesized that the Pkinase gene family, regulated by transcription factors (TFs) and miRNAs, might play key roles in salt stress tolerance. These findings inferred comprehensive information on possible structure and function of Pkinase gene family in cotton under salt stress.

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.

We identified 672, 374, and 379 CYPs proteins encoded by the CYPs genes in Gossypium hirsutum, Gossypium raimondii, and Gossypium arboreum, respectively. The genes were found to be distributed in all 26 chromosomes of the tetraploid cotton, with chrA05, chrA12, and their homeolog chromosomes harboring the highest number of genes. The physiochemical properties of the proteins encoded by the CYP450 genes varied in terms of their protein lengths, molecular weight, isoelectric points (pI), and even grand hydropathy values (GRAVY). However, over 99% of the cotton proteins had GRAVY values below 0, which indicated that the majority of the proteins encoded by the CYP450 genes were hydrophilic in nature, a common property of proteins encoded by stress-responsive genes. Moreover, through the RNA interference (RNAi) technique, the expression levels of Gh_D07G1197 and Gh_A13G2057 were suppressed, and the silenced plants showed a higher concentration of hydrogen peroxide (H2O2) with a significant reduction in the concentration levels of glutathione (GSH), ascorbate peroxidase (APX), and proline compared to the wild types under drought and salt stress conditions. Furthermore, the stress-responsive genes 1-Pyrroline–5-Carboxylate Synthetase (GhP5CS), superoxide dismutase (GhSOD), and myeloblastosis (GhMYB) were downregulated in VIGS plants, but showed upregulation in the leaf tissues of the wild types under drought and salt stress conditions. In addition, CYP450-silenced cotton plants exhibited a high level of oxidative injury due to high levels of oxidant enzymes, in addition to negative effects on CMS, ELWL, RLWC, and chlorophyll content The results provide the basic foundation for future exploration of the proteins encoded by the CYP450 genes in order to understand the physiological and biochemical mechanisms in enhancing drought and salt stress tolerance in plants.