ABSTRACT. This study evaluated the efficiency of genome-wide selection (GWS) based on regularized quantile regression (RQR) to obtain genomic growth curves based on genomic estimated breeding values (GEBV) of individuals with different probability distributions. The data were simulated and composed of 2,025 individuals from two generations and 435 markers randomly distributed across five chromosomes. The simulated phenotypes presented symmetrical, skewed, positive, and negative distributions. Data were analyzed using RQR considering nine quantiles (0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9) and traditional methods of genomic selection (specifically, RR-BLUP, BLASSO, BayesA, and BayesB). In general, RQR-based estimation of the GEBV was efficient-at least for a quantile model, the results obtained were more accurate than those obtained by the other evaluated methodologies. Specifically, in the symmetrical-distribution scenario, the highest accuracy values were obtained for the parameters with the models RQR0.4, RQR0.3, and RQR0.4. For positive skewness, the models RQR0.2, RQR0.3, and RQR0.1 presented higher accuracy values, whereas for negative skewness, the best model was RQR0.9. Finally, the GEBV vectors obtained by RQR facilitated the construction of genomic growth curves at different levels of interest (quantiles), illustrating the weight-age relationship. ABSTRACT genomewide genome wide GWS (GWS (RQR (GEBV distributions 2025 2 025 2,02 43 chromosomes symmetrical skewed 0.1, 01 0 1 (0.1 02 0.2 03 3 0.3 04 4 0.4 05 5 0.5 06 6 0.6 07 7 0.7 08 8 0.8 0.9 09 9 specifically, specifically (specifically RRBLUP, RRBLUP RR BLUP, BLUP RR-BLUP BLASSO BayesA BayesB. BayesB . BayesB) general RQRbased efficientat efficient methodologies Specifically symmetricaldistribution distribution scenario RQR04 RQR0 RQR0.4 RQR03 RQR0.3 skewness RQR02 RQR0.2 RQR01 RQR0. RQR09 RQR0.9 Finally quantiles, , (quantiles) weightage weight age relationship 202 2,0 0.1 (0. 0. (quantiles 20 2, (0 (

ABSTRACT: Quantile Random Forest (QRF) is a non-parametric methodology that combines the advantages of Random Forest (RF) and Quantile Regression (QR). Specifically, this approach can explore non-linear functions, determining the probability distribution of a response variable and extracting information from different quantiles instead of just predicting the mean. This evaluated the performance of the QRF in the genomic prediction for complex traits (epistasis and dominance). In addition, compare the accuracies obtained with those derived from the G-BLUP. The simulation created an F2 population with 1,000 individuals and genotyped for 4,010 SNP markers. Besides, twelve traits were simulated from a model considering additive and non-additive effects, QTL (Quantitative trait loci) numbers ranging from eight to 120, and heritability of 0.3, 0.5, or 0.8. For training and validation, the 5-fold cross-validation approach was used. For each fold, the accuracies of all the proposed models were calculated: QRF in five different quantiles and three G-BLUP models (additive effect, additive and epistatic effects, additive and dominant effects). Finally, the predictive performance of these methodologies was compared. In all scenarios, the QRF accuracies were equal to or greater than the methodologies evaluated and proved to be an alternative tool to predict genetic values in complex traits. ABSTRACT (QRF nonparametric non parametric RF (RF QR. QR . (QR) Specifically nonlinear linear functions mean epistasis dominance. dominance dominance) addition GBLUP. GBLUP G BLUP. BLUP F 1000 1 000 1,00 4010 4 010 4,01 markers Besides nonadditive effects Quantitative loci 120 03 0 3 0.3 05 5 0.5 08 8 0.8 validation 5fold fold crossvalidation cross used calculated effect effects. effects) Finally compared scenarios (QR 100 00 1,0 401 01 4,0 12 0. 10 1, 40 4,

RESUMO: Quantile Random Forest (QRF) é uma metodologia não paramétrica, que combina as vantagens do Random Forest (RF) e da Regressão Quantílica (QR). Especificamente, essa abordagem pode explorar funções não lineares, determinando a distribuição de probabilidade de uma variável resposta e extraindo informações de diferentes quantis em vez de apenas prever a média. O objetivo deste trabalho foi avaliar o desempenho do QRF em predizer o valor genético genômico para características com arquitetura genética não aditiva (epistasia e dominância). Adicionalmente, as acurácias obtidas foram comparadas com aquelas advindas do G-BLUP. A simulação criou uma população F2 com 1.000 indivíduos genotipados para 4.010 marcadores SNP. Além disso, doze características foram simuladas a partir de um modelo considerando efeitos aditivos e não aditivos, com número de QTL (Quantitative trait loci) variando de oito a 120 e herdabilidade de 0,3, 0,5 ou 0,8. Para treinamento e validação foi usada a abordagem da validação cruzada 5-fold. Para cada um dos folds foram calculadas as acurácias de todos os modelos propostos: QRF em cinco quantis diferentes e três modelos do G-BLUP (com efeito aditivo, aditivo e epistático, aditivo e dominante). Por fim, o desempenho preditivo dessas metodologias foi comparado. Em todos os cenários, as acurácias do QRF foram iguais ou superiores às metodologias avaliadas e mostrou ser uma ferramenta alternativa para predizer valores genéticos em características complexas. RESUMO (QRF paramétrica RF (RF QR. QR . (QR) Especificamente lineares média epistasia dominância. dominância dominância) Adicionalmente GBLUP. GBLUP G BLUP. BLUP F 1000 1 000 1.00 4010 4 010 4.01 SNP disso Quantitative loci 12 03 0 3 0,3 05 5 0, 08 8 0,8 5fold. 5fold fold fold. 5-fold propostos epistático dominante. dominante dominante) fim comparado cenários complexas (QR 100 00 1.0 401 01 4.0 10 1. 40 4.

ABSTRACT Identifying Coffea arabica cultivars that are more efficient in the use of nitrogen is an important strategy and a necessity in the context of environmental and economic impacts attributed to excessive nitrogen fertilization. Although Coffea arabica breeding data have a multi-trait structure, they are often analyzed under a single trait structure. Thus, the objectives of this study were to use a Bayesian multitrait model, to estimate heritability in the broad sense, and to select arabica coffee cultivars with better genetic potential (desirable agronomic traits) in nitrogen-restricted cultivation. The experiment was carried out in a greenhouse with 20 arabica coffee cultivars grown in a nutrient solution with low-nitrogen content (1.5 mM). The experimental design used was in randomized blocks with three replications. Six agromorphological traits of the arabica coffee breeding program and five nutritional efficiency indices were used. The Markov Chain Monte Carlo algorithm was used to estimate genetic parameters and genetic values. The agromorphological traits were considered highly heritable, with a credibility interval (95% probability): H2 = 0.9538 – 5.89E-01. The Bayesian multitrait model presents an adequate strategy for the genetic improvement of arabica coffee grown in low-nitrogen concentrations. Coffee arabica cultivars Icatu Precoce 3282, Icatu Vermelho IAC 4045, Acaiá Cerrado MG 1474, Tupi IAC 1669-33, Catucaí 785/15, Caturra Vermelho and Obatã IAC 1669/20 demonstrated greater potential for cultivation in low-nitrogen concentration. fertilization multi structure Thus sense desirable nitrogenrestricted restricted 2 lownitrogen low 1.5 15 1 5 (1. mM. mM . mM) replications values heritable 95% 95 (95 probability probability) H 09538 0 9538 0.953 5.89E01. 589E01 E 5.89E 01. 89E 01 5.89E-01 concentrations 3282 4045 1474 166933, 166933 1669 33, 33 1669-33 78515 785 785/15 166920 1669/2 concentration 1. (1 9 (9 0953 953 0.95 89E01 5.89E01 589E0 589E 5.89E-0 328 404 147 16693 166 3 1669-3 7851 78 785/1 16692 1669/ ( 095 0.9 89E0 5.89E0 5.89E- 32 40 14 16 1669- 7 785/ 09 0. 4

Abstract This study was conducted to test the significance of adding row and column factors in the frequentist and Bayesian models used in the evaluation of a population of Passiflora edulis, as well as selecting promising genotypes to form the next generation. The following parameters were evaluated: number of fruits, yield, fruit weight, transverse fruit diameter, longitudinal fruit diameter, pulp percentage, skin thickness and total soluble solids. For the Bayesian model, two priors were considered, namely, inverse gamma and a priori distribution with extended parameters. The model with a priori distribution with extended parameters showed lower root mean square error and higher correlation coefficient between observed and predicted values than the inverse gamma model. Furthermore, for a selection intensity of 37%, the mixed and Bayesian models selected practically the same progenies in both experiments. The use of the 5-fold cross-validation technique indicated that both tested models were efficient. edulis generation evaluated fruits yield weight diameter percentage solids considered namely Furthermore 37 37% experiments 5fold fold 5 crossvalidation cross validation efficient 3

ABSTRACT Among the multi-trait models selected to study several traits and environments jointly, the Bayesian framework has been a preferred tool when constructing a more complex and biologically realistic model. In most cases, non-informative prior distributions are adopted in studies using the Bayesian approach. However, the Bayesian approach presents more accurate estimates when informative prior distributions are used. The present study was developed to evaluate the efficiency and applicability of multi-trait multi-environment (MTME) models within a Bayesian framework utilizing a strategy for eliciting informative prior distribution using previous data on rice. The study involved data pertaining to rice (Oryza sativa L.) genotypes in three environments and five crop seasons (2010/2011 until 2014/2015) for the following traits: grain yield (GY), flowering in days (FLOR) and plant height (PH). Variance components, genetic and non-genetic parameters were estimated using the Bayesian method. In general, the informative prior distribution in Bayesian MTME models provided higher estimates of individual narrow-sense heritability and variance components, as well as minor lengths for the highest probability density interval (HPD), compared to their respective non-informative prior distribution analyses. More informative prior distributions make it possible to detect genetic correlations between traits, which cannot be achieved with non-informative prior distributions. Therefore, this mechanism presented to update knowledge for an elicitation of an informative prior distribution can be efficiently applied in rice breeding programs.

RESUMO: O processo de desenvolvimento de uma nova cultivar de trigo requer tempo entre a obtenção da população base e a seleção da linhagem mais promissora. Estimar parâmetros genéticos com mais precisão nas primeiras gerações com vistas a antecipar a seleção significa avanços importantes para os programas de melhoramento de trigo. Assim, o presente estudo estima os parâmetros genéticos de populações F2 de trigo tropical e o ganho genético da seleção via abordagem Bayesiana. Para tanto, os autores avaliaram a produtividade de grãos por parcela de 34 populações F2 de trigo tropical. A abordagem Bayesiana proporcionou um ajuste adequado ao modelo, estimando parâmetros genéticos dentro do espaço paramétrico. A herdabilidade (h2) foi de 0,51. Dentre as selecionadas, 11 populações F2 obtiveram desempenho superior às cultivares controle, com ganho genético de seleção de 7,80%. As seguintes populações foram as mais promissoras: Tbio Sossego/CD 1303, CD 1303/Tbio Ponteiro, BRS 254/CD 1303, Tbio Duque/Tbio Aton e Tbio Aton/CD 1303. A inferência Bayesiana pode ser usada para melhorar significativamente programas de melhoramento de trigo tropical.

ABSTRACT: The development process of a new wheat cultivar requires time between obtaining the base population and selecting the most promising line. Estimating genetic parameters more accurately in early generations with a view to anticipating selection means important advances for wheat breeding programs. Thus, the present study estimated the genetic parameters of F2 populations of tropical wheat and the genetic gain from selection via the Bayesian approach. To this end, the authors assessed the grain yield per plot of 34 F2 populations of tropical wheat. The Bayesian approach provided an adequate fit to the model, estimating genetic parameters within the parametric space. Heritability (h2) was 0.51. Among those selected, 11 F2 populations performed better than the control cultivars, with genetic gain of 7.80%. The following populations were the most promising: TbioSossego/CD 1303, CD 1303/TbioPonteiro, BRS 254/CD 1303, Tbio Duque/Tbio Aton, and Tbio Aton/CD 1303. Bayesian inference can be used to significantly improve tropical wheat breeding programs.

Resumo O objetivo deste trabalho foi avaliar um procedimento para elicitação de distribuição a priori informativa, comparada à distribuição a priori não informativa, em tamanho amostral reduzido, com uso de 14 caracteres de três genótipos de linhaça (Linum usitatissimum), em sete épocas de semeadura. Os valores dos hiperparâmetros regulam a informatividade da distribuição a priori; portanto, para cada época, foram calculados os hiperparâmetros a serem utilizados na próxima época. As duas distribuições a priori, não informativa e informativa, foram comparadas pelo comprimento dos intervalos de credibilidade e pela variância da distribuição a posteriori. Em geral, quando a distribuição a priori informativa é adotada, os parâmetros genéticos apresentam menor comprimento do intervalo de credibilidade e estimativas mais precisas. O mecanismo de elicitação a priori informativa com uso de informações prévias de programas de melhoramento é eficiente para estimativa de parâmetros genéticos, incluindo herdabilidade e variância genética, mesmo quando o tamanho da amostra é pequeno. Na avaliação genética, o uso da distribuição a priori informativa é melhor do que o da distribuição não informativa para tamanho amostral pequeno. Em geral, os resultados das distribuições a priori informativas indicam que os valores genéticos da primeira época de semeadura são maiores para os seguintes caracteres: comprimento do ciclo, altura da planta, e número de cápsulas não granuladas e de ramos produtivos.

Abstract The objective of this work was to evaluate a procedure for the elicitation of informative prior distribution, compared with non-informative prior distribution, in a small sample size, using 14 traits of three linseed (Linum usitatissimum) genotypes in seven sowing seasons. The values of the hyperparameters regulate the informativeness of the prior distribution; therefore, for each season, the hyperparameters to be used in the next season were calculated. The two prior distributions, non-informative and informative, were compared by the length of the credible interval and variance of the posterior distribution. In general, when the informative prior distribution is adopted, the genetic parameters present a shorter length of the credible interval and more precise estimates. The mechanism for informative prior elicitation using previous information from breeding programs is efficient for the estimation of genetic parameters, including heritability and genetic variance, even when the sample size is small. In genetic evaluation, the use of informative prior distribution is better than that of non-informative distribution for a small sample size. In general, the results of the informative prior distributions are indicative that the genetic values of the first sowing season are greater for the following traits: cycle length, plant height, and number of non-grained capsules and of productive branches.

Abstract Bayesian approaches applied in association studies select regions of single-nucleotide polymorphisms, indicating genes with important effects. The Bayesian methods differ in terms of the distribution assumed for the marker effects. Here, we used the window posterior probability of association to detect potential regions. The present study evaluated the efficiency of these methods in identifying regions located close to genes. Data were simulated in six scenarios. Considering the lack of dominance, BayesA was more efficient in the scenario with three QTLs. For scenarios with 10 or 100 QTLs, BayesCπ and BayesDπ were more efficient according to the false positive rate and detection power. Considering the presence of dominance, all methods were similar in the scenario with three QTLs, except in terms of accuracy. BayesDπ was superior in the scenario with 10 QTLs, while BRR was more efficient in the scenario with 100 QTLs.

ABSTRACT: The development of efficient methods for genome–wide association studies (GWAS) between quantitative trait loci (QTL) and genetic values is extremely important to animal and plant breeding programs. Bayesian approaches that aim to select regions of single nucleotide polymorphisms (SNPs) proved to be efficient, indicating genes with important effects. Among the selection criteria for SNPs or regions, selection criterion by percentage of variance can be explained by genomic regions (%var), selection of tag SNPs, and selection based on the window posterior probability of association (WPPA). To also detect potentially associated regions, we proposed measuring posterior probability of the interval PPint), which aims to select regions based on the markers of greatest effects. Therefore, the objective of this work was to evaluate these approaches, in terms of efficiency in selecting and identifying markers or regions located within or close to genes associated with traits. This study also aimed to compare these methodologies with single–marker analyses. To accomplish this, simulated data were used in six scenarios, with SNPs allocated in non–overlapping genomic regions. Considering traits with oligogenic inheritance, WPPA criterion followed by %var and PPint criteria were shown to be superior, presenting higher values of detection power, capturing higher percentages of genetic variance and larger areas. For traits with polygenic inheritance, PPint and WPPA criteria were considered superior. Single–marker analyses identified SNPs associated only in oligogenic inheritance scenarios and was lower than the other criteria.

ABSTRACT Methods for genetic improvement of semi–perennial species, such as passion fruit, often involve large areas, unbalanced data, and lack of observations. Some strategies can be applied to solve these problems. In this work, different models and approaches were tested to improve the precision of estimates of genetic evaluation models for several characteristics of the passion fruit. A randomized block design (RBD) model was compared to a posteriori correction, adding two factors to the model (post–hoc blocking Row–Col). These models were also combined with the frequentist and Bayesian approaches to identify which combination yields the most accurate results. These approaches are part of a strategic plan in a perennial plant breeding program to select promising genitors of passion to compose the next selection cycle. For Bayesian, we tested two priors, defining different values for the distribution parameters of effect variances of the model. We also performed a cross–validation test to choose a priori values and compare the frequentist and Bayesian approaches using the root mean square error (RMSE) and the correlation between the predicted and observed values, called Predictive capacity of the model (PC). The model with the post–hoc blocking Row–Col design captured the spatial variability for productivity and number of fruits, directly affecting the experimental precision. Both approaches applied to the models showed a similar performance, with predictive capacity and selective efficiency leading to the selection of the same individuals.

ABSTRACT: The Fisher’s infinitesimal model is traditionally used in quantitative genetics and genomic selection, and it attributes most genetic variance to additive variance. Recently, the dominance maximization model was proposed and it prioritizes the dominance variance based on alternative parameterizations. In this model, the additive effects at the locus level are introduced into the model after the dominance variance is maximized. In this study, the new parameterizations of additive and dominance effects on quantitative genetics and genomic selection were evaluated and compared with the parameterizations traditionally applied using the genomic best linear unbiased prediction method. As the parametric relative magnitude of the additive and dominance effects vary with allelic frequencies of populations, we considered different minor allele frequencies to compare the relative magnitudes. We also proposed and evaluated two indices that combine the additive and dominance variances estimated by both models. The dominance maximization model, along with the two indices, offers alternatives to improve the estimates of additive and dominance variances and their respective proportions and can be successfully used in genetic evaluation.

ABSTRACT: The principal component regression (PCR) and the independent component regression (ICR) are dimensionality reduction methods and extremely important in genomic prediction. These methods require the choice of the number of components to be inserted into the model. For PCR, there are formal criteria; however, for ICR, the adopted criterion chooses the number of independent components (ICs) associated to greater accuracy and requires high computational time. In this study, seven criteria based on the number of principal components (PCs) and methods of variable selection to guide this choice in ICR are proposed and evaluated in simulated and real data. For both datasets, the most efficient criterion and that drastically reduced computational time determined that the number of ICs should be equal to the number of PCs to reach a higher accuracy value. In addition, the criteria did not recover the simulated heritability and generated biased genomic values.

ABSTRACT Genomic selection (GS) emphasizes the simultaneous prediction of the genetic effects of thousands of scattered markers over the genome. Several statistical methodologies have been used in GS for the prediction of genetic merit. In general, such methodologies require certain assumptions about the data, such as the normality of the distribution of phenotypic values. To circumvent the non-normality of phenotypic values, the literature suggests the use of Bayesian Generalized Linear Regression (GBLASSO). Another alternative is the models based on machine learning, represented by methodologies such as Artificial Neural Networks (ANN), Decision Trees (DT) and related possible refinements such as Bagging, Random Forest and Boosting. This study aimed to use DT and its refinements for predicting resistance to orange rust in Arabica coffee. Additionally, DT and its refinements were used to identify the importance of markers related to the characteristic of interest. The results were compared with those from GBLASSO and ANN. Data on coffee rust resistance of 245 Arabica coffee plants genotyped for 137 markers were used. The DT refinements presented equal or inferior values of Apparent Error Rate compared to those obtained by DT, GBLASSO, and ANN. Moreover, DT refinements were able to identify important markers for the characteristic of interest. Out of 14 of the most important markers analyzed in each methodology, 9.3 markers on average were in regions of quantitative trait loci (QTLs) related to resistance to disease listed in the literature.

ABSTRACT Genetic diversity studies are performed based on information on a set of traits measured in a group of genotypes, considering one or more environments. The pattern recognition methods allow classifying genotypes from a set of important agronomic information. Thus, this study aimed to present and compare pattern recognition methods to inquire about the similarity of environments and genotypes in flood-irrigated rice for the recommendation of cultivars. The experim>ents were performed in the municipalities of Leopoldina, Lambari, and Janaúba, state of Minas Gerais, Brazil. To evaluate the pattern of similarity, 25 rice genotypes in three environments belonging to the flood-irrigated rice breeding program were used. Among these genotypes, five cultivars were used as an experimental control for the grain yield, the height of the plant, flowering, panicle length, grains filled by panicles, percentages of grains filled by panicles, in the 2012/2013 agricultural year. The methods us>ed were mixtures of multivariate normal distributions and density-based clustering algorithm. It was observed, therefore, that the genotypes are distributed in three distinct groups, in which there are intragroup homogeneity and intergroup heterogeneity for the agronomic traits of the flooded rice culture. The methods >used to assess the dissimilarity of environments using pattern recognition methods were efficient in classifying flooded rice irrigated environments.

ABSTRACT: Genome-wide selection (GWS) is currently a technique of great importance in plant breeding, since it improves efficiency of genetic evaluations by increasing genetic gains. The process is based on genomic estimated breeding values (GEBVs) obtained through phenotypic and dense marker genomic information. In this context, GEBVs of N individuals are calculated through appropriate models, which estimate the effect of each marker on phenotypes, allowing the early identification of genetically superior individuals. However, GWS leads to statistical challenges, due to high dimensionality and multicollinearity problems. These challenges require the use of statistical methods to approach the regularization of the estimation process. Therefore, we aimed to propose a method denominated as triple categorical regression (TCR) and compare it with the genomic best linear unbiased predictor (G-BLUP) and Bayesian least absolute shrinkage and selection operator (BLASSO) methods that have been widely applied to GWS. The methods were evaluated in simulated populations considering four different scenarios. Additionally, a modification of the G-BLUP method was proposed based on the TCR-estimated (TCR/G-BLUP) results. All methods were applied to real data of cassava (Manihot esculenta) with to increase efficiency of a current breeding program. The methods were compared through independent validation and efficiency measures, such as prediction accuracy, bias, and recovered genomic heritability. The TCR method was suitable to estimate variance components and heritability, and the TCR/G-BLUP method provided efficient GEBV predictions. Thus, the proposed methods provide new insights for GWS.