Canola (Brassica napus L.) is an important oilseed crop grown worldwide. In recent years, canola cultivation has expanded rapidly due to its high grain prices and demand for healthy vegetable oil, stockfeed and biodiesel markets . Weed infestation, however, is a major constraint limiting canola production . In Australia, total yield loss due to weeds in canola and pulses has been estimated at 122,048 tonnes, resulting in a revenue loss of $54 million . Annual ryegrass (Lolium rigidum Gaudin) has been the most widespread weed in winter crops, occurring in 86% of canola crops in south-eastern Australia .
The primary method of weed control is the application of herbicides. However, the prolonged and widespread chemical use has been increased especially after the introduction of herbicide-tolerant cultivars to triazine, imidazolinone and glyphosate. This has led to an increase in the evolution of herbicide resistance . Canola seems to be particularly vulnerable for developing to herbicide resistance as there are limited options of commercial herbicide available to control broad-leaf weeds. In recent years, the canola industry is increasingly becoming dependent on herbicide tolerant varieties including genetically modified herbicide tolerant varieties, which are meant to provide control options for the major weeds of that crop, such as annual ryegrass and wild radish (Raphanus raphanistrum). Evaluation of the herbicides against different weed species showed that 8 of the top 15 are likely to be utilised in canola production, including imazamox and imazethapyr for the Clearfield® HR canola, glyphosate for the Roundup Ready® canola, and atrazine and simazine for the triazine tolerant lines are resistant to herbicides . In recent years, some countries are imposing restriction on the usage of certain herbicides such as glyphosate for weed control; this practice necessitates the development of alternative and sustainable options for weed management. In addition to agronomic interventions that can influence weed management including seeding rate, row spacing, row orientation and fertilizer  ,crop interference is worth investigating as a tool for weed management  - . Considerable genotypic variation for weed competition exists in crop plants including canola, although some species are considered more competitive than others    . For example, studies have shown canola to be less competitive on weeds than wheat and barley   . Vigorous hybrid canola varieties were found to be more competitive than open-pollinated varieties, largely attributed to greater hybrid vigour traits such as plant height and early-season crop biomass accumulation     . The competitive ability of a crop to weeds can be measured either on the basis of the ability of crop to maintain growth and seed yield in the presence of weed, or on the basis of the ability of crop to suppress growth and seed production of weed species . Weed competitiveness in canola has been evaluated mainly under laboratory and field conditions and to a limited extent under glasshouse conditions   . Under field conditions, it is difficult to achieve precise and uniform plant densities across a trial site, and this may influence the differential responses obtained  . Secondly, field conditions can compromise the outcomes through environmental variance   . Lastly, screening large numbers of genotypes for weed-crop competition under field conditions is labour and space intensive .
The objective of this study is to evaluate the suppression ability of different canola genotypes against ryegrass. Obtaining reliable estimates of weed competitive ability and understanding the canola ideotype for interference traits are important for designing sustainable weed control strategies with low herbicide use for improving canola productivity and profitability.
2. Materials and Methods
2.1. Canola Genotype and Weed Populations
Previously, Asaduzzaman et al.  utilised 70 genotypes of canola to investigate genotypic variation for allelopathy and weed competitiveness. In this study, a set of 26 diverse Brassica genotypes (Table 1), including a subset of canola genotypes utilised by Asaduzzaman et al.  was characterised for their competitive ability against annual ryegrass cv. Wimmera under glasshouse and field conditions. This rye grass cultivar is well-suited to dry and low fertile soils and extensively used for productive, nutritious pasture crop. Seeds of canola genotypes were accessed from the National Brassica Germplasm Improvement Program (Wagga Wagga, Australia).
2.2. Evaluation of Canola-Weed Suppressive Ability under Glasshouse Conditions
The glasshouse experiment was conducted at the Wagga Wagga Agricultural Institute, NSW, Australia. The experiment was arranged in split-plot design with three replicates; main plots were 26 canola genotypes and subplots were the weed and weed-free treatments. Each experimental unit had four pots. The target canola density was five plants per pot in each treatment and ten annual ryegrass plants per pot in the weedy treatment. The 200 millimetre in diameter pots were filled with commercial garden mix containing: compost, manure, pine fines, sand and loam (3:3:2:1:1: by volume, 2 kg/pot). Fertiliser, 50 mL∙pot−1 (22.23 mL∙L−1 of water) of Thrive™ (NPK = 25:5:8.8), was added once per week. When canola plants started flowering, magnesium sulphate (2 g∙L−1 of water) was added once per week. Imidacloprid (15 g∙L−1) was used against aphids as needed, and Prothioconazole and Tebuconazole (210 g∙L−1) were used to control downy mildew disease as needed. The experiment was conducted under a 25/15˚C day/night temperature regime with a 16-h photoperiod. Pots were watered daily to field capacity, avoiding any confounding effect due to moisture stress. At 50 d after sowing, shoots of canola were cut at the soil level to estimate stem and leaf biomass. Numbers of leaves on the main stem were counted and then leaves were separated from the stem and placed in separate paper bags. Both leaves and stems were dried at 70˚C for 48 h, and weighed. The shoots of 10 annual ryegrass plants were harvested from each pot, dried at 70˚C for 48 h and weighed. Plant height was measured from the soil level to the top of the
Table 1. Phenology and country of origin of canola accessions.
canopy of the main (primary) stem.
2.3. Evaluation of Canola-Weed Suppressive Ability in the Field
The field experiments were conducted at the Wagga Wagga Agricultural Institute, Australia (35˚30'07"S; 147˚21'06"0E) in a duplex Red Kandosol of pH 5.3. The field had a history of naturally high annual ryegrass population. Herbicides including glyphosate (450 g∙L−1) and oxyfluorfen (450 g∙L−1) were used to control annual ryegrass before sowing. Canola was sown on June 1, 2017 (due to unavailability of sufficient soil moisture during the main canola growing season i.e. late April to mid-May) with sowing rate of 1400 seeds plot−1 and with 120 kg∙ha−1 of Croplift™ fertiliser treated with fungicide (Flutriafol, 500 g∙L−1). Urea (46% N) was applied twice at 75 kg∙ha−1, at sowing time and in the canola vegetative phase. Foliar applications of Prothioconazole® and Tebuconazole® (210 g∙L−1) against the fungal disease blackleg caused by Leptosphaeria maculans and Pirimicarb® (500 g∙kg−1) against aphids were applied as needed.
The experimental design was a complete randomised block with 26 genotypes (Sturt-TT was not included) and with three replicates in a spatially optimised layout of 5 ranges × 15 rows. Plot size was 1.6 m wide by 8 m long with 6 rows at a row spacing of 22 cm. On September 7, 2017 (100 days of sowing), canola and annual ryegrass densities were recorded in three random quadrats of 1 m by 1 m per plot. Canola biomass was determined from ten plants in the middle row of the plot whereas weed biomass was determined from three quadrats of 1 m by 1 m per plot. Both canola and weed biomass were cut at early canola flowering on September 20, 2017 and dried at 70˚C for 48 h. Plant height of canola was recorded twice in each plot by measuring the height of 10 randomly-selected plants per plot at early flowering and at crop maturity stages. The time of 50% canola flowering of each genotype was recorded (mid-September to mid-October). Plots were harvested at maturity with a small-plot header (Kingaroy, Australia) and grain yield was expressed in g plot−1.
3. Statistical Analysis
Glasshouse trial: Data were analysed using R software  and the ASReml package . Graphics were produced in the lattice package . The typical ASReml mixed model employed had fixed effects for genotype, weed (weedy v weed-free), and their interaction. Random effects were included (where significantly different from zero) for main plots and subplots.
Field trial: A Shapiro-Wilk test was used to examine the data for normality and scatter plot of residual versus fitted value was used to assess the homoscedasticity. On confirming data normality, the ASReml package in R was used to perform a linear mixed effect model. The observed data for crop shoot biomass (g plot−1) and crop height (cm plant−1) were normally distributed, whereas the data for yield (g plot−1) was transformed to obtain normal distribution. The Box-Cox transformation using the parameter lambda equal to 0.25 was used to transform yield data. The aim of the analysis was to examine both the genotype effects and weed biomass while accounting for the block stratification and spatial layout of the plots (including ryegrass density across plots). The mean with 95% confidence interval for yield plot−1, average crop height plot−1 at flowering time and crop shoot biomass plot−1 were presented graphically.
Correlation analysis: The Performance Analytics package in R was used for understanding their inter-relationships, including to describe their correlations in graphical format. The ggplot2 package in R was used to investigate the annual ryegrass response to 26 canola genotypes under the glasshouse and field conditions.
4.1. Genetic Variation for Weedcompetitive Ability under Glasshouse Conditions
The analysis revealed highly significant differences between genotypes and between weed treatments for all traits evaluated, whereas the interaction between the genotype and weed treatments was only significant in stem biomass and plant height (Table 2). Several genotypes revealed a strong ability to interfere with ryegrass growth under glasshouse conditions. Significant genotypic effects were found on shoot biomass of ryegrass that ranged from 1.05 to 2.28 g plant−1 (Figure 1(a)). Two winter varieties (cv. Akela and cv. Gross-Luesewitze) and a semi-winter variety (cv. Ningyou 7) were the most ryegrass suppressive genotypes, resulting in a lower annual ryegrass biomass of 1.05 to 1.26 g plant−1, followed by three Australian spring varieties (Av-Garnet, Av-Opal and Tarcoola-141). However, the other three Australian spring varieties (Sturt-TT, Ag-Spectrum and Lantern) were least suppressive, having higher levels of ryegrass biomass (2.15 to 2.28 g plant−1). These results suggested that growth habit of canola genotypes does not have any consistent relationship with ryegrass suppression.
Weed treatment influenced shoot biomass of canola genotypes as compared with the weed-free treatment (Figure 1(c)). However, the crop height was less affected by the presence of ryegrass (Figure 1(b)). This explains why crop biomass had higher negative correlation with weed biomass relative to the low correlation between crop height and ryegrass biomass. To identify competitive traits suitable for genetic selection, we sought correlation relationships between ryegrass and canola phenotypes. Leaf and shoot biomass showed a significant negative correlation (r = −0.50 to −0.76) with biomass of annual ryegrass (Figure 2). In contrast crop heights and leaf number had positive relationships with ryegrass biomass (Figure 2).
4.2. Validation of Genetic Variation for Weed Competitive Ability under Field Conditions
Our glasshouse experiment revealed that vigorous canola genotypes having
Table 2. P-values for fixed effects (genotype, weed treatment and genotype x weed interaction for each six traits measured under the glasshouse condition.
*n/a, Not applicable; only weedy treatment are included in the analysis of weed biomass.
Figure 1. Weed competitiveness of the 26 canola genotypes against ryegrass. Predicted means from the fitted model described in Table 1 are presented: (a) weed biomass (g plant−1), (b) crop height (cm plant−1), (c) crop biomass (g plant−1) and (d) leaf biomass (g plant−1), as predicted.
Figure 2. Correlation and correlation coefficient (r) of five traits (dry shoot biomass, dry leaf biomass, dry stem biomass, plant height and leaf number) for 26 canola genotypes and of biomass for annual ryegrass under the glasshouse condition.
higher leaf and shoot biomass control ryegrass better as compared to less vigorous genotypes. In order to validate the competitiveness of vigorous varieties against ryegrass, we conducted a field experiment using 26 canola genotypes. Results reconfirmed that plant vigour (shoot biomass) is a genetic trait and significantly affects ryegrass (shoot biomass), but the interaction between genotypes and weed biomass was non-significant (Table 3). This indicates that genotype and/or weed treatment is not influenced with growing environment (ryegrass influenced the crop biomass at a constant rate across all genotypes). As observed under glasshouse experiment, there were significant genotypic differences in plant height; however it did not significantly affect ryegrass (shoot biomass). The interaction between genotypes and weed biomass was also significant
Table 3. P-values for fixed effects (genotype) and covariate (weed biomass g plot−1) and interaction between genotype and weed biomass (g plot−1), crop biomass (g plot−1), crop height (cm plant−1) at flowering (growth stage 65, height 1) and physiological maturity (growth stage 80, height 2), measured under the field conditions.
*n/a, Not applicable, covariate included in the analysis.
(Table 3). The main effect of genotypes on grain yield was highly significant, whereas the effects of weed biomass and flowering time were not significant. The significant genetic effects of crop biomass and plant height may have accounted the variation in yield (Table 3).
4.3. Correlation between Weed Biomass, and Canola Development and Productivity Traits under Field Conditions
Crop biomass and grain yield showed significantly negative correlations with ryegrass biomass. As anticipated, crop biomass was strongly correlated (r = 0.55) with grain yield (Figure 3). Four canola genotypes: Ag-Outback, Av-Opal, Skipton and X6-06-3275-3 showed the higher dry biomass, indicating that these genotypes have strong ability for ryegrass competitiveness. These genotypes produced higher grain yield as compared to cvs. Ag-Spectrum, Beluga, CB-Argyle, Lantern and Tapidor that had the lower biomass and grain yield (Figure 4(a) and Figure 4(c)). Early flowering resulted in the increase of yield and showed a negative correlation (Figure 3). Crop height measured at both flowering and physiological stages was positively correlated with yield, with correlation coefficient 0.74 and 0.67, respectively. Five canola genotypes, PAK85388-502, Rivette, Av-Garnet, Tarcoola-141 and X6-06-3275-3, were taller and those genotypes had the higher yields (Figure 4(b) and Figure 4(c)). However, ryegrass biomass was poorly correlated with crop height at the flowering and physiological maturity (Figure 3). Winter varieties were the shortest compared with semi-winter and spring varieties which could be attributed to their vernalisation requirement.
4.4. Correlation between Weed Suppression Ability of Canola Genotypes under the Glasshouse and Field Conditions
A positive relationship was observed between ryegrass and crop biomass data
Figure 3. Correlation, correlation coefficient and normal distribution of five traits for 25 canola genotypes and of biomass for annual ryegrass grown under the field condition.
collected under glasshouse and field conditions (Figure 5). Our results showed that five canola genotypes: Akela, PAK85388-502, Av-Opal, Av-Garnet and Tarcoola-141 had lower weed biomass under both conditions, whereas Ag-Spectrum had the highest weed biomass. There was higher biomass of ryegrassin Lantern under the glasshouse condition but it was moderate under the field condition and this may have attributed to the variable density of ryegrass and require further verification.
Weed biomass, crop height, crop biomass and yield are the common proxy traits used to screen canola genotypes for the weed-crop competition under field conditions. However, field-based screening has some limitations which influence the reliability and accuracy of phenotypic estimations. Screening different germplasm for weed-crop competition under glasshouse conditions can provide
Figure 4. The 25 canola genotypes means with 95% confidence interval for: (a) crop biomass (g plot−1) as predicted from the fitted model: canola genotypes * weed biomass (g plot−1), (b) crop height at flowering stage (cm plant−1) as predicted from the fitted model: canola genotype * weed biomass (g plot−1), (c) grain yield (g plot−1) as predicted from the fitted model: canola genotypes + weed biomass (g plot−1) + crop biomass (g plot−1) + crop height at physiological stage (cm plant−1) + crop height at flowering stage (cm plant−1) + time to 50% flowering.
some advantages over field conditions. In our study we compared the competitiveability of 26 Brassica genotypes to suppress ryegrass under glasshouse and field conditions. For most traits, genotypes and weed treatments had a high significant effect. The interaction was only significant on crop height in both conditions. Five genotypes: Akela, Tarcoola-141, Av-Opal, Av-Garnet and PAK85388-502 consistently showed the ability to suppress weed ryegrass in both glasshouse and field conditions. Majority of these varieties had also the highest grain yields in the field conditions except in Akela.
A good understanding of morphology and phenology of both weed and crop could lead to an improvement in weed-crop competition via the use of genetic and genomic tools. In cereals, extensive leaf display, leaf area index, long flag leaves and good ground cover have been associated with superior competitive ability   . Morphological traits associated with the interception of radiation by leaves including plant height, leaf size, number and leaf area index, have been implicated in competiveness for light . Our study showed a
Figure 5. Correlation between two set of data of annual ryegrass biomass grown with 25 canola genotypes under the glasshouse and field conditions.
significant correlation between crop biomass and ryegrass suppression in both glasshouse and field conditions. These results indicate that plant vigour, manifested in higher canola biomass is associated with suppression of annual ryegrass. This is in agreement with previous studies where diverse canola accessions were compared under field conditions and showed that early-season crop biomass accumulation (early vigour) was associated with weed suppression   . In the glasshouse experiment, we found that leaf biomass was significantly correlated with weed suppression (r = −0.50), whereas crop stem biomass did not show any correlation (r = −0.05). Number of leaves also was not correlated with weed suppression. These results indicate that larger leaves are likely to provide shade and thus interfere with light inception to ryegrass. Further studies to understand the association between the area and angle of leaves and the suppression of weed growth are required for the germplasm used in our study. This would identify effective phenotypic traits for weed competition in a breeding program. Under field conditions, increase in crop biomass was associated with reductions in ryegrass biomass and with increased crop yields. This may reveal that canola biomass and particularly leaf biomass related to weed suppression and consequently improvement in yield.
The glasshouse and field experiments showed a low correlation between crop height and weed suppression but no correlation between crop stem biomass and weed suppression. Such lack of relationship has also been found in screening 111 rice genotypes against Echinochloa crus-galli in the field . However, this result is not in agreement with previous studies in canola. The glasshouse experiment (Figure 1(b)) indicated that crop height was less affected by the presence of weed compared with those in weed free conditions, whereas crop biomass was highly affected by weed (Figure 1(c)). Similarly, in the field experiment, crop height was not significantly influenced by the weed, whereas the crop biomass was. Most winter varieties found in our study are good weed suppressors in the glasshouse experiment (Figure 1(a)) and only Akela was superior in weed suppression in the field experiment (Figure 4). A spring genotype, Tarcoola-141 showed better ability in weed suppression under the glasshouse and field conditions with higher yield in the field. Other two spring genotypes, Av-Opal and Av-Garnet and a semi-winter PAK85388-502 demonstrated good weed suppression under both conditions as it was found in other studies  . Three spring genotypes including Ag-Spectrum, Lantern and Sturt-TT were found less to be weed suppressive under glasshouse and field conditions and had low grain yield in the presence of weeds.
Our study suggests that: 1) vigorous canola varieties can provide competitiveness to ryegrass under glasshouse and field conditions and 2) glasshouse conditions can be used to evaluate weed suppression ability of a large number of canola accessions while maintaining uniform density of annual ryegrass. We conclude that vigorous genotypes with highest weed suppression can be exploited for weed control in canola. Further research is required to develop structured populations between highest and lowest weed suppression genotypes can be developed to understand genetics underlying this valuable trait. Molecular markers can be developed for marker assisted selection leading to an acceleration of improved varietal development pipeline.
Authors are thankful to NSW Department of Primary Industries for supporting this research.
 AOF (2017) Australian Oilseeds Federation. Crop Report.
 Lemerle, D., Blackshaw, R., Smith, A., Potter, T. and Marcroft, S. (2001) Comparative Survey of Weeds Surviving in Triazine-Tolerance and Conventional Canola Crops in South-Eastern Australia. Plant Protection Quarterly, 16, 37-40.
 Llewellyn, R.R.D., Ouzman, J., Walker, S., Mayfield, A. and Clarke, M. (2016) Impact of Weeds on Australian Grain Production: The Cost of Weeds to Australian Grain Growers and the Adoption of Weed Management and Tillage Practices. Report to GRDC and CSIRO, Australia.
 Broster, J., Pratley, J., Ip, R., Ang, L. and Seng, K. (2019) A Quarter of a Century of Monitoring Herbicide Resistance in Lolium rigidum in Australia. Crop and Pasture Science, 70, 283-293.
 Heap, I. (2017) The International Survey of Herbicide Resistant Weeds.
 Asaduzzaman, M., An, M., Pratley, J.E., Luckett, D.J. and Lemerle, D. (2014) Canola (Brassica napus) Germplasm Shows Variable Allelopathic Effects against Annual Ryegrass (Lolium rigidum). Plant Soil, 380, 47-56.
 Asaduzzaman, M., Luckett, D., Cowley, R., An, M. and Pratley, J. (2014) Canola Cultivar Performance in Weed-Infested Field Plots Confirms Allelopathy Ranking from in Vitro Testing. Biocontrol Science and Technology, 24, 1394-1411.
 Bais, H.P., Weir, T.L., Perry, L.G., Gilroy, S. and Vivanco, J.M. (2006) The Role of Root Exudates in Rhizosphere Interactions with Plants and Other Organisms. Annual Review of Plant Biology, 57, 233-266.
 Lemerle, D., Verbeek, B., Cousens, R.D. and Coombes, N.E. (1996) The Potential for Selecting Wheat Varieties Strongly Competitive against Weeds. Weed Research, 36, 505-513.
 Asaduzzaman, M., Luckett, D.J., An, M., Pratley, J.E. and Lemerle, D. (2014) Management of Paterson’s Curse (Echium plantagineum) through Canola Interference. In: Nineteenth Australasian Weeds Conference, Hobart, 162-165.
 Wu, H., Pratley, J.E., Lemerle, D. and Haig, T. (2000) Laboratory Screening for Allelopathic Potential of Wheat (Triticum aestivum) Accessions against Annual Ryegrass (Lolium rigidum). Australian Journal of Agricultural Research, 51, 259-266.
 Daugovish, O., Thill, D. and Shafii, B. (2002) Competition between Wild Oat (Avena fatua) and Yellow Mustard (Sinapis alba) or Canola (Brassica napus). Weed Science, 50, 587-594.
 Lemerle, D., Luckett, D., Lockley, P., Koetz, E. and Wu, H. (2014) Competitive Ability of Australian Canola (Brassica napus) Genotypes for Weed Management. Crop and Pasture Science, 65, 1300-1310.
 Raman, H., Shamaya, N. and Pratley, J. (2018) Genetic Variation for Weed Competition and Allelopathy in Rapeseed (Brassica napus L.). In: Radhakrishnan, R., Ed., Biological Approaches for Controlling Weeds, IntechOpen, London, United Kingdom.
 Lutman, P.J.W., Dixon, F.L. and Risiott, R. (1994) The Response of Four Spring-Sown Combination Arable Crops to Weed Competition. Weed Research, 34, 137-146.
 Melander, B. (1993) Modelling the Effects of Elymus repens (L.) Gould Competition on Yield of Cereal, Peas and Oilseed Rape. Weed Research, 33, 99-108.
 Lemerle, D.V. and Coombes, N.E. (1995) Losses in Grain Yield of Winter Crops from Lolium rigidum (Gaud.) Depend on Crop Species, Cultivar and Season. Weed Research, 35, 503-509.
 Harker, N.C., Blackshaw, R. and O’Donovan, J. (2003) Seeding Rate, Herbicide Timing and Competitive Hybrids Contribute to Integrated Weed Management in Canola (Brassica napus). Canadian Journal of Plant Science, 83, 433-440.
 Harker, N., O’Donovan, J., Blackshaw, R., Johnson, E., Holm, F. and Clayton, G. (2011) Environmental Effects on the Relative Competitive Ability of Canola and Small-Grain Cereals in a Direct-Seeded System. Weed Science, 59, 404-415.
 Lemerle, D., Lockley, P., Koetz, E. and Wu, H. (2014) Competitive Ability of Australian Canola (Brassica napus) Genotypes for Weed Management. Crop and Pasture Science, 65, 1300-1310.
 Coleman, R., Gill, G. and Rebetzke, G.J. (2001) Identification of Quantitative Trait Loci for Traits Conferring Weed Competitiveness in Wheat (Triticum aestivum L.). Australian Journal of Agricultural and Research, 52, 1235-1246.
 Asaduzzaman, M., An, M., Pratley, J.E., Luckett, D.J. and Lemerle, D. (2014) Laboratory Bioassay for Canola (Brassica napus) Allelopathy. Journal of Crop Science and Biotechnology, 17, 267-272.
 Zand, E. and Beckie, J. (2002) Competitive Ability of Hybrid and Open-Pollinated Canola (Brassica napus) with Wild Oat (Avena fatua). Canadian Journal of Plant Science, 82, 473-480.
 Worthington, M. and Reberg-Horton, C. (2013) Breeding Cereal Crops for Enhanced Weed Suppression: Optimizing Allelopathy and Competitive Ability. Journal of Chemical Ecology, 39, 213-231.
 Johnson, D.E., Dingkuhn, M., Jone, M.P. and Mahamane, M.C. (1998) The Influence of Rice Plant Type on the Effect of Weed Competition on Oryza sativa and Oryza glaberrima. Weed Research, 38, 207-216.
 Team, R.D.C. (2014) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.
 Butler, D.G., Cullis, B.R., Gilmour, A.R. and Gogel, B.J. (2009) ASReml-R Reference Manual. Release 3.0. Technical Report, Queensland Department of Primary Industries, Australia.
 Donald, C.M. and Hamblin, J. (1976) Biological Yield and Harvest Index of Cereal as Agronomic and Plant Breeding Criteria. Advances in Agronomy, 28, 361-405.
 Olofsdotter, M., Navarez, D., Rebulanan, M. and Streibig, J. (1999) Weed-Suppressing Rice Cultivars: Does Allelopathy Play a Role? Weed Research, 39, 441-454.