Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium.
ABSTRACT: BACKGROUND:Cotton (Gossypium spp.) is the most important world-wide fiber crop but salt stress limits cotton production in coastal and other areas. Growth regulation factors (GRFs) play regulatory roles in response to salt stress, but their roles have not been studied in cotton under salt stress. RESULTS:We identified 19 GRF genes in G. raimondii, 18 in G. arboreum, 34 in G. hirsutum and 45 in G. barbadense, respectively. These GRF genes were phylogenetically analyzed leading to the recognition of seven GRF clades. GRF genes from diploid cottons (G. raimondii and G. arboreum) were largely retained in allopolyploid cotton, with subsequent gene expansion in G. barbadense relative to G. hirsutum. Most G. hirsutum GRF (GhGRF) genes are preferentially expressed in young and growing tissues. To explore their possible role in salt stress, we used qRT-PCR to study expression responses to NaCl treatment, showing that five GhGRF genes were down-regulated in leaves. RNA-seq experiments showed that seven GhGRF genes exhibited decreased expression in leaves under NaCl treatment, three of which (GhGRF3, GhGRF4, and GhGRF16) were identified by both RNA-seq and qRT-PCR. We also identified six and three GRF genes that exhibit decreased expression under salt stress in G. arboreum and G. barbadense, respectively. Consistent with its lack of leaf withering or yellowing under the salt treatment conditions, G. arboreum had better salt tolerance than G. hirsutum and G. barbadense. Our results suggest that GRF genes are involved in salt stress responses in Gossypium. CONCLUSION:In summary, we identified candidate GRF genes that were involved in salt stress responses in cotton.
Project description:Nucleotide binding site (NBS) genes encode a large family of disease resistance (R) proteins in plants. The availability of genomic data of the two diploid cotton species, Gossypium arboreum and Gossypium raimondii, and the two allotetraploid cotton species, Gossypium hirsutum (TM-1) and Gossypium barbadense allow for a more comprehensive and systematic comparative study of NBS-encoding genes to elucidate the mechanisms of cotton disease resistance.Based on the genome assembly data, 246, 365, 588 and 682 NBS-encoding genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The distribution of NBS-encoding genes among the chromosomes was nonrandom and uneven, and was tended to form clusters. Gene structure analysis showed that G. arboreum and G. hirsutum possessed a greater proportion of CN, CNL, and N genes and a lower proportion of NL, TN and TNL genes compared to that of G. raimondii and G. barbadense, while the percentages of RN and RNL genes remained relatively unchanged. The percentage changes among them were largest for TNL genes, about 7 times. Exon statistics showed that the average exon numbers per NBS gene in G. raimondii and G. barbadense were all greater than that in G. arboretum and G. hirsutum. Phylogenetic analysis revealed that the TIR-NBS genes of G. barbadense were closely related with that of G. raimondii. Sequence similarity analysis showed that diploid cotton G. arboreum possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. hirsutum, while diploid G. raimondii possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. barbadense. The synteny analysis showed that more NBS genes in G. raimondii and G. arboreum were syntenic with that in G. barbadense and G. hirsutum, respectively.The structural architectures, amino acid sequence similarities and synteny of NBS-encoding genes between G. arboreum and G. hirsutum, and between G. raimondii and G. barbadense were the highest among comparisons between the diploid and allotetraploid genomes, indicating that G. hirsutum inherited more NBS-encoding genes from G. arboreum, while G. barbadense inherited more NBS-encoding genes from G. raimondii. This asymmetric evolution of NBS-encoding genes may help to explain why G. raimondii and G. barbadense are more resistant to Verticillium wilt, whereas G. arboreum and G. hirsutum are more susceptible to Verticillium wilt. The disease resistances of the allotetraploid cotton were related to their NBS-encoding genes especially in regard from which diploid progenitor they were derived, and the TNL genes may have a significant role in disease resistance to Verticillium wilt in G. raimondii and G. barbadense.
Project description:Acetyl-CoA carboxylase is an important enzyme, which catalyzes acetyl-CoA's carboxylation to produce malonyl-CoA and to serve as a committed step for de novo fatty acid biosynthesis in plastids. In this study, 24 putative cotton BCCP genes were identified based on the lately published genome data in Gossypium. Among them, 4, 4, 8, and 8 BCCP homologs were identified in Gossypium raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. These genes were divided into two classes based on a phylogenetic analysis. In each class, these homologs were relatively conserved in gene structure and motifs. The chromosomal distribution pattern revealed that all the BCCP genes were distributed equally on corresponding chromosomes or scaffold in the four cotton species. Segmental duplication was a predominant duplication event in both of G. hirsutum and G. barbadense. The analysis of the expression profile showed that 8 GhBCCP genes expressed in all the tested tissues with changed expression levels, and GhBCCP genes belonging to class II were predominantly expressed in developing ovules. Meanwhile, the expression analysis for the 16 cotton BCCP genes from G. raimondii, G. arboreum and G. hirsutum showed that they were induced or suppressed by cold or salt stress, and their expression patterns varied among different tissues. These findings will help to determine the functional and evolutionary characteristics of the BCCP genes in Gossypium species.
Project description:Intron length polymorphisms (ILPs), a type of gene-based functional marker, could themselves be related to the particular traits. Here, we developed a genome-wide cotton ILPs based on orthologs annotation from two sequenced diploid species, A-genome Gossypium arboreum and D-genome G. raimondii. We identified 10,180 putative ILP markers from 5,021 orthologous genes. Among these, 535 ILP markers from 9 gene families related to stress were selected for experimental verification. Polymorphic rates were 72.71% between G. arboreum and G. raimondii and 36.45% between G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124. Furthermore, 14 polymorphic ILP markers were detected in 264 G. hirsutum accessions. Coupled with previous simple sequence repeats (SSRs) evaluations and salt tolerance assays from the same individuals, we found a total of 25 marker-trait associations involved in nine ILPs. The nine genes, temporally named as C1 to C9, showed the various expressions in different organs and tissues, and five genes (C3, C4, C5, C7 and C9) were significantly upregulated after salt treatment. We verified that the five genes play important roles in salt tolerance. Particularly, silencing of C4 (encodes WRKY DNA-binding protein) and C9 (encodes Mitogen-activated protein kinase) can significantly enhance cotton susceptibility to salt stress.
Project description:Plant P-type H?-ATPase (P-ATPase) is a membrane protein existing in the plasma membrane that plays an important role in the transmembrane transport of plant cells. To understand the variety and quantity of P-ATPase proteins in different cotton species, we combined four databases from two diploid cotton species (Gossypium raimondii and G. arboreum) and two tetraploid cotton species (G. hirsutum and G. barbadense) to screen the P-ATPase gene family and resolved the evolutionary relationships between the former cotton species. We identified 53, 51, 99 and 98 P-ATPase genes from G. arboretum, G. raimondii, G. barbadense and G. hirsutum, respectively. The structural and phylogenetic analyses revealed that the gene structure was consistent between P-ATPase genes, with a close evolutionary relationship. The expression analysis of P-ATPase genes showed that many P-ATPase genes were highly expressed in various tissues and at different fiber developmental stages in G. hirsutum, suggesting that they have potential functions during growth and fiber development in cotton.
Project description:The protein phosphatase (PP2C) gene family, known to participate in cellular processes, is one of the momentous and conserved plant-specific gene families that regulate signal transduction in eukaryotic organisms. Recently, PP2Cs were identified in Arabidopsis and various other crop species, but analysis of PP2C in cotton is yet to be reported. In the current research, we found 87 (Gossypium arboreum), 147 (Gossypium barbadense), 181 (Gossypium hirsutum), and 99 (Gossypium raimondii) PP2C-encoding genes in total from the cotton genome. Herein, we provide a comprehensive analysis of the PP2C gene family in cotton, such as gene structure organization, gene duplications, expression profiling, chromosomal mapping, protein motif organization, and phylogenetic relationships of each species. Phylogenetic analysis further categorized PP2C genes into 12 subgroups based on conserved domain composition analysis. Moreover, we observed a strong signature of purifying selection among duplicated pairs (i.e., segmental and dispersed) of Gossypium hirsutum. We also observed the tissue-specific response of GhPP2C genes in organ and fiber development by comparing the RNA-sequence (RNA-seq) data reported on different organs. The qRT-PCR validation of 30 GhPP2C genes suggested their critical role in cotton by exposure to heat, cold, drought, and salt stress treatments. Hence, our findings provide an overview of the PP2C gene family in cotton based on various bioinformatic tools that demonstrated their critical role in organ and fiber development, and abiotic stress tolerance, thereby contributing to the genetic improvement of cotton for the resistant cultivar.
Project description:The Dof (DNA-binding one zinc finger) transcription factor family is a representative of plant-specific classes of transcription factors. In this study, we performed a genome-wide screening and characterization of the Dof gene family within two tetraploid species Gossypium barbadense, Gossypium hirsutum, and two diploid species Gossypium arboreum, Gossypium raimondii. 115, 116, 55 and 56 Dof genes were identified respectively and all of the genes contain a sequence encoding the Dof DNA-binding domain. Those genes were unevenly distributed across 13/26 chromosomes of the cotton. Genome comparison revealed that segmental duplication may have played crucial roles in the expansion of the cotton Dof gene family, and tandem duplication also played a minor role. Analysis of RNA-Seq data indicated that cotton Dof gene expression levels varied across different tissues and in response to different abiotic stress. Overall, our results could provide valuable information for better understanding the evolution of cotton Dof genes, and lays a foundation for future investigation in cotton.
Project description:BACKGROUND:Plant Na+/H+ antiporters (NHXs) are membrane-localized proteins that maintain cellular Na+/K+ and pH homeostasis. Considerable evidence highlighted the critical roles of NHX family in plant development and salt response; however, NHXs in cotton are rarely studied. RESULTS:The comprehensive and systematic comparative study of NHXs in three Gossypium species was performed. We identified 12, 12, and 23 putative NHX proteins from G. arboreum, G. raimondii, and G. hirsutum, respectively. Phylogenetic study revealed that repeated polyploidization of Gossypium spp. contributed to the expansion of NHX family. Gene structure analysis showed that cotton NHXs contain many introns, which will lead to alternative splicing and help plants to adapt to high salt concentrations in soil. The expression changes of NHXs indicate the possible differences in the roles of distinct NHXs in salt response. GhNHX1 was proved to be located in the vacuolar system and intensively induced by salt stress in cotton. Silencing of GhNHX1 resulted in enhanced sensitivity of cotton seedlings to high salt concentrations, which suggests that GhNHX1 positively regulates cotton tolerance to salt stress. CONCLUSION:We characterized the gene structure, phylogenetic relationship, chromosomal location, and expression pattern of NHX genes from G. arboreum, G. raimondii, and G. hirsutum. Our findings indicated that the cotton NHX genes are regulated meticulously and differently at the transcription level with possible alternative splicing. The tolerance of plants to salt stress may rely on the expression level of a particular NHX, rather than the number of NHXs in the genome. This study could provide significant insights into the function of plant NHXs, as well as propose promising candidate genes for breeding salt-resistant cotton cultivars.
Project description: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 (AD? 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.
Project description:Metacaspases (MCs) are cysteine proteases that are important for programmed cell death (PCD) in plants. In this study, we identified 89 MC genes in the genomes of four Gossypium species (Gossypium raimondii, Gossypium barbadense, Gossypium hirsutum, and Gossypium arboreum), and classified them as type-I or type-II genes. All of the type-I and type-II MC genes contain a sequence encoding the peptidase C14 domain. During developmentally regulated PCD, type-II MC genes may play an important role related to fiber elongation, while type-I genes may affect the thickening of the secondary wall. Additionally, 13 genes were observed to be differentially expressed between two cotton lines with differing fiber strengths, and four genes (GhMC02, GhMC04, GhMC07, and GhMC08) were predominantly expressed in cotton fibers at 5-30 days post-anthesis (DPA). During environmentally induced PCD, the expression levels of four genes were affected in the root, stem, and leaf tissues within 6 h of an abiotic stress treatment. In general, the MC gene family affects the development of cotton fibers, including fiber elongation and fiber thickening while four prominent fiber- expressed genes were identified. The effects of the abiotic stress and hormone treatments imply that the cotton MC gene family may be important for fiber development. The data presented herein may form the foundation for future investigations of the MC gene family in Gossypium species.
Project description:The protein kinase (PK, kinome) family is one of the largest families in plants and regulates almost all aspects of plant processes, including plant development and stress responses. Despite their important functions, comprehensive functional classification, evolutionary analysis and expression patterns of the cotton PK gene family has yet to be performed on PK genes. In this study, we identified the cotton kinomes in the Gossypium raimondii, Gossypium arboretum, Gossypium hirsutum and Gossypium barbadense genomes and classified them into 7 groups and 122-24 subfamilies using software HMMER v3.0 scanning and neighbor-joining (NJ) phylogenetic analysis. Some conserved exon-intron structures were identified not only in cotton species but also in primitive plants, ferns and moss, suggesting the significant function and ancient origination of these PK genes. Collinearity analysis revealed that 16.6 million years ago (Mya) cotton-specific whole genome duplication (WGD) events may have played a partial role in the expansion of the cotton kinomes, whereas tandem duplication (TD) events mainly contributed to the expansion of the cotton RLK group. Synteny analysis revealed that tetraploidization of G. hirsutum and G. barbadense contributed to the expansion of G. hirsutum and G. barbadense PKs. Global expression analysis of cotton PKs revealed stress-specific and fiber development-related expression patterns, suggesting that many cotton PKs might be involved in the regulation of the stress response and fiber development processes. This study provides foundational information for further studies on the evolution and molecular function of cotton PKs.