Computational Methods for Mapping, Assembly and Quantification for Coding and Non-coding Transcripts

Splicing inhibition mediated by reduced splicing factors and helicases is associated with the cellular response of lung cancer cells to cisplatin

Lung cancer's mortality is predominantly linked to post-chemotherapy recurrence, driven by the reactivation of dormant cancer cells. Despite the critical role of these reactivated cells in cancer recurrence and metastasis, the molecular mechanisms governing their therapeutic selection remain poorly understood. In this study, we conducted an integrative analysis by combining PacBio single molecule real-time (SMRT) sequencing with short reads Illumina RNA-seq. Our study revealed that cisplatin-induced dormant and reactivated cancer cells exhibited a noteworthy reduction in gene transcripts and alternative splicing events. Particularly, the differential alternative splicing events were found to be overlapping with the differentially expression genes and enriched in genes related to cell cycle and cell division. Utilizing ENCORI database and correlation analysis, we identified key splicing factors, including SRSF7, SRSF3, PRPF8, and HNRNPC, as well as RNA helicase such as EIF4A3, DDX39A, DDX11, and BRIP1, which were associated with the observed reduction in alternative splicing and subsequent decrease in gene expression. Our study demonstrated that lung cancer cells reduce gene transcripts through diminished alternative splicing events mediated by specific splicing factors and RNA helicase in response to the chemotherapeutic stress. These findings provide insights into the molecular mechanisms underlying the therapeutic selection and reactivation of dormant cancer cells. This discovery opens a potential avenue for the development of therapeutic strategies aimed at preventing cancer recurrence following chemotherapy.

Bioinformatics for Prohormone and Neuropeptide Discovery

Neuropeptides and peptide hormones are signaling molecules produced via complex posttranslational modifications of precursor proteins known as prohormones. Neuropeptides activate specific receptors and are associated with the regulation of physiological systems and behaviors. The identification of prohormones-and the neuropeptides created by these prohormones-from genomic assemblies has become essential to support the annotation and use of the rapidly growing number of sequenced genomes. Here we describe a well-validated methodology for identifying the prohormone complement from genomic assemblies that employs widely available public toolsets and databases. The uncovered prohormone sequences can then be screened for putative neuropeptides to enable accurate proteomic discovery and validation.

Exosomes from cisplatin-induced dormant cancer cells facilitate the formation of premetastatic niche in bone marrow through activating glycolysis of BMSCs

Introduction

Lung cancer is the leading cause of cancer-related deaths worldwide. Chemotherapy kills most cancer cells; however, residual cells enter a dormant state. The dormant cancer cells can be reactivated under specific circumstances. The "premetastatic niche" that is suitable for colonization of cancer cells is formed before the arrival of cancer cells. Tumor-derived exosomes are the main mediators of tumorigenesis. We are aiming to elucidate the roles of exosomes from cisplatin-induced dormant lung cancer cells in the formation of premetastatic niches in bone marrow.

Methods

We performed differential proteomics in dormant A549 cell- and A549 cell-derived exosomes. Non-targeted metabolomics and RNA sequencing were performed to explore the molecular and metabolic reprogramming of bone marrow stromal cells (BMSCs). The growth and metastasis of A549 cells in vivo were monitored by bioluminescence imaging.

Results

We found that Insulin-like growth factor 2 (IGF-2) and Insulin-like growth factor binding protein 2 (IGFBP2) were upregulated in dormant A549 cell-derived exosomes. BMSCs that took up exosomes from dormant A549 cells showed enhanced glycolysis and promoted the growth and metastasis of A549 cells possibly through Insulin-like growth factor 1 receptor (IGF-1R)-induced metabolic reprogramming. Inhibition of the production of lactate and IGF-1R signaling can suppress the growth and metastasis of A549 cells from bone marrow.

Discussion

Overall, we demonstrated that BMSCs formed a premetastatic niche upon taking up exosomes from cisplatin-induced dormant lung cancer cells. BMSCs promoted lung cancer cell growth and metastasis through the reverse Warburg effect.

Cell-cell contact-driven EphB1 cis- and trans- signalings regulate cancer stem cells enrichment after chemotherapy

Reactivation of chemotherapy-induced dormant cancer cells is the main cause of relapse and metastasis. The molecular mechanisms underlying remain to be elucidated. In this study, we introduced a cellular model that mimics the process of cisplatin responsiveness in NSCLC patients. We found that during the process of dormancy and reactivation induced by cisplatin, NSCLC cells underwent sequential EMT-MET with enrichment of cancer stem cells. The ATAC-seq combined with motif analysis revealed that OCT4-SOX2-TCF-NANOG motifs were associated with the enrichment of cancer stem cells induced by chemotherapy. Gene expression profiling suggested a dynamic regulatory mechanism during the process of enrichment of cancer stem cells, where Nanog showed upregulation in the dormant state and SOX2 showed upregulation in the reactivated state. Further, we showed that EphB1 and p-EphB1 showed dynamic expression in the process of cancer cell dormancy and reactivation, where the expression profiles of EphB1 and p-EphB1 showed negatively correlated. In the dormant EMT cells which showed disrupted cell-cell contacts, ligand-independent EphB1 promoted entry of lung cancer cells into dormancy through activating p-p38 and downregulating E-cadherin. On the contrary, in the state of MET, in which cell-cell adhesion was recovered, interactions of EphB1 and ligand EphrinB2 in trans promoted the stemness of cancer cells through upregulating Nanog and Sox2. In conclusion, lung cancer stem cells were enriched during the process of cellular response to chemotherapy. EphB1 cis- and trans- signalings function in the dormant and reactivated state of lung cancer cells respectively. It may provide a therapeutic strategy that target the evolution process of cancer cells induced by chemotherapy.

aTAP: automated transcriptome analysis platform for processing RNA-seq data by de novo assembly

RNA-seq is a sequencing technique that uses next-generation sequencing (NGS) to explore and study the entire transcriptome of a biological sample. NGS-based analyses are mostly performed via command-line interfaces, which is an obstacle for molecular biologists and researchers. Therefore, the higher throughputs from NGS can only be accessed with the help of bioinformatics and computer science expertise. As the cost of sequencing is continuously falling, the use of RNA-seq seems certain to increase. To minimize the problems encountered by biologists and researchers in RNA-seq data analysis, we propose an automated platform with a web application that integrates various bioinformatics pipelines. The platform is intended to enable academic users to more easily analyze transcriptome datasets. Our automated Transcriptome Analysis Platform (aTAP) offers comprehensive bioinformatics workflows, including quality control of raw reads, trimming of low-quality reads, de novo transcriptome assembly, transcript expression quantification, differential expression analysis, and transcript annotation. aTAP has a user-friendly graphical interface, allowing researchers to interact with and visualize results in the web browser. This project offers an alternative way to analyze transcriptome data, by integrating efficient and well-known tools, that is simpler and more accessible to research communities. aTAP is freely available to academic users at https://atap.psu.ac.th/.

The effects of sequencing depth on the assembly of coding and noncoding transcripts in the human genome

Investigating the functions and activities of genes requires proper annotation of the transcribed units. However, transcript assembly efforts have produced a surprisingly large variation in the number of transcripts, and especially so for noncoding transcripts. This heterogeneity in assembled transcript sets might be partially explained by sequencing depth. Here, we used real and simulated short-read sequencing data as well as long-read data to systematically investigate the impact of sequencing depths on the accuracy of assembled transcripts. We assembled and analyzed transcripts from 671 human short-read data sets and four long-read data sets. At the first level, there is a positive correlation between the number of reads and the number of recovered transcripts. However, the effect of the sequencing depth varied based on cell or tissue type, the type of read and the nature and expression levels of the transcripts. The detection of coding transcripts saturated rapidly with both short and long-reads, however, there was no sign of early saturation for noncoding transcripts at any sequencing depth. Increasing long-read sequencing depth specifically benefited transcripts containing transposable elements. Finally, we show how single-cell RNA-seq can be guided by transcripts assembled from bulk long-read samples, and demonstrate that noncoding transcripts are expressed at similar levels to coding transcripts but are expressed in fewer cells. This study highlights the impact of sequencing depth on transcript assembly.

Comparative transcriptome analysis of differentially expressed genes and pathways in Procambarus clarkii (Louisiana crawfish) at different acute temperature stress

Procambarus clarkii is an important economic species in China, and exhibit heat and cold tolerance in the main culture regions. To understand the mechanisms, we analyzed the hepatopancreas transcriptome of P. clarkii treated at 10 °C, 25 °C, and 30 °C, then 2092 DEGs and 6929 DEGs were found in 30 °C stress group and 10 °C stress group, respectively. KEGG pathway enrichment results showed that immune pathway is the main stress pathway for 10 °C treatment and metabolic pathway is the main response pathway for 30 °C treatment, which implies low temperature stress induces the damage of the immune system and increases the susceptibility of bacteria while the body response to high temperature stress through metabolic adjustment. In addition, flow cytometry proved that both high and low temperature stress caused different degrees of apoptosis of hemocytes, and dynamic transcription heat map analysis also identified the differential expression of HSPs family genes and apoptosis pathway genes under different heat stresses. This indicates that preventing damaged protein misfolding and accelerating cell apoptosis are necessary mechanisms for P. clarkii to cope with high and low temperature stress. Our research has deepened our understanding of the complex molecular mechanisms of P. clarkii in response to acute temperature stress, and provided a potential strategy for aquatic animals to relieve environmental duress.

Celastrol inhibits the proliferation and migration of MCF-7 cells through the leptin-triggered PI3K/AKT pathway

Leptin is the pivotal modulator in the onset and progression of breast cancer and obesity. Celastrol, which is extracted from the roots of Tripterygium wilfordi plants, exerts various anticancer bioactivities and has recently emerged as a candidate to treat obesity by improving leptin sensitivity. However, the relationship between leptin and celastrol in the treatment of breast cancer is unknown. Here, the growth and migration of MCF-7 cells induced by leptin were tested to demonstrate the antineoplastic activity of celastrol. Transcriptomic analysis and western blotting were conducted to explore the biological roles of leptin in treating breast cancer with celastrol. The present findings showed that celastrol remarkably reversed leptin-triggered cell proliferation and migration in MCF-7 cells. Fifty-two mRNAs with fivefold higher counts and 149 mRNAs with fivefold lower counts were identified in the celastrol-treated MCF-7 cells. According to the GO and KEGG analyses, the effects of celastrol on MCF-7 cells forced lipid metabolism and the endocrine system. Moreover, leptin treatment induced phosphorylation of leptin receptor and PI3K/AKT in MCF-7 cells, whereas pretreatment with celastrol partly abrogated leptin activation. The binding of celastrol to the leptin receptor was also confirmed by molecular docking. The antitumor effect of celastrol is proposed to be mediated by its binding to the leptin receptor and controlled downregulation of the PI3K/AKT pathway.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Noncoding RNAs in apoptosis: identification and function

Apoptosis is a vital cellular process that is critical for the maintenance of homeostasis in health and disease. The derailment of apoptotic mechanisms has severe consequences such as abnormal development, cancer, and neurodegenerative diseases. Thus, there exist complex regulatory mechanisms in eukaryotes to preserve the balance between cell growth and cell death. Initially, protein-coding genes were prioritized in the search for such regulatory macromolecules involved in the regulation of apoptosis. However, recent genome annotations and transcriptomics studies have uncovered a plethora of regulatory noncoding RNAs that have the ability to modulate not only apoptosis but also many other biochemical processes in eukaryotes. In this review article, we will cover a brief summary of apoptosis and detection methods followed by an extensive discussion on microRNAs, circular RNAs, and long noncoding RNAs in apoptosis.

Transposable element sequence fragments incorporated into coding and noncoding transcripts modulate the transcriptome of human pluripotent stem cells

Transposable elements (TEs) occupy nearly 40% of mammalian genomes and, whilst most are fragmentary and no longer capable of transposition, they can nevertheless contribute to cell function. TEs within genes transcribed by RNA polymerase II can be copied as parts of primary transcripts; however, their full contribution to mature transcript sequences remains unresolved. Here, using long and short read (LR and SR) RNA sequencing data, we show that 26% of coding and 65% of noncoding transcripts in human pluripotent stem cells (hPSCs) contain TE-derived sequences. Different TE families are incorporated into RNAs in unique patterns, with consequences to transcript structure and function. The presence of TE sequences within a transcript is correlated with TE-type specific changes in its subcellular distribution, alterations in steady-state levels and half-life, and differential association with RNA Binding Proteins (RBPs). We identify hPSC-specific incorporation of endogenous retroviruses (ERVs) and LINE:L1 into protein-coding mRNAs, which generate TE sequence-derived peptides. Finally, single cell RNA-seq reveals that hPSCs express ERV-containing transcripts, whilst differentiating subpopulations lack ERVs and express SINE and LINE-containing transcripts. Overall, our comprehensive analysis demonstrates that the incorporation of TE sequences into the RNAs of hPSCs is more widespread and has a greater impact than previously appreciated.

SAUTE: sequence assembly using target enrichment

Background

Illumina is the dominant sequencing technology at this time. Short length, short insert size, some systematic biases, and low-level carryover contamination in Illumina reads continue to make assembly of repeated regions a challenging problem. Some applications also require finding multiple well supported variants for assembled regions.

Results

To facilitate assembly of repeat regions and to report multiple well supported variants when a user can provide target sequences to assist the assembly, we propose SAUTE and SAUTE_PROT assemblers. Both assemblers use de Bruijn graph on reads. Targets can be transcripts or proteins for RNA-seq reads and transcripts, proteins, or genomic regions for genomic reads. Target sequences are nucleotide and protein sequences for SAUTE and SAUTE_PROT, respectively.

Conclusions

For RNA-seq, comparisons with Trinity, rnaSPAdes, SPAligner, and SPAdes assembly of reads aligned to target proteins by DIAMOND show that SAUTE_PROT finds more coding sequences that translate to benchmark proteins. Using AMRFinderPlus calls, we find SAUTE has higher sensitivity and precision than SPAdes, plasmidSPAdes, SPAligner, and SPAdes assembly of reads aligned to target regions by HISAT2. It also has better sensitivity than SKESA but worse precision.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04174-9.

Single cells and transposable element heterogeneity in stem cells and development

Recent innovations in single cell sequencing-based technologies are shining a light on the heterogeneity of cellular populations in unprecedented detail. However, several cellular aspects are currently underutilized in single cell studies. One aspect is the expression and activity of transposable elements (TEs). TEs are selfish sequences of DNA that can replicate, and have been wildly successful in colonizing genomes. However, most TEs are mutated, fragmentary and incapable of transposition, yet they are actively bound by multiple transcription factors, host complex patterns of chromatin modifications, and are expressed in mRNAs as part of the transcriptome in both normal and diseased states. The contribution of TEs to development and cellular function remains unclear, and the routine inclusion of TEs in single cell sequencing analyses will potentially lead to insight into stem cells, development and human disease.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

Chromatin and Epigenetic Rearrangements in Embryonic Stem Cell Fate Transitions

A major event in embryonic development is the rearrangement of epigenetic information as the somatic genome is reprogrammed for a new round of organismal development. Epigenetic data are held in chemical modifications on DNA and histones, and there are dramatic and dynamic changes in these marks during embryogenesis. However, the mechanisms behind this intricate process and how it is regulating and responding to embryonic development remain unclear. As embryos develop from totipotency to pluripotency, they pass through several distinct stages that can be captured permanently or transiently in vitro. Pluripotent naïve cells resemble the early epiblast, primed cells resemble the late epiblast, and blastomere-like cells have been isolated, although fully totipotent cells remain elusive. Experiments using these in vitro model systems have led to insights into chromatin changes in embryonic development, which has informed exploration of pre-implantation embryos. Intriguingly, human and mouse cells rely on different signaling and epigenetic pathways, and it remains a mystery why this variation exists. In this review, we will summarize the chromatin rearrangements in early embryonic development, drawing from genomic data from in vitro cell lines, and human and mouse embryos.

Bioinformatics as a Tool for the Structural and Evolutionary Analysis of Proteins

This chapter deals with the topic of bioinformatics, computational, mathematics, and statistics tools applied to biology, essential for the analysis and characterization of biological molecules, in particular proteins, which play an important role in all cellular and evolutionary processes of the organisms. In recent decades, with the next generation sequencing technologies and bioinformatics, it has facilitated the collection and analysis of a large amount of genomic, transcriptomic, proteomic, and metabolomic data from different organisms that have allowed predictions on the regulation of expression, transcription, translation, structure, and mechanisms of action of proteins as well as homology, mutations, and evolutionary processes that generate structural and functional changes over time. Although the information in the databases is greater every day, all bioinformatics tools continue to be constantly modified to improve performance that leads to more accurate predictions regarding protein functionality, which is why bioinformatics research remains a great challenge.

Transcriptome reveals the important role of metabolic imbalances, immune disorders and apoptosis in the treatment of Procambarus clarkii at super high temperature

Temperature is an important environmental factor in the living environment of crustaceans. Changes in temperature can affect their normal growth and metabolism and even cause bacterial disease. Currently, the potential anti-reverse molecular reaction mechanism of crustaceans during high-temperature conditions has not yet been fully understood. Therefore, in this study, we characterised the transcriptome of Procambarus clarkii using RNA sequencing and performed a comparison between super-high-temperature treated samples and controls. After assembly and annotation, 81,097 unigenes with an average length of 069 bp and 358 differentially expressed genes (DEGs) were identified. Among these DEGs, 264 were differentially upregulated and 94 were differentially downregulated. To obtain comprehensive gene function information, we queried seven databases, namely, Nr, Nt, Pfam, KOG, Swiss-Prot, KEGG, and GO to annotate gene functions. Transcriptome analysis revealed that the identified DEGs have significant effects on immune-related pathways, including lysosomal and phagosomal pathways, and that super-high-temperature conditions can cause disease in P. clarkii. Some significantly downregulated genes are involved in oxidative phosphorylation and the PPAR signalling pathway; this suggests a metabolic imbalance in P. clarkia during extreme temperature conditions. In addition, elevated temperature changed the expression patterns of key apoptosis genes XIAP, CASP2, CASP2, CASP8, and CYTC, thereby confirming that high-temperature conditions caused immune disorders, metabolic imbalance, and, finally, triggered apoptosis. Our results provide a useful foundation for understanding the molecular mechanisms underlying the responses of P. clarkii during high-temperature conditions.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Seq-ing answers: Current data integration approaches to uncover mechanisms of transcriptional regulation

Advancements in the field of next generation sequencing lead to the generation of ever-more data, with the challenge often being how to combine and reconcile results from different OMICs studies such as genome, epigenome and transcriptome. Here we provide an overview of the standard processing pipelines for ChIP-seq and RNA-seq as well as common downstream analyses. We describe popular multi-omics data integration approaches used to identify target genes and co-factors, and we discuss how machine learning techniques may predict transcriptional regulators and gene expression.

Transcriptome visualization and data availability at the Saccharomyces Genome Database

The Saccharomyces Genome Database (SGD; www.yeastgenome.org) maintains the official annotation of all genes in the Saccharomyces cerevisiae reference genome and aims to elucidate the function of these genes and their products by integrating manually curated experimental data. Technological advances have allowed researchers to profile RNA expression and identify transcripts at high resolution. These data can be configured in web-based genome browser applications for display to the general public. Accordingly, SGD has incorporated published transcript isoform data in our instance of JBrowse, a genome visualization platform. This resource will help clarify S. cerevisiae biological processes by furthering studies of transcriptional regulation, untranslated regions, genome engineering, and expression quantification in S. cerevisiae.

Computational Oncology in the Multi-Omics Era: State of the Art

Cancer is the quintessential complex disease. As technologies evolve faster each day, we are able to quantify the different layers of biological elements that contribute to the emergence and development of malignancies. In this multi-omics context, the use of integrative approaches is mandatory in order to gain further insights on oncological phenomena, and to move forward toward the precision medicine paradigm. In this review, we will focus on computational oncology as an integrative discipline that incorporates knowledge from the mathematical, physical, and computational fields to further the biomedical understanding of cancer. We will discuss the current roles of computation in oncology in the context of multi-omic technologies, which include: data acquisition and processing; data management in the clinical and research settings; classification, diagnosis, and prognosis; and the development of models in the research setting, including their use for therapeutic target identification. We will discuss the machine learning and network approaches as two of the most promising emerging paradigms, in computational oncology. These approaches provide a foundation on how to integrate different layers of biological description into coherent frameworks that allow advances both in the basic and clinical settings.

Full-length transcriptome sequencing provides insights into the evolution of apocarotenoid biosynthesis in Crocus sativus

Crocus sativus, containing remarkably amounts of crocin, picrocrocin and safranal, is the source of saffron with tremendous medicinal, economic and cultural importance. Here, we present a high-quality full-length transcriptome of the sterile triploid C. sativus, using the PacBio SMRT sequencing technology. This yields 31,755 high-confidence predictions of protein-coding genes, with 50.1% forming paralogous gene pairs. Analysis on distribution of Ks values suggests that the current genome of C. sativus is probably a product resulting from at least two rounds of whole-genome duplication (WGD) events occurred at ~28 and ~114 million years ago (Mya), respectively. We provide evidence demonstrating that the recent β WGD event confers a major impact on family expansion of secondary metabolite genes, possibly leading to an enhanced accumulation of three distinct compounds: crocin, picrocrocin and safranal. Phylogenetic analysis unravels that the founding member (CCD2) of CCD enzymes necessary for the biosynthesis of apocarotenoids in C. sativus might be evolved from the CCD1 family via the β WGD event. Based on the gene expression profiling, CCD2 is found to be expressed at an extremely high level in the stigma. These findings may shed lights on further genomic refinement of the characteristic biosynthesis pathways and promote germplasm utilization for the improvement of saffron quality.