Identification of genes associated with fatty acid biosynthesis based on 214 safflower core germplasm.

BACKGROUND: Safflower (Carthamus tinctorius L.) is an oilseed crop with substantial medicinal and economic value. However, the methods for constructing safflower core germplasm resources are limited, and the molecular mechanisms of lipid biosynthesis in safflower seeds are not well understood. RESULTS: In this study, 11 oil-related quantitative traits and 50 pairs of InDel markers were used to assess the diversity of a collection of 605 safflower germplasms. The original safflower germplasm exhibited rich phenotypic diversity, with high variation for most of the phenotypic traits under investigation. Similarly, high genetic diversity was evaluated in the original germplasm, in which the mean Shannon's information index (I), observed heterozygosity (H0), and expected heterozygosity (He) were 0.553, 0.182, and 0.374, respectively. Four subgroups with strong genetic structures were identified and a core germplasm of 214 cultivars was constructed, which is well represented in the original germplasm. Meanwhile, differential expression analysis of the transcriptomes of high and low linoleic acid safflower varieties at two stages of seed development identified a total of 47 genes associated with lipid biosynthesis. High expression of the genes KAS II and SAD enhanced the synthesis and accumulation of oleic acid, while FAD genes like FAD2 (Chr8G0104100), FAD3, FAD7 and FAD8 promoted the consumption of oleic acid conversion. The coordinated regulation of these multiple genes ensures the high accumulation of oleic acid in safflower seed oil. CONCLUSIONS: Based on these findings, a core germplasm of 214 cultivars was constructed and 47 candidate genes related to unsaturated fatty acid biosynthesis and lipid accumulation were identified. These results not only provide guidance for further studies to elucidate the molecular basis of oil lipid accumulation in safflower seeds, but also contribute to safflower cultivar improvements.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L.) response to repeated drought stresses.

BACKGROUND: Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions. RESULTS: We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes. CONCLUSION: Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L) response to repeated drought stresses

BACKGROUND: Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions. RESULTS: We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes. CONCLUSION: Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

Transcriptome analysis and molecular mechanism of linseed (Linum usitatissimum L.) drought tolerance under repeated drought using single-molecule long-read sequencing.

BACKGROUND: Oil flax (linseed, Linum usitatissimum L.) is one of the most important oil crops., However, the increases in drought resulting from climate change have dramatically reduces linseed yield and quality, but very little is known about how linseed coordinates the expression of drought resistance gene in response to different level of drought stress (DS) on the genome-wide level. RESULTS: To explore the linseed transcriptional response of DS and repeated drought (RD) stress, we determined the drought tolerance of different linseed varieties. Then we performed full-length transcriptome sequencing of drought-resistant variety (Z141) and drought-sensitive variety (NY-17) under DS and RD stress at the seedling stage using single-molecule real-time sequencing and RNA-sequencing. Gene Ontology (GO) and reduce and visualize GO (REVIGO) enrichment analysis showed that upregulated genes of Z141 were enriched in more functional pathways related to plant drought tolerance than those of NY-17 were under DS. In addition, 4436 linseed transcription factors were identified, and 1190 were responsive to stress treatments. Moreover, protein-protein interaction (PPI) network analysis showed that the proline biosynthesis pathway interacts with stress response genes through RAD50 (DNA repair protein 50) interacting protein 1 (RIN-1). Finally, proline biosynthesis and DNA repair structural gene expression patterns were verified by RT- PCR. CONCLUSIONS: The drought tolerance of Z141 may be related to its upregulation of drought tolerance genes under DS. Proline may play an important role in linseed drought tolerance by maintaining cell osmotic and protecting DNA from ROS damage. In summary, this study provides a new perspective to understand the drought adaptability of linseed.

Genomic Comparison and Population Diversity Analysis Provide Insights into the Domestication and Improvement of Flax.

Flax has been cultivated for its oil and fiber for thousands of years. However, it remains unclear how the modifications of agronomic traits occurred on the genetic level during flax cultivation. In this study, we conducted genome-wide variation analyses on multiple accessions of oil-use, fiber-use, landraces, and pale flax to identify the genomic variations during flax cultivation. Our findings indicate that, during flax domestication, genes relevant to flowering, dehiscence, oil production, and plant architecture were preferentially selected. Furthermore, regardless of origins, the improvement of the modern oil-use flax preceded that of the fiber-use flax, although the dual selection on oil-use and fiber-use characteristics might have occurred in the early flax domestication. We also found that the expansion of MYB46/MYB83 genes may have contributed to the unique secondary cell wall biosynthesis in flax and the directional selections on MYB46/MYB83 may have shaped the morphological profile of the current oil-use and fiber-use flax.

Sequencing of Chinese castor lines reveals genetic signatures of selection and yield-associated loci.

Oil produced by castor (Ricinus communis) has broad industrial applications. However, knowledge on the genetic diversity, especially genetic alterations that occurred during domestication and subsequent traits selection, of this oil crop is limited. Here, our population genomics analyses show that the Chinese castors have developed a geographic pattern, classified into the southern-, the middle-, and the northern-China groups. We detect a number of candidate genomic loci that are associated with the selection signals during the geographical differentiation and domestication. Using genome-wide association analysis, we identify candidate genes associated with nine agronomically important traits. One of the candidate genes encoding a glycosyltransferase related to cellulose and lignin biosynthesis is associated with both capsule dehiscence and endocarp thickness. We hypothesize that the abundance of cellulose or lignin in endocarp is an important factor for capsule dehiscence. Our results provide foundation for castor breeding and genetic study.

A survey of transcriptome complexity using PacBio single-molecule real-time analysis combined with Illumina RNA sequencing for a better understanding of ricinoleic acid biosynthesis in Ricinus communis.

BACKGROUND: Ricinus communis is a highly economically valuable oil crop plant from the spurge family, Euphorbiaceae. However, the available reference genomes are incomplete and to date studies on ricinoleic acid biosynthesis at the transcriptional level are limited. RESULTS: In this study, we combined PacBio single-molecule long read isoform and Illumina RNA sequencing to identify the alternative splicing (AS) events, novel isoforms, fusion genes, long non-coding RNAs (lncRNAs) and alternative polyadenylation (APA) sites to unveil the transcriptomic complexity of castor beans and identify critical genes related to ricinoleic acid biosynthesis. Here, we identified 11,285 AS-variants distributed in 21,448 novel genes and detected 520 fusion genes, 320 lncRNAs and 9511 (APA-sites). Furthermore, a total of 6067, 5983 and 4058 differentially expressed genes between developing beans of the R. communis lines 349 and 1115 with extremely different oil content were identified at 7, 14 and 21 days after flowering, respectively. Specifically, 14, 18 and 11 DEGs were annotated encoding key enzymes related to ricinoleic acid biosynthesis reflecting the higher castor oil content of 1115 compared than 349. Quantitative real-time RT-PCR further validated fifteen of these DEGs at three-time points. CONCLUSION: Our results significantly improved the existed gene models of R. communis, and a putative model of key genes was built to show the differences between strains 349 and 1115, illustrating the molecular mechanism of castor oil biosynthesis. A multi-transcriptome database and candidate genes were provided to further improve the level of ricinoleic acid in transgenic crops.

Transcriptomic Analysis for Different Sex Types of Ricinus communis L. during Development from Apical Buds to Inflorescences by Digital Gene Expression Profiling.

The castor plant (Ricinus communis L.) is a versatile industrial oilseed crop with a diversity of sex patterns, its hybrid breeding for improving yield and high purity is still hampered by genetic instability of female and poor knowledge of sex expression mechanisms. To obtain some hints involved in sex expression and provide the basis for further insight into the molecular mechanisms of castor plant sex determination, we performed DGE analysis to investigate differences between the transcriptomes of apices and racemes derived from female (JXBM0705P) and monoecious (JXBM0705M) lines. A total of 18 DGE libraries were constructed from the apices and racemes of a wild monoecious line and its isogenic female derivative at three stages of apex development, in triplicate. Approximately 5.7 million clean tags per library were generated and mapped to the reference castor genome. Transcriptomic analysis showed that identical dynamic changes of gene expression were indicated in monoecious and female apical bud during its development from vegetation to reproduction, with more genes expressed at the raceme formation and infant raceme stages compare to the early leaf bud stage. More than 3000 of differentially expressed genes (DEGs) were detected in Ricinus apices at three developmental stages between two different sex types. A number of DEGs involved in hormone response and biosynthesis, such as auxin response and transport, transcription factors, signal transduction, histone demethylation/methylation, programmed cell death, and pollination, putatively associated with sex expression and reproduction were discovered, and the selected DEGs showed consistent expression between qRT-PCR validation and the DGE patterns. Most of those DEGs were suppressed at the early leaf stage in buds of the mutant, but then activated at the following transition stage (5-7-leaf stage) of buds in the mutant, and ultimately, the number of up-regulated DEGs was equal to that of down-regulation in the small raceme of the mutant. In this study, a large number of DEGs and some suggestions involved in sex expression and reproduction were discovered using DGE analysis, which provides large information and valuable hints for next insights into the molecular mechanism of sex determination. It is useful for other further studies in Ricinus.

The phytoremediation potential of bioenergy crop Ricinus communis for DDTs and cadmium co-contaminated soil.

Cadmium (Cd) and dichlorodiphenyltrichloroethane (DDT) or its metabolite residues are frequently detected in agricultural soils and food, posing a threat to human health. The objective of this study was to compare the ability of 23 genotypes of Ricinus communis in mobilizing and uptake of Cd and DDTs (p,p'-DDT, o,p'-DDT, p,p'-DDD and p,p'-DDE) in the co-contaminated soil. The plant genotypes varied largely in the uptake and accumulation of DDTs and Cd, with mean concentrations of 0.37, 0.43 and 70.51 for DDTs, and 1.22, 2.27 and 37.63 mg kg(-1) dw for Cd in leaf, stem and root, respectively. The total uptake of DDTs and Cd varied from 83.1 to 267.8 and 66.0 to 155.1 µg per pot, respectively. These results indicate that R. communis has great potential for removing DDTs and Cd from contaminated soils attributed to its fast growth, high biomass, strong absorption and accumulation for both DDTs and Cd.