Genetic Parameters, Prediction of Gains and Intraspeciﬁc Hybrid Selection of Paspalum notatum Flügge for Forage Using REML/BLUP

: Genetic improvement of native forage species is a sustainable alternative for maximizing livestock production. Paspalum notatum Flügge is the most important forage grass in the native grasslands of southern Brazil, with substantial potential available for further genetic improvement. The objective of this study was to quantify a range of genetic parameters and predict yield gains in a population of P. notatum intraspeciﬁc hybrids. Results indicated intraspeciﬁc hybrids of P. notatum had high magnitudes of heritability in the broad and average sense of genotype, plus high selective accuracy and genetic variation for all forage characteristics evaluated. This indicated REML/BLUP can contribute useful information for plant selection in future plant breeding programs. The genetic material studied showed high genetic variability for forage production. Analysis indicated hybrids 336, 332, 437, 132 and male parent ‘30N’ should be included in new crosses to increase the dry matter production of P. notatum . Parents need to be selected from different groups in order to maximize genetic variability and heterosis. In addition, these parents must be included in diallel crosses. The results obtained in this study provide important information for the future breeding of improved P. notatum cultivars for commercialization.


Introduction
Native forage species in pastoral agriculture systems help to provide productive stability and conservation of natural resources and can reduce costs and risks associated with livestock production, culminating in increased sustainability of the system [1].Paspalum notatum Flügge is the main constituent of natural pastures in South America [2] and is a native perennial grass [3].In addition to being cultivated as fodder and cover all over the world [4,5], it has high forage yield [6,7] and is of wide economic importance [8].Furthermore, when intercropped with legumes, the dry matter production of P. notatum was similar to when fertilized with 240 kg N ha −1 [9].
The first studies to identify the mode of reproduction were carried out by Burton [10].Later, Forbes and Burton [11] performed the artificial chromosomal duplication, with colchicine, of a sexual diploid ecotype.A sexual tetraploid was obtained in Argentina through chromosomal duplication directly in the culture medium in plants [12].In Brazil, Weiler et al. [13] were successful in chromosomal duplication from sexual tetraploids.There was the possibility of carrying out crosses to create improved hybrids for commercial Agronomy 2022, 12, 1654 2 of 13 release.Some commercial cultivars were released previously in the southeastern United States and Argentina [4,14].However, in Brazil the only cultivar available is 'Pensacola', which is diploid and has lower forage production compared with tetraploid ecotypes [8,9].It is also less efficient in the use of nitrogen (N) compared with tetraploid ecotypes [10].
Beef cattle grazing systems in Brazil predominantly use forages which have apomixis as a form of reproduction [15].Apomixis refers to a form of asexual reproduction through seeds, which produces progenies genetically identical to the mother plant [16].Among the forages most used in agricultural production in this region, the generas Panicum, Urochloa, Cenchrus and Paspalum are biologically important [17].Paspalum notatum has two cytotypes associated with different modes of reproduction [18], the diploid (2n = 2x = 20), which is sexual and cross-pollinating [10], and the tetraploid (2n = 4x = 40), which is apomictic for apospory [2,10,19].In P. notatum, apomixis is gametophytic of the aposporic type, where embryo sacs originate from cells of the nucellus, and includes processes of apomeiosis, parthenogenesis and pseudogamy [20][21][22][23][24].The most abundant cytotype is the apomictic tetraploid [2], and, consequently, in Rio Grande do Sul/Brazil, P. notatum ecotypes are generally tetraploid and reproduce apomictically [11].Given the mode of reproduction, the exploitation of heterosis becomes a key element in breeding superior hybrids.Therefore, identifying progenitors with superior performance and with dominant allele frequencies for the characteristics of interest is fundamental for exploring hybrid vigor [25] because heterosis results from the accumulation of favorable dominant alleles [26].
The plant breeding process is time consuming, expensive and laborious due to the prolonged period of experimental evaluation and selection of superior genotypes before commercialization [27,28].Therefore, the adoption of more robust and efficient statistical methodologies is essential to guide the process of genetic gain, especially in perennial species [29].Furthermore, in plant breeding, estimates of genetic parameters that produce reliable predictions and information on genetic values are crucial for the success of the program [30].Restricted/residual maximum likelihood/best linear unbiased prediction (REML/BLUP) is considered the standard evaluation method for perennial species.This method offers precision and the possibility to model both fixed and random effects [25,[31][32][33], which is more informative than analysis of variance [34].The variance components of genetic parameters are estimated via REML and genotypic values are predicted via BLUP [31].Effectively, the BLUP presents a favorable characteristic for the shrinkage of the estimators towards the mean and reduces its variance and increases its predictive precision [35].The REML/BLUP procedure has been widely used in annual allogamy breeding crops such as Zea mays (L.) [36][37][38] and autogamy such as Glycine max (L.) [39,40].The procedure has also been applied to selection for perennial crops e.g., Prunus persica (L.) [41], but the approach is less common.Recently, the BLUP procedure was used to determine the specific combining ability of Urochola decumbens [42].In P. notatum, some genetic parameters for forage characteristics have previously been estimated via REML/BLUP [25].
In this study the objective was to estimate a group of genetic parameters and predict yield gains for a population of P. notatum intraspecific hybrids.In addition, two clustering methods were used in order to verify the variability within the tested population.
A randomized complete block design with four replicates was established at the UFRGS (Federal University of Rio Grande do Sul) Experiment Station.A total of 84 genotypes of P. notatum were evaluated.Clones were transplanted into the field with an on-thesquare spacing of 1.0 m within and between rows.All genotypes tested were transplanted simultaneously.Three genotypes, 'C44X' [12], 'Q4188' and 'Q4205' [46], were tetraploid sexual genotypes sourced from the Instituto de Botánica del Nordeste (IBONE-UNNE), Corrientes, Argentina.Seven tetraploid apomitic genotypes ('30N', '36N', '48N', '70N', '83N', '95N' and 'V4') were sourced from collections originally made in South America by the United States Department of Agriculture (USDA).The ploidy level of male parents was determined according to Fachineto et al. [47] and the reproduction mode according to Machado et al. [48].The determination of the ploidy level was performed from the gametic chromosome number, using pollen mother cells analyzed in young inflorescences collected.
In addition, 74 intraspecific hybrids were evaluated (Table 1).These were a result of crosses between female sexual plants (from IBONE-UNNE) and apomitic male genotypes (from the USDA germplasm bank).Hybrids were created according the methodology described by Burton [10] and later adapted by Weiler et al. [49].The genotypes from USDA used in these crosses were selected in a previous evaluation for forage production [47].In the mother plants, before the anthesis, rooted culms bearing inflorescences were collected and placed in a jar with water.These culms were placed in an artificial fog chamber with a high level of air humidity to avoid the anthers from dehiscence.In the morning of the next day, the flowers were emasculated, using sharp pointed tweezers to remove the anthers.Non-flowering spikelets were eliminated from the inflorescence.Fresh pollen from the apomictic parent was collected in paper envelopes and later dusted on emasculated inflorescences of the sexual plants.After pollination, the inflorescences were bagged and labeled.Seeds were collected at least 21 days after pollination.The progeny from each cross was referred as a family, and a number was given to identify each hybrid (Figure 1).
Seeds from the parents and the F1 generation were incubated on Germest paper lined petri dishes, for germination under controlled temperature and day length in a germination chamber: 8 h of light at 30 • C and 16 h of darkness at 20 • C. Afterwards, the seedlings were kept in honeycomb trays until they had five fully expanded leaves.Seedlings with four tillers were then selected and transplanted into pots filled with substrate Carolina Soil™, a commercial substrate, composed of peat, vermiculite, organic residue and limestone.The evaluations were carried out in two growing seasons: (15 March and 26 April 2013) and (12 November 2013; 17 December 2013; 9 January 2014; 2 February 2014).The plants were cut to a 5 cm residual height when they reached an average height of 20 cm to quantify dry matter yield.Samples were sorted into morphological components (leaf blades, stem and inflorescences) then dried in an oven at 60 • C, until constant weight.
to quantify dry matter yield.Samples were sorted into morphological components (leaf blades, stem and inflorescences) then dried in an oven at 60 °C, until constant weight.Plant height (PH, cm), tiller population density (TPD, tillers plant −1 ), accumulated total dry mass (ATDM, kg DM plant −1 ), accumulated leaf dry mass (ALDM, kg DM plant −1 ), accumulated stem dry mass (ASDM, kg DM plant −1 ) and accumulated inflorescence dry mass (AIDM, kg DM plant −1 ) were measured.Non-destructive observations were taken prior to cutting at each date.PH was measured from the soil surface to the average bend of the leaves and then the total number of tillers on each plant were counted to determine TPD.The ATDM, ALDM, ASDM and AIDM yields are the total accumulated DM of each component summed across the entire evaluation period.The leaf:stem ratio (LSR) was calculated from ALDM and ASDM.
Estimates of variance components and genetic parameters were obtained using SEL-EGEN-REML/BLUP software [30] following the REML/BLUP procedure.The genetic statistical model used considered a randomized block design in one location and one year, according to the model below: where y is the data vector, r is the vector of replicate effects (assumed to be fixed), g is the vector of genotypic effects (assumed to be random), e is the error vector (random) and X and Z are the incidence matrices.The mixed model equations for the prediction of r and g are equivalent to: The significance of random effects was obtained through deviance analysis (ANADEV) by the restricted maximum likelihood method, via LRT (likelihood-ratio test).Deviances were obtained following the method of Resende [50], using the model with and without the respective effects, subtracting the deviance obtained in the complete model from the model without the effect and compared with the chi-square (χ²) value with a degree of freedom.Plant height (PH, cm), tiller population density (TPD, tillers plant −1 ), accumulated total dry mass (ATDM, kg DM plant −1 ), accumulated leaf dry mass (ALDM, kg DM plant −1 ), accumulated stem dry mass (ASDM, kg DM plant −1 ) and accumulated inflorescence dry mass (AIDM, kg DM plant −1 ) were measured.Non-destructive observations were taken prior to cutting at each date.PH was measured from the soil surface to the average bend of the leaves and then the total number of tillers on each plant were counted to determine TPD.The ATDM, ALDM, ASDM and AIDM yields are the total accumulated DM of each component summed across the entire evaluation period.The leaf:stem ratio (LSR) was calculated from ALDM and ASDM.
Estimates of variance components and genetic parameters were obtained using SELEGEN-REML/BLUP software [30] following the REML/BLUP procedure.The genetic statistical model used considered a randomized block design in one location and one year, according to the model below: where y is the data vector, r is the vector of replicate effects (assumed to be fixed), g is the vector of genotypic effects (assumed to be random), e is the error vector (random) and X and Z are the incidence matrices.
The mixed model equations for the prediction of r and g are equivalent to: The significance of random effects was obtained through deviance analysis (ANADEV) by the restricted maximum likelihood method, via LRT (likelihood-ratio test).Deviances were obtained following the method of Resende [50], using the model with and without the respective effects, subtracting the deviance obtained in the complete model from the model without the effect and compared with the chi-square (χ 2 ) value with a degree of freedom.The block factor, considered fixed effect, was tested by Snedecor's F test.The genetic parameters via REML for phenotypic (σ 2 p ), genotypic (σ 2 g ) and environmental (σ 2 e ) variance by the EM (expectation-maximization) algorithm are specified by the formulas: where the number of random elements (genotypes) is N g , the sum of the elements of the diagonal matrix (matrix operator) is tr, the total number of data is N, the number of independent linear columns is X and C 22 is determined by the formula: Heritability in the broad sense (h 2 g ) and selective accuracy (r gg ), genetic variation coefficient (CV g ), residual variation coefficient (CV res ) and relative variation coefficient (CV r ) were estimated as follows: ) CV r (%) = CV g CV res (11) Genetic divergence among the hybrids was evaluated by Tocher's clustering method and Unweighted Pair Group Method using Arithmetic averages (UPGMA), by the matrix of genetic distances of Mahalanobis [50].The predicted values were obtained from the variance and covariance matrix of these genetic values as follows: The block factor, considered fixed effect, was tested by Snedecor's F test.The genetic parameters via REML for phenotypic (σ ), genotypic (σ ) and environmental (σ ) variance by the EM (expectation-maximization) algorithm are specified by the formulas: σ   + σ   /  (4) where the number of random elements (genotypes) is  , the sum of the elements of the diagonal matrix (matrix operator) is tr, the total number of data is N, the number of independent linear columns is X and C 22 is determined by the formula: Heritability in the broad sense (h ) and selective accuracy (r ), genetic variation coefficient (CVg), residual variation coefficient (CVres) and relative variation coefficient (CVr) were estimated as follows: Genetic divergence among the hybrids was evaluated by Tocher's clustering method and Unweighted Pair Group Method using Arithmetic averages (UPGMA), by the matrix of genetic distances of Mahalanobis [50].The predicted values were obtained from the variance and covariance matrix of these genetic values as follows:  = ᵟ′Gᵟ, where  is the Mahalanobis distance between genotypes i and i′; G is the matrix of genotypic variance and covariance; ᵟ [d1, d2, … dj], being dj = Yij − Yi′j; and Yij where mean of the i-th genotype in relation to the j-th variable.

Results
In all characteristics studied, the LRT for genotypic effects was significant (p < 0.01) (Table 2).The highest estimates of genetic (σ ) and phenotypic (σ ) variance were for ALDM and TPD.The results indicated there was potential for selection gains within the studied germplasm.For all characteristics there was high genetic variability (CVg), which ranged from a minimum of 27.8 (PH) to a maximum of 78.7 (ASDM).Additionally, high heritability values (h ) were identified within the population studied.High values of these genetic parameters led to high average heritability of genotypes (h ) and high selective accuracy (r ) with all values close to 1.0.This indicated strong genetic control in the studied characteristics in addition to the potential for selection among intraspecific hybrids.The relative variation coefficient (CVr) of all characteristics studied exceeded 1.0.This was The block factor, considered fixed effect, was tested by Snedecor's F test.The genetic parameters via REML for phenotypic (σ ), genotypic (σ ) and environmental (σ ) variance by the EM (expectation-maximization) algorithm are specified by the formulas: σ   + σ   /  (4) where the number of random elements (genotypes) is  , the sum of the elements of the diagonal matrix (matrix operator) is tr, the total number of data is N, the number of independent linear columns is X and C 22 is determined by the formula: Heritability in the broad sense (h ) and selective accuracy (r ), genetic variation coefficient (CVg), residual variation coefficient (CVres) and relative variation coefficient (CVr) were estimated as follows: (%) = 100 σ  (10) Genetic divergence among the hybrids was evaluated by Tocher's clustering method and Unweighted Pair Group Method using Arithmetic averages (UPGMA), by the matrix of genetic distances of Mahalanobis [50].The predicted values were obtained from the variance and covariance matrix of these genetic values as follows:  = ᵟ′Gᵟ, where  is the Mahalanobis distance between genotypes i and i′; G is the matrix of genotypic variance and covariance; ᵟ [d1, d2, … dj], being dj = Yij − Yi′j; and Yij where mean of the i-th genotype in relation to the j-th variable.

Results
In all characteristics studied, the LRT for genotypic effects was significant (p < 0.01) (Table 2).The highest estimates of genetic (σ ) and phenotypic (σ ) variance were for ALDM and TPD.The results indicated there was potential for selection gains within the studied germplasm.For all characteristics there was high genetic variability (CVg), which ranged from a minimum of 27.8 (PH) to a maximum of 78.7 (ASDM).Additionally, high heritability values (h ) were identified within the population studied.High values of these genetic parameters led to high average heritability of genotypes (h ) and high selective accuracy (r ) with all values close to 1.0.This indicated strong genetic control in the studied characteristics in addition to the potential for selection among intraspecific hybrids.The relative variation coefficient (CVr) of all characteristics studied exceeded 1.0.This was , where D 2 ii is the Mahalanobis distance between genotypes i and i ; G is the matrix of genotypic variance and covariance; σ   + σ   /  (4) of random elements (genotypes) is  , the sum of the elements of the atrix operator) is tr, the total number of data is N, the number of indemns is X and C 22 is determined by the formula: the broad sense (h ) and selective accuracy (r ), genetic variation codual variation coefficient (CVres) and relative variation coefficient (CVr) follows: (%) = 100 σ  (10) ence among the hybrids was evaluated by Tocher's clustering method air Group Method using Arithmetic averages (UPGMA), by the matrix s of Mahalanobis [50].The predicted values were obtained from the iance matrix of these genetic values as follows:  = ᵟ′Gᵟ, where  distance between genotypes i and i′; G is the matrix of genotypic varice; ᵟ [d1, d2, … dj], being dj = Yij − Yi′j; and Yij where mean of the i-th n to the j-th variable.
ere performed using the SELEGEN-REML/BLUP genetic-statistical lication [30] and GENES [51] to obtain multivariate analyzes.ristics studied, the LRT for genotypic effects was significant (p < 0.01) est estimates of genetic (σ ) and phenotypic (σ ) variance were for he results indicated there was potential for selection gains within the .For all characteristics there was high genetic variability (CVg), which imum of 27.8 (PH) to a maximum of 78.7 (ASDM).Additionally, high h ) were identified within the population studied.High values of these led to high average heritability of genotypes (h ) and high selective all values close to 1.0.This indicated strong genetic control in the studin addition to the potential for selection among intraspecific hybrids.on coefficient (CVr) of all characteristics studied exceeded 1.0.This was , being d j = Y ij − Y i j ; and Y ij where mean of the i-th genotype in relation to the j-th variable.

Results
In all characteristics studied, the LRT for genotypic effects was significant (p < 0.01) (Table 2).The highest estimates of genetic ( σ2 g ) and phenotypic ( σ2 p ) variance were for ALDM and TPD.The results indicated there was potential for selection gains within the studied germplasm.For all characteristics there was high genetic variability (CV g ), which ranged from a minimum of 27.8 (PH) to a maximum of 78.7 (ASDM).Additionally, high heritability values (h 2 g ) were identified within the population studied.High values of these genetic parameters led to high average heritability of genotypes (h 2 mc ) and high selective accuracy (r gg ) with all values close to 1.0.This indicated strong genetic control in the studied characteristics in addition to the potential for selection among intraspecific hybrids.The relative variation coefficient (CV r ) of all characteristics studied exceeded 1.0.This was a particularly important result because heritability can vary under different environmental conditions, between years and across evaluated characteristics.Therefore, this parameter is of great importance to assist in decision making within the breeding program.Based on the BLUP methodology 15 superior intraspecific hybrids were then selected to quantify forage characteristics.These hybrids represented 17.8% of the total germplasm evaluated (Table 3) and were chosen because their predicted breeding values were higher than the grand mean for all parameters evaluated.For ATDM, ALDM, ASDM, AIDM, LSR, PH, TPD and LS, new averages (X new ) were 458.70 kg DM plant −1 , 270.87 kg DM plant −1 , 144.7855 kg DM plant −1 , 67.6287 kg DM plant −1 , 5.98:1, 27.37 cm and 207.72 tillers per plant −1 , respectively.The genetic value for ATDM and ALDM revealed that intraspecific hybrids 336, 332, 437, 132 and male parent '30N' were superior to other genotypes and hybrids.Genetic gains ranged from 201.1 to 399.51 kg DM plant −1 for ATDM and 223.4 to 347.1 kg DM plant −1 for ALDM, which raised the population mean.
The LSR of hybrids 1026, 525, 225 and female parents 'Q4188 and 'Q4205 showed superior values within the selected population (Table 3).Genetic gains ranged from 2.33 to 4.49.For TPD, genetic gains ranged from 81.6 to 145.8 tillers plant −1 .Hybrids 137, 216, 132, 332 and male parent '48N', were superior to other evaluated material.Based on the most important forage characteristics, the aforementioned intraspecific hybrids were identified as potential parents in new crosses aiming at greater genetic gain for forage production of P. notatum.
The Tocher optimization method identified eight distinct groups of intraspecific hybrids (Table 4).Group I contained the highest concentration of hybrids (64.3%).This demonstrated that these hybrids are more related and have less genetic variation among them.The highest average ATDM was 514.1 kg DM plant −1 (Group V), 501.9 kg DM plant −1 (Group VII), 327.3 kg DM plant −1 (Group VI) and 305.9 kg DM plant −1 (Group II).Group VI had the highest average ALDM (238.1 kg DM plant −1 ) and Group V had the highest average TPD (192.7 tillers plant −1 ).The constituent genotypes of Groups V, VI, VII and II can be included as parents in future crosses.The observed cophenetic correlation coefficient demonstrated an adequate relationship between the distance matrix and the generated dendrogram (Figure 2).The UPGMA clustering method found 10 homogeneous groups which were heterogeneous to the others identified.These results differed to those found by the Tocher method (Table 3).Both methods demonstrated the presence of high genetic variability within the intraspecific hybrids studied (Table 3; Figure 2).The first two groups contained 71% of the hybrids studied, with Group I containing 46% of all hybrids, followed by Group II (Figure 2) with a further 25%.For the ATDM characteristic, groups VI, VII and VIII had the highest average yields of 514.1, 501.9 and 356.7 kg DM plant −1 , respectively.Groups VI (320.1 kg DM plant −1 ), VII (263.4 kg DM plant −1 ) and IV (20.65 kg DM plant −1 ) had superior ALDM compared with the other groups identified.LSR was higher for Groups V (6.3), III (6.1) and IV (3.9).can be included as parents in future crosses.The observed cophenetic correlation coefficient demonstrated an adequate relationship between the distance matrix and the generated dendrogram (Figure 2).The UPGMA clustering method found 10 homogeneous groups which were heterogeneous to the others identified.These results differed to those found by the Tocher method (Table 3).Both methods demonstrated the presence of high genetic variability within the intraspecific hybrids studied (Table 3; Figure 2).The first two groups contained 71% of the hybrids studied, with Group I containing 46% of all hybrids, followed by Group II (Figure 2) with a further 25%.For the ATDM characteristic, groups VI, VII and VIII had the highest average yields of 514.1, 501.9 and 356.7 kg DM plant −1 , respectively.Groups VI (320.1 kg DM plant −1 ), VII (263.4 kg DM plant −1 ) and IV (20.65 kg DM plant −1 ) had superior ALDM compared with the other groups identified.LSR was higher for Groups V (6.3), III (6.1) and IV (3.9).Cophenetic correlation index = 0.76.Group I (blue); Group II (light blue); Group III (orange); Group IV (purple); Group V (light purple); Group VI (green); Group VII (orange); Group VIII (light green); Group IX (red) and Group X (light red).

Discussion
The significance for genetic effects (Table 2) indicated the existence of genetic variability and the possibility of yield gains through targeted selection criteria [52,53].Genetic variance ( σ2 g ) was positive and non-zero values demonstrate greater genetic variation and chances of production gains with selection and relative coefficient of variation values (CV r ) greater than 1.0 indicated a benefit in selection [55].Heritability values (h 2 g ) were >0.50, which is considered high, mainly because it indicates how much of the genotypic variation is of genetic origin [31].Therefore, it is clear that the variance is a particular case of covariance of a variable with itself.Covariance plays a fundamental role in selection as the association and similarity between parents and children is of interest to the breeder and guarantees progress in selection and can be measured by covariance.The narrow sense heritability for forage and seed yield was previously also identified as high [25].However, in this study the correlation values were lower.These values indicated a high probability of success in the selection of hybrids with increased potential forage yields.Selective accuracy (r gg ) is the correlation between the true genotypic value of the genotypes and the estimated/predicted value [56].High values of rgg (close to 1) indicate the experimental design was adequate and the results obtained were reliable (Table 2).The quality of the genotypic assessment should be based on rgg because they provide more accurate inferences of the genetic value of the hybrids evaluated [53].It is important to note that selective accuracy (r gg ) is not associated with the genetic merit of the genotypes studied [57].The set of information generated through the estimation of genetic parameters within the population assists the breeder in selection decisions [30].Estimating genetic parameters via REML and prediction of genotypic values via BLUP offers a more robust and efficient statistical methodology in genetic improvement [26,29,34].Thus, the data obtained here demonstrate high potential to produce genetic gains in the forage production traits studied in future breeding selections.
Breeding programs need to prioritize crossing genotypes with high average values for desirable traits [58].The formation of several groups (Table 4) demonstrated genetic divergence between the genotypes studied.This means there will be numerous alternatives for crosses among the groups identified [59].These authors emphasized that the analysis of genetic divergence simplifies and facilitates the use of germplasm in future crosses.Similar to other studies of genetic divergence, the Tocher method used here also showed a trend towards the formation of a larger group and genotypes isolated in other groups [60][61][62].Multivariate analyses, such as Tocher (Table 4) and UPMG (Figure 2), are powerful tools to assist in the characterization [63] and discrimination of genotypes [64,65] and should be extensively used in genetic diversity studies [66][67][68].Understanding genetic diversity is one of the fundamental steps within genetic improvement programs, with the multivariate analysis approach being the most commonly used [69].Generally, when there is a large amount of data, such as morphological, physiological, biochemical and molecular data, multivariate analyses are used with emphasis on the selection and preservation of genotypes with potential for use within the breeding program [69,70].The use of dendrograms (Figure 2) to graphically describe the clustering method requires the cophenetic correlation coefficient to exceed 0.70 [63].Crosses between elite genotypes with complementary characteristics is desirable to obtain hybrids with enough genetic variability to outperform the parents [71].Crosses between hybrids assigned to the same group are not interesting, mainly because of low variability and the non-exploitation of the existing potential diversity (Table 4; Figure 2).Similarly, the formation of groups with many genotypes can be a limiting factor in the choice of parents for breeding programs because of the proximity between the genotypes within the group [72].Thus, a combination of desirable characteristics and low genetic similarity is required to increase the probability of exploitation of heterosis [8].
This study has shown intraspecific hybrids of P. notatum had high magnitudes of heritability in the broad and average sense of genotype, plus high selective accuracy and genetic variation for all forage characteristics evaluated.This indicated there is potential to select superior hybrids using the REML/BLUP method in future plant breeding forages programs.The genetic material studied included hybrids with high genetic variability for forage production.In the next stage of the breeding program, the selected hybrids (Table 3; Figure 2) can be included in new crosses with female parents that have high genetic value.These parents will be selected from different groups to maximize genetic variability and heterosis.In addition, these must be included in diallel crosses mainly, aiming to select the best parents capable of generating new productive populations, mainly aiming at the general combining ability, which refers to the ability of a parent to produce progenies with a given performance when crossed with a number of other parents, or even in relation to specific combination ability, which refers to the performance of a specific combination.This will allow the results to be confirmed by BLUP and multivariate analysis (Tocher and UPGMA) because these analyses are predictive in nature.Based on diallel crosses, parameters such as heterosis, heterobeltiosis and combining ability (general and specific) can be estimated.The advances obtained with the study of a promising species such as P. notatum will future favor the availability of commercial cultivars for the purpose of use in intensive livestock production systems and in the recovery of degraded areas.In addition, it will also contribute to the conservation of natural grassland areas in South America.

Figure 2 .
Figure 2. Dendrogram of genetic dissimilarity among 74 intraspecific hybrids and eight parents of P. notatum Flügge, obtained by the UPGMA method, based on the average genetic Euclidean distance matrix considering the yield and morphological characteristics evaluated.Cophenetic correlation index = 0.76.Group I (blue); Group II (light blue); Group III (orange); Group IV (purple); Group V (light purple); Group VI (green); Group VII (orange); Group VIII (light green); Group IX (red) and Group X (light red).

Table 3 .
Predicted genotypic effect (g), genotypic value (u + g), genetic gains (gain) and new average (X new ) for characteristics studied in intraspecific hybrids of P. notatum Flügge.

Table 4 .
Composition of groups formed by the Tocher optimization method for evaluated intraspecific hybrids of P. notatum Flügge, based on the average genetic Euclidean distance matrix.

Table 4 .
Composition of groups formed by the Tocher optimization method for evaluated intraspecific hybrids of P. notatum Flügge, based on the average genetic Euclidean distance matrix.