Genetic diversity of a Silybum marianum (L.) Gaertn. germplasm collection revealed by DNA Diversity Array Technology (DArTseq).

Silybum marianum (L.) Gaertn. is a multipurpose crop native to the Mediterranean and middle east regions and mainly known for the hepatoprotective properties of fruit-derived silymarin. Despite growing interest in milk thistle as a versatile crop with medicinal value, its potential in agroindustry is hindered by incomplete domestication and limited genomic knowledge, impeding the development of competitive breeding programs. The present study aimed to evaluate genetic diversity in a panel of S. marianum accessions (n = 31), previously characterized for morphological and phytochemical traits, using 5,178 polymorphic DArTseq SNP markers. The genetic structure investigated using both parametric and non-parametric approaches (e.g. PCA, AWclust, Admixture), revealed three distinctive groups reflecting geographical origins. Indeed, Pop1 grouped accessions from Central Europe and UK, Pop3 consisted mainly of accessions of Italian origin, and Pop2 included accessions from different geographical areas. Interestingly, Italian genotypes showed a divergent phenotypic distribution, particularly in fruit oleic and linoleic acid content, compared to the other two groups. Genetic differentiation among the three groups, investigated by computing pairwise fixation index (FST), confirmed a greater differentiation of Pop3 compared to other subpopulations, also based on other diversity indices (e.g. private alleles, heterozygosity). Finally, 22 markers were declared as putatively under natural selection, of which seven significantly affected some important phenotypic traits such as oleic, arachidonic, behenic and linoleic acid content. These findings suggest that these markers, and overall, the seven SNP markers identified within Pop3, could be exploited in specific breeding programs, potentially aimed at diversifying the use of milk thistle. Indeed, incorporating genetic material from Pop3 haplotypes carrying the selected loci into milk thistle breeding populations might be the basis for developing milk thistle lines with higher levels of oleic, arachidonic, and behenic acids, and lower levels of linoleic acid, paving new avenues for enhancing the nutritional and agronomic characteristics of milk thistle.

Chromosome-level genome assembly of milk thistle (Silybum marianum (L.) Gaertn.).

Silybum marianum (L.) Gaertn., commonly known as milk thistle, is a medicinal plant belonging to the Asteraceae family. This plant has been recognized for its medicinal properties for over 2,000 years. However, the genome of this plant remains largely undiscovered, having no reference genome at a chromosomal level. Here, we assembled the chromosome-level genome of S. marianum, allowing for the annotation of 53,552 genes and the identification of transposable elements comprising 58% of the genome. The genome assembly from this study showed 99.1% completeness as determined by BUSCO assessment, while the previous assembly (ASM154182v1) showed 36.7%. Functional annotation of the predicted genes showed 50,329 genes (94% of total genes) with known protein functions in public databases. Comparative genome analysis among Asteraceae plants revealed a striking conservation of collinearity between S. marianum and C. cardunculus. The genomic information generated from this study will be a valuable resource for milk thistle breeding and for use by the larger research community.

Phenylpropanoids in Silybum marianum cultures treated with cyclodextrins coated with magnetic nanoparticles.

The glucose oligosaccharide-derived cyclodextrins (CDs) are used for improving bioactive compound production in plant cell cultures because, in addition to their elicitation activity, CDs promote product removal from cells. However, despite these advantages, the industrial application of CDs is hampered by their high market price. A strategy to overcome this constraint was recently tested, in which reusable CD polymers coated with magnetic Fe3O4 nanoparticles were harnessed in Vitis vinifera cell cultures to produce t-resveratrol (t-R). In this study, we applied hydroxypropyl-ß-CDs (HPCD) and HPCDs coated with magnetic nanoparticles (HPCD-EPI-MN) in methyl jasmonate (MJ)-treated transgenic Silybum marianum cultures ectopically expressing either a stilbene synthase gene (STS) or a chalcone synthase gene (CHS), and compared their effects on the yields of t-R and naringenin (Ng), respectively. HPCD-EPI-MN at 15 g/L stimulated the accumulation of metabolites in the culture medium of the corresponding transgenic cell lines, with up to 4 mg/L of t-R and 3 mg/L of Ng released after 3 days. Similar amounts were produced in cultures treated with HPCD. Concentrations higher than 15 g/L of HPCD-EPI-MN and prolonged incubation periods negatively affected cell growth and viability in both transgenic cell lines. Reutilization of HPCD-EPI-MN was possible in three elicitation cycles (72 h each), after which the polymer retained 25-30% of its initial efficiency, indicating good stability and reusability. Due to their capacity to adsorb metabolites and their recyclability, the application of magnetic CD polymers may reduce the costs of establishing efficient secondary metabolite production systems on a commercial scale. KEY POINTS: • Long-term transgenic S. marianum suspensions stably produce transgene products • t-R and Ng accumulated extracellularly in cultures elicited with HPCD and HPCD-EPI-MN • The recyclability of HPCD-EPI-MN for metabolite production was proven.

Advances in biochemistry and the biotechnological production of taxifolin and its derivatives.

Taxifolin (dihydroquercetin) and its derivatives are medicinally important flavanonols with a wide distribution in plants. These compounds have been isolated from various plants, such as milk thistle, onions, french maritime, and tamarind. In general, they are commercially generated in semisynthetic forms. Taxifolin and related compounds are biosynthesized via the phenylpropanoid pathway, and most of the biosynthetic steps have been functionally characterized. The knowledge gained through the detailed investigation of their biosynthesis has provided the foundation for the reconstruction of biosynthetic pathways. Plant- and microbial-based platforms are utilized for the expression of such pathways for generating taxifolin-related compounds, either by whole-cell biotransformation or through reconfiguration of the genetic circuits. In this review, we summarize recent advances in the biotechnological production of taxifolin and its derivatives.

Effect of gamma rays and colchicine on silymarin production in cell suspension cultures of Silybum marianum: A transcriptomic study of key genes involved in the biosynthetic pathway.

The aim of this study was to investigate secondary metabolite production in Silybum marianum L. cell suspension cultures obtained from seeds treated with gamma rays (200 and 600 Gy) and 0.05% colchicine. The effects of these treatments on callus induction, growth, viability, and silymarin production were studied, along with the changes in the transcriptome and DNA sequence of chalcone synthase (CHS) genes. The effect of gamma radiation (200 and 600 Gy) on silymarin production in S. marianum dry seeds was also studied using HPLC-UV. All three treatments induced high callus biomass production from leaf segments. The viability of the cell suspension cultures was over 90%. The flavonolignan content measured in the extracellular culture medium of the S. marianum cell suspension was highest after treatment with 600 Gy, followed by 0.05% colchicine, and finally, 200 Gy, after a growth period of 12 days. In general, an increased expression of CHS1, CHS2, and CHS3 genes, accompanied by an increase of silymarin content, was observed in response to all the studied treatments, although the effect was greatest on CHS2 expression. Bioinformatics analysis confirmed that the three CHS2 clones exhibited the highest genetic variation, both in relation to each other and to the CHS1 and CHS3 clones. Based on the results, S. marianum plants obtained from seeds previously exposed to 600 and 200 Gy as well as colchicine constitute a renewable resource with the potential to obtain large amounts of silymarin.

Selection and validation of reference genes for RT-qPCR normalization in different tissues of milk thistle (Silybum marianum, Gaert.).

Quantitative reverse transcription PCR is a sensitive technique for evaluating transcriptional profiles in different experimental datasets. To obtain a reliable quantification of the transcripts level, data normalization with stable reference genes is required. Stable reference genes are identified after analysis of their transcripts profile in every new experiment and species of interest. In Silybum marianum, a widely cultivated officinal plant, only few gene expression studies exist, and reference genes for RT-qPCR studies in the diverse plant tissues have never been investigated before. In this work, the expression stability of 10 candidate reference genes was evaluated in leaves, roots, stems and fruits of S. marianum grown under physiological environmental condition. The stability values for each candidate reference gene were calculated by four canonical statistical algorithms GeNorm, NormFinder, Bestkeeper and ΔCt method in different subsets of samples, then they were ranked with RefFinder from the most to the least suitable for normalization. Best combinations of reference genes are finally proposed for different experimental data sets, including all tissues, vegetative, and reproductive tissues separately. Three target genes putatively involved in important biosynthetic pathway leading to key metabolites in the fruits of milk thistle, such as silymarin and fatty acids, were analyzed with the chosen panels of reference genes, in comparison to the ones used in previous papers. To the best of our knowledge, this is the first report on a reliable and systematic identification and validation of the reference genes for RT-qPCR normalization to study gene expression in S. marianum.

Efficient Biosynthesis of (2S)-Eriodictyol from (2S)-Naringenin in Saccharomyces cerevisiae through a Combination of Promoter Adjustment and Directed Evolution.

The compound (2S)-eriodictyol is an important flavonoid that can be derived from (2S)-naringenin through flavonoid 3'-hydroxylase (F3'H) catalyzation. F3'H is a cytochrome P450 enzyme that requires a cytochrome P450 reductase (CPR) to function. However, P450s have limited applications in industrial scale biosynthesis, owing to their low activity. Here, an efficient SmF3'H and a matched SmCPR were identified from Silybum marianum. To improve the efficiency of SmF3'H, we established a high-throughput detection method for (2S)-eriodictyol, in which the promoter combination of SmF3'H and SmCPR were optimized in Saccharomyces cerevisiae. The results revealed that SmF3'H/SmCPR should be expressed by using promoters with similar and strong expression levels. Furthermore, directed evolution was applied to further improve the efficiency of SmF3'H/SmCPR. With the optimized promoter and mutated combinations SmF3'HD285N/SmCPRI453V, the (2S)-eriodictyol titer was improved to 3.3 g/L, the highest titer in currently available reports. These results indicated that S. cerevisiae is an ideal platform for functional expression of flavonoid related P450 enzymes.

Gene Expression Analysis and Metabolite Profiling of Silymarin Biosynthesis during Milk Thistle (Silybum marianum (L.) Gaertn.) Fruit Ripening.

Mature fruits (i.e., achenes) of milk thistle (Silybum marianum (L.) Gaertn., Asteraceae) accumulate high amounts of silymarin (SILM), a complex mixture of bioactive flavonolignans deriving from taxifolin. Their biological activities in relation with human health promotion and disease prevention are well described. However, the conditions of their biosynthesis in planta are still obscure. To fill this gap, fruit development stages were first precisely defined to study the accumulation kinetics of SILM constituents during fruit ripening. The accumulation profiles of the SILM components during fruit maturation were determined using the LC-MS analysis of these defined developmental phases. The kinetics of phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS) and peroxidase (POX) activities suggest in situ biosynthesis of SILM from l-Phenylalanine during fruit maturation rather than a transport of precursors to the achene. In particular, in contrast to laccase activity, POX activity was associated with the accumulation of silymarin, thus indicating a possible preferential involvement of peroxidase(s) in the oxidative coupling step leading to flavonolignans. Reference genes have been identified, selected and validated to allow accurate gene expression profiling of candidate biosynthetic genes (PAL, CAD, CHS, F3H, F3'H and POX) related to SILM accumulation. Gene expression profiles were correlated with SILM accumulation kinetic and preferential location in pericarp during S. marianum fruit maturation, reaching maximum biosynthesis when desiccation occurs, thus reinforcing the hypothesis of an in situ biosynthesis. This observation led us to consider the involvement of abscisic acid (ABA), a key phytohormone in the control of fruit ripening process. ABA accumulation timing and location during milk thistle fruit ripening appeared in line with a potential regulation of the SLIM accumulation. A possible transcriptional regulation of SILM biosynthesis by ABA was supported by the presence of ABA-responsive cis-acting elements in the promoter regions of the SILM biosynthetic genes studied. These results pave the way for a better understanding of the biosynthetic regulation of SILM during the maturation of S. marianum fruit and offer important insights to better control the production of these medicinally important compounds.

A New Strategy to Increase Production of Genoprotective Bioactive Molecules from Cotyledon-Derived Silybum marianum L. Callus.

There is a need to enhance the production of bioactive secondary metabolites and to establish new production systems, e.g., for liver-protective compounds of Silybum marianum seeds. Quantifying and identifying the produced phytochemicals, and examining their protective effects against genotoxic agents, is of great interest. This study established a protocol for the qualitative and quantitative production of hepatoprotective compounds in cotyledon-derived Silybum marianum callus through optimized supplementation of the MS medium with the growth regulators 2,4-D, benzylaminopurine, myoinositol, and asparagine. High-performance liquid chromatography (HPLC) coupled with electrospray ionisation mass spectrometry (ESI-MS) allowed for identification and quantification of the produced compounds. None of the growth medium combinations supported a detectable production of silymarin. Instead, the generated calli accumulated phenolic acids, in particular chlorogenic acid and dicaffeoylquinic acid, as revealed by HPLC and mass spectrometric analysis. 4-Nitro-o-phenylenediamine (NPD) was employed in the AMES-test with Salmonella typhimurium strain TA98 because it is a potent mutagen for this strain. Results revealed that callus extract had a high anti-genotoxic activity with respect to standard silymarin but more evident with respect seed extract. The callus produced chlorogenic acid and dicaffeoylquinic acid, which revealed higher bioactivity than silymarin. Both compounds were not formed or could not be detected in the seeds of Silybum marianum Egyptian ecotype.

Molecular insights of genetic variation in milk thistle (Silybum marianum [L.] Gaertn.) populations collected from southwest Iran.

Milk thistle (Silybum marianum) is among the world's popular medicinal plants. Start Codon Targeted (SCoT) marker system was utilized to investigate the genetic variability of 80 S. marianum genotypes from eight populations in Iran. SCoT marker produced 255 amplicons and 84.03% polymorphism was generated. The SCoT marker system's polymorphism information content value was 0.43. The primers' resolving power values were between 4.18 and 7.84. The percentage of polymorphic bands was between 33.3 and 100%. The Nei's gene diversity (h) was 0.19-1.30 with an average 0.72. The Shannon's index (I) ranged from 0.29 to 1.38 with an average value of 0.83. The average gene flow (0.37) demonstrated a high genetic variation among the studied populations. The variation of 42% was displayed by the molecular variance analysis among the populations while a recorded variation of 58% was made within the populations. Current investigation suggested that SCoT marker system could effectively evaluate milk thistle genotypes genetic diversity.

Spatial organization of silybin biosynthesis in milk thistle [Silybum marianum (L.) Gaertn].

Silymarin is a collection of compounds extracted from the medicinal herb milk thistle, among which silybin is the major flavonolignan. However, the biosynthesis pathway of silybin remains unclear. In this study, biomimetic reactions demonstrated that silybin can be synthesized from coniferyl alcohol and taxifolin by the action of peroxidase. The concentration profiles of silybin and its precursors and RNA-Seq analysis of gene expression revealed that the amount of taxifolin and the activity of peroxidase serve as the limiting factors in silybin biosynthesis. Hierarchical clustering of the expression profile of genes of the flavonoid biosynthesis pathway distinguished flowers from other organs. RNA-Seq revealed five candidates for the peroxidase involved in silybin production, among which APX1 (ascorbate peroxidase 1) showed a distinct peroxidase activity and the capacity to synthesize silybin. The spatial organization of silybin biosynthesis in milk thistle was elucidated, which could help our understanding of the biosynthesis of silybin and other flavonolignans.

Silybin content and overexpression of chalcone synthase genes in Silybum marianum L. plants under abiotic elicitation.

Silymarin, a Silybum marianum seed extract containing a mixture of flavonolignans including silybin, is being used as an antihepatotoxic therapy for liver diseases. In this study, the enhancing effect of gamma irradiation on plant growth parameters of S. marianum under salt stress was investigated. The effect of gamma irradiation, either as a single elicitor or coupled with salinity, on chalcone synthase (CHS) gene expression and silybin A + B yield was also evaluated. The silybin A + B content in S. marianum fruits was estimated by liquid chromatography-mass spectrometry (LC-MS/MS). An increase in silybin content was accompanied by up-regulation of the CHS1, CHS2 and CHS3 genes, which are involved in the silybin biosynthetic pathway. The highest silybin A + B production (0.77 g/100 g plant DW) and transcript levels of the three studied genes (100.2-, 91.9-, and 24.3-fold increase, respectively) were obtained with 100GY gamma irradiation and 4000 ppm salty water. The CHS2 and CHS3 genes were partially sequenced and submitted to the NCBI database under the accession numbers KT252908.1 and KT252909.1, respectively. Developing new approaches to stimulate silybin biosynthetic pathways could be a useful tool to potentiate the use of plants as renewable resources of medicinal compounds.

Gene expression and flavonolignan production in fruits and cell cultures of Silybum marianum.

The hepatoprotectant flavonolignan silymarin (Sm) is synthesized through 4-coumaroyl-CoA, which enters both the flavonoid and the monolignol pathway giving the two immediate precursors taxifolin (Tx) and coniferyl alcohol (CA), respectively. Sm formation occurs via oxidative radicalization of Tx and CA and is accumulated at high levels at final stages of maturation of Silybum marianum fruits. By contrast, Sm production is severely reduced in cell cultures of this species, although suspensions are able to excrete Sm compounds into the medium upon elicitation with methyl jasmonate (MeJA) or cyclodextrins (CD). Knowledge of gene expression is important to understand Sm dynamics and to develop strategies aimed at increasing production by means of cell cultures but, to date, only one gene of the pathway (chalcone synthase, SmCHS) has been cloned. Therefore, to elucidate the relationship between expression of Sm pathway genes and production of these metabolites, four cDNA fragments of genes putatively involved in flavonolignan biosynthesis, chalcone isomerase, flavanone 3-hydroxylase, flavonol 3'-hydroxylase and cinnamyl alcohol dehydrogenase, were isolated from Sm producing S. marianum fruits and their expression, together with that of the SmCHS, were studied both in fruits at different maturation stages and in elicited cell suspensions. Combined results at both transcript expression and metabolite levels at three different stages of fruit maturation revealed that the formation of the flavonoid moiety precedes flavonolignan biosynthesis, being Sm accumulation associated to expression of the monolignol pathway. There was not detectable accumulation of transcripts in cell suspensions, however, elicitation with MeJA or CD notably induced expression of the studied fragments. These results indicate that the five genes expressed during maturation of S. marianum fruits may contribute to observed increases in flavonolignan accumulation upon treatment of cell cultures with elicitors.

Chalcone synthase genes from milk thistle (Silybum marianum): isolation and expression analysis.

Silymarin is a flavonoid compound derived from milk thistle (Silybum marianum) seeds which has several pharmacological applications. Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoids; thereby, the identification of CHS encoding genes in milk thistle plant can be of great importance. In the current research, fragments of CHS genes were amplified using degenerate primers based on the conserved parts of Asteraceae CHS genes, and then cloned and sequenced. Analysis of the resultant nucleotide and deduced amino acid sequences led to the identification of two different members of CHS gene family,SmCHS1 and SmCHS2. Third member, full-length cDNA (SmCHS3) was isolated by rapid amplification of cDNA ends (RACE), whose open reading frame contained 1239 bp including exon 1 (190 bp) and exon 2 (1049 bp), encoding 63 and 349 amino acids, respectively. In silico analysis of SmCHS3 sequence contains all the conserved CHS sites and shares high homology with CHS proteins from other plants.Real-time PCR analysis indicated that SmCHS1 and SmCHS3 had the highest transcript level in petals in the early flowering stage and in the stem of five upper leaves, followed by five upper leaves in the mid-flowering stage which are most probably involved in anthocyanin and silymarin biosynthesis.

[Establishment of culture system of Silybum marianum hairy roots and determination of silybin].

This research uses six Agrobacterium rhizogenes R1601, R15384, R1000, A4, R1025 and R1 to infect silymarin explants to induce hairy roots and silibin. All of the six A. rhizogenes can induce Silybum marianum to generate hairy roots and the A. rhizogene A4 shows comparatively high infection on the plant. This research determines the condition to induce silymarin hairy roots by the factors of infection time, pre-culturing, co-culturing and pH value. The fact that MS liquid medium fits the proliferation of silymarin hairy roots is determined. Through PCR molecular identification, it can be seen that the DNA plasmids in the A. rhizogenes are successfully integrated into the genome of transformed roots. Using liquid chromatography, it is determined that the silibin content in silymarin hairy roots is 2.5 times that in the plant In this research, the silymarin hairy roots culturing system is established, which lays a foundation for the study of culturing silymarin hairy roots and producing silibin.

An integrated proteomic approach to decipher the effect of methyl jasmonate elicitation on the proteome of Silybum marianum L. hairy roots.

Jasmonate and its methyl derivative, methyl jasmonate (MeJA), are naturally occurring compounds that mediate several plant physiological processes in response to pathogen attack, wounding, and ozone. Exogenous application of jasmonates triggers defense responses that resemble those initiated by pathogen infection and also modulates the production of certain secondary metabolites in a variety of plant species. In this study, we treated the hairy root cultures of Silybum marianum L. with 100 µM MeJA and then measured the content of Silymarin (SLM). We observed that the SLM content increased significantly after 48 h of MeJA treatment and remained constant for 120 h. However, MeJA treatment caused a significant growth reduction after 96 h incubation. The activity of lipoxygenase as a key enzyme in the jasmonate biosynthesis pathway and anti-oxidative enzymes; peroxidase and ascorbate peroxidase was also significantly increased after MeJA treatment. To elucidate the global effect of jasmonate on gene expression of S. marianum, we employed high resolution two-dimensional gel electrophoresis coupled with tandem mass spectrometry. Out of 670 reproducibly detected protein spots which were analyzed on each given gel, 32 spots were up- or down regulated upon MeJA treatment. Of them, ten proteins such as ER binding protein, glutamine synthetase, pathogenesis-related protein, caffeoyl CoA O-methyltransferase, and profilin-1 could be identified by mass spectrometry analysis. The possible implications of the identified proteins on physiological outcome of MeJA application in S. marianum hairy root culture will be discussed.

Genetic properties of milk thistle ecotypes from Iran for morphological and flavonolignans characters.

The aim of present study was to investigate the genetic variation within and between 32 milk thistle ecotypes collected from northern (23 accessions) and southern (9 accessions) regions of Iran along with two introduced varieties, CN seeds and Budakalaszi, for morphological and flavonolignans properties. The two collections were assessed at separate field experiments. MANOVA for all the morphological traits showed significant difference between ecotypes. Univariate ANOVA verified these differences for most of the traits in the northern ecotypes (first collection) while for southern ecotypes no significant differences were obtained for the studied traits except seed yield. Among and within ecotypes genotypic coefficient of variation indicated higher level of variation among ecotypes than within ecotypes. In both of the experiments, there was a large genetic variation for silybin and silymarin quality and quantity. Cluster analysis of 34 accessions was performed for morphological traits and silymarin and silybin characteristics, separately. The resulting dendrogram based on silybin and silymarin characteristics revealed that the native accessions such as Dezfoul, Fereydounkenar and Nour, had highest flavonolignans and they were better than the foreign varieties. Also, there was no clear relationship between clustering based on morphological traits and flavonolignan compounds.

Variability pattern and correlation studies in Silybum marianum Gaertn.

Phenotypic and genotypic coefficient of variability, heritability in broad sense and genetic advance were determined investigating the characters of 15 accessions of Silybum marianum. Seed yield/plant and number of capsules/plant had highest estimates of genotypic variation, heritability and genetic advance which suggest that direct selection for these traits is suitable for the improvement of this crop. Number of capsules/plant had a significant positive correlation with the number of branches/plant and leaf length (r=0.3398, 0.7547), whereas seed yield/plant had a positive significant correlation with leaf length, stem diameter, capsule diameter and silymarin content (r=0.6830, 0.3140, 0.3484, 0.2925).