CottonGen: a genomics, genetics and breeding database for cotton research

Identification of cotton PIP5K genes and role of GhPIP5K9a in primary root development

Phosphatidylinositol 4 phosphate 5-kinase (PIP5K) is crucial for the phosphatidylinositol (PI) signaling pathway. It plays a significant role in plant growth and development, as well as stress response. However, its effects on cotton are unknown. This study identified PIP5K genes from four cotton species and conducted bioinformatic analyses, with a particular emphasis on the functions of GhPIP5K9a in primary roots. The results showed that cotton PIP5Ks were classified into four subgroups. Analysis of gene structure and motif composition showed obvious conservation within each subgroup. Synteny analysis suggested that the PIP5K gene family experienced significant expansion due to both whole-genome duplication (WGD) and segmental duplication. Transcriptomic data analysis revealed that the majority of GhPIP5K genes had the either low or undetectable levels of expression. Moreover, GhPIP5K9a is highly expressed in the root and was located in plasmalemma. Suppression of GhPIP5K9a transcripts resulted in longer primary roots, longer primary root cells and increased auxin polar transport-related genes expression, and decreased abscisic acid (ABA) content, indicating that GhPIP5K9a negatively regulates cotton primary root growth. This study lays the foundation for further exploration of the role of the PIP5K genes in cotton.

Characteristic of GEX1 genes reveals the essential roles for reproduction in cotton

GEX1 (gamete expressed 1) proteins are critical membrane proteins conserved among flowering plants that are involved in the nuclear fusion and embryonic development. Herein, we identified the 32 GEX1 proteins from representative land plants. In cotton, GEX1 genes expressed in various tissues across all stages of the life cycle, especially in pollen. Subcellular localization indicated the position of GhGEX1 protein was localized in the endoplasmic reticulum. Experimental research has demonstrated that GhGEX1 has the potential to improve the partial abortion phenotype in Arabidopsis. CRISPR/Cas9-mediated knockout of GhGEX1 exhibited the seed abortion. Paraffin section of the ovule revealed that the polar nuclear fusion of ghgex1 plants remains at a standstill when the wild type has developed into a normal embryo. Comparative transcriptome analysis showed that the DEGs of reproductive-related processes and membrane-related processes were repressed in the pollen of knockout lines. The predicted protein interactions showed that GhGEX1 probably functioned through interactions with proteins related to reproduction and membrane. From all these investigations, it was possible to conclude that the GEX1 proteins are evolutionarily conserved in flowering plants and elucidated the pivotal roles during fertilization and early embryonic development in cotton.

Genotype and phenotype data standardization, utilization and integration in the big data era for agricultural sciences

Large-scale genotype and phenotype data have been increasingly generated to identify genetic markers, understand gene function and evolution and facilitate genomic selection. These datasets hold immense value for both current and future studies, as they are vital for crop breeding, yield improvement and overall agricultural sustainability. However, integrating these datasets from heterogeneous sources presents significant challenges and hinders their effective utilization. We established the Genotype-Phenotype Working Group in November 2021 as a part of the AgBioData Consortium (https://www.agbiodata.org) to review current data types and resources that support archiving, analysis and visualization of genotype and phenotype data to understand the needs and challenges of the plant genomic research community. For 2021-22, we identified different types of datasets and examined metadata annotations related to experimental design/methods/sample collection, etc. Furthermore, we thoroughly reviewed publicly funded repositories for raw and processed data as well as secondary databases and knowledgebases that enable the integration of heterogeneous data in the context of the genome browser, pathway networks and tissue-specific gene expression. Based on our survey, we recommend a need for (i) additional infrastructural support for archiving many new data types, (ii) development of community standards for data annotation and formatting, (iii) resources for biocuration and (iv) analysis and visualization tools to connect genotype data with phenotype data to enhance knowledge synthesis and to foster translational research. Although this paper only covers the data and resources relevant to the plant research community, we expect that similar issues and needs are shared by researchers working on animals. Database URL: https://www.agbiodata.org.

Glyceraldehyde-3-phosphate dehydrogenase Gh_GAPDH9 is associated with drought resistance in Gossypium hirsutum

Background

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is the central enzyme of glycolysis and plays important regulatory roles in plant growth and development and responses to adverse stress conditions. However, studies on the characteristics and functions of cotton GAPDH family genes are still lacking.

Methods

In this study, genome-wide identification of the cotton GAPDH gene family was performed, and the phylogeny, gene structures, promoter progenitors and expression profiles of upland cotton GAPDH gene family members were explored by bioinformatics analysis to highlight potential functions. The functions of GhGAPDH9 in response to drought stress were initially validated based on RNA-seq, qRT‒PCR, VIGS techniques and overexpression laying a foundation for further studies on the functions of GAPDH genes.

Results

This study is the first systematic analysis of the cotton GAPDH gene family, which contains a total of 84 GAPDH genes, among which upland cotton contains 27 members. Quantitative, phylogenetic and covariance analyses of the genes revealed that the GAPDH gene family has been conserved during the evolution of cotton. Promoter analysis revealed that most cis-acting elements were related to MeJA and ABA. Based on the identified promoter cis-acting elements and RNA-seq data, it was hypothesized that Gh_GAPDH9, Gh_GAPDH11, Gh_GAPDH19 and Gh_GAPDH21 are involved in the response of cotton to abiotic stress. The expression levels of the Gh_GAPDH9 gene in two drought-resistant and two drought-sensitive materials were analyzed by qRT‒PCR and found to be high early in the treatment period in the drought-resistant material. The silencing of Gh_GAPDH9 based on virus-induced gene silencing (VIGS) technology resulted in significant leaf wilting or whole-plant dieback in silenced plants after drought stress compared to the control. The content of-malondialdehyde (MDA) in cotton leaves was significantly increased, and the content of proline (Pro) and chlorophyll (Chl) was reduced. In addition, the leaf wilting and dryness of transgenic lines under drought stress were lower than those of wild-type Arabidopsis, indicating that Gh_GAPDH9 is a positive regulator of drought resistance. In conclusion, our results demonstrate that GAPDH genes play an important role in the response of cotton to abiotic stresses and provide preliminary validation of the function of the Gh_GAPDH9 gene under drought stress. These findings provide an important theoretical basis for further studies on the function of the Gh_GAPDH9 gene and the molecular mechanism of the drought response in cotton.

Decoding the guardians of cotton resilience: A comprehensive exploration of the βCA genes and its role in Verticillium dahliae resistance

Plant Carbonic anhydrases (Cas) have been shown to be stress-responsive enzymes that may play a role in adapting to adverse conditions. Cotton is a significant economic crop in China, with upland cotton (Gossypium hirsutum) being the most widely cultivated species. We conducted genome-wide identification of the βCA gene in six cotton species and preliminary analysis of the βCA gene in upland cotton. In total, 73 βCA genes from six cotton species were identified, with phylogenetic analysis dividing them into five subgroups. GHβCA proteins were predominantly localized in the chloroplast and cytoplasm. The genes exhibited conserved motifs, with motifs 1, 2, and 3 being prominent. GHβCA genes were unevenly distributed across chromosomes and were associated with stress-responsive cis-regulatory elements, including those responding to light, MeJA, salicylic acid, abscisic acid, cell cycle regulation, and defence/stress. Expression analysis indicated that GHβCA6, GHβCA7, GHβCA10, GHβCA15, and GHβCA16 were highly expressed under various abiotic stress conditions, whereas GHβCA3, GHβCA9, GHβCA10, and GHβCA18 had higher expression patterns under Verticillium dahliae infection at different time intervals. In Gossypium thurberi, GthβCA1, GthβCA2, and GthβCA4 showed elevated expression across stress conditions and tissues. Silencing GHβCA10 through VIGS increased Verticillium wilt severity and reduced lignin deposition compared to non-silenced plants. GHβCA10 is crucial for cotton's defense against Verticillium dahliae. Further research is needed to understand the underlying mechanisms and develop strategies to enhance resistance against Verticillium wilt.

Using association mapping and local interval haplotype association analysis to improve the cotton drought stress response

Drought stress has a serious impact on the growth and development of cotton. To explore the relevant molecular mechanism of the drought stress response in cotton, gene mapping based on the QTL interval mapped by simplified genome BSA-seq of the drought-resistance-related RIL population was performed. A QTL region spanning 2.02 Mb on chromosome D07 was selected, and 201 resource materials were genotyped using 9 KASP markers in the interval. After local interval haplotype association analysis, the overlap of the 110 kb peak region confirmed the reliability of this region, and at the same time, the role of GhGF14-30, the only gene in the overlapping region, was modeled in the response of cotton to drought stress. qRTPCR analysis of the materials and population parents proved that this gene plays a role in the drought stress response in cotton. Virus-induced gene silencing proved the importance of this gene in drought-sensitive materials, and drought-resistance-related marker genes also proved that the GhGF14-30 gene may play an important role in the ABA and SOS signaling pathways. This study provides a basis for mining drought stress response functional genes in cotton and lays the foundation for the molecular mechanism of the GhGF14-30 gene in response to drought stress in cotton.

Analysis of the histidine kinase gene family and the role of GhHK8 in response to drought tolerance in cotton

As an important member of the two-component system (TCS), histidine kinases (HKs) play important roles in various plant developmental processes and signal transduction in response to a wide range of biotic and abiotic stresses. So far, the HK gene family has not been investigated in Gossypium. In this study, a total of 177 HK gene family members were identified in cotton. They were further divided into seven groups, and the protein characteristics, genetic relationship, gene structure, chromosome location, collinearity, and cis-elements identification were comprehensively analyzed. Whole genome duplication (WGD) / segmental duplication may be the reason why the number of HK genes doubled in tetraploid Gossypium species. Expression analysis revealed that most cotton HK genes were mainly expressed in the reproductive organs and the fiber at initial stage. Gene expression analysis revealed that HK family genes are involved in cotton abiotic stress, especially drought stress and salt stress. In addition, gene interaction networks showed that HKs were involved in the regulation of cotton abiotic stress, especially drought stress. VIGS experiments have shown that GhHK8 is a negative regulatory factor in response to drought stress. Our systematic analysis provided insights into the characteristics of the HK genes in cotton and laid a foundation for further exploring their potential in drought stress resistance in cotton.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development

Introduction

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels.

Methods

Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes.

Results

Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers.

Discussion

Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

Alternative Splicing during Fiber Development in G. hirsutum

Cotton is a valuable cash crop in many countries. Cotton fiber is a trichome that develops from a single epidermal cell and serves as an excellent model for understanding cell differentiation and other life processes. Alternative splicing (AS) of genes is a common post-transcriptional regulatory process in plants that is essential for plant growth and development. The process of AS during cotton fiber formation, on the other hand, is mainly unknown. A substantial number of multi-exon genes were discovered to be alternatively spliced during cotton fiber formation in this study, accounting for 23.31% of the total number of genes in Gossypium hirsutum. Retention intron (RI) is not necessarily the most common AS type, indicating that AS genes and processes during fiber development are very temporal and tissue-specific. When compared to fiber samples, AS is more prevalent at the fiber initiation stages and in the ovule, indicating that development stages and tissues use different AS strategies. Genes involved in fiber development have gone through stage-specific AS, demonstrating that AS regulates cotton fiber development. Furthermore, AS can be regulated by trans-regulation elements such as splicing factor and cis-regulation elements such as gene length, exon numbers, and GC content, particularly at exon-intron junction sites. Our findings also suggest that increased DNA methylation may aid in the efficiency of AS, and that gene body methylation is key in AS control. Finally, our research will provide useful information about the roles of AS during the cotton fiber development process.

GhSTP18, a member of sugar transport proteins family, negatively regulates salt stress in cotton

The sugar transporter protein (STP) family has been shown to play important roles in plant growth, development, and stress response. However, it has not been studied in cotton compared to other major crops. In this study, we identified 90 STP genes from four cotton species, performed bioinformatic analysis, and focused on the role of GhSTP18 in salt stress. According to our results, cotton STP proteins were divided into four subgroups according to the phylogenetic tree. A synteny analysis suggested that whole-genome duplication (WGD) and segmental duplication were key drivers in the expansion of the STP gene family. The transcriptomic data analysis showed that 29 GhSTP genes exhibited sink-specific expression. Quantitative real time-polymerase chain reaction (qRT-PCR) analyses revealed that expression of GhSTP18 was induced by salt treatment, heat treatment, cold treatment, and drought treatment, and continuously increased during a salt stress time course. Notably, GhSTP18 encodes a plasma membrane-localized galactose transporter. Suppression of GhSTP18 transcription by a virus-induced gene silencing (VIGS) assay reduced sensitivity to salt stress in cotton, indicating that GhSTP18 negatively regulates plant salt tolerance. These results provide an important reference and resource for further studying and deploying STP genes for cotton improvement.

Genome-wide identification of xylan glucuronosyltransferase family in cotton and function characterization of GhGUX5 in regulating Verticillium wilt resistance

Xylan glucuronosyltransferase (GUX) is widely involved in a variety of physiological processes in plants, including plant development, growth and the defense response to pathogens. However, the function of GUX regulators in Verticillium dahliae (V. dahliae) infection has not been considered previously in cotton. Overall, 119 GUX genes were identified from multiple species and were phylogenetically categorized into seven classes. Duplication event analysis indicated that GUXs in Gossypium hirsutum primarily originated from segmental duplication. GhGUXs promoter analysis indicated cis-regulatory elements capable of reacting to several different stresses. RNA-Seq data and qRT-PCR analysis both indicated that most GhGUXs were associated with V. dahliae infection. Gene interaction network analysis showed that GhGUX5 interacted with 11 proteins, and the relative expression of these 11 proteins changed significantly following V. dahliae infection. In addition, silencing and overexpression of GhGUX5 results to enhance and reduce plant's susceptibility to V. dahliae. Further study showed that TRV: GhGUX5 silenced cotton plants exhibited a decrease in the degree of lignification, total lignin content, gene expression levels involved in lignin biosynthesis, and enzyme activity compared with TRV: 00. The above results indicate that GhGUX5 enhances Verticillium wilt resistance through the lignin biosynthesis pathway.

MDSi: Multi-omics Database for Setaria italica

BACKGROUND

Foxtail millet (Setaria italica) harbors the small diploid genome (~ 450 Mb) and shows the high inbreeding rate and close relationship to several major foods, feed, fuel and bioenergy grasses. Previously, we created a mini foxtail millet, xiaomi, with an Arabidopsis-like life cycle. The de novo assembled genome data with high-quality and an efficient Agrobacterium-mediated genetic transformation system made xiaomi an ideal C₄ model system. The mini foxtail millet has been widely shared in the research community and as a result there is a growing need for a user-friendly portal and intuitive interface to perform exploratory analysis of the data.

RESULTS

Here, we built a Multi-omics Database for Setaria italica (MDSi, http://sky.sxau.edu.cn/MDSi.htm ), that contains xiaomi genome of 161,844 annotations, 34,436 protein-coding genes and their expression information in 29 different tissues of xiaomi (6) and JG21 (23) samples that can be showed as an Electronic Fluorescent Pictograph (xEFP) in-situ. Moreover, the whole-genome resequencing (WGS) data of 398 germplasms, including 360 foxtail millets and 38 green foxtails and the corresponding metabolic data were available in MDSi. The SNPs and Indels of these germplasms were called in advance and can be searched and compared in an interactive manner. Common tools including BLAST, GBrowse, JBrowse, map viewer, and data downloads were implemented in MDSi.

CONCLUSION

The MDSi constructed in this study integrated and visualized data from three levels of genomics, transcriptomics and metabolomics, and also provides information on the variation of hundreds of germplasm resources that can satisfies the mainstream requirements and supports the corresponding research community.

Identification of fertility restoration candidate genes from a restorer line R186 for Gossypium harknessii cytoplasmic male sterile cotton

BACKGROUND

The utilization of heterosis based on three-line system is an effective strategy in crop breeding. However, cloning and mechanism elucidation of restorer genes for cytoplasmic male sterility (CMS) in upland cotton have yet been realized.

RESULTS

This research is based on CMS line 2074A with the cytoplasm from Gossypium harknessii (D_2-2) and restorer line R186. The offspring of 2074A × R186 were used to conduct genetic analysis. The fertility mechanism of 2074A can be speculated to be governed by multiple genes, since neither the single gene model nor the double genes model could be used. The bulked segregant analysis (BSA) for (2074A × R186) F₂ determined the genetic interval of restorer genes on a region of 4.30 Mb on chromosome D05 that contains 77 annotated genes. Four genes were identified as candidates for fertility restoration using the RNA-seq data of 2074A, 2074B, and R186. There are a number of large effect variants in the four genes between 2074A and R186 that could cause amino acid changes. Evolutionary analysis and identity analysis revealed that GH_D05G3183, GH_D05G3384, and GH_D05G3490 have high identity with their homologs in D_2-2, respectively. Tissue differential expression analysis revealed that the genes GH_D05G3183, GH_D05G3384, and GH_D05G3490 were highly expressed in the buds of the line R186. The predicted results demonstrated that GH_D05G3183, GH_D05G3384 and GH_D05G3490 might interact with GH_A02G1295 to regulate orf610a in mitochondria.

CONCLUSION

Our study uncovered candidate genes for fertility restoration in the restorer line R186 and predicted the possible mechanism for restoring the male fertility in 2074A. This research provided valuable insight into the nucleoplasmic interactions.

Functional characterization of the GhNRT2.1e gene reveals its significant role in improving nitrogen use efficiency in Gossypium hirsutum

Background

Nitrate is the primary type of nitrogen available to plants, which is absorbed and transported by nitrate transporter 2 (NRT2) at low nitrate conditions.

Methods

Genome-wide identification of NRT2 genes in G. hirsutum was performed. Gene expression patterns were revealed using RNA-seq and qRT-PCR. Gene functions were characterized using overexpression in A. thaliana and silencing in G. hirsutum. Protein interactions were verified by yeast two-hybrid and luciferase complementation imaging (LCI) assays.

Results

We identified 14, 14, seven, and seven NRT2 proteins in G. hirsutum, G. barbadense, G. raimondii, and G. arboreum. Most NRT2 proteins were predicted in the plasma membrane. The NRT2 genes were classified into four distinct groups through evolutionary relationships, with members of the same group similar in conserved motifs and gene structure. The promoter regions of NRT2 genes included many elements related to growth regulation, phytohormones, and abiotic stresses. Tissue expression pattern results revealed that most GhNRT2 genes were specifically expressed in roots. Under low nitrate conditions, GhNRT2 genes exhibited different expression levels, with GhNRT2.1e being the most up-regulated. Arabidopsis plants overexpressing GhNRT2.1e exhibited increased biomass, nitrogen and nitrate accumulation, nitrogen uptake and utilization efficiency, nitrogen-metabolizing enzyme activity, and amino acid content under low nitrate conditions. In addition, GhNRT2.1e-silenced plants exhibited suppressed nitrate uptake and accumulation, hampered plant growth, affected nitrogen metabolism processes, and reduced tolerance to low nitrate. The results showed that GhNRT2.1e could promote nitrate uptake and transport under low nitrate conditions, thus effectively increasing nitrogen use efficiency (NUE). We found that GhNRT2.1e interacts with GhNAR2.1 by yeast two-hybrid and LCI assays.

Discussion

Our research lays the foundation to increase NUE and cultivate new cotton varieties with efficient nitrogen use.

A mutant cotton fatty acid desaturase 2-1d allele causes protein mistargeting and altered seed oil composition

BACKGROUND

Cotton (Gossypium sp.) has been cultivated for centuries for its spinnable fibers, but its seed oil also possesses untapped economic potential if, improvements could be made to its oleic acid content.

RESULTS

Previous studies, including those from our laboratory, identified pima accessions containing approximately doubled levels of seed oil oleic acid, compared to standard upland cottonseed oil. Here, the molecular properties of a fatty acid desaturase encoded by a mutant allele identified by genome sequencing in an earlier analysis were analyzed. The mutant sequence is predicted to encode a C-terminally truncated protein lacking nine residues, including a predicted endoplasmic reticulum membrane retrieval motif. We determined that the mutation was caused by a relatively recent movement of a Ty1/copia type retrotransposon that is not found associated with this desaturase gene in other sequenced cotton genomes. The mutant desaturase, along with its repaired isozyme and the wild-type A-subgenome homoeologous protein were expressed in transgenic yeast and stably transformed Arabidopsis plants. All full-length enzymes efficiently converted oleic acid to linoleic acid. The mutant desaturase protein produced only trace amounts of linoleic acid, and only when strongly overexpressed in yeast cells, indicating that the missing C-terminal amino acid residues are not strictly required for enzyme activity, yet are necessary for proper subcellular targeting to the endoplasmic reticulum membrane.

CONCLUSION

These results provide the biochemical underpinning that links a genetic lesion present in a limited group of South American pima cotton accessions and their rare seed oil oleic acid traits. Markers developed to the mutant desaturase allele are currently being used in breeding programs designed to introduce this trait into agronomic upland cotton varieties.

Genome-wide identification and expression analysis of the HD2 protein family and its response to drought and salt stress in Gossypium species

Histone deacetylase 2 (HD2) proteins play an important role in the regulation of gene expression. This helps with the growth and development of plants and also plays a crucial role in responses to biotic and abiotic stress es. HD2s comprise a C₂H₂-type Zn²⁺ finger at their C-terminal and an HD2 label, deacetylation and phosphorylation sites, and NLS motifs at their N-terminal. In this study, a total of 27 HD2 members were identified, using Hidden Markov model profiles, in two diploid cotton genomes (Gossypium raimondii and Gossypium arboretum) and two tetraploid cotton genomes (Gossypium hirsutum and Gossypium barbadense). These cotton HD2 members were classified into 10 major phylogenetic groups (I-X), of which group III was found to be the largest with 13 cotton HD2 members. An evolutionary investigation showed that the expansion of HD2 members primarily occurred as a result of segmental duplication in paralogous gene pairs. Further qRT-PCR validation of nine putative genes using RNA-Seq data suggested that GhHDT3D.2 exhibits significantly higher levels of expression at 12h, 24h, 48h, and 72h of exposure to both drought and salt stress conditions compared to a control measure at 0h. Furthermore, gene ontology, pathways, and co-expression network study of GhHDT3D.2 gene affirmed their significance in drought and salt stress responses.

Flame resistant cotton lines generated by synergistic epistasis in a MAGIC population

Textiles made from cotton fibers are flammable and thus often include flame retardant additives for consumer safety. Transgressive segregation in multi-parent populations facilitates new combinations of alleles of genes and can result in traits that are superior to those of any of the parents. A screen of 257 recombinant inbred lines from a multi-parent advanced generation intercross (MAGIC) population for naturally enhance flame retardance (FR) was conducted. All eleven parents, like all conventional white fiber cotton cultivars produce flammable fabric. MAGIC recombinant inbred lines (RILs) that produced fibers with significantly lower heat release capacities (HRC) as measured by microscale combustion calorimetry (MCC) were identified and the stability of the phenotypes of the outliers were confirmed when the RILs were grown at an additional location. Of the textiles fabricated from the five superior RILs, four exhibited the novel characteristic of inherent flame resistance. When exposed to open flame by standard 45° incline flammability testing, these four fabrics self-extinguished. To determine the genetic architecture of this novel trait, linkage, epistatic and multi-locus genome wide association studies (GWAS) were conducted with 473k SNPs identified by whole genome sequencing (WGS). Transcriptomes of developing fiber cells from select RILs were sequenced (RNAseq). Together, these data provide insight into the genetic mechanism of the unexpected emergence of flame-resistant cotton by transgressive segregation in a breeding program. The incorporation of this trait into global cotton germplasm by breeding has the potential to greatly reduce the costs and impacts of flame-retardant chemicals.

Whole-Genome Resequencing Deciphers New Insight Into Genetic Diversity and Signatures of Resistance in Cultivated Cotton Gossypium hirsutum

Cotton is an important crop that produces fiber and cottonseed oil for the textile and oil industry. However, cotton leaf curl virus disease (CLCuD) stress is limiting its yield in several Asian countries. In this study, we have sequenced Mac7 accession, a Gossypium hirsutum resistance source against several biotic stresses. By aligning with the Gossypium hirsutum (AD1) 'TM-1' genome, a total of 4.7 and 1.2 million SNPs and InDels were identified in the Mac7 genome. The gene ontology and metabolic pathway enrichment indicated SNPs and InDels role in nucleotide bindings, secondary metabolite synthesis, and plant-pathogen interaction pathways. Furthermore, the RNA-seq data in different tissues and qPCR expression profiling under CLCuD provided individual gene roles in resistant and susceptible accessions. Interestingly, the differential NLR genes demonstrated higher expression in resistant plants rather than in susceptible plants expression. The current resequencing results may provide primary data to identify DNA resistance markers which will be helpful in marker-assisted breeding for development of Mac7-derived resistance lines.

Dual Domestication, Diversity, and Differential Introgression in Old World Cotton Diploids

Domestication in the cotton genus is remarkable in that it has occurred independently four different times at two different ploidy levels. Relatively little is known about genome evolution and domestication in the cultivated diploid species Gossypium herbaceum and Gossypium arboreum, due to the absence of wild representatives for the latter species, their ancient domestication, and their joint history of human-mediated dispersal and interspecific gene flow. Using in-depth resequencing of a broad sampling from both species, we provide support for their independent domestication, as opposed to a progenitor-derivative relationship, showing that diversity (mean π = 6 × 10-3) within species is similar, and that divergence between species is modest (FST = 0.413). Individual accessions were homozygous for ancestral single-nucleotide polymorphisms at over half of variable sites, while fixed, derived sites were at modest frequencies. Notably, two chromosomes with a paucity of fixed, derived sites (i.e., chromosomes 7 and 10) were also strongly implicated as having experienced high levels of introgression. Collectively, these data demonstrate variable permeability to introgression among chromosomes, which we propose is due to divergent selection under domestication and/or the phenomenon of F2 breakdown in interspecific crosses. Our analyses provide insight into the evolutionary forces that shape diversity and divergence in the diploid cultivated species and establish a foundation for understanding the contribution of introgression and/or strong parallel selection to the extensive morphological similarities shared between species.

GhALKBH10 negatively regulates salt tolerance in cotton

The alpha-ketoglutarate-dependent dioxygenase (AlkB) gene family plays an essential role in regulating plant development and stress response. However, the AlkB gene family is still not well understood in cotton. In this study, 40 AlkB genes in cotton and Arabidopsis are identified and classified into three classes based on phylogenetic analysis. Their protein motifs and exon/intron structures are highly conserved. Chromosomal localization and synteny analysis suggested that segmental or whole-genome duplication and polyploidization events contributed to the expansion of the cotton AlkB gene family. Furthermore, the AlkB genes showed dynamic spatiotemporal expression patterns and diverse responses to abiotic stresses. Among them, GhALKBH10 was down-regulated under various abiotic stresses and its subcellular expression was localized in cytoplasm and nucleus. Silencing GhALKBH10 in cotton increased antioxidant capacity and reduced cytoplasmic Na+ concentration, thereby improved the plant tolerance to salinity. Conversely, overexpression (OE) of GhALKBH10 in Arabidopsis markedly weakened the plant tolerance to salinity. The global m6A levels measured in VIGS and OE transgenic lines showed that they were significantly higher in TRV: GhALKBH10 plants (VIGS) than in TRV: 00 plants but significantly lower in OE plants than wild-type plants under salt stress, which could be considered as a potential m6A demethylase in cotton. Our results suggest that the GhALKBH10 gene negatively regulates salt tolerance in plants, which provides information of the cotton AlkB family and an understanding of GhALKBH10 function under salt condition as well as a new gene for salt-tolerant cotton breeding.

Transcriptomic analysis reveals the key role of histone deacetylation via mediating different phytohormone signalings in fiber initiation of cotton

Background

Histone deacetylation is one of the most important epigenetic modifications and plays diverse roles in plant development. However, the detailed functions and mechanisms of histone deacetylation in fiber development of cotton are still unclear. HDAC inhibitors (HDACi) have been commonly used to study the molecular mechanism underlying histone deacetylation or to facilitate disease therapy in humans through hindering the histone deacetylase catalytic activity. Trichostatin A (TSA)—the most widely used HDACi has been extensively employed to determine the role of histone deacetylation on different developmental stages of plants.

Results

Through in vitro culture of ovules, we observed that exogenous application of TSA was able to inhibit the fiber initiation development. Subsequently, we performed a transcriptomic analysis to reveal the underlying mechanisms. The data showed that TSA treatment resulted in 4209 differentially expressed genes, which were mostly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, and carbon metabolism pathways. The phytohormone signal transduction pathways harbor the most differentially expressed genes. Deeper studies showed that some genes promoting auxin, Gibberellic Acid (GA) signaling were down-regulated, while some genes facilitating Abscisic Acid (ABA) and inhibiting Jasmonic Acid (JA) signaling were up-regulated after the TSA treatments. Further analysis of plant hormone contents proved that TSA significantly promoted the accumulation of ABA, JA and GA₃.

Conclusions

Collectively, histone deacetylation can regulate some key genes involved in different phytohormone pathways, and consequently promoting the auxin, GA, and JA signaling, whereas repressing the ABA synthesis and signaling to improve the fiber cell initiation. Moreover, the genes associated with energy metabolism, phenylpropanoid, and glutathione metabolism were also regulated by histone deacetylation. The above results provided novel clues to illuminate the underlying mechanisms of epigenetic modifications as well as related different phytohormones in fiber cell differentiation, which is also very valuable for the molecular breeding of higher quality cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00840-4.

Systematic analysis of Histidine photosphoto transfer gene family in cotton and functional characterization in response to salt and around tolerance

BACKGROUND

Phosphorylation regulated by the two-component system (TCS) is a very important approach signal transduction in most of living organisms. Histidine phosphotransfer (HP) is one of the important members of the TCS system. Members of the HP gene family have implications in plant stresses tolerance and have been deeply studied in several crops. However, upland cotton is still lacking with complete systematic examination of the HP gene family.

RESULTS

A total of 103 HP gene family members were identified. Multiple sequence alignment and phylogeny of HPs distributed them into 7 clades that contain the highly conserved amino acid residue "XHQXKGSSXS", similar to the Arabidopsis HP protein. Gene duplication relationship showed the expansion of HP gene family being subjected with whole-genome duplication (WGD) in cotton. Varying expression profiles of HPs illustrates their multiple roles under altering environments particularly the abiotic stresses. Analysis is of transcriptome data signifies the important roles played by HP genes against abiotic stresses. Moreover, protein regulatory network analysis and VIGS mediated functional approaches of two HP genes (GhHP23 and GhHP27) supports their predictor roles in salt and drought stress tolerance.

CONCLUSIONS

This study provides new bases for systematic examination of HP genes in upland cotton, which formulated the genetic makeup for their future survey and examination of their potential use in cotton production.

Comparative Genomics and Functional Studies of Putative m(6)A Methyltransferase (METTL) Genes in Cotton

N6-methyladenosine (m6A) RNA modification plays important regulatory roles in plant development and adapting to the environment, which requires methyltransferases to achieve the methylation process. However, there has been no research regarding m6A RNA methyltransferases in cotton. Here, a systematic analysis of the m6A methyltransferase (METTL) gene family was performed on twelve cotton species, resulting in six METTLs identified in five allotetraploid cottons, respectively, and three to four METTLs in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that METTL genes from cottons, Arabidopsis thaliana, and Homo sapiens could be classified into three clades (METTL3, METTL14, and METTL-like clades). Cis-element analysis predicated the possible functions of METTL genes in G. hirsutum. RNA-seq data revealed that GhMETTL14 (GH_A07G0817/GH_D07G0819) and GhMETTL3 (GH_A12G2586/GH_D12G2605) had high expressions in root, stem, leaf, torus, petal, stamen, pistil, and calycle tissues. GhMETTL14 also had the highest expression in 20 and 25 dpa fiber cells, implying a potential role at the cell wall thickening stage. Suppressing GhMETTL3 and GhMETTL14 by VIGS caused growth arrest and even death in G. hirsutum, along with decreased m6A abundance from the leaf tissues of VIGS plants. Overexpression of GhMETTL3 and GhMETTL14 produced distinct differentially expressed genes (DEGs) in A. thaliana, indicating their possible divergent functions after gene duplication. Overall, GhMETTLs play indispensable but divergent roles during the growth of cotton plants, which provides the basis for the systematic investigation of m6A in subsequent studies to improve the agronomic traits in cotton.

Systematical Characterization of the Cotton Di19 Gene Family and the Role of GhDi19-3 and GhDi19-4 as Two Negative Regulators in Response to Salt Stress

Drought-induced 19 (Di19) protein is a Cys2/His2 (C2H2) type zinc-finger protein, which plays a crucial role in plant development and in response to abiotic stress. This study systematically investigated the characteristics of the GhDi19 gene family, including the member number, gene structure, chromosomal distribution, promoter cis-elements, and expression profiles. Transcriptomic analysis indicated that some GhDi19s were up-regulated under heat and salt stress. Particularly, two nuclear localized proteins, GhDi19-3 and GhDi19-4, were identified as being in potential salt stress responsive roles. GhDi19-3 and GhDi19-4 decreased sensitivity under salt stress through virus-induced gene silencing (VIGS), and showed significantly lower levels of H2O2, malondialdehyde (MDA), and peroxidase (POD) as well as significantly increased superoxide dismutase (SOD) activity. This suggested that their abilities were improved to effectively reduce the reactive oxygen species (ROS) damage. Furthermore, certain calcium signaling and abscisic acid (ABA)-responsive gene expression levels showed up- and down-regulation changes in target gene-silenced plants, suggesting that GhDi19-3 and GhDi19-4 were involved in calcium signaling and ABA signaling pathways in response to salt stress. In conclusion, GhDi19-3 and GhDi19-4, two negative transcription factors, were found to be responsive to salt stress through calcium signaling and ABA signaling pathways.

Protein S-Acyl Transferase GhPAT27 Was Associated with Verticillium wilt Resistance in Cotton

Protein palmitoylation is an ability of the frame of the cell marker protein is one of the most notable reversible changes after translation. However, studies on protein palmitoylation in cotton have not yet been performed. In our current research, the PAT gene family was systematically identified and bioinformatically analyzed in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, and 211 PAT genes were authenticated and classified into six subfamilies. Sixty-nine PAT genes were identified in upland cotton, mainly at the ends of its the 26 chromosomes of upland cotton. The majority of these genes are located in the nucleus of the plant. Gene structure analysis revealed that each member encodes a protein that which contains at least one DHHC structural domain. Cis-acting element analysis indicated that GhPATs genes are mainly involved in hormone production, light response and stress response. Gene expression pattern analysis indicated that most GhPATs genes were differentially expressed upon induction by pathogenic bacteria, drought, salt, hot and cold stresses, and some GhPATs could be induced by multiple abiotic stresses simultaneously. GhPATs genes showed different expression patterns in tissue-specific assays and were found to be preferentially expressed in roots, followed by expression in stems and leaves. Virus-induced gene silencing (VIGS) experiments showed that cotton was significantly less resistant to Verticillium dahliae when GhPAT27 was silenced. We conclude that the GhPAT27 gene, which mediates S-palmitoylation acetylation, may be involved in the regulation of upland cotton resistance to Verticillium wilt (VW). Overall, this work has provided a fundamental framework for understanding the latent capabilities of GhPATs and a solid foundation for molecular breeding and plant pathogen resistance in cotton.

CottonMD: a multi-omics database for cotton biological study

Cotton is an important economic crop, and many loci for important traits have been identified, but it remains challenging and time-consuming to identify candidate or causal genes/variants and clarify their roles in phenotype formation and regulation. Here, we first collected and integrated the multi-omics datasets including 25 genomes, transcriptomes in 76 tissue samples, epigenome data of five species and metabolome data of 768 metabolites from four tissues, and genetic variation, trait and transcriptome datasets from 4180 cotton accessions. Then, a cotton multi-omics database (CottonMD, http://yanglab.hzau.edu.cn/CottonMD/) was constructed. In CottonMD, multiple statistical methods were applied to identify the associations between variations and phenotypes, and many easy-to-use analysis tools were provided to help researchers quickly acquire the related omics information and perform multi-omics data analysis. Two case studies demonstrated the power of CottonMD for identifying and analyzing the candidate genes, as well as the great potential of integrating multi-omics data for cotton genetic breeding and functional genomics research.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons

Wild relatives of domesticated plants provide a rich resource for crop improvement and a valuable comparative perspective for understanding genomic, physiological, and agricultural traits. Here, we provide high-quality reference genomes of one early domesticated form of the economically most important cotton species, Gossypium hirsutum, and two other wild species, to clarify evolutionary relationships and understand the genomic changes that characterize these species and their close relatives. We document abundant gene resources involved in adaptation to environmental challenges, highlighting the potential for introgression of favorable genes into domesticated cotton and for increasing resilience to climate variability. Our study complements other recent genomic analyses in the cotton genus and provides a valuable foundation for breeding improved cotton varieties.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)₆, Ge] and Gossypium stephensii [(AD)₇, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)₁, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates—including phenotypic differentiation, genetic isolation, and genetic convergence—that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

COTTONOMICS: a comprehensive cotton multi-omics database

The rapid advancement of sequencing technology, including next-generation sequencing (NGS), has greatly improved sequencing efficiency and decreased cost. Consequently, huge amounts of genomic, transcriptomic and epigenetic data concerning cotton species have been generated and released. These large-scale data provide immense opportunities for the study of cotton genomic structure and evolution, population genetic diversity and genome-wide mining of excellent genes for important traits. However, the complexity of NGS data also causes distress, as it cannot be utilized easily. Here, we presented the cotton omics data platform COTTONOMICS (http://cotton.zju.edu.cn/), an easily accessible web database that integrates 32.5 TB of omics data including seven assembled genomes, resequencing data from 1180 allotetraploid cotton accessions and RNA-sequencing (RNA-seq), small RNA-sequencing (smRNA-seq), Chromatin Immunoprecipitation sequencing (ChIP-seq), DNase hypersensitive sites sequencing (DNase-seq) and Bisulfite sequencing (BS-seq). COTTONOMICS allows users to employ various search scenarios and retrieve information concerning the cotton genomes, genomic variation (Single nucleotide polymorphisms (SNPs) and Insertion and Deletion (InDels)), gene expression, smRNA expression, epigenetic regulation and quantitative trait locus (QTLs). The user-friendly web interface offers a variety of modules for storing, retrieving, analyzing and visualizing cotton multi-omics data to diverse ends, thereby enabling users to decipher cotton population genetics and identify potential novel genes that influence agronomically beneficial traits.

Database URL: http://cotton.zju.edu.cn

Retrieving a disrupted gene encoding phospholipase A for fibre enhancement in allotetraploid cultivated cotton

After polyploidization originated from one interspecific hybridization event in Gossypium, Gossypium barbadense evolved to produce extra-long staple fibres than Gossypium hirsutum (Upland cotton), which produces a higher fibre yield. The genomic diversity between G. barbadense and G. hirsutum thus provides a genetic basis for fibre trait variation. Recently, rapid accumulation of gene disruption or deleterious mutation was reported in allotetraploid cotton genomes, with unknown impacts on fibre traits. Here, we identified gene disruptions in allotetraploid G. hirsutum (18.14%) and G. barbadense (17.38%) through comparison with their presumed diploid progenitors. Relative to conserved genes, these disrupted genes exhibited faster evolution rate, lower expression level and altered gene co-expression networks. Within a module regulating fibre elongation, a hub gene experienced gene disruption in G. hirsutum after polyploidization, with a 2-bp deletion in the coding region of GhNPLA1D introducing early termination of translation. This deletion was observed in all of the 34 G. hirsutum landraces and 36 G. hirsutum cultivars, but not in 96% of 57 G. barbadense accessions. Retrieving the disrupted gene GhNPLA1D using its homoeolog GhNPLA1A achieved longer fibre length in G. hirsutum. Further enzyme activity and lipids analysis confirmed that GhNPLA1A encodes a typical phospholipase A and promotes cotton fibre elongation via elevating intracellular levels of linolenic acid and 34:3 phosphatidylinositol. Our work opens a strategy for identifying disrupted genes and retrieving their functions in ways that can provide valuable resources for accelerating fibre trait enhancement in cotton breeding.

AAQSP increases mapping resolution of stable QTLs through applying NGS-BSA in multiple genetic backgrounds

In this study, we present AAQSP as an extension of existing NGS-BSA applications for identifying stable QTLs at high resolution. GhPAP16 and GhIQD14 fine mapped on chromosome D09 of upland cotton are identified as important candidate genes for lint percentage (LP). Bulked segregant analysis combined with next generation sequencing (NGS-BSA) allows rapid identification of genome sequence differences responsible for phenotypic variation. The NGS-BSA approach applied to crops mainly depends on comparing two bulked DNA samples of individuals from an F₂ population. Since some F₂ individuals still maintain high heterozygosity, heterosis will exert complications in pursuing NGS-BSA in such populations. In addition, the genetic background influences the stability of gene expression in crops, so some QTLs mapped in one segregating population may not be widely applied in crop improvement. The AAQSP (Association Analysis of QTL-seq on Semi-homologous Populations) reported in our study combines the optimized scheme of constructing BSA bulks with NGS-BSA analysis in two (or more) different parental genetic backgrounds for isolating the stable QTLs. With application of AAQSP strategy and construction of a high-density linkage map, we have successfully identified a QTL significantly related to lint percentage (LP) in cultivated upland cotton, followed by map-based cloning to dissect two candidate genes, GhPAP16 and GhIQD14. This study demonstrated that AAQSP can efficiently identify stable QTLs for complex traits of interest, and thus accelerate the genetic improvement of upland cotton and other crop plants.

Functional divergence of GLP genes between G. barbadense and G. hirsutum in response to Verticillium dahliae infection

Germin-like proteins (GLPs) play important roles in plant disease resistance but are rarely reported in cotton. We compared the expression of GLPs in Verticillium dahliae inoculate G. hirsutum (susceptible) and G. barbadense (resistant) and enriched 11 differentially expressed GLPs. 2741 GLP proteins identified from 53 species determined that GLP probably originated from algae and could be classified into 7 clades according to phylogenetic analysis, among which Clade I is likely the most ancient. Cotton GLP (two allopolyploids and two diploids) genes within a shared clade were highly conserved. Intriguingly, clade VII genes were mainly located in gene clusters that derived from the expansion of LTR transposons. Clade VII members expressed mainly in root which is the first battle against Verticillium dahlia and could be induced more intensely in G. barbadense than G. hirsutum. The GLP genes are resistant to Verticillium dahliae, which can be further investigated against Verticillium wilt.

Structure and character analysis of cotton response regulator genes family reveals that GhRR7 responses to draught stress

Background

Cytokinin signal transduction is mediated by a two-component system (TCS). Two-component systems are utilized in plant responses to hormones as well as to biotic and abiotic environmental stimuli. In plants, response regulatory genes (RRs) are one of the main members of the two-component system (TCS).

Method

From the aspects of gene structure, evolution mode, expression type, regulatory network and gene function, the evolution process and role of RR genes in the evolution of the cotton genome were analyzed.

Result

A total of 284 RR genes in four cotton species were identified. Including 1049 orthologous/paralogous gene pairs were identified, most of which were whole genome duplication (WGD). The RR genes promoter elements contain phytohormone responses and abiotic or biotic stress-related cis-elements. Expression analysis showed that RR genes family may be negatively regulate and involved in salt stress and drought stress in plants. Protein regulatory network analysis showed that RR family proteins are involved in regulating the DNA-binding transcription factor activity (COG5641) pathway and HP kinase pathways. VIGS analysis showed that the GhRR7 gene may be in the same regulatory pathway as GhAHP5 and GhPHYB, ultimately negatively regulating cotton drought stress by regulating POD, SOD, CAT, H₂O₂ and other reactive oxygen removal systems.

Conclusion

This study is the first to gain insight into RR gene members in cotton. Our research lays the foundation for discovering the genes related to drought and salt tolerance and creating new cotton germplasm materials for drought and salt tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-022-00394-2.

Parental legacy versus regulatory innovation in salt stress responsiveness of allopolyploid cotton (Gossypium) species

Polyploidy provides an opportunity for evolutionary innovation and species diversification, especially under stressful conditions. In allopolyploids, the conditional dynamics of homoeologous gene expression can be either inherited from ancestral states pre-existing in the parental diploids or novel upon polyploidization, the latter potentially permitting a wider range of phenotypic responses to stresses. To gain insight into regulatory mechanisms underlying the diversity of salt resistance in Gossypium species, we compared global transcriptomic responses to modest salinity stress in two allotetraploid (AD-genome) cotton species, Gossypium hirsutum and G. mustelinum, relative to their model diploid progenitors (A-genome and D-genome). Multivariate and pairwise analyses of salt-responsive changes revealed a profound alteration of gene expression for about one third of the transcriptome. Transcriptional responses and associated functional implications of salt acclimation varied across species, as did species-specific coexpression modules among species and ploidy levels. Salt responsiveness in both allopolyploids was strongly biased toward the D-genome progenitor. A much lower level of transgressive downregulation was observed in the more salt-tolerant G. mustelinum than in the less tolerant G. hirsutum. By disentangling inherited effects from evolved responses, we show that expression biases that are not conditional upon salt stress approximately equally reflect parental legacy and regulatory novelty upon allopolyploidization, whereas stress-responsive biases are predominantly novel, or evolved, in allopolyploids. Overall, our work suggests that allopolyploid cottons acquired a wide range of stress response flexibility relative to their diploid ancestors, most likely mediated by complex suites of duplicated genes and regulatory factors.

Function deficiency of GhOMT1 causes anthocyanidins over-accumulation and diversifies fibre colours in cotton (Gossypium hirsutum)

Naturally coloured cotton (NCC) fibres need little or no dyeing process in textile industry to low-carbon emission and are environment-friendly. Proanthocyanidins (PAs) and their derivatives were considered as the main components causing fibre coloration and made NCCs very popular and healthy, but the monotonous fibre colours greatly limit the wide application of NCCs. Here a G. hirsutum empurpled mutant (HS2) caused by T-DNA insertion is found to enhance the anthocyanidins biosynthesis and accumulate anthocyanidins in the whole plant. HPLC and LC/MS-ESI analysis confirmed the anthocyanidins methylation and peonidin, petunidin and malvidin formation are blocked. The deficiency of GhOMT1 in HS2 was associated with the activation of the anthocyanidin biosynthesis and the altered components of anthocyanidins. The transcripts of key genes in anthocyanidin biosynthesis pathway are significantly up-regulated in HS2, while transcripts of the genes for transport and decoration were at similar levels as in WT. To investigate the potential mechanism of GhOMT1 deficiency in cotton fibre coloration, HS2 mutant was crossed with NCCs. Surprisingly, offsprings of HS2 and NCCs enhanced PAs biosynthesis and increased PAs levels in their fibres from the accumulated anthocyanidins through up-regulated GhANR and GhLAR. As expected, multiple novel lines with improved fibre colours including orange red and navy blue were produced in their generations. Based on this work, a new strategy for breeding diversified NCCs was brought out by promoting PA biosynthesis. This work will help shed light on mechanisms of PA biosynthesis and bring out potential molecular breeding strategy to increase PA levels in NCCs.

Gb_ANR-47 Enhances the Resistance of Gossypium barbadense to Fusarium oxysporum f. sp. vasinfectum (FOV) by Regulating the Content of Proanthocyanidins

Anthocyanidin reductase (ANR) is an important regulator of flavonoid metabolism, and proanthocyanidins, the secondary metabolites of flavonoids, play an important role in the response of plants to pathogenic stress. Therefore, in this study, the expression analysis of the ANR gene family of Gossypium barbadense after inoculation with Fusarium oxysporum f. sp. vasinfectum (FOV) was performed at different time points. It was found that Gb_ANR-47 showed significant differences in the disease-resistant cultivar 06-146 and the susceptible cultivar Xinhai 14, as well as in the highest root expression. It was found that the expression of Gb_ANR-47 in the resistant cultivar was significantly higher than that in the susceptible cultivar by MeJA and SA, and different amounts of methyl jasmonate (MeJA) and salicylic acid (SA) response elements were found in the promoter region of Gb_ANR-47. After silencing GbANR-47 in 06-146 material by VIGS technology, its resistance to FOV decreased significantly. The disease severity index (DSI) was significantly increased, and the anthocyanin content was significantly decreased in silenced plants, compared to controls. Our findings suggest that GbANR-47 is a positive regulator of FOV resistance in Gossypium barbadense. The research results provide an important theoretical basis for in-depth analysis of the molecular mechanism of GbANR-47 and improving the anti-FOV of Gossypium barbadense.

A genome-wide comparative evolutionary analysis of zinc finger-BED transcription factor genes in land plants

Zinc finger (Zf)-BED proteins are a novel superfamily of transcription factors that controls numerous activities in plants including growth, development, and cellular responses to biotic and abiotic stresses. Despite their important roles in gene regulation, little is known about the specific functions of Zf-BEDs in land plants. The current study identified a total of 750 Zf-BED-encoding genes in 35 land plant species including mosses, bryophytes, lycophytes, gymnosperms, and angiosperms. The gene family size was somewhat proportional to genome size. All identified genes were categorized into 22 classes based on their specific domain architectures. Of these, class I (Zf-BED_DUF-domain_Dimer_Tnp_hAT) was the most common in the majority of the land plants. However, some classes were family-specific, while the others were species-specific, demonstrating diversity at different classification levels. In addition, several novel functional domains were also predicated including WRKY and nucleotide-binding site (NBS). Comparative genomics, transcriptomics, and proteomics provided insights into the evolutionary history, duplication, divergence, gene gain and loss, species relationship, expression profiling, and structural diversity of Zf-BEDs in land plants. The comprehensive study of Zf-BEDs in Gossypium sp., (cotton) also demonstrated a clear footprint of polyploidization. Overall, this comprehensive evolutionary study of Zf-BEDs in land plants highlighted significant diversity among plant species.

Genome-wide identification and functional analysis of class III peroxidases in Gossypium hirsutum

Class III peroxidase (PRX) genes play essential roles in various processes, such as auxin catabolism, removal of H₂O₂, crosslinking cell wall components, and response to biotic and abiotic stresses. In this study, we identified 166, 78 and 89 PRX genes from G. hirsutum, G. arboretum and G. raimondii, respectively. These PRX genes were classified into seven subfamilies based on phylogenetic tree analysis and the classification of PRX genes in Arabidopsis. Segmental duplication and purifying selection were the major factors driving the evolution of GhPRXs. GO and KEGG enrichment analysis revealed that GhPRX genes were mainly associated with responding to oxidative stresses, peroxidase activities and phenylpropanoid biosynthesis pathways. Transcriptome data analysis showed that GhPRX genes expression were significantly different in microspore development between the sterility line-JinA and the maintainer line MB177. We confirmed the up-regulation of GhPRX107 and down-regulation of GhPRX128 in the sterile line compared to its maintainer line using qRT-PCR, suggesting their roles in pollen fertility. In addition, silencing GhPRX107 in cotton showed a significant decrease of the reactive oxygen species (ROS) levels of microsporocyte stage anthers compared to control. Overexpressing GhPRX107 in Arabidopsis significantly increased the ROS levels of anthers compared to wild type. In conclusion, we identified GhPRX107 as a determinant of ROS levels in anther. This work sets a foundation for PRX studies in pollen development.

Expression Patterns Divergence of Reciprocal F1 Hybrids Between Gossypium hirsutum and Gossypium barbadense Reveals Overdominance Mediating Interspecific Biomass Heterosis

Hybrid breeding has provided an impetus to the process and achievement of a higher yield and quality of crops. Interspecific hybridization is critical for resolving parental genetic diversity bottleneck problems. The reciprocal interspecific hybrids and their parents (Gossypium hirsutum and Gossypium barbadense) have been applied in this study to elucidate the transcription regulatory mechanism of early biomass heterosis. Phenotypically, the seed biomass, plant height over parent heterosis, leaf area over parent heterosis, and fresh and dry biomass were found to be significantly higher in hybrids than in parents. Analysis of leaf areas revealed that the one-leaf stage exhibits the most significant performance in initial vegetative growth vigor and larger leaves in hybrids, increasing the synthesis of photosynthesis compounds and enhancing photosynthesis compound synthesis. Comparative transcriptome analysis showed that transgressive down-regulation (TDR) is the main gene expression pattern in the hybrids (G. hirsutum × G. barbadense, HB), and it was found that the genes of photosystem I and Adenosine triphosphate (ATP)-binding may promote early growth vigor. Transgressive up-regulation (TUR) is the major primary gene expression pattern in the hybrids (G. barbadense × G. hirsutum, BH), and photosystem II-related genes mediated the performance of early biomass heterosis. The above results demonstrated that overdominance mediates biomass heterosis in interspecific hybrid cotton and the supervisory mechanism divergence of hybrids with different females. Photosynthesis and other metabolic process are jointly involved in controlling early biomass heterosis in interspecific hybrid cotton. The expression pattern data of transcriptome sequencing were supported using the qRT-PCR analysis. Our findings could be useful in theoretical and practical studies of early interspecific biomass heterosis, and the results provide potential resources for the theoretical and applied research on early interspecific biomass heterosis.

Enhanced seedling growth by 3-n-pentadecylphenolethanolamide is mediated by fatty acid amide hydrolases in upland cotton (Gossypium hirsutum L.)

Fatty acid amide hydrolase (FAAH) is a conserved amidase that is known to modulate the levels of endogenous N-acylethanolamines (NAEs) in both plants and animals. The activity of FAAH is enhanced in vitro by synthetic phenoxyacylethanolamides resulting in greater hydrolysis of NAEs. Previously, 3-n-pentadecylphenolethanolamide (PDP-EA) was shown to exert positive effects on the development of Arabidopsis seedlings by enhancing Arabidopsis FAAH (AtFAAH) activity. However, there is little information regarding FAAH activity and the impact of PDP-EA in the development of seedlings of other plant species. Here, we examined the effects of PDP-EA on growth of upland cotton (Gossypium hirsutum L. cv Coker 312) seedlings including two lines of transgenic seedlings overexpressing AtFAAH. Independent transgenic events showed accelerated true-leaf emergence compared with non-transgenic controls. Exogenous applications of PDP-EA led to increases in overall seedling growth in AtFAAH transgenic lines. These enhanced-growth phenotypes coincided with elevated FAAH activities toward NAEs and NAE oxylipins. Conversely, the endogenous contents of NAEs and NAE-oxylipin species, especially linoleoylethanolamide and 9-hydroxy linoleoylethanolamide, were lower in PDP-EA treated seedlings than in controls. Further, transcripts for endogenous cotton FAAH genes were increased following PDP-EA exposure. Collectively, our data corroborate that the enhancement of FAAH enzyme activity by PDP-EA stimulates NAE-hydrolysis and that this results in enhanced growth in seedlings of a perennial crop species, extending the role of NAE metabolism in seedling development beyond the model annual plant species, Arabidopsis thaliana.

Detection of Stable Elite Haplotypes and Potential Candidate Genes of Boll Weight Across Multiple Environments via GWAS in Upland Cotton

Boll weight (BW) is a key determinant of yield component traits in cotton, and understanding the genetic mechanism of BW could contribute to the progress of cotton fiber yield. Although many yield-related quantitative trait loci (QTLs) responsible for BW have been determined, knowledge of the genes controlling cotton yield remains limited. Here, association mapping based on 25,169 single-nucleotide polymorphisms (SNPs) and 2,315 insertions/deletions (InDels) was conducted to identify high-quality QTLs responsible for BW in a global collection of 290 diverse accessions, and BW was measured in nine different environments. A total of 19 significant markers were detected, and 225 candidate genes within a 400 kb region (± 200 kb surrounding each locus) were predicted. Of them, two major QTLs with highly phenotypic variation explanation on chromosomes A08 and D13 were identified among multiple environments. Furthermore, we found that two novel candidate genes (Ghir_A08G009110 and Ghir_D13G023010) were associated with BW and that Ghir_D13G023010 was involved in artificial selection during cotton breeding by population genetic analysis. The transcription level analyses showed that these two genes were significantly differentially expressed between high-BW accession and low-BW accession during the ovule development stage. Thus, these results reveal valuable information for clarifying the genetic basics of the control of BW, which are useful for increasing yield by molecular marker-assisted selection (MAS) breeding in cotton.

Potential Roles of 1-Aminocyclopropane-1-carboxylic Acid Synthase Genes in the Response of Gossypium Species to Abiotic Stress by Genome-Wide Identification and Expression Analysis

Ethylene plays a pivotal role in plant stress resistance and 1-aminocyclopropane-1-carboxylic acid synthase (ACS) is the rate-limiting enzyme in ethylene biosynthesis. Upland cotton (Gossypium hirsutum L.) is the most important natural fiber crop, but the function of ACS in response to abiotic stress has rarely been reported in this plant. We identified 18 GaACS, 18 GrACS, and 35 GhACS genes in Gossypiumarboreum, Gossypium raimondii and Gossypiumhirsutum, respectively, that were classified as types I, II, III, or IV. Collinearity analysis showed that the GhACS genes were expanded from diploid cotton by the whole-genome-duplication. Multiple alignments showed that the C-terminals of the GhACS proteins were conserved, whereas the N-terminals of GhACS10 and GhACS12 were different from the N-terminals of AtACS10 and AtACS12, probably diverging during evolution. Most type II ACS genes were hardly expressed, whereas GhACS10/GhACS12 were expressed in many tissues and in response to abiotic stress; for example, they were highly and hardly expressed at the early stages of cold and heat exposure, respectively. The GhACS genes showed different expression profiles in response to cold, heat, drought, and salt stress by quantitative PCR analysis, which indicate the potential roles of them when encountering the various adverse conditions, and provide insights into GhACS functions in cotton’s adaptation to abiotic stress.

SmGDB: genome database of Salvia miltiorrhiza, an important TCM Plant

BACKGROUND

Salvia miltiorrhiza is an important traditional Chinese medicinal (TCM) plant and a model plant in the genetic study of TCM. A series of omics related to Danshen have been published. Integrating, managing, storing, and sharing data has become an urgent problem to be solved in S. miltiorrhiza genetic studies.

OBJECTIVES

The genome database is the link for the exchange, acquisition, and use of different omics data between data producers and users, maximizing value and utilization of data.

METHODS

The genome database included DSS3 genome and five RNA-Seq data. The back-end performs data search and retrieval through the LAMP (Linux, Apache, MySQL, PHP) framework.

RESULTS

Here, we present SmGDB (S. miltiorrhiza genome database; http://8.140.162.85/ ), which houses the latest version of genome sequence and annotation data for S. miltiorrhiza, combining three unpublished RNA-Seq data from our group and two released RNA-Seq data. We also identified a novel gene cluster including seven CYP71D genes involved in the tanshinone synthesis pathway based on genome sequences and expression data. Besides, SmGDB provides user-friendly web interfaces for querying and browsing gene annotation, structure, location, and expression profiles for concerned genes. Popular bioinformatics tools such as 'BLAST', 'Search', 'Heatmap', 'JBrowse', etc., were also provided in SmGDB.

CONCLUSIONS

SmGDB will provide utility for characterizing the structure of the S. miltiorrhiza genome and better understanding gene functions and biological processes underlying complex secondary metabolism in Danshen.

Mining for allelic gold: finding genetic variation in photosynthetic traits in crops and wild relatives

Improvement of photosynthetic traits in crops to increase yield potential and crop resilience has recently become a major breeding target. Synthetic biology and genetic technologies offer unparalleled opportunities to create new genetics for photosynthetic traits driven by existing fundamental knowledge. However, large 'gene bank' collections of germplasm comprising historical collections of crop species and their relatives offer a wealth of opportunities to find novel allelic variation in the key steps of photosynthesis, to identify new mechanisms and to accelerate genetic progress in crop breeding programmes. Here we explore the available genetic resources in food and fibre crops, strategies to selectively target allelic variation in genes underpinning key photosynthetic processes, and deployment of this variation via gene editing in modern elite material.

Identification and expression analysis of arabinogalactan protein genes in cotton reveal the function of GhAGP15 in Verticillium dahliae resistance

Verticillium wilt, primarily caused by the fungal pathogen Verticillium dahliae, is a serious disease in cotton. Arabinogalactan proteins (AGPs), a class of hydroxyproline-rich glycoproteins, have been widely implicated in plant growth and environmental adaptation. The purpose of this study is to identify and characterize AGP members in cotton plants and explore their roles in responding to environmental stressors. In total, 65 GhAGP members were identified in upland cotton (Gossypium hirsutum), along with 43, 35, and 37 AGP members that were also identified in G. barbadense, G. arboreum, and G. raimondii, respectively. According to gene structure and protein domains analysis, GhAGP genes in upland cotton are highly conserved. Meanwhile, tandem duplication events have occurred frequently throughout cotton's evolutionary history. Expression analysis showed that GhAGP genes were widely expressed during growth and development and in response to abiotic stressors. Many cis-elements related to hormonal responses and environmental stressors were detected in GhAGP promoter regions. GhAGP genes participate in responding to cold, drought, and salt stress, and were sensitive to ET signaling. Furthermore, the expression level of GhAGP15 was elevated during V. dahliae infection and resistance against V. dahliae in upland cotton was significantly weakened by silencing GhAGP15 using a virus-induced gene silencing (VIGS) approach. Our results further suggest that the function of GhAGP15 in V. dahliae resistance might be involved in regulation of the JA, SA, and reactive oxygen species (ROS) pathways. The comprehensive analysis of AGP genes in cotton performed in this study provides a basic framework for further functional research of these genes.

Analysis of a tetraploid cotton line Mac7 transcriptome reveals mechanisms underlying resistance against the whitefly Bemisia tabaci

Whitefly inflicts both direct and indirect losses to cotton crop. Whitefly resistant cotton germplasm is a high priority and considered among the best possible solutions to mitigate this issue. In this study, we evaluated cotton leaf curl disease (CLCuD) resistant cotton line Mac7 under whitefly stress. Furthermore, we utilized the already available transcriptome data of Mac7 concerning whitefly stress to elucidate associated mechanisms and identify functionally important genes in cotton. In transcriptomic data analysis, differentially expressed genes (DEGs) were found involved in complex relay pathways, activated on whitefly exposure. The response implicates signalling through resistance genes (R-genes), MAPK, ROS, VQs or RLKs, transcription factors, which leads to the activation of defence responses including, Ca²⁺messengers, phytohormonal cross-talk, gossypol, flavonoids, PhasiRNA and susceptibility genes (S-genes). The qRT-PCR assay of 10 functionally important genes also showed their involvement in differential responses at 24 and 48 h post whitefly infestation. Briefly, our study helps in understanding the resistant nature of Mac7 under whitefly stress.

Analysis of the PRA1 Genes in Cotton Identifies the Role of GhPRA1.B1-1A in Verticillium dahliae Resistance

Verticillium wilt in cotton (Gossypium hirsutum) is primarily caused by Verticillium dahliae. Previous data suggest that prenylated RAB acceptors (PRAs) play essential roles in environmental plant adaptation, although the potential roles of PRA1 in cotton are unclear. Therefore, in this study, PRA1 family members were identified in G. hirsutum, and their roles in biotic and abiotic stresses were analyzed. Thirty-seven GhPRA1 family members were identified in upland cotton, which were divided into eight groups. Gene structure and domain analyses revealed that the sequences of GhPRA1 members in each group were highly conserved. Many environmental stress-related and hormone-response cis-acting elements were identified in the GhPRA1 promoter regions, indicating that they may respond to biotic and abiotic stresses. Expression analysis revealed that GhPRA1 members were widely expressed in upland cotton. The GhPRA1 genes responded to abiotic stress: drought, cold, salt, and heat stress. GhPRA1.B1-1A expression increased after V. dahliae infection. Furthermore, the functional role of GhPRA1.B1-1A was confirmed by overexpression in Arabidopsis thaliana, which enhanced the resistance to V. dahliae. In contrast, V. dahliae resistance was significantly weakened via virus-induced gene silencing of GhPRA1.B1-1A in upland cotton. Simultaneously, reactive oxygen species accumulation; the H2O2, salicylic acid, and jasmonic acid contents; and callose deposition were significantly decreased in cotton plants with GhPRA1.B1-1A silencing. These findings contribute to a better understanding of the biological roles of GhPRA1 proteins and provide candidate genes for cotton breeders for breeding V. dahliae-resistant cultivars.

Transcriptome Analysis Reveals Differences in Anthocyanin Accumulation in Cotton (Gossypium hirsutum L.) Induced by Red and Blue Light

Many factors, including illumination, affect anthocyanin biosynthesis and accumulation in plants. light quality is the key factor affecting the process of photoinduced anthocyanin biosynthesis and accumulation. We observed that the red color of the Upland cotton accession Huiyuan with the R1 mutation turned to normal green color under light-emitting diodes (LEDs), which inspired us to investigate the effect of red and blue lights on the biosynthesis and accumulation of anthocyanins. We found that both red and blue lights elevated accumulation of anthocyanins. Comparative transcriptomic analyses, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and GSEA, revealed that genes differentially expressed under different light conditions were enriched with the pathways of circadian rhythm, phenylpropanoid biosynthesis, anthocyanin biosynthesis, and flavone and flavonol biosynthesis. Not surprisingly, all the major structural genes related to biosynthesis of anthocyanins, including the key regulatory MYB transcription factor (GhPAP1D) and anthocyanin transporter (GhGSTF12), were induced by red or blue light treatment. However, LARs and MATEs related to biosynthesis of proanthocyanidins were more significantly up-regulated by red light radiation than by blue light radiation. Vice versa, the accumulation of anthocyanins under red light was not as high as that under blue light. In addition, we demonstrated a potential role of GhHY5, a key regulator in plant circadian rhythms, in regulation of anthocyanin accumulation, which could be achieved via interaction with GhPAP1D. Together, these results indicate different effect of red and blue lights on biosynthesis and accumulation of anthocyanins and a potential module including GhHY5 and GhPAP1D in regulation of anthocyanin accumulation in cotton. These results also suggest that the substrates responsible the synthesis of anthocyanins under blue light is diverted to biosynthesis of proanthocyanidin under red light.

Comprehensive identification and expression analysis of CRY gene family in Gossypium

BACKGROUND

The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance.

RESULTS

In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments.

CONCLUSION

We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.