{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,5]],"date-time":"2026-09-05T11:35:51Z","timestamp":1788608151905,"version":"build-2803163510"},"reference-count":43,"publisher":"Oxford University Press (OUP)","issue":"1","license":[{"start":{"date-parts":[[2026,3,26]],"date-time":"2026-03-26T00:00:00Z","timestamp":1774483200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100007172","name":"Agriculture and Food Research Initiative Competitive","doi-asserted-by":"publisher","award":["2024-67015-42255"],"award-info":[{"award-number":["2024-67015-42255"]}],"id":[{"id":"10.13039\/501100007172","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100005825","name":"US Department of Agriculture\u2019s National Institute of Food and Agriculture","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100005825","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,5,6]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>The landscape of animal and plant breeding is rapidly evolving. Companies now invest millions in collecting extensive data on health, disease, welfare, and fertility, which are key fitness traits relevant for sustainable production. Likewise, genomic-based predictions or GWAS-based searches for causal variants in categorically scored psychiatric disorders and other diseases are key areas in clinical genetics. The categorical nature (yes\/no, scores, etc.) of these traits violates the normality assumptions of standard linear models; threshold or threshold-linear models represent conceptually attractive alternatives. In Bayesian settings, these models can be fitted either via Gibbs sampling or by maximizing the posterior density (maximum a posteriori, MAP). While Gibbs sampling is flexible and can handle multivariate models with several categorical and continuous traits, its computational and memory demands limit its applicability in large-scale applications and genetic and genomic predictions. MAP methods are computationally efficient and suitable for large-scale datasets. However, existing MAP methods are restricted to a single categorical trait and many continuous ones, making them insufficient for the new generation of data. The lack of proper methods for fitting highly dimensional threshold-linear models for quantitative genetics has persisted for over 2 decades. Failing to model categorical traits appropriately can lead to reduced genetic gain and deterioration of the traits over time. This study develops theoretically sound and computationally efficient framework for jointly analyzing multiple categorical and Gaussian traits, specifically multi-trait threshold-linear models, including or not pedigree and genomic information. The proposed MAP multiple-trait threshold-linear model methods are based on Newton-Raphson and Expectation-Maximization schemes. Using a simulated dataset and treating Gibbs sampling as the benchmark, we found that the breeding values estimated by our methods agree nearly perfectly with those obtained from Gibbs sampling, while greatly reducing computational cost and allowing scalability to very large datasets. As expected, the Newton-Raphson iteration outperformed the Expectation-Maximization algorithm in terms of computational time. Our results demonstrate that routine genetic evaluations incorporating multiple categorical traits are now feasible using the presented methodology. Solving this decades-old gap enables accurate, large-scale genetic and genomic predictions of key categorical fitness traits under multi-trait models, unlocking the full potential of genomic selection for complex traits in animal and plant populations and genomic-based predictions in clinical genetics.<\/jats:p>","DOI":"10.1093\/genetics\/iyag086","type":"journal-article","created":{"date-parts":[[2026,3,26]],"date-time":"2026-03-26T13:52:09Z","timestamp":1774533129000},"source":"Crossref","is-referenced-by-count":1,"title":["Methods for joint genetic prediction of multiple ordinal categorical and continuous traits"],"prefix":"10.1093","volume":"233","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-5374-0710","authenticated-orcid":false,"given":"Matias","family":"Bermann","sequence":"first","affiliation":[{"name":"Department of Animal and Dairy Science, University of Georgia , Athens, GA 30602 ,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8893-7620","authenticated-orcid":false,"given":"Andres","family":"Legarra","sequence":"additional","affiliation":[{"name":"Department of Animal and Dairy Science, University of Georgia , Athens, GA 30602 ,","place":["United States"]},{"name":"Council on Dairy Cattle Breeding , Bowie, MD 20716 ,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-0382-1897","authenticated-orcid":false,"given":"Ignacy","family":"Misztal","sequence":"additional","affiliation":[{"name":"Department of Animal and Dairy Science, University of Georgia , Athens, GA 30602 ,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-3140-1002","authenticated-orcid":false,"given":"Daniela","family":"Lourenco","sequence":"additional","affiliation":[{"name":"Department of Animal and Dairy Science, University of Georgia , Athens, GA 30602 ,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,3,26]]},"reference":[{"key":"2026050609302207300_iyag086-B1","doi-asserted-by":"publisher","first-page":"743","DOI":"10.3168\/jds.2009-2730","article-title":"Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score","volume":"93","author":"Aguilar","year":"2010","journal-title":"J Dairy Sci"},{"key":"2026050609302207300_iyag086-B2","doi-asserted-by":"publisher","first-page":"51","DOI":"10.1186\/s12711-019-0469-3","article-title":"Frequentist p-values for large-scale-single step genome-wide association, with an application to birth weight in American Angus cattle","volume":"20","author":"Aguilar","year":"2019","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B3","doi-asserted-by":"publisher","first-page":"4702","DOI":"10.1038\/s41467-023-40330-w","article-title":"Multi-PGS enhances polygenic prediction by combining 937 polygenic scores","volume":"14","author":"Albi\u00f1ana","year":"2023","journal-title":"Nat Commun"},{"key":"2026050609302207300_iyag086-B4","doi-asserted-by":"publisher","first-page":"1518","DOI":"10.3168\/jds.2022-22629","article-title":"Invited review: reliability computation from the animal model era to the single-step genomic model era","volume":"106","author":"Ben Zaabza","year":"2023","journal-title":"J Dairy Sci"},{"key":"2026050609302207300_iyag086-B5","doi-asserted-by":"publisher","first-page":"65","DOI":"10.1137\/141000671","article-title":"Julia: a fresh approach to numerical computing","volume":"59","author":"Bezanson","year":"2017","journal-title":"SIAM Rev"},{"key":"2026050609302207300_iyag086-B6","doi-asserted-by":"publisher","first-page":"880","DOI":"10.1038\/nrg2898","article-title":"Predicting genetic predisposition in humans: the promise of whole-genome markers","volume":"11","author":"de los Campos","year":"2010","journal-title":"Nat Rev Genet"},{"key":"2026050609302207300_iyag086-B7","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1111\/j.2517-6161.1977.tb01600.x","article-title":"Maximum likelihood from incomplete data via the EM algorithm","volume":"39","author":"Dempster","year":"1977","journal-title":"J R Stat Soc Ser B"},{"key":"2026050609302207300_iyag086-B8","doi-asserted-by":"publisher","first-page":"212","DOI":"10.1093\/genetics\/35.2.212","article-title":"Heritability of threshold characters","volume":"35","author":"Dempster","year":"1950","journal-title":"Genetics"},{"key":"2026050609302207300_iyag086-B9","doi-asserted-by":"publisher","first-page":"457","DOI":"10.1093\/biomet\/65.3.457","article-title":"Assessing the accuracy of the maximum likelihood estimator: observed versus expected fisher information","volume":"65","author":"Efron","year":"1978","journal-title":"Biometrika"},{"key":"2026050609302207300_iyag086-B10","volume-title":"Introduction to quantitative genetics.","author":"Falconer","year":"1996","edition":"4th edn"},{"key":"2026050609302207300_iyag086-B11","doi-asserted-by":"publisher","first-page":"401","DOI":"10.1186\/1297-9686-15-3-401","article-title":"Prediction of genetic merit from data on binary and quantitative variates with an application to calving difficulty, birth weight and pelvic opening","volume":"15","author":"Foulley","year":"1983","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B12","doi-asserted-by":"publisher","first-page":"197","DOI":"10.1186\/1297-9686-19-2-197","article-title":"Empirical Bayes estimation of parameters for n polygenic binary traits","volume":"19","author":"Foulley","year":"1987","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B13","doi-asserted-by":"publisher","first-page":"1079","DOI":"10.2527\/jas1982.5451079x","article-title":"Theory and analysis of threshold characters","volume":"54","author":"Gianola","year":"1982","journal-title":"J Anim Sci"},{"key":"2026050609302207300_iyag086-B14","doi-asserted-by":"publisher","first-page":"201","DOI":"10.1186\/1297-9686-15-2-201","article-title":"Sire evaluation for ordered categorical data with a threshold model","volume":"15","author":"Gianola","year":"1983","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B15","doi-asserted-by":"publisher","first-page":"393","DOI":"10.2307\/2531393","article-title":"A mixed-model procedure for analyzing ordered categorical data","volume":"40","author":"Harville","year":"1984","journal-title":"Biometrics"},{"key":"2026050609302207300_iyag086-B16","first-page":"309","article-title":"Estimation of genetic parameters","volume":"21","author":"Henderson","year":"1950","journal-title":"Ann Math Statist"},{"key":"2026050609302207300_iyag086-B17","doi-asserted-by":"publisher","first-page":"skae307","DOI":"10.1093\/jas\/skae307","article-title":"Transforming estimated breeding values from observed to probability scale: how to make categorical data analyses more efficient","volume":"102","author":"Hidalgo","year":"2024","journal-title":"J Anim Sci"},{"key":"2026050609302207300_iyag086-B18","doi-asserted-by":"publisher","first-page":"1609","DOI":"10.2527\/1995.7361609x","article-title":"Multiple-trait genetic evaluation for one polychotomous trait and several continuous traits with missing data and unequal models","volume":"73","author":"Hoeschele","year":"1995","journal-title":"J Anim Sci"},{"key":"2026050609302207300_iyag086-B19","doi-asserted-by":"publisher","first-page":"299","DOI":"10.1186\/1297-9686-18-3-299","article-title":"Genetic evaluation for multiple binary responses","volume":"18","author":"H\u00f6schele","year":"1986","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B20","doi-asserted-by":"publisher","first-page":"183","DOI":"10.1016\/0301-6226(93)90001-X","article-title":"Bivariate analysis for one continuous and one discrete trait with unequal design matrices","volume":"33","author":"Janss","year":"1993","journal-title":"Livest Prod Sci"},{"key":"2026050609302207300_iyag086-B21","doi-asserted-by":"publisher","first-page":"60","DOI":"10.1186\/s12711-020-00580-4","article-title":"Why breed disease-resilient livestock, and how?","volume":"52","author":"Knap","year":"2020","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B22","doi-asserted-by":"publisher","first-page":"159","DOI":"10.1186\/1297-9686-35-2-159","article-title":"Multivariate Bayesian analysis of Gaussian, right censored Gaussian, ordered categorical and binary traits using Gibbs sampling","volume":"35","author":"Korsgaard","year":"2003","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B23","doi-asserted-by":"publisher","first-page":"2872","DOI":"10.1038\/s41588-025-02370-4","article-title":"Liability threshold model-based disease risk prediction based on electronic health record phenotypes","volume":"57","author":"Lee","year":"2025","journal-title":"Nat Genet"},{"key":"2026050609302207300_iyag086-B24","doi-asserted-by":"publisher","first-page":"619","DOI":"10.1111\/j.2517-6161.1996.tb02105.x","article-title":"Hierarchical generalized linear models","volume":"58","author":"Lee","year":"1996","journal-title":"J R Stat Soc Ser B"},{"key":"2026050609302207300_iyag086-B25","doi-asserted-by":"publisher","first-page":"52","DOI":"10.1186\/s12711-023-00816-z","article-title":"Using social interaction models for genetic analysis of skin damage in gilts","volume":"55","author":"Leite","year":"2023","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B26","doi-asserted-by":"publisher","first-page":"59","DOI":"10.1186\/s12711-024-00925-3","article-title":"Marker effect p-values for single-step GWAS with the algorithm for proven and young in large genotyped populations","volume":"56","author":"Leite","year":"2024","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B27","doi-asserted-by":"publisher","first-page":"226","DOI":"10.1111\/j.2517-6161.1982.tb01203.x","article-title":"Finding the observed information matrix when using the EM algorithm","volume":"44","author":"Louis","year":"1982","journal-title":"J R Stat Soc Series B Stat Methodol"},{"key":"2026050609302207300_iyag086-B28","author":"Lourenco","year":"2022"},{"key":"2026050609302207300_iyag086-B29","doi-asserted-by":"publisher","first-page":"1045","DOI":"10.1016\/j.biopsych.2020.01.006","article-title":"Maternal effects as causes of risk for obsessive\u2013compulsive disorder","volume":"87","author":"Mahjani","year":"2020","journal-title":"Biol Psychiatry"},{"key":"2026050609302207300_iyag086-B30","doi-asserted-by":"publisher","first-page":"17","DOI":"10.35566\/jbds\/v1n1\/p2","article-title":"Moments calculation for the doubly truncated multivariate normal density","volume":"1","author":"Manjunath","year":"2021","journal-title":"J Behav Data Sci"},{"key":"2026050609302207300_iyag086-B31","doi-asserted-by":"publisher","first-page":"3062","DOI":"10.3168\/jds.2023-23668","article-title":"A comprehensive meta-analysis of genetic parameters for resilience and productivity indicator traits in Holstein cattle","volume":"107","author":"Maskal","year":"2024","journal-title":"J Dairy Sci"},{"key":"2026050609302207300_iyag086-B32","doi-asserted-by":"publisher","first-page":"267","DOI":"10.1093\/biomet\/80.2.267","article-title":"Maximum likelihood estimation via the ECM algorithm: a general framework","volume":"80","author":"Meng","year":"1993","journal-title":"Biometrika"},{"key":"2026050609302207300_iyag086-B33","doi-asserted-by":"publisher","first-page":"223","DOI":"10.1111\/1467-9868.00382","article-title":"The evaluation of general non-centred orthant probabilities","volume":"65","author":"Miwa","year":"2003","journal-title":"J R Stat Soc Ser B"},{"key":"2026050609302207300_iyag086-B34","doi-asserted-by":"publisher","first-page":"291","DOI":"10.1534\/g3.114.016188","article-title":"Threshold models for genome-enabled prediction of ordinal categorical traits in plant breeding","volume":"5","author":"Montesinos-L\u00f3pez","year":"2014","journal-title":"G3 (Bethesda)"},{"key":"2026050609302207300_iyag086-B35","author":"Quaas","year":"1994"},{"key":"2026050609302207300_iyag086-B36","doi-asserted-by":"publisher","first-page":"229","DOI":"10.1051\/gse:19950303","article-title":"Bayesian inference in threshold models using Gibbs sampling","volume":"27","author":"Sorensen","year":"1995","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B37","doi-asserted-by":"publisher","first-page":"223","DOI":"10.1111\/j.2517-6161.1961.tb00408.x","article-title":"The moment generating function of the truncated multi-normal distribution","volume":"23","author":"Tallis","year":"1961","journal-title":"J R Stat Soc Ser B"},{"key":"2026050609302207300_iyag086-B38","doi-asserted-by":"crossref","first-page":"1428","DOI":"10.3168\/jds.S0022-0302(98)75707-8","article-title":"Generalized linear mixed models in dairy cattle breeding","volume":"81","author":"Tempelman","year":"1998","journal-title":"J Dairy Sci"},{"key":"2026050609302207300_iyag086-B39","doi-asserted-by":"publisher","first-page":"3366","DOI":"10.3168\/jds.S0022-0302(03)73940-X","article-title":"Implementation of a sire-maternal grandsire model for evaluation of calving ease in the United States","volume":"86","author":"Van Tassell","year":"2003","journal-title":"J Dairy Sci"},{"key":"2026050609302207300_iyag086-B40","doi-asserted-by":"publisher","first-page":"335","DOI":"10.1111\/j.1467-9469.2007.00585.x","article-title":"Simple and globally convergent methods for accelerating the convergence of any EM algorithm","volume":"35","author":"Varadhan","year":"2008","journal-title":"Scand J Stat"},{"key":"2026050609302207300_iyag086-B41","doi-asserted-by":"publisher","first-page":"165","DOI":"10.1159\/000446931","article-title":"Concepts and misconceptions about the polygenic additive model applied to disease","volume":"80","author":"Visscher","year":"2015","journal-title":"Hum Hered"},{"key":"2026050609302207300_iyag086-B42","doi-asserted-by":"publisher","first-page":"117","DOI":"10.1186\/1297-9686-29-2-117","article-title":"Bayesian analysis of calving ease scores and birth weights","volume":"29","author":"Wang","year":"1997","journal-title":"Genet Sel Evol"},{"key":"2026050609302207300_iyag086-B43","doi-asserted-by":"publisher","first-page":"506","DOI":"10.1093\/genetics\/19.6.506","article-title":"An analysis of variability in number of digits in an inbred strain of Guinea pigs","volume":"19","author":"Wright","year":"1934","journal-title":"Genetics"}],"container-title":["GENETICS"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/genetics\/advance-article-pdf\/doi\/10.1093\/genetics\/iyag086\/67551540\/iyag086.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/genetics\/article-pdf\/233\/1\/iyag086\/67551540\/iyag086.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/genetics\/article-pdf\/233\/1\/iyag086\/67551540\/iyag086.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,5,6]],"date-time":"2026-05-06T13:30:27Z","timestamp":1778074227000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/genetics\/article\/doi\/10.1093\/genetics\/iyag086\/8541764"}},"subtitle":[],"editor":[{"given":"M","family":"Calus","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2026,3,26]]},"references-count":43,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2026,3,26]]},"published-print":{"date-parts":[[2026,5,6]]}},"URL":"https:\/\/doi.org\/10.1093\/genetics\/iyag086","relation":{},"ISSN":["1943-2631"],"issn-type":[{"value":"1943-2631","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,5]]},"published":{"date-parts":[[2026,3,26]]},"article-number":"iyag086"}}