Progress report for GW21-218
Project Information
Cirsium arvense (L.) Scop. (Canada thistle) is an aggressive perennial weed which threatens sustainable crop production. Organic producers throughout Montana and across the Western region have expressed the need for the development of integrated Canada thistle management programs which include mechanical, cultural, and biological control tactics. The fungal pathogen, Puccinia punctiformis (thistle rust), is a highly selective Canada thistle parasite that has the potential to aid in the suppression of Canada thistle when integrated with other agronomic management tactics. This research and education project will (1) explore methods to integrate thistle rust into tillage practices and cropping rotations to reduce the spread and impact of Canada thistle in organic settings, and (2) convey the generated information to organic and conventional farmers to facilitate the development of integrated weed management practices.
In this study, we evaluate in certified organic farms and experimental plots: (1) the impact of conventional and reduced tillage practices on the establishment and spread of thistle rust in Canada thistle populations, and (2) Canada thistle growth and spread within competitive cropping rotations. Research results have been used to inform producers who wish to adopt tillage and rotational tactics which maximize disease establishment and minimize the impact of Canada thistle. The ecological principles developed in this study will be useful to both organic and conventional growers, and our findings have been shared with agricultural stakeholders through extension meetings, grower conferences, and virtual content. The ultimate contribution from this project is the sustainable suppression of Canada thistle infestations within agroecosystems.
The overall goal of this study is to assess the integration of the Canada thistle pathogen, Puccinia punctiformis (thistle rust) to manage Canada thistle in organic agroecosystems.
- Research Objective 1: Assess how thistle rust responds to conventional and reduced tillage practices.
- Research Objective 2: Compare the impacts of crop competition in rust infected Canada thistle and non-infected Canada thistle.
- Education Objective 1: Extend knowledge of thistle rust and its integration into organic weed management.
Cooperators
- - Technical Advisor - Producer
Research
Research Objective 1: Assess how thistle rust responds to conventional and reduced tillage practices.
An integration of thistle rust as a Canada thistle biocontrol in agricultural systems requires a clear understanding of the pathogen's response to management tactics. Thistle rust spores remain dormant within Canada thistle debris on the soil surface until they come in contact with a new host. Soil tillage has the potential to bury, transport, or destroy thistle rust spores which are responsible for transmission of the pathogen. This objective is designed to examine the potential for tillage to impact the spread and establishment of the thistle rust pathogen.
Cooperator Farm
During the summer of 2020, producer Ole Norgaard cooperated in the establishment of a field study on his organically certified farm near Shonkin, Montana. The farm has pre-existing Canada thistle infestations which contain varying degrees of naturally established thistle rust infection. The preliminary study explored tillage impacts on thistle rust establishment and spread. At Mr. Norgaard's farm, four discrete Canada thistle patches were identified, flagged and GPS mapped. Two of the patches contained varying degrees of symptomatic rust infected Canada thistle, while the other two patches showed no signs of infection. Tillage occurred in one patch containing thistle rust infection and one patch with no signs of infection, using the cooperator’s tillage method of choice. The remaining two patches were left untilled throughout the rest of the growing season. Prior to tillage, vegetation growth data was collected along 20-meter transects, placed through the center of each thistle patch. Ten 1-m2 sample frames were placed at even intervals along the transects. Vegetation density data was collected for symptomatically infected thistle stems and total thistle stems within each frame. Additionally, percent vegetation cover was estimated within frames, including % symptomatically infected Canada thistle, % asymptomatic Canada thistle, % crop, % other, and % bare ground.
Mr. Norgaard decided to plant alfalfa as a perennial forage crop in the fall of 2020. During the 2021 growing season, sampling methods were repeated in within the same four discrete thistle patches that were monitored in 2020. Mr. Norgaard did not till the soil or swathe the field because the alfalfa crop was in the early to mid-vegetative stage of growth. Low moisture conditions in the summer of 2021 led to poor establishment of Mr. Norgaard's alfalfa crop, and an expansion of the Canada thistle populations within the fields. As a result, Mr. Norgaard decided to remove the fields from organic certification in the fall of 2021. The fields which hosted this study were sprayed with 123.5 oz./hectare of glyphosate in the fall, followed by seeding of a forage mixture.
As a response to the new management, the methods for the 2022 sampling period were altered to best accommodate the project objectives. The new objectives at the Norgaard Farm focused on re-evaluation of the original four thistle patches within the crop fields, with a new evaluation of thistle patches in a grazed pasture that borders Mr. Norgaard's crop fields. The goal was to determine if there were any differences in thistle density and the density of thistle stems with rust symptoms between the two habitat types at the Norgaard Farm. The crop fields were revisited in June 2022 where all four experimental thistle patches were sampled under previously established protocols. Three additional thistle patches, located within grazed pastures along Mr. Norgaard's crop fields, were also sampled using the original protocols. Canada thistle stem density and density of thistle stems with rust symptoms data were analyzed in RStudio using a generalized linear model with the gaussian distribution. Comparisons were made between thistle growth and rust symptoms from 2022 in the crop field habitat and the pasture habitat, and analysis was conducted in RStudio. The margin of each thistle patch was also mapped with an Emlid Reach RS2 GNSS receiver on June of 2022, and point vectors were taken at the location of all observations of symptomatically infected Canada thistle within each patch. The data was analyzed in QGIS on the WGS 84 Pseudo-Mercator coordinate reference system. Patch areas from each patch polygon was derived in QGIS, and RStudio was used to evaluate the difference in mean patch areas between the crop field habitat and the pasture habitat.
Montana State University (MSU) Fort Ellis Research Farm
Integration of thistle rust into tillage is being evaluated at the MSU Ft. Ellis Research Farm in Bozeman, Montana. The research fields at the Ft. Ellis research farm are approximately 0.34 acres (0.14 hectares) each, with a mean annual precipitation of 56 cm and a mean annual air temperature of 6°C (USDA NRCS, 2019).
During the summers of 2020 and 2021, discrete thistle patches with pre-existing rust infections were surveyed and mapped within four organically certified fields. Each field was uniformly seeded with an organic green manure barley crop during the spring of 2020, at a rate of 197-lbs/hectare. In July, prior to Canada thistle anthesis, a census of symptomatically infected Canada thistle was taken from each thistle patch in all four replications. Additionally, asymptomatic Canada thistle stem density was counted within 1-m2 sample frames per thistle patch, using a simple random sampling method. After thistle assessment, all four fields were randomly split into a reduced tillage mowing termination and a conventional tillage disc termination, with nearly equal representation of thistle rust in each treatment. The mowing termination was performed with a 72” wide flail mower set to leave 4” stubble. Disc cultivation was set at an 8” soil depth, and three passes were made per conventional tillage treatment. In the fall of 2020, a green manure mixture of organic winter peas and winter triticale was seeded into all four fields at a rate of 81-lbs/acre winter pea and 75-lbs/acre winter triticale. Sampling protocols were repeated in July of the following summer, prior to Canada thistle anthesis. Reduced and conventional tillage events were conducted within assigned treatment areas immediately after sampling, followed by a secondary tillage treatment that took place in September. In May of 2022, the fields were prepared for seeding with respective reduced and conventional tillage treatments. Each field was sown with spring wheat seed, immediately after seedbed preparation. The spring wheat crop will be grown to full maturity, and the final sampling event will take place prior to harvest. All sampling protocols will be repeated, with the addition of a final biomass and crop yield assessment.
Canada thistle stem density and the density of thistle rust infected Canada thistle stems was evaluated in RStudio. Data was analyzed between the years 2020 and 2021, and between the disc and mowed treatments. A Levene test was used to check the assumptions of equal variance in the data. As a response to the Levene tests, Welch two sample t-tests were used to compare the effects of tillage and years on thistle stem density and density of rust infected thistle stems.
Research Objective 2: Compare the impacts of cropping management between rust-infected Canada thistle and non-infected Canada thistle.
The experiments under this objective consist of a greenhouse study and a field study that assess the effects of thistle rust on Canada thistle competitiveness and spread. The greenhouse study is designed to look directly at the growth responses of thistle rust inoculated Canada thistle when it is under stress from crop competition, and within a controlled environment. The field study is designed to evaluate the performance of the biocontrol agent in a setting that more closely represents a "real world" cropping system.
Montana State University Greenhouse
The effects of crop competition on non-inoculated (control) and thistle rust inoculated Canada thistle was evaluated at the MSU greenhouse in Bozeman, MT. Ten-inch diameter, eight-inch-deep greenhouse pots were sown with 1-gram (±0.2-grams) cuttings of adventitious Canada thistle rhizome, sourced from wild thistle populations in Gallatin County and Hill County Montana. Once the rhizome produced basal thistle rosettes, the pots were randomly assigned to thistle rust inoculated or control treatments. Treatments were segregated into separate greenhouses to prevent unintentional rust infection in the non-inoculation treatment. Each pot in the inoculated treatment received five-mL thistle rust spores topically applied to Canada thistle rosettes. The spores were sourced from symptomatically infected thistle stems in Gallatin County, and prepared by grinding symptomatic thistle leaves into a dry powder. Inoculations occurred at an average temperature of 22°C, and were misted once a day with deionized water for 72 hours after inoculation. Temperature and water applications were held constant for the duration of the study and supplemental lighting was used during the winter months. Nutrients were added bi-weekly in the form of a water-soluble 20-20-20 NPK fertilizer. Individual greenhouse pots were randomly rotated bi-weekly to prevent microclimate impacts. A competition treatment was nested into the control and thistle rust inoculated pots to simulate the impacts of a multi-phase crop rotation on inoculated thistle. The levels of competition included Canada thistle grown by itself (monoculture), Canada thistle grown with a crop, and a crop grown by itself (see figure 1).

A 4-phase crop rotation was used to assess the performance of thistle rust inoculated Canada thistle in response to competition from a multi-phase crop rotation consisting of phase 1.) fallow with 1-gram (±0.2-grams) of Canada thistle rhizome per pot, phase 2.) spring wheat seeded at 100 kg/hectare, phase 3.) forage peas seeded at 89 kg/hectare and, phase 4.) safflower seeded at 33 kg/hectare. Canada thistle was grown within its original pot throughout the four phases. Harvests occurred at the crop maturity stage within each phase (approximately every 3 months), where all above-ground thistle and crop vegetation were cut at soil-level and dried for 72 hours at 40.5°C, and weighed. After each harvest, the next crop phase was seeded at recommended seeding rates, and thistle was allowed to regrow. Once Canada thistle rosettes re-established, the inoculated treatment received a repeated application of thistle rust spores following the original protocols. Each trial was terminated after harvest of phase 4. At the termination of each trial the thistle rhizome was removed from pots, cleaned of soil and debris, dried, and weighed.
A relative competition index was used to evaluate the performance of above-ground Canada thistle growth among treatments. The thistle monoculture level of competition acted as an index for maximum potential growth, where Canada thistle had no interspecific competition to impact its growth. In contrast, the polyculture level represented a system of competition, where Canada thistle development was directly influenced by the growth of crop species throughout the 4-phase rotation. The relative competition index is a measure of plant performance which contrasts monoculture with polyculture dry biomass where Pmonoculture is the average Canada thistle performance (biomass) in monoculture and Ppolyculture is the average thistle performance (biomass) in polyculture.
RCI = ((Pmonoculture - Ppolyculture) ÷ (Pmonoculture)) x 100%
Relative competition from above-ground biomass of Canada thistle within three competitive phases of the crop rotation in trials one and two was calculated in RStudio. A linear mixed effects model was fit with relative competition as the response to inoculation treatments (inoculated & control), rotational phase (wheat, pea, safflower), and trial (one & two) as covariates.
Rhizome performance was compared for both treatments using a relative growth rate (RGR) that calculated the difference of rhizome weight between the beginning and end of each trial. The pots containing Canada thistle were all initially sown with 1-gram (±0.2-grams) pieces of rhizome. Relative growth of Canada thistle rhizome within the control, inoculated, monoculture, and polyculture treatments were evaluated after the four phases of growth. In the following RGR equation W1 represents initial rhizome weight at planting, W2 represents dry rhizome weight after 4 phases of growth, and t is the total days of growth for each trial.
RRGtreatment = (lnW2 -lnW1) ÷(t)
The relative growth rate of Canada thistle rhizome (grams/day), grown in monoculture and polyculture within trials one and two, was evaluated in RStudio. Linear mixed effects model was fit with relative growth rate as the response to inoculation treatments (inoculated & control), levels of competition (monoculture and polyculture), and trial (one & two) as covariates.
Montana State University Post Agronomy Farm:
The impacts of the thistle rust biocontrol on Canada thistle growth is being evaluated at the MSU Post Agronomy Farm in Bozeman, Montana. The research fields at the site are each 0.02 hectares with a mean annual precipitation of 48 cm and a mean annual temperature of 7°C (USDA NRCS, 2019). This experiment combines biological, mechanical, and cultural controls for Canada thistle infestations within agroecosystems. During the spring of 2020, 1-gram cuttings of thistle rhizome were planted in two separate fallow fields. Soil samples were taken from each field for nutrient comparisons, and Delmhorst model GB-1 gypsum blocks were installed 15 centimeters below the soil surface for soil moisture monitoring. Each field contains 10 thistle rust inoculated and 10 non-inoculated Canada thistle. To reduce potential for unintentional thistle rust inoculation of the non-inoculated control treatment, the controls were intentionally planted in the predominant upwind direction, 12 meters away from the inoculated treatments. During the summer of 2020, both fields were maintained in fallow so that the planted rhizome could establish. Thistle growth data was collected in the fall within 1-m2 sample frames centered over the original rhizome planting locations. Data collection included thistle stem density, density of thistle stems showing rust symptoms, % cover asymptomatic thistle, % cover symptomatic thistle, % cover crop, % cover bare ground, and % cover other. Fall inoculations were performed during an optimal climate range, where temperatures were forecasted between 8 and 25°C. A maximum of five thistle rosettes were inoculated within 1-m2 sample frames centered over the original rhizome planting locations. Each inoculation consisted of a 5-mL application of thistle rust spores, applied directly to thistle rosettes. The spores were sourced from symptomatically infected thistle stems in Gallatin County, and prepared by grinding symptomatic thistle leaves into a dry powder. Thistle was allowed to senesce naturally in each field, which were left no-till to avoid disturbance of patch development.
In the spring of 2021, 80-lbs/acre of forage barley was drill seeded into each no-till field with an in-furrow application of 13-36-0 fertilizer at a rate of 100-lbs/acre. Thistle was allowed to freely grow through the summer within the barley crop. Data collection occurred at the mealy ripening stage of the barley kernels. An Emlid Reach RS2 GNSS receiver was used to map the spread of thistle patch boundaries, and collect total patch area data for each treatment. Sampling of three additional transects was conducted, west and east from the original planting locations, to account for the spread of thistle beyond the original planting locations. Each transect was 20-m long, with ten 1-m2 samples taken at every 0.5-m interval. Data collection followed the same protocols that were established in 2020. Each field was swathed with a gas-powered trimmer, hand raked, and cleared of above-ground biomass following data collection. The fields were left no-till, and thistle rhizomes were allowed to re-establish through the fall. Once fall climate conditions reached an optimal range, a maximum of five thistle rosettes were inoculated within all 1-m2 sample frames along the four transects in the inoculated treatment. Inoculations were only applied within sample frames containing live thistle rosettes, following protocols established in 2020.
Analysis was conducted in RStudio, comparing thistle stem density between the inoculated and non-inoculated treatments over 2020 and 2021. A linear mixed effects model was fit with thistle stem density as a response to treatment levels (inoculated & non-inoculated), levels of growing season (2020 & 2021), and levels of block (A & B) as covariates. The spatial spread of the planted Canada thistle at the Post Agronomy Farm was derived in Arc GIS as total patch areas (m2), and analyzed in RStudio using a generalized linear model.
In the spring of 2022, spring wheat was drill seeded in the two no-till fields at a rate of 85-lbs/acre. The sampling protocols that were established in 2021 will be repeated at the harvest-ready stage of wheat in 2022. Above ground thistle biomass, crop biomass, and grain yield will be collected from all 1-m2 samples as a conclusion to the study.

Research Objective 1: Assess how thistle rust responds to conventional and reduced tillage practices.
Cooperator Farm
The mean thistle density per 1-m2 was 38.9 thistle stems in the crop field, and 34.8 thistle stems in the pasture. There is no evidence for a difference in mean thistle stem density per 1-m2 between habitat types at the Norgaard Farm in 2022 (generalized linear model, gaussian distribution, p-value= 0.770). The mean density of Canada thistle stems with rust symptoms per 1-m2 was 0.6 thistle stems with rust symptoms in the crop field and 1.5 thistle stems with rust symptoms in the pasture. There was no evidence for a significant difference in the mean density of Canada thistle stems with rust symptoms per 1-m2 between habitat types (generalized linear model, gaussian distribution, p-value= 0.134). There is strong evidence for a difference in mean patch area (m2) between the habitats at the Norgaard Farm in 2022 (generalized linear model, gaussian distribution, p-value<0.0001). The average patch area in the herbicide treated crop field was 3989.9- m2, and 331.0- m2 in the grazed pasture.
Habitat Type |
Management Type |
Average Thistle Patch Area (m2) |
Average Thistle Stem Density (1-m2) |
Average Thistle Stems with Rust Symptoms (1-m2) |
Average Thistle Stems with Rust Symptoms (Patch) |
Crop Field |
Herbicide |
3989.8 |
38.9 |
0.6 |
61.8 |
Pasture |
Grazing |
331.0 |
34.8 |
1.5 |
16.6 |
Mean Canada thistle stem densities (generalized linear model, p-value= 0.406) and stem densities with rust symptoms (generalized linear model, p-value= 0.12) were not impacted by the habitat type or management type at the Norgaard Farm. Mr. Norgaard's Fall 2021 application of glyphosate did not appear to have any significant impact on the overall thistle performance indicating that these thistle patches were resilient to a singular application of herbicide, likely attributed to well established below-ground energy stores. The differences in patch area between habitats may be attributed to past tillage events. Tillage within the crop field can change the soil profile by reducing compaction, leading to less obstruction to thistle rhizome spread. The pasture area consisted of perennial grasses and forbs that are rotationally grazed by 4 heads of cattle. The combination of vegetative competition, soil stability and minimized grazing has appeared to prevent large spatial spread of thistle within the pasture habitat. Thistle rust symptoms were generally consolidated to dense infection areas within patches, but showed no overall difference in densities per 1-m2 between habitat or management types. The observed grouping of symptoms within patches could mean that the pathogen is not effectively spread by grazers or by a combined legacy of tillage and herbicide application.

Montana State University (MSU) Fort Ellis Research Farm
The mean thistle stem density per 1-m2 within the standard tillage (disc) treatment was 56 thistle stems, and 55 thistle stems in the reduced tillage (mow) treatment. To check the assumption of equal variance, the Levene test was performed. According to the Lavene test, the variance between the standard tillage (disc) and reduced tillage (mow) treatments is significantly different (p-value <0.0001). It was concluded that thistle stem density per 1-m2 is not different between the standard tillage (disc) and reduced tillage (mow) treatments for 2020 and 2021 (Welch two sample t-test, p-value= 0.7432). A 95% confidence interval estimates that the true difference in mean thistle stem density per 1-m2 is between -8.7 thistle stems per 1-m2 and 6.2 thistle stems per 1-m2.
The mean thistle stem density per 1-m2 in 2020 was 69 thistle stems, and 45 thistle stems in 2021. The Levene test was performed to test the assumption of equal variance between years, resulting in a difference in variance between 2020 and 2021 (p-value<0.0001). There is very strong evidence for a difference in thistle stem density per 1-m2 between 2020 and 2021 (Welch two sample t-test, p-value<0.0001). A 95% confidence interval estimates that the true difference in mean thistle stem density per 1-m2 between 2020 and 2021 is between 17.2 stems per 1-m2 and 31.8 stems per 1-m2 (figure 4).
Canada thistle stem density per 1-m2 was reduced by an average of 54% between 2020 and 2021 within the study. However, the Ft. Ellis Research Farm experienced a drought from 2019 to 2021. The five-month precipitation average (March-July) was 22-inches of rain in 2019, 13-inches of rain in 2020, and 17-inches of rain in 2021 (wunderground.com). The soil moisture profile at the farm was highly reduced by the time vegetation data was collected in 2021.

The mean density of Canada thistle stems with rust symptoms within the standard tillage (disc) treatment was 9.5 thistle stems, and 8.2 thistle stems in the reduced tillage (mow) treatment. To check the assumption of equal variance, the Levene test was performed. According to the Lavene test, the variance between the standard tillage (disc) and reduced tillage (mow) treatments is significantly different (p-value <0.0001). Therefore, a Welch two sample t-test was used to evaluate the difference in means between the standard tillage (disc) and reduced tillage (mow) treatments. There is no evidence that density of thistle stems with rust symptoms per patch is different between the standard tillage (disc) and reduced tillage (mow) treatments for 2020 and 2021 (Welch two sample t-test, p-value= 0.319). A 95% confidence interval estimates that the true difference in mean density of thistle stems with rust symptoms per patch is between -4.0 thistle stems and 1.3 thistle stems (figure 5).
The difference in mean density of Canada thistle stems with rust symptoms per discrete thistle patch between 2020 and 2021 was also evaluated in RStudio. The mean density of thistle stems with rust symptoms per patch in 2020 was 6.1 symptomatic thistle stems, and 10.9 symptomatic thistle stems in 2021. The Levene test was performed to test the assumption of equal variance between years, resulting in a significant difference in variance between 2020 and 2021 (p-value<0.0001). A Welch two sample t-test was used to evaluate the difference in means Canada thistle stem density per patch between 2020 and 2021. There is very strong evidence for a difference in density of Canada thistle stems with rust symptoms per discrete thistle patch between 2020 and 2021 (Welch two sample t-test, p-value<0.0001). A 95% confidence interval estimates that the true difference in mean density of Canada thistle stems with rust symptoms per discrete thistle patch between 2020 and 2021 is between 6.1 symptomatic stems and 10.9 symptomatic stems.

Research Objective 2: Compare the impacts of cropping management between rust infected Canada thistle and non-infected Canada thistle.
Montana State University Greenhouse
The mean percentage of Canada thistle biomass loss within the control treatment of trails one and two was 49% loss the wheat phase, 65% loss in the pea phase, and 12% loss in the safflower phase. The mean percentage of Canada thistle biomass loss within the inoculated treatment was 57% loss in the wheat phase 74% loss in the pea phase, and 36% loss in the safflower phase. There was strong evidence for a difference in mean percent thistle biomass loss in competition between the wheat, pea and safflower phases of trials one and two (Type III analysis of variance with Satterthwaite's method, p-value=0.0087). There is no evidence for a difference for a difference in mean percent thistle biomass loss in competition between the inoculated and control treatments in trials one and two (Type III analysis of variance with Satterthwaite's method, p-value=0.1634).
The difference in competition observed between the three competitive crop phases reflects the different competitive abilities of the crops (figure 6). Wheat, pea, and safflower have different biological traits that determine vigor, rooting depth, canopy development, and resource consumption. The relative competition results in this greenhouse study validate biological differences between each crop. However, these results may suggest that competitive effects on Canada thistle were delayed by one phase of crop growth. It is generally expected that wheat is a strong competitor against weeds, peas are weak in a competitive system, and safflower is a moderate competitor. Our results show that average relative competition was 53% thistle biomass loss in wheat, 70% thistle biomass loss in peas, and 24% biomass loss in safflower. This may be due to a suppression of rhizome development of thistle when grown in polyculture with a strong competitive crop, that isn't expressed in above ground biomass until the following phase of thistle growth. In contrast, a weak competitor would facilitate an increase in thistle's rhizome development, that is expressed in above ground biomass during the following phase of thistle growth.

There is strong evidence for different thistle rhizome relative growth rates between the thistle grown in monoculture and thistle grown in polyculture (Type III analysis of variance with Satterthwaite's method, p-value=0.0062). There is also moderate evidence for different rhizome relative growth rates between thistle that was rust inoculated and thistle that was the non-inoculated control (Type III analysis of variance with Satterthwaite's method, p-value=0.0679) (figure 7).
The relative growth rate data from greenhouse trials one and two show that the integration of crop competition with the thistle rust biocontrol can reduce thistle rhizome biomass. Thistle rhizome biomass was lowered for rhizome that were grown in polyculture with three phases of competitive annual crops, as expected. The addition of thistle rust inoculations to thistle grown in polyculture further decreased the relative growth rate of thistle rhizome. The integration of these two management tools appears to be an effective way to reduce Canada thistle rhizome in a greenhouse setting. These results will be finalized by the end of the August of 2022, upon the completion of the final experimental greenhouse trial.

Montana State University Post Agronomy Farm:
The mean thistle stem density in 2020 was 8.2 thistle stems in the inoculated treatment and 7.4 thistle in the non-inoculated treatment. In 2021, the mean thistle stem density per 1-m2 was 6.9 thistle stems in the inoculated treatment and 7.6 thistle stems in the non-inoculated treatment. 10% of the twenty inoculated thistle produced thistle stems with rust symptoms in 2020, and 15% produced thistle stems with rust symptoms in 2021 (figure 8).
There was no significant difference in thistle stem density per 1-m2 between the inoculated treatment and non-inoculated treatment over the two years of thistle growth (Type III analysis of variance with Satterthwaite's method, p-value= 0.8651). There was also no difference in thistle stem density per 1-m2 between the two years, 2020 and 2021 (Type III analysis of variance with Satterthwaite's method, p-value= 0.4129).

After one year of thistle growth, the 40 original rhizome plantings had merged into discrete patchworks of thistle, with a maximum spread of approximately 6-meters away from the original planting location. The mean patch areas in 2021 were 82.2-m2 for the rust inoculated treatments and 92.0-m2 for the non-inoculated control treatments. There is no evidence for a difference in mean thistle patch area between the inoculated and non-inoculated treatments at the Post Agronomy Farm (generalized linear model, gaussian distribution, p-value= 0.142).

Research Outcomes
Education and Outreach
Participation Summary:
Newsletters
The Montana State University Extension webpage on Integrated Pest Management (IPM) is a public access resource that provides research-based education for citizens of the region. In 2021, the Cropland Weeds Laboratory published a newsletter on the MSU Extension IPM webpage that was focused on the management of Canada thistle with the thistle rust biocontrol. This newsletter discussed basic biology of Canada thistle and the thistle rust pathogen, how to identify thistle rust symptoms, how to collect and process spores, and how to inoculate Canada thistle with the biocontrol agent.
Presentations
The Montana Organic Association (MOA) is a collaboration of producers, businesses, and agencies whose goal is to promote and advocate for organic production in Montana. MOA hosts an annual conference at the end of each year where the public is given access to research presentations. Content from this study was presented in the 2021 MOA annual conference, in a virtual format. The presentation summarized the methods and preliminary results from our study in integration of thistle rust into organic management for Canada thistle. Approximately 40 participants attended the live virtual presentation, which was later uploaded to the MOA YouTube channel.
MOA 2021 Annual Conference Presentation
As part of the Montana State University Land Resources and Environmental Sciences graduate student requirements, this research was presented in a department seminar. Approximately 20 students, faculty and researchers attended the presentation. The fifteen-minute presentation summarized the biology of Canada thistle and the thistle rust biocontrol agent, our research methods, and preliminary results.
During the summer of 2022, the content from this study was also presented at the Montana State University Central Agriculture Research Center (CARC) organic field day, and at the Post Agronomy Farm field day. Approximately 20 producers, professionals, researchers and students attended the CARC organic field day. Approximately 50 producers, professionals, researchers and students attended the Post Agronomy field day.
Social Media
An Instagram account (@montanaweedecology) was created in the spring of 2021, to act an educational platform for a public audience. A total of 11 posts and 8 stories have been made public over the last year. All posts and stories have are focused on the work being conducted by the cropland weed ecology laboratory, with highlights on Canada thistle, thistle rust, and experimental management. The Instagram page has gained 583 public followers since its creation.