DEVELOPMENT OF THE CROPPOPSTRUCTURE GENOMIC PLATFORM FOR GENETIC DIVERSITY ANALYSIS AND SUPPORT OF BREEDING DECISIONS
Main Article Content
Authors
D. Yerbolatov
Institute of Plant Biology and Biotechnology, Almaty, Kazakhstan
A. Amalova
Institute of Plant Biology and Biotechnology, Almaty, Kazakhstan
G. Yuliya
Institute of Plant Biology and Biotechnology, Almaty, Kazakhstan
Y. Turuspekov
Institute of Plant Biology and Biotechnology, Almaty, Kazakhstan
Abstract
Understanding genetic diversity and population structure is essential for efficiently using crop germplasm collections in breeding, yet a complete analysis requires many separate programs and command-line skills. This study aimed to develop CropPopStructure, a web-based genomic platform that performs a full population genetic analysis of a germplasm collection from a VCF file and a population list, and to validate it on data with a known structure. The platform combines a Nextflow pipeline based on PLINK, VCFtools, and R packages with an R Shiny web interface accessible from a browser. It validates the input data and runs eleven analysis modules: genetic diversity indices, principal component analysis (PCA), discriminant analysis of principal components (DAPC), sNMF ancestry estimation, genomic kinship, a neighbor-joining tree, analysis of molecular variance (AMOVA), pairwise FST, private SNPs, SNP density, and linkage disequilibrium (LD) decay. HapMap files can be converted to VCF with TASSEL within the platform. The platform was tested on a simulated collection of 98 homozygous lines from four populations genotyped with 35,000 SNPs and generated under a hierarchical coalescent model. PCA separated all four populations (PC1 = 6.32%, PC2 = 1.96% of variance), and the BIC criterion indicated two main genetic groups corresponding to the model's two lineages. AMOVA attributed 21.9 % of molecular variance to differences among populations (ΦPT = 0.219; p < 0.001). Pairwise FST values were 0.037–0.039 between sister populations and 0.124–0.125 between lineages, reproducing the simulated hierarchy. CropPopStructure allows breeders and geneticists without bioinformatics training to obtain reproducible results that support forming core collections, selecting genetically distant parents, and controlling population structure in association studies.
Keywords
population structure, genetic diversity, SNP genotyping, germplasm collection, breeding, Nextflow, Shiny, simulation
Article Details
References
Glaszmann J.C., Kilian B., Upadhyaya H.D., Varshney R.K. Accessing genetic diversity for crop improvement // Current Opinion in Plant Biology. 2010. Vol. 13, No. 2. P. 167–173.
McCouch S., Baute G.J., Bradeen J. et al. Agriculture: Feeding the future // Nature. 2013. Vol. 499. P. 23–24.
Varshney R.K., Graner A., Sorrells M.E. Genomics-assisted breeding for crop improvement // Trends in Plant Science. 2005. Vol. 10, No. 12. P. 621–630.
Flint-Garcia S.A., Thornsberry J.M., Buckler E.S. Structure of linkage disequilibrium in plants // Annual Review of Plant Biology. 2003. Vol. 54. P. 357–374.
Yu J., Pressoir G., Briggs W.H. et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness // Nature Genetics. 2006. Vol. 38, No. 2. P. 203–208.
Price A.L., Patterson N.J., Plenge R.M. et al. Principal components analysis corrects for stratification in genome-wide association studies // Nature Genetics. 2006. Vol. 38, No. 8. P. 904–909.
Elshire R.J., Glaubitz J.C., Sun Q. et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species // PLoS ONE. 2011. Vol. 6, No. 5. e19379.
Wang S., Wong D., Forrest K. et al. Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array // Plant Biotechnology Journal. 2014. Vol. 12, No. 6. P. 787–796.
Comadran J., Kilian B., Russell J. et al. Natural variation in a homolog of Antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley // Nature Genetics. 2012. Vol. 44, No. 12. P. 1388–1392.
Pritchard J.K., Stephens M., Donnelly P. Inference of population structure using multilocus genotype data // Genetics. 2000. Vol. 155, No. 2. P. 945–959.
Alexander D.H., Novembre J., Lange K. Fast model-based estimation of ancestry in unrelated individuals // Genome Research. 2009. Vol. 19, No. 9. P. 1655–1664.
Frichot E., Mathieu F., Trouillon T., Bouchard G., François O. Fast and efficient estimation of individual ancestry coefficients // Genetics. 2014. Vol. 196, No. 4. P. 973–983.
Patterson N., Price A.L., Reich D. Population structure and eigenanalysis // PLoS Genetics. 2006. Vol. 2, No. 12. e190.
Jombart T., Devillard S., Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations // BMC Genetics. 2010. Vol. 11. P. 94.
Excoffier L., Smouse P.E., Quattro J.M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data // Genetics. 1992. Vol. 131, No. 2. P. 479–491.
Weir B.S., Cockerham C.C. Estimating F-statistics for the analysis of population structure // Evolution. 1984. Vol. 38, No. 6. P. 1358–1370.
Purcell S., Neale B., Todd-Brown K. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses // American Journal of Human Genetics. 2007. Vol. 81, No. 3. P. 559–575.
Danecek P., Auton A., Abecasis G. et al. The variant call format and VCFtools // Bioinformatics. 2011. Vol. 27, No. 15. P. 2156–2158.
Jombart T. adegenet: a R package for the multivariate analysis of genetic markers // Bioinformatics. 2008. Vol. 24, No. 11. P. 1403–1405.
Bradbury P.J., Zhang Z., Kroon D.E. et al. TASSEL: software for association mapping of complex traits in diverse samples // Bioinformatics. 2007. Vol. 23, No. 19. P. 2633–2635.
Afgan E., Baker D., Batut B. et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update // Nucleic Acids Research. 2018. Vol. 46, No. W1. P. W537–W544.
Di Tommaso P., Chatzou M., Floden E.W. et al. Nextflow enables reproducible computational workflows // Nature Biotechnology. 2017. Vol. 35, No. 4. P. 316–319.
Köster J., Rahmann S. Snakemake — a scalable bioinformatics workflow engine // Bioinformatics. 2012. Vol. 28, No. 19. P. 2520–2522.
Ewels P.A., Peltzer A., Fillinger S. et al. The nf-core framework for community-curated bioinformatics pipelines // Nature Biotechnology. 2020. Vol. 38, No. 3. P. 276–278.
Wilkinson M.D., Dumontier M., Aalbersberg I.J. et al. The FAIR Guiding Principles for scientific data management and stewardship // Scientific Data. 2016. Vol. 3. 160018.
Hoban S., Bertorelle G., Gaggiotti O.E. Computer simulations: tools for population and evolutionary genetics // Nature Reviews Genetics. 2012. Vol. 13, No. 2. P. 110–122.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2024. URL: URL
Chang W., Cheng J., Allaire J.J. et al. shiny: Web Application Framework for R. R package. URL: URL
Chang C.C., Chow C.C., Tellier L.C. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets // GigaScience. 2015. Vol. 4. P. 7.
Grüning B., Dale R., Sjödin A. et al. Bioconda: sustainable and comprehensive software distribution for the life sciences // Nature Methods. 2018. Vol. 15, No. 7. P. 475–476.
Wickham H. ggplot2: Elegant Graphics for Data Analysis. 2nd ed. New York: Springer, 2016. 260 p.
Botstein D., White R.L., Skolnick M., Davis R.W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms // American Journal of Human Genetics. 1980. Vol. 32, No. 3. P. 314–331.
Wright S. The genetical structure of populations // Annals of Eugenics. 1951. Vol. 15, No. 4. P. 323–354.
Schwarz G. Estimating the dimension of a model // The Annals of Statistics. 1978. Vol. 6, No. 2. P. 461–464.
Jombart T., Ahmed I. adegenet 1.3-1: new tools for the analysis of genome-wide SNP data // Bioinformatics. 2011. Vol. 27, No. 21. P. 3070–3071.
Frichot E., François O. LEA: An R package for landscape and ecological association studies // Methods in Ecology and Evolution. 2015. Vol. 6, No. 8. P. 925–929.
Endelman J.B. Ridge regression and other kernels for genomic selection with R package rrBLUP // The Plant Genome. 2011. Vol. 4, No. 3. P. 250–255.
VanRaden P.M. Efficient methods to compute genomic predictions // Journal of Dairy Science. 2008. Vol. 91, No. 11. P. 4414–4423.
Gu Z., Eils R., Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data // Bioinformatics. 2016. Vol. 32, No. 18. P. 2847–2849.
Saitou N., Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees // Molecular Biology and Evolution. 1987. Vol. 4, No. 4. P. 406–425.
Paradis E., Schliep K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R // Bioinformatics. 2019. Vol. 35, No. 3. P. 526–528.
Nei M. Genetic distance between populations // The American Naturalist. 1972. Vol. 106, No. 949. P. 283–292.
Pembleton L.W., Cogan N.O.I., Forster J.W. StAMPP: an R package for calculation of genetic differentiation and structure of mixed-ploidy level populations // Molecular Ecology Resources. 2013. Vol. 13, No. 5. P. 946–952.
Paradis E. pegas: an R package for population genetics with an integrated–modular approach // Bioinformatics. 2010. Vol. 26, No. 3. P. 419–420.
Hill W.G., Robertson A. Linkage disequilibrium in finite populations // Theoretical and Applied Genetics. 1968. Vol. 38, No. 6. P. 226–231.
Kelleher J., Etheridge A.M., McVean G. Efficient coalescent simulation and genealogical analysis for large sample sizes // PLoS Computational Biology. 2016. Vol. 12, No. 5. e1004842.
Baumdicker F., Bisschop G., Goldstein D. et al. Efficient ancestry and mutation simulation with msprime 1.0 // Genetics. 2022. Vol. 220, No. 3. iyab229.
Bhatia G., Patterson N., Sankararaman S., Price A.L. Estimating and interpreting FST: the impact of rare variants // Genome Research. 2013. Vol. 23, No. 9. P. 1514–1521.
Meirmans P.G., Hedrick P.W. Assessing population structure: FST and related measures // Molecular Ecology Resources. 2011. Vol. 11, No. 1. P. 5–18.
Evanno G., Regnaut S., Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study // Molecular Ecology. 2005. Vol. 14, No. 8. P. 2611–2620.
Puechmaille S.J. The program structure does not reliably recover the correct population structure when sampling is uneven: subsampling and new estimators alleviate the problem // Molecular Ecology Resources. 2016. Vol. 16, No. 3. P. 608–627.
Lawson D.J., van Dorp L., Falush D. A tutorial on how not to over-interpret STRUCTURE and ADMIXTURE bar plots // Nature Communications. 2018. Vol. 9. P. 3258.
Kurtzer G.M., Sochat V., Bauer M.W. Singularity: Scientific containers for mobility of compute // PLoS ONE. 2017. Vol. 12, No. 5. e0177459.