A selected set of 136 rapeseed accessions were genotyped using the Brassica 60K Illumina Infinium@ SNP array (Illumina) as described previously (Liu et al., 2016 (link)). After quality control (minor allele frequency>0.05 and missing data<20%) (Liu and Muse, 2005 (link)), a total of 21426 SNPs were selected for subsequent GWAS analysis. SNP density map was drawn using the CMplot package (Yin et al., 2021 (link)) in R3.3.3 (Team RC, 2013 ) to gauge the marker coverage across the B. napus genome. Principal component analysis was performed using Tassel v5.0 software (Bradbury et al., 2007 (link)). Only the first two principal components were plotted by the ggplot2 package (Wickham, 2011 (link)) in R3.3.3 (Team RC, 2013 ). MEGA-X software (Kumar et al., 2018 (link)) was used to draw the phylogenetic tree and show the genetic relationship among GWAS accessions. Both principal component analysis and phylogenetic tree are used to reveal the distribution of population structure. LD (linkage disequilibrium) attenuation was estimated to demonstrate the degree of genetic linkage, and it is drawn by PopLDdecay software (max decay distant = 2500, break = 6000bp, bin1 = 1000, bin2 = 6000) as described previously (Zhang et al., 2019 (link)).
Free full text: Click here