Genetic diversity, population structure in a historical panel of Brazilian soybean cultivars.

Soybean [Glycine max (L.) Merrill] is one of the most widely grown legumes in the world, with Brazil being its largest producer and exporter. Breeding programs in Brazil have resulted from multiple cycles of selection and recombination starting from a small number of USA cultivar ancestors in the 1950s and 1960s years. This process has led to the successful adaptation of this crop to tropical conditions, a phenomenon known as tropicalization. Many studies describe a narrow genetic background in Brazilian soybean cultivars. Various factors can affect the genetic diversity in species, especially in cultivated crops, such as the reproduction type, artificial selection, and the number and sources of variability in the breeding programs. In turns, the genetic diversity can affect the linkage disequilibrium blocks (LD) patterns and, consequently, molecular breeding strategies for selection of target loci for agronomic traits. We used high-throughput genotyping with SoySNP50K Illumina SNP markers to assess a collection of 370 Brazilian soybean accessions covering more than 60 years of soybean breeding in Brazil. Our goal was to investigate population structure and genetic diversity in the Brazilian germplasm, detect patterns of LD blocks, and identify regions presenting signals of selective swaps linked with quantitative trait loci (QTLs) of agronomic interest. Population structure analysis revealed two major groups among all genotypes, primarily differentiated by the year of release, separating old and new cultivars (before and after 2000´s years), and by growth habit (stem termination type-SST). The group I comprises about 75% of the panel and includes cultivars release before 2000`s years, including the oldest cultivars released in Brazil, most of which exhibit a determinate growth habit and maturity groups VI and VII. Group II includes only 83 materials, but shows higher levels of diversity than group I, representing most recent introductions in Brazilian germplasm. Further analysis of substructure within Group I, identified seven subgroups with no clear trend for segregation based on maturity group, STT or year of release. Instead, these subgroups were based on the contribution of key donors of disease resistance and adaptability, as soybean cultivation expanded from the South to Central region of Brazil. This finding is consistent with the history of soybean expansion in Brazil. We identified 123 genomic regions under selection among the groups of Brazilian cultivars associated with 440 quantitative trait loci (QTLs), revealing regions fixed across the breeding process associated with yield, disease resistance, water efficiency use, and others.