Next Article in Journal
Deciphering Reserve Mobilization, Antioxidant Potential, and Expression Analysis of Starch Synthesis in Sorghum Seedlings under Salt Stress
Next Article in Special Issue
Reference-Guided De Novo Genome Assembly to Dissect a QTL Region for Submergence Tolerance Derived from Ciherang-Sub1
Previous Article in Journal
Radical Scavenging Actions and Immunomodulatory Activity of Aronia melanocarpa Propylene Glycol Extracts
Previous Article in Special Issue
Opportunities for Improving Waterlogging Tolerance in Cereal Crops—Physiological Traits and Genetic Mechanisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Formation of Potential Heterotic Groups of Oat Using Variation at Microsatellite Loci

by
Michaela Havrlentová
1,2,*,
Katarína Ondreičková
1,
Peter Hozlár
1,
Veronika Gregusová
2,
Daniel Mihálik
1,2 and
Ján Kraic
1,2
1
National Agricultural and Food Center, Research Institute of Plant Production, Bratislavská Cesta 122, 92168 Piešťany, Slovakia
2
Department of Biotechnology, Faculty of Natural Sciences, University of Ss. Cyril and Methodius, Námestie Jozefa Herdu 2, 91701 Trnava, Slovakia
*
Author to whom correspondence should be addressed.
Plants 2021, 10(11), 2462; https://doi.org/10.3390/plants10112462
Submission received: 30 September 2021 / Revised: 11 November 2021 / Accepted: 12 November 2021 / Published: 15 November 2021
(This article belongs to the Special Issue Cereals Genetic Resources and Improvement)

Abstract

:
An evaluation of polymorphism at the microsatellite loci was applied in distinguishing 85 oat (Avena sativa L.) genotypes selected from the collection of genetic resources. The set of genotypes included oats with white, yellow, and brown seeds as well as a subgroup of naked oat (Avena sativa var. nuda Koern). Variation at these loci was used to form potential heterotic groups potentially used in the oat breeding program. Seven from 20 analyzed microsatellite loci revealed polymorphism. Altogether, 35 microsatellite alleles were detected (2–10 per locus). Polymorphic patterns completely differentiated all genotypes within the subgroups of white, brown, and naked oats, respectively. Only within the greatest subgroup of yellow genotypes, four pairs of genotypes remained unseparated. Genetic differentiation between the oat subgroups allowed the formation of seven potential heterotic groups using the STRUCTURE analysis. The overall value of the fixation index (Fst) suggested a high genetic differentiation between the subgroups and validated a heterotic grouping. This approach can be implemented as a simple predictor of heterosis in parental crosses prior to extensive field testing or development and implementation of more accurate genomic selection.

1. Introduction

Cultivated oat (Avena sativa L.) is an economically important crop, ranking sixth in world cereal production after wheat, rice, maize, barley, and sorghum. Oat is used as green fodder, straw, hay, or silage with good balanced feed components for livestock. Oat grains are also part of humans’ diet, either directly as raw food (flakes, milk) or as raw material and ingredients in food production. They possess unique and important nutritive characteristics, particularly high contents of lipids, proteins, and micronutrients [1]. Compared with other cereals, oat grains are rich in antioxidants (e.g., α-tocotrienol, α-tocopherol, and avenanthramides) total dietary fiber, and the water-soluble dietary fiber, the β-D-glucans [2,3,4]. Generally, oat grains have multifunctional use in human nutrition, animal feeding, and in the production of health care and cosmetic products [5].
Plant breeders frequently use genetically similar parents in crosses to generate progenies to create new cultivars. Such an approach is also common among oat breeders. Elite but genetically very similar parents have the potential to generate advanced offspring without an undesirably changed genetic background. However, if breeders need to substantially extend genetic diversity and introduce new genes, they must select and cross parents genetically as different as possible. Genetically divergent parents can generate advanced hybrids with higher heterozygosity and manifested heterotic effects. Genetic diversity as well as unrelatedness is usually determined by pedigree analysis, morphological, physiological, agronomical traits, and biometric analysis. These approaches are limited mainly by low variation in these parameters and have been enriched gradually after the introduction of molecular, especially DNA markers. Molecular markers can effectively contribute to the formation of heterotic groups. Subsequently, parental genotypes can be selected from distinct groups with the expectation of greater heterosis effect in offspring. However, the reliability of molecular markers to confirm the relationship between genetic diversity and heterosis in offspring is neither completely clear nor universal. It depends on the type of molecular markers used, the plant species, the observed phenotypic traits, and other factors. Melchinger [6] accumulated results for numerous crops and concluded that DNA markers are well-suited for assessment of the genetic diversity, grouping of genotypes, and selecting parents to establish the base populations. Nevertheless, he [6] expected that the heterotic response between these groups cannot be simply predicted from genetic distances detected by DNA markers, but it needs to be evaluated and confirmed in real field trials. Later, the relationship between the heterosis effect and the genetic distance between parental components determined by several types of molecular markers was evaluated [7]. Among three types of molecular markers (SSR, AFLP, RAPD), the parental components for heterosis crosses were the most effectively selected by microsatellite markers (SSR) and confirmed by calculated Jaccard, Kluczyński, Nei, and Rogers coefficients, respectively [7]. Other marker types were less useful (AFLP) or declared no clear relationship between genetic distance and the heterotic effect (RAPD) [7]. The number of such studies performed in oats is very limited. However, the crossing between genotypes belonging to genetically different groups may result in high heterosis usable in breeding improvement of oat [8,9]. A genetic diversity analysis in oats themselves was performed using various molecular markers [10,11,12,13,14,15]. Relatively high polymorphism and ease of use favor SSRs located in non-coding (SSRs) as well as in coding (EST-SSRs) DNA sequences [16,17,18]. Oat breeders are also interested in the effect of heterosis obtained by the crossing of thoroughly selected parent components. An analysis of genetic dissimilarity between selected parents can be a tool to achieve this. Therefore, the aim of this study was to determine the extent of genetic diversity in a set of oats (Avena sativa L.) genotypes from a maintained collection of genetic resources using the microsatellite polymorphism for selection of genetically different genotypes and formation of potential heterotic groups for the oat breeding program.

2. Results

2.1. Informativeness of SSR Markers

Twenty pairs of primers for the polymorphism analysis at the microsatellite loci were tested. Only seven (AM1, AM14, AM22, AM83, AM87, AM102, AM115) revealed polymorphism within the set of 85 analyzed oat genotypes. Altogether, 35 alleles were detected, 2–10 alleles per locus, and on average 5 alleles per locus. The most polymorphic one was the locus AM1 (10 alleles) where the highest number of different microsatellite patterns (31) were found. A heterozygote status was detected at several loci.
The genetic variation at all polymorphic microsatellite loci was declared by several parameters (Table 1). The heterozygosity indices (i.e., expected heterozygosity) showed the probability that an individual genotype would be heterozygous at a given locus. Its values ranged from 0.327 to 0.5. The polymorphic information content (PIC) values showed that all seven analyzed markers (loci) had approximately the same predictive value. Although the PIC values were not high (below 0.5), together with the Heterozygosity index (H), the values suggested that primers designed for the analysis at these microsatellite loci have discriminatory competence.
The marker index (MI) is the product of the effective multiplex ratio (E) and Havp and should reveal the distinguishing power of markers or techniques. It predicts the relative utility of various marker systems [19]. The MI proved to be useful in comparing the efficacy of different marker systems in some studies [20,21]. Prevost and Wilkinson [22] argued that there is little or no correlation between MI and the ability of primers to distinguish genotypes, and the practical value of MI for this purpose is limited. Only one type of marker (microsatellite) was used in our work, so the MI had no informative value in this case. Its values were low due to very low Havp values.
The parameter resolving power (R) compared the diagnostic effectiveness of used primers within the analyzed set of oats. The most effective was the marker locus AM1 with the highest R value (4.1176). According to the equation of Prevost and Wilkinson [22] 0.15x + 1.78 = R (where x is the number of genotypes identified), this marker alone can identify 15 oat genotypes. Other markers had lower R values (Table 3). Even this marker is useful for discrimination of this set of oats only in combination with many others. Theoretically, at least six markers with the R value higher than 4.2 would be required for complete discrimination of all 85 evaluated oat genotypes.
The Discriminating power (D) described by Tessier et al. [23] can be a good estimator of the efficiency of individual markers or marker combinations to describe the probability that two randomly chosen genotypes have different patterns. However, the efficiency of a given marker does not depend only on the number of patterns it generates. The higher the D value (closest to 1) the lower the probability of confusion between microsatellite profiles within the analyzed oat genotypes. Both, the mean (0.8538) and individual (0.7248–0.9578) D values were high to very high in our experiment (Table 1).
Generally, the polymorphism revealed at the genomic microsatellite loci in the oat genome tends to be high. It is characterized by heterozygosity and a high number of alleles per locus, which usually results in high PIC values, diversity indices, and a high ability to differentiate genotypes [16,17,24,25]. Many of the developed EST-SSR markers have similar robustness as microsatellites (SSR) markers [26,27,28]. The effectiveness of oat genotypes distinguishing relates to the extent of the genetic diversity within the analyzed set of genotypes. This, in turn, usually depends on the geographical origin, pedigrees, and type of germplasm (cultivars, breeding lines, landraces, wild relatives). The second factor is the appropriate choice of microsatellite marker, either genomic (SSR) or expression sequence tags (EST-SSRs). The set of 85 oats used in this study contained considerable genetic variation and the parameters of all seven used polymorphic microsatellite markers had high informativeness. This allowed the subsequent analysis of genetic diversity and population structure within a set of oat genotypes.

2.2. Genetic Diversity

The relationships between subgroups of oat genotypes differing in seed parameters (glume color and hull/naked) faithfully demonstrated the principal component analysis (PCA). This type of statistical analysis is suitable for the discrimination of genotypes into groups and subgroups based on selected parameters. The classification of oat genotypes by PCA into subgroups according to seed parameters (white, yellow, brown, and naked) was created in this case. Four subgroups, white, yellow, brown, and naked oats, contained a different number of genotypes. A relatively high variation in analyzed microsatellite loci was revealed within each of these four subgroups. Nevertheless, it should be noted that separate subgroups overlapped and no boundaries and no separation between them were evident (Figure 1). Therefore, the relationships between the aggregating of genotypes based on variations at the microsatellite loci and the observed seed parameters could not be identified.
Although the PCA did not show the formation of separate subgroups, it revealed the differentiation ability of microsatellite markers. This analysis apparently demonstrated the suitability for differentiating genotypes of a relatively large group as a whole and individuals within four subgroups (white, yellow, brown, and nude). All genotypes within the subgroups of white, brown, and naked oats were completely differentiated from each other. Only four pairs remained unseparated within the subgroup of yellow genotypes: Hecht-Hron (PS-100), Auron-Senator, Ardo-Cyril, and Expander-Expo, respectively. Pairwise comparisons based on the PCA analysis revealed that only the brown oats were not statistically different from other subgroups of oats (Table 2)
A high differentiation competence of the microsatellite markers is better demonstrated though a cluster analysis that also indicates the possible grouping of genotypes into smaller groups/subgroups (Figure 2). It revealed the tendency for the grouping of naked and brown oats (in the middle of the dendrogram), but partially also white genotypes (in the upper half of the dendrogram). Yellow oats were scattered throughout the dendrogram.

2.3. Population Structure

The estimation of the population structure was performed by the STRUCTURE software using the program Structure Harvester [29]. The ∆K determining the best fitted value was at K = 2 (Figure 3). The genotypes were differentiated at K = 2 only into two subgroups that correspond to the two main clusters formed also by the hierarchical cluster analysis (Figure 2). There were nineteen genotypes in one subgroup (cluster at the bottom of the dendrogram, Figure 2). All the other genotypes were in the second large subgroup. However, there was no rational reason for the separation of genotypes into only two subgroups. They have different geographical origins, the period of origin, pedigrees, agronomic parameters, chemical composition of seeds, and other traits. In addition, the molecular markers used were genomic microsatellites, so it was neither reasonable to separate the genotypes into two subgroups, nor to look for associations with the observed seed parameters. Therefore, K = 2 was not used in the STRUCTURE analysis. The aim of the study was to create more subsets of oats that are closer to the potential heterotic groups. To achieve a more reliable and useful grouping it may be appropriate to test the K-values that are not the best fitted to ∆K. The ∆K helps in identifying the correct number of clusters in most situations, but it should not be used exclusively [30,31]. The ΔK method should be used with caution and a meaningful genetic structure should be considered when selecting an appropriate K-value [32]. Such an approach may allow the detection of additional substructure layer with more closely related accessions [33]. Considering these references, it was hypothesized that the formation of seven potential heterotic groups can be more applicable. The analysis conducted for K = 7, where the second highest peak was observed in ∆K plot (Figure 3), permitted a more detailed genetic grouping to seven clusters.
The alignment of individual clusters containing genotypes grouped to seven subgroups is marked in the Figure 4 (colors of clusters are the same as in the Figure 5). It reflected the microsatellite allele frequencies between subgroups (Table 3). The most distant were subgroups 3 and 7, and the most similar were 1 and 4. The average distances (expected heterozygosity) between individual genotypes within seven individual clusters indicated great variation within each of them (values are inside the Figure 4). Their values were quite comparable and showed relatively high average distances.
Seed parameters were no longer decisive for genotype classification in the presented output of the STRUCTURE analysis (Figure 5). The variation at the microsatellite loci was crucial; therefore, the STRUCTURE analysis did not present the commonly used form of clustering according to the affiliation with individual clusters. However, the grouping of genotypes according to their seed parameters (white, yellow, brown, naked) was intentionally applied. Such a presentation showed a high genetic variation within all four subgroups of oats (white, yellow, brown, and naked).
Oat genotypes with white and yellow seeds were distributed across all seven clusters. White oats had a high proportion in clusters 1 and 4 (red and yellow colors in Figure 4 and Figure 5), 24.3% and 29.7%, respectively. Yellow oats had a proportion of 28.8% in cluster 6 (turquoise) and brown oats of 43.9% proportion in cluster 4 (yellow). Naked oats were grouped mainly in cluster 7 (orange) with a proportion 42.1%. According to the frequency divergence of the analyzed microsatellite alleles, it was possible to identify the most genetically distant genotypes within individual white, yellow, brown, and naked subgroups. Among the white oat genotypes, the most distant were Edit versus Pendek, among the yellow oats Euro, Expander, Expo versus Consul, and among the naked oats genotypes Detvan versus PS-106. Within the hulled genotypes such can be Edmund, Flamingstern, and Unisignum versus Consul and Edit. For the selection of the most genetically different genotypes within the whole set of 85 oats, it would be appropriate to select genotypes from clusters 3 and 7. The classification created in this way made it possible to search the most different genotypes with identical seed parameters as well as between genotypes regardless of this parameter.
The calculated fixation indices (Fst) from the STRUCTURE analysis helped to know how different subgroups can be from each other. The overall fixation index Fst was 0.367, suggesting high genetic differentiation between the subgroups in the frame of the assayed germplasm and validated heterotic grouping suitable for development of advanced offspring. In addition, the Fst values of seven individual subgroups were: Fst1 = 0.181, Fst2 = 0.412, Fst3 = 0.498, Fst4 = 0.434, Fst5 = 0.214, Fst6 = 0.490, and Fst7 = 0.343, indicating high variation found within subgroups, especially in those included in clusters 2, 3, 4, and 6. This can suggest that the separation of genotypes according to genetic diversity at microsatellite loci can continue into even more heterotic groups. This is also supported by the ∆K plot (Figure 3) where another interesting peak was observed for K = 9.

3. Discussion

Melchinger and Gumber [34] in their concept of heterosis suggested approaches to identify heterotic groups. With many germplasm accessions, it is not feasible in most crops to form diallel crosses and produce sufficient F1 seeds for field testing in different environments. Their suggestion was to identify heterotic groups, and the first step was to group the germplasm (i.e., genotypes) according to genetic similarity. Molecular markers have been marked as extremely powerful tools for grouping of germplasm for a long time [35]. This theory was later tested in practical hybrid breeding of different plant species. Several authors have reported in different crops no or low correlations between heterosis and the genetic distance defined by molecular markers not linked to the required trait and higher heterosis effect for a yield reported in intra-group than for inter-group hybrids [36,37,38,39]. Others stated that some significant correlations between genetically different heterotic groups formed by molecular markers and yield in hybrids were found only in particular environments and across environments [40]. However, many other studies have confirmed the theory published by Melchinger [35] and Melchinger and Gumber [34]. Geng et al. [41] found that the crossing of parents from different heterotic groups, formed by SSR and SNP markers resulted in a heterotic effect in F1 progenies. They suggested that the genetic distance between parents determined by molecular markers can be helpful in heterosis prediction for some traits, and categorization for heterotic groups and parental selection would be beneficial in cotton hybrid breeding. Similar statements were presented in sorghum and maize [7,42,43], rye [44], barley [45], wheat [46], rice [47], and in other crops. Wang et al. [48] concluded that genomic selection predicted higher accuracy of genotypic value of genotypes for hybrid breeding but only if a higher number of field trials in more testing locations and years are used. According to them, genomic selection can be extended by marker-assisted selection using a higher number of molecular markers distributed over the genome with sufficient density. Therefore, a much simpler genetic distance analysis using different molecular markers, especially SSRs and SNPs, can be easily implemented as a simple predictor of hybrid performance in parental crosses before the implementation of such more accurate but more complicated genomic selection procedures [48,49,50].

4. Materials and Methods

4.1. Plant Material

The set of 85 oats (Avena sativa L.) contained registered cultivars and breeding lines from the period of 1952–2007 originating from 18 countries (including former Czechoslovakia and Soviet Union, Table S1). The oat genotypes were with hulled or naked (hull-less) grains (Avena sativa var. nuda Koern) and differed in the color of the glume (white, yellow, brown). Seeds of all oat genotypes were obtained from the Avena sp. collection of genetic resources maintained in the Gene Bank of the Slovak Republic (Piešťany, Slovakia).

4.2. Molecular Analysis

The total DNA was extracted from young leaves using the protocol of Dellaporta et al. [51]. DNA from each oat genotype represented bulk DNA prepared from equivalent amounts of DNA from 1 to 2 g of young leaves. Twenty pairs of microsatellite-derived primers were used for analysis (Table S2). Polymerase chain reactions were performed in 15 μL volumes containing 25 ng of DNA, 1 × PCR buffer (50 mmol/L KCl, 10 mmol/L TRIS-HCl, pH 8.3), 1.5 mmol/L MgCl2, 0.1 μmol/L of both primers, 0.1 μmol/L each of dNTPs, and 0.8 U of Taq-DNA polymerase. The amplification conditions included an initial denaturation 3 min at 94 °C, followed by 30 cycles of denaturation, 1 min at 94 °C, annealing 1 min at annealing temperature (Table S2), extension 1 min at 72 °C, and the final extension 8 min at 72 °C. The equivalent volume of loading buffer (95% formamide, 10 mmol/L NaOH, 0.05% bromophenol blue) was added to each sample. Samples were denatured for 3 min at 99 °C and 5 μL of each sample was loaded into 6% polyacrylamide gel containing 7 mol/L of urea. Gels in the electrophoretic unit (SEQ 3341, Scie-Plas Ltd., Waterbeach Cambridge, United Kingdom) were run in 0.5 × TBE buffer at the constant power of 45 W for 3.5–5 h, depending on the size of amplified fragments. The microsatellite DNA was stained by the silver staining method [52]. The sizes of the microsatellite alleles were determined using 10 bp, 25 bp, 50 bp, and 100 bp DNA Ladders (Invitrogen Thermo Fisher Scientific, Waltham, MA, USA).

4.3. Data Analysis

The selected polymorphic indices characterizing used SSR markers and loci were calculated using the iMEC Online Marker Efficiency Calculator [21].
The principal component analysis (PCA), using the Euclidean distance measure and the cluster analysis using the neighbor joining clustering and Jaccard similarity index were performed using the Paleontological Statistics (PAST) software version 3.19 [53]. This software also used the Euclidean similarity index and 9999 permutations for the Pairwise PERMANOVA analysis using score (eigenvalues) from the first 10 significant principal components based on the Scree plot in PCA.
The population structure was evaluated through the STRUCTURE v. 2.3.4 software [54] using the default setting of the admixture model for the ancestry of individuals and correlated allele frequencies. The models were tested for K-values ranging from one to fifteen with ten independent runs each. The Burn-in and Markov Chain Monte Carlo iterations were set to 100,000. The number of clusters was chosen by plotting the LnP(D) values against ∆K values with the K value selected according to the Evanno test [30]. The tree in STRUCTURE was estimated using the program NEIGHBOR by Mary Kuhner and John Yamato, implementing Saitou and Nei’s neighbor joining method [55]. The plot was produced using DRAWTREE as part of his PHYLIP phylogeny package [56,57,58,59].

5. Conclusions

The analysis of genetic distance using genomic microsatellite markers revealed a relatively high degree of polymorphism and heterozygosity at given loci. This can be used for the almost perfect mutual differentiation of genotypes as well as for the classification of genotypes into several subgroups. Microsatellite polymorphism dominated in the way these subgroups were discriminated, but signs of grouping by seed parameters were also identifiable. High genetic differentiation between the subgroups supported a heterotic grouping. Moreover, the high heterozygosity found within subgroups suggests that microsatellite polymorphism can also be used in the formation of other heterotic groups. This approach can be implemented as a simple predictor of heterosis in parental crosses prior to extensive field testing or the development and implementation of more accurate genomic selection.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10112462/s1, Table S1: List of oat genotypes included in SSR analyses, Table S2: Primer sequences for oat microsatellite loci.

Author Contributions

Conceptualization, J.K. and M.H.; methodology, J.K., M.H. and D.M.; plant material, P.H.; laboratory experiments, M.H. and J.K.; data analysis, M.H., J.K. and K.O.; writing (draft preparation, review and editing), J.K., M.H., K.O. and V.G.; graphical visualization, K.O.; supervision, J.K.; project administration and funding acquisition, M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovak Research and Development Agency, grant number APVV-18-0154 and by the Operational Programme Research and Development co-financed by the European Regional Development Fund, grant number ITMS 26210120039.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bityutskii, N.P.; Loskutov, I.; Yakkonen, K.; Konarev, A.; Shelenga, T.; Khoreva, V.; Blinova, E.; Ryumin, A. Screening of Avena sativa cultivars for iron, zinc, manganese, protein and oil content and fatty acid composition in whole grains. Cereal Res. Commun. 2020, 48, 87–94. [Google Scholar] [CrossRef]
  2. Peterson, D.M. Oat antioxidants. J. Cereal Sci. 2001, 33, 115–129. [Google Scholar] [CrossRef]
  3. Havrlentová, M.; Kraic, J. Content of beta-D-glucan in cereal grains. J. Food Nutr. Res. 2006, 45, 97–103. [Google Scholar]
  4. Redaelli, R.; Del Frate, V.; Bellato, S.; Terracciano, G.; Ciccoritti, R.; Germeier, C.U.; De Stefanis, E.; Sgrulletta, D. Genetic and environmental variability in total and soluble β-glucan in European oat genotypes. J. Cereal Sci. 2013, 57, 193–199. [Google Scholar] [CrossRef]
  5. Ahmad, M.; Gul-Zaffar Dar, Z.A.; Habib, M. A review of oat (Avena sativa L.) as a dual-purpose crop. Sci. Res. Essays 2014, 9, 52–59. [Google Scholar] [CrossRef] [Green Version]
  6. Melchinger, A.E. Genetic diversity and heterosis. In The Genetics and Exploitation of Heterosis in Crops; Coors, J.G., Pandey, S., Eds.; American Society of Agronomy Crop Science Society of America: Madison, WI, USA, 1999; pp. 99–118. [Google Scholar]
  7. Tomkowiak, A.; Bocianowski, J.; Kwiatek, M.; Kowalczewski, P.Ł. Dependence of the heterosis effect on genetic distance, determined using various molecular markers. Open Life Sci. 2020, 15, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Ahmed, S.; Roy, A.K.; Majumdar, A.B. Genetic diversity and variability analysis in oat (Avena sativa L.). Range Manag. Agrofor. 2011, 32, 96–99. [Google Scholar]
  9. Krishna, A.; Ahmed, S.; Pandey, H.C.; Kumar, V. Correlation, path and diversity analysis of oat (Avena sativa L.) genotypes for grain and fodder yield. J. Plant Sci. Res. 2014, 1, 110. [Google Scholar]
  10. O’Donoughue, L.S.; Souza, E.; Tanksley, S.D.; Sorrells, M.E. Relationships among North American oat cultivars based on restriction fragment length polymorphisms. Crop Sci. 1994, 34, 1251–1258. [Google Scholar] [CrossRef]
  11. Tinker, N.A.; Kilian, A.; Wight, C.P.; Heller-Uszynska, K.; Wenzl, P.; Rines, H.W.; Bjørnstad, Å.; Howarth, C.J.; Jannink, J.-L.; Anderson, J.M.; et al. New DArT markers for oat provide enhanced map coverage and global germplasm characterization. BMC Genom. 2009, 10, 39. [Google Scholar] [CrossRef] [Green Version]
  12. Boczkowska, M.; Nowosielski, J.; Nowosielska, D.; Podyma, W. Assessing genetic diversity in 23 early Polish oat cultivars based on molecular and morphological studies. Genet. Res. Crop Evol. 2014, 61, 927–941. [Google Scholar] [CrossRef] [Green Version]
  13. Cömertpay, G.; Habyarimana, E.; Baloch, F.S.; Güngör, H.; Dokuyucu, T.; Akkaya, A.; Dumlupinar, Z. Geographical description and molecular characterization of genetic structure and diversity using 6K SNP array in Turkish oat germplasm. Can. J. Plant Sci. 2019, 99, 12–21. [Google Scholar] [CrossRef]
  14. Jan, S.F.; Khan, M.R.; Iqbal, A.; Khan, F.U.; Ali, S. Genetic diversity in exotic oat germplasm and resistance against barley yellow dwarf virus. Saudi J. Biol. Sci. 2020, 27, 2622–2631. [Google Scholar] [CrossRef]
  15. Yan, H.; Zhou, P.; Peng, Y.; Bekele, W.A.; Ren, C.; Tinker, N.A.; Peng, Y. Genetic diversity and genome-wide association analysis in Chinese hulless oat germplasm. Theor. Appl. Genet. 2020, 133, 3365–3380. [Google Scholar] [CrossRef]
  16. Montilla-Bascón, G.; Sánchez-Martín, J.; Rispail, N.; Mur, L.; Langdon, T.; Griffiths, I.; Howarth, C.; Prats, E. Genetic diversity and population structure among oat cultivars and landraces. Plant Mol. Biol. Rep. 2013, 31, 1305–1314. [Google Scholar] [CrossRef] [Green Version]
  17. Dumlupinar, Z.; Brown, R.; Campbell, R.; Jellen, E.N.; Anderson, J.; Bonman, J.M.; Carson, M.; Chao, S.; Obert, D.; Jackson, E. The art of attrition: Development of robust oat microsatellites. Plant Breed. 2020, 135, 323–334. [Google Scholar] [CrossRef]
  18. Leišová-Svobodová, L.; Michel, S.; Tamm, I.; Chourová, M.; Janovska, D.; Grausgruber, H. Diversity and pre-breeding prospects for local adaptation in oat genetic resources. Sustainability 2019, 11, 6950. [Google Scholar] [CrossRef] [Green Version]
  19. Powell, W.; Morgante, M.; Andre, C.; Hanafey, M.; Vogel, J.; Tingey, S.; Rafalski, A. The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Mol. Breed. 1996, 2, 225–238. [Google Scholar] [CrossRef]
  20. Milbourne, D.; Meyer, R.; Bradshaw, J.E.; Baird, E.; Bonar, N.; Provan, J.; Powell, W.; Waugh, R. Comparison of PCR-based marker systems for the analysis of genetic relationships in cultivated potato. Mol. Breed. 1997, 3, 127–136. [Google Scholar] [CrossRef]
  21. Amiryousefi, A.; Hyvönen, J.; Poczai, P. iMEC: Online marker efficiency calculator. Appl. Plant Sci. 2018, 6, e1159. [Google Scholar] [CrossRef]
  22. Prevost, A.; Wilkinson, M.J. A new system of comparing PCR primers applied to ISSR fingerprint of potato cultivars. Theor. Appl. Genet. 1999, 98, 107–112. [Google Scholar] [CrossRef]
  23. Tessier, C.J.; David, P.; This, J.; Boursiquot, M.; Charrier, A. Optimization of the choice of molecular markers for varietal identification in Vitis vinifera L. Theor. Appl. Genet. 1999, 98, 171–177. [Google Scholar] [CrossRef]
  24. Nikoloudakis, N.; Bladenopoulos, K.; Katsiotis, A. Structural patterns and genetic diversity among oat (Avena) landraces assessed by microsatellite markers and morphological analysis. Genet. Res. Crop Evol. 2016, 63, 801–811. [Google Scholar] [CrossRef]
  25. Kapoor, R.; Choudhary, K. Genetic diversity analysis of fodder oats (Avena sativa L.) germplasm by microsatellite markers. J. Agric. Sci. Technol. 2017, 19, 1369–1379. [Google Scholar]
  26. Bräutigam, M.; Lindlöf, A.; Zakhrabekova, S.; Gharti-Chhetri, G.; Olsson, B.; Olsson, O. Generation and analysis of 9792 EST sequences from cold acclimated oat, Avena sativa. BMC Plant Biol. 2005, 5, 18. [Google Scholar] [CrossRef]
  27. Becher, R. EST-derived microsatellites as a rich source of molecular markers for oats. Plant Breed. 2007, 126, 274–278. [Google Scholar] [CrossRef]
  28. Isabel, L.P.; Park, J.-R.; Lee, G.S.; Park, G.-H.; Kim, K.-M. Development of EST-SSR markers and analysis of genetic relationship it’s resources in hexaploid oats. J. Crop Sci. Biotechnol. 2019, 22, 243–251. [Google Scholar] [CrossRef]
  29. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  30. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  31. Wegary, D.; Teklewold, A.; Prasanna, B.M.; Ertiro, B.T.; Alachiotis, N.; Negera, D.; Awas, G.; Abakemal, D.; Ogugo, V.; Gowda, M.; et al. Molecular diversity and selective sweeps in maize inbred lines adapted to African highlands. Sci. Rep. 2019, 9, 13490. [Google Scholar] [CrossRef] [Green Version]
  32. Cullingham, C.; Miller, J.; Peery, R.; Dupuis, J.; Malenfant, R.; Gorrell, J.; Janes, J. Confidently identifying the correct K value using the ΔK method: When does K = 2? Mol. Ecol. 2020, 29, 862–869. [Google Scholar] [CrossRef] [PubMed]
  33. Roncallo, P.F.; Beaufort, V.; Larsen, A.O.; Dreisigacker, S.; Echenique, V. Genetic diversity and linkage disequilibrium using SNP (KASP) and AFLP markers in a worldwide durum wheat (Triticum turgidum L. var durum) collection. PLoS ONE 2019, 14, e0218562. [Google Scholar] [CrossRef] [PubMed]
  34. Melchinger, A.E.; Gumber, R.K. Overview of Heterosis and Heterotic Groups in Agronomic Crops. In Concepts and Breeding of Heterosis in Crop Plants, 1st ed.; Larnkey, K.R., Staub, J.E., Eds.; Crop Science Society of America: Madison, WI, USA, 1998; pp. 29–44. [Google Scholar]
  35. Melchinger, A.E. Use of RFLP markers for analysis of genetic relationship among breeding materials and prediction of hybrid performance. In International Crop Science I; Buxton, D.R., Shibles, R., Forsberg, R.A., Blad, B.L., Asay, K.H., Paulsen, G.M., Wilson, R.F., Eds.; Crop Science Society of America: Madison, WI, USA, 1993; pp. 621–628. [Google Scholar]
  36. Kwon, S.-J.; Ha, W.-G.; Hwang, H.-G.; Yang, S.-J.; Choi, H.-C.; Moon, H.-P.; Ahn, S.-N. Relationship between heterosis and genetic divergence in Tongil-type rice. Plant Breed. 2002, 121, 487–492. [Google Scholar] [CrossRef]
  37. Tams, S.H.; Bauer, E.; Oettler, G.; Melchinger, A.E.; Schön, C.-C. Prospects for hybrid breeding in winter triticale: II. Relationship between parental genetic distance and specific combining ability. Plant Breed. 2006, 125, 331–336. [Google Scholar] [CrossRef]
  38. Dermail, A.; Suriharn, B.; Chankaew, S.; Sanitchon, J.; Lertrat, K. Hybrid prediction based on SSR-genetic distance, heterosis and combining ability on agronomic traits and yields in sweet and waxy corn. Sci. Hortic. 2020, 259, 108817. [Google Scholar] [CrossRef]
  39. Kumar, P.; Sarma, D.; Hazarika, M. Characterization and heterotic grouping of traditional Assam rice (Oryza sativa L.). bioRxiv 2021. [Google Scholar] [CrossRef]
  40. Xie, F.; He, Z.; Esguerra, M.Q.; Qiu, F.; Ramanathan, V. Determination of heterotic groups for tropical Indica hybrid rice germlasm. Theor. Appl. Genet. 2014, 127, 407–417. [Google Scholar] [CrossRef]
  41. Geng, X.; Qu, Y.; Jia, Y.; He, S.; Pan, Z.; Wang, L.; Du, X. Assessment of heterosis based on parental genetic distance estimated with SSR and SNP markers in upland cotton (Gossypium hirsutum L.). BMC Genom. 2021, 22, 123. [Google Scholar] [CrossRef]
  42. Amelework, B.; Abakemal, D.; Shimelis, H.; Laing, M. Application of microsatellites in genetic diversity analysis and heterotic grouping of sorghum and maize. In Microsatellite Markers; Abdurakhmonov, I.Y., Ed.; InTech: Rijeka, Slovenia, 2016; pp. 117–138. [Google Scholar]
  43. Van Inghelandt, D.; Melchinger, A.E.; Lebreton, C.; Stich, B. Population structure and genetic diversity in a commercial maize breeding program assessed with SSR and SNP markers. Theor. Appl. Genet. 2010, 120, 1289–1299. [Google Scholar] [CrossRef] [Green Version]
  44. Fischer, S.; Melchinger, A.E.; Korzun, V.; Wilde, P.; Schmiedchen, B.; Möhring, J.; Piepho, H.-P.; Dhillon, B.S.; Würschum, T.; Reif, J.C. Molecular marker assisted broadening of the Central European heterotic groups in rye with Eastern European germplasm. Theor. Appl. Genet. 2010, 120, 291–299. [Google Scholar] [CrossRef] [PubMed]
  45. Sommer, L.; Spiller, M.; Stiewe, G.; Pillen, K.; Reif, J.C.; Schulthess, A.W. Proof of concept to unmask the breeding value of genetic resources of barley (Hordeum vulgare) with a hybrid strategy. Plant Breed. 2019, 139, 536–549. [Google Scholar] [CrossRef]
  46. Al-Ashkar, I.; Alotaibi, M.; Refay, Y.; Ghazy, A.; Zakri, A.; Al-Doss, A. Selection criteria for high-yielding and early-flowering bread wheat hybrids under heat stress. PLoS ONE 2020, 15, e0236351. [Google Scholar] [CrossRef] [PubMed]
  47. Sruthi, K.; Divya, B.; Senguttuvel, P.; Revathi, P.; Kemparaju, K.B.; Koteswararao, P.; Sundaram, R.M.; Singh, V.J.; Ranjith Kumar, E.; Bhowmick, P.K.; et al. Evaluation of genetic diversity of parental lines for development of heterotic groups in hybrid rice (Oryza sativa L.). J. Plant Biochem. Biotechnol. 2020, 29, 236–252. [Google Scholar] [CrossRef]
  48. Wang, Y.; Mette, M.F.; Miedaner, T.; Gottwald, M.; Wilde, P.; Reif, J.C.; Zhao, Y. The accuracy of prediction of molecular marker assisted selection and genomic selection in elite hybrid rye populations surpasses the accuracy of marker-assisted selection and is equally augmented by multiple field evaluation locations and test years. BMC Genom. 2014, 15, 556. [Google Scholar] [CrossRef] [Green Version]
  49. Vendelbo, N.M.; Sarup, P.; Orabi, J.; Kristensen, P.S.; Jahoor, A. Genetic structure of a germplasm for hybrid breeding in rye (Secale cereale L.). PLoS ONE 2020, 15, e0239541. [Google Scholar] [CrossRef]
  50. Auinger, H.J.; Schönleben, M.; Lehermeier, C.; Schmidt, M.; Korzun, V.; Geiger, H.H.; Piepho, H.-P.; Gordillo, A.; Wilde, P.; Bauer, E.; et al. Model training across multiple breeding cycles significantly improves genomic prediction accuracy in rye (Secale cereale L.). Theor. Appl. Genet. 2016, 129, 2043–2053. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Dellaporta, S.L.; Wood, J.; Hicks, J.B. A plant DNA minipreparation, version II. Plant Mol. Biol. Rep. 1993, 4, 19–21. [Google Scholar] [CrossRef]
  52. Bassam, B.J.; Caetano-Anolles, G.; Gresshoff, P.M. Fast and sensitive silver staining of DNA in polyacrylamide gels. Anal. Biochem. 1991, 196, 80–83. [Google Scholar] [CrossRef]
  53. Hammer, Ø.; Harper, D.A.T. Paleontological Data Analysis; Blackwell Publishing Ltd.: Oxford, UK, 2006; 351p. [Google Scholar]
  54. Earl, D.A.; vonHoldt, B.M. Structure Harvester: A website and program for visualizing Structure output and implementing the Evanno method. Conserv. Gene. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  55. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef]
  56. Felsenstein, J. PHYLIP—Phylogeny Inference Package (Version 3.2). Cladistics 1989, 5, 164–166. [Google Scholar]
  57. Li, C.D.; Rossnagel, B.G.; Scoles, G.J. The development of oat microsatellite markers and their use in identifying relationships among Avena species and oat cultivars. Theor. Appl. Genet. 2000, 101, 1259–1268. [Google Scholar] [CrossRef]
  58. Pal, N.; Sandhu, J.S.; Domier, L.L.; Kolb, F.L. Development and characterization of microsatellite and RFLP-derived PCR markers in oat. Crop Sci. 2002, 42, 912–918. [Google Scholar] [CrossRef]
  59. Portyanko, V.A.; Hoffman, D.L.; Lee, M.; Holland, J.B. A linkage map of hexaploid oat based on grass anchor DNA clones and its relationship to other oat maps. Genome 2001, 44, 249–265. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The principal component analysis (PCA) of microsatellite data and separation of oats with different seed parameters (white, yellow, brown, naked).
Figure 1. The principal component analysis (PCA) of microsatellite data and separation of oats with different seed parameters (white, yellow, brown, naked).
Plants 10 02462 g001
Figure 2. Cluster analysis of oat genotypes according to microsatellite variation.
Figure 2. Cluster analysis of oat genotypes according to microsatellite variation.
Plants 10 02462 g002
Figure 3. The ln(K) and ∆K values based on the dataset of used microsatellites. Hypothesized number of populations ranged from 1 to 14.
Figure 3. The ln(K) and ∆K values based on the dataset of used microsatellites. Hypothesized number of populations ranged from 1 to 14.
Plants 10 02462 g003
Figure 4. Seven clusters of grouped 85 oats according to differences in expected heterozygosity (colors of numbers correlate with colors of clusters in the Figure 5).
Figure 4. Seven clusters of grouped 85 oats according to differences in expected heterozygosity (colors of numbers correlate with colors of clusters in the Figure 5).
Plants 10 02462 g004
Figure 5. Classification of 85 oat genotypes based on microsatellite variation and seed parameters. Each bar represented individual genotype numbered in accordance with Table S1.
Figure 5. Classification of 85 oat genotypes based on microsatellite variation and seed parameters. Each bar represented individual genotype numbered in accordance with Table S1.
Plants 10 02462 g005
Table 1. Parameters of polymorphism in analyzed SSR loci.
Table 1. Parameters of polymorphism in analyzed SSR loci.
LocusnHPICEHavpMIDR
AM1100.32700.37202.05880.00040.00080.95784.1176
AM1480.37790.35402.02350.00060.00110.93632.5176
AM2230.49870.30111.57650.00200.00310.72481.8353
AM8320.50000.30041.00000.00290.00290.75150.1412
AM8740.37500.35511.00000.00110.00110.93810.9882
AM10240.37500.35511.00000.00110.00110.93810.9412
AM11540.49920.30082.08240.00150.00310.72970.2118
Mean50.42180.33411.53450.00140.00190.85381.5361
n—number of alleles; H—heterozygosity index; PIC—polymorphism information content; E—effective multiplex ratio; Havp—mean heterozygosity; MI—marker index; D—discriminating power; R—resolving power.
Table 2. p-values from Pairwise PERMANOVA between all pairs of subgroups (white, yellow, brown, naked) calculated from the values of the first 10 significant principal components from PCA. Significant comparisons at p ≤ 0.05 are shown in bold.
Table 2. p-values from Pairwise PERMANOVA between all pairs of subgroups (white, yellow, brown, naked) calculated from the values of the first 10 significant principal components from PCA. Significant comparisons at p ≤ 0.05 are shown in bold.
WhiteYellowBrownNaked
White-
Yellow0.0414-
Brown0.52600.1204-
Naked0.00780.00140.1103-
Table 3. Divergences in allele frequency among oat subgroups differentiated into 7 clusters.
Table 3. Divergences in allele frequency among oat subgroups differentiated into 7 clusters.
Cluster 1Cluster 2Cluster 3Cluster 4Cluster 5Cluster 6Cluster 7
Cluster 1-
Cluster 20.228-
Cluster 30.2020.232-
Cluster 40.0550.2290.200-
Cluster 50.0770.2020.2010.063-
Cluster 60.1110.2630.1780.1110.081-
Cluster 70.1880.2450.2740.2100.1470.179-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Havrlentová, M.; Ondreičková, K.; Hozlár, P.; Gregusová, V.; Mihálik, D.; Kraic, J. Formation of Potential Heterotic Groups of Oat Using Variation at Microsatellite Loci. Plants 2021, 10, 2462. https://doi.org/10.3390/plants10112462

AMA Style

Havrlentová M, Ondreičková K, Hozlár P, Gregusová V, Mihálik D, Kraic J. Formation of Potential Heterotic Groups of Oat Using Variation at Microsatellite Loci. Plants. 2021; 10(11):2462. https://doi.org/10.3390/plants10112462

Chicago/Turabian Style

Havrlentová, Michaela, Katarína Ondreičková, Peter Hozlár, Veronika Gregusová, Daniel Mihálik, and Ján Kraic. 2021. "Formation of Potential Heterotic Groups of Oat Using Variation at Microsatellite Loci" Plants 10, no. 11: 2462. https://doi.org/10.3390/plants10112462

APA Style

Havrlentová, M., Ondreičková, K., Hozlár, P., Gregusová, V., Mihálik, D., & Kraic, J. (2021). Formation of Potential Heterotic Groups of Oat Using Variation at Microsatellite Loci. Plants, 10(11), 2462. https://doi.org/10.3390/plants10112462

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop