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.

Multi-environment Evaluations Across Ecological Regions Reveal That the Kernel Oil Content of Cottonseed Is Equally Determined by Genotype and Environment.

Cotton is the fifth-largest oil crop in the world. A high kernel oil content (KOC) and high stability are important cottonseed attributes for food security. In this study, the phenotype of KOC and the genotype-by-environment interaction factors were collectively dissected using 250 recombinant inbred lines, their parental cultivars sGK156 and 901-001, and CCRI70 across multi-environments. ANOVA and correlation analysis showed that both genotype and environment contributed significantly to KOC accumulation. Analyses of additive main effect multiplicative interaction and genotype-by-environment interaction biplot models presented the effects of genotype, environment, and genotype by environment on KOC performance and the stability of the experimental materials. Interaction network analysis revealed that meteorological and geographical factors explained 38% of the total KOC variance, with average daily rainfall contributing the largest positive impact and cumulative rainfall having the largest negative impact on KOC accumulation. This study provides insight into KOC accumulation and could direct selection strategies for improved KOC and field management of cottonseed in the future.

Identification and analysis of oil candidate genes reveals the molecular basis of cottonseed oil accumulation in Gossypium hirsutum L.

KEY MESSAGE: Based on the integration of QTL-mapping and regulatory network analyses, five high-confidence stable QTL regions, six candidate genes and two microRNAs that potentially affect the cottonseed oil content were discovered. Cottonseed oil is increasingly becoming a promising target for edible oil with its high content of unsaturated fatty acids. In this study, a recombinant inbred line (RIL) cotton population was constructed to detect quantitative trait loci (QTLs) for the cottonseed oil content. A total of 39 QTLs were detected across eight different environments, of which five QTLs were stable. Forty-three candidate genes potentially involved in carbon metabolism, fatty acid synthesis and triacylglycerol biosynthesis processes were further obtained in the stable QTL regions. Transcriptome analysis showed that nineteen of these candidate genes expressed during  the developing cottonseed ovules and may affect the cottonseed oil content. Besides, transcription factor (TF) and microRNA (miRNA) co-regulatory network analyses based on the nineteen candidate genes suggested that six genes, two core miRNAs (ghr-miR2949b and ghr-miR2949c), and one TF GhHSL1 were considered to be closely associated with the cottonseed oil content. Moreover, four vital genes were validated by quantitative real-time PCR (qRT-PCR). These results provide insights into the oil accumulation mechanism in developing cottonseed ovules through the construction of a detailed oil accumulation model.