Minnesota (U.S.A.) population or North American forage cultivar codes, river/wet meadow name and location or forage cultivar name and germplasm source; GPS coordinates for site of collection (wild populations) or germplasm bank identifier number for
Abstract
The spread of invasive plant species in natural habitats has become a worldwide problem with negative impacts. Phalaris arundinacea, an important forage and ornamental crop, is widespread worldwide. In recent years there has been a massive spread of P. arundinacea across North America and Canada. Production of Phalaris seed as a forage crop occurs in northern Minnesota; seeds are sold throughout the world, particularly in central Europe. We tested genetic similarities and differences between populations in the US (Minnesota) and the Czech Republic using ISSRs to determine potential gene flow for this forage crop. The cultivated forage and wild genotypes were dispersed into two groups that overlapped. At least four sets of wild US genotypes are dissimilar to European counterparts and potentially native to N. America. Future work to prove the ancestry of each accession will be necessary. Nonetheless, the sale of forage cultivars related to or derived from European types causes genetic mixing with N. American types. Part of this intercontinental gene flow is exacerbated by the production of Phalaris forage seed in Minnesota. The implications of these findings for management of invasive crops native to both continents are significant for forage producers, managers, and breeders.
Keywords
- reed canarygrass
- invasive species
- forage cultivars
- ornamental cultivars
- ISSRs
1. Introduction
In Europe, the standard forage cultivar is ‘Palaton’ (from the US), while other important ones include: ‘Luba son. Motycka’ (Poland); ‘Motterwizer’ (Denmark); ‘Peti’, ‘Szarvasi 50’, ‘Szarvasi 60’, ‘Keszthelyi 52’ (Hungary); ‘Lara’ (Norway); ‘Vantage’, ‘Venture’ (US); ‘Bellevue’, ‘Rival’ (Canada); ‘Chrastava’ (Czech Republic) [23]. Current breeding is focused on improving of its yield potential as a fodder crop as well as for wastewater treatment plants and, more recently, biomass production.
Despite unverified assertions that “reed canarygrass is native to the northern half of the United States…” and “native to the temperate portions of Europe, Asia, and North America” [32], invasion biologists and ecologists have consistently postulated that
Original and introduced
Previous work by our labs [18, 31] analyzed phenotypic and genotypic markers in genotypes obtained from wild populations growing along the six main rivers within the Czech Republic (Berounka, Dyje, Labe, Lužnice, Orlice, Vltava) and commercial cultivars (forage, ornamental types) grown in the Czech Republic to serve as a foundation for Central European reed canarygrass diversity. ISSRs or inter-simple sequence repeats, for the first time ever, showed distinct genetic differences between ornamental cultivars and wild
The objective of the present study was to extend the focus on assessment of genetic structure to wild
2. Materials and methods
2.1. Genotypes
A total of 16 wild
Population or forage cultivar codes | River/wet meadow name and location or forage cultivar name and germplasm source | GPS1 coordinates for site of collection (wild populations) or germplasm bank identifier number; [citations] |
---|---|---|
2.1.2 | St. Croix River, South of Bayport, MN; by the Bayport Marina | Lat.2 45°0′32.8710″ N Long.2 −92°46′40.4286″ W |
6.1.3 | St. Croix River, St. Croix State Park; along the river by boat launch and swimming areas | Lat. 45°57.012′ N Long. −92°34.044′ W |
8.I.A.1; 8.I.C.3; 8.I.G.3; 8.II.A.2; 8.II.F.3 | Wet Meadow, Chanhassen, MN; Horticulture Research Center′s “Rice Paddy” wetlands | Lat. 44°51′43.3296″ N Long. −93°35′59.4126″ W |
9.3.1 | Mississippi River, Reno, MN; along the dead arms, S from the dam of the “big” lake | Lat. 43°36.128′ N Long. 91°16.151′ W |
14.2.1 | Mississippi River, Red Wing, MN; along the river banks in a wooded area | Lat. 44°35′03.9444′′ N Long. 92°38′39.6918″ W |
21.5.1 | Mississippi River, between Little Falls and Rice, MN; in open areas between wooded banks | Lat. 45°49.597′ N Long. 94°21.262′ W |
34.3.1 | Mississippi River, near the headwaters; W of Bear Den Landing, Mississippi Headwaters State Forest | Lat. 47°26.012′ N Long. 95°07.748′ W |
38.1.B.3 | Minnesota River, Blakeley, MN; W of Belle Plaine. MN; in open wet meadows | Lat. 44°36′47.1708′′ N Long. 93°51′32.8320′′ W |
38.2.3 | Minnesota River, Blakeley, MN; W of Belle Plaine. MN; in open wet meadows | Lat. 44°36′43.7214′′ N Long. 93°51′35.2620′′ W |
46.1.1 | Minnesota River, SE of Montevideo, MN at the confluence of Highways 212/15 | Lat. N 44°54′09.8′′ N Long. 95°41′07.9′′ W |
50.1.1 | Des Moines River, S of Petersburg, MN at the border with the State of Iowa | Lat. 43°31′33.2′′ N Long. 94°55′07.4′′ W |
54.3.2 | Des Moines River, SW of Dovray, MN; adjacent to Highway 8 | Lat. 44°00′09.1′′ N Long. 95°35′00.3′′ W |
56.2.2 | Roseau River, in the Red Lake State Wildlife Mgt. Area, W of Mulligan Lake, adjacent to the Red Lake Indian Reservation; Co. Rd. 704, at headwaters (source) of the river | Lat. 48°32.774′ N Long. 95°19.204′ W |
58.1.3 | Roseau River, N of Roseau, MN; wet meadows near Hwy. 3 | Lat. 48°54.504′ N Long. 95°49.778′ W |
58.2.2 | Roseau River, N of Roseau, MN; wet meadows near Hwy. 3 | Lat. 48°54.546′ N Long. 95°49.711′ W |
58.3.1 | Roseau River, N of Roseau, MN; wet meadows near Hwy. 3 | Lat. 48°54.562′ N Long. 95°49.635′ W |
58.IV.A.1 | Roseau River, N of Roseau, MN; wet meadows near Hwy. 3; transect in cultivated field | Lat. 48°54.699′ N Long. 95°52.130′ W |
58.IV.H.3 | Roseau River, N of Roseau, MN; wet meadows near Hwy. 3; transect in cultivated field | Lat. 48°54.753′ N Long. 95°52.084′ W |
61.1.2 | Roseau River, Caribou, MN; Hwy. 4 near confluence with State Ditch; S of the Canadian Border | Lat. 48°59.006′ N Long. 96°26.951′ W |
63.4.3 | Red River, S of McCauleville, MN and SW of Kent, MN | Lat. 46°26′43.0″ N Long. 96°42′57.9″ W |
74.1.2 | Red River, Oslo, MN; S of Big Woods, County Ditch 38 | Lat. 48°18′40.3″ N Long. 97°07′24.4″ W |
VEN | ‘Venture’ (Minnesota); derived from crossing ‘Vantage’ × ‘Flare’); low alkaloid variety; does not contain any tryptamine-carboline alkaloids; USDA-GRIN; https://npgsweb.ars-grin.gov | PI3 531089 [14, 42] |
PAL | ‘Palaton’ (Minnesota); derived from ‘Flare’, ‘Vantage’ and ‘Rise’); low alkaloid variety; does not contain any tryptamine-carboline alkaloids; USDA-GRIN; https://npgsweb.ars-grin.gov | PI 531088 [14] |
AUB | ‘Auburn’ (Alabama); landrace, most likely derived from native N. American germplasm; USDA-GRIN; https://npgsweb.ars-grin.gov | PI 422031 [42] |
IOR | ‘Ioreed’ (Iowa); high levels of alkaloids; USDA-GRIN; https://npgsweb.ars-grin.gov | PI 422030 [42] |
365 | 367 (British Columbia, Canada); USDA-GRIN; https://npgsweb.ars-grin.gov | PI 387929 |
PHA |
|
PI 241065 |
PN-609 | Unknown origin; USDA-GRIN; https://npgsweb.ars-grin.gov | PI 371754 |
GRO | ‘Grove’ (Ontario, Canada); USDA-GRIN; https://npgsweb.ars-grin.gov | PI 357645 [42] |
MN-76 | MN-76 (Minnesota) 4-clone double cross hybrid; low alkaloid variety; does not contain any tryptamine-carboline alkaloids; USDA-GRIN; https://npgsweb.ars-grin.gov | PI 578797 [42] |
CANA | ‘Cana’ (California); USDA-GRIN; https://npgsweb.ars-grin.gov | PI 578795 |
VAN | ‘Vantage’ (Iowa); high alkaloid content; USDA-GRIN; https://npgsweb.ars-grin.gov | PI 578794 [14, 42] |
MCRC1 | NCRC-1 (Minnesota); USDA-GRIN; https://npgsweb.ars-grin.gov | PI 578793 |
SUP | ‘Superior’ (Oregon); most likely derived from native N. American germplasm; USDA-GRIN; https://npgsweb.ars-grin.gov | PI 578792 [14, 42] |
Table 1.
1GPS, global positioning system.
2Lat., latitude; Long., longitude.
3PI, plant introduction; USDA-GRIN, U.S. Dept. of Agriculture, Germplasm Resources Information Network.
We included data from our previous paper [18] and that of Kávová’s M.S. thesis [31] for comparative purposes, namely 110 European genotypes from Czech wild populations (1 site/river; 1–9 genotypes/collection site/river) collected in 2011 along the six main rivers of the Czech Republic (Berounka, Dyje, Labe, Lužnice, Orlice, and Vltava). Similar to the Minnesota wild populations, five of the Czech rivers empty into the North Sea basin while the Dyje River flow into the Black Sea basin [18]. Additional wild population samples were made at the OSEVA PRO, Ltd., Grassland Research Station (Rožnov-Zubří, CZ); commercial forage, ornamental cultivars either bred and/or grown in the Czech Republic were also included. These all were grown and previously analyzed in our previous study [18] and ISSR molecular data from these were used herein to compare with the results found with the Minnesota and N. American types. Genotypic codes for all Czech germplasm consisted of the following: BE-1, 2, 3 (Berounka); DY1, 2, 3 (Dyje); LA1 (Labe); LU1, 2, 3 (Lužnice); OR1, 3 (Orlice); VL1, 3 (Vltava); CHR (‘Chrastava’; forage cultivar); Z13, Z77, Z83, Z124, Z125 (OSEVA PRO, Ltd, Grassland Research Station, Rožnov-Zubří, CZ); ZP/COV1, 17 (Gardening Pelikán, Spálené Poříčí), AT/P6, 7 (‘Picta’), AT/T2, 6 (‘Tricolor’), F/L1, 4 (‘Luteopicta’), F/Pa3, 4 (
2.2. Genetic analyses
Genetic variability was assessed using ISSR markers. This molecular technique is often used in studies focused on genetic variation of plant populations and plant germplasm and we verified its suitability and stability in analyses of Phalaris genotypes. ISSR is also marker system with high detectable extent of genetic variation/diversity and also with the ability to detect the genetic diversity among individual accessions.
2.3. DNA extraction and ISSR analyses
DNA extraction from leaf samples and subsequent ISSR analyses of all Minnesota and N. American samples followed the protocols delineated by Kávová [31]. Four primers from the University of British Columbia were used to generate scorable ISSR markers: UBC 810—[GA]8T, UBC 825—[AC]8T, UBC 881—G3[TGGGG]2TG, and UBC 890—VHV[GT]7 [31]; these have been used in our subsequent studies for
2.4. Statistical analyses
Genetic distance matrices were created with Nei and Li’s [44] metrics. PCoA (principal coordinate analysis) and UPGMA (unweighted pair group method with arithmetic mean) cluster analyses were calculated with MVSP, version 3.1 (Multi-Variate Statistical Package; Kovach Computing Services U.K.) and DARwin, version 5.0.158 (Dissimilarity Analysis and Representation for windows; CIRAD, F) software. Genetic structure was calculated using STRUCTURE version 2.3.4, a Bayesian clustering algorithm (Admixture Model; correlated allele frequencies; K = 2, K = 4, K = 6, and K = 10 groupings; 100,000 burnin repetitions) [43, 45]. STRUCTURE groupings refer to relationship patterns. After plotting, the K = 2 grouping had the necessary decrease in slope and increase in variance, diagnostic of the true K value, with the greatest number of genotypes/grouping; all other groupings were eliminated [30]. Only results from the K = 2 grouping will be shown.
3. Results
The four ISSR primers generated 76 scorable bands (56.6% were polymorphic). The UPGMA cluster analysis showed three distinct grouping of genotypes, all of which separated at a genetic distance of 0.0 (
Figure 1
). The first grouping consisted of strictly forage cultivars from Iowa (PAL, VEN), Minnesota (MN-76), California (CANA) and Missouri (AUB) (
Table 1
), all of which differed significantly (

Figure 1.
UPGMA, based on ISSR markers, for wild Minnesota populations and N. American comparative forage cultivars of
Principal coordinate analysis (PCoA) of inter-simple sequence repeat (ISSR) markers in reed canarygrass

Figure 2.
Principal coordinate analysis (PCoA) of inter-simple sequence repeat (ISSR) markers in reed canarygrass
When the wild and cultivated US genotypes were comparatively analyzed for PCoA together with the Czech/European genotypes [2] this resulted into forming two primary clusters ( Figure 3 ). Cluster I (lower circle) included all samples from wild Czech (European) populations along rivers and the forage ‘Chrastava’ as established for European genotypes by Anderson et al. [18]; this cluster was enriched with all samples of US origin. All US genotypes were clustered into a small, oval sub cluster of Cluster I, on the border of the European wild genotypes and showing high similarity in ISSR marker pattern ( Figure 3 ). Cluster II (upper oval) is represented by European horticultural and forage cultivars and genotypes from The Nursery of Genetic Resources, OSEVA PRO, Ltd., Grassland Research Station (Rožnov-Zubří, Czech Republic) with both variegated and nonvariegated leaf types [18].

Figure 3.
Principal coordinate analysis (PCoA) of inter-simple sequence repeat (ISSR) markers in reed canarygrass
Assessing the genetic structure of analyzed Czech cultivated and wild genotypes showed classification of genotypes according to Q1/Q2 values (membership probabilities in the C [rows or genotypes] × K [columns or clusters] matrix for a single cluster analysis); K = 2 had the best stratification in STRUCTURE ( Figure 4 ). One group, ‘PN-609’, contains several forage cultivars and a few wild genotypes from Site 8. Whereas the larger group, ‘54.3.2’ contains the remaining genotypes from all rivers, wet meadows and any remaining forage cultivars.

Figure 4.
Genetic structure analysis of the US and N. American reed canarygrass collection using STRUCTURE software package (Admixture Model, allele frequencies correlated, K = 6, length of burnin period: 100,000). Key: the population code is located left from the corresponding color bars with two groups of accessions: black— ‘PN-609’; grey—‘54.3.2’.
UPGMA analyses of both the US and Czech populations, based on ISSRs, showed distinct groupings of reed canarygrass genotypes ( Figure 5 ). The first group was a small set of 6 genotypes, ZPCOV, collected at The Nursery of Genetic Resources, OSEVA PRO, Ltd., Grassland Research Station (Rožnov-Zubří, Czech Republic). The second grouping was a large series of sub clusters divided as follows. The most distant genotypes from the ZPCOV cluster were primarily horticultural cultivars from the Czech Republic along with one sole US genotype from the wet meadow in Chanhassen, MN (8.I.A.1; Figure 5 and Table 1 ). The next cluster was divided into two groups of: (a) 6 Czech genotypes (3 wild from the Vltava River and 3 GFP/GNP or ‘Picta’). Next were two sub clusters which bifurcated at a genetic distance of ~ 0.2 ( Figure 5 ). One formed a small grouping of 18 genotypes, namely Czech accessions and MN genotypes while the other was a large grouping of all remaining wild and cultivated US and CZ genotypes.

Figure 5.
UPGMA analysis, based on ISSR data, of combined Czech, Minnesota and N. American reed canarygrass samples analyzed (
4. Discussion
There were two overlapping groupings for the forage cultivars and wild reed canarygrass
Based on the UPGMA analysis of US cultivated and wild types of reed canarygrass (
Figure 1
), potentially the 4 wild populations from the Mississippi (34.3.1, 38.2.3), Minnesota (46.1.1), and Red (63.4.3) rivers are the least related to the N. American forage cultivars SUP (‘Superior’), VAN (‘Vantage’) and MCRC-1 and may be native American genotypes. These MN wild populations also differed from the Czech wild populations (
Figure 4
). Casler et al. [42] and Jakubowski et al. [38, 39] used 15 SSR molecular markers to distinguish among N. American native and exotic (European)
In the STRUCTURE analysis of the US reed canarygrass collected along Minnesota rivers and in wet meadows, along with the North American cultivars, the cultivars were distributed throughout both groups (
Figure 4A
and
B
). This was unexpected and surprising since, for instance, the Red and Roseau Rivers running through northern Minnesota do not flow to the Gulf of Mexico and the Atlantic Ocean via the Mississippi River, but instead flow to Manitoba, Canada into Lake Winnipeg and would have limited opportunities for gene exchange. Additionally, since reed canarygrass is native in Minnesota, there could have been divergent evolution within isolated rivers creating distinct populations but this was not found to be the case. This could be due to wind pollination, which may allow for gene flow (pollen) between rivers. Also likely could be the small sample sizes collected along all rivers and/or the choice of genetic markers that, even though they are polymorphic among the populations and cultivars, may not be able to discriminate among
In the PCoA and STRUCTURE analyses of both the N. American and European sample sets from river habitats and forage cultivars, no clear differentiation among groupings was evident (
Figure 5
). The pattern of genetic markers in the European (Czech) genotypes from alluvial habitats was inclusive of all US wild and cultivated forage genotypes. Both groups of genotypes (wild/cultivated) overlapped and, in contrast with our previous analysis of the European genotypes [18], it was not possible to distinguish between wild and cultivated genotypes with precision. One reason for this may be the low genetic variation and differentiation among genotypes and their high genetic similarity. What is surprising is the very low extent of genetic variability among US genotypes, which formed one “dense” group in this pooled analysis, with low levels of genetic dissimilarity. This fact may also explain poor differentiation between US wild genotypes and US cultivated forage, because of their low genetic dissimilarity. Another reason may be the small sample sizes tested herein. Future work will be devoted to conducting a more thorough sampling of wild and commercial
Since all grasses, including reed canarygrass, are anemophilous (wind-pollinated), it would be easy for genetic mixing to occur in adjacent plantings of cultivated and wild types. Likewise, as most forage cultivars bred and/or produced in Minnesota and N. America are closely related to or derived from European types [42], this also could be a reason why the wild and forage types overlapped in their genetic similarity (
Figures 1
and
2
). The numerous influxes of exotic, European types and cross-pollination effects (either occurring naturally or by hand pollination by plant breeders), combined with migration have mixed the gene pools [42]. For instance, while ‘Rival’ has both European and Scandinavian ancestors, ‘Ioreed’ is a hybrid mixture with the European nuclear haplotype but N. American cytoplasmic haplotype [42]. However, maintaining the integrity of N. American
5. Conclusion
In Minnesota populations of
References
- 1.
Baldini RM. Revision of genus Phalaris arundinacea . Webbia. 1995;49 :265‐329 - 2.
Lavergne S, Molofsky J. Reed canary grass ( Phalaris arundinacea ) as a biological model in the study of plant invasions. Critical Reviews in Plant Sciences. 2004;23 :415-429 - 3.
Lewandowski I, Scurlock JMO, Lindvall E, Christou M. The development and current status of perennial rhizomatous grasses as energy crops in the US and Europe. Biomass & Bioenergy. 2003; 25 :335-361 - 4.
Nilsson D, Hansson PA. Influence of various machinery combinations, fuel proportions and storage capacities on cost of co-handling of straw and reed canary grass to district heating plants. Biomass & Bioenergy. 2001; 20 :247-260 - 5.
Rice J, Pinkerton B. Reed canarygrass survival under cyclic inundation. Journal of Soil and Water Conservation. 1993; 48 :132-135 - 6.
Sheaffer CC, Rosen CJ, Gupta SC. Reed canarygrass forage yield and nutrient uptake on a year-round wastewater application site. Journal of Agronomy and Crop Science. 2008; 194 :465-469 - 7.
Chekol T, Vough LR, Chaney RL. 2002. Plant-soil-contaminant specificity affects phytoremediation of organic contaminants. International Journal of Phytoremediation. 2002; 4 :17-26 - 8.
Lasat MM, Norvell WA, Kochian LV. 1997. Potential for phytoextraction of Cs-137 from a contaminated soil. Plant Soil. 1997; 195 :99-106 - 9.
Olsen FJ, Chong SK. Reclamation of acid coal refuse. Landscape Urban Planning.1991; 20 : 309-313 - 10.
Samecka-Cymerman A, Kempers AJ. Concentrations of heavy metals and plant nutrients in water, sediments and aquatic macrophytes of anthropogenic lakes (former open cut brown coal mines) differing in stage of acidification. Science of the Total Environment. 2001; 281 :87-98 - 11.
Vymazal J. Constructed wetlands for wastewater treatment in the Czech Republic. Water Science & Technology. 2001; 44 :369-374 - 12.
Vymazal J. Vegetation development in subsurface flow constructed wetlands in the Czech Republic. Ecological Engineering. 2013; 61 :575-581 - 13.
Zhu T, Sikora FJ. Ammonium and nitrate removal in vegetated and unvegetated gravelbed microcosm wetlands. Water Science & Technology. 1995; 32 :219-228 - 14.
Sheaffer CC, Marten GC, Rabas DL, Martin NP, Miller DW. Reed Canarygrass. Minnesota Agricultural Experiment Station Bulletin No. 595-1990 (Item No. AD-SB-5533). St. Paul, MN; 1990 - 15.
Alway FJ. Early trials and use of reed canary grass as a forage plant. Journal. American Society of Agronomy. 1931; 23 :64-66 - 16.
Casler MD, Undersander DJ. Selection for establishment capacity in reed canary grass. Crop Science. 2006; 46 :1277-1285 - 17.
Cagaš B. Reed canary grass ‘Chrastava’. Czech Journal of Plant Breeding and Genetics. 2008; 44 (1):41-42 - 18.
Anderson NO, Kávová T, Bastlová D, Čurn V, Kubátová B, Edwards KR, Januš J, and Květ J. Phenotypic and genotypic variation in Czech forage, ornamental and wild populations of reed canarygrass 56 (5):2421-2435 - 19.
Juntila O, Landgraff L, Nilsen AJ. Germination of Phalaris seeds: Seed problems. Acta Horticulturae. 1978;82 :163-166 - 20.
Conchou O, Pautou G. 1987. Modes of colonization of an heterogenous alluvial area on the edge of the Garonne River by Phalaris arundinacea . Regulated Rivers. 1987;1 :37-48 - 21.
Anonymous. Ribbon grass. New England Farmer. 1834; 13 :41, 129 - 22.
Harris A. Ribbon grass. New England Farmer. 1835; 14 :125 - 23.
Strasil Z. Ekonomická analýza vybraných energetických rostlin urených prospalování. In: Sbor. Technika a technologie pro nepotravináské vyuití pody a jejíudrování v klidu., Brno. 2000; pp. 17‐22 - 24.
Dostál J. Nová kvtena CSSR 2. Prague: Academia Praha; 1989 - 25.
Weber E. Invasive Plants of the World. CAB International, Wallingford, UK: CABI Publishing; 2003 - 26.
Čížková H, Rychterová J, Hamadejová L, Suchý K, Filipová M, Květ J, Anderson NO. Biomass production in permanent wet grasslands dominated with Phalaris arundinacea: Case study of the Třeboň basin biosphere reserve, Czech Republic. In: Vymazal V, editor, The Role of Natural and Constructed Wetlands in Nutrient Cycling and Retention on the Landscape. Springer International, Cham (ZG), Switzerland. 2015. DOI 10.1007/978-3-319-08177-9_1 - 27.
Gyulai G, Mester Z, Kiss J, Szeman L, Idnurm A, Heszky L. Somaclonal breeding of reed canary grass ( Phalaris arundinacea L.). Grass and Forage Science. 2003;58 :210-214 - 28.
Lavergne S, Molofsky J. Increased genetic variation and evolutionary potential drive the success of an invasive grass. Proceedings of the National Academy of Sciences of the United States of America. 2007; 104 :3883-3888 - 29.
Morrison SL, Molofsky J. Effects of genotypes, soil moisture, and competition on the growth of an invasive grass, Phalaris arundinacea (reed canary grass). Canadian Journal of Botany. 1998;76 :1939-1946 - 30.
Nelson MF, Anderson NO, Casler MD, Jakubowski AR. Population genetic structure of N. American and European Phalaris arundinacea L. as inferred from inter-simple sequence repeat markers. Biological Invasions. 2014;15 :353-363. DOI: 10.1007/s10530-013-0525-9 - 31.
Kávová T. Genetická variabilita v populacích chrastice rákosovité ( Phalaris arundinacea L). MS Thesis. University of South Bohemia, Ceske Budejovice, CZ. 2013 - 32.
Schoth HA. Reed Canary Grass. Washington, DC: United States Department of Agriculture. Farmers’ Bull. No. 1602; 1938 - 33.
Galatowitsch SM, Anderson NO, Ascher P. Invasiveness in wetland plants in temperate North America. Wetlands. 1999; 19 :733-755 - 34.
Saltonstall K. Cryptic invasion by a non-native genotype of the common reed, Phragmites australis , into North America. Proceedings of the National Academy of Sciences of the United States of America. 2002;99 :2445-2449 - 35.
Dore WG, McNeill J. Grasses of Ontario. Ottawa, Canada: Agriculture Canada, Research Branch, Biosystematics Research Institute; 1980 - 36.
Merigliano MF, Lesica P. The native status of reed canarygrass ( Phalaris arundinacea L.) in the inland northwest, USA. Natural Areas Journal. 1998;18 :223-230 - 37.
Lavoie C, Dufresne C, Delisle FF. The spread of reed canary grass ( Phalaris arundinacea ) in Québec: A spatio-temporal perspective. Écoscience. 2005;12 (3):355-375 - 38.
Jakubowski AR, Casler MD, Jackson RD. Has selection for improved agronomic traits made reed canarygrass invasive? PLoS ONE. 2011; 6 (10):e25757. DOI: 10.1371/journal.pone.0025757 - 39.
Jakubowski AR, Casler MD, Jackson RD. Genetic evidence suggests a widespread distribution of native North American populations of reed canarygrass. Biological Invasions. 2013; 15 :261-268 - 40.
Jakubowski AR, Jackson RD, Casler MD. The history of reed canarygrass in North America: Persistence of natives among invading Eurasian populations. Crop Science. 2014; 54 (1):210-219 - 41.
Maurer DA, Zedler JB. Differential invasion of a wetland grass explained by tests of nutrients and light availability on establishment and clonal growth. Oecologia. 2002; 131 : 279-288 - 42.
Casler MD, Phillips MM, Krohn AL. DNA polymorphisms reveal geographic races of reed canarygrass. Crop Science. 2009; 49 :2139-2148 - 43.
Nelson MF, Anderson NO. How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms. Ecology and Evolution. 2013; 3 (10):3455-3470 - 44.
Nei M, Li W-H. 1979. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proceedings of the National Academy of Sciences of the United States of America. 1979; 76 (10):5269-5273 - 45.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000; 155 :945-959 - 46.
Densmore F. How Indians Use Wild Plants for Food, Medicine, and Crafts - 47.
Kindscher K, Noguera E. Cultural use of plants from the Baker Wetlands - 48.
Lerman N. Okanogan (Salish) ethnology - 49.
Turner NJ, Bouchard R, Kennedy DID. Ethnobotany of the Okanagan-Colville Indians of British Columbia and Washington. Victoria, BC: British Columbia Provincial Museum; 1980