Interacting Environmental Factors Affect Targeted Milk Thistle Metabolomic Profile and other Growth Components: A Review

 

Muhammad Bilal Hafeez

College of Agronomy, Northwest A&F University, Yangling, China

 

METADATA

 

Paper history

Received: 28 February 2025

Revised: 30 June 2025

Accepted: 20 November 2025

Published online: 26 October 2025

 

Corresponding author

Email: bilalhafeez@nwafu.edu.cn

(Muhammad Bilal Hafeez)

 

Keywords

Cancer

Hepatitis

Biotic stress

Abiotic stress

Metabolites

 

Citation

Hafeez MB (2025). Interacting vnvironmental factors affect targeted milk thistle metabolomic profile and other growth components: A review. Innovations in STEAM: Research & Education 3: 25030202. https://doi.org/10.63793/ISRE/0027

ABSTRACT

 

Background: Milk thistle (Silybum marianum L. Gaertn.) is a medicinally important herb of the family Asteraceae. Its achenes contain the active compound silymarin, which has gained significant attention in the pharmaceutical industry for its hepatoprotective properties, including protection against hepatotoxic agents and stimulation of liver regeneration. However, biosynthesis and accumulation of active ingredients are strongly influenced by environmental variability.

Objective: To provide an inclusive overview of the physiological and phenotypic variations in milk thistle under different abiotic and biotic stresses, with a particular focus on silymarin synthesis and accumulation.

Methodology: A literature-based review was conducted, compiling information from available sources regarding germination, growth behavior, and secondary metabolite production of milk thistle under variable environmental conditions.

Results: Milk thistle not only serves as a medicinally valuable plant but also behaves as a noxious weed. Its germination, growth, and metabolite accumulation, especially silymarin, are highly affected by environmental factors. Despite available research, the full potential of this plant under diverse environmental conditions remains underexplored.

Conclusion: Milk thistle can adapt and grow under diverse conditions. However, stresses such as salinity, temperature, and rainfall adversely affect its growth and development, particularly silymarin production. Understanding the physiological responses and secondary metabolite production of milk thistle under different environmental stresses is crucial for optimizing its medicinal use and managing its weed potential in agriculture.

 


INTRODUCTION

 

Cells, organs, tissues, and metabolic functions at different developmental stages respond differently to environmental conditions. Environmental stresses pose serious challenges to agriculture by increasing consumption demands, limiting land availability, and reducing plant-derived medicinal product yields. Abiotic stresses exert considerable influence on the synthesis of secondary metabolites (Jaleel et al. 2007, Zahra et al. 2022). Milk thistle (Silybum marianum L. Gaertn.) is distributed across several countries but is specifically indigenous to Mediterranean regions. It grows at diverse altitudes, ranging from 700 to 1100 m, and thrives in sub-mountainous to coastal areas (Morazzoni and Bombardelli 1995). It can tolerate a wide range of pH but grows best at 5.5 to 7.6 (Andrzejewska and Sadowska 2008). As a dietary supplement, it is ranked among the top ten and is widely used for liver- and bile-related diseases (Kurkin 2003). Its achenes contain 20–35% fatty oil (Ramasamy and Agarwal 2008). Its oil is rich in vitamins (El-Mallah et al. 2003). Medicinally, milk thistle is used to treat gallbladder and various liver diseases (Abenavoli et al. 2010). It also hinders cholesterol biosynthesis, reduces certain cancer risks, and inhibits leukotriene production. Smith et al. (2008) reported that silymarin sales reached approximately 16.6 million USD in 2018 due to the presence of bioactive compounds. Its medicinal importance lies in the active compound silymarin, an isomeric mixture of flavonolignans including silychristin, silydianin, isosilybin, and silybin (Afshar et al. 2014). Silybin, a major component of silymarin, is in high demand due to its anti-carcinogenic properties. Silymarin stabilizes cell membranes, prevents hepatotoxic damage and stimulates liver regeneration (Fraschini et al. 2002).

The production of secondary metabolites is genetically and environmentally regulated, varying across plant families. Such metabolites enable plants to cope with severe environmental stresses and are used as therapeutic agents. Metabolic pathways and active substances are also severely affected by environmental stresses (Bohnert et al. 1995). Their biosynthesis and accumulation depend strongly on soil properties (Selmar and Kleinwachter 2013). Silymarin content in achenes is influenced by both genotypic variation and environmental conditions (Ghavami and Ramin 2008). Interestingly, silybin content was reported to be higher in cultivated ecotypes, whereas isosilybin, silydianin, and silychristin levels were higher in native ecotypes (Radjabian et al. 2008). Milk thistle can adapt and grow under diverse conditions. However, stresses such as salinity, temperature, and rainfall adversely affect its growth and development, particularly silymarin production (Fig. 1). However, no comprehensive review is present on the effect of biotic and abiotic stresses on milk thistle production and metabolites synthesis. Being one of the most important medicinal plants for treating liver diseases in humans, understanding its ecophysiological behavior is crucial for promoting large-scale cultivation. This review applies a nonlinear regression model to describe milk thistle’s responses under different stresses and to highlight the challenges faced in its cultivation and utilization. Such insights can guide future research and support sustainable production

 

Fig. 1: Illustration of abiotic and biotic stresses affecting milk thistle

.

ABIOTIC STRESSES

 

Milk thistle under salinity stress

 

Germination of plants faces a life-threatening challenge in salt marshes and saline desert areas, leading to the mortality of germinating plants. However, different plant species have their own salinity tolerance mechanisms (Brady and Weil 1996). In general, satisfactory achene germination of milk thistle was recorded up to 6 dS/m salinity stress. A 50% reduction in achene germination and seedling emergence was reported at a salinity level of 9 dS/m. Significant reductions in the number of leaves per plant, main capitulum per plant, achene weight per capitulum, achene weight per plant, and 1000 achene weights were observed at 9 dS/m salinity stress. However, at the 15 dS/m salinity level, plants still produced achenes, but the yield was one-third compared to the control. At low salinity (< 9 dS/m), milk thistle shows limited growth and no effect on grain yield compared to control plants, which is why it is categorized as a facultative halophyte (Ghavami and Ramin 2007). Sedghi et al. (2010) recorded a severe reduction of growth attributes of milk thistle seedlings under salinity, including plumule and radicle length, plumule fresh and dry weight, and germination percentage with increasing salinity. Maximum reduction was observed at 10 dS/m. Kashmir et al. (2016) documented that salinity levels up to 100 mM had a non-significant effect on germination and growth-related parameters of milk thistle, but concentrations higher than Xiao-fang et al. (2000) also found that germination percentage decreased with increasing salinity. Similarly, Ghanbari et al. (2013) reported that shoot and root growth were negatively affected under salinity in milk thistle. Moreover, Solouki et al. (2015) also reported that increasing Na+ concentration decreased germination time, germination percentage, seedling number, coefficient of germination time, radicle length, vigor index, seedling length, plumule fresh weight, and radicle fresh weight. Safikhan et al. (2018) reported that salt stress, especially salinity levels of 8 and 12 dS/m, decreased growth as well as other biochemical attributes, including chlorophyll content, carbohydrates, enzymatic activity, and proline concentration. Hydrogen peroxide (H2O2) concentration also increased, indicating stress severity. It was concluded that salt stress, especially under 8 and 12 dS/m, decreased growth characteristics and chlorophyll content, while proline, carbohydrates, enzymatic activity, and H₂O₂ concentration increased in milk thistle leaves (Fig. 2). Little information has been published regarding the correlation of active substances, yield components, and grain yield under salinity stress in milk thistle (Omidbaigi and Nobakht 2001). Maximum oil content was observed up to 6–9 dS/m salinity level; however, further increase in salinity levels gradually decreased oil content. The positive effect of salinity stress on achene silymarin and silybin content has also been reported (Ghavami and Ramin 2008). Similarly, Zahra et al. (2021a; 2021b; 2022) also reported that silymarin content was enhanced under salinity stress, while severely deteriorating all growth and yield parameters.

In crux, salinity stress adversely affects the growth and yield attributes. Moreover, salinity stress increased the production of silymarin and silybin, which are medicinally important phytochemicals present in the achene. However, further research is required to explore its metabolic shifts under a saline environment.

 

Milk thistle and drought stress

 

Worldwide, plant growth and development are severely affected by drought stress, especially in arid and semi-arid regions (Afshar et al. 2014). Deliri et al. (2010) worked on different milk thistle ecotypes and observed that ecotypic differences are highly significant in relation to drought stress. They emphasized that decreases in chlorophyll content, dry weight, root volume, and root tolerance index, along with an increase in electrolyte leakage, are related to drought stress severity. Moreover, Afshar et al. (2015) also noted that silymarin content increases in drought-affected milk thistle achenes. Furthermore, they elaborated that the amount of silybin increased under water stress, which is a more biologically active compound compared to others. Zahir et al. (2014) found enhanced accumulation of total flavonoids and phenolic content in drought-affected milk thistle. Malekzade et al. (2011) proved that milk thistle oil increased under drought stress. A high content of unsaturated fatty acids accumulated under severe drought stress. Ghassemi-Golezani et al. (2017) reported that under water stress, harvest index, 1000 achene weight, achene yield per plant, number of achenes per plant, and plant biomass decreased. Furthermore, oil percentage and yield also decreased; however, flavonoid content increased in water-stressed milk thistle.

Essential oil levels are significantly reduced under acute water stress. Afshar et al. (2016) evaluated that relative water content, stem diameter, leaf dry weight, and leaf area remained unaffected under moderate drought and were only affected under severe drought stress. They observed a decrease of about 19 and 44% in photosynthesis under moderate and severe drought stress, respectively. Moreover, Zahir et al. (2014) elaborated that water deficiency inhibited shoot and root growth; however, total phenolic content, total protein, antioxidant enzymes, and flavonoids increased under drought stress (Table 1). The potential use of drought stress is to enhance the production of active compounds, especially phenolic compounds (Bettaieb et al. 2009). It has been observed that under drought stress, the total flavonoid and phenolic content increase in milk thistle (Zahir et al. 2014). A significant increase in silymarin accumulation and synthesis was observed in milk thistle achenes under drought stress. So, severe and moderate drought stress enhanced silymarin by 4 and 17% respectively, compared to the control. Under drought conditions, silymarin, silychristin, isosilybin, and silybin also increased, but decreased silydianin content (Afshar et al. 2015). A decrease in grain yield of milk thistle was also reported under drought, so the enhanced concentration of silymarin is not economically beneficial according to Afshar et al. (2014). In conclusion, drought stress causes a severe reduction of all the growth and yield-related traits but enhances its medicinally important secondary metabolite production. The production of these metabolites is stress stress-relieving strategy, but their metabolic profile characterization under mild stress may play a plausible role in uplifting its economic benefits.

 

Milk thistle under temperature stress

 

Temperature is an important abiotic factor that influences plant growth and development. Rahman et al. (2016) documented that temperature changes have regulatory effects on plant height, number of flowers per plant, number of achenes per plant, and crop yield per hectare. Milk thistle achene germination, germination percentage, and the number of seedlings at 15 °C were higher compared to 25 or 35°C. Germination percentage was about 95 and 70% under 15 and 35°C, respectively (Ghavami and Ramin 2007). Kashmir et al. (2016) reported that 25°C (optimum temperature) resulted in higher growth and germination rates; however, lower (15°C) and higher temperatures (40°C) resulted in poor germination and growth. Pourreza and Bahrani (2012) reported that temperatures ranging from

Table 1: Changes in silymarin content under different stresses

 

Stress

Change in silymarin contents

References

Salinity stress

Silymarin ↑; Silybin ↑

Ghavami and Ramin (2008)

Drought stress

Silymarin ↑

Afshar et al. (2014); Afshar et al. (2015);

Silymarin + silybin A & B ↑

Shawky (2015)

Silymarin ↓

Afshar (2015)

Oil ↑

Malekzade et al. (2011)

Silymarin ↑

Zahir et al. (2014)

Silymarin ↓

Afshar (2014)

Density

Silymarin ↑

Azizi et al. (2018); Katar et al. (2013)

Heavy metal

Silymarin is not affected

Rio-Celestino et al. (2006)

Herbicides Metribuzin bentazon

Silymarin ↑

Zheljazkov (2006)

Silymarin ↓

Zheljazkov (2006)

Higher population

Not affected —

Omer et al. (1993)

 

 

Fig. 2: Milk thistle necrosis under salinity stress

 

 

21–27°C were effective in enhancing germination percentage. Heidari et al. (2014) used three varieties of milk thistle to confirm cardinal temperatures related to germination response. They concluded that varietal and temperature differences are of prime importance, especially regarding germination rate, reciprocal time to 50% germination, germination uniformity, germination percentage, and time to 5, 10, 50, 90, and 95% germination. They suggested that the optimum temperature for milk thistle growth is 28–29.5°C.

 

Milk thistle under heavy metal stress

 

Milk thistle often faces heavy metal stress due to its cosmopolitan nature. Khatamipour et al. (2011) reported that cadmium toxicity affected germination rate, germination percentage, seedling growth, fresh and dry weight of shoot and root, and shoot and root length of milk thistle. They also concluded that all concentrations of cadmium (Cd) slightly increased the shoot/root ratio and proline content. Moreover, results indicated that roots were more affected by Cd than shoots. Several researchers reported that milk thistle can grow well in contaminated soils with heavy metals such as zinc (Zn), manganese (Mn), copper (Cu), lead (Pb), chromium (Cr), and Cd (Zheljazkov and Nikolov 1996), and even tolerates the radioactive element cesium (Cs). Zheljazkov and Nikolov (1996) reported that Zn accumulated mainly in leaves and stems, while Mn, Cu, Pb, and Cd accumulated in leaves and roots. Achene yield decreased by 16% under heavy metal stress compared to the control. It was noted that the species can accumulate zinc and lead and can also relocate them to the harvestable parts. For this reason, Del Rio-Celestino et al. (2006) suggested that milk thistle is not a hyperaccumulator. However, silymarin content remained unaffected under heavy metal stress (Zheljazkov and Nikolov 1996). According to Ikram et al. (2025), arsenic (As) stress increased silymarin production up to 80%, and they suggested that increasing its production plays a pivotal role in neutralizing stress and initiating tolerance mechanisms.

 

BIOTIC STRESSES

 

Effect of insect attack on milk thistle

 

Milk thistle is susceptible to insect attack. For instance, Goeden (1971) noticed that an assemblage of phytophagous insects fed or reproduced on milk thistle plants, but apparently no deleterious effect was observed on the root, stem, or reproductive parts of milk thistle. Rhinocyllus conicus (weevil) attacks thistle genera Onopordum, Carduus, Cirsium, and Silybum (Fig. 3) (Goeden and Ricker 1974). R. conicus larvae were also found in the achene tissues and achene heads of milk thistle (Coombs et al. 1996). Clarke and Walter (1993) observed that Nezara viridula infects milk thistle in Queensland, Australia. Abdel-Moniem (2002) reported the presence of the achene head weevil (Larinus latus Herbst) on milk thistle. They noted that weevil achene larvae have an injurious effect on the flower head. A single larva can destroy all the achenes of a flower head ranging from 2 to 3 cm in diameter. In Greece and Iran, Aphis fabae cirsiiacanthoidis and Dysaphis lappae cynarae are well-known aphids (Fig. 3) that attack milk thistle plants (Kavallieratos et al. 2007; Rezwani 2008). Khan et al. (2009) observed that caterpillars of Spodoptera sp. damage leaves at the end of flowering. Snails are pests recorded frequently in wet weather conditions. Abdel-Moniem (2002) reported a reduction in achene heads by L. latus. Dodd (1989) pointed out that weevils have restricted oviposition and low-density larvae per capitulum, with little effect on prolonged flowering of milk thistle. Scientists in Israel focused on dense plant occurrence near ant nests and achene dispersal by ants. Ants move the achene into their nest and remove the oily body (elaiosome) to feed their

Fig. 3: Rhinocyllus conicus attack on milk thistle

larvae, which increases milk thistle vigor and germination (Gabay et al. 1994).

 

Pest and disease attack

 

Pest and disease attacks on plants not only affect growth but also yield. Like other plants, milk thistle is also infected by various pests and microbes. Septoria silybi is a fungus that interferes with photosynthesis and causes leaf lesions (Moscow and Lindow 1989). Roche (1991) observed that S. silybi infects milk thistle plants during daylight when there is a high humidity inoculation period, but rare infestation was observed when light was excluded. The reason behind this infestation is related to the requirement for open stomata for pathogen penetration in milk thistle leaves. Berner et al. (2002) observed that the rust fungus Puccinia punctiformis is a pathogen of Canada thistle but often affects milk thistle. El-Elimat et al. (2014) isolated Aspergillus iizukae from the leaves of milk thistle. Souissi et al. (2005) suggested that Microbotryum silybum (a smut fungus) is a naturally occurring pathogen of Silybum marianum (Tamouridou et al. 2018). Moscow and Lindow (1989) observed S. silybi infection in milk thistle plants over several years in central California. Saccardo (1884) and Oudemans (1923) reported that S. silybi is the only pathogen on the sole host of milk thistle. Moscow and Lindow (1989) conducted a detailed experiment on S. silybi-infected milk thistle and reported that it can survive under dry periods when rain and dew are inadequate. A very low inoculum of S. silybi spores is enough to cause considerable infection, leading the leaves to become necrotic. Under high inoculum, severe disease was observed with numerous necrotic leaves that reduced plant growth and eventually killed the plant (Jamali 2015). Puccinia cruchetiana, P. tyrimni, P. mariana and P. laschii also cause infestation in milk thistle (Brandenburger 1985). Kovαčikovα and Kubνnek (1986) noted that milk thistle is severely infected by the Fusarium genus. Cwalina-Ambroziak et al. (2012) reported approximately six species of this genus that infect milk thistle. Šafrαnkovα et al. (2015) observed mildew, Golovinomyces orontii on milk thistle during the vegetation period. Besides, gray mold (Botrytis cinerea) was observed during the rainy season. At the achene ripening stage, vast infestation was observed on stems, leaves, and anthodia. Additionally, they observed the presence of Fusarium and Rhizoctonia sp. on milk thistle roots and Rhizoglyphus sp. infection on roots and root collars. Milk thistle is also a host for cucumber mosaic virus (Souissi et al. 2005) and tomato spotted virus (Chatzivasiliou et al. 2001). Furthermore, Chatzivasiliou et al. (2001) observed that it is also a host for TSWV (tomato spotted wilt virus).

 

Weed attack and milk thistle productivity

 

One of the most important limiting factors in the production of milk thistle is the lack of weed control (Topalov et al. 1983). Zheljazkov et al. (2006) noted 16 species, of which the most abundant were green foxtail (Setaria viridis L. Beauv.), bermudagrass (Cynodon dactylon L. Pers.), and redroot pigweed (Amaranthus retroflexus L.). The most observed perennials were monocotyledonous johnsongrass (Sorghum halepense L. Pers.), Canada thistle (Cirsium arvense L.), motherwort (Leonurus cardiaca L.), and bindweed (Convolvulus arvensis L.). Other weed species included large crabgrass (Digitaria sanguinalis L.), prostrate pigweed (Amaranthus blitoides S. Wats.), velvetleaf (Abutilon theophrasti Medik.), common cocklebur (Xanthium strumarium L.), black nightshade (Solanum nigrum L.), prostrate knotweed (Polygonum aviculare L.), wild buckwheat (P. convolvulus), jimsonweed (Datura stramonium L.), and common lambsquarters (Chenopodium album L.). They also documented that the highest infestation was found in untreated milk thistle plants.

 

Animal attack

 

The achene bank for milk thistle is very limited (Sofer-Arad et al. 2007), and continuous grazing might control this species within a few years (Fig. 4). The density of milk thistle is severely affected by cattle grazing, including rotational and continuous grazing (De Bruijn and Brok 2006). However, Sofer-Arad et al. (2007) observed that cattle grazing may be associated with higher thistle frequency in the mid-eastern rangelands. Spines in milk thistle deter cattle, thus hampering grazing (Danin and Yom-Tov 1990). Grazing milk thistle is toxic for cattle due to lethal and high content of nitrates (Clark County Noxious Weed Program; CCNWP 2015). Campbell et al. (1979) reported that goats can limit milk thistle biomass and reduce achene production. Goats will graze on milk thistle, but less than 1% of achenes pass through their digestive tract (Sindel 1991). Vinograd et al. (2011) reported that in Israel, milk thistle is a dominant weed and cannot be grazed by sheep and goats.

 

 

Fig. 4: Animal attack on milk thistle

Milk thistle and population density

 

Population density and row spacing between plants have a significant effect on growth, yield, and active compounds of milk thistle. Austin et al. (1988) maintained eight plants per pot (18 cm diameter pot) and found the highest shoot yield after 6 weeks of plantation, while a decrease in shoot yield was observed with increased density of plants per pot. Gabucci et al. (2002) noted that higher population density decreased achene yield, number of blooms per plant, and bloom diameter. Belitz and Sams (2007) observed that achene yield decreased when population density increased, showing a negative correlation between yield, mature achene counts, bloom diameter, number of blooms per plant, and population density. Contrarily, Duran Katar et al. (2013) reported that higher population density increased achene yield and silymarin content in milk thistle. They found a higher yield (83.13 kg ha⁻Ή) and silymarin (1.413 kg ha⁻Ή) at a sowing density of 40,000 plants ha⁻Ή. Omidbaigi et al. (2003) concluded that 50 Χ 30 cm is the most suitable density for milk thistle. Recently, Azizi et al. (2018) observed the tallest plants, the highest grain, and biological yield at 8 plants m⁻² density. Moreover, they noted that population density had no impact on silymarin concentration. However, Omer et al. (1993) noted that narrow row spacing of approximately 25 cm increased achene yield but decreased flavonolignan and oil content compared with 50 cm spacing in milk thistle. They also found that row spacing greater than 25 cm significantly increased silymarin, isosilybin, silychristin, and silybin concentrations.

 

Anthropogenic activities

 

Milk thistle is a noxious weed that is harmful to economic and environmental resources; therefore, plants are targeted for eradication. It has pappi-bearing achenes that are easily pollinated even before harvest, thus emerging as a weed for the next crops. In North America, it is classified as a noxious weed in Washington (category A), Oregon (category B), and Texas (category S2) (Plant Protection and Quarantine 2002), but no case has been reported from Canada (USDA-ARS 2005). A limiting factor for milk thistle production is weed control. Milk thistle is very sensitive to herbicides used for other crops (Topalov et al. 1983). Parsons (1973) noted that it is very easy to eradicate milk thistle plants with several herbicides; however, large flowering and rosette plants are difficult to kill. Shimi et al. (2006) observed that clopyralid (0.24 kg ha⁻Ή) can control 94% of milk thistle growth. In cereals, 2,4-D ester and MCPA, 2,4-D amine can be used to control milk thistle (Department of Primary Industries, Water and Environment, 2008). Zand et al. (2007) noted that bromoxynil plus MCPA at 560 g ha⁻Ή, metsulfuron plus sulfosulfuron at 36 g ha⁻Ή, and chlorsulfuron at 10.5 g ha⁻Ή suppressed milk thistle achene production. So, these herbicides also affect the milk thistle production.

 

CONCLUSION

 

Milk thistle (Silybum marianum L. Gaertn.) is a valuable medicinal plant recognized for its hepatoprotective, antioxidant, and anti-inflammatory properties, primarily attributed to its bioactive compound silymarin. Despite its therapeutic potential, milk thistle faces several abiotic and biotic stresses, including drought, salinity, heavy metals, insect pests, diseases, weeds, and grazing pressure, which significantly affect its growth, yield, and active constituents. Research findings indicate that appropriate agronomic practices, such as optimized population density, nutrient management, and protective measures against pests and weeds, can improve its productivity and phytochemical composition. Furthermore, its classification as both a medicinal resource and a noxious weed highlights the dual challenges in its management. Overall, milk thistle represents a promising plant species with significant pharmaceutical and ecological importance, but sustainable cultivation strategies are essential to maximize its benefits while minimizing its invasive potential.

 

DATA AVAILABILITY

 

Not applicable to this paper

 

ETHICS APPROVAL

 

Not applicable to this paper.

 

FUNDING SOURCE

 

Not applicable to this paper.

 

REFERENCES

 

Abdel-Moniem ASH (2002) The seed-head weevil Larinus latus Herbst (Coleoptera: Curculionidae) as a new record in Egypt on the milk thistle Silybum marianum (L.) (Asteraceae: Compositae). Archives of Phytopathology and Plant Protection 35: 157–160. https://doi.org/10.1080/03235400214210.

Abenavoli L, Capasso R, Milic N, Capasso F (2010) Milk thistle in liver diseases: past, present, future. Phytotherapy Research 24: 1423–1432. https://doi.org/10.1002/ptr.3207.

Abenavoli L, Izzo AA, Milić N, Cicala C, Santini A, Capasso R (2018) Milk thistle (Silybum marianum): A concise overview on its chemistry, pharmacological, and nutraceutical uses in liver diseases. Phytotherapy Research 32: 2202–2213. https://doi.org/10.1002/ptr.6171.

Afshar RK, Chaichi MR, Ansari Jovini M, Jahanzad E, Hashemi M (2015) Accumulation of silymarin in milk thistle seeds under drought stress. Planta 242: 539–543. https://doi.org/10.1007/s00425-015-2265-9.

Afshar RK, Chaichi MR, Asareh MH, Hashemi M, Liaghat A (2014) Interactive effect of deficit irrigation and soil organic amendments on seed yield and flavonolignan production of milk thistle (Silybum marianum L. Gaertn.). Industrial Crops and Products 58: 166–172. https://doi.org/10.1016/j.indcrop.2014.03.043.

Afshar RK, Hashemi M, DaCosta M, Spargo J, Sadeghpour A (2016) Biochar application and drought stress effects on physiological characteristics of Silybum marianum. Communications in Soil Science and Plant Analysis 47: 743–752. https://doi.org/10.1080/00103624.2016.1146752.

Andrzejewska J, Sadowska K (2008) Effect of cultivation conditions on the variability and interrelation of yield and raw material quality in milk thistle (Silybum marianum L. Gaertn.). Acta Scientiarum Polonorum Agricultura 7: 3–11.

Austin MP, Fresco LFM, Nicholls AO, Groves RH, Kaye PE (1988) Competition and relative yield: Estimation and interpretation at different densities and under various nutrient concentrations using Silybum marianum and Cirsium vulgare. Journal of Ecology 76: 157–171. https://doi.org/10.2307/2260460.

Azizi K, Nazari Alam J, Faizian M, Heydari R (2018) Effect of climate and density on different populations of milk thistle (Silybum marianum L.). Iranian Journal of Field Crop Science 49: 65–75. https://doi.org/10.22059/ijfcs.2017.226910.654273.

Belitz AR, Sams CE (2007) Effect of population density on growth, yield, and flavonolignan content in milk thistle (Silybum marianum). Acta Horticulturae 756: 251–258. https://doi.org/10.17660/ActaHortic.2007.756.26.

Berner DK, Paxson LK, Bruckart WL, Luster DG, McMahon MB, Michael JL (2002) First report of Silybum marianum as a host of Puccinia punctiformis. Plant Disease 86: 1271. https://doi.org/10.1094/PDIS.2002.86.11.1271B.

Bettaieb I, Zakhama N, Wannes WA, Kchouk ME, Marzouk B (2009) Water deficit effects on Salvia officinalis fatty acids and essential oils composition. Scientia Horticulturae 120: 271–275. https://doi.org/10.1016/j.scienta.2008.10.016.

Bohnert HJ, Nelson DE, Jensen RG (1995) Adaptations to environmental stresses. Plant Cell 7: 1099–1111. https://doi.org/10.2307/3870060.

Brady NC, Weil RR (1996) The Nature and Properties of Soils, 11th edn. Upper Saddle River, NJ: Prentice Hall, Inc.

Brandenburger W (1985) Parasitische Pilze an Gefδsspflanzen in Europa. Stuttgart: Spektrum Akademischer Verlag.

Campbell MH, Holst PJ, Auld BA, Medd RW (1979) Control of three pasture weeds using goats. In: Proceedings of the 7th Asian-Pacific Weed Science Society Conference, Sydney, Australia, 26–30 November 1979. Asian-Pacific Weed Science Society, Sydney, 201–205. Available at: https://catalogue.nla.gov.au/catalog/2937787.

Chatzivasiliou EK, Bourakas I, Drossos E, Eleftherohorinos I, Jenser G, Peter D, Katis (2001) Weeds in greenhouses and tobacco fields are differentially infected by tomato spotted wilt virus and infested by its vector species. Plant Disease 85: 40–46. https://doi.org/10.1094/PDIS.2001.85.1.40.

Clark County Noxious Weed Program (2015) Milk thistle (Silybum marianum). Clark County Noxious Weed Program. Available at: https://clark.wa.gov/sites/default/files/dept/files/publicworks/Vegetation/Milk%20Thistle%20Control%20Sheet.pdf.

Clarke AR, Walter GH (1993) Variegated thistle (Silybum marianum (L.)): a non-crop host plant of Nezara viridula (Hemiptera: Pentatomidae) in Southeastern Queensland. Journal of the Australian Entomological Society 32: 81–83. https://doi.org/10.1111/j.1440-6055.1993.tb00549.x.

Coombs E, Turner C, Rees N (1995) Milk thistle Silybum marianum. In: Rees N, Quimby P, Piper G, Coombs E, Turner C, Spencer N, Knutson L (eds) Biological Control of Weeds in the West. Western Society of Weed Science in cooperation with USDA-ARS, Montana Department of Agriculture and Montana State University.

Cwalina-Ambroziak B, Wierzbowska J, Damszel M, Bowszys T (2012) The effect of mineral fertilization on achenes yield and fungal communities isolated from the stems of milk thistle Silybum marianum. Acta Scientiarum Polonorum Hortorum Cultus 11: 157–168.

Danin Y, Yom-Tov Y (1990) Ant nests as primary habitats of Silybum marianum (Compositae). Plant Systematics and Evolution 169: 209–217. https://doi.org/10.1007/BF00937675.

De Bruijn SL, Bork EW (2006) Biological control of Canada thistle in temperate pastures using high-density rotational cattle grazing. Biological Control 36: 305–315. https://doi.org/10.1016/j.biocontrol.2005.10.007.

Del Rio-Celestino MD, Font R, Moreno-Rojas R, De Haro-Bailon A (2006) Uptake of lead and zinc by wild plants growing on contaminated soils. Industrial Crops and Products 24: 230–237. https://doi.org/10.1016/j.indcrop.2006.06.013.

Deliri R, Shokrpour M, Asghari A, Esfandiari E, Sharifi RS (2010) Assessment of milk thistle ecotypes for drought resistance in a hydroponic system. Journal of Science and Technology of Greenhouse Culture 1: 9–18.

Department of Primary Industries, Water and Environment (2008) Variegated Thistle (Silybum marianum L.) Control Guide. Hobart: Department of Primary Industries, Water and Environment.

Dodd J (1989) Phenology and seed production of variegated thistle, Silybum marianum (L.) Gaertn., in Australia in relation to mechanical and biological control. Weed Research 29: 255–263. https://doi.org/10.1111/j.1365-3180.1989.tb00910.x.

Duran Katar K, Arslan Y, Subasi I (2013) Effect of different plant densities on growth and yield of milk thistle [Silybum marianum (L.) Gaertn.] grown under ecological conditions of Ankara, Turkey. Research on Crops 14: 304–310.

El-Elimat T, Raja HA, Graf TN, Faeth SH, Cech NB, Oberlies NH (2014) Flavonolignans from Aspergillus iizukae, a fungal endophyte of milk thistle (Silybum marianum). Journal of Natural Products 77(2): 193–199. https://doi.org/10.1021/np400955q.

El-Mallah MH, El-Shami SM, Hassanein MM (2003) Detailed studies on some lipids of Silybum marianum seed oil. Asian Journal of Weed Technology 54: 397–402. https://doi.org/10.3989/gya.2003.v54.i4.227.

Fraschini F, Demartini G, Esposti D (2002) Pharmacology of silymarin. Clinical Drug Investigation 22: 51–65. https://doi.org/10.2165/00044011-200222010-00007.

Gabay R, Plitmann U, Danin A (1994) Factors affecting the dominance of Silybum marianum L. (Asteraceae) in its specific habitats. Flora 189: 201–206. https://doi.org/10.1016/S0367-2530(17)30594-7.

Gabucci L, Curioni A, Garcia M, Urrutia ME (2002) Seeds production in milk thistle crop. Acta Horticulturae 569: 121–128. https://doi.org/10.17660/ActaHortic.2002.569.20.

Ghanbari F, Zarie B, Pour-Aboughadarch A, Kordai S (2013) Investigation of 5-aminolevolinic acid (ALA) effects on seed germination and seedling growth of Silybum marianum under salinity stress. International Journal of Biosciences 3(9): 95–101. https://doi.org/10.12692/ijb/3.9.95-101.

Ghassemi-Golezani K, Ghassemi S, Yaghoubian I (2017) Improving oil and flavonoid contents of milk thistle under water stress by salicylic acid. Advances in Horticultural Science 31: 19–23. https://doi.org/10.13128/ahs-20721.

Ghavami N, Ramin AA (2008). Grain yield and active substanc­es of milk thistle as affected by soil salinity. Communications in Soil Science and Plant Analysis 39: 2608–2618. https://doi.org/10.1080/00103620802358672

Ghavami, N., & Ramin, A. A. (2007). Salinity and temperature effects on seed germination of milk thistle. Communications in Soil Science and Plant Analysis 38: 2681–2691. https://doi.org/10.1080/00103620701662901.

Goeden RD (1971) The phytophagous insect fauna of milk thistle in Southern California. Journal of Economic Entomology 64: 1101–1104. https://doi.org/10.1093/jee/64.5.1101.

Goeden RD, Ricker DW (1974) Imported seed weevils attack Italian and milk thistles in Southern California. California Agriculture 28: 8–9.

Heidari Z, Kamkar B, Sinaki JM (2014) Determination of cardinal temperatures of milk thistle (Silybum marianum L.) germination. Advances in Plants and Agriculture Research 1: 172–178. https://doi.org/10.15406/apar.2014.01.00027.

Ikram M, Pirzada AS, Zulfat M, Khan A and Saqib MN (2025) Silybum marianum (milk thistle): as antioxidant, alleviating pollution-induced oxidative stress in the presence of heavy metals contamination. Global Research Journal of Natural Science and Technology 3: 25–33. https://doi.org/10.53762/grjnst.03.03.02.

Jaleel CA, Manivannan P, Sankar B (2007) Water deficit stress mitigation by calcium chloride in Catharanthus roseus: Effects on oxidative stress, proline metabolism and indole alkaloid accumulation. Colloids and Surfaces B Biointerfaces 60: 110–116. https://doi.org/10.1016/j.colsurfb.2007.06.006.

Jamali S (2015) First report of Septoria silybi associated with leaf blotch of Silybum marianum from Iran. Plant Science Today 2: 21–23. https://doi.org/10.14719/pst.2015.2.1.82.

Kashmir S, Khan MA, Shad AA, Marwat KB, Khan H (2016) Temperature and salinity affect the germination and growth of Silybum marianum Gaertn and Avena fatua L. Pakistan Journal of Botany 48: 469–476.

Kavallieratos NG, Tomanović Z, Sarlis GP, Vayias BJ, Zikić V, Emmanouel NE (2007) Aphids (Hemiptera: Aphidoidea) on cultivated and self-sown plants in Greece. Biology 62: 335–344. https://doi.org/10.2478/s11756-007-0056-x.

Khalil-ur-Rahman, Jahan N, Ahmed Z, Shakeel A, Aslam HKW, Sakandar HA, Raza MS (2016) Behavioral pattern of milk thistle (Silybum marianum L.) agronomic parameters to variable temperature and land zones of Punjab. Journal of Global Innovations in Agricultural and Social Sciences 4: 15–22.

Khan MA, Blackshaw RE, Marwat KB (2009) Biology of milk thistle (Silybum marianum) and the management options for growers in north-western Pakistan. Weed Biology and Management 9: 99–105. https://doi.org/10.1111/j.1445-6664.2009.00326.x.

Khatamipour M, Piri E, Esmaeilian Y, Tavassoli A (2011) Toxic effect of cadmium on germination, seedling growth and proline content of milk thistle (Silybum marianum). Annals of Biological Research 2: 527–532.

Kovαčikovα E, Kubνnek J (1986) Fading of the milk thistle (Silybum marianum L. Gaertn.). Proceedings of the Institute of Scientific and Technical Information for Agriculture, Plant Protection 22: 217–222.

Kurkin VA (2003) Saint-Mary thistle: A source of medicine. Pharmaceutical Chemistry Journal 37: 189–202. https://doi.org/10.1023/A:1024782728074.

Malekzade M, Mirmazloum SI, Angoran HR, Alirezalu A, Radacsi P (2011) Climate change: challenges and opportunities in agriculture. In: AGRISAFE Final Conference Proceedings, Budapest, Hungary, 21–23 March 2011. 187–190. Available at: https://www.cabidigitallibrary.org/doi/pdf/10.5555/20133092387.

Morazzoni P, Bombardelli E (1995) Silybum marianum (Carduus marianus). Fitoterapia 66: 3–42.

Moscow D, Lindow SE (1989) Infection of milk thistle (Silybum marianum) leaves by Septoria silybi. Phytopathology 79: 1085–1090. https://doi.org/10.1094/Phyto-79-1085.

Omer EA, Refaat AM, Ahmed SS, Kamel A, Hammouda FM (1993) Effect of spacing and fertilization on the yield and active constituents of milk thistle, Silybum marianum. Journal of Herbs Spices and Medicinal Plants 1: 17–23. https://doi.org/10.1300/J044v01n04_04.

Omidbaigi R, Karimzadeh G, Koshki MH (2003) A study on the influence of sowing date and plant density on the productivity of Silybum marianum (L.) Gaertn. and characteristics correlation. Iranian Journal of Science and Technology 1: 203–212.

Omidbaigi R, Nobakht A (2001) Nitrogen fertilizer affecting growth, seed yield and active substances of milk thistle (Silybum marianum). Pakistan Journal of Biological Sciences 11: 1345–1349. https://doi.org/10.3923/pjbs.2001.1345.1349.

Oudemans CAJA (1923) Enumeratio Systemica Fungorum, Vol. 4. The Hague: Apud Martinum Nijhoff.

Parsons WT (1973) Noxious Weeds of Victoria. Melbourne: Inkata Press.

Plant Protection and Quarantine (2002) Federal Noxious Weed List. USDA Animal and Plant Health Inspection Service. Available at: http://plants.usda.gov/cgi-bin/federal_noxious.cgi.

Pourreza J, Bahrani A (2012) Estimating cardinal temperatures of milk thistle (Silybum marianum) seed germination. American-Eurasian Journal of Agricultural and Environmental Sciences 12: 1030–1034.

Radjabian T, Rezazadeh SH, Fallah-Huseini H (2008) Analysis of silymarin components in the seed extracts of some milk thistle ecotypes from Iran by HPLC. Iranian Journal of Science and Technology 32: 141–146.

Ramasamy K, Agarwal R (2008) Multitargeted therapy of cancer by silymarin. Cancer Letters 269: 352–362. https://doi.org/10.1016/j.canlet.2008.03.053.

Rezwani A (2008) A new species, and a new record of a subspecies belonging to the genus Dysaphis (Homoptera: Aphididae) from Iran. Journal of the Entomological Society of Iran 27: 23–26.

Roche C (1991) Milk thistle (Silybum marianum (L.) Gaertn.). Pacific Northwest Extension Publication no. PNW382. Washington State University/Oregon State University.

Saccardo PA (1884) Sylloge Fungorum, Vol. 3. Berlin: Friedlander & Sohn.

Safikhan S, Khoshbakht K, Chaichi MR, Amini A, Motesharezadeh B (2018) Role of chitosan on the growth, physiological parameters and enzymatic activity of milk thistle (Silybum marianum L. Gaertn.) in a pot experiment. Journal of Applied Research on Medicinal and Aromatic Plants 10: 49–58. https://doi.org/10.1016/j.jarmap.2018.06.002.

Šafrαnkovα I, Kolαčkovα P, Růžičkovα G (2015) The pathogens of milk thistle (Silybum marianum [L.] Gaertn.) in the Czech Republic. Journal of Advances in Agriculture 5: 538–545. https://doi.org/10.24297/jaa.v5i1.4264.

Saller R, Meier R, Brignoli R (2001) The use of silymarin in the treatment of liver diseases. Drugs 61: 2035–2063. https://doi.org/10.2165/00003495-200161140-00003.

Sedghi, M, Gholipouri a, Sharifi RS (2008) γ-Tocopherolaccumulation and floraldifferentiation of medicinal pumpkin (Cucurbita pepo L.) in response to plant growth regulators. Notulae Botanicae Horti Agrobotanici Cluj-Napoca 36: 80–84. https://doi.org/10.15835/nbha361101.

Selmar D, Kleinwachter K (2013) Stress enhances the synthesis of secondary plant products: The impact of stress-related over-reduction on the accumulation of natural products. Plant and Cell Physiology 54: 817–826. https://doi.org/10.1093/pcp/pct054.

Shimi P, Poorazar R, Jamali M, Bagherani-Torshiz N (2006) Evaluating clopyralid as a broadleaf herbicide in canola fields of Iran. Pakistan Journal of Weed Science Research 12: 307–311.

Sindel BM (1991) A review of the ecology and control of thistles in Australia. Weed Research 31: 189–201. https://doi.org/10.1111/j.1365-3180.1991.tb01758.x.

Smith T, Gillespie M, Eckl V, Knepper J, Reynolds CM (2019) Herbal supplement sales in US increase by 9.4% in 2018. HerbalGram 123: 62–73.

Sofer-Arad S, Kigel H, Henkin Z and Kol M (2007) Reproduction and dispersion of annual thistles in natural pastures. In: Proceedings of the 15th Meeting Israeli Association for Range Science, March 14–15 2007, Ramat Hanadiv.

Solouki H, Rezvani Moghaddam P, Phiroozjaee GK (2015) Investigating the seed germination characteristics of milk thistle (Silybum marianum L.) affected by magnetic field, sodium chloride and hydro-priming. Advances in Plants and Agriculture Research 2: 1–6. https://doi.org/10.15406/apar.2015.02.00070.

Souissi T, Berner DK, Smallwood EL (2005) First report of smut caused by Microbotryum silybum on inovythistle. Plant Disease 89: 1242. https://doi.org/10.1094/PD-89-1242A.

Tamouridou AA, Pantazi XE, Alexandridis T, Lagopodi A, Kontouris G, Moshou D (2018) Spectral identification of disease in weeds using multilayer perceptron with automatic relevance determination. Sensors 18: 2770. https://doi.org/10.3390/s18092770.

Topalov V, Dechev I and Pehlivanov M (1983) Crop production of blessed thistle. In: Topalov VD, Dechev I and Pehlivanov M (eds) Crop Production. Sofia: Zemizdat Press, 421–424.

USDA-ARS Agricultural Research Service (2005) Noxious weeds in the US and Canada. Available at: http://invader.dbs.umt.edu/noxious_weeds/.

Vinograd A, Zaady E, Kigel H (2011) Dynamics of plant populations in abandoned sheep corrals: processes and management. In: Proceedings of the 19th Meeting Israeli Association for Range Science, Ramat Hanadiv, Israel, 14–15 March 2011, 48–49.

Xiao-Fang S, Qingsong Z, Youlinag L (2000) Regulation of salt tolerance of cotton plants at seedling emergence stage by soaking seeds in pix (DPC) and CaCl₂ solutions. Jiangsu Journal of Agricultural Sciences 16: 204–207.


Zahir A, Abbasi BH, Adil M, Anjum S, Zia M, Ihsan-Ul-Haq (2014) Synergistic effects of drought stress and photoperiods on phenology and secondary metabolism of Silybum marianum. Applied Biochemistry and Biotechnology 174: 693–707. https://doi.org/10.1007/s12010-014-1098-5.

Zahra N, Wahid A, Hafeez MB, et al. (2021a) Plant growth promoters mediated quality and yield attributes of milk thistle (Silybum marianum L.) ecotypes under salinity stress. Scientific Reports 11: 23200. https://doi.org/10.1038/s41598-021-02435-4.

Zahra N, Wahid A, Hafeez MB, Shaukat K, Shahzad S, Shah T, Alyemeni MN (2022) Plant growth promoters alleviate oxidative damages and improve the growth of milk thistle (Silybum marianum L.) under salinity stress. Journal of Plant Growth Regulation 41: 3091–3116. https://doi.org/10.1007/s00344-021-10498-w.

Zahra N, Wahid A, Shaukat K, Hafeez MB, Batool A, Hasanuzzaman M (2021b) Oxidative stress tolerance potential of milk thistle ecotypes after supplementation of different plant growth-promoting agents under salinity. Plant Physiology and Biochemistry 166: 53–65. https://doi.org/10.1016/j.plaphy.2021.05.042.

Zand E, Baghestani MA, Soufizadeh S, Eskandari A, PourAzar R, Veysi M et al. (2007) Evaluation of some newly registered herbicides for weed control in wheat (Triticum aestivum L.) in Iran. Crop Protection 26: 1349–1358. https://doi.org/10.1016/j.cropro.2006.10.011.

Zheljazkov V, Nikolov S (1996) Accumulation of Cd, Pb, Cu, Mn and Zn by Silybum marianum L. on polluted soils. Acta Horticulturae 426: 297. https://doi.org/10.17660/ActaHortic.1996.426.35.

Zheljazkov VD, Zhalnov I, Nedkov NK (2006) Herbicides for weed control in blessed thistle (Silybum marianum). Weed Technology 20: 1030–1034. https://doi.org/10.1614/WT-05-135.1.