A global survey of the gene network and key genes for oil accumulation in cultivated tetraploid cottons

Pan-genome analysis of GT64 gene family and expression response to Verticillium wilt in cotton

BACKGROUND

The GT64 subfamily, belonging to the glycosyltransferase family, plays a critical function in plant adaptation to stress conditions and the modulation of plant growth, development, and organogenesis processes. However, a comprehensive identification and systematic analysis of GT64 in cotton are still lacking.

RESULTS

This study used bioinformatics techniques to conduct a detailed investigation on the GT64 gene family members of eight cotton species for the first time. A total of 39 GT64 genes were detected, which could be classified into five subfamilies according to the phylogenetic tree. Among them, six genes were found in upland cotton. Furthermore, investigated the precise chromosomal positions of these genes and visually represented their gene structure details. Moreover, forecasted cis-regulatory elements in GhGT64s and ascertained the duplication type of the GT64 in the eight cotton species. Evaluation of the Ka/Ks ratio for similar gene pairs among the eight cotton species provided insights into the selective pressures acting on these homologous genes. Additionally, analyzed the expression profiles of the GT64 gene family. Overexpressing GhGT64_4 in tobacco improved its disease resistance. Subsequently, VIGS experiments conducted in cotton demonstrated reduced disease resistance upon silencing of the GhGT64_4, may indicate its involvement in affecting lignin and jasmonic acid biosynthesis pathways, thus impacting cotton resistance. Weighted Gene Co-expression Network Analysis (WGCNA) revealed an early immune response against Verticillium dahliae in G. barbadense compared to G. hirsutum. Quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR) analysis indicated that some GT64 genes might play a role under various biotic and abiotic stress conditions.

CONCLUSIONS

These discoveries enhance our knowledge of GT64 family members and lay the groundwork for future investigations into the disease resistance mechanisms of this gene in cotton.

Genome-Wide and Expression Pattern Analysis of the HIT4 Gene Family Uncovers the Involvement of GHHIT4_4 in Response to Verticillium Wilt in Gossypium hirsutum

Chromatin remodelers are essential for regulating plant growth, development, and responses to environmental stresses. HIT4 (HEAT-INTOLERANT 4) is a novel stress-induced chromatin remodeling factor that has been less studied in abiotic stress and stress resistance, particularly in cotton. In this study, we conducted a comprehensive analysis of the members of the HIT4 gene family in Gossypium hirsutum using bioinformatics methods, including phylogenetic relationships, gene organization, transcription profiles, phylogenetic connections, selection pressure, and stress response. A total of 18 HIT4 genes were identified in four cotton species, with six HIT4 gene members in upland cotton. Based on the evolutionary relationships shown in the phylogenetic tree, the 18 HIT4 protein sequences were classified into four distinct subgroups. Furthermore, we conducted chromosome mapping to determine the genomic locations of these genes and visually represented the structural characteristics of HIT4 in G. hirsutum. In addition, we predicted the regulatory elements in HIT4 in G. hirsutum and conducted an analysis of repetitive sequences and gene collinearity among HIT4 in four cotton species. Moreover, we calculated the Ka/Ks ratio for homologous genes to assess the selection pressure acting on HIT4. Using RNA-seq, we explored the expression patterns of HIT4 genes in G. hirsutum and Gossypium barbadense. Through weighted gene co-expression network analysis (WGCNA), we found that GHHIT4_4 belonged to the MEblue module, which was mainly enriched in pathways such as DNA replication, phagosome, pentose and glucuronate interconversions, steroid biosynthesis, and starch and sucrose metabolism. This module may regulate the mechanism of upland cotton resistance to Verticillium wilt through DNA replication, phagosome, and various metabolic pathways. In addition, we performed heterologous overexpression of GH_D11G0591 (GHHIT4_4) in tobacco, and the results showed a significant reduction in disease index compared to the wild type, with higher expression levels of disease resistance genes in the transgenic tobacco. After conducting a VIGS (virus-induced gene silencing) experiment in cotton, the results indicated that silencing GHHIT4_4 had a significant impact, the resistance to Verticillium wilt weakened, and the internode length of the plants significantly decreased by 30.7% while the number of true leaves increased by 41.5%. qRT-PCR analysis indicated that GHHIT4_4 mainly enhanced cotton resistance to Verticillium wilt by indirectly regulating the PAL, 4CL, and CHI genes. The subcellular localization results revealed that GHHIT4_4 was predominantly distributed in the mitochondria and nucleus. This study offers preliminary evidence for the involvement of the GHHIT4_4 in cotton resistance to Verticillium wilt and lays the foundation for further research on the disease resistance mechanism of this gene in cotton.

GhL1L1 regulates the contents of unsaturated fatty acids by activating the expression of GhFAD2 genes in cotton

Edible oils with high unsaturated fatty acids, particularly oleic acid, are beneficial to human health. Cotton is one of the top five oil crops in the world, but the mechanism of high-quality oil synthesis and regulatory networks in cotton are largely unclear. Here, we identified Leafy cotyledon1-like 1 (GhL1L1), a NF-YB subfamily gene that is specifically expressed during somatic embryogenesis and seed maturation in cotton. Overexpression of GhL1L1 regulates the contents of unsaturated fatty acids in cotton, especially in the seeds, which is associated with altered expression of the cotton fatty acid biosynthesis-related genes. GhL1L1 synergistically enhanced the expression of GhFAD2-1A by binding to the G-box in its promoter, leading to an increase in the content of linoleic acid. Furthermore, this activation could be enhanced by GhNF-YC2 and GhNF-YA1 by form a transcriptional complex. Collectively, these results contribute to provide new insights into the molecular mechanism of oil biosynthesis in cotton and can facilitate genetic manipulation of cotton varieties with enhanced oil content.

Genomic and co-expression network analyses reveal candidate genes for oil accumulation based on an introgression population in Upland cotton (Gossypium hirsutum)

KEY MESSAGE

Integrated QTL mapping and WGCNA condense the potential gene regulatory network involved in oil accumulation. A glycosyl hydrolases gene (GhHSD1) for oil biosynthesis was confirmed in Arabidopsis, which will provide useful knowledge to understand the functional mechanism of oil biosynthesis in cotton. Cotton is an economical source of edible oil for the food industry. The genetic mechanism that regulates oil biosynthesis in cottonseeds is essential for the genetic enhancement of oil content (OC). To explore the functional genomics of OC, this study utilized an interspecific backcross inbred line population to dissect the quantitative trait locus (QTL) interlinked with OC. In total, nine OC QTLs were identified, four of which were novel, and each QTL explained 3.62-34.73% of the phenotypic variation of OC. The comprehensive transcript profiling of developing cottonseeds revealed 3,646 core genes differentially expressed in both inbred parents. Functional enrichment analysis determined 43 genes were annotated with oil biosynthesis processes. Implementation of weighted gene co-expression network analysis showed that 803 differential genes had a significant correlation with the OC phenotype. Further integrated analysis identified seven important genes located in OC QTLs. Of which, the GhHSD1 gene located in stable QTL qOC-Dt3-1 exhibited the highest functional linkages with the other network genes. Phylogenetic analysis showed significant evolutionary differences in the HSD1 sequences between oilseed- and starch- crops. Furthermore, the overexpression of GhHSD1 in Arabidopsis yielded almost 6.78% higher seed oil. This study not only uncovers important genetic loci for oil accumulation in cottonseed, but also provides a set of new candidate genes that potentially influence the oil biosynthesis pathway in cottonseed.

Hub Genes in Stable QTLs Orchestrate the Accumulation of Cottonseed Oil in Upland Cotton via Catalyzing Key Steps of Lipid-Related Pathways

Upland cotton is the fifth-largest oil crop in the world, with an average supply of nearly 20% of vegetable oil production. Cottonseed oil is also an ideal alternative raw material to be efficiently converted into biodiesel. However, the improvement in kernel oil content (KOC) of cottonseed has not received sufficient attention from researchers for a long time, due to the fact that the main product of cotton planting is fiber. Previous studies have tagged QTLs and identified individual candidate genes that regulate KOC of cottonseed. The regulatory mechanism of oil metabolism and accumulation of cottonseed are still elusive. In the current study, two high-density genetic maps (HDGMs), which were constructed based on a recombinant inbred line (RIL) population consisting of 231 individuals, were used to identify KOC QTLs. A total of forty-three stable QTLs were detected via these two HDGM strategies. Bioinformatic analysis of all the genes harbored in the marker intervals of the stable QTLs revealed that a total of fifty-one genes were involved in the pathways related to lipid biosynthesis. Functional analysis via coexpression network and RNA-seq revealed that the hub genes in the co-expression network that also catalyze the key steps of fatty acid synthesis, lipid metabolism and oil body formation pathways (ACX4, LACS4, KCR1, and SQD1) could jointly orchestrate oil accumulation in cottonseed. This study will strengthen our understanding of oil metabolism and accumulation in cottonseed and contribute to KOC improvement in cottonseed in the future, enhancing the security and stability of worldwide food supply.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1

BACKGROUND

Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear.

RESULTS

In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR).

CONCLUSIONS

Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

Genetic and Morpho-Physiological Differences among Transgenic and No-Transgenic Cotton Cultivars

Three carbon-chain extension genes associated with fatty acid synthesis in upland cotton (Gossypium hirsutum), namely GhKAR, GhHAD, and GhENR, play important roles in oil accumulation in cotton seeds. In the present study, these three genes were cloned and characterized. The expression patterns of GhKAR, GhHAD, and GhENR in the high seed oil content cultivars 10H1014 and 10H1041 differed somewhat compared with those of 10H1007 and 2074B with low seed oil content at different stages of seed development. GhKAR showed all three cultivars showed higher transcript levels than that of 2074B at 10-, 40-, and 45-days post anthesis (DPA). The expression pattern of GhHAD showed a lower transcript level than that of 2074B at both 10 and 30 DPA but a higher transcript level than that of 2074B at 40 DPA. GhENR showed a lower transcript level than that of 2074B at both 15 and 30 DPA. The highest transcript levels of GhKAR and GhENR were detected at 15 DPA in 10H1007, 10H1014, and 10H1041 compared with 2074B. From 5 to 45 DPA cotton seed, the oil content accumulated continuously in the developing seed. Oil accumulation reached a peak between 40 DPA and 45 DPA and slightly decreased in mature seed. In addition, GhKAR and GhENR showed different expression patterns in fiber and ovule development processes, in which they showed high expression levels at 20 DPA during the fiber elongation stage, but their expression level peaked at 15 DPA during ovule development processes. These two genes showed the lowest expression levels at the late seed maturation stage, while GhHAD showed a peak of 10 DPA in fiber development. Compared to 2074B, the oil contents of GhKAR and GhENR overexpression lines increased 1.05~1.08 folds. These results indicated that GhHAD, GhENR, and GhKAR were involved in both seed oil synthesis and fiber elongation with dual biological functions in cotton.

Control of arbuscule development by a transcriptional negative feedback loop in Medicago

Most terrestrial plants establish a symbiosis with arbuscular mycorrhizal fungi (AMF), which provide them with lipids and sugars in exchange for phosphorus and nitrogen. Nutrient exchange must be dynamically controlled to maintain a mutually beneficial relationship between the two symbiotic partners. The WRI5a and its homologues play a conserved role in lipid supply to AMF. Here, we demonstrate that the AP2/ERF transcription factor MtERM1 binds directly to AW-box and AW-box-like cis-elements in the promoters of MtSTR2 and MtSTR, which are required for host lipid efflux and arbuscule development. The EAR domain-containing transcription factor MtERF12 is also directly activated by MtERM1/MtWRI5a to negatively regulate arbuscule development, and the TOPLESS co-repressor is further recruited by MtERF12 through EAR motif to oppose MtERM1/MtWRI5a function, thereby suppressing arbuscule development. We therefore reveal an ERM1/WRI5a-ERF12-TOPLESS negative feedback loop that enables plants to flexibly control nutrient exchange and ensure a mutually beneficial symbiosis.

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.

Molecular Characterization of the Acyl-CoA-Binding Protein Genes Reveals Their Significant Roles in Oil Accumulation and Abiotic Stress Response in Cotton

Members of the acyl-CoA-binding protein (ACBP) gene family play vital roles in diverse processes related to lipid metabolism, growth and development, and environmental response. Plant ACBP genes have been well-studied in a variety of species including Arabidopsis, soybean, rice and maize. However, the identification and functions of ACBP genes in cotton remain to be elucidated. In this study, a total of 11 GaACBP, 12 GrACBP, 20 GbACBP, and 19 GhACBP genes were identified in the genomes of Gossypium arboreum, Gossypium raimondii, Gossypium babardense, and Gossypium hirsutum, respectively, and grouped into four clades. Forty-nine duplicated gene pairs were identified in Gossypium ACBP genes, and almost all of which have undergone purifying selection during the long evolutionary process. In addition, expression analyses showed that most of the GhACBP genes were highly expressed in the developing embryos. Furthermore, GhACBP1 and GhACBP2 were induced by salt and drought stress based on a real-time quantitative PCR (RT-qPCR) assay, indicating that these genes may play an important role in salt- and drought-stress tolerance. This study will provide a basic resource for further functional analysis of the ACBP gene family in cotton.

Identification of candidate genes in cotton associated with specific seed traits and their initial functional characterization in Arabidopsis

Oilseed crops are used to produce vegetable oil to satisfy the requirements of humans and livestock. Cotton (Gossypium spp.) is of great economic value because it is used as both an important textile commodity and a nutrient-rich resource. Cottonseed oil is rich in polyunsaturated fatty acids and does not contain trans fatty acids; hence, it is considered a healthy vegetable oil. However, research on the genetic basis for cottonseed protein content, oil production, and fatty acid composition is lacking. Here, we investigated the protein content, oil content, and fatty acid composition in terms of oleic acid (C18:1) and linoleic acid (C18:2) in mature cottonseeds from 318 Gossypium hirsutum accessions. Moreover, we examined the dynamic change of protein content and lipid composition including palmitic acid (C16:0), stearic acid (C18:0), oleic acid (C18:1), linoleic acid (C18:2), and linolenic acid (C18:3) in developing seeds from 258 accessions at 10 and 20 days post-anthesis. Then, we conducted a genome-wide association study and identified 152 trait-associated loci and 64 candidate genes responsible for protein and oil-related contents in mature cottonseeds and ovules. Finally, six candidate genes were experimentally validated to be involved in the regulation of fatty acid biosynthesis through heterologous expression in Arabidopsis. These results comprise a solid foundation for expanding our understanding of lipid biosynthesis in cotton, which will help breeders manipulate protein and oil contents to make it a fully developed 'fiber, food, and oil crop'.

Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton

Cotton is not only the most important fiber crop but also the fifth most important oilseed crop in the world because of its oil-rich seeds as a byproduct of fiber production. By comparative transcriptome analysis between two germplasms with diverse oil accumulation, we reveal pieces of the gene expression network involved in the process of oil synthesis in cottonseeds. Approximately, 197.16 Gb of raw data from 30 RNA sequencing samples with 3 biological replicates were generated. Comparison of the high-oil and low-oil transcriptomes enabled the identification of 7682 differentially expressed genes (DEGs). Based on gene expression profiles relevant to triacylglycerol (TAG) biosynthesis, we proposed that the Kennedy pathway (diacylglycerol acyltransferase-catalyzed diacylglycerol to TAG) is the main pathway for oil production, rather than the phospholipid diacylglycerol acyltransferase-mediated pathway. Using weighted gene co-expression network analysis, 5312 DEGs were obtained and classified into 14 co-expression modules, including the MEblack module containing 10 genes involved in lipid metabolism. Among the DEGs in the MEblack module, GhCYSD1 was identified as a potential key player in oil biosynthesis. The overexpression of GhCYSD1 in yeast resulted in increased oil content and altered fatty acid composition. This study may not only shed more light on the underlying molecular mechanism of oil accumulation in cottonseed oil, but also provide a set of new gene for potential enhancement of oil content in cottonseeds.

Linkage and association analyses reveal that hub genes in energy-flow and lipid biosynthesis pathways form a cluster in upland cotton

Upland cotton is an important allotetraploid crop that provides both natural fiber for the textile industry and edible vegetable oil for the food or feed industry. To better understand the genetic mechanism that regulates the biosynthesis of storage oil in cottonseed, we identified the genes harbored in the major quantitative trait loci/nucleotides (QTLs/QTNs) of kernel oil content (KOC) in cottonseed via both multiple linkage analyses and genome-wide association studies (GWAS). In ‘CCRI70′ RILs, six stable QTLs were simultaneously identified by linkage analysis of CHIP and SLAF-seq strategies. In ‘0-153′ RILs, eight stable QTLs were detected by consensus linkage analysis integrating multiple strategies. In the natural panel, thirteen and eight loci were associated across multiple environments with two algorithms of GWAS. Within the confidence interval of a major common QTL on chromosome 3, six genes were identified as participating in the interaction network highly correlated with cottonseed KOC. Further observations of gene differential expression showed that four of the genes, LtnD, PGK, LPLAT1, and PAH2, formed hub genes and two of them, FER and RAV1, formed the key genes in the interaction network. Sequence variations in the coding regions of LtnD, FER, PGK, LPLAT1, and PAH2 genes may support their regulatory effects on oil accumulation in mature cottonseed. Taken together, clustering of the hub genes in the lipid biosynthesis interaction network provides new insights to understanding the mechanism of fatty acid biosynthesis and TAG assembly and to further genetic improvement projects for the KOC in cottonseeds.

Genetics, Breeding and Genetic Engineering to Improve Cottonseed Oil and Protein: A Review

Upland cotton (Gossypium hirsutum) is the world's leading fiber crop and one of the most important oilseed crops. Genetic improvement of cotton has primarily focused on fiber yield and quality. However, there is an increased interest and demand for enhanced cottonseed traits, including protein, oil, fatty acids, and amino acids for broad food, feed and biofuel applications. As a byproduct of cotton production, cottonseed is an important source of edible oil in many countries and could also be a vital source of protein for human consumption. The focus of cotton breeding on high yield and better fiber quality has substantially reduced the natural genetic variation available for effective cottonseed quality improvement within Upland cotton. However, genetic variation in cottonseed oil and protein content exists within the genus of Gossypium and cultivated cotton. A plethora of genes and quantitative trait loci (QTLs) (associated with cottonseed oil, fatty acids, protein and amino acids) have been identified, providing important information for genetic improvement of cottonseed quality. Genetic engineering in cotton through RNA interference and insertions of additional genes of other genetic sources, in addition to the more recent development of genome editing technology has achieved considerable progress in altering the relative levels of protein, oil, fatty acid profile, and amino acids composition in cottonseed for enhanced nutritional value and expanded industrial applications. The objective of this review is to summarize and discuss the cottonseed oil biosynthetic pathway and major genes involved, genetic basis of cottonseed oil and protein content, genetic engineering, genome editing through CRISPR/Cas9, and QTLs associated with quantity and quality enhancement of cottonseed oil and protein.