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.

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.

Genome-Wide Characterization and Analysis of CIPK Gene Family in Two Cultivated Allopolyploid Cotton Species: Sequence Variation, Association with Seed Oil Content, and the Role of GhCIPK6.

Calcineurin B-like protein-interacting protein kinases (CIPKs), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we identified 80 and 78 CIPK genes in the Gossypium hirsutum and G. barbadense, respectively. The phylogenetic and gene structure analysis divided the cotton CIPK genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of Gossypium CIPK gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that GhCIPK genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen GhCIPK genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four GhCIPKs were found to be significantly associated with oil content in one of the three field tests. Although most GhCIPK genes were not associated with natural variations in cotton oil content, the overexpression of the GhCIPK6 gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the GhCIPK genes. In brief, these results enhance our understanding of the roles of CIPK genes in oil synthesis and stress responses.

Genome-wide association study of the oil content in upland cotton (Gossypium hirsutum L.) and identification of GhPRXR1, a candidate gene for a stable QTLqOC-Dt5-1.

Cottonseed oil is one of the most important renewable resources for edible oil and biodiesel. To detect QTLs associated with cottonseed oil content (OC) and identify candidate genes that regulate oil biosynthesis, a panel of upland cotton germplasm lines was selected among those previously used to perform GWASs in China. In the present study, 13 QTLs associated with 53 common SNPs on 13 chromosomes were identified in multiple environments based on 15,369 polymorphic SNPs using the Cotton63 KSNP array. Of these, the OC QTL qOC-Dt5-1 delineated by nine SNPs occurred in a confidence interval of 4 SSRs with previously reported OC QTLs. A combined transcriptome and qRT-PCR analysis revealed that a peroxidase gene (GhPRXR1) was predominantly expressed during the middle-late stage (20-35 days post anthesis) of ovule development. The overexpression of GhPRXR1 in yeast significantly increased the OC by 20.01-37.25 %. Suppression of GhPRXR1 gene expression in the virus-induced gene-silenced cotton reduced the OC by 18.11%. Our results contribute to identifying more OC QTLs and verifying a candidate gene that influences cottonseed oil biosynthesis.

A genome-wide analysis of the phospholipid: diacylglycerol acyltransferase gene family in Gossypium.

BACKGROUND: Cotton (Gossypium spp.) is the most important natural fiber crop worldwide, and cottonseed oil is its most important byproduct. Phospholipid: diacylglycerol acyltransferase (PDAT) is important in TAG biosynthesis, as it catalyzes the transfer of a fatty acyl moiety from the sn-2 position of a phospholipid to the sn-3 position of sn-1, 2-diacylglyerol to form triacylglycerol (TAG) and a lysophospholipid. However, little is known about the genes encoding PDATs involved in cottonseed oil biosynthesis. RESULTS: A comprehensive genome-wide analysis of G. hirsutum, G. barbadense, G. arboreum, and G. raimondii herein identified 12, 11, 6 and 6 PDATs, respectively. These genes were divided into 3 subfamilies, and a PDAT-like subfamily was identified in comparison with dicotyledonous Arabidopsis. All GhPDATs contained one or two LCAT domains at the C-terminus, while most GhPDATs contained a preserved single transmembrane region at the N-terminus. A chromosomal distribution analysis showed that the 12 GhPDAT genes in G. hirsutum were distributed in 10 chromosomes. However, none of the GhPDATs was co-localized with quantitative trait loci (QTL) for cottonseed oil content, suggesting that their sequence variations are not genetically associated with the natural variation in cottonseed oil content. Most GhPDATs were expressed during the cottonseed oil accumulation stage. Ectopic expression of GhPDAT1d increased Arabidopsis seed oil content. CONCLUSIONS: Our comprehensive genome-wide analysis of the cotton PDAT gene family provides a foundation for further studies into the use of PDAT genes to increase cottonseed oil content through biotechnology.

A genome-wide analysis of the lysophosphatidate acyltransferase (LPAAT) gene family in cotton: organization, expression, sequence variation, and association with seed oil content and fiber quality.

BACKGROUND: Lysophosphatidic acid acyltransferase (LPAAT) encoded by a multigene family is a rate-limiting enzyme in the Kennedy pathway in higher plants. Cotton is the most important natural fiber crop and one of the most important oilseed crops. However, little is known on genes coding for LPAATs involved in oil biosynthesis with regard to its genome organization, diversity, expression, natural genetic variation, and association with fiber development and oil content in cotton. RESULTS: In this study, a comprehensive genome-wide analysis in four Gossypium species with genome sequences, i.e., tetraploid G. hirsutum- AD1 and G. barbadense- AD2 and its possible ancestral diploids G. raimondii- D5 and G. arboreum- A2, identified 13, 10, 8, and 9 LPAAT genes, respectively, that were divided into four subfamilies. RNA-seq analyses of the LPAAT genes in the widely grown G. hirsutum suggest their differential expression at the transcriptional level in developing cottonseeds and fibers. Although 10 LPAAT genes were co-localised with quantitative trait loci (QTL) for cottonseed oil or protein content within a 25-cM region, only one single strand conformation polymorphic (SSCP) marker developed from a synonymous single nucleotide polymorphism (SNP) of the At-Gh13LPAAT5 gene was significantly correlated with cottonseed oil and protein contents in one of the three field tests. Moreover, transformed yeasts using the At-Gh13LPAAT5 gene with the two sequences for the SNP led to similar results, i.e., a 25-31% increase in palmitic acid and oleic acid, and a 16-29% increase in total triacylglycerol (TAG). CONCLUSIONS: The results in this study demonstrated that the natural variation in the LPAAT genes to improving cottonseed oil content and fiber quality is limited; therefore, traditional cross breeding should not expect much progress in improving cottonseed oil content or fiber quality through a marker-assisted selection for the LPAAT genes. However, enhancing the expression of one of the LPAAT genes such as At-Gh13LPAAT5 can significantly increase the production of total TAG and other fatty acids, providing an incentive for further studies into the use of LPAAT genes to increase cottonseed oil content through biotechnology.

Breeding Potential of Introgression Lines Developed from Interspecific Crossing between Upland Cotton (Gossypium hirsutum) and Gossypium barbadense: Heterosis, Combining Ability and Genetic Effects.

Upland cotton (Gossypium hirstum L.), which produces more than 95% of the world natural cotton fibers, has a narrow genetic base which hinders progress in cotton breeding. Introducing germplasm from exotic sources especially from another cultivated tetraploid G. barbadense L. can broaden the genetic base of Upland cotton. However, the breeding potential of introgression lines (ILs) in Upland cotton with G. barbadense germplasm integration has not been well addressed. This study involved six ILs developed from an interspecific crossing and backcrossing between Upland cotton and G. barbadense and represented one of the first studies to investigate breeding potentials of a set of ILs using a full diallel analysis. High mid-parent heterosis was detected in several hybrids between ILs and a commercial cultivar, which also out-yielded the high-yielding cultivar parent in F1, F2 and F3 generations. A further analysis indicated that general ability (GCA) variance was predominant for all the traits, while specific combining ability (SCA) variance was either non-existent or much lower than GCA. The estimated GCA effects and predicted additive effects for parents in each trait were positively correlated (at P<0.01). Furthermore, GCA and additive effects for each trait were also positively correlated among generations (at P<0.05), suggesting that F2 and F3 generations can be used as a proxy to F1 in analyzing combining abilities and estimating genetic parameters. In addition, differences between reciprocal crosses in F1 and F2 were not significant for yield, yield components and fiber quality traits. But maternal effects appeared to be present for seed oil and protein contents in F3. This study identified introgression lines as good general combiners for yield and fiber quality improvement and hybrids with high heterotic vigor in yield, and therefore provided useful information for further utilization of introgression lines in cotton breeding.