09.12.2012 Views

RAPD-based assessment of genetic diversity among annual ...

RAPD-based assessment of genetic diversity among annual ...

RAPD-based assessment of genetic diversity among annual ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

EurAsian Journal <strong>of</strong> BioSciences<br />

Eurasia J Biosci 5, 37-47 (2011)<br />

DOI:10.5053/ejobios.2011.5.0.4<br />

<strong>RAPD</strong>-<strong>based</strong> <strong>assessment</strong> <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> <strong>among</strong><br />

<strong>annual</strong> caraway (Carum carvi) populations<br />

Bochra Laribi1*, Nejia Zoghlami2, Myriam Lamine2, Karima Kouki1, Abdelwahed Ghorbel2, Abdelaziz Mougou1 1Department <strong>of</strong> Agronomy and Vegetal Biotechnology, National Agronomic Institute<br />

<strong>of</strong> Tunisia, 1082 Tunis, Tunisia<br />

2Laboratory <strong>of</strong> Plant Molecular Biology, Center <strong>of</strong> Biotechnology <strong>of</strong> Borj Cedria,<br />

BP 901, 2050 Hammam-Lif, Tunisia<br />

*Corresponding Author: bochra_laribi@yahoo.fr<br />

Abstract<br />

For the first time, <strong>genetic</strong> variability and differentiation <strong>among</strong> five <strong>annual</strong> caraway (Carum<br />

carvi) populations originating from Tunisia, Germany and Egypt were examined. Random<br />

Amplified Polymorphic DNA (<strong>RAPD</strong>) marker data were obtained and analysed with respect to<br />

<strong>genetic</strong> <strong>diversity</strong>, population structure and gene flow. Fourteen primers generated a total <strong>of</strong><br />

136 discernible and reproducible bands across the analyzed populations, out <strong>of</strong> which 56 were<br />

polymorphic. The UPGMA cluster analysis permitted the discrimination <strong>of</strong> all the genotypes and<br />

their sorting into 3 main groups. Tunisian caraway populations diverged significantly from<br />

German and Egyptian ones. Population clustering was made dependently from geographic<br />

origin. This has been further explained at the DNA level as we were able to select a set <strong>of</strong><br />

<strong>RAPD</strong> fingerprints unique to each <strong>of</strong> the studied populations. Furthermore, dimensional graph<br />

derived from factorial analysis <strong>of</strong> <strong>RAPD</strong> frequency data, allowed significant grouping <strong>of</strong> the<br />

genotypes into five sub-plots, representing each one population. Shannon's index values<br />

showed that variation ranks between rather than within populations. These results indicated<br />

that considerable <strong>genetic</strong> differences <strong>among</strong> C. carvi populations were registered.<br />

Keywords: Carum carvi, differentiation, <strong>diversity</strong>, population.<br />

Laribi B, Zoghlami N, Lamine M, Kouki K, Ghorbel A, Mougou A (2011) <strong>RAPD</strong>-<strong>based</strong><br />

<strong>assessment</strong> <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> <strong>among</strong> <strong>annual</strong> caraway (Carum carvi) populations. Eurasia J<br />

Biosci 5: 37-47.<br />

DOI:10.5053/ejobios.2011.5.0.5<br />

INTRODUCTION<br />

Caraway (Carum carvi L.) is one <strong>of</strong> the<br />

oldest aromatic and medicinal plants <strong>of</strong> the<br />

Apiaceae family which is native to Asia,<br />

Europe and North Africa (Nemeth 1998). It is<br />

mainly cultivated in the Netherlands, Finland,<br />

Hungary, Morocco, Iran, India and Russia<br />

(Hornok 1986, Toxopeus and Bowmeester<br />

1993). Besides it is being used as a<br />

condiment, caraway has been utilized for a<br />

long time in traditional medicine especially for<br />

the treatment <strong>of</strong> digestive disorders and was<br />

commonly used in phytomedicine as diuretic<br />

(Lahlou et al. 2007), anti-hyperglycaemic<br />

(Eddouks et al. 2004, Ene et al. 2007,<br />

Tahraoui et al. 2007), anti-hyperchlolesterol<br />

(Lemhadri et al. 2006) and anti-cancerous<br />

(Naderi-Kalali et al. 2005, Kamaleeswari et al.<br />

2006).<br />

Little is known about caraway <strong>genetic</strong><br />

<strong>diversity</strong> and population structure despite its<br />

importance as an aromatic and medicinal<br />

plant. Only a few reports <strong>based</strong> on the species<br />

<strong>of</strong> Lithuania are available and relevant to the<br />

<strong>diversity</strong> and ex situ stability <strong>of</strong><br />

morphological, biochemical and productivity<br />

traits <strong>of</strong> C. carvi cenopopulations growing<br />

wild in natural habitats such as grasslands<br />

(Petraityte 2003, Petraityte and Dastikaite<br />

2007, Dabkevicius et al. 2008). Thus, with<br />

the decline <strong>of</strong> these lasts, mainly by<br />

anthropogenic activities, the C. carvi<br />

bio<strong>diversity</strong> and <strong>genetic</strong> resources become<br />

sparse (Petraityte 2003, Petraityte and<br />

Dastikaite 2007).<br />

Received: February 2011<br />

Received in revised form: May 2011<br />

Accepted: June 2011<br />

Printed: June 2011<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011 37


EurAsian Journal <strong>of</strong> BioSciences<br />

During the last few years, the<br />

characterization and evaluation <strong>of</strong> <strong>genetic</strong><br />

<strong>diversity</strong> and relationships within and between<br />

species and populations were performed<br />

generally by molecular techniques that<br />

substituted the classic ones such as<br />

morphological and physiological characters<br />

(Dolezalova et al. 2003).<br />

Indeed, the nucleotide sequences <strong>of</strong> the<br />

internal transcribed spacer region were used<br />

as molecular markers in order to establish<br />

phylo<strong>genetic</strong> relationships <strong>of</strong> some Apiaceae<br />

genera like Chamaesciadium C. A. Meyer<br />

(Papini 2006) and Carum L. (Papini et al.<br />

2007). It have been demonstrated that this<br />

last genus is polyphyletic. Indeed, maximum<br />

parsimony with bootstrap resampling,<br />

maximum likelihood and Bayesian inference<br />

analyses resulted in three distinct clades for<br />

this genus (Papini et al. 2007).<br />

In addition, the <strong>RAPD</strong> remains one <strong>of</strong> the<br />

most extensively used molecular techniques<br />

due to its simplicity, low cost and high speed.<br />

Thus, <strong>RAPD</strong> markers have been successfully<br />

used in many crops in providing a convenient<br />

and rapid <strong>assessment</strong> <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong><br />

<strong>among</strong> different genotypes (Williams et al.<br />

1990, Rafalski and Tingey 1993, Ragot and<br />

Hoisington 1993). In general, <strong>RAPD</strong> can<br />

provide valuable data in the analysis <strong>of</strong><br />

population <strong>genetic</strong> structure including <strong>genetic</strong><br />

<strong>diversity</strong> within and <strong>among</strong> populations,<br />

population subdivision, and degree <strong>of</strong><br />

inbreeding and individual relatedness (Lunch<br />

and Milligan 1994).<br />

In case <strong>of</strong> aromatic and medicinal plants,<br />

<strong>RAPD</strong> markers have been scarcely used in<br />

<strong>genetic</strong> <strong>diversity</strong> studies, e.g. in Origanum<br />

majorana L. (Klocke et al. 2002), Cunila<br />

galioides Benth (Fracaro et al. 2005),<br />

Curcuma zedoaria (Christm.) Rosc. (Islam et<br />

al. 2007, Syankumar and Sasikuma 2007),<br />

Matricaria chamomilla L. (Solouki et al. 2008),<br />

Satureja hortensis L. (Hadian et al. 2008) and<br />

Foeniculum vulgare Mill. (Zahid et al. 2009).<br />

However, molecular data regarding Tunisian<br />

caraway <strong>genetic</strong> resources still have been<br />

missing. Thus, in this study we present the<br />

first survey <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> <strong>among</strong> five C.<br />

carvi populations; three from Tunisia, one<br />

38<br />

Laribi et al.<br />

from Germany and the other from Egypt,<br />

detected with molecular <strong>RAPD</strong> markers. This<br />

work aims at investigating the extent <strong>of</strong> gene<br />

flow within and between these populations <strong>of</strong><br />

C. carvi as a first step to provide baseline<br />

information for the development <strong>of</strong> strategies<br />

for the conservation and breeding <strong>of</strong> this<br />

aromatic and medicinal plant.<br />

MATERIAL AND METHODS<br />

Plant material and growth conditions<br />

Five Carum carvi L. populations <strong>of</strong><br />

cultivated origin were used in this study.<br />

Three populations were collected from the<br />

field in the principal production regions <strong>of</strong><br />

caraway in Tunisia which are Haouaria,<br />

Menzel Temime and Souassi. The two others<br />

were imported, one from Germany and the<br />

other from Egypt (Table 1).<br />

Seeds were sown in pots <strong>of</strong> 1.5 L volume<br />

filled by agricultural soil containing 0.22,<br />

0.34, 0.05 and 0.08 mequiv. (100 g -1) <strong>of</strong> dry<br />

soil <strong>of</strong> Na +, K +, Ca 2+, Fe 2+, respectively. The<br />

pots were kept in a greenhouse at day light<br />

(photoperiod varying from 13 to 16 h) and at<br />

a temperature varying between 18 and 20°C<br />

during the day and between 10 and 12°C<br />

during the night at the Centre <strong>of</strong><br />

Biotechnology <strong>of</strong> Borj Cedria (36°41'47 N,<br />

10°29'30 E, 26 m above sea level). Preirrigation<br />

was done immediately after sowing<br />

and the moisture content <strong>of</strong> the soil was<br />

controlled by regular irrigation intervals<br />

according to weather conditions for uniform<br />

emergence and seedlings establishment which<br />

took place approximately four weeks after<br />

sowing. After that, fresh leaves were sampled<br />

separately from 10 plants per population on<br />

the basis <strong>of</strong> one plant per pot, and were<br />

immediately ground to a powder in liquid<br />

nitrogen and then stored at -80°C.<br />

DNA extraction<br />

Total foliar DNA was extracted following<br />

the protocol <strong>of</strong> Bowers et al. (1993) later<br />

modified by This et al. (1997) and Zoghlami et<br />

al. (2001). Extracted DNA was quantified by<br />

visual comparison with a molecular size<br />

marker (1 kb ladder) on ethidim bromide<br />

stained agarose gels.<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011


EurAsian Journal <strong>of</strong> BioSciences<br />

DNA amplification<br />

Seventeen universal decamer oligonucleotides,<br />

purchased from the University <strong>of</strong><br />

British Colombia were used for the<br />

polymerase chain reaction (PCR) amplification<br />

(Table 2). They were tested on three<br />

populations for their ability to produce<br />

polymorphic, unambiguous and stable <strong>RAPD</strong><br />

markers (Table 2).<br />

DNA amplification reactions and electrophoresis<br />

were performed according to Zoghlami<br />

et al. (2007). Amplifications were repeated at<br />

least twice and only reproducible products<br />

were taken into account for further data<br />

analysis.<br />

Data analysis<br />

Polymorphic DNA bands were screened<br />

and photographed by a Biometra Bio-doc-IITM<br />

system (Biometra, Göttingen, Germany).<br />

Bands were carefully selected from replicated<br />

amplifications in order to ensure the absence<br />

<strong>of</strong> artefacts and were designated by their<br />

primer code and their size in base pairs.<br />

Polymorphic DNA bands were recorded as<br />

discrete variables for the presence (1) or for<br />

the absence (0) <strong>of</strong> a similar band.<br />

For each primer, the number <strong>of</strong> bands and<br />

the polymorphic ones were calculated as well<br />

as the percentage <strong>of</strong> polymorphic bands<br />

(PPB). The latter was determined as the<br />

percentage <strong>of</strong> polymorphic bands over the<br />

total number <strong>of</strong> the yielded bands. The<br />

number <strong>of</strong> <strong>RAPD</strong> banding pr<strong>of</strong>iles (pr<strong>of</strong>iles<br />

generated by all the populations per primer)<br />

has also been calculated. The ability <strong>of</strong> the<br />

most informative primers to differentiate<br />

between populations was assessed by the<br />

estimation <strong>of</strong> their resolving power (Rp)<br />

(Prevost and Wilkinson 1999). The Rp has<br />

been described to correlate strongly with the<br />

ability to distinguish between the populations<br />

according to the following Gilbert et al.<br />

(1999) formula:<br />

Rp= Ib<br />

Where Ib= 1- ( ) where p is the<br />

population proportion containing the I band.<br />

<strong>RAPD</strong> bands were transformed into a<br />

binary character matrix. Populations were<br />

then clustered using the Unweighed Pair<br />

Group Method with Arithmetic averaging<br />

Laribi et al.<br />

algorithm (UPGMA, Sneath and Sokal 1973)<br />

using Darwin s<strong>of</strong>tware (Version 5.0.148,<br />

http://darwin.cirad.fr/darwin). The origin <strong>of</strong><br />

populations under investigation was taken<br />

into account in order to examine their<br />

potential effect on the <strong>genetic</strong> clustering. We<br />

also investigated <strong>genetic</strong> variation <strong>among</strong><br />

populations using Factorial Correspondence<br />

Analysis (FCA).<br />

Population differentiation analysis has been<br />

assessed <strong>among</strong> and within the studied<br />

populations according to Nei (1978) using the<br />

computer program POPGENE 1.32 (Yeh et al.<br />

1999). Shannon's index <strong>of</strong> genotypic <strong>diversity</strong><br />

(H0) for <strong>RAPD</strong> was estimated to quantify the<br />

degree <strong>of</strong> within-population <strong>diversity</strong> following<br />

the calculation:<br />

H0=- piln pi ; where pi was the frequency<br />

<strong>of</strong> presence or absence <strong>of</strong> a <strong>RAPD</strong> band in a<br />

population, as described by Yeh et al. (1995).<br />

The average <strong>diversity</strong> over all populations<br />

(Hpop) was calculated as:<br />

1<br />

H H<br />

pop= o<br />

n ; where n was the number <strong>of</strong><br />

populations Yeh et al. (1995).<br />

The mean <strong>diversity</strong> at species level (Hsp) was calculated as:<br />

Hsp=- psln ps; where ps was the frequency<br />

<strong>of</strong> presence or absence <strong>of</strong> a <strong>RAPD</strong><br />

band in all the populations.<br />

The component <strong>of</strong> within-population<br />

<strong>diversity</strong> was estimated as Hpop/Hsp, and that<br />

<strong>of</strong> between-population <strong>diversity</strong> as 1-Hpop/Hsp. The gene flow (Nm, number <strong>of</strong> migrants per<br />

generation) (McDermott and McDonald 1993)<br />

was approximated as:<br />

Nm= 0.5*(1-GST)/GST; where GST (analogous<br />

to FST) (Wright 1951) was the <strong>diversity</strong><br />

coefficient <strong>among</strong> populations.<br />

RESULTS<br />

Genetic polymorphism<br />

A total <strong>of</strong> 17 primers were screened for<br />

their ability to generate consistently amplified<br />

band patterns and to assess polymorphism in<br />

the tested populations (Table 2). One hundred<br />

and thirty six (136) bands were obtained,<br />

which is an average <strong>of</strong> 9.7 bands per primer:<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011 39


EurAsian Journal <strong>of</strong> BioSciences<br />

Table 1. Origin <strong>of</strong> C. carvi populations included in<br />

this study with their respective geographic<br />

coordinates and altitude.<br />

80 bands were common in all populations and<br />

56 were polymorphic (Table 2).<br />

Except primers UBC 204, UBC 235 and<br />

UBC 248, all primers generated multiple<br />

banding patterns with 1 to 8 polymorphic<br />

amplified DNA bands ranging in size from 250<br />

to 3000 pb. The polymorphic markers yielded<br />

44 different electrophoretic banding pr<strong>of</strong>iles.<br />

An example <strong>of</strong> <strong>RAPD</strong> banding pr<strong>of</strong>ile<br />

generated by the primers UBC 162, UBC 213<br />

and UBC 225 is presented in Fig.1. Hence, we<br />

may assume that a large <strong>genetic</strong> <strong>diversity</strong> at<br />

the DNA level characterises the five <strong>annual</strong><br />

caraway populations.<br />

Estimation <strong>of</strong> Rp values exhibited a<br />

collective rate <strong>of</strong> 30 and varied from 0.4 for<br />

UBC 211 and UBC 231 primers to 5.2 for<br />

UBC 264 one with a mean <strong>of</strong> 2.14 (Table 2).<br />

UBC 261 and UBC 264 primers were the most<br />

informative primers as displaying the highest<br />

Rp rates. Thus, these latter could be referred<br />

to as the most efficient loci in assessing<br />

<strong>genetic</strong> <strong>diversity</strong> within C. carvi.<br />

The distinctiveness <strong>of</strong> the clusters<br />

identified in the UPGMA derived dendrogram<br />

exhibited 3 distinct groups (A, B and C) (Fig.<br />

2). Hence, <strong>genetic</strong> divergence was dependent<br />

from the geographic origin and all tested<br />

populations were differentiated.<br />

German and Egyptian genotypes housed in<br />

Group A were referred to as an out-put <strong>of</strong><br />

analysis as diverging significantly from<br />

Tunisian populations. This has been further<br />

explained at the DNA level as we were able to<br />

select a set <strong>of</strong> <strong>RAPD</strong> fingerprints unique to<br />

each <strong>of</strong> the studied populations (Table 2). For<br />

example, we counted 5 unique <strong>RAPD</strong> bands<br />

for the Tunisian population generated by the<br />

primers UBC 226, UBC 225, UBC 243, UBC<br />

261 and UBC 264 and sizing respectively<br />

700, 835, 1700, 1100 and 975 base pairs.<br />

40<br />

Laribi et al.<br />

The two-dimensional FCA graph allowed<br />

also significant grouping <strong>of</strong> the populations<br />

that were plotted into five sub-plots,<br />

representing each one population (Fig. 3).<br />

Population analysis<br />

Primers varied in their ability to detect<br />

variation at both within and between<br />

populations, with <strong>among</strong> populations H0 ranging from 8.5% (primer UBC-226) to<br />

100% (primers UBC-162 and UBC-231)<br />

(Table 2). On average over all primers, the<br />

population Haouaria (1) exhibited the lowest<br />

level <strong>of</strong> within population <strong>genetic</strong> <strong>diversity</strong><br />

(mean H0 0.012), while the populations<br />

Menzel Temime (2), Germany (4) and Egypt<br />

(5) displayed similar estimates <strong>of</strong> mean H0 (0.022-0.031) (Table 2). However, population<br />

Souassi (3) was the most variable (mean H0 0.107) (Table 2).<br />

The results derived from Shannon's index<br />

showed that variation within (Hpop/Hsp) and<br />

between (1- Hpop/Hsp) populations was 12.9<br />

and 87.1%, respectively (Table 2). Thus, high<br />

level <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> was detected<br />

between populations rather than within<br />

populations.<br />

A matrix <strong>of</strong> pair-wise GST values, the<br />

significance and the effective number <strong>of</strong><br />

migrants (Nm) between populations are<br />

presented in Table 3. Values <strong>of</strong> GST ranged<br />

from 0.288 (between populations Haouaria<br />

and Menzel Temime) to 0.765 (between<br />

populations Souassi and Germany), averaging<br />

0.671. This indicates that all populations may<br />

be considered different from each other, with<br />

population Germany (4) being the most<br />

different from the others and populations<br />

Haouaria (1) and Menzel Temime (2) being the<br />

most similar.<br />

The significant differentiation between<br />

populations was also revealed in the estimates<br />

<strong>of</strong> gene flow (Nm) (Table 3). Values <strong>of</strong> Nm ranged from a moderate value <strong>of</strong> 0.465 to<br />

0.749, averaging 0.585 which indicated low<br />

gene flow between the five studied<br />

populations <strong>of</strong> C. carvi.<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011


EurAsian Journal <strong>of</strong> BioSciences<br />

In this study, we fingerprinted a set <strong>of</strong> five<br />

<strong>annual</strong> caraway populations originating from<br />

Tunisia, Germany and Egypt by means <strong>of</strong><br />

<strong>RAPD</strong> markers in order to obtain molecular<br />

data on their <strong>genetic</strong> background. Thus, we<br />

have demonstrated the reliability <strong>of</strong> <strong>RAPD</strong><br />

analysis to detect DNA polymorphisms and<br />

relationships within C. carvi populations.<br />

The selected primers were highly<br />

discriminating since they were characterised<br />

by relatively high collective Rp rate (30.00), a<br />

high number <strong>of</strong> polymorphic markers and<br />

electrophoretic banding patterns. The primers<br />

generated 56 polymorphic bands out <strong>of</strong> 136<br />

with a mean <strong>of</strong> 4. Only three primers did not<br />

allow DNA amplification (Table 2). In <strong>RAPD</strong><br />

literature, the presence <strong>of</strong> primers does not<br />

allow amplification to occur (Caetano-Annoles<br />

1994), whereas others primers yielding faint<br />

banding pr<strong>of</strong>iles (Moreno et al. 1995, Ortiz et<br />

al. 1997) were reported. In addition,<br />

according to the results forwarded by Devos<br />

and Gale (1992), Penner et al. (1993) and<br />

This et al. (1997), some primers are more<br />

Laribi et al.<br />

Table 2. Genetic <strong>diversity</strong> analysis through <strong>RAPD</strong> markers and Shannon's index <strong>among</strong> five <strong>annual</strong> C. carvi<br />

populations originated from Tunisia, Germany and Egypt.<br />

*Population labels are given in Table 1.<br />

Table 3. Pair-wise GST values (below diagonal) and<br />

Nm (above diagonal) between the five<br />

<strong>annual</strong> C. carvi populations.<br />

*Population labels are given in Table 1.<br />

DISCUSSION<br />

Fig. 1. An example <strong>of</strong> <strong>RAPD</strong> banding pr<strong>of</strong>ile <strong>of</strong> five<br />

C. carvi genotypes generated using the<br />

primers UBC 162, UBC 213 and UBC 225.<br />

Lane 1: Haouaria, lane 2: Menzel Temime,<br />

lane 3: Souassi, lane 4: Germany, lane 5:<br />

Egypt; lane-M: Molecular size marker (1Kb<br />

DNA ladder). Size <strong>of</strong> DNA fragments is<br />

presented on right.<br />

Fig. 2. UPGMA cluster analysis within the five<br />

<strong>annual</strong> C. carvi populations established by<br />

means <strong>of</strong> 56 <strong>RAPD</strong> markers.<br />

The letters (from A to C) stand for the group<br />

individualized (1: Haouaria, 2: Menzel Temime,<br />

3: Souassi, 4: Germany, 5: Egypt).<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011 41


EurAsian Journal <strong>of</strong> BioSciences<br />

Fig. 3. FCA Two-dimensional plot <strong>of</strong> C. carvi<br />

populations. (1: Haouaria, 2: Menzel Temime,<br />

3: Souassi, 4: Germany, 5: Egypt)<br />

efficient than others in producing stable and<br />

reproducible pr<strong>of</strong>iles.<br />

The <strong>genetic</strong> divergence <strong>of</strong> the populations<br />

under investigation was confirmed at the DNA<br />

level. In fact, the UPGMA cluster analysis<br />

permitted the discrimination <strong>of</strong> all the<br />

genotypes and their sorting into 3 main<br />

groups (Fig. 2). German and Egyptian caraway<br />

populations (Group A, Fig. 2) diverged<br />

significantly from Tunisian ones. Thus,<br />

population clustering was made dependently<br />

from the original region. In the same way, the<br />

occurrence <strong>of</strong> <strong>RAPD</strong> fingerprints unique to<br />

each <strong>of</strong> the studied populations (Table 2)<br />

suggests that <strong>RAPD</strong> markers may constitute<br />

rapid molecular tools for assigning one<br />

genotype to its origin. On the other hand,<br />

<strong>RAPD</strong> markers are suitable to perform <strong>genetic</strong><br />

variation analyses at both intra and interpopulation<br />

levels (Syankumar and Sasikuma<br />

2007, Hadian et al. 2008, Zahid et al. 2009).<br />

According to Shannon's index, high levels<br />

<strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> were detected between C.<br />

carvi populations rather than within<br />

populations (Table 2). The existence <strong>of</strong> low<br />

<strong>genetic</strong> <strong>diversity</strong> within the species has been<br />

mostly attributed to self pollination, unless<br />

other environmental pressures are influencing<br />

<strong>genetic</strong> <strong>diversity</strong> (Archak et al. 2002). In our<br />

case, despite being an out-crossing and<br />

insect-pollinated plant (Nemeth et al. 1999,<br />

Nemeth and Szekely 2000), C. carvi exhibited<br />

low levels <strong>of</strong> <strong>genetic</strong> variation within<br />

populations. In general, C. carvi<br />

cenopopulations are grown wild in natural<br />

habitats such as grasslands and are<br />

characterized by rich <strong>genetic</strong> <strong>diversity</strong><br />

42<br />

Laribi et al.<br />

strongly influenced by ecological conditions<br />

(Petraityte 2003, Petraityte and Dastikaite<br />

2007). Therefore, C. carvi confined to seminatural<br />

grasslands responds to decreases in<br />

fragment size by forming smaller populations<br />

(Kiviniemi 2008). This may have directly<br />

reduced <strong>genetic</strong> <strong>diversity</strong>, and also subjected<br />

the remaining populations to <strong>genetic</strong> drift<br />

(Nunney and Elam 1994). Thus, C. carvi<br />

population persistence is dependent on<br />

continuous management that maintains<br />

habitat quality and creates suitable conditions<br />

for regeneration, and hence opportunities for<br />

positive population growth (Kiviniemi 2008).<br />

Furthermore, another possible explanation<br />

for the low <strong>genetic</strong> <strong>diversity</strong> <strong>of</strong> C. carvi is the<br />

influence <strong>of</strong> human activity. In fact, the recent<br />

rapid development <strong>of</strong> tourism and urbanization<br />

has led to low population densities and<br />

restricted distribution <strong>of</strong> C. carvi<br />

cenopopulations, resulting in loss and<br />

reduction <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong>. As a result, the<br />

disappearance <strong>of</strong> C. carvi cenopopulations<br />

makes the <strong>genetic</strong> resources <strong>of</strong> the species<br />

poor (Petraityte 2003, Petraityte and<br />

Dastikaite 2007) and conventional caraway<br />

breeding is hampered mainly by the limited<br />

available <strong>genetic</strong> variation <strong>of</strong> the caraway<br />

gene pool (Nemeth 1998).<br />

Detecting the patterns <strong>of</strong> variation, in<br />

terms <strong>of</strong> polymorphism and differentiation<br />

<strong>among</strong> C. carvi populations using <strong>RAPD</strong><br />

markers, is similar to the finding obtained by<br />

Fu et al. (2003) for Changium smyrnioides<br />

Wolff (Apiaceae) who also reported that<br />

variation was be partitioned between rather<br />

than within populations. However, these<br />

results are not in accordance with those <strong>of</strong><br />

Barbosa et al. (2010) in a medicinal plant<br />

Palicourea coriacea (Cham.) K. Schum.<br />

(Rubiaceae) for which, 23% <strong>of</strong> <strong>genetic</strong><br />

variability is found <strong>among</strong> populations and<br />

77% <strong>of</strong> <strong>genetic</strong> variability is within<br />

populations.<br />

GST (0.671) and Nm (Nm 0.585, Nm


EurAsian Journal <strong>of</strong> BioSciences<br />

structure <strong>of</strong> C. carvi populations. However,<br />

pollen flow may be inhibited by the distance<br />

between these populations. Indeed, the<br />

caraway species have not been able by<br />

natural means, e.g. by wind or by domestic<br />

cattle, to exploit all suitable habitat sites<br />

within grassland fragments (Kiviniemi 2008).<br />

Furthermore, seed dispersal capacity <strong>of</strong> this<br />

plant may be limited because its seeds lack<br />

morphological adaptations for dispersal.<br />

Consequently, gene flow <strong>among</strong> populations<br />

was expected to be low and <strong>genetic</strong><br />

differentiation to be high.<br />

The same results were obtained by<br />

Aegisdóttir et al. (2009) who analyzed<br />

molecular <strong>diversity</strong> <strong>of</strong> a rare alpine plant<br />

species Campanula thyrsoides L. using five<br />

polymorphic microsatellite loci. They<br />

demonstrated that this species showed high<br />

<strong>genetic</strong> <strong>diversity</strong> and considerable population<br />

differentiation.<br />

However, these results are in contrast to<br />

those reported by Sebastien et al. (2010) who<br />

found low levels <strong>of</strong> differentiation in a<br />

medicinal plant Tylophora rotundifolia Buch.-<br />

Ham. ex Wight using <strong>RAPD</strong> markers, which<br />

they attributed to a common <strong>genetic</strong> pool,<br />

lack <strong>of</strong> differential selection pressure,<br />

diminished mutation rate and short life span <strong>of</strong><br />

the species.<br />

As expected, the groupings <strong>of</strong> populations<br />

were largely in accordance with their<br />

geographical localities, and the representation<br />

<strong>of</strong> <strong>genetic</strong> relatedness <strong>among</strong> individuals<br />

obtained with the two-dimensional graph<br />

derived from factorial analysis <strong>of</strong> <strong>RAPD</strong><br />

REFERENCES<br />

Laribi et al.<br />

frequency data (Fig. 3) also indicated the low<br />

level <strong>of</strong> gene flow <strong>among</strong> C. carvi populations.<br />

CONCLUSION<br />

To conclude, this is the first report on the<br />

<strong>assessment</strong> <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> and<br />

population differentiation analysis in five<br />

<strong>annual</strong> caraway (Carum carvi) populations<br />

originating from Tunisia, Germany and Egypt,<br />

using <strong>RAPD</strong> markers. Though the number <strong>of</strong><br />

populations used in this work does represent<br />

the total <strong>diversity</strong> <strong>of</strong> Tunisian caraway<br />

populations; nor those <strong>of</strong> the Egyptian and the<br />

German populations. We successfully<br />

provided deep insights on the <strong>genetic</strong><br />

background <strong>of</strong> the above mentioned<br />

populations. Considerable <strong>genetic</strong> <strong>diversity</strong><br />

has been detected within the species.<br />

Nevertheless, using more informative markers<br />

such as SSRs and including more populations<br />

would allow detecting the <strong>genetic</strong> background<br />

resources <strong>of</strong> our local caraway. Therefore, to<br />

prevent further substantial loss <strong>of</strong> <strong>genetic</strong><br />

<strong>diversity</strong>, we should conserve as many<br />

populations as possible. In this context, we<br />

suggest that an appropriate strategy for<br />

sampling may be formulated when ex situ<br />

conservation is required. Besides, the best in<br />

situ conservation strategy for C. carvi<br />

population is to preserve its natural habitats.<br />

We further aim to include wild germplasm in<br />

order to select cultivars <strong>of</strong> economically<br />

importance and to transfer important <strong>genetic</strong><br />

traits from wild to cultivated species by<br />

marker-assisted-selection.<br />

Aegisdottir H, Kuss P, Stöcklin J (2009) Isolated populations <strong>of</strong> a rare alpine plant show high<br />

<strong>genetic</strong> <strong>diversity</strong> and considerable population differentiation. Annals <strong>of</strong> Botany 104: 1313-<br />

1322.<br />

Archak S, Karihaloo JL, Jain A (2002) <strong>RAPD</strong> markers reveal narrowing <strong>genetic</strong> base <strong>of</strong> Indian<br />

tomato cultivars. Current Science 82: 1139-1143.<br />

Barbosa TCS, Sibov ST, Telles MPC, Soares TN (2010) Genetic characterization <strong>of</strong> natural<br />

populations <strong>of</strong> the medicinal plant Palicourea coriacea (Rubiaceae) with molecular markers.<br />

Genetics and Molecular Research 9: 695-704.<br />

Bowers JE, Bandmann EB, Meredith CP (1993) DNA fingerprinting and characterization <strong>of</strong> some<br />

wine grape cultivars. American Journal <strong>of</strong> Enology and Viticulture 44: 266-274.<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011 43


EurAsian Journal <strong>of</strong> BioSciences<br />

44<br />

Laribi et al.<br />

Caetano-Annoles G (1994) A versatile and universal tool for genome analysis. Plant Molecular<br />

Biology 25: 1011-1026.<br />

Dabkevicius Z, Gelvonauskis B, Leistrumaite A (2008) Investigation <strong>of</strong> <strong>genetic</strong> resources <strong>of</strong><br />

cultivated plants in Lithuania. Biologija 54: 51-55.<br />

Devos KM, Gale MD (1992) The use <strong>of</strong> random amplified polymorphic DNA markers in wheat.<br />

Theoretical and Applied Genetics 84: 567-572.<br />

Dolezalova I, Lebeda A, Dziechciarkova M, Kristkova E, Astley D, Wiel C (2003) Relationships<br />

<strong>among</strong> morphological characters, isozymes polymorphism and DNA variability-the impact on<br />

lactuca germplasm taxonomy. Czech Journal <strong>of</strong> Genetics and Plant Breeding 39: 59-67.<br />

Eddouks M, Lemhadri A, Michel JB (2004) Caraway and caper: potential anti-hyperglycaemic<br />

plants in diabetic rats. Journal <strong>of</strong> Ethnopharmacology 94: 143-148.<br />

Ene AC, Nwankwo EA, Samdi LM (2007) Alloxan-induced diabetes in rats and the effects <strong>of</strong><br />

black caraway (Carum carvi L.) oil on their body weight. Research Journal <strong>of</strong> Medicine and<br />

Medical Sciences 2: 48-52.<br />

Fracaro F, Zacaria J, Echeverrigaray S (2005) <strong>RAPD</strong> <strong>based</strong> <strong>genetic</strong> relationships between<br />

populations <strong>of</strong> three chemotypes <strong>of</strong> Cunila galioides Benth. Biochemical and Systematics<br />

Ecology 33: 409-417.<br />

Fu C, Yingxiong QIU, Kong H (2003) <strong>RAPD</strong> analysis for <strong>genetic</strong> <strong>diversity</strong> in Changium<br />

smyrnioides (Apiaceae), an endangered plant. Botanical Bulletin <strong>of</strong> Academia Sinica 44: 13-<br />

18.<br />

Gilbert JE, Lewis RV, Wilkinson MJ (1999) Developing an appropriate strategy to assess <strong>genetic</strong><br />

variability in plant germplasm collections. Theoretical and Applied Genetics 98: 1125-1131.<br />

Hadian J, Tabatabaei SMF, Naghavi MR, Jamzad Z, Ramak-Masoumi T (2008) Genetic <strong>diversity</strong><br />

<strong>of</strong> Iranian accessions <strong>of</strong> Satureja hortensis L. <strong>based</strong> on horticultural traits and <strong>RAPD</strong> markers.<br />

Scientia Horticulturae 115: 196-202.<br />

Hornok L (1986) Effect <strong>of</strong> environmental factors on growth, yield and on the active principles<br />

<strong>of</strong> some spice plants. Acta Horticulturae 188: 69-176.<br />

Islam MA, Meister A, Schubert V, Kloppstech K, Esch E (2007) Genetic <strong>diversity</strong> and<br />

cyto<strong>genetic</strong> analyses in Curcuma zedoaria (Christm.) Roscoe from Bangladesh. Genetic<br />

Resources and Crop Evolution 54: 149-156.<br />

Kamaleeswari M, Kumaraswami Deeptha K, Murugan Sengottuvelan M, Nalini N (2006) Effect<br />

<strong>of</strong> dietary caraway (Carum carvi L.) on aberrant crypt foci development, fecal steroids, and<br />

intestinal alkaline phosphatase activities in 1,2-dimethylhydrazine-induced colon<br />

carcinogenesis. Toxicology and Applied Pharmacology 214: 290-296.<br />

Kiviniemi K (2008) Effects <strong>of</strong> fragment size and isolation on the occurrence <strong>of</strong> four short-lived<br />

plants in semi-natural grasslands. Acta Oecologica 33: 56-65.<br />

Klocke E, Langbehn J, Grewe C, Pank F (2002) DNA fingerprinting by <strong>RAPD</strong> on Origanum<br />

majorana. Journal <strong>of</strong> Herbs, Spices and Medicinal Plants 9: 171-176.<br />

Lahlou S, Tahraoui A, Zafar Israili Z, Lyoussi B (2007) Diuretic activity <strong>of</strong> the aqueous extracts<br />

<strong>of</strong> Carum carvi and Tanacetum vulgare in normal rats. Journal <strong>of</strong> Ethnopharmacology 110:<br />

458-463.<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011


EurAsian Journal <strong>of</strong> BioSciences<br />

Laribi et al.<br />

Lemhadri A, Hajji L, Michel J, Eddouks M (2006) Cholesterol and triglycerides lowering activities<br />

<strong>of</strong> caraway fruits in normal and streptozotocin diabetic rats. Journal <strong>of</strong> Ethnopharmacology<br />

106: 321-326.<br />

Lunch M, Milligan BG (1994) Analysis <strong>of</strong> population <strong>genetic</strong> structure with <strong>RAPD</strong> markers.<br />

Molecular Ecology 3: 91-99.<br />

McDermott JM, McDonald BA (1993) Gene flow in plant pathosystems. Annual Review <strong>of</strong><br />

Phytopathology 31: 353-373.<br />

Moreno S, Gogorcena Y, Ortiz JM (1995) The use <strong>of</strong> <strong>RAPD</strong> markers for identification <strong>of</strong><br />

cultivated grapevine (Vitis vinifera L.). Scientia Horticulturae 62: 237-243.<br />

Naderi-Kalali B, Allameh A, Rasaee MJ, Bach HJ, Behechti A, Doods K, Kettrup A, Schramm<br />

KW (2005) Suppressive effects <strong>of</strong> caraway (Carum carvi) extracts on 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin-dependent<br />

gene expression <strong>of</strong> cytochrome P450 1A1 in the rat H4IIE cells.<br />

Toxicology in Vitro 19: 373-377.<br />

Nei M (1978) Estimation <strong>of</strong> average heterozygosity and <strong>genetic</strong> distance from a small number<br />

<strong>of</strong> individuals. Genetics 89: 583-590.<br />

Nemeth E (1998) Caraway: The Genus Carum. Harwood Academic Publishers, Amsterdam.<br />

Nemeth E, Bernath J, Petheo F (1999) Study on flowering dynamic and fertilization properties<br />

<strong>of</strong> caraway and fennel. Acta Horticulturae 502: 77-83.<br />

Nemeth E, Szekely G (2000) Floral biology <strong>of</strong> medicinal plants. I. Apiaceae species. International<br />

Journal Horticultural Science 6: 133-136.<br />

Nunney L, Elam DR (1994) Estimating the effective population size <strong>of</strong> conserved populations.<br />

Conservation Biology 8: 175-184.<br />

Ortiz A, Renaud R, Calzada I, Ritter E (1997) Analysis <strong>of</strong> Palm cultivars with <strong>RAPD</strong> markers.<br />

Journal <strong>of</strong> Horticultural Science 72: 1-9.<br />

Papini A (2006) The systematic position <strong>of</strong> Chamaesciadium C. A. Meyer (Umbelliferae) on the<br />

basis <strong>of</strong> nuclear ITS sequence. Flora Mediterranea 16: 5-15.<br />

Papini A, Banci F, Nardi E (2007) Molecular evidence that genus Carum L. (Apiaceae) is<br />

polyphyletic. Genetics and Molecular Biology 30(2): 475-482.<br />

Penner GA, Bush A, Wise R, Kim W, Domier L, Kasha K, Laroche A, Scoles G, Molnar SJ, Fedak<br />

G (1993) Reproducibility <strong>of</strong> random amplified polymorphic DNA (<strong>RAPD</strong>) analysis <strong>among</strong><br />

laboratories. PCR Methods and Applications 2: 341-345.<br />

Petraityte N (2003) Ex-situ stability <strong>of</strong> morpho-biochemical properties <strong>of</strong> common caraway<br />

(Carum carvi L.). Ekologija (Vilnius) 1: 43-46.<br />

Petraityte N, Dastikaite A (2007) Agrobiological <strong>assessment</strong> <strong>of</strong> wild Carum carvi L.<br />

cenopopulation bio<strong>diversity</strong> ex situ. Biologija 53: 74-79.<br />

Prevost A, Wilkinson MJ (1999) A new system <strong>of</strong> comparing PCR primers applied to ISSR<br />

fingerprinting <strong>of</strong> potato cultivars. Theoretical and Applied Genetics 98:107-112.<br />

Rafalski JA, Tingey SV (1993) Genetic diagnosis in plant breeding: <strong>RAPD</strong>s microsatellites and<br />

machines. Trends in Genetics 9: 275-280.<br />

Ragot M, Hoisington DA (1993) Molecular markers for plant breeding: comparisons <strong>of</strong> RFLP and<br />

<strong>RAPD</strong> genotyping costs. Theoretical and Applied Genetics 86: 975-984.<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011 45


EurAsian Journal <strong>of</strong> BioSciences<br />

46<br />

Laribi et al.<br />

Sebastian VA, D'Cruz L, Subramanian RB, Braganza VJ (2010) Assessment <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong><br />

within and <strong>among</strong> populations <strong>of</strong> Tylophora rotundifolia using <strong>RAPD</strong> markers. Gene conserve<br />

9: 94-117.<br />

Sneath PHA, Sokal RR (1973) Numerical Taxonomy. The Principle and Practice <strong>of</strong> Numerical<br />

Classification. W.H. Freeman and Co., San Francisco.<br />

Solouki M, Mehdikhani H, Zeinali H, Emamjomeh AA (2008) Study <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> in<br />

Chamomile (Matricaria chamomilla) <strong>based</strong> on morphological traits and molecular markers.<br />

Scientia Horticulturae 117: 281-287.<br />

Syamkumar S, Sasikumar B (2007) Molecular marker <strong>based</strong> <strong>genetic</strong> <strong>diversity</strong> analysis <strong>of</strong><br />

Curcuma species from India. Scientia Horticulturae 112: 235-241.<br />

Tahraoui A, El-Hilaly J, Israili ZH, Lyoussi B (2007) Ethnopharmacological survey <strong>of</strong> plants used<br />

in the traditional treatment <strong>of</strong> hypertension and diabetes in south-eastern Morocco (Errachidia<br />

province). Journal <strong>of</strong> Ethnopharmacology 110: 105-117.<br />

This P, Cuisset C, Boursiquot JM (1997) Development <strong>of</strong> stable <strong>RAPD</strong> markers for the<br />

identification <strong>of</strong> grapevine rootstocks and the analysis <strong>of</strong> <strong>genetic</strong> relationships. American<br />

Journal <strong>of</strong> Enology and Viticulture 48: 492-501.<br />

Toxopeus H, Bouwmeester HJ (1993) Improvement <strong>of</strong> caraway essential oil and carvone<br />

production in The Netherlands. Industrial Crops and Products 1: 295-301.<br />

Williams JGK, Kubelik AR, Livak KJ, Rafalski JA (1990) DNA polymorphisms amplified as<br />

arbitrary primers are useful <strong>genetic</strong> markers. Nucleic Acids Research 18: 6531-6535.<br />

Wright S (1951) The <strong>genetic</strong>al structure <strong>of</strong> populations. Annals <strong>of</strong> Eugenics 15: 323-354.<br />

Yeh FC, Chong DKX, Yang RC (1995) <strong>RAPD</strong> variation within and <strong>among</strong> natural populations <strong>of</strong><br />

trembling aspen (Populus tremuloides) from Alberta. Journal <strong>of</strong> Heredity 86: 454-460.<br />

Yeh FC, Yang R, Boyle T (1999) POPGENE Version 1.32. Micros<strong>of</strong>t Window-Based Freeware<br />

for Population Genetic Analysis. University <strong>of</strong> Alberta, Edmonton.<br />

Zahid NY, Abbasi NA, Hafiz IA, Ahmad Z (2009) Genetic <strong>diversity</strong> <strong>of</strong> indigenous fennel<br />

(Foeniculum vulgare Mill.) germoplasm in Pakistan assessed by <strong>RAPD</strong> markers. Pakistan<br />

Journal <strong>of</strong> Botany 41: 1759-1767.<br />

Zoghlami N, Mliki A, Ghorbel A (2001) Evaluation <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> <strong>among</strong> Tunisian<br />

grapevines by <strong>RAPD</strong> markers. Vitis 40: 31-37.<br />

Zoghlami N, Chrita I, Bouamama B, Gargouri M, Zemni H, Ghorbel A, Mliki A (2007) Molecular<br />

<strong>based</strong> <strong>assessment</strong> <strong>of</strong> <strong>genetic</strong> <strong>diversity</strong> within Barbary fig (Opuntia ficus indica (L.) Mill.) in<br />

Tunisia. Scientia Horticulturae 113: 134-141.<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011


EurAsian Journal <strong>of</strong> BioSciences<br />

Laribi et al.<br />

Frenk Kimyonu (Carum carvi) Popülasyonlarinda Genetik Çesitliligin <strong>RAPD</strong> Temelli Incelenmesi<br />

Özet<br />

Tunus, Almanya ve Misir kaynakli yillik Frenk kimyonuna (Carum carvi) ait bes popülasyon<br />

arasindaki genetik çesitlilik ve farklilasma, ilk kez incelendi. Genetik çesitlilik, popülasyon yapisi ve<br />

gen akisi ile alakali Random Amplified Polymorphic DNA (<strong>RAPD</strong>) markör verileri elde edildi ve<br />

incelendi. Analiz edilen popülasyonlarda; 14 primer, toplam 136 ayird edilebilir ve tekrarlanabilir<br />

bant üretti. Bunlardan 56 tanesi polimorfik idi. UPGMA yigin analizi bütün genotipin ayirt<br />

edilmesini ve 3 ana grupta toplanmasini sagladi. Tunus kimyon popülasyonlari, anlamli ölçüde<br />

Alman ve Misir popülasyonlarindan ayrildi. Popülasyon gruplanmasi, cografik kaynaktan bagimsiz<br />

olarak yapildi. Çalisilan popülasyonlara özel <strong>RAPD</strong> izleri seçebildigimiz için, DNA seviyesinde de bu<br />

durumu göstermis olduk. Ek olarak, <strong>RAPD</strong> frekans verilerinin analizinden elde edilen boyutsal<br />

grafik; genotipleri, her biri popülasyonu temsil eden bes alt dal halinde gruplandirmamizi sagladi.<br />

Shannon indeks degerleri; farkliligin popülasyonlarin kendi içinde degil, popülasyonlar arasinda<br />

gerçeklestigini gösterdi. Bu sonuçlar, C. carvi populasyonlari arasinda önemli genetik farkliliklar<br />

oldugunu göstermistir.<br />

Anahtar Kelimeler: Carum, çesitlilik, farklilik, popülasyon.<br />

©EurAsian Journal <strong>of</strong> BioSciences, 2011 47

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!