Simultaneously accounting for population structure, genotype by environment interaction, and spatial variation in marker-trait associations in sugarcane

Few association mapping studies have simultaneously accounted for population structure, genotype by environment interaction (GEI), and spatial variation. In this sugarcane association mapping study we tested models accounting for these factors and identified the impact that each model component had on the list of markers declared as being significantly associated with traits. About 480 genotypes were evaluated for cane yield and sugar content at three sites and scored with DArT markers. A mixed model was applied in analysis of the data to simultaneously account for the impacts of population structure, GEI, and spatial variation within a trial. Two forms of the DArT marker data were used in the analysis: the standard discrete data (0, 1) and a continuous DArT score, which is related to the marker dosage. A large number of markers were significantly associated with cane yield and sugar content. However, failure to account for population structure, GEI, and (or) spatial variation produced both type I and type II errors, which on the one hand substantially inflated the number of significant markers identified (especially true for failing to account for GEI) and on the other hand resulted in failure to detect markers that could be associated with cane yield or sugar content (especially when failing to account for population structure). We concluded that association mapping based on trials from one site or analysis that failed to account for GEI would produce many trial-specific associated markers that would have low value in breeding programs.