Skip to main content
Canna~Fangled Abstracts

Origin, early expansion, domestication and anthropogenic diffusion of Cannabis, with emphasis on Europe and the Iberian Peninsula

By June 16, 2024January 5th, 2025No Comments


https://doi.org/10.1016/j.ppees.2022.125670Get rights and content
Under a Creative Commons license
open access

Highlights

  • Oligo-Miocene origin (NE Tibetan Plateau) and Mio-Pliocene preanthropic expansion.
  • Extant subspecies (C. sativa subsp. sativaC. sativa subsp. indica) of Pleistocene origin.
  • Early Neolithic (ca. 12 kyr BP) domestication in a single (E Asia) or more regions?
  • Wild forms present in Europe during the Pleistocene. European domestication?
  • Postglacial wild forms in the Iberian Peninsula (IP) during post-glacial times (18.5 kyr BP)?
  • IP widespread cultivation and retting only since the Middle Ages? Insufficient evidence.

Abstract

Cannabis is among the oldest human domesticates and has been subjected to intensive artificial (human-mediated) selection throughout history to create a wide array of varieties and biotypes for diverse uses, including fiber, food, biofuel, medicine and drugs. This paper briefly reviews the available literature on the taxonomy, evolutionary origin and domestication of this plant, as well as its worldwide dispersal, in both its wild and cultivated forms. Emphasis is placed on Europe and especially on the Iberian Peninsula. Today, it is accepted that Cannabis is a monospecific genus with two subspecies, C. sativa subsp. sativa and C. sativa subsp. indica, originating in Europe and Asia, respectively, by allopatric differentiation after geographic isolation fostered by Pleistocene glacial-interglacial cycles. Palynological and phylogeographic evidence situates the Cannabis ancestor on the NE Tibetan Plateau during the mid-Oligocene. The timing and place of domestication is still a matter of debate between contrasting views that defend single or multiple Neolithic domestication centers situated in different parts of the Eurasian supercontinent, notably central/southeastern China and the Caucasus region. Recent meta-analyses have suggested that wild Cannabis may have already been spread across Europe in the Pleistocene, and its domestication could have occurred during the European Copper/Bronze ages. According to the available reviews and meta-analyses, pre-anthropic dispersal of Cannabis into the Iberian Peninsula seems to have occurred only in postglacial times, and the earlier signs of cultivation date to the Early Medieval Ages. However, the palynological and archeological evidence used to date is insufficient for a sound assessment, and the development of thorough Iberian databases to address further meta-analysis is essential for more robust conclusions. Some clues are provided for these achievements to be fulfilled.

Keywords

Cannabis
Taxonomy
Biogeography
Evolution
Domestication
Dispersal
Pollen

1. Introduction

Since ancient times, Cannabis has been extensively utilized by humans for a variety of uses, such as textiles, paper, food, medicine, biofuel or recreational drugs. Cannabis has been an integral part of human life since its domestication and continues into the present (e.g., Clarke and Merlin, 2013Fike, 2016Gray et al., 2016aGrotenhermen and Müller-Vahl, 2016). However, research on the uses and applications of this plant significantly decreased and almost ceased in the mid-20th century when the species and its varieties were declared illegal in most Western countries, regardless of the concentrations of psychoactive compounds (Duvall, 2014Warf, 2014). Therefore, one of the most ancient crops, long-valued for its multiple uses, became largely ignored during the last 50 years. Thus, compared with other crop plants, Cannabis has not fully benefited from modern scientific technologies, which has created a major knowledge gap that remains to be filled (Gray et al., 2016b).
The last decade has witnessed a revival of interest in Cannabis research, especially in relation to the genetic and phytochemical features of the different varieties and to its evolutionary origin, domestication and further geographic diffusion (e.g., Clarke and Merlin, 2013Small, 2015Gray et al., 2016bVergara et al., 2016). For example, the use of modern methods of molecular analysis has helped clarify the taxonomy of the Cannabis complex and its genetic modifications over time as a result of artificial (human-mediated) selection during the domestication process (Clarke and Merlin, 2016). Palynological and archeological records have been used mostly to identify the geographical center of origin of Cannabis, its center of domestication and further human-driven diffusion patterns throughout the world. Modern global and regional databases and associated handling facilities, notably geographic information system (GIS) tools, have been instrumental in this type of research (McPartland, 2020).
This paper briefly reviews the newly available information on the center of origin and the center of domestication of Cannabis and its further worldwide dispersal, with a focus on Europe, where spatiotemporal diffusion patterns remain controversial, and the Iberian Peninsula, which represents a knowledge gap regarding the tempo and mode of Cannabis arrival and diffusion. Previous reviews on the subject date from roughly a decade ago (Clarke and Merlin, 2013Duvall, 2014Warf, 2014Small, 2015). A review including the research developed on Cannabis during the last decade, which has been fundamental for setting the present standards of knowledge, is lacking. This is especially true for the meta-analyses and molecular phylogenetic investigations developed during the last five years. The present review represents an updated outlook on the evolution, domestication and worldwide diffusion of Cannabis that includes the latest developments on this topic.
The review begins with a brief section regarding the current state of Cannabis taxonomy, ecology, phenology and human uses. The next section discusses the evolutionary center of origin and the time of appearance of the original wild Cannabis, followed by similar geographic and chronological considerations on Cannabis domestication and diffusion of its cultivated forms. The next section focuses on Europe, where domestication and diffusion patterns are still under discussion. The last section highlights the scarcity of information in regard to the Iberian Peninsula, in comparison with most European regions, evaluates the potential causes for this paucity and suggests how further research could contribute to bridging this knowledge gap. Time units used: Ma, million years before present; (k)yr BP, (kilo)years before present; CE/BCE, Common Era/Before Common Era.

2. Present-day status

2.1. Species

The genus Cannabis is in the family Cannabaceae, with 11 genera and approximately 170 species. Discussions on the existence of one (C. sativa) or two species (C. sativa and C. indica) started during the time of Linnaeus and Lamarck and have continued until very recently. Some morphological differences exist between the two taxa, as C. sativa is taller with a fibrous stalk, whereas C. indica is shorter with a woody stalk. There are also phytochemical divergences, as manifested in the tetrahydrocannabinol/cannabidiol (THC/CBD) ratio, which is higher in C. indica. The discontinuous geographic ranges – the sativa lineage in Europe and the indica lineage in Asia – would also support the occurrence of two different taxa within Cannabis. The question is whether these differences qualify for defining two separate species or two subspecific taxa. Using quantitative criteria based on key molecular DNA sequences (DNA barcode), McPartland (2018) concluded that the genetic differences are in the rank of subspecies and that the proper nomenclature is C. sativa subsp. sativa and C. sativa subsp. indica. Therefore, Cannabis is now considered a monospecific genus (Barcaccia et al., 2020Kovalchuck et al., 2020).

2.2. Biotypes

A number of new varieties have been obtained from these two original subspecies by artificial selection, which have different morphologies, phytochemical compositions, geographic distributions and uses (Clarke and Merlin, 2016Lynch et al., 2016Rahn et al., 2016). According to Small (2015), the consequences of artificial selection make it impossible to determine if unaltered primeval or ancestral populations still exist. This author recommends that C. sativa be recognized as a single species with a narcotic subspecies (drug) and a non-narcotic subspecies (hemp), each with domesticated and ruderal varieties. A similar approach is the definition of four major biotypes of cultivated Cannabis (NLH, BLH, NLD and BLD) with morphological and chemical differences submitted to different uses, either fiber/oil or drug production (see Table 1 for more information on these biotypes). There is also a hybrid between NLD and BLD, known as sinsemilla, which is highly psychoactive (high THC and low CBD) and, hence, is used mainly as a drug. The sinsemilla cultivars were developed in the New World and diffused worldwide. In contrast to the former biotypes, this hybrid is only cultivated, with no feral escapes or ruderal populations. The current distribution of these biotypes is shown in Fig. 1.

Table 1. The four biotypes of cultivated Cannabis as defined by Clarke and Merlin (2016), following the criteria of Small (2015). See Fig. 1 for the present geographical distribution of these biotypes.

Biotype Name Psychoactive THC CBD Use Origin Diffusion
NLH Narrow-leaf hemp Rarely Low High Fiber/oil Europe New World
BLH Broad-leaf hemp Mildly Low/Moderate High Fiber/oil East Asia Europe, New World
NLD Narrow-leaf drug Very High Low/absent Drug South Asia Africa, Europe, Middle East, New World
BLD Broad-leaf drug Moderately Moderate/high Moderate/high Drug Afghanistan Europe, New World
Fig. 1

  1. Download: Download high-res image (148KB)
  2. Download: Download full-size image

Fig. 1. Worldwide distribution of the different Cannabis biotypes developed by humans (Table 1). Hemp biotypes are in green and drug biotypes in red. NLH, narrow-leaf hemp; BLH, broad-leaf hemp; NLD, narrow-leaf drug; BLD, broad-leaf drug; NLD/BLD, hybrid between NLD and BLD; PA?, putative cannabis ancestor. Redrawn from Fig. 2 of Clarke and Merlin (2016).

2.3. Ecology and phenology

Cannabis sativa is a sun-loving (heliotropic) species that requires well-drained and nitrogen-rich soils, warmth and moisture. Therefore, most natural populations are found seasonally across accommodating northern temperate latitudes. This plant grows well along exposed riverbanks, lakesides, margins of agricultural lands and other areas disturbed by humans. Cannabis plants are annual and usually diecious, as determined by X and Y chromosomes, and anemophylous (wind-pollinated). The annual cycle extends from spring (germination) to summer (fast juvenile growing) and autumn (flowering). Male plants, which are slightly taller than female plants, die shortly before pollination. Female plants ripen viable seeds just before the arrival of winter killing frosts. Seed dissemination is carried out mostly by wind or feeding birds. During germination, seeds are surrounded by bracts with hairs that produce a resinous blend of cannabinoids and aromatic compounds as secondary metabolites, which are believed to protect seeds against pests and pathogens. Cannabinols may be psychoactive (THC) or not psychoactive (CBD) for humans (Small, 2015Clarke and Merlin, 2016).

2.4. Uses

Almost all parts of the Cannabis plants are utilized for a variety of uses (Table 2). For example, stem bark and fiber are used for cordage, woven textiles, building materials, paper, animal bedding and fuel. Seeds and seed oils are used for human food, animal feed, industrial feedstock and fuel. Female flowers and seeds are used for medicine or recreational drugs. All parts of the plant, primarily bark, seeds and female flowers, are used in ritual and social activities (healing and life cycle rituals, inebriation). Cannabis populations may also be used for environmentally related activities, such as the control of soil erosion and to increase CO2 sequestration. Esthetic and educational applications include botanical gardens and the iconic character of the plant as a symbol of a very ancient crop deeply rooted in human culture (Gray et al., 2016b). An extensive database on the traditional uses of Cannabis (CANNUSE) is available at http://cannusedb.csic.es/. See Balant et al. (2021) for detailed explanations on this database.

Table 2. Examples of the variety of uses of the different parts of the C. sativa plants. Modified from Clarke and Merlin (2016).

Plant parts Use category Material type or benefits
Stem bark Cordage Long cellulose fibers
Stem fiber Cordage, woven textiles, building materials Long cellulose fibers, concrete reinforcement
Wood/bark Paper, building materials, animal bedding, fuel Long/short cellulose fibers, chip board, concrete matrix, heat, light
Female flowers/seeds Medicinal Herbal remedies, pharmaceuticals, nutraceuticals
Female flowers and associated resin glands Recreational drugs Marijuana, hashish
Seeds (oil) Human food, industrial feedstock, fuel Proteins/essential fatty acids, paint/plastic manufacture
Seeds (cake), foliage Animal feed Proteins and essential fatty acids
All parts Ritual and social Healing, life cycle rituals, inebriation
Plant-people interplay Esthetic Intrinsic beauty of the plant
Genus Educational Iconic example of an economic plant and its ancient human relationships

3. Time and place of origin

Early attempts to identify the place of origin of wild Cannabis prior to human contact were based on the geographic distribution of its wild, cultivated and ruderal populations, combined with the known ecological requirements and reproductive strategies. However, the distribution of this plant and its biotypes/varieties is closely associated with human settlements and trade routes, and therefore, the original native range is obscured (Clarke and Merlin, 2013). In spite of this, a broad area referred to as central Asia (presently China) was proposed as the center of origin of Cannabis (Schultes, 1969Merlin, 1972). Another, less generally accepted, possibility mentioned was south Asia (presently India). Regarding timing, accepting that the central Asian steppes were colonized by humans by 35,000 years ago (Wells, 2002Finlayson, 2005), it has been assumed that wild Cannabis could have originated earlier. These hypotheses, however, were based on circumstantial evidence, and robust empirical evidence was lacking. This empirical evidence was provided by the fossil record and the use of time-calibrated molecular DNA phylogenies.

3.1. Fossil record

The macrofossil record of Cannabis is relatively scarce and consists of only a few leaf and fruit/seed impressions with ages ranging between the Oligocene and the late Miocene (McPartland et al., 2019). However, microfossils, specifically pollen, are abundant and widespread and have commonly been utilized to reliably reconstruct the history of Cannabis. Nevertheless, the identification of Cannabis pollen deserves special attention because of its similarity with other members of the family Cannabaceae, especially Humulus (hop), a sister genus that bears contrasting ecological requirements and cultural connotations. Therefore, inaccurate identification might lead to erroneous conclusions (Rull and Vegas-Vilarrúbia, 2014). This is why different authors have used broader taxonomic categories for this pollen type, such as Cannabis-type, Cannabis/Humulus or Cannabaceae. In addition, it is unclear whether pollen from wild and cultivated Cannabis may be distinguished morphologically.
Several morphological details have been suggested to be useful to differentiate Cannabis and Humulus pollen. Godwin (1967) emphasized several differential characteristics of the pore complex. Further statistical studies revealed that Cannabis pollen is generally larger than Humulus, but this character alone was not sufficient to allow reliable separation (Whittington and Edwards, 1989Whittington and Gordon, 1987). Pollen size also seems to be an unreliable parameter to separate wild from cultivated Cannabis (review in McPartland et al., 2018). The combination of pore complexes and size seems to provide a more reliable, yet not universally accepted, identification criterion to differentiate between Cannabis and Humulus pollen (Fleming and Clarke, 1998Mercuri et al., 2002).
However, even in the case of conclusive Cannabis pollen identifications, comparisons with studies referring to this pollen type as Cannabis-type, Cannabis/Humulus or Cannabaceae remain problematic, which may be a handicap for the development of meta-analyses aimed at reconstructing past biogeographic and cultural patterns. Recently, some meta-analyses have been conducted using different criteria, such as considering the entire Cannabis/Humulus complex or taking into account only those studies that explicitly identified Cannabis pollen (Clarke and Merlin, 2013Long et al., 2017). These approaches tend to overestimate or underestimate the actual pollen record of Cannabis.
A different approach, called here the assemblage approach, has recently been proposed by McPartland et al. (2018). These authors noted that wild C. sativa is typical of open temperate steppe habitats dominated by grasses, chenopods and Artemisia, whereas Humulus is a vine plant that requires trees to climb and is common in temperate deciduous forests dominated by alder (Alnus), willow (Salix) and poplar (Populus). A third assemblage corresponds to cultivated Cannabis, which is usually found together with cultivated cereals such as Avena (oats), Hordeum (barley), Secale (rye) and Triticum (wheat), as well weeds such as Centaurea (cornflowers) or Scleranthus (knawels) species. Using these phytosociological affinities, McPartland et al. (2018) attributed the Cannabis-like pollen types recorded in the literature to wild Cannabis if this pollen occurred together with steppe assemblages; to cultivated Cannabis when it was part of crop assemblages; or to Humulus if the dominant pollen assemblage corresponded to temperate deciduous forests. Based on these premises, these authors developed a more complicated identification algorithm that also considers the relationship between arboreal (AP) and non-arboreal (NAP) pollen (Fig. 2).
Fig. 2

  1. Download: Download high-res image (104KB)
  2. Download: Download full-size image

Fig. 2. Algorithm developed to differentiate between Cannabis and Humulus pollen in sedimentary records, using the assemblage approach. As explained in the text, crop pollen usually includes cereal and weed pollen. AP, arboreal pollen; NAP, non-arboreal pollen; PAC, steppe assemblage (Poaceae-Artemisia-Chenopodiaceae); ASP, forest assemblage (Alnus-Salix-Populus). Original figure based on Fig. 1 of McPartland et al. (2018).

Using these criteria, the oldest known pollen compatible with Cannabis was found in 19.6 Ma-old (early Miocene) rocks from the NE Tibetan Plateau (presently Ningxia, China), which was proposed as the center of origin of Cannabis (McPartland et al., 2019). Interestingly, this proposal roughly coincides with the former hypothesis based on indirect biogeographic evidence.

3.2. Calibrated DNA phylogenies

The use of DNA molecular phylogenies calibrated with fossils of related genera such as HumulusCeltisMorus and Ficus (Fig. 3) allowed to estimate the age of divergence of Cannabis and Humulus to 27.8 Ma (mid-Oligocene). Using the same DNA phylogeny and the associated molecular clock, the divergence between C. indica (or C. sativa subsp. indica) and C. sativa (C. sativa subsp. sativa) would have occurred in the Middle Pleistocene approximately 1 Ma (McPartland, 2018). These authors noted that there is a gap of ca. 8 million years between the age of origin estimated by the molecular clock and the first fossil pollen encountered. In spite of this, they favored the mid-Oligocene age for the origin of Cannabis on the NE Tibetan Plateau (Fig. 4), assuming that, as demonstrated by the presence of Artemisia and other steppe elements, the region was covered by this type of vegetation, which would have been particularly well suited for the development of Cannabis. A previous study using Bayesian calibration estimated the divergence between Cannabis and Humulus to have occurred 21 Ma (Zerega et al., 2005), which is closer to the age (19.6 Ma) of the first fossil pollen evidence mentioned above (McPartland et al., 2019).
Fig. 3

  1. Download: Download high-res image (121KB)
  2. Download: Download full-size image

Fig. 3. Time-calibrated molecular DNA phylogenetic tree used for the estimation of the divergence rates within the Cannabaceae and some related families (molecular clock). Red dots are the nodes used for calibrating the phylogenetic tree using HumulusCeltisMorus and Ficus fossils of known age. Numbers are ages in million years before present (Ma). Redrawn from Fig. 2 of McPartland (2018).

Fig. 4

  1. Download: Download high-res image (168KB)
  2. Download: Download full-size image

Fig. 4. Center of origin (red dot) and pre-anthropic dispersal of Cannabis, based on fossil pollen data and the use of the assemblage algorithm explained above (Fig. 2). Dot colors indicate the age of first occurrences of Cannabis pollen (see legend), in million years before present (bold numbers) and in thousand years before present (normal numbers). Original figure based on raw data from McPartland et al. (2019).

3.3. Pre-anthropic expansion

According to pollen and seed fossil records, Cannabis would have experienced some expansion from its center of origin to Europe and East Asia well before the evolutionary appearance of the genus Homo (Fig. 4). The suggested dispersal agents are water (hydrochory) and animals (zoochory) (McPartland, 2018McPartland et al., 2019). The first expansions occurred to the west (eastern Europe) and the east (north-eastern China) during the Miocene-Pliocene. Most parts of the Asian continent were colonized by Cannabis during the Pleistocene (the last 2.6 Ma), before the onset of the Neolithic, when humans domesticated the first plants. During the Pleistocene, glacial-interglacial recurrence could have contributed to Cannabis diversification without human intervention. It has been suggested that Cannabis underwent recurrent range contractions (glacials) and expansions (interglacials) that facilitated allopatric processes, possibly leading to the differentiation between the European (C. sativa subsp. sativa) and Asian (C. sativa subsp. indica) subspecies, which would have diverged nearly 1 Ma (Clarke and Merlin, 2013McPartland, 2020). The first is considered the putative hemp ancestor (PHA), and the second is the putative drug ancestor (PDA) (Clarke and Merlin, 2013).

4. Domestication and diffusion

In addition to fossils and DNA phylogenies, archeological evidence is of paramount importance to reconstruct domestication and anthropogenic diffusion trends within Cannabis. The main types of evidence of Cannabis from archeological sites are pollen, seeds, fibers, fiber/seed impressions, carbonized remains, phytoliths and chemical remains. During the historical period, written and graphic documents are also of fundamental help (Clarke and Merlin, 2013). Pollen identification has some additional clues related to its abundance in sediments, as Cannabis produces much more pollen than Humulus, which is usually underrepresented (Lewis et al., 1983). This is especially useful in cases of very high percentages of this pollen type, which are difficult to explain unless the sediments come from a former hemp-retting site. Indeed, when flowering male hemp plants are soaked in a retting pond to separate the fibers from the stalk, large quantities of pollen settle into pond sediments. In these situations, the percentages of Cannabis pollen in sediments may reach 80–90% of the total, but percentages over 15% or 25% have been considered sufficient to infer hemp retting (Peglar, 1993Mercuri et al., 2002Lavrieux et al., 2013Demske et al., 2016). The anemophylous pollen of Cannabis may be transported long distances; therefore, the finding of a few grains or their scattered occurrence throughout a stratigraphic section is not necessarily evidence for the local presence of the parent plant. Recent developments in molecular DNA analytical methods have increased the probability of identifying Cannabis sativa in lake sediments and comparing these results with the abundance of Cannabis/Humulus pollen. For example, in lake sediments from the French Alps, Giguet-Covex et al. (2019) found a good agreement between DNA and pollen records for the period 1500–1000 yr BP but significant inconsistencies for the periods 2000–1500 yr BP and 500 yr BP-present. During these periods, the Cannabis/Humulus pollen attained values of 10–15% and DNA was undetectable, which could be attributed to the dominance of Humulus pollen in the pollen record.

4.1. Domestication center(s)

The debate regarding the Cannabis center of domestication has paralleled the taxonomic controversy on Cannabis species that, as seen before, was already active by the time of Linnaeus and Lamarck. See, for example, Clarke and Merlin, 2013Clarke and Merlin, 2016 and Small (2015) for more details on this long-standing discussion. Whether single or multiple centers of Cannabis domestication occurred has profound implications for the tempo and mode of Cannabis diffusion in Eurasia, where the plant evolved and was domesticated.
Some authors have proposed that Cannabis was domesticated in Central Asia (presently Mongolia and northern China) by 12,000 BCE (ca. 10,000 yr BP), which would place Cannabis among the oldest human domesticates (Warf, 2014). A recent genome-wide phylogeographic study supports a unique center of domestication in east Asia (China), from which all biotypes would have emerged and dispersed throughout the world (Ren et al., 2021). According to this analysis, early domesticated ancestors of hemp and drug types diverged from wild Cannabis ca. 12,000 yr BP, which indicates that the species had already been domesticated by early Neolithic times. Among the defenders of the single center hypothesis, Small (2015) suggested that Cannabis was domesticated in the northern Caspian Sea and that the four different domesticated groups were transported to other parts of the world during the last millennium. These four groups coincide with the abovementioned biotypes NLH (Europe), BLH (E Asia), NLD (South-Central Asia) and BLD (Afghanistan). Other authors propose two centers of domestication, one in the Caucasus (between the Caspian and the Black Sea), where the ancestor of hemp fiber domesticated biotypes (PHA) evolved, and another in south-eastern China, the place where the precursor of narcotic domesticates (PDA) originated (Fig. 5). The idea of at least two domestication centers situated in Europe and Asia is supported by archaeobotanical studies using pollen, seeds and fibers (Long et al., 2017McPartland et al., 2018).
Fig. 5

  1. Download: Download high-res image (318KB)
  2. Download: Download full-size image

Fig. 5. The first two phases of human-mediated Cannabis diffusion. As in Fig. 1, hemp biotypes are indicated by green letters and drug biotypes are in red letters. Original figure based on Maps 10 and 11 from Clarke and Merlin (2013).

4.2. Anthropogenic diffusion

The human-mediated diffusion of Cannabis outside Eurasia occurred relatively recently and is well documented historically. After a thorough review of the available archeological and historical evidence, Clarke and Merlin (2013) subdivided anthropogenic diffusion of Cannabis dispersal beyond Eurasia into six phases: 1) primary dispersal across Eurasia (10,000–2000 yr BP), 2) spread into Africa and SE Asia (2000–500 yr BP), 3) diffusion into the Americas from Europe (1545–1800 CE), 4) diffusion to the Americas from Europe and Asia (1800–1945), 5) expansion after World War II (1945–1990), and 6) proliferation of industrial hemp (1990 CE to present).
In the first phase, Paleolithic nomadic peoples – who were expanding their range during the postglacial ice retreat – could have contributed, either consciously or accidentally, to the dispersal of the wild hemp (PHA) and drug (PDA) precursors from their corresponding distribution areas located in the Caucasus region and southeastern China, respectively (Fig. 5). This was the onset of differentiation of the four major cannabis biotypes. NLH evolved from the PHA via an intermediate form known as the NLH ancestor (NLHA), whereas NLD and BLD originated from the PDA through the intermediate ancestral form called the NLDA. BLH could also have evolved from isolated PDA populations in north-east China. The second phase coincided with the expansion of the Arab Empire into Africa and the Indian Empire into SE Asia. These cultures used Cannabis mainly as a drug and introduced the NLD biotype in their dominions, primarily through adventurers and traders. After the initial introductions, migrants and traders spread the NLD biotype across much of Africa and Asia, whereas the hemp biotype (NLH) remained restricted to Europe until the colonization of the Americas.
Diffusion into the New World (phase 3) occurred during European colonization and consisted mainly of the production of hemp (NLH) for cordage, cloth and seed (Fig. 6). It is speculated whether Vikings could have carried hemp to North America in pre-Columbian times, but there is no evidence for a sound assessment. In addition to domestic uses, hemp cultivation was stimulated by European governments of the colonizing countries to provide their sailing ships with ropes and sails. During this phase (1545–1800 CE), drug biotypes remained in the Old World, where the African expansion of NLD continued. It was not until the middle 1800s CE (phase 4) that Asian hemp (BLH) was introduced in North America through the Pacific Ocean (Fig. 6) to replace the inferior European hemp (NLH) in terms of fiber quality. The first psychoactive biotypes (NLD) were introduced from India into the Americas (mainly Latin America) by indentured laborers following the abolition of slavery in 1834.
Fig. 6

  1. Download: Download high-res image (279KB)
  2. Download: Download full-size image

Fig. 6. Phases three to six of anthropogenic cannabis expansion. As in Fig. 1, hemp biotypes are indicated by green letters and drug biotypes are in red letters. Original figure based on Maps 12, 13 and 14 from Clarke and Merlin (2013).

After World War II (1945–1949), Cannabis cultivation and use were prohibited at the international level, which signified the end of legal diffusion of this plant and initiated the fifth spreading phase (Fig. 6). The illegal traffic of marijuana (NLD) for recreational use did not stop but has experienced a significant increase since the 1960s. Most producers of clandestine Cannabis were from Africa, Asia and Latin America, which inadvertently exported seeds to the main consuming centers from Europe and North America, where local cultivation attempts began. In the 1970 s, Asian varieties of BLD were introduced into North America and western Europe, where the cultivation of NLD × BLD hybrids flourished.
During the last three decades (phase 6), Cannabis prosecution has generated new cultivation and dispersal modes to avoid legal detection. At the same time, renewed economic interest in hemp fiber has been sparked, which has promoted the legal industrial cultivation of hemp. Medical applications have also increased by artificial selection toward varieties with increased THC and other cannabinoids. New varieties have been created by hybridization of NLH and BLH biotypes aimed at expanding hemp cultivation to equatorial areas where European NLH does not grow well. Regarding drugs, the production of NLD/BLD hybrids has also increased and disseminated worldwide, and vegetative indoor cultivation of these hybrids has proliferated. As a result, the global genetic diversity of both hemp and drug forms experienced a significant reduction in comparison with the 1970s and 1980s, when the different biotypes were cultivated by traditional farmers in isolated geographic habitats.

5. Cannabis in Europe

According to the single center of domestication hypothesis, European Cannabis would have been carried by humans from Asia already in a cultivated form. However, the finding of Pleistocene Cannabis pollen across Europe suggests that this plant would had arrived to this continent before its domestication, and therefore, the plant could have also been domesticated in Europe. Using the assemblage approach for the identification of C. sativa pollen explained above (Fig. 2), McPartland et al. (2018) performed a meta-analysis on nearly 480 sites and demonstrated that this pollen was widespread across Europe during the Pleistocene prior to human agency (Fig. 7). According to the same authors, no evident signs of Cannabis cultivation have been found either in the pollen or in the archeological record for most of the European Neolithic. The first palynological evidence for potential Cannabis cultivation appeared in present-day Bulgaria during the late Neolithic-Copper age and expanded during the Bronze and Iron ages (Fig. 7). Cannabis cultivation expanded across the entire continent between 2000 and 800 yr BP during the Roman Empire and the Early Middle Ages. These results were supported by a further meta-analysis of archeological evidence (textiles, cordage, fiber, seeds, pottery impressions and phytoliths) from almost 140 sites (McPartland and Hegman, 2018). These authors concluded that Cannabis dispersed from Asia to Europe in its wild form during the Pleistocene and was domesticated in situ during the Copper and Bronze ages. This challenges the single-site domestication hypothesis, according to which Cannabis would have been domesticated in Asia and expanded worldwide as a cultivated plant.
Fig. 7

  1. Download: Download high-res image (304KB)
  2. Download: Download full-size image

Fig. 7. Pleistocene (upper panel) and Holocene (lower panel) European pollen records consistent with Cannabis, according to the assemblage identification approach (Fig. 2). Note that, in the Pleistocene panel, post-glacial ages have been differentiated from the rest of Late Pleistocene ages. PHA is the precursor of the European Cannabis according to Clarke and Merlin (2013); compare with Fig. 5, phase 1. Ages in million years before present (bold numbers) and in thousand years before present (normal numbers). Original figure based on raw data from McPartland et al. (2018).

The possibility of an European domestication of Cannabis is consistent with the former view of Clarke and Merlin (2013), who proposed that the precursor of the European NLH biotype – the wild PHA, related to C. sativa subsp. sativa – would have been centered in the Caucasus region and expanded to western Eurasia and eastern Europe during the Pleistocene (compare Fig. 7 with Fig. 5, phase 1). Chronologically, this possibility also seems plausible, as C. sativa subsp. sativa would have diverged from its Asian sister C. sativa subsp. indica in the Middle Pleistocene, approximately 1 Ma (McPartland, 2018), and the European records of wild Cannabis are younger than this divergence. In addition, the ages of these wild European records show consistent westward and southward dispersal patterns during the Middle and Late Pleistocene, with younger postglacial records situated in the west and south (Fig. 7). Notably, the European domestication of Cannabis (7–5 kyr BP) estimated from pollen and archeological records (McPartland et al., 2018McPartland and Hegman, 2018) occurred after the European Neolithic and much later than Asian domestication (12 kyr BP), as estimated by calibrated phylogenetic trees (Ren et al., 2021).

6. The Iberian Peninsula

The Iberian Peninsula (hereafter IP) is a key biogeographic region due to its high biodiversity and endemism levels, along with its transitional character between contrasting biogeographical regions, which results in a very peculiar biota in the European context. The IP is one of the main centers of Mediterranean plant diversity, together with the Anatolian, the Balkan and the Italian peninsulas (Médail and Quetzel, 1997). These peninsulas have had a fundamental role not only as biodiversity cradles but also as glacial refugia, especially during the Last Glacial Maximum (LGM) (Hewitt, 1999). It has been estimated that roughly a quarter of the LGM refugia have been located in the IP and the Balearic Islands (Médail and Diadema, 2009).
The present Iberian vascular flora is especially rich – 6276 species (739 non-native), distributed in 1278 genera and 189 families – and shows a high degree of specific endemism (~ 23% of the native flora). This represents more than 50% of the European flora. The main Iberian centers of diversity and endemism are situated in the mountain ranges and the richest groups are the Compositae, the Leguminosae and the Gramineae, whereas the Gymnosperms and the Pteridophytes are the less represented (Aedo et al., 2017). Traditionally, the IP has been subdivided into two major biogeographical domains characterized by contrasting bioclimatic features: the Mediterranean and the Eurosiberian regions. The Mediterranean region is characterized by the occurrence of Mediterranean macrobioclimates, with at least two consecutive arid months during the warmest season (summer), and dominates most of the IP (Fig. 8). The Eurosiberian region is restricted to the N and NW sectors of the peninsula (~ 20% of the total surface) and is characterized by wet Temperate macrobioclimates, without a period of two or more consecutive with summer aridity (Rivas-Martínez et al., 2017). The Mediterranean region is dominated by species distributed across southern Europe and northern Africa, whereas the Eurosiberian region is characterized by species distributed across central and northern Europe (Aedo et al., 2017).
Fig. 8

  1. Download: Download high-res image (160KB)
  2. Download: Download full-size image

Fig. 8. Map of the Iberian Peninsula indicating the location of the sites included in the available reviews and meta-analyses (Clarke and Merlin, 2013McPartland and Hegman, 2018McPartland et al., 2018). The Eurosiberian bioclimatic region is in green and the Mediterranean region in yellow. The numbers in brackets above the site names are the age of the first appearance of the Cannabis/Humulus pollen, in kyr BP. Original figure based on raw data from the following references: Alcúdia (Burjachs et al., 1994); Andorra (Ejarque et al., 2010); Antas (Pantaleón-Cano et al., 2003); Algendar (Yll et al., 1997); Banyoles (Pérez-Obiol and Julià, 1994); Coll del Moro (Alonso and Juan, 1994); Estanya (Riera et al., 2004); La Roya (Allen et al., 1996); Montcortès (Rull et al., 2011); Rascafría (Franco-Mújica et al., 1998); Somolinos (Currás, 2012); Totana (López, 1988).

Land use is also strongly influenced by bioclimatic and biogeographic patterns. In the Eurosiberian region, land use is comparable to the rest of Atlantic Europe, where the traditional agrarian economy is based on cattle raising for milk and meat, combined with small-scale crops such as maize, potatoes and other vegetables, as well as fruits such as apples and chestnuts. The resulting landscape is a mosaic of small fields, woodlands and pastures, combined with larger areas of heathlands, forests and grasslands, especially in higher elevations. The Mediterranean Iberia is dominated by a totally different landscape resulting from land-use practices adapted to the occurrence of the climatic summer drought. In these conditions, crops not needing irrigation are favored. Annual crops, especially cereals, largely dominate the central IP, which is characterized by huge medium-elevation flat terrains. In the southern half, where climates are not so cold, typical Mediterranean crops such as vineyards and olive grows dominate the landscape. In the Mediterranean region, extensive crops are the norm and animal husbandry occupies a subordinate position, except in some low-mountain areas, where extensive parkland landscapes dominate (Loidi, 2017). The Mediterranean biome and its climatic features are better suited than its Eurosuberian counterpart for the growth of Cannabis (Section 2.3).
The oldest human remains found in the Iberian Peninsula date from approximately 1.2 Ma (Early Pleistocene) and correspond to Homo antecessor, an autochthonous species related to H. erectus (Bermúdez de Castro et al., 1997). The Neanderthals entered the IP by 200 kyr BP (Early Paleolithic) and the modern humans did the same by 40 kyr BP. The Neolithization of the IP began at approximately 5500 BCE (~ 7500 yr BP), when the former hunter-gatherer nomadic societies were replaced by more stable agricultural-based cultures (García-Martínez de Lagrán, 2015). During the Bronze Age (2700–800 BCE) and the Iron Age (800–70 BCE), a diversity of cultures developed on the IP, which were influenced by Indo-European migrations, from the north (first millennium BCE) and the Phoenicians, arriving by the Mediterranean coasts during the 11th century BCE. The autochthonous Iberian culture fully developed during the 7th-6th centuries BCE but colonizations and invasions from the Mediterranean area continued with the Greeks (8th century BCE), the Carthaginians (5th century BCE) and the Romans, who arrived by 220 BCE and occupied the IP for almost 600 years. After a transitional phase in which the IP was occupied by the Germanic Visigoths, coming from the north, the next long-lasting invasion was that of the Muslims, who entered the IP from Africa and dominated Iberia for almost eight centuries (711–1492 CE). The Middle Ages (5th–15th centuries CE) signified the progressive expansion of Christian cultures and the unification of monarchies, followed by the consolidation of the modern geographical and political configuration. After the Middle Ages, Spain and Portugal, the major countries of the IP, shifted from colonized to colonizer countries – favored by the expansion of their naval industry, which largely benefited from the development of local the hemp industry – and expanded their dominions all over the world, with emphasis on the American continent, a situation that remained until the 19th–20th centuries (Marín, 1995).
This brief summary shows that the IP is a strategic spot also in cultural and historical terms, mainly due to its geographic location, which has promoted the interaction of a diversity of cultures over prehistoric and historic times, favored by numerous and varied land and sea connections. Therefore, it might be expected that the IP would also be a relevant place for the study of a plant such as Cannabis, which has been subjected to intensive and extensive natural and human-mediated dissemination. It should be noted, however, that the IP is far from any proposed center of origin for Cannabis, and the wild form of this plant seems to have reached the peninsula only in postglacial times (after 18.5 kyr BP), after a long journey from the Caucasus region (Fig. 4Fig. 5Fig. 6). Similarly, the former continental-scale meta-analysis of fossil pollen and archeological evidence (McPartland et al., 2018McPartland and Hegman, 2018) indicates that the IP was among the last European areas – along with the British Isles and the Scandinavian region – where Cannabis was cultivated, and this occurred during the Roman Empire and the Early Middle Ages (2000–800 yr BP).
However, these continental-wide analyses include only a few IP sites in their compilations. For example, Clarke and Merlin (2013) mentioned a single, yet representative, site (Lake Estanya; Fig. 8) where Cannabis was submitted to intensive cultivation and retting and used it to propose the arrival of Cannabis to the IP by 600 CE. McPartland et al. (2018) retrieved their Iberian sites from the European Pollen Database (http://www.europeanpollendatabase.net/index.php) and other internet facilities and considered 10 IP sites, including two from the Balearic Islands (Fig. 8). In this way, although the incoming of Cannabis to the IP may be placed within a rather general European context, the precise timing and pathways of arrival, as well as the internal and external dispersal trends, remain largely unknown. Due to the peculiar climatic, biogeographic and cultural features of the IP mentioned above, a thorough survey focused on this specific region to reconstruct Cannabis history is worth attempting.

6.1. Relevant case studies

The potential of the IP for providing useful information in this sense is high, as demonstrated by some studies that have recorded relevant historical developments in the hemp industry based on pollen analyses of lake sediments. Among them, the most significant examples were found in lakes Estanya and Montcortès (Fig. 8), where hemp was exploited since the Early Middle Ages. The advantage of these sites is that the amount of Cannabis pollen deposited in the sediments is consistent with local cultivation and/or retting rather than with dispersal from regional and long-distance sources. Indeed, some aerobiological studies have demonstrated that relevant amounts of north-African Cannabis pollen can be transported to the southern IP in a few days (Cabezudo et al., 1997). Similarly, a recent study of modern sedimentation in Lake Montcortès (Fig. 8) recorded significant percentages of Cannabis pollen, but the parent plant was absent from the region (Rull et al., 2017). This long-distance dispersal ability implies that only the presence of this pollen is not sufficient to infer the local occurrence of Cannabis populations. However, this is not the case for Estanya and Montocortès, as we discuss at following.
In the sediments of Lake Estanya, the first appearance of Cannabaceae pollen was recorded ca. 600 CE (Riera et al., 2004) during the Early Middle Ages (Fig. 9). From then, the record was continuous until the 1990s. Hemp pollen was accompanied by other cultivated plants, such as Olea (olive tree), Secale (rye) and other undifferentiated cereal pollen. Hemp percentages remained relatively stable (< 10%) until the 14th century, when they underwent a significant increase of up to 25% around the middle of the 18th century (Modern Age). This increase was interpreted in terms of hemp retting in the lake, coinciding with a local increase in the cultivation of this plant (Riera et al., 2004). This was supported by the study of proxies for water quality and non-pollen palynomorphs (NPP) (Riera et al., 2006). This phase coincided with the maximum hemp production in Spain due to the high demand from the Spanish navy (Riera et al., 2006). After these dates, hemp pollen decreased abruptly to values below 10% during the 20th century. The authors attributed this hemp crisis to a general decrease in cultivation due to the depopulation of the area during the first half of the 20th century.
Fig. 9

  1. Download: Download high-res image (175KB)
  2. Download: Download full-size image

Fig. 9. Pollen diagram of selected cultivated plants from Lake Estanya record, which was dated using radiocarbon and known historical events. Dashed lines are the percentage boundaries (15% and 25%) proposed by different authors for inferring retting practices (see text). The red arrow indicates the first appearance of Cannabaceae pollen. The boundaries of the cultural phases are approximate. Original figure based on raw data from Riera et al. (2004).

The first pollen records of Cannabis pollen in Lake Montcortès sediments also occurred at approximately 600 CE and coincided with the disappearance of cereals (Fig. 10), which indicates a shift in local cultivation practices from cereals to hemp (Rull et al., 2021). Further increases in wild grasses and weeds such as Artemisia (mugworts) and Plantago (plantains) are consistent with the expansion of pastures. Hemp pollen experienced three main phases of abundance, separated by two phases of scarcity. The first two phases (ca. 600–750 CE and ca. 850–250 CE) were interpreted in terms of low-intensity cultivation/retting to cover local needs for fiber. The third phase, however, was characterized by significantly high abundances (50–60%) and was difficult to explain in terms of only local consumption (Rull et al., 2021). As in Lake Estanya, this phase (ca 1550–1850 CE) was coeval with the maximum development of the Spanish navy and, as a consequence, of hemp cultivation across the entire country. The same abrupt decrease in hemp pollen was recorded at the end of the 19th century, which has been related to the dismantling of the royal navy, the onset of hemp importation from other countries, the substitution of hemp fiber by other materials such as cotton and synthetic fibers, and the decrease in human pressure (Rull and Vegas-Vilarrúbia, 2014Trapote et al., 2018). The further increase in hemp pollen in the late 20th century may have been due to the renewed interest in hemp, likely favored by EU subsides (Trapote et al., 2018).
Fig. 10

  1. Download: Download high-res image (219KB)
  2. Download: Download full-size image

Fig. 10. Pollen diagram of selected cultivated plants and weeds from Lake Montcortès sediments, which were dated up to ca. 3000 yr BP using varves. Dashed lines are the percentage boundaries (15% and 25%) proposed by different authors for inferring retting practices (see text). The red arrow indicates the first appearance of Cannabis pollen. Original figure based on raw data from Rull et al. (2021).

The combination of the meta-analysis discussed above (Fig. 7) and the case studies of Estanya and Montcortès may suggest that wild Cannabis reached the IP during the postglacial period (18.5–15 kyr BP) and that cultivated Cannabis entered much later, by 600 CE. Notably, both entries would have proceeded from the north-eastern sector. However, in the present state of knowledge, it is still premature to confirm these assessments. The lag of sufficient localities also hinders knowing what happened on the IP with Cannabis during the large gap between post-glacial times and the Middle Ages. The development of a thorough database for the IP, as a basis for further meta-analyses, is essential to understand when and how wild and cultivated Cannabis reached the IP, as well as what happened since those times.

6.2. Insights for future research

In addition to the information available from the above reviews and meta-analyses, other sources of information should be accessed. For example, many other sites are available in the compilation by Carrión (2012) that have not been included in the former studies. This compilation gathered almost all pollen records available for the IP by the time of publication and is now being updated with new studies developed during the last decade. There are also other individual sites with Cannabis/Humulus pollen that need more bibliographic research. A number of these studies are not easy to locate, as they are available only in dissertations and local journals. Finally, some studies do not include Cannabis/Humulus pollen in the diagrams due to its scarcity, but the authors have data in their counting sheets and they can be recovered. Therefore, personal contact with palynologists working on the IP is also needed. All these information sources, along with others that may be located further, should be included in a thorough IP-wide study.

7. Conclusions

  • Currently, it is accepted that Cannabis is a monospecific genus with two subspecies: C. sativa subsp. sativa and C. sativa subsp. indica. However, artificial (human-made) selection has created a multitude of varieties for a diversity of uses, which have been grouped into five biotypes, namely narrow-leaf hemp (NLH), broad-leaf hemp (BLH), narrow-leaf drug (NLD) and broad-leaf drug (BLD), along with a hybrid between the latter two (NLD/BLD).
  • According to DNA phylogenetics calibrated with fossils (molecular clock), Cannabis diverged from its sister genus Humulus during the mid-Oligocene (27.8 Ma). Other similar studies using Bayesian calibration estimate the age of divergence between these two genera in 21 Ma (Early Miocene).
  • The oldest known fossil pollen record compatible with Cannabis was found in 19.6 Ma-old (Early Miocene) rocks from the north-eastern Tibetan Plateau (presently China), which has been proposed as the center of origin of Cannabis.
  • Fossil evidence also shows that wild Cannabis expanded from its center of origin to Europe and eastern Asia between the Miocene and the Pliocene, before the appearance of the genus Homo.
  • During the Pleistocene (the last 2.6 Ma), glacial-interglacial cycles would have fostered recurrent contractions and expansion of Cannabis range leading to the geographical isolation of European and Asian populations. This would have facilitated the differentiation of the extant subspecies: C. sativa subsp. sativa (Europe) and C. sativa subsp. indica (Asia). The first is considered the putative hemp ancestor (PHA), and the second is the putative drug ancestor (PDA).
  • The time and place of Cannabis domestication is still debated between those who propose a single domestication center in central or eastern Asia and those who believe in the possibility of two domestication centers, situated in Asia (PDA) and the Caucasus region (PHA).
  • Recent DNA phylogenetic studies suggest that the early ancestors of hemp and drug types diverged from wild Cannabis in eastern Asia ca. 12 kyr BP (early Neolithic), which situate this plant among the earliest human domesticates.
  • Cannabis was transported outside Eurasia in its cultivated forms and reached Africa (NLD) only after 2000 yr BP. The diffusion to the Americas (NLH and NLD) did not occur until 1545 CE and, by 1945 CE, the Cannabis biotypes were already widespread worldwide.
  • Recent meta-analyses using fossil pollen suggest that wild Cannabis was already widespread across Europe during the Pleistocene and, therefore, the plant could also have been domesticated in this continent. This domestication would have occurred during the late Neolithic-Copper Age, with further expansions during the Bronze-Iron ages.
  • The same meta-analyses indicate that Cannabis could have reached the Iberian Peninsula in its wild form in post-glacial times (18.5 kyr BP) and that cultivation and hemp retting date from the Middle Ages (ca. 600 CE). However, these conclusions are based on only a few localities and consider mostly presence/absence evidence.
  • It is proposed that further meta-analyses in the Iberian Peninsula are based on thorough databases and consider not only the presence but also the abundance of Cannabis pollen, which may be relatively abundant in sites without local Cannabis crops because of its anemophilous character. Some clues for further research are provided

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgements

No specific funding was received for the development of this paper. The comments by John McPartland and two anonymous reviewers contributed to improve the manuscript.

References

Cited by (14)

  • Historical biogeography of Cannabis in the Iberian Peninsula: A probabilistic approach using palynological evidence

    2023, Perspectives in Plant Ecology, Evolution and Systematics
    Citation Excerpt :

    The fossil record of Cannabis is largely based on pollen, which is the most extensively used evidence to reconstruct the historical biogeography of this plant and to trace its cultural relationships. However, the morphological similarity between Cannabis and Humulus pollen has been a handicap for the correct identification of these genera in past records, which is a critical point as the species of these genera have very different ecological requirements and cultural connotations (Portillo et al., 2020; Rull, 2022). This is why different authors refer to the fossil/subfossil representatives of this pollen type using unspecific names such as Cannabis/Humulus (thereafter C/H), Cannabis-type, Humulus-type or Cannabaceae, among others (McPartland et al., 2018).

  • Bacterial phylogenetic markers in lake sediments provide direct evidence for historical hemp retting

    2022, Quaternary Science Reviews
    Citation Excerpt :

    Hemp pollen remained above 20% during the whole retting phases, suggesting that this is the pollen threshold to infer hemp retting in this lake, which falls within the range of previous studies from other areas (15–25%). The results of this study can be useful to interpret paleoecological sequences from other similar lakes in the way toward a more regional view of the introduction, spreading, uses and associated cultural connotations of Cannabis in the Iberian Peninsula within European and Mediterranean contexts (Rull, 2022). VR: conceptualization, methodology, formal analysis, investigation, data curation, writing, visualization.

  • Cannabis de-domestication and invasion risk

    2022, Biological Conservation
    Citation Excerpt :

    This has greatly influenced the taxonomy of the species (Kudo et al., 2009). Cannabis is now generally accepted to be a monotypic genus with genetic differences at the rank of subspecies C. sativa subsp. sativa and C. sativa subsp. indica, although the systematics and nomenclature of the plant is highly debated (Barcaccia et al., 2020; Long et al., 2017; McPartland, 2018; McPartland and Small, 2020; Rull, 2022; Sawler et al., 2015; Serrano-Serrano Martha et al., 2021). A recent whole-genome resequencing study by Ren et al. (2021) identified well-separated genetic clusters which can be divided into four groups: basal cannabis (feral plants collected from China and the US [likely originating from Chinese landraces]), hemp-type (distributed globally), drug-type (distributed globally), and drug-type feral (collected in China, India, and Pakistan).

  • A Bibliometric Analysis of Cannabis-Related Research from 2010 To 2022

    2024, Palestinian Medical and Pharmaceutical Journal

View all citing articles on Scopus


Leave a Reply