Echinochloa
(L.) P. Beauv. includes some of the noxious weeds, causing a serious yield loss when they are dominant in the fields. Identification of the
Echinochloa
is very difficult because many interspecific and intraspecific forms of the species are found. However, it is important to identify the species exactly and to know the genetic diversity of the species for effective weed management. This study was conducted to identify and summarize the
Echinochloa
species by comparing the genetic variation and relationship among Korean
Echinochloa
species using SSR. The genetic diversity of 107 individuals, including seven species were assessed using five SSR markers. UPGMA dendrogram generated two clades (I and II) and clade II divided again into two subclades (II-1 and II-2) whereas the model based genetic structure proposed four subpopulations. The two subpopulations were corresponded to clades I and II-1 and the other two were arranged to clade II-2 of the UPGMA dendrogram. We have concluded that
E. colona
and
E. glabrescens
might have not distributed in Korea. The biological varieties,
praticola
and
echinata
, of
E. crus-galli
should be treated as
E. crus-galli
. Korean
Echinochloa
should be summarized with four species, i.e.,
E. oryzicola, E. crus-galli, E. esculenta
, and
E. oryzoides
.
Introduction
The genus
Echinochloa
(L.) P. Beauv. is a very difficult genus due to many interspecific and intraspecific forms found in the nature. This genus is still taxonomically controversial due to the absence of reliable morphological characters, resulting in poor taxonomic understanding of the genus. It has approximately 30- 40 species, mainly distributed in the tropical and warm temperate regions, usually associated with wet or damp places, of the world (
Clayton and Renvoize, 1999
).
This genus includes some of noxious weeds, such as
E. crusgalli
(L.) P. Beauv.,
E. oryzicola
Vasinger, and
E. colona
(L.) Link in agricultural areas of the world. The dominance of these species causes serious yield loss if they are not appropriately controlled in the fields (
Hill et al., 1985
). The two species,
E. crus-galli
and
E. oryzicola
are problematic weeds in Korean agricultural fields as well. The Korean agronomists are confused about identifying this difficult genus using morphological characters because E. crus-galli not only exhibits very diverse morphological and ecological types, but also includes several varieties such as var.
echinata
(Willd.) A. Chev. and var.
praticola
Ohwi. However, these two varieties,
echinata
and
praticola
, have been treated as synonyms of var.
crus-galli
depending on taxonomists (The Plant List, 2014).
It has been known that eight morphological taxa of the genus,
E. crus-galli
var.
crus-galli
, var.
echinata
, var.
praticola, E. oryzicola, E. colona, E. glabrescens
Munro ex Hook.f.,
E. oryzoides
(Ard.) Fritsch, and
E. esculenta
(A.Braun) H.Scholz, are distributed throughout Korea (
Lee et al., 2013
).
Echinochloa crus-galli
grows in the upland fields and also in paddy fields whereas
E. oryzicola
is frequently found in paddy fields. These two species are mainly found in the whole country. Existence of
E. colona
and
E. glabrescens
in Korea, however, should be reexamined. The result of a previous study showed that two species
E. colona
and
E. glabrescens
collected from Korea was clustered with
E. crus-galli
. They were not clustered with
E. colona
collected from China and/or formed independent clade (in press).
The continuous use of the herbicides which have the same mode of action has led to select the herbicide resistant species in the fields. Since the first herbicide resistant species,
Monochoria vaginalis
var.
plantaginea
(Roxb.) Solms, was discovered in Korea in 1998, the herbicide resistant
E. crusgalli
and
E. oryzicola
have also been reported and spreaded in Korea (
Park et al., 2010
). Since then, herbicide resistant
Echinochloa
spp. is getting more attention to both agronomists and farmers for effective control of these noxious weeds. Therefore, it is important to identify the species exactly and to know the genetic diversity of the species for designing the rational strategies for weed management, especially for the herbicide resistant weeds.
Simple sequence repeats (SSR) have many advantages over DNA sequencing for inferring phylogeny due to their neutral but faster evolution that may lead to more informative characteristics. Although the usefulness of SSR for resolving phylogenetic relationships among the closely related species was suggested (Takezaki and Nei, 1996), a few SSR phylogenies have been reconstructed (
Petren et al., 1999
;
Ochieng et al., 2007
). A better knowledge of the distribution of genetic variation within the
Echinochloa
would help its taxonomic classification and designing a rational strategy for herbicide test. The objective of this study was to identify and summarize the
Echinochloa
species distributed in Korea by comparing the genetic variation and relationship among Korean
Echinochloa
species using SSR fingerprinting.
Materials and Methods
- Sampling, DNA extraction and quantification
total of 107
Echinochloa
accessions was subjected to the fingerprinting study. Most of the Korean
Echinochloa
species were collected by authors from August to September of 2011 throughout the South Korea and the foreign accessions were the same accessions used for previous study (in press). The plant materials were identified based on the descriptions of
Chen and Peterson (2006)
,
Michael (2007)
and
Park et al. (2011)
. Plant materials are presented in
Table 1
.
The 107 accessions used for genetic diversity analyses.
*E. sp C or G: considered as E. colona or E. glabrescens when they were collected.
A total genomic DNA was isolated from the green leaves of plant materials using a Genomic DNA Isolation Kit (NucleoGen, Germany) according to the instructions provided. The isolated DNA concentration and relative purity were checked using a Nanodrop ND-1000 (Dupont Agricultural Genomic Laboratory) and adjusted to 25 ng/μL for PCR amplification.
- SSR fingerprinting and Data analysis
Five SSR markers developed by
Danquah et al. (2002a)
were subjected to study the genetic variation of the Korean
Echinochloa
species. Forward primers of the five SSR loci were labeled with blue (FAM), green (NED), or yellow (HEX) fluorescent tags (AB-PEC, Foster City, CA). PCR reaction mixtures were prepared according to
Lee et al. (2005)
. The amplification was performed using a T100 Thermal Cycler (BIO-RAD, USA) that was set to run at 94℃ for 3min for initial denaturation followed by 30 cycles of 94℃ for 30 s, 55℃ for 1min, 72℃ for 90 s, and final extension at 72℃ for 5 min. The amplification product of 1 μL was combined with 10 μL of Hi-Di formamide and 0.5 μL of an internal size standard, Genescan-500 ROX (6-carbon-rhodamine) molecular size standards (35-500 bp). The samples were denatured at 94℃ for 3 min and analyzed with an ABI 3130xl Genetic Analyzer (Applied Biosystems / Hitachi, Inc., Foster City, CA, USA). Automated sizing and labeled alleles were determined and visualized relative to an internal Rox-labeled size standard using Genemapper 4.0 for automated data output.
Basic statistics at each SSR locus, including the number of alleles (NA), major allele frequency, gene diversity (GD), heterozygosity, and polymorphism information content (PIC), were calculated using the genetic analysis package PowerMarker ver. 3.25 (
Liu and Muse, 2005
). Genetic distances between each pair of the accessions were measured by calculating the shared allele frequencies using PowerMarker ver. 3.25 (
Liu and Muse, 2005
). The unweighted pair-group method with arithmetic averaging (UPGMA) was used to construct a phylogram from a distance matrix using MEGA ver. 5.2 (
Tamura et al., 2011
). The model-based clustering analysis was performed using all the 107 individuals and five independent runs of STUCTURE for each
K
value (the number of subpopulations) from two to 10 without prior population information. Runs were carried out by setting for 100,000 iterations with 100,000 burn-in and assuming an admixture model with correlated allele frequencies (
Falush et al., 2003
). For each
K
value, the runs showing the highest posterior probability of the data were considered. The true value of
K
was detected by an ad hoc quantity based on the second order rate of change of the likelihood function with respect to
K
(
Evanno et al., 2005
). An individual having more than 90% of its genome fraction value was assigned to a group.
Results and Discussion
- SSR polymorphism
The 107 samples examined for SSR polymorphism represent seven species of Korean and foreign origins. The
E
. sp. in
Table 1
were the taxa which were considered as either
E. colona
or
E. glabrescens
when they were collected. One of SSR markers, EC4, was not used to assess the genetic diversity of
Echinochloa
species because it did not show polymorphism across all the samples employed in this study. A total of 19 alleles, ranging from three (EC2) to six (EC3 and EC5), was detected from 107
Echinochloa
accessions using four SSR markers with an average of 4.75 alleles per loci. The major allele frequency per locus varied from 0.5047 (EC5) to 0.8551 (EC2). The highest gene diversity (0.647) and Polymorphism Information Content (PIC; 0.594) were detected for the EC5, while the lowest gene diversity (0.250) and PIC (0.223) were detected using EC2 (
Table 2
). It was not congruent with the result of
Danquah et al.(2002b)
which revealed the least in EC1 (0.549) and the greatest in EC2 (0.715). Although this study included seven different species from several different origins, the average gene diversity, 0.469, seemed relatively lower than others comparing to those (0.556 and 0.66, respectively) of
Danquah et al. (2002b)
and
Xu et al. (2004)
. The higher genetic diversity of
Danquah et al. (2002b)
might be due to inclusion of samples collected from broader geographic areas, such as the rice fields of Bangladesh, India, Côte D’Ivoire and Philippines. The low genetic diversity of the genus
Echinochloa
may mean that the genus will have not become a problem for agriculture in a certain changing environment. It could be, however, a big problem if they persist genes that can make the domesticated plants resistant to the intruder. Therefore, it is necessary to monitor the genetic diversity of this genus periodically in Korea.
Major allele frequency, number of genotypes, sample size, number of alleles, gene diversity, heterozygosity, and polymorphism information content of 107Echinochloaaccessions.
1Major allele frequency2No. of genotype 3No. of alleles 4Gene diversity 5Polymorphism information content
- Distance based phylogeny and population genetic structure
A genetic distance-based analysis was performed by calculating the shared allele frequencies among the 107 accessions and an unrooted and rooted UPGMA phylograms (
Fig. 1A
) were computed using MEGA5 software (
Tamura et al., 2011
). In the UPGMA phylogram, all
Echinochloa
accessions clustered into two main clades (I and II). The first group (I) included nine accessions which are
E. colona
, all coming from China and Taiwan, and its progenitor,
E. frumentacea
. Independent lineage of
E. colona-frumentacea
within the genus
Echinochloa
was also supported by many studies, including sequencing data derived from cpDNA and nrDNA in a previous study (in press) (
Yamaguchi et al., 2005
;
Aoki and Yamaguchi, 2008
). The second clade is composed of the remaining accessions, constituted by two clades (II-1 and II-2) again. II-1 is constituted by
E. crus-pavonis
and
E. oryzicola
. The rest of the accessions, a crus-galli complex, are arranged to the II-2 (
E. crus-galli
var.
crus-galli
, var.
praticola
, var.
echinata, E. oryzoides
, and
E. esculenta
). The close relationship between
E. crus-pavonis
and a crus-galli complex was supported by the maximum parsimony tree derived from cpDNA sequences, but not by that from nrDNA sequences (
Aoki and Yamaguchi, 2008
). We assume that the clustering in this clade II might be in relation to the number of chromosomes, 2n = 4× = 36 or 2n = 6× = 54; clade I is composed of the species having the number of chromosomes, 2n = 4× = 36, whereas clade II is consisted of the species with the hexaploids, 2n = 6× = 54. One interesting point is that
E. oryzicola
(E55) collected from Choonchon, Gangwondo, showed the close relationship with a crus-galli complex, but not with the other
E. oryzicola
accessions. Therefore, it is necessary to check the number of chromosomes of
E. oryzicola
(E55) collected from Choonchon, Gangwondo. If E55 has the number of chromosomes, 2n = 54, it will be a first report about the hexaploid
E. oryzicola
.
Yamaguchi et al. (2005)
called
E. crus-galli, E. esculenta
, and
E. oryzoide
s as a crus-galli complex and considered to belong to the same species with more supporting studies such as isozyme analysis (
Nakayama et al., 1999
) and polymerase chain reactionrestiction fragment length polymorphism (PCR-RFLP) technique (
Yasuda et al., 2002
).
The model based population genetic structure of 107
Echinochloa
individuals was conducted using STRUCTURE v2.3.3.(
Pritchard et al., 2000
). The true value of
K
was four and the real structure showed a clear peak of 107
Echinochloa
individuals was set at
K
= 4 (
Fig. 1B
). What means that the 107
Echinochloa
individuals should be divided into four subpopulations. In total, nine (C) and 14 (A) individuals were clearly assigned to two of the four subpopulations, whereas the rest of 84 individuals were divided into two different subpopulations (B and D). Subpopulations C and A were corresponded mainly to the clades I and II-1 of UPGMA dendrogram and subpopulations B and D were corresponded to the clade II-2.
A. UPGMA dendrogram generated from 107 Echinochloa individuals using four SSR markers. Clade I was consisted of E. colona and its progenitor, E. frumentacea. Clade II was divided into two subclades, II-1 and II-2. Subclade II-1 was composed of mainly E. oryzicola and E. crus-pavonis whereas the rest the individuals, E. crus-galli, E. oryzoides, and E. esculenta, were arranged to the subclade II-2. B. Model-based populations of 107 Echinochloa individuals. Each accession is divided into a number of hypothetical sub-populations based on the proportional membership coefficients totaling 1 at K = 4. The two of subpopulations, C and A, were corresponded to the I and II-1 and the other two, B and D, were arranged to the clade II-2 of the UPGMA dendrogram.
The three clusterings, including subclade II-1 and II-2, of UPGMA phylogram might be more reasonable than four subgroups of the population structure because of no distinct morphological characters for each subpopulation found between the two subpopulations, B and D. Therefore, we have concluded from the results of the previous study (in press) and this study as follows: First, we could not find
E. colona
and
E. glabrescens
in Korea for our studies. Vouchers (KHB 1045294 and 1086605, KHB 1264632 and 1264634) loaned from Korea National Arboretum were turned out
E. crus-galli
, but neither
E. glabrescens
nor
E. colona
, respectively. Samples collected as those two species were
E. crus-galli
as well. Therefore, their natural distribution in Korea should be more surveyed. Second, Biological varieties
praticola
and
echinata
of
E. crusgalli
should be treated as
E. crus-galli
because many integrated forms can be easily found in the nature. In addition, some of the characters traditionally used for distinguishing taxa, such as awn length of var.
echinata
, are affected by the amount of moisture available in the habitats (
Michael, 2007
). Third, the Korean
Echinochloa
should be summarized with four species, i.e.,
E. oryzicola, E. crus-galli, E. esculenta
, and
E. oryzoides
. The most abundant species is
E. crus-galli
. We also consider treating the latter three species, regarded as a crus-galli complex, as one species with care. Further detailed studies, however, should be warranted for better understanding of the
Echinochloa
in Korea.
Acknowledgements
This research was supported by a project from the Cooperative Research Program for Agricultural Science & Technology Development (Project No. PJ008548) of the RDA.
Aoki D.
,
Yamaguchi H.
2008
Genetic relationship between Echinochloa crus-galli and Echinochloa oryzicola accessions inferred from internal transcribed spacer and chloroplast DNA sequences
Weed Biol. Manage.
8
233 -
242
DOI : 10.1111/j.1445-6664.2008.00303.x
Chen S.
,
Peterson P.M.
2006
Flora of China
Science Press and Missouri Botanical Garden Press
Beijing and St. Louis, Missouri, USA
Poaceae (Gramineae)
515 -
518
Clayton W.D.
,
Renvoize S.A.
1999
Genera Graminum, Grasses of the World
Royal Botanic Gardens
Kew, London, UK
280 -
281
Danquah E.Y.
,
Hanley S.J.
,
Brookes R.C.
,
Aldam C.
,
Karp A.
2002
Isolation and characterization of microsatellites in Echinochloa (L.) Beauv. spp
Mol. Ecol. Notes
2
(1)
54 -
56
DOI : 10.1046/j.1471-8286.2002.00144.x
Danquah E.Y.
,
Johnson D.E.
,
Riches C.
,
Arnold G.M.
,
Karp A.
2002
Genetic diversity in Echinochloa spp. collected from different geographic origins and within rice elds in Côte d’Ivoire
Weed Res
42
(5)
394 -
405
DOI : 10.1046/j.1365-3180.2002.00300.x
Evanno G.
,
Regnaut S.
,
Goudet J.
2005
Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study
Mol. Biol.
14
(8)
2611 -
2620
Falush D.
,
Stephens M.
,
Pritchard J.K.
2003
Inference of population structure using multilocus genotype data. Linked loci and correlated allele frequencies
Genetics
164
(4)
1567 -
1587
Hill J.E.
,
Strange M.L.
,
Bayer D.E.
,
Williams J.F.
1985
Integrated weed management in California
Proc. West. Soc. Weed Sci.
Reno. NV
100 -
104
Lee J.
,
Baek H.-J.
,
Yoon M.-S.
,
Kim C.-Y.
,
Cho E.-G.
2005
Taxonomic review and genetic diversity of cowpea species and related taxa
Korean J. Breed. Sci.
37
(4)
187 -
191
Lee J.
,
Kim C.-S.
,
Lee I.-Y.
2013
Identification of Echinochloa oryzicola (Vasinger) Vasinger and E. oryzoides (Ard.) Fritsch in Korea
Korean J. Pla. Tax.
43
56 -
62
DOI : 10.11110/kjpt.2013.43.1.56
Michael P.W.
,
Barkworth M.E.
,
Capels K.M.
,
Long S.
,
Piep M.B.
,
(Eds.)
2007
Flora of North America, North of Mexico
Oxford University Press
New York, USA
(Eds.)
Echinochloa P. Beauv.
390 -
403
Nakayama Y.
,
Umemoto S.
,
Yamaguchi H.
1999
Identification of polyploid groups in the genus Echinochloa by isozyme analysis
J. Weed Sci. Tech.
44
205 -
217
DOI : 10.3719/weed.44.205
Ochieng J.W.
,
Steane D.A.
,
Ladiges P.Y.
,
Baverstock P.R.
,
Henry R.J.
,
et al.
2007
Microsatellites retain phylogenetic signals across genera in eucalypts (Myrtaceae)
Genet. Mol. Biol.
et al.
30
1125 -
1134
DOI : 10.1590/S1415-47572007000600016
Park S.H.
,
Lee Y.M.
,
Chung S.Y.
,
Chang G.S.
,
Kang W.C.
,
et al.
2011
Illustrated grasses of Korean
(Revised and enlarged edition)
Korea National Arboretum
Pocheon, Kyonggido, Korea
et al.
448 -
457
Park T.-S.
,
Ku B.-I.
,
Kang S.-K.
,
Choi M.-K.
,
Park H.-K.
,
et al.
2010
Response of the resistant biotype of Echinochloa oryzoides to ACCase and ALS inhibitors, and effect of alternative herbicides
Korean J. Weed Sci.
et al.
30
291 -
299
DOI : 10.5660/KJWS.2010.30.3.291
Petren K.
,
Grant B.R.
,
Grant P.R.
1999
A phylogeny of Darwin’s finches based on microsatellite DNA length variation
Proc. Roy. Soc., London, UK
266
321 -
329
DOI : 10.1098/rspb.1999.0641
Pritchard J.K.
,
Stephens M.
,
Donnelly P
2000
Inference of population structure using multilocus genotype data
Genetics
155
945 -
959
Tamura K.
,
Peterson D.
,
Peterson N.
,
Stecher G.
,
Nei M.
,
et al.
2011
MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods
Mol. Biol. Evol.
et al.
28
(10)
2731 -
2739
DOI : 10.1093/molbev/msr121
Xu Y.
,
Beachell H.
,
McCouch S.R
2004
A marker-based approach to broadening the genetic base of rice in the USA
Crop Sci
44
(6)
1947 -
1959
DOI : 10.2135/cropsci2004.1947
Yamaguchi H.
,
Utano A.
,
Yasuda K.
,
Yano A.
,
Soejima A.
2005
A molecular phylogeny of wild and cultivated Echinochloa in East Asia inferred from non-coding region sequences of trnT-L-F
Weed Biol. Manage.
5
(4)
210 -
218
DOI : 10.1111/j.1445-6664.2005.00185.x
Yasuda K.
,
Yano A.
,
Nakayama Y.
,
Yamaguchi H.
2002
Molecular identification of Echinochloa oryzicola Vasing. and E. crus-galli (L.) Beauv. using a polymerase chain reaction-restriction fragment length polymorphism technique
Weed Biol. Manage.
2
(1)
11 -
17
DOI : 10.1046/j.1445-6664.2002.00041.x