# Increasing plant group productivity through latent genetic variation for cooperation

Samuel E Wuest, Nuno D Pires, Shan Luo, François Vasseur, Julie Messier, Ueli Grossniklaus, Pascal A Niklaus

► To cite this version:

Samuel E Wuest, Nuno D Pires, Shan Luo, François Vasseur, Julie Messier, et al. Increasing plant group productivity through latent genetic variation for cooperation. PLoS Biology, 2022, 20 (11), pp.e3001842. 10.1371/journal.pbio.3001842. hal-04309707

HAL Id: hal-04309707

https://hal.inrae.fr/hal-04309707

Submitted on 27 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

METHODS AND RESOURCES

# Increasing plant group productivity through latent genetic variation for cooperation

Samuel E. Wuest<sup>1,2,3*</sup>, Nuno D. Pires<sup>1</sup>, Shan Luo<sup>4</sup>, Francois Vasseur<sup>5</sup>, Julie Messier<sup>6</sup>, Ueli Grossniklaus<sup>1</sup>, Pascal A. Niklaus<sup>2</sup>

1 Department of Plant and Microbial Biology & Zurich-Basel Plant Science Center, University of Zurich, Zurich, Switzerland, 2 Department of Evolutionary Biology and Environmental Studies & Zurich-Basel Plant Science Center, University of Zurich, Zurich, Switzerland, 3 Group Breeding Research, Division Plant Breeding, Agroscope, Wädenswil, Switzerland, 4 Lancaster Environment Centre, Lancaster University, Lancaster, United Kingdom, 5 CEFE, Univ Montpellier, CNRS, EPHE, IRD, Montpellier, France, 6 Department of Biology, University of Waterloo, Waterloo, Canada

* samuel.wuest@agroscope.admin.ch

## Abstract

Historic yield advances in the major crops have, to a large extent, been achieved by selection for improved productivity of groups of plant individuals such as high-density stands. Research suggests that such improved group productivity depends on “cooperative” traits (e.g., erect leaves, short stems) that—while beneficial to the group—decrease individual fitness under competition. This poses a problem for some traditional breeding approaches, especially when selection occurs at the level of individuals, because “selfish” traits will be selected for and reduce yield in high-density monocultures. One approach, therefore, has been to select individuals based on ideotypes with traits expected to promote group productivity. However, this approach is limited to architectural and physiological traits whose effects on growth and competition are relatively easy to anticipate. Here, we developed a general and simple method for the discovery of alleles promoting cooperation in plant stands. Our method is based on the game-theoretical premise that alleles increasing cooperation benefit the monoculture group but are disadvantageous to the individual when facing noncooperative neighbors. Testing the approach using the model plant *Arabidopsis thaliana*, we found a major effect locus where the rarer allele was associated with increased cooperation and productivity in high-density stands. The allele likely affects a pleiotropic gene, since we find that it is also associated with reduced root competition but higher resistance against disease. Thus, even though cooperation is considered evolutionarily unstable except under special circumstances, conflicting selective forces acting on a pleiotropic gene might maintain latent genetic variation for cooperation in nature. Such variation, once identified in a crop, could rapidly be leveraged in modern breeding programs and provide efficient routes to increase yields.

## Introduction

Crop breeding is currently undergoing fundamental transformations. Speed breeding and genomic prediction can shorten generation times and increase effective population sizes,

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Wuest SE, Pires ND, Luo S, Vasseur F, Messier J, Grossniklaus U, et al. (2022) Increasing plant group productivity through latent genetic variation for cooperation. PLoS Biol 20(11): e3001842. https://doi.org/10.1371/journal.pbio.3001842

**Academic Editor:** Pamela C. Ronald, University of California, Davis, UNITED STATES

**Received:** March 30, 2022

**Accepted:** September 20, 2022

**Published:** November 29, 2022

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pbio.3001842

**Copyright:** © 2022 Wuest et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The datasets described and a basic analysis script are available through the Zenodo data repository (DOI:10.5281/zenodo.6983283).
**Funding:** This work was supported by the University of Zurich, Agroscope and an Ambizione Fellowship (PZ00P3_148223) of the Swiss National Science Foundation (to SEW), and an Advanced Grant of the European Research Council (to UG). FV acknowledges funding from the French Agency for Research (ANR grant ANR-17-CE02-0018-01, 'AraBreed'), and the Agreenskills fellowship programme (grant agreement n° 3215), which has received funding from the EU's Seventh Framework Programme under the agreement N° FP7-609398. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

**Abbreviations:** G-I, group versus individual; GWAS, genome-wide association study; NIL, near-isogenic line; SNP, single nucleotide polymorphism.

leveraging rates of phenotypic change to unprecedented levels [1,2]. At the same time, large-scale, high-throughput phenotyping platforms have become available and allow for the simultaneous quantification of multiple traits in ever larger greenhouse and field settings [3]. Yet, some studies suggest that current rates of yield increase are insufficient to meet growing demands from increasing human and animal populations under a changing climate [4,5].

Historically, the highest rates of yield increase were achieved in the middle of the 20th century, at the beginning of the “Green Revolution.” Combining breeding with improved management, yield potentials of major crops, such as wheat and tropical rice, approximately doubled in only a few plant generations [6–8]. In retrospect, these gains in yield potential have been unusually large. In contrast to most classical breeding that operates through selection on polygenic variation, they were largely realized by capitalizing on single genes, notably by the introgression of discrete but pleiotropic dwarfing alleles with major effects on plant form and function [9,10]. These alleles resulted in smaller and less bushy individuals, which diverted less resources to competition and allocated more into reproduction, i.e., into grain yield. The resulting more “cooperative” genotypes could be cultivated in high-density stands, which resulted in an increased crop yield per unit land area, with little productivity increase at the individual level. In other words, a trade-off between individual fitness and group-level yield was exploited to promote yield in well-fertilized high-density stands [8,11–17]. This trade-off, however, also causes the traits that benefit the group to be selected against in mixed or segregating populations containing larger, more competitive individuals, where cooperative individuals have lower relative fitness [14,18].

Thus, a key lesson from the “Green Revolution” is that selection for competitive individuals must be avoided to achieve the highest-possible yields. One way to achieve this in breeding is to anticipate “ideotypes” with a suite of phenotypic traits that reduce competition and thereby promote cooperation. Traits desirable for the breeder are, for example, a short stature, erect leaves, and a compact root system [12,13,19]. However, a limitation of ideotype breeding is that the specific trait variation relevant for high crop yield in high-density groups often remains elusive to the human observer, especially for belowground traits. In principle, this can be avoided by shifting the focus on selection for quantitative genetic variation for cooperative traits at the group level [20–22], but this is often impractical (especially early in the breeding cycle), requiring more land and resources for phenotyping. Further, when genetic variation for the trait of interest does not exist in the breeding population, it first must be discovered in a wider crop population or in wild relatives (e.g., secondary gene pool).

Cooperative strategies in natural plant communities may be selected for under certain circumstances [e.g., spatial clustering of closely related individuals, or in cases where individuals recognized close relatives (kin recognition) and reduce competitive responses] [23,24], but when competing individuals are unrelated, the evolutionary stable strategy generally is one that maximizes individual performance (and resource use) to the detriment of the group as a whole [25–28]. However, if genetic variation for cooperation exists in nature, for whatever reason, it could potentially be leveraged in breeding [29]. The difficulty, however, is to discover relevant alleles in a diverse set of plants, in particular if they do not correspond to classical ideotypes. Here, we present a methodological framework for the discovery of such alleles. We aimed for such a framework to be as general and unbiased as possible in order to detect yield or biomass gains (or any other performance measure) that emerge from any type of cooperation, including for unknown resources and through unidentified traits. Specifically, we designed competition experiments and corresponding analytical tools that allow to rank plant genotypes on a scale ranging from “competitive” to “cooperative.” We then applied these methods in a proof-of-concept experiment with a population of *Arabidopsis thaliana* genotypes and produced a genetic map of a group versus individual (G-I) performance trade-off.
Variation at one genomic region was strongly associated with variation along this trade-off, whereby the allele of minor frequency promoted cooperation. However, the identified genetic variants were not associated with individual, or monoculture-level, productivity. They only became evident when mapping the G-I performance trade-off itself. In separate experiments, we then showed that the associated genetic polymorphisms allowed predicting productivity responses of genotypes along plant—density gradients. We argue that the method we present here has a large application potential, not least because of recent advances in technology including genome-wide association studies (GWAS) and large-scale phenotyping that make these methods widely applicable.

## Results

### Quantifying a plant group-versus-individual (G-I) performance trade-off by growing genotypes in different social contexts

The core principle of our method is to measure the performance of a set of genotypes in different “social contexts” [30], allowing for the quantification and subsequent genetic mapping of variation in a G-I trade-off (Weiner [30] discusses this trade-off in detail). To achieve this, we grew individuals of a set of genotypes (1) in monoculture, i.e., we confronted them with their own social strategy; and (2) in competition with a set of “tester” genotypes that represent the range of social strategies present in the population. According to game theory [28], the most cooperative genotype of the set will perform best with similarly cooperative neighbors (i.e., as a group in monoculture) but will lose when facing selfish, highly competitive neighbors (i.e., as an individual in mixtures). Conversely, the most competitive genotype will perform worst when grown with similarly competitive neighbors but will fare well when facing a cooperative neighbor. Comparing the performance of a genotype in monoculture relative to the performance with the tester set therefore allows to quantify the genotypes’ G-I performance trade-off (Fig 1) and thereby rank all genotypes on a continuum between cooperative and competitive. In principle, this approach is applicable to any species, but for simplicity we tested it with an association panel of 98 natural *A. thaliana* genotypes—a subset of the RegMap population [31]. Aboveground dry matter production (shoot biomass) served as the measure of performance (Materials and methods), but any other target characteristic, such as agricultural yield, could also be used. Each of the 98 focal genotypes was grown in a pot that contained another congenotypic individual (genotype monoculture). We further grew individuals of all genotypes in full factorial combination with an individual of each of ten tester genotypes (genotype mixtures; Figs 1A and 1B and S1). This design was replicated in two blocks.

As expected, competitive interactions among individuals were strong, with large negative effects of average tester shoot biomass (average across all pots) on shoot biomass of the focal genotypes (ANOVA $F_{1,960} = 88.23; P < 0.001$). To evaluate the G-I performance trade-off of genotypes, we related the monoculture shoot biomass of the target genotype to the average individual biomass of the same genotype grown in competition with tester genotypes (average individual performance in mixed stand; Fig 1C). Not surprisingly, across genotypes, group and individual performance were highly positively associated, with more vigorous genotypes producing more biomass both in monoculture groups and as individuals subject to competition by testers. This relationship was also nonlinear (second degree polynomial; $F_{1,95} = 8.4, P = 0.005$ for quadratic term), possibly because space became progressively limited as plant size increased [30,32]. However, the exact nature of this deviation from linearity is unimportant since we treated the overall relationship as heuristic. We then used the distance of each genotype’s data from this empirical relationship to locate each genotype on an orthogonal axis that quantified the G-I trade-off (Materials and methods and Fig 1C). In other words, this
<figure><img src="image_4.png" /><figcaption>**Fig 1.** A general method for the genetic dissection of the group-vs-individual performance (G-I) trade-off. (A) Experimental design of the competition experiment. G1, G2, ... G98: focal genotypes sampled from the RegMap panel of natural *A. thaliana* genotypes. T1, T2, ... T10: tester genotypes representing a range in plant size to reflect a large portion of the genetic variation present within *A. thaliana*. (B) Experimental setup. (C) Relationship between a genotype's mean performance as an individual across all mixtures with tester genotypes (average individual shoot biomass) and its group performance in monoculture (average combined shoot biomass of the two individuals in a pot). Black line: second degree polynomial with 99% confidence intervals. The inset outlines three of many possible genetic effects a hypothetical allele could have on a genotype's strategy: variation along the mixture-vs-monoculture regression describes variation in vigor, and variation perpendicular to the regression describes the G-I performance trade-off. Red and blue dots show genotypes carrying different alleles at position 15'294'955 on chromosome 3 (see below), whereby red alleles confer a more cooperative strategy. Data available at https://zenodo.org/record/6983283, file competition. csv.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001842.g001
procedure transforms the separate values for group performance in genotype monoculture and mean individual performance in genotype mixtures into two orthogonal metrics: the position along the overall relationship reflects general genotypic vigor (e.g., increased productivity due to better adaptation to the specific growth conditions), and the position perpendicular to the general relationship reflects a G-I trade-off value that characterizes the communal properties of the focal genotype (inset Fig 1C). Specifically, the G-I value is positive for more cooperative genotypes, which have relatively lower individual performances in mixtures (competitive environment) but a higher performance in monoculture (cooperative environment).

## Genome-wide association mapping identifies allelic variation associated with G-I trade-off

Next, we performed genome-wide association tests for the genotypic G-I trade-off value. Genome-wide polymorphism data of our population were available through the RegMap panel, and single nucleotide polymorphism (SNP) information was available for 214,000 sites. The G-I trade-off value was significantly associated with a major effect locus on chromosome three (Fig 2A and 2B). The rarer allele (C-allele) was found in 18% of the RegMap population and was associated with lower individual but higher group performance, i.e., with increased cooperation (Fig 1C). The SNP with the strongest association resides in the center of a transposon-rich region and explained approximately 25% of the variation in the genotypic G-I trade-off values (Fig 2C). Direct mapping of untransformed data, i.e., of variation in either individual or monoculture group biomass alone, did not reveal any significant association with this locus (S2A and S2B Fig). This was expected, because such an analysis fails to separate general vigor from the G-I trade-off value that we used here to quantify cooperative properties. We did not find any other associations, likely because monoculture productivity is a complex polygenic trait and depends on a broad range of underlying processes. A more detailed genomic analysis based on a subset of 68 genotypes for which genome-wide resequencing data are available [33] revealed association signals across many polymorphisms in a region of approximately 150 kb around the identified RegMap SNP, all in high linkage disequilibrium (Fig 2B).

## Cooperative allele improves monoculture productivity at high density

To experimentally validate our finding of enhanced cooperation (i.e., of improved group performance in lines carrying the minor allele at SNP 15'294'955), we performed a stratified sampling of genotypes for plants of different sizes carrying either allele and grew these in monoculture stands of increasing planting density. Cooperative plant traits typically improve group performance especially at high planting densities [15,34,35]. Genotypes carrying the cooperation-associated allele at SNP 15'294'955 indeed exhibited superior productivity at the highest sowing density (+15% biomass, average across all genotypes; Fig 3A and 3B; ANOVA F₁,₁₀.₆ = 7.5, *P* = 0.02), despite slightly lower individual performance across mixtures in the competition experiment (−4% biomass in mixtures; Fig 3B). As anticipated, group performance gains were more pronounced at higher densities where plants experienced more severe space limitations relative to those at lower densities (Fig 3A; ANOVA F₁,₁₄.₉ = 7.0, *P* = 0.019 for allele × ground area per individual). These results demonstrate that the molecular method presented here is able to predict group-level features that cannot be deduced from individual-level properties and that these allow improving monoculture stand productivity.

## Reduced root allocation is associated with enhanced cooperation

We next explored the relationship between the allele for cooperation and phenotypic traits to identify potential mechanism underlying enhanced cooperation. We measured two traits that
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Allelic variation at a major effect locus affects the group-vs.-individual performance (G-I) trade-off in *A. thaliana*. (A) Manhattan plots of genome-wide association tests for variation in the G-I trade-off, based on the 250k SNP chip data. The genotypic G-I trade-off value is the distance from the overall trend between group and individual performance in monoculture and mixtures, respectively (inset). (B) Zoom in on a segment of chromosome 3, showing Manhattan plots of either an association analysis using SNP chip polymorphisms (top), or, for a subset of 68 genotypes, genome-wide resequencing polymorphisms (bottom). Models of protein-coding genes are drawn as boxes below, on either + (upper) or–(lower) strand. (C) Association of variation at SNP 15'294'955 and the G-I trade-off. Error bars denote means ± SEM. Data available at https://zenodo.org/record/6983283, file competition.csv.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001842.g002

characterize growth and resource-acquisition strategies of the genotypes in our panel. Specifically, we measured rosette diameter, because rosette lateral size affects light harvesting, and monoculture root-to-shoot ratio as an indicator of relative investment into belowground nutrient acquisition (Materials and methods). We further included two publicly available phenotypic traits into our analysis [36], namely flowering time in the field and vegetative growth rate. Genotypes carrying the allele for cooperation had significantly lower root-to-shoot ratios (ANOVA F₁,₉₅ = 5.13, *P* = 0.026; effect size: 19% lower root-to-shoot ratio), and these were significantly negatively associated with their G-I trade-off value (ANOVA F₁,₉₅ = 18.4, *P* < 0.001; Fig 4A and 4B). None of the other traits examined showed such an association (Figs 4A and S3). We confirmed that this pattern of lower root-to-shoot ratios in cooperative genotypes holds across environments by conducting an independent experiment, where we explored this effect in both monocultures and in isolated individual plants, and in a different soil type (Figs
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Genotypes carrying the allele for cooperation exhibit superior monoculture performances at high density. (A) Changes in monoculture biomass of genotypes carrying either the cooperation-associated allele (C-allele, red) or the alternative allele (A-allele, blue) across a realized planting density gradient. Lines show linear regression estimates ± SEM. Arrowheads show per plant areas at the sown target densities, which were not always realized due to seedling mortality. (B) Comparison of genotype's mean individual shoot biomass in mixtures versus monoculture biomass at densities of 25 plants per pot (2.56 cm² per plant). Horizontal lines: mean values across all genotypes carrying either allele. Red and blue: cooperation-associated allele (C-allele) and alternative allele (A-allele) at SNP Chr 3 15'294'955, respectively. Note the different scales of the left and the right y-axes. Data available at https://zenodo.org/record/6983283, file densitygrad.csv.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001842.g003

4C and S3). Overall, our findings indicate that relative root growth is associated with the identified genetic variation and that reduced root allocation may be part of a strategy associated with enhanced cooperation [37].

**The cooperative allele is geographically widespread and associated with pathogen resistance**

Evolutionary game theory predicts that an allele that promotes cooperation will be selected against in a natural plant population, except under special circumstances [38]. We were thus surprised that the allele for cooperation that we identified is found over a wide geographic range and at remarkably high frequency (Fig 5A; minor allele frequency 18%). Genes often have multiple functions and we therefore expected that, in natural populations, conflicting selective forces acting on such pleiotropic genes (or on genes in tight linkage) might underpin the persistence of alleles for cooperation [39]. Examining genes in the identified genomic region, we found
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Altered allocation to roots but not growth or life history is associated with increased levels of cooperation. (A) Association of allelic variation at SNP Chr3:15'294'955 with variation in traits related to different plant strategies. (B and C) Relationship between the G-I performance trade-off and plant root-to-shoot ratio in monocultures of the competition experiment shown in Fig 1 (B) or monocultures of an independent experiment (C; see also Materials and Methods and S3 Fig). Bars and regression lines show means ± SEM; * ANOVA *P* < 0.05; n.s.: not significant. Data available at https://zenodo.org/record/6983283, files competition.csv and sand.csv.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001842.g004
<figure><img src="image_8.png" /></figure>

B Emoy2

Noco2

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

Fig 5. The cooperation-associated allele exhibits a wide geographic distribution and is correlated with increased race-specific pathogen resistance. (A) Occurrence of natural *A. thaliana* genotypes carrying the allele for cooperation (red 'C') or the alternative allele (blue 'A') across sampling sites in Europe. R package "maps," version 3.4.0, was used to draw the map, which uses the database from the Natural Earth data project, and all data are in the public domain
(https://cran.r-project.org/web/packages/maps/readme/README.html) (B) Association of Chr 3 SNP 15'294'955 with resistance or susceptibility against strains of *Hyalonperonospora arabidopsidis* (Emoy2 and Noco2), based on published data [41]. Numbers indicate genotype counts. n.s. = no significant difference. (C) Association of Noco2 resistance levels with the G-I trade-off. s = susceptible, i+r = intermediately and fully resistant. Data available in Supplementary File of Nemri and colleagues [41].

https://doi.org/10.1371/journal.pbio.3001842.g005

*AtMIN7*, a documented regulator of both growth and defense. The *AtMIN7* protein targets pathogen effectors that suppress the plant immune response [40]. Importantly, variation at the *AtMIN7* gene is associated with resistance against *Hyalonperonospora arabidopsidis*, an obligate pathogen of *A. thaliana* causing downy mildew [41]. Analyzing this published data on resistance against different strains of *H. arabidopsidis* (incl. Noco2 and Emoy2), we found that the cooperation-associated allele indeed was highly significantly associated with partial or full resistance against strain Noco2 (Fig 5B, Fisher's exact test; *P* < 0.001). The resistance level against Noco2 also explained significant amounts of variation in the G-I trade-off value of our genotypes (ANOVA F<sub>2,79</sub> = 3.57, *P* = 0.03; Fig 5C). We therefore refer to this naturally occurring genetic variation as latent variation for cooperation, since contributions to pathogen resistance rather than cooperation has likely maintained the minor allele in *A. thaliana* populations.

## Lines carrying the cooperator allele have a competitive advantage in conditions of high disease pressure

Further evidence for the idea that the allele for cooperation is advantageous in environments with pathogens emerged from an additional experiment in which we quantified benefits of individual alleles (payoffs) in well-defined competitive interactions modeled according to the classic prisoner's dilemma. For the experiment, we isolated four different pairs of near-iso- genic lines (NILs) from Shahdara Bayreuth recombinant inbred lines [42], i.e., pairs of nearly identical genotypes, in which allelic variation was restricted to the region around the focal SNP (see Materials and methods for details). It is predicted that NILs homozygous for the Shah- dara (Sha) allele exhibit a more selfish and competitive, and lines homozygous for the Bayreuth (Bay) allele a more cooperative strategy. To separate effects of the Sha and Bay alleles from the general genetic background of the respective NILs, we replicated this design with four different pairs of NILs, i.e., in four different genetic backgrounds. In accordance with classic social dilemmas in evolutionary game theory (e.g., prisoner's dilemma) [28], we assumed the com- petitive genotypes would receive the highest payoff when growing with a cooperative neighbor (Fig 6A). However, such payoffs are difficult to quantify directly, and, therefore, we refer to payoffs as the reduction in aboveground biomass of plants competing relative to plants grown individually (e.g., a small reduction in biomass under competition = high payoff). We thus ran replicated experiments matching focal plants with neighbors so that all four possible combina- tions of Sha and Bay alleles were realized. All genotypes were also grown individually without neighbor, and the relative reduction of individual performance in competition was used to quantify the cost of growing with a specific neighbor.

Unexpectedly, the experiment was invaded by a parasite that affected genotypes and allele carriers differently and influenced competitive outcomes. Over the course of the experiment, some plants showed early leaf blotching, and, at harvest, many showed senesced and wilted leaves (Figs 6B and S5). This is a phenotype typical of soil-borne diseases. The symptoms were dependent on genotype and were mainly seen in NIL-background 33RV113 and individuals carrying the competitor (Sha) allele (Fig 6C). At the same time, we detected abundant spores of the obligate plant parasite *Olpidium brassicae* in roots of all sampled pots (n = 18; S5 Fig). We took advantage of this to test the hypothesis that the allele increasing cooperation also
<figure><img src="image_12.png" /><figcaption>**Fig 6.** A trade-off between competitive dominance and disease resistance expressed in different near-isogenic lines. (A) Expected payoff matrix of interactions between plants ("cooperate" vs. "compete," as determined by carrying Bay vs. Sha alleles), with Temptation (T) > Reward (R) > Punishment (P) > Sucker's Payoff (S) (color coded). (B) Experimental outline: benefits refer to the cost reduction as measured by the average aboveground biomass of genotypes in competition compared to the biomass of the plant grown alone. A strong reduction refers to a high costs and therefore low benefit, a weak reduction to a high benefit. (C) Pot containing a pair of near-isogenic plants isolated from the RIL parent 33RV113. The plant carrying the competitive Sha allele exhibits strong signs of leaf wilting (arrows). (D) Summary of the competition experiment, separated by genetic background (33RV085, 33RV142, 33RV192, 33RV192), showing disease symptoms (number of plants exhibiting leaf blotching; root browning score) and payoff matrices for all genetic backgrounds and allele combinations. The values in payoff matrix are the average biomass reduction ± SEM in brackets. Data available at https://zenodo.org/record/6983283, file `Payoff_longFormat.csv`.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001842.g006

conferred pathogen resistance. As a proxy for root disease expression, we visually scored root yellowing/browning, a phenotype caused by *O. brassicae* [43]. We found that this score significantly depended on genetic background (ordinal logistic regression *p*-value < 0.001 for NIL pair) and was always lower in genotypes carrying the cooperative (Bay) allele (ordinal logistic regression *p*-value = 0.001; Fig 6C).

As predicted, plants carrying the allele promoting cooperation were competitively inferior when matched with a competitive neighbor, but only in the two genetic backgrounds that showed low disease expression (Fig 6D). In the two genetic backgrounds with high disease
expression, the competitive relationship was very different: It was most pronounced in the most susceptible genetic background, the genotype carrying the allele promoting cooperation was competitively superior.

Together, our results suggest that the cooperative allele is associated with increased resistance against several diseases, which may confer a competitive advantage under high disease pressure. This effect may be caused by antagonistic pleiotropy or by different alleles at closely linked genes. Whether or not such a trade-off between competition and disease resistance affects selection in nature will have to be tested. However, it could explain why a cooperative allele is maintained in a natural population and over such a wide geographic range.

## Discussion

Here we used specifically designed competition experiments, combined with genetic association mapping, to identify (1) plant genotypes that perform well when grown in high-density stands; and (2) potential genes and alleles that underlie these cooperative growth strategies. Interestingly, in *A. thaliana*, we detected a specific allele that is associated with increased growth in high-density stands, representing a strategy we here refer to as “cooperation.” The importance of cooperation among individuals for agronomic production has long been recognized. Animal breeding has focused more explicitly on cooperation and social strategies [20,44]. Domesticated animals are often kept at relatively high densities, which can trigger aggressive behavior among individuals, resulting in high levels of stress, injuries, and increased mortality, all of which are detrimental to animal welfare and productivity. A breeding objective thus has been to obtain relatively docile, cooperative animals that limit adverse interactions when kept at high densities [45]. Conceptually similar considerations apply to plant communities. However, in plants, the traits that underlie cooperation are less easily recognized. In general, such traits may be related to resource sharing in the sense that resources remain accessible to the individuals that most efficiently convert these into biomass or yield. In other words, any excessive investment into competition in the presence of neighbors undermines what is, from an agronomic perspective, an efficient use of resources at the group level [25,26,46–50]. In plants, cooperative growth strategies thus may be characterized by a relative insensitivity to the presence of neighbors, or even a reduction in the development of competitive structures. Indeed, improved group functioning was observed in crop experiments where genotypes exhibited more vertical leaf angles [51], where leaf area [52] or branch number were reduced [53], or where belowground competition was experimentally prevented [46,50]. Cooperative responses may also be plastic and expressed predominantly in monocultures; for example, kin recognition studies have shown that plants competing with close relatives reduced root allocation relative to when they competed with strangers [54]. Common to all these experiments is the finding that restricting competitive interference between neighbors, either experimentally or by using genotypes with corresponding traits, is beneficial for group yield, especially at high planting densities. In our study, we observed a lower relative investment into root growth in cooperative genotypes, which is in line with this consideration.

Cooperation, although beneficial for group yield, is typically also associated with individual fitness costs in noncooperative environments, i.e., under strong competition [28–30]. Thus, addressing this G-I trade-off is key to identifying and using cooperative genotypes. In breeding, different strategies have been developed that achieve this goal. First, cooperation can be selected for by anticipating corresponding ideotypes possessing traits expected to foster cooperative strategies. This has been very successfully accomplished during the Green Revolution [6,12,13,55] and led to the recognition of the relevance of a G-I trade-off in crop improvement [14,18]. Second, cooperative strategies can be directly selected for at the group level [17,22,56]. For example,
this is often done in field trials with maize, whereby identical hybrid genotypes are planted and evaluated in high-density rows and subsequently in larger field plots. Such a selection at the group level also promotes traits that benefit group yield, for example, reduced tassel size, more vertical leaves, and reduced tillering [17,29,57]. Another idea is to incorporate kin selection schemes into the breeding process (e.g., in the early cycles of breeding programs where individual genotypes compete with cosegregants), which could also select for crops that exhibit cooperative traits [24,56,58]. Though it is currently unclear how widespread kin recognition is in plants [59], kin recognition studies in wild populations are another avenue that could promote the identification of cooperative traits in plants [22]. Finally, genomic predictions are increasingly adopted to predict a genotype's merit purely based on its genome, whereby prediction models are trained on yield variation of high-density stands of pure varieties. These approaches have arguably contributed to historical and contemporary yield advances and implicitly or explicitly include elements that selected to some extent for cooperative traits [22]. However, they also rely on genetic variation in the relevant traits to be present in the breeding populations. The method we present here—while motivated by the same broader topic—differs in important respects. First, it directly assesses the G-I trade-off rather than group yield. By exposing genotypes to two “social environments” that represent different ends of the cooperation-competition continuum, genotypes expressing cooperative traits can more easily be distinguished. This contrasts with other methods that either focus only on groups (e.g., row or plot-level selection) or use an environment with mixed neighborhoods (e.g., typical setting in quantitative genetic analyses of indirect genetic effects; [20,60,61]). Second, our method systematically decomposes group yield into two effects: general vigor and cooperation. This separation is critical because alleles promoting cooperation may be found at low frequency and in genotypes with low vigor. A prominent example is the recessive semidwarf allele *sd1* in rice. It was originally found in the Taiwanese Dee-Geo-Woo-Gen variety, which showed only average yields in the tropics, but was instrumental for boosting tropical rice yields [9, 10, 18]. Third, the fact that we use an empirical, heuristic relationship to correct for general vigor implicitly accounts for nonlinear relationships between genetics and productivity, which is typically difficult to address with other methods.

In combination with genetic mapping, our method provides an “agnostic” approach to discover latent genetic variation for cooperation. In our proof-of-concept study, we identified a genetic variant of large effect. This was surprising, since traits underlying individual growth and yield or environmental adaptations, such as drought tolerance, are often highly polygenic [62–64]. On the other hand, there is now increasing evidence that different plant-plant interactions can be affected by few genetic variants of large effects ([65,66]; this study), maybe resembling the pattern more generally found in plant-parasite interactions [67]. Regardless of the reason, the existence of such alleles is important given that many biotechnological methods, including genome editing, typically rely on the availability of traits that are influenced by a single or few genes. Therefore, we propose that biotechnological and breeding research should extend their focus beyond the performance of individuals, often grown as single plants in pots on “conveyor belt systems,” towards the performance of groups and the relevant competitive interactions that drive it.

In this work, we focused—primarily for the sake of simplicity—on a model population of *A. thaliana* and on aboveground biomass (which is a target trait for some crops, such as forage grasses). This is a first case study using the proposed methodology, and a larger-scale systematic search may reveal alleles with comparable effects in crops or their wild relatives. Here, we found that the allele for cooperation among plants remained relatively frequent (18%) in natural populations, likely because it conferred resistance to diseases to the individuals carrying this allele, and thereby higher fitness under high disease pressure. This finding raises hopes that crops and wild plants may also hold hidden cooperative traits that have persisted despite their
disadvantages in competitive situations. Once identified, such latent variation in cooperation
could rapidly be co-opted in marker-assisted breeding programs or through biotechnology. At
a more fundamental level, the finding that large-effect genetic variants for cooperation are
maintained in natural populations leads to the intriguing thought that social traits could arise as
evolutionary exaptations, i.e., by co-option of an existing trait unrelated to cooperation [68].

**Materials and methods**

**Plant material**

The natural *A. thaliana* genotypes used (S1 Data) are a subset of the RegMap population [31] for which a comprehensive list of traits has been collected [36]. For pairs of NILs isolated from four Bay-Sha RILs, we confirmed homozygosity for a genetic polymorphism upstream of the focal SNP (primer F: 5'-TGAGAGAGAGCTGATGATGGATG-3'; primer R: 5'-CGCCTTGA TTGACACAGATTC-3'; approx. 100 bp deletion in Sha sequence at position 15'333'951) and a genetic polymorphism downstream of the focal SNP by the use of PCR markers (markers primer F: 5'-GCAAGAGGGAGCTAAAGAACAG-3'; primer R: 5'-GCCCTTATCGCCAT GAACTG-3'; approx. 50 bp deletion in Sha sequence at position 14'914'025). These deletions in the Sha sequence were predicted by the Polymorph tool (http://polymorph.weigelworld.org). All eight genomes were also resequenced on the DNBseq platform (BGI, Hongkong) with a minimum of three Gb sequences per genome. The reads were aligned to the *Arabidopsis* Col-0 reference sequence with the use of the bwa software (version 0.7.16a) [69]; read sorting and variant calling was performed using samtools (version 1.5) [70]. Whole-genome reconstruction was performed as described previously [65], following the approach developed by Xie and colleagues [71]. Genome reconstructions were in agreement with the published marker maps of the corresponding RILs ([42]; S5A Fig). We confirmed genetic differentiation of each of the four NIL pairs at the locus of interest. Two pairs were otherwise completely isogenic; the two other pairs showed some residual genetic differences on other chromosomes (S5B Fig).

**Experimental design**

**Competition experiment:** Pairs of individual plants were grown in small pots (6 × 6 × 5.5 cm) in a factorial design in which the 97 genotypes of the panel were each grown together with one of ten tester genotypes (Bay-0, C24, Col-0, Cvi-0, Ler-1, Sav-0, Sf-2, Shahdara, St-0, Uk-1). These tester genotypes were chosen because they constitute the parents of different publicly available recombinant inbred line populations, i.e., a genetic resource that allows for certain genetic experiments (for example, isolating NILs from heterogeneous inbred families; see Fig 6). The tester genotypes were a subset of the panel, but tester genotypes not part of the original panel would have worked equally well. Each genotype was further grown in a monoculture of two individuals. Each genotype composition was replicated twice, in separate blocks. In the second block, not enough seeds were available for line LP-2-6; in the second block, we therefore replaced this genotype by Kn-0. We thus effectively had 98 genotypes grown with the ten tester genotypes. Overall, the experiment consisted of 2,134 pots with two plants each. Each tester line was also grown as individual plant, once per block. Pots containing single plants, and pots in which one plant died at the seedling stage, were excluded from data analyses.

**Density gradient:** To test for decreased self-inhibition in plants carrying the cooperation-associated allele, genotype monocultures with different individual densities (3 × 3, 4 × 4, or 5 × 5 plants on 8 × 8 cm area) were constructed. Six genotypes that carried the cooperation-associated allele (Bor-4, Est-1, Mt-0, Ra-0, Sav-0, Wa-1) but varied in average individual performances were paired with seven genotypes (An-1, Br-0, Can-0, Kondara, Nfa-10, Shahdara, St-0) that carried the alternative allele but had similar average individual performances. By
matching genotypes by size, we were able to separate the size dependence of self-inhibition
from the allele effect.

**Preparation of plant material and growth conditions**

Competition experiment: Seeds of all genotypes were sown directly onto soil (four parts Ein-
heitserde ED73, Gebrüder Patzer, Germany; one-part quartz sand) in February 2016. For ger-
mination, the pots of size 6 × 6 × 5.5 cm from a given block were randomly placed into trays
covered with plastic lids. The two plants per pot were established 3 to 4 cm apart by sowing 5
to 20 seeds per position and, once seeds had germinated, removing surplus seedlings. Blocks 1
and 2 were sown on February 17 and 18, 2016, respectively. Plants were grown with a photope-
riod of 14 hours, providing additional light when necessary. Daytime and nighttime tempera-
tures were maintained around 20 to 25°C and 16 to 20°C, respectively. Trays were randomly
rearranged within the greenhouse every 3 to 5 days. After 5 to 5.5 weeks, pots were transferred
from trays onto three tables with automated watering and randomly rearranged every week.
Flowering shoots of individual plants were tied to wooden sticks when they grew taller than
approx. 10 cm. All plants were harvested on April 14 (Block 1) and April 15 (Block 2), 2016,
i.e., approx. 8 weeks after sowing. Each plant was cut below the rosette and individually dried
at 65°C for 4 to 5 days and then stored at room temperature until weighing. Roots from a pot
were washed on a metal sieve, and total root mass determined after drying at 65°C for 4 days.
Flowering time was determined by checking every 2 to 3 days whether flowering bolts were
present that exceeded 0.5 cm height.

Density gradient: Monocultures were sown in pots of 9 × 9 × 10 cm (inner pot diameter approximately 8 × 8 cm) at densities of either 9, 16, or 25 plants per pot. Soil and growth conditions were as described above for the competition experiment. Because we observed some seedling mortality early in the experiment, realized planting density was determined from photographs taken 27 days after sowing, i.e., at a time when competition still was very limited. Aboveground biomass was harvested, dried, and weighed 54 days after sowing.

Independent biomass allocation measurements: For an independent assessment of root-to-
shoot biomass ratios in the studied natural genotypes, 80 genotypes that were used in the main
competition experiment were grown for 43 days either as single plants or as monoculture
(consisting of four plants per pot) and in pots of 7 × 7 × 8 cm size on a mixture of one-part
ED73 and four-parts quartz sand. The measurements were performed as described above.

Competition experiments with NIL pairs: We used four NIL pairs derived from 33RV085, 33RV113, 33RV142, and 33RV192. Each pair consisted of one line carrying the Bay allele and one line carrying the Sha allele. For each pair, we grew 16 monocultures per line and 32 mixtures of the two lines, all consisting of two plants per pot. We additionally grew each line as single plants per pot, also replicated 16 times. Seeds were directly sown into 6 × 6 × 5.5 cm pots filled with a 4:1 parts ED73 soil:quartz sand. Leaf blotching and excessive senescence were noticed midway through the experiment and scored 31 and 50 days after sowing. Scoring was done blindly, i.e., without information about pot or genotype identity. Shoot and root biomass were harvested 57 days after sowing, roots washed, and all biomass dried and weighed as described above. Root browning scores were determined visually by categorizing root browning into three ordinary levels (S5 Fig). Roots were stained using a 1% cotton blue in 80% lactic acid solution and examined using a Leica microscope for a visual inspection of potential root pathogens.

**Statistical analyses**

All statistical analyses were performed using the statistical software R version 3.4.1 (http://r-project.org). The block-adjusted individual performance of genotypes across mixtures were
estimated using least square means from a model including just block and genotype. Monoculture biomass per individual (i.e., total average monoculture biomass divided by two) was then fitted as function of linear and quadratic forms of individual biomass, using the R-function "lm." The G-I trade-off value was determined as orthogonal distance by determining the point in the quadratic heuristic that was closest to the respective point by nonlinear minimization using the R-function "nlm." The GWAS analyses were performed with easyGWAS (https://easygwas.ethz.ch) [72], using the EMMAX algorithm [73] and using SNPs from the 250k SNP chip (http://bergelson.uchicago.edu/) or the 1001 genomes project (http://1001genomes.org/) [33]. SNPs with a minor allele frequency below 5% were removed.

To test for the dependency of the G-I trade-off value on relative root allocation, we linearly regressed the trade-off value against root-to-shoot ratio and tested for significance of the slope term. We did not transform the axes because the relationship was sufficiently linear. Note that analysis on a log-log scale, as frequently used in allometric analyses, was not appropriate here because the G-I trade-off value was centered around zero (i.e., not systematically positive).

For the density gradient experiment, productivity was modeled in dependence of the fixed terms area per individual, allele, plus their interaction (we chose area per individual (available space) instead of density (the inverse of available space) because the relationship with density was nonlinear). The corresponding random terms were genotype, and the interaction between genotype and area per individual. Individual performance in mixture of the previous experiment was included as a covariate in the model, because it was used in the design of the experiment for stratified sampling of genotypes according to size. The realized densities deviated from sown densities because of a relatively high initial mortality. Therefore, we instead used densities determined from photographs of each pot that were made midway through the experiment. Two pots were removed from the analysis because realized densities were much higher than planted densities, probably because they accidentally had not been thinned to the intended densities.

In the NIL pair competition experiment, the payoff for a given combination of strategies was determined as shoot mass of a line grown as individual plant minus its shoot mass under competition, i.e., the payoff quantified the shot biomass reduction induced by the presence of particular neighbors. The standard error of the payoff $\sigma_{payoff}$ was calculated by error propagation as $\sqrt{(\sigma_{single plant}^2 + \sigma_{individual | neighbour}^2)}$.

Root browning scores were analyzed using an ordinal logit regression as implemented in the polr-function of the R package MASS, with ordinal root scores as independent variable and community type (single plant or monoculture), genetic background (NIL pair), and allele (Sha or Bay) as explanatory variables. Mixtures were excluded from the analysis of root browning scores.

## Supporting information

**S1 Fig. Experimental setup of the competition experiment.** (A-C) Photos show the experiment at sowing (A), midway through the experiment (B), and at harvest day (C). (TIFF)

**S2 Fig. Association tests for variation in (A) average individual performance across mixtures or (B) average monoculture performance.** Data available at https://zenodo.org/record/6983283, file competition.csv. (TIFF)

**S3 Fig. Associations of SNP Chr3:15'294'955 with phenotypic variation in traits related to plant strategies.** (A and B) Shoot-to-root ratios for genotypes grown in monocultures (A) or as individual plants (B) in an independent experiment and on sand-rich soil are shown. (C) Published data of genotypic means in flowering time in the field [36]. Bars show mean ± SEM.
** = ANOVA *p* < 0.01; * = ANOVA *p* < 0.05; n.s., not significant. Data available at https://zenodo.org/record/6983283, filesand.csv, and at https://arapheno.1001genomes.org/phenotype/86/.
(TIFF)

**S4 Fig. Decreased self-inhibition of genotypes carrying the cooperation-associated allele along a density gradient.** Lines represent reaction norms or genotypes carrying different alleles at Chr 3 SNP 15'294'955; red lines: cooperator allele carriers; blue lines: competitor allele carriers. Data available at https://zenodo.org/record/6983283, file densitygrad.csv.
(TIFF)

**S5 Fig. A trade-off between competitive dominance and disease resistance.** (A) Genome reconstructions of two near-isogenic lines (NILs) from the 33RV085 parental RIL. Viterbi paths of each NIL are shown in red (Bay-allele at Chr3:15'294'955) and blue (Sha-allele), together with published marker data of the parent (dashed line). (B) Genetic differences identified in whole-genome comparisons between the pairs of near-isogenic lines. All NIL pairs differ from each other at Chr3:15'294'955 (light blue: "locus of interest"). (C) Earliest sign of shoot-level disease expression in line 33RV113 carrying the Sha-allele. (D) Microscopic image of dried root samples containing resting spores of the obligate parasite *Olpidium brassicae* (arrows). Scale bar = 50 µm. (E) Visual scoring scheme of root browning for washed and dried root samples from the competition experiment. Scoring was performed fully blinded to the sample ID or genotype.
(TIFF)

**S1 Data. List of *A. thaliana* genotypes used in the study, their estimated productivities across mixtures and monocultures, and measured trait values.**
(CSV)

**Acknowledgments**

We thank Bernhard Schmid (UZH) and Andrea Patocchi (Agroscope) for support and helpful discussions, Matthias Lutz (Agroscope) for help in identifying *Olpidium brassicae* in root samples, Cyrille Violle (CEFE) for helpful comments on the manuscript, Matthias Philipp, Daniel Trujillo, and Mariela Soto Araya for help with sowing and harvesting the competition experiment, and Matthias Furler for technical support in the greenhouse.

**Author Contributions**

**Conceptualization:** Samuel E. Wuest, Nuno D. Pires, Ueli Grossniklaus, Pascal A. Niklaus.

**Data curation:** Samuel E. Wuest.

**Formal analysis:** Samuel E. Wuest, Pascal A. Niklaus.

**Funding acquisition:** Samuel E. Wuest, Ueli Grossniklaus.

**Investigation:** Samuel E. Wuest, Nuno D. Pires, Shan Luo, Pascal A. Niklaus.

**Methodology:** Samuel E. Wuest, Pascal A. Niklaus.

**Resources:** Nuno D. Pires, Francois Vasseur, Julie Messier, Ueli Grossniklaus.

**Visualization:** Samuel E. Wuest, Pascal A. Niklaus.

**Writing – original draft:** Samuel E. Wuest, Pascal A. Niklaus.
**Writing - review & editing:** Samuel E. Wuest, Nuno D. Pires, Shan Luo, Francois Vasseur, Julie Messier, Ueli Grossniklaus, Pascal A. Niklaus.

## References

1. Watson A, Ghosh S, Williams MJ, Cuddy WS, Simmonds J, Rey MD, et al. Speed breeding is a powerful tool to accelerate crop research and breeding. Nat Plants. 2018; 4:23–29. https://doi.org/10.1038/s41477-017-0083-8 PMID: 29292376

2. Crossa J, Pérez-Rodríguez P, Cuevas J, Montesinos-López O, Jarquín D, de los Campos G, et al. Genomic Selection in Plant Breeding: Methods, Models, and Perspectives. Trends Plant Sci. 2017; 22:961–975. https://doi.org/10.1016/j.tplants.2017.08.011 PMID: 28965742

3. Furbank RT, Tester M. Phenomics—technologies to relieve the phenotyping bottleneck. Trends Plant Sci. 2011; 16:635–644. https://doi.org/10.1016/j.tplants.2011.09.005 PMID: 22074787

4. Ray DK, Ramankutty N, Mueller ND, West PC, Foley JA. Recent patterns of crop yield growth and stagnation. Nat Commun. 2012; 3. https://doi.org/10.1038/ncomms2296 PMID: 23250423

5. Ray DK, Mueller ND, West PC, Foley JA. Yield Trends Are Insufficient to Double Global Crop Production by 2050. PLoS ONE. 2013; 8:e66428. https://doi.org/10.1371/journal.pone.0066428 PMID: 23840465

6. Borlaug NE. Wheat Breeding and its Impact on World Food Supply. Third International Wheat Genetics Symposium. 1968. pp. 1–36.

7. Khush GS. Green revolution: the way forward. Nat Rev Genet. 2001; 2:815–821. https://doi.org/10.1038/35093585 PMID: 11584298

8. De Datta SK, Tauro AC, Balaoing SN. Effect of Plant Type and Nitrogen Level on the Growth Characteristics and Grain Yield of Indica Rice in the Tropics1. Agron J. 1968; 60:643–647. https://doi.org/10.2134/agronj1968.00021962006000060017x

9. Hedden P. The genes of the Green Revolution. Trends Genet. 2003; 19:5–9. https://doi.org/10.1016/s0168-9525(02)00009-4 PMID: 12493241

10. Spielmeyer W, Ellis MH, Chandler PM. Semidwarf (sd-1), “green revolution” rice, contains a defective gibberellin 20-oxidase gene. Proc Natl Acad Sci. 2002; 99:9043–9048. https://doi.org/10.1073/pnas.132266399 PMID: 12077303

11. Donald CM. Competitive plants, communal plants, and yield in wheat crops. In: Evans LT, Peacock WJ, editors. Wheat Science—Today and Tomorrow. Cambridge University Press; 1981. p. 223–247.

12. Donald CM. The breeding of crop ideotypes. Euphytica. 1968; 17:385–403. https://doi.org/10.1007/BF00056241

13. Jennings PR. Plant Type as a Rice Breeding Objective. Crop Sci. 1964; 4:13–15. https://doi.org/10.2135/cropsci1964.0011183x000400010005x

14. Jennings PR, De Jesus JJ. Studies on competition in rice I Competition in mixtures of varieties. Evolution. 1968; 22:119–124. https://doi.org/10.1111/j.1558-5646.1968.tb03455.x

15. Duvick DN, Smith JSC, Cooper M. Long-Term Selection in a Commercial Hybrid Maize Breeding Program. In: Janick J, editor. Plant Breeding Reviews Part 2 Long Term Selection: Crops, Animals and Bacteria, Vol 24. New York: John Wiley & Sons; 2004. p. 109–151. https://doi.org/10.1002/9780470650288.ch4

16. Tsunoda S. A Developmental Analysis of Yielding Ability in Varieties of Field Crops: II. The assimilation-system of plants as affected by the form, direction and arrangement of single leaves. Japanese J Breed. 1959; 9:237–244. https://doi.org/10.1270/jsbbs1951.9.237

17. Weiner J. Looking in the Wrong Direction for Higher-Yielding Crop Genotypes. Trends Plant Sci. 2019; 19:S1360–1385. https://doi.org/10.1016/j.tplants.2019.07.001 PMID: 31358472

18. Jennings PR, Herrera RM. Studies on competition in rice II: Competition in segregating populations. Evolution. 1968; 22:332–336. https://doi.org/10.1111/j.1558-5646.1968.tb05901.x

19. Peng S, Khush GS, Virk P, Tang Q, Zou Y. Progress in ideotype breeding to increase rice yield potential. F Crop Res. 2008; 108:32–38. https://doi.org/10.1016/j.fcr.2008.04.001

20. Wade MJ, Bijma P, Ellen ED, Muir W. Group selection and social evolution in domesticated animals. Evol Appl. 2010; 3:453–465. https://doi.org/10.1111/j.1752-4571.2010.00147.x PMID: 25567938

21. Goodnight C. The influence of environmental variation on group and individual selection in a cress. Evolution. 1985; 39:545–558. https://doi.org/10.1111/j.1558-5646.1985.tb00394.x PMID: 28561972

22. Murphy GP, Swanton CJ, Van Acker RC, Dudley SA. Kin recognition, multilevel selection and altruism in crop sustainability. J Ecol. 2017; 105:930–934. https://doi.org/10.1111/1365-2745.12787
23. Dudley SA. Plant cooperation. AoB Plants. 2015; 7:plv113. https://doi.org/10.1093/aobpla/plv113
PMID: 26410832

24. File AL, Murphy GP, Dudley SA. Fitness consequences of plants growing with siblings: Reconciling kin selection, niche partitioning and competitive ability. Proc R Soc B Biol Sci. 2012; 279:209–218. https://doi.org/10.1098/rspb.2011.1995 PMID: 22072602

25. Falster DS, Westoby M. Plant height and evolutionary games. Trends Ecol Evol. 2003; 18:337–343.
https://doi.org/10.1016/S0169-5347(03)00061-2

26. Anten NPR, Vermeulen PJ. Tragedies and Crops: Understanding Natural Selection To Improve Crop-
ping Systems. Trends Ecol Evol. 2016; 31:429–439. https://doi.org/10.1016/j.tree.2016.02.010 PMID:
27012675

27. Rankin DJ, Bargum K, Kokko H. The tragedy of the commons in evolutionary biology. Trends Ecol Evol.
2007; 22:643–51. https://doi.org/10.1016/j.tree.2007.07.009 PMID: 17981363

28. Axelrod R, Hamilton WD. The evolution of cooperation. Science. 1981:1390–1396. https://doi.org/10.
1126/science.7466396 PMID: 7466396

29. Denison RF. Darwinian Agriculture. Princeton, NJ: Princeton Univ. Press; 2012. https://doi.org/10.
23943/princeton/9780691139500.001.0001

30. Weiner J, Du YL, Zhang C, Qin XL, Li FM. Evolutionary agroecology: individual fitness and population
yield in wheat (*Triticum aestivum*). Ecology. 2017; 98:2261–2266. https://doi.org/10.1002/ecy.1934
PMID: 28783218

31. Horton MW, Hancock AM, Huang YS, Toomajian C, Atwell S, Auton A, et al. Genome-wide patterns of
genetic variation in worldwide *Arabidopsis thaliana* accessions from the RegMap panel. Nat Genet.
2012; 44:212–216. https://doi.org/10.1038/ng.1042 PMID: 22231484

32. Laird RA, Aarssen LW. Size inequality and the tragedy of the commons phenomenon in plant competi-
tion. Plant Ecol. 2005; 179:127–131. https://doi.org/10.1007/s11258-004-6737-4

33. Alonso-Blanco C, Andrade J, Becker C, Bemm F, Bergelson J, Borgwardt KM, et al. 1,135 genomes
reveal the global pattern of polymorphism in *Arabidopsis thaliana*. Cell. 2016; 166:481–491. https://doi.
org/10.1016/j.cell.2016.05.063 PMID: 27293186

34. York LM, Galindo-Castañeda T, Schussler JR, Lynch JP. Evolution of US maize (*Zea mays* L.) root
architectural and anatomical phenes over the past 100 years corresponds to increased tolerance of
nitrogen stress. J Exp Bot. 2015; 66:2347–2358. https://doi.org/10.1093/jxb/erv074 PMID: 25795737

35. Zhu YH, Weiner J, Yu MX, Li FM. Evolutionary agroecology: Trends in root architecture during wheat
breeding. Evol Appl. 2018; 00:1–11. https://doi.org/10.1111/eva.12749 PMID: 30976306

36. Atwell S, Huang YS, Vilhjálmsson BJ, Willems G, Horton M, Li Y, et al. Genome-wide association study
of 107 phenotypes in *Arabidopsis thaliana* inbred lines. Nature. 2010; 465:627–631. https://doi.org/10.
1038/nature08800 PMID: 20336072

37. Biedrzycki ML, Jilany TA, Dudley SA, Bais HP. Root exudates mediate kin recognition in plants. Com-
mun Integr Biol. 2010; 3:28–35. https://doi.org/10.4161/cib.3.1.10118 PMID: 20539778

38. Nowak MA. Five rules for the evolution of cooperation. Science. 2006; 314:1560–1563. https://doi.org/
10.1126/science.1133755 PMID: 17158317

39. Todesco M, Balasubramanian S, Hu TT, Traw MB, Horton M, Epple P, et al. Natural allelic variation
underlying a major fitness trade-off in *Arabidopsis thaliana*. Nature. 2010; 465:632–636. https://doi.org/
10.1038/nature09083 PMID: 20520716

40. Nomura K, DebRoy S, Lee YH, Pumplin N, Jones J, He SY. A bacterial virulence protein suppresses
host innate immunity to cause plant disease. Science. 2006; 313:220–223. https://doi.org/10.1126/
science.1129523 PMID: 16840699

41. Nemri A, Atwell S, Tarone AM, Huang YS, Zhao K, Studholme DJ, et al. Genome-wide survey of *Arabi-
dopsis* natural variation in downy mildew resistance using combined association and linkage mapping.
Proc Natl Acad Sci. 2010; 107:10302–10307. https://doi.org/10.1073/pnas.0913160107 PMID:
20479233

42. Loudet O, Chaillou S, Camilleri C, Bouchez D, Daniel-Vedele F. Bay-0 x Shahdara recombinant inbred
line population: A powerful tool for the genetic dissection of complex traits in *Arabidopsis*. Theor Appl
Genet. 2002; 104:1173–1184. https://doi.org/10.1007/s00122-001-0825-9 PMID: 12582628

43. Van Der Meer JHH. The rhizoctonia- and oplidium disease of cauliflower seedlings. Tijdschr Over Plan-
tenziekten. 1926; 32:235–242. https://doi.org/10.1007/bf02812706

44. Muir WM. Group selection for adaptation to multiple-hen cages: Selection program and direct
responses. Poult Sci. 1996; 75:447–458. https://doi.org/10.3382/ps.0750447 PMID: 8786932

45. Friedrich J, Brand B, Schwerin M. Genetics of cattle temperament and its impact on livestock production
and breeding—a review. Arch Anim Breed. 2015; 58:13–21. https://doi.org/10.5194/aab-58-13-2015
46. Gersani M, Brown JS, O'Brien EE, Maina GM, Abramsky Z. Tragedy of the commons as a result of root competition. J Ecol. 2001; 89:660–669. https://doi.org/10.1046/j.0022-0477.2001.00609.x

47. Cabal C, Martínez-García R, De Castro Aguilar A, Valladares F, Pacala SW. The exploitative segregation of plant roots. Science. 2020; 370:1197–1199. https://doi.org/10.1126/science.aba9877 PMID: 33273098

48. Zhang DY, Sun GJ, Jiang XH. Donald's ideotype and growth redundancy: A game theoretical analysis. F Crop Res. 1999; 61:179–187. https://doi.org/10.1016/S0378-4290(98)00156-7

49. Zhu L, Zhang DY. Donald's Ideotype and Growth Redundancy: A Pot Experimental Test Using an Old and a Modern Spring Wheat Cultivar. PLoS ONE. 2013; 8:e70006. https://doi.org/10.1371/journal.pone.0070006 PMID: 23936133

50. Zhu Y-H, Weiner J, Li F-M. Root proliferation in response to neighbouring roots in wheat (Triticum aestivum). Basic Appl Ecol. 2019; 39:10–14. https://doi.org/10.1016/j.baae.2019.07.001

51. Tian J, Wang C, Xia J, Wu L, Xu G, Wu W, et al. Teosinte ligule allele narrows plant architecture and enhances high-density maize yields. Science. 2019; 365:658–664. https://doi.org/10.1126/science.aax5482 PMID: 31416957

52. Srinivasan V, Kumar P, Long SP. Decreasing, not increasing, leaf area will raise crop yields under global atmospheric change. Glob Chang Biol. 2017; 23:1626–163. https://doi.org/10.1111/gcb.13526 PMID: 27860122

53. Kokubun M. Design and examination of soybean ideotypes. Japan Agric Res Q. 1987; 21:237–243.

54. Dudley SA, File AL. Kin recognition in an annual plant. Biol Lett. 2007; 3:435–438. https://doi.org/10.1098/rsbl.2007.0232 PMID: 17567552

55. Vogel OA, Craddock JC, Muir CE, Everson EH, Rohde CR. Semidwarf Growth Habit in Winter Wheat Improvement for the Pacific Northwest1. Agron J. 1956; 48:76–78. https://doi.org/10.2134/agronj1956.00021962004800020009x

56. Biernaskie JM. Kin selection theory and the design of cooperative crops. Evol Appl. 2022; 00:1–10. https://doi.org/10.1111/eva.13418 PMID: 36330299

57. Denison RF. Evolutionary tradeoffs as opportunities to improve yield potential. F Crop Res. 2015; 182:3–8. https://doi.org/10.1016/j.fcr.2015.04.004

58. Montazeaud G, Rousset F, Fort F, Violle C, Fréville H, Gandon S. Farming plant cooperation in crops. Proceedings Biol Sci. 2020; 287:20191290. https://doi.org/10.1098/rspb.2019.1290 PMID: 31964305

59. Milla R, Forero DM, Escudero A, Iriondo JM. Growing with siblings: A common ground for cooperation or for fiercer competition among plants? Proc R Soc B Biol Sci. 2009; 276. https://doi.org/10.1098/rspb.2009.0369 PMID: 19403541

60. Costa E Silva J, Potts BM, Gilmour AR, Kerr RJ. Genetic-based interactions among tree neighbors: Identification of the most influential neighbors, and estimation of correlations among direct and indirect genetic effects for leaf disease and growth in Eucalyptus globulus. Heredity. 2017; 119:125–135. https://doi.org/10.1038/hdy.2017.25 PMID: 28561806

61. Camerlink I, Ursinus WW, Bijma P, Kemp B, Bolhuis JE. Indirect Genetic Effects for Growth Rate in Domestic Pigs Alter Aggressive and Manipulative Biting Behaviour. Behav Genet. 2015; 45:117–126. https://doi.org/10.1007/s10519-014-9671-9 PMID: 25227986

62. Lynch M, Walsh B. Genetics and Analysis of Quantitative Traits. Sunderland, Massachusetts: Sinauer; 1998. https://doi.org/10.1086/318209

63. Boyle EA, Li YI, Pritchard JK. An Expanded View of Complex Traits: From Polygenic to Omnigenic. Cell. 2017; 169:1177–1186. https://doi.org/10.1016/j.cell.2017.05.038 PMID: 28622505

64. Wieters B, Steige KA, He F, Koch EM, Ramos-Onsins SE, Gu H, et al. Polygenic adaptation of rosette growth in *Arabidopsis thaliana*. PLoS Genet. 2021; 17:e1008748. https://doi.org/10.1371/journal.pgen.1008748 PMID: 33493157

65. Wuest SE, Niklaus PA. A plant biodiversity effect resolved to a single chromosomal region. Nat Ecol Evol. 2018; 2:1933–1939. https://doi.org/10.1038/s41559-018-0708-y PMID: 30397303

66. Montazeaud G, Flutre T, Ballini E, Morel J-B, David J, Girodolle J, et al. From cultivar mixtures to allelic mixtures: opposite effects of allelic richness between genotypes and genotype richness in wheat. New Phytol. 2022; 233:2573–2584. https://doi.org/10.1111/nph.17915 PMID: 35081666

67. Wilfert L, Schmid-Hempel P. The genetic architecture of susceptibility to parasites. BMC Evol Biol. 2008; 8:187. https://doi.org/10.1186/1471-2148-8-187 PMID: 18590517

68. Gould SJ, Vrba ES. Exaptation—a Missing Term in the Science of Form. Paleobiology. 1982; 8:4–15. https://doi.org/10.1017/S0094837300004310

69. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25:1754–1760. https://doi.org/10.1093/bioinformatics/btp324 PMID: 19451168
70. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25:2078–2079. https://doi.org/10.1093/bioinformatics/btp352
PMID: 19505943

71. Xie W, Feng Q, Yu H, Huang X, Zhao Q, Xing Y, Yu S, Han B, Zhang Q (2010) Parent-independent genotyping for constructing an ultrahigh-density linkage map based on population sequencing. Proc Natl Acad Sci 107: 10578–10583 https://doi.org/10.1073/pnas.100593110

72. Grimm DG, Roqueiro D, Salomé PA, Kleeberger S, Greshake B, Zhu W, et al. easyGWAS: A Cloud-Based Platform for Comparing the Results of Genome-Wide Association Studies. Plant Cell. 2016; 29:5–19. https://doi.org/10.1105/tpc.16.00551 PMID: 27986896

73. Kang HM, Sul JH, Service SK, Zaitlen NA, Kong SY, Freimer NB, et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010; 42:348–354. https://doi.org/10.1038/ng.548 PMID: 20208533