Next Article in Journal
Phosphorus Deficiency Alters Nutrient Accumulation Patterns and Grain Nutritional Quality in Rice
Next Article in Special Issue
Development of a Molecular Breeding Strategy for the Integration of Transgenic Traits in Outcrossing Perennial Grasses
Previous Article in Journal
Commercial Humates in Agriculture: Real Substance or Smoke and Mirrors?
Previous Article in Special Issue
Evidence for Heterosis in Italian Ryegrass (Lolium multiflorum Lam.) Based on Inbreeding Depression in F2 Generation Offspring from Biparental Crosses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Trait Loci (QTL) Identification in the Progeny of a Polycross

1
Jouffray-Drillaud La Litière, Saint Sauvant 86600, France
2
GIE GRASS La Litière, Saint Sauvant, 86600 France
3
INRA, UR4, Unité de Recherche Pluridisciplinaire Prairies et Plantes Fourragères, RD150, Le Chêne, CS80006, Lusignan 86600, France
*
Author to whom correspondence should be addressed.
Agronomy 2016, 6(4), 51; https://doi.org/10.3390/agronomy6040051
Submission received: 23 September 2016 / Revised: 18 October 2016 / Accepted: 20 October 2016 / Published: 28 October 2016

Abstract

:
Connected multiparental crosses are valuable for detecting quantitative trait loci (QTL) with multiple alleles. The objective of this study was to show that the progeny of a polycross can be considered as connected mutiparental crosses and used for QTL identification. This is particularly relevant in outbreeding species showing strong inbreeding depression and for which synthetic varieties are created. A total of 191 genotypes from a polycross with six parents were phenotyped for plant height (PH) and plant growth rate (PGR) and genotyped with 82 codominant markers. Markers allowed the identification of the male parent for each sibling and so the 191 genotypes were divided into 15 full-sib families. The number of genotypes per full-sib family varied from 2 to 28. A consensus map of 491 cM was built and QTL were detected with MCQTL-software dedicated to QTL detection in connected mapping populations. Two major QTL for PH and PGR in spring were identified on linkage groups 3 and 4. These explained from 12% to 22% of phenotypic variance. The additive effects reached 12.4 mm for PH and 0.11 mm/C°d for PGR. This study shows that the progenies of polycrosses can be used to detect QTL.

1. Introduction

Molecular markers are being used increasingly in plant breeding either to construct new genotypes having favourable alleles or to better estimate the breeding values of genotypes [1,2]. The first option implies the detection of the loci involved in the variation of quantitative traits (QTL). This detection is usually realised either in populations derived from bi-parental crosses between contrasting parents or in populations for which the lineage is unknown (association studies). The first strategy is interesting because linkage disequilibrium (LD) is long in the populations created, allowing the use of a limited number of markers to cover the genome. Nevertheless, this strategy can be limited by difficulties in creating the populations by the low number of alleles surveyed and by the accuracy of the location of the QTL. The second strategy (association studies) is interesting because the populations already exist and can be highly diverse. Nevertheless, it can be limited by short LD, forcing breeders to focus on just a few genes (candidate gene approach) or to use a very large number of markers [3]. In particular, LD is often very short in outbreeding species [3].
Synthetic varieties are produced in outbreeding species for which crosses cannot be controlled at large scale to produce hybrid varieties (e.g., in the majority of perennial forage and tree species). Synthetic varieties are created by inter-crossing a number of selected plants (from four to several hundred parents (i.e., making a polycross)) and then multiplying up the resulting population over three or four generations (Syn 3–4), without selection, to obtain sufficient seed for commercialisation. The genetic diversity remaining in synthetic varieties avoids inbreeding depression and also allows better stability of performance in variable environments (spatial and temporal variations) [4]. When a low number of parents is used, these can be selected for their specific combining ability in addition to their general combining ability [5]. To apply marker assisted selection (MAS) in such species, QTL can be detected in bi-parental populations and then the favourable alleles must be introgressed into more diverse populations, such as elite varieties. This process is time consuming and the effect of the QTL can change in a new genetic background. Another way is to detect QTL directly in diverse populations such as in the progeny of a polycross or in a given generation of multiplication of a synthetic variety (Syn 2, 3 or 4). Linkage disequilibrium decreases with successive generations of multiplication [6]. Depending on the genetic architecture of selected traits, the knowledge of candidate genes and the budgetary allocation for genotyping (number of markers necessary to cover the genome depending on the LD) breeders can choose the best population for MAS. For example, a Syn 4 has been used to detect an association between a candidate gene and leaf elongation rate in perennial ryegrass [7].
Perennial ryegrass (Lolium perenne L.) is the most commonly sown forage and turf-grass species in temperate climates and it is considered to be a model for genomics in forage grasses [8]. Varieties are synthetics, due to biological constraints. In this species, LD decreases very rapidly (r2 < 0.2 over less than 1 kb) in natural populations but also in synthetic varieties, except for in varieties produced from a very small number of parental plants [9,10]. Leaf length and leaf elongation rate are important traits affecting: (i) vegetative yield [11,12,13,14,15]; (ii) intake rate by dairy cows [16]; and (iii) plant survival under light competition conditions [17]. At the plant level, leaf length can be estimated by stretched plant height (PH) and leaf elongation rate by plant growth rate (PGR).
The objective of this study was to demonstrate that the progeny from a polycross can be used advantageously for the detection of QTL with interval mapping because (i) the variability within the progeny can be sufficient to avoid inbreeding depression in the following step of selection to produce a variety; (ii) LD is high, allowing the use of a moderate number of markers to cover the genome; and (iii) polycrosses are classically made in breeding programs for species in which synthetics are produced. Our approach was to identify QTL in the progeny of a polycross of perennial ryegrass including six parents, for PH and for PGR, the increase in plant height over a certain time after a defoliation event.

2. Results

2.1. Paternity Identification

A total of 82 codominant molecular markers allowed the identification of the male parent for the 191 genotypes (Table 1). The number of genotypes per full-sib family ranged from 2 to 28, clearly showing that the crosses within the polycross were not in panmixia. That is, Nemo A, Nemo B, Nemo D, and Nemo G tend to cross preferentially together, except for Nemo B and Nemo D. Also, Nemo C and Nemo H tend to cross preferentially together but there was a deficit of genotypes in all other crosses involving either Nemo C or Nemo H. In addition, few genotypes were obtained from Nemo B × Nemo D. This non-random mating could be at least partly due to differences in flowering date. Indeed, Nemo C and Nemo H showed heading dates three to four days earlier than the other potential parents (Table 1). Another reason for low mating between Nemo A × Nemo C, Nemo B × Nemo D, and Nemo D × Nemo H could have arisen from the spatial distance between parents in the polycross (Figure 1). These pairs were not nearby in any of the four repetitions of the polycross.

2.2. Phenotypic Analyses

The 191 genotypes were measured for PH in spring and autumn in a spaced-plant trial with four replicates. Six measurements of PH in April were used to estimate PGR (Figure 2). After the defoliation of all plants on 24 March, PH increased linearly up to a maximum reached after three to six weeks, depending on the plant. This means that at this date (24 March) all plants had re-started their growth after winter. The fast regrowth indicates all the plants were certainly induced for flowering [18].
For each trait, the distribution of the error term of the variance analysis was not significantly different from a normal distribution. The broad sense heritability was medium for PH and low for PGR (Table 2).
All traits showed large genetic variability (i.e., the variability of the mean per genotype (Table 3)). On average, PH reached 205 mm after six weeks of regrowth in spring and 203 mm after 11 weeks of regrowth in autumn.
All PH measurements done in spring were highly inter-correlated and were also highly correlated with the PGR estimated in spring (Table 4). The high correlation between PH in March and in April indicates the majority of the plants had already started their regrowth in March, after winter. As expected from the estimation of PGR, PGR was very highly correlated with the first three PH measurements in April. The value of PGR was also highly correlated with the last two measurements in April, although slightly less than with the first three measurements. This indicates the rate of plant regrowth in spring had a greater influence on PH after six weeks of regrowth than the duration of the linear phase of regrowth. The PH in autumn was only poorly correlated with PH and PGR in spring.

2.3. Genetic Map

A consensus map was built based on the recombination information of the six parents included in the polycross (Figure 3). The map included seven linkage groups (LG) as expected in perennial ryegrass. The LG were numbered as in the map built by the International Lolium Genome Initiative [19]. It covered a distance of 491.4 cM with an average distance between two consecutive markers of 5.1 cM. The maximum distance between two markers was 36.9 cM on LG 2. LG 1 was surprisingly short, indicating a lack of markers on this LG.

2.4. QTL Identification

Two loci (LG 3 and 4) were identified with significant effects for PH and PGR in spring for all parents (Table 5). These explained between 12% and 22% of the variance of traits. The locus on LG 3 was detected only at the beginning of spring. No QTL was detected for PH in autumn. Significant effects were identified in all parents with Nemo C, showing the highest number of significant effects. The most favourable alleles were identified in Nemo C, Nemo D, and Nemo G on LG3 and Nemo B, Nemo C, and Nemo H on LG4.

3. Discussion

In this study, we successfully carried out genetic mapping and QTL detection on the progeny of a polycross, considered as a sum of connected populations. In the progeny, molecular markers allowed the identification of the male parent for each sibling. This approach has previously been used in red clover and alfalfa allowing increase of the selection gain [20]. As expected, we found the progeny of the polycross with six parents consisted of 15 full-sib families. However, the number of plants per full-sib family was highly variable with some families comprising fewer than five plants. This is a limitation for the estimation of a trait mean for very small families. To avoid this problem in future studies, it would be best to genotype a larger number of plants per half-sib family than the one needed for phenotyping and then to select seedlings based on their genotype so as to have equal numbers of plants per full-sib family. Another point that could be improved is the spatial design of the polycross, which should allow all crosses between the constituents of the polycross [21].
As expected in perennial ryegrass, we were able to build a genetic map with seven linkage groups, but the map was shorter than expected with only 491 cM, instead of between 700 and 800 cM [19,22,23,24,25]. The belonging of markers to LG was in agreement with previous studies [22,26,27,28]. The order of the markers on the LG was generally identical to the one in [27,28] except for LG3 and some punctual changes. For example, on LG2, the upper part of the LG, including G01_040 on the map in the present study, was at the opposite extremity in the map of [28]. The markers used in the present study cover the entire map of [28], but the size of the map was shrunk. For example, on LG1 the distance between LMgSSR04-09F and LpSSR085 was 33.5 cM in [28] and only 6.3 cM in this study. This could be due to the limited number of plants within each half-sib family and/or the number of markers in this study. In future, the number of markers used for this kind of study should be increased. This is now feasible at reasonable cost with marker platforms such as Genotyping by Sequencing (GBS) [29,30].
Values of QTL with strong effects were detected on LG 3 and 4 for PH and PGR in the spring. These regions have already been identified as QTL involved in traits related to leaf growth, such as leaf length, leaf elongation rate, and plant height [22,24,31,32,33,34,35]. In particular, the QTL on LG4 nearby LpSSR011 and G05_014 is located on the map of [22] nearby a QTL for leaf elongation rate, lamina length, and plant height. Moreover this QTL was near LpSSR082 which is close to the gene Gibberellic Acid Insensitive [9], which is involved in leaf elongation rate and leaf length [7]. The QTL on LG3 was not far from a QTL for leaf elongation rate in [32] (marker pps0164 in common) and in [34] (marker 25ca1 in common), but the very low number of common markers between studies does not allow a real accurate comparison. In the literature cited above, QTL were identified on all seven LG with the percentages of phenotypic variance explained by the QTL ranging from 5% to 43%. The fact that we detected only two strong QTL could arise from (i) the very small size of some full-sib families; (ii) a relatively low number of plants surveyed in comparison to the number of parents [36]; (iii) a partial genetic map; and (iv) the change of the effect of a QTL with genetic background. Nevertheless, the results obtained in this study are very encouraging for future work with larger populations and higher numbers of markers. It would improve the robustness of QTL to phenotype the plants on a multi-local trial, but due to the cost it could be a good option to keep the multi-local trials for advanced elite material tested on swards.
This study shows it was possible to detect QTL in the progeny of a polycross. This finding should be useful for breeding synthetic varieties because it allows maintenance of enough diversity to avoid inbreeding, while still selecting for the best alleles from different parents. Practically, from the progeny of a polycross, we propose to perform 1/ a selection based on molecular markers and on phenotypes (cycle 1) then 2/ a cycle of selection only on molecular markers (cycle 2) and 3/ a selection on seedlings based on molecular markers followed by a selection based on phenotypes on the rest of plants (cycle 3) (Figure 4). The selected plants could then be used to create a variety or to create new polycrosses for further cycles of selection. This breeding scheme could be done on several polycrosses at the same time and new polycrosses could be created with selected plants from different polycrosses. The selection on molecular markers should be performed for both increasing the frequencies of favourable alleles and for maintaining diversity 1/ on the rest of the genome and if possible 2/ at the QTL locations by selecting different favourable alleles. In our case, genotypes bearing the favourable alleles from Nemo C, Nemo D, and Nemo G for the QTL on LG 3 and the favourable alleles from Nemo B, Nemo C, and Nemo H for the QTL on LG 4 were selected. In addition, we were careful to represent all six parents in the selected genotypes in order to avoid inbreeding depression. Moreover, the LD in the progeny allows the genome to be covered with a moderate number of markers. Connected populations have been used successfully for the detection of QTL in inbreeding [24,37] and outbreeding species [38], but it is the first time the progeny of a polycross has been used.

4. Materials and Methods

The plants were obtained from Gie GRASS (a French plant breeding company) and come from their elite perennial ryegrass breeding programme. Plant material consisted of six half-sib (HS) families derived from a polycross including six parental plants: Nemo A–D, Nemo G, and Nemo H. The design of the polycross is presented in Figure 1. Four repetitions of lines (10 clones per line) of the parents were used to minimise any effects of spatial proximity. For each parent, all 40 clones were harvested and threshed together. A total of 42 genotypes per HS family were used. Four clones of each genotype were produced by splitting sets of three to five tillers from a mother plant grown in a spaced-plant nursery. Clones were planted in November 2008 in a spaced-plant trial in a randomised block design with four blocks, at Saint-Sauvant (France) (46°23′1′′ N 0°05′ E). Non-surveyed plants of perennial ryegrass were planted around the trial to avoid border effects. The distance between plants was 75 cm in all directions.
Plants were cut at 5 cm height and fertilised with 50 U of nitrogen on 24 March 2009 and again on 16 July 2009. Stretched plant height (PH) was measured, on 18 March (PH0309), then every seven days from 1 April to 28 April (PH0409_1 to PH0409_5) and on 29 September after 11 weeks of regrowth (PH0909). The value of PH was measured using a modified HerboMETRE® (ARVALIS-Institut du Végétal, Paris, France) which made it possible to stretch the leaves and to record plant height from the ground to the top of the longest leaf.
For each genotype, 50 mg of fresh leaves were harvested and DNA was extracted as in [24] (i.e., CTAB (Cetyl Trimethyl Ammonium Bromide) followed by chloroform/octanol (24/1) purification). A total of 283 SSR (Simple Sequence Repeat) and STS (Sequence-Tagged Site) were tested for their polymorphism within the six parents of the polycross. A set of 82 markers evenly spread across the genome were selected to genotype the 252 genotypes (Table 6) as described in [24] (i.e., separation of PCR (Polymerase Chain Reaction) products on a 6.5% acrylamide gel with a LI-COR DNA sequencer 4200 and coding with SAGA Generation 2 software).
To estimate the maximum plant growth rate (PGR) in spring 2009, a non-equilateral hyperbolic function was fitted between PH from 1 April to 28 April (five measurements) and thermal time (base 0 °C), using the NLIN procedure of SAS [46] as in [22] to estimate maximum leaf elongation rate.
For both PH and PGR, an analysis of variance was carried out using the GLM procedure in the SAS software package [46] with the following model: Yij = μ + Gi + Bj + Eij, where μ is the global mean, Gi is the effect of genotype i, Bj is the effect of replication (four replications), and Eij is the error term. The Kolmogorov–Smirnov test in the UNIVARIATE procedure was used to test the normality of the error term Eij in order to detect potential aberrant data. The adjusted mean for each genotype, used for QTL analyses, was computed with LSMEANS in the GLM procedure. The genotype and error variances were estimated with the VARCOMP procedure in the SAS software package. Broad-sense heritabilities (H2) were calculated as H2 = σ2g/(σ2g + σ2e) [5], where σ2g is the genetic variance and σ2e is the error variance. Pearson correlation coefficients were calculated with the Statistica software [47] for pair-wise combinations of mean data for each genotype.
Molecular markers allowed the detection of 32 genotypes for which the male parent was not present in the polycross and 21 plants obtained by self-fertilisation. These plants were discarded in the analyses. Moreover, eight genotypes died in the field. Finally, 191 genotypes were used in the analyses.
For each of the six parents included in the polycross, one genetic map was built with the Joinmap software [48]. The Haldane distance was used. To increase the number of recombination events counted for each parent, the molecular data for Nemo A, Nemo B, and Nemo H from [24] were used to construct the maps of these genotypes and to determine the phase between markers belonging to the same linkage group. For each map, the linkage groups were defined with a Logarithm of Odds (LOD score) higher than 3. For each parent, the phases between markers belonging to the same linkage group were determined with Joinmap. The maps were then combined with the “combine maps” option of Joinmap. The linkage map was drawn using MapChart 2.1 [49].
The QTL were detected with the software MCQTL with the Outbred module [50] as in [24]. The multipopulation connected analysis was used [38] to detect QTL in the six parents within a single model, taking into account additive effects. For each parent, the additive effect was calculated as half the difference between the averages of the two allelic classes.

Acknowledgments

The authors thank the National Association for Research and Technology (ANRT) for funding. We thank the technicians of the experimental station of Jouffray-Drillaud and the GIE Grass for their help and advice on the trial. We thank the employees of the laboratory of biotechnologies of the INRA of Lusignan for helpful advice. We are grateful to Brigitte Mangin and Sylvain Jasson for their help on MCQTL.

Author Contributions

P.B., V.B., S.F. and B.J., conceived and designed the experiment; L.P., J.G., and S.F. carried out the experiments; L.P., S.F., and P.B. analysed the data; L.P., B.J., and P.B. wrote the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hill, W.G. Quantitative genetics in the genomic era. Curr. Genom. 2012, 13, 196–206. [Google Scholar] [CrossRef] [PubMed]
  2. Pérez-de-Castro, A.M.; Vilanova, S.; Cañizares, J.; Pascual, L.; Blanca, J.M.; Diez, M.J.; Prohens, J.; Pico, B. Application of genomic tools in plant breeding. Curr. Genom. 2012, 13, 179–195. [Google Scholar] [CrossRef] [PubMed]
  3. Rafalski, A. Application of single nucleotide polymorphisms in crop genetics. Curr. Opin. Plant Biol. 2002, 54, 94–100. [Google Scholar] [CrossRef]
  4. Smith, S.E. Breeding synthetic cultivars. In Encyclopedia of Plant and Crop Science; Goodman, R.M., Ed.; Rutgers University: New Brunswick, NJ, USA, 2004; pp. 205–206. [Google Scholar]
  5. Gallais, A. Théorie De La Sélection en Amélioration Des Plantes; Masson: Paris, France; Milan, Italy; Barcelon, Spain; Mexico, NM, USA, 1990; pp. 12–26. [Google Scholar]
  6. Flint-Garcia, S.A.; Thornsberry, J.M.; Buckler, E.S. Structure of linkage disequilibrium in plants. Annu. Rev. Plant Biol. 2003, 54, 357–374. [Google Scholar] [CrossRef] [PubMed]
  7. Auzanneau, J.; Huyghe, C.; Escobar-Gutierrez, A.J.; Julier, B.; Gastal, F.; Barre, P. Association study between the gibberellic acid insensitive gene and leaf length in a Lolium perenne L. Synthetic variety. BMC Plant Biol. 2011, 11, 183–196. [Google Scholar] [CrossRef] [PubMed]
  8. Wilkins, P.; Humphreys, M. Progress in breeding perennial forage grasses for temperate agriculture. J. Agric. Sci. 2003, 140, 129–150. [Google Scholar] [CrossRef]
  9. Auzanneau, J.; Huyghe, C.; Julier, B.; Barre, P. Linkage disequilibrium in synthetic varieties of perennial ryegrass. Theor. Appl. Genet. 2007, 115, 837–847. [Google Scholar] [CrossRef] [PubMed]
  10. Brazauskas, G.; Lenk, I.; Pedersen, M.G.; Studer, B.; Lubberstedt, T. Genetic variation, population structure, and linkage disequilibrium in european elite germplasm of perennial ryegrass. Plant Sci. 2011, 181, 412–420. [Google Scholar] [CrossRef] [PubMed]
  11. Rhodes, I. The relationship between productivity and some components of canopy structure in ryegrass (Lolium spp.). II. Yield, canopy structure and ligth interception. J. Agric. Sci. 1971, 77, 283–292. [Google Scholar] [CrossRef]
  12. Rhodes, I. The relationship between productivity and some components of canopy structure in ryegrass (Lolium spp.) III. Spaced plant characters, their heritabilities and relationship to sward yield. J. Agric. Sci. 1973, 80, 171–176. [Google Scholar] [CrossRef]
  13. Rhodes, I.; Mee, S. Changes in dry matter yield associated with selection for canopy characters in ryegrass. Grass Forage Sci. 1980, 35, 35–39. [Google Scholar] [CrossRef]
  14. Horst, G.L.; Nelson, C.J.; Asay, K.H. Relationship of leaf elongation to forage yield of tall fescue genotypes. Crop Sci. 1978, 18, 715–719. [Google Scholar] [CrossRef]
  15. Hazard, L.; Ghesquière, M. Productivity under contrasting cutting regimes of perennial ryegrass selected for short and long leaves. Euphytica 1997, 95, 295–299. [Google Scholar] [CrossRef]
  16. Barre, P.; Emile, J.C.; Betin, M.; Surault, F.; Ghesquière, M.; Hazard, L. Morphological characteristics of perennial ryegrass leaves that influence short-term intake in dairy cows. Agron. J. 2006, 98, 978–985. [Google Scholar] [CrossRef]
  17. Hazard, L.; Ghesquière, M. Evidence from the use of isozyme markers of competition in swards between short-leaved and long-leaved perennial ryegrass. Grass Forage Sci. 1995, 50, 241–248. [Google Scholar] [CrossRef]
  18. Parsons, A.J.; Robson, M.J. Seasonal changes in the physiology of s24 perennial ryegrass (Lolium perenne L.). I. Response of leaf extension to temperature during the transition from vegetative to reproductive growth. Ann. Bot. 1980, 46, 435–444. [Google Scholar]
  19. Jones, E.S.; Mahoney, N.L.; Hayward, M.D.; Armstead, I.P.; Jones, J.; Humphreys, M.O.; King, I.P.; Kishida, T.; Yamada, T.; Balfourier, F.; et al. An enhanced molecular marker based genetic map of perennial ryegrass (Lolium perenne) reveals comparative relationships with other poaceae genome. Genome 2002, 45, 282–295. [Google Scholar] [CrossRef] [PubMed]
  20. Riday, H. Paternity testing: A non-linkage based marker-assisted selection scheme for outbred forage species. Crop Sci. 2011, 51, 631–641. [Google Scholar] [CrossRef]
  21. Gillet, M. Un plan systématique de polycross pour les nombres de base pairs. Ann. Amélior. Plantes 1963, 13, 269–276. [Google Scholar]
  22. Barre, P.; Moreau, L.; Mi, F.; Turner, L.; Gastal, F.; Julier, B.; Ghesquiere, M. Quantitative trait loci for leaf length in perennial ryegrass (Lolium perenne L.). Grass Forage Sci. 2009, 64, 310–321. [Google Scholar] [CrossRef]
  23. Jensen, L.B.; Muylle, H.; Arens, P.; Andersen, C.H.; Holm, P.B.; Ghesquière, M.; Julier, B.; Lübberstedt, T.; Nielsen, K.K.; De Riek, J.; et al. Development and mapping of a public reference set of ssr markers in Lolium perenne L. Mol. Ecol. Notes 2005, 5, 951–957. [Google Scholar] [CrossRef]
  24. Pauly, L.; Flajoulot, S.; Garon, J.; Julier, B.; Beguier, V.; Barre, P. Detection of favorable alleles for plant height and crown rust tolerance in three connected populations of perennial ryegrass (Lolium perenne L.). Theor. Appl. Genet. 2012, 124, 1139–1153. [Google Scholar] [CrossRef] [PubMed]
  25. Faville, M.; Vecchies, A.; Schreiber, M.; Drayton, M.; Hughes, L.; Jones, E.; Guthridge, K.; Smith, K.; Sawbridge, T.; Spangenberg, G.; et al. Functionally associated molecular genetic marker map construction in perennial ryegrass (Lolium perenne L.). Theor. Appl. Genet. 2004, 110, 12–32. [Google Scholar] [CrossRef] [PubMed]
  26. Gill, G.; Wilcox, P.; Whittaker, D.; Winz, R.; Bickerstaff, P.; Echt, C.; Kent, J.; Humphreys, M.; Elborough, K.; Gardner, R. A framework linkage map of perennial ryegrass based on ssr markers. Genome 2006, 49, 354–364. [Google Scholar] [CrossRef] [PubMed]
  27. Studer, B.; Byrne, S.; Nielsen, R.O.; Panitz, F.; Bendixen, C.; Islam, M.S.; Pfeifer, M.; Lubberstedt, T.; Asp, T. A transcriptome map of perennial ryegrass (Lolium perenne L.). BMC Genom. 2012, 13, 140–153. [Google Scholar] [CrossRef] [PubMed]
  28. Studer, B.; Kolliker, R.; Muylle, H.; Asp, T.; Frei, U.; Roldan-Ruiz, I.; Barre, P.; Tomaszewski, C.; Meally, H.; Barth, S.; et al. Est-derived ssr markers used as anchor loci for the construction of a consensus linkage map in ryegrass (Lolium spp.). BMC Plant Biol. 2010, 10, 177–187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Byrne, S.; Czaban, A.; Studer, B.; Panitz, F.; Bendixen, C.; Asp, T. Genome wide allele frequency fingerprints (gwaffs) of populations via genotyping by sequencing. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed]
  30. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 2011, 6. [Google Scholar] [CrossRef] [PubMed]
  31. Kobayashi, S.; Humphreys, M.O.; Tase, K.; Sanada, Y.; Yamada, T. Molecular marker dissection of ryegrass plant development and its response to growth environments and foliage cuts. Crop Sci. 2011, 51, 600–611. [Google Scholar] [CrossRef]
  32. Sartie, A.M.; Matthew, C.; Easton, H.S.; Faville, M.J. Phenotypic and qtl analyses of herbage production-related traits in perennial ryegrass (Lolium perenne L.). Euphytica 2011, 182, 295–315. [Google Scholar] [CrossRef]
  33. Yamada, T.; Jones, E.S.; Cogan, N.O.I.; Vecchies, A.C.; Nomura, T.; Hisano, H.; Shimamoto, Y.; Smith, K.F.; Hayward, M.D.; Forster, J.W. QTL analysis of morphological, developmental, and winter hardiness-associated traits in perennial ryegrass. Crop Sci. 2004, 44, 925–935. [Google Scholar] [CrossRef]
  34. Turner, L.; Cairns, A.; Armstead, I.; Thomas, H.; Humphreys, M.W.; Humphreys, M.O. Does fructan have a functional role in physiological traits? Investigation by quantitative trait locus mapping. New Phytol. 2008, 179, 765–777. [Google Scholar] [CrossRef] [PubMed]
  35. Studer, B.; Jensen, L.B.; Hentrup, S.; Brazauskas, G.; Kölliker, R.; Lübberstedt, T. Genetic characterisation of seed yield and fertility traits in perennial ryegrass (Lolium perenne L.). Theor. Appl. Genet. 2008, 117, 781–791. [Google Scholar] [CrossRef] [PubMed]
  36. Muranty, H. Power of tests for quantitative trait loci detection using full-sib families in different schemes. Heredity 1996, 76, 156–165. [Google Scholar] [CrossRef]
  37. Blanc, G.; Charcosset, A.; Mangin, B.; Gallais, A.; Moreau, L. Connected populations for detecting quantitative trait loci and testing for epistasis: An application in maize. Theor. Appl. Genet. 2006, 113, 206–224. [Google Scholar] [CrossRef] [PubMed]
  38. Billotte, N.; Jourjon, M.F.; Marseillac, N.; Berger, A.; Flori, A.; Asmady, H.; Adon, B.; Singh, R.; Nouy, B.; Potier, F.; et al. QTL detection by multi-parent linkage mapping in oil palm (Elaeis guineensis Jacq.). Theor. Appl. Genet. 2010, 120, 1673–1687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Hirata, M.; Cai, H.W.; Inoue, M.; Yuyama, N.; Miura, Y.; Komatsu, T.; Takamizo, T.; Fujimori, M. Development of simple sequence repeat (SSR) markers and construction of an ssr-based linkage map in italian ryegrass (Lolium multiflorum Lam.). Theor. Appl. Genet. 2006, 113, 270–279. [Google Scholar] [CrossRef] [PubMed]
  40. Lauvergeat, V.; Barre, P.; Bonnet, M.; Ghesquière, M. Sixty simple sequence repeat markers for use in the festuca—Lolium complex of grasses. Mol. Ecol. Notes 2005, 5, 401–405. [Google Scholar] [CrossRef]
  41. Studer, B.; Asp, T.; Frei, U.; Hentrup, S.; Meally, H.; Guillard, A.; Barth, S.; Muylle, H.; Roldan-Ruiz, I.; Barre, P.; et al. Expressed sequence tag-derived microsatellite markers of perennial ryegrass (Lolium perenne L.). Mol. Breed. 2008, 21, 533–548. [Google Scholar] [CrossRef]
  42. Jones, E.S.; Dupal, M.P.; Kölliker, R.; Drayton, M.C.; Forster, J.W. Development and characterisation of simple sequence repeat (SSR) markers for perennial ryegrass (Lolium perenne L.). Theor. Appl. Genet. 2001, 102, 405–415. [Google Scholar] [CrossRef]
  43. Kubik, C.; Sawkins, M.; Meyer, W.A.; Gaut, B.S. Genetic diversity in seven perennial ryegrass (Lolium perenne L.) cultivars based on ssr markers. Crop Sci. 2001, 41, 1565–1572. [Google Scholar] [CrossRef]
  44. Saha, M.C.; Rouf Mian, M.A.; Eujayl, I.; Zwonitzer, J.C.; Wang, L.; May, G.D. Tall fescue est-ssr markers with transferability across several grass species. Theor. Appl. Genet. 2004, 109, 783–791. [Google Scholar] [CrossRef] [PubMed]
  45. Lem, P.; Lallemand, J. Grass consensus sts markers: An efficient approach for detecting polymorphism in Lolium. Theor. Appl. Genet. 2003, 107, 1113–1122. [Google Scholar] [CrossRef] [PubMed]
  46. SAS, Version 9.3; Sas Institute Inc.: Cary, NC, USA, 2000.
  47. StatSoft. Statistica Pour Windows. Available online: Http://www.Statsoft.Fr (accessed on 20 September 2016).
  48. Van Ooijen, J.W.; Voorrips, R.E. Joinmap 3.0 Software for the Calculation of Genetic Linkage Maps; Plant Research International: Wageningen, The Netherlands, 2001; pp. 1–51. [Google Scholar]
  49. Voorrips, R.E. Mapchart: Software for the graphical presentation of linkage maps and qtls. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed]
  50. Jourjon, M.F.; Jasson, S.; Marcel, J.; Ngom, B.; Mangin, B. Mcqtl: Multi-allelic qtl mapping in multi-cross design. Bioinformatics 2005, 21, 128–130. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Design of the polycross including Nemo A–D, Nemo G, and Nemo H organised in four repetitions. Nemo E and Nemo F were discarded before crossing. In each repetition, each parent is represented by a line of 10 clones (dots). The polycross was planted in a field of rye (in yellow).
Figure 1. Design of the polycross including Nemo A–D, Nemo G, and Nemo H organised in four repetitions. Nemo E and Nemo F were discarded before crossing. In each repetition, each parent is represented by a line of 10 clones (dots). The polycross was planted in a field of rye (in yellow).
Agronomy 06 00051 g001
Figure 2. Example of plant regrowth in spring for the genotype 23 (four replicates). A non-equilateral hyperbolic function was fitted between stretched plant height (PH) and thermal time (mm/°Cd, base 0 °C).
Figure 2. Example of plant regrowth in spring for the genotype 23 (four replicates). A non-equilateral hyperbolic function was fitted between stretched plant height (PH) and thermal time (mm/°Cd, base 0 °C).
Agronomy 06 00051 g002
Figure 3. Consensus map of the six parents of the polycross. Seven linkage groups (LG) are presented.
Figure 3. Consensus map of the six parents of the polycross. Seven linkage groups (LG) are presented.
Agronomy 06 00051 g003
Figure 4. Breeding scheme including marker assisted selection (MAS) with QTL identification in the progenies of polycrosses.
Figure 4. Breeding scheme including marker assisted selection (MAS) with QTL identification in the progenies of polycrosses.
Agronomy 06 00051 g004
Table 1. Number of genotypes per family and heading dates of the parents.
Table 1. Number of genotypes per family and heading dates of the parents.
Nemo ANemo BNemo CNemo DNemo GTotal Number of Genotypes per HS Family Heading Dates 1
Nemo A 84157
Nemo B20 55157
Nemo C53 48153
Nemo D27510 67156
Nemo G2322216 70157
Nemo H95289758154
1 Calendar day counting from 1 January; HS: half-sib.
Table 2. Genetic variance, error variance, and broad sense heritability (H2 = σ2g/(σ2g + σ2e)) for stretched plant heights (PH in mm) in spring (0309 and 0409) and autumn (0909) and maximum plant growth rate (PGR in mm/°Cd, base 0 °C) in spring.
Table 2. Genetic variance, error variance, and broad sense heritability (H2 = σ2g/(σ2g + σ2e)) for stretched plant heights (PH in mm) in spring (0309 and 0409) and autumn (0909) and maximum plant growth rate (PGR in mm/°Cd, base 0 °C) in spring.
Traitsσ2gσ2eH2
PH03094143820.52
PH0409_11721640.51
PH0409_24514510.50
PH0409_33783600.51
PH0409_45125520.48
PH0409_55956650.47
PGR0.0390.0830.32
PH09094263940.52
Table 3. Distribution of the mean per genotype for stretched plant heights (PH in mm) and maximum plant growth rate (PGR in mm/ Cd, base 0 °C).
Table 3. Distribution of the mean per genotype for stretched plant heights (PH in mm) and maximum plant growth rate (PGR in mm/ Cd, base 0 °C).
TraitsAverageMinimumMaximumStandard Deviation
PH03091297119022
PH0409_11146614715
PH0409_21345318524
PH0409_31379623722
PH0409_419811426625
PH0409_520512328228
PGR0.950.481.660.25
PH090920311627124
Table 4. Pearson correlation coefficients for pair-wise combinations of mean data for each genotype for stretched plant heights (PH in mm) and maximum plant growth rate (PGR in mm/°Cd, base 0 °C).
Table 4. Pearson correlation coefficients for pair-wise combinations of mean data for each genotype for stretched plant heights (PH in mm) and maximum plant growth rate (PGR in mm/°Cd, base 0 °C).
TraitsPH0309PH0409_1PH0409_2PH0409_3PH0409_4PH0409_5PGR
PH0409_10.83 ***
PH0409_20.89 ***0.90 ***
PH0409_30.83 ***0.83 ***0.88 ***
PH0409_40.81 ***0.83 ***0.87 ***0.88 ***
PH0409_50.77 ***0.82 ***0.83 ***0.87 ***0.88 ***
PGR0.81 ***0.86 ***0.90 ***0.78 ***0.73 ***0.68 ***
PH09090.18 *0.11 ns0.15 *0.13 ns0.20 **0.17 *0.13 ns
*** significant at 0.001; ** significant at 0.01; * significant at 0.05; ns: not significant at 0.05.
Table 5. Quantitative traits (QTL) detected by the multipopulation connected analysis with additive effects. The software MCQTL with the Outbred module was used. Significant effects are indicated in bold.
Table 5. Quantitative traits (QTL) detected by the multipopulation connected analysis with additive effects. The software MCQTL with the Outbred module was used. Significant effects are indicated in bold.
Additive Effects
TraitsLGMarkersMap Location (cM)Confidence Intervalr2 (%)Global r2 (%)Nemo ANemo BNemo CNemo DNemo GNemo H
PH03093pps000729.328.1–41.212.028.1−1.8 14.55.24.6−9.51.9
4LpSSR01130.826.4–35.920.1 −3.87.410.7−1.23.410.8
PH0409_13pps000729.327.0–30.614.028.4−1.70.95.84.8−5.51.0
4G01-07328.827.9–34.419.2 −3.34.15.70.14.06.0
PH0409_23pps000729.327.2–30.615.627.7−1.23.08.88.3−9.73.1
4G01-07527.826.2–32.617.3 −4.64.59.61.55.610.6
PH0409_34pps004028.225.8–35.117.3 −6.64.37.41.26.78.7
PH0409_44LpSSR01131.726.2–35.016.1 −5.311.87.80.84.09.1
PH0409_54LpSSR01130.827.5–33.022.1 −6.612.411.7−2.68.411.6
PGR4G05-01433.819.2–36.212.3 −0.050.080.11−0.040.020.06
1 Within each parent, additive effects having the same sign are in phase.
Table 6. Marker description: names, references, and linkage group (LG).
Table 6. Marker description: names, references, and linkage group (LG).
Marker NamesReferenceLGMarker NamesReferenceLG
25ca1[26]3LmgSSR01-08H[39]2
B1A2[40]3LmgSSR02-10B 6
B1C8 7LmgSSR03-04E 3
B1C9 3LmgSSR04-09F 1
B3B8 3LmgSSR10-12D 7
B3-C5 7LmgSSR15-08E 4
B4D7op 1LmgSSR15-09F 7
DLF008[23]7LpSSR006[23]4
DLF027 1LpSSR011 4
G01-001[41]7LpSSR020 7
G01-037 4LpSSR023 4
G01-040 2LpSSR058 6
G01-043 2LpSSR066 7
G01-047 6LpSSR085 1
G01-053 3LpSSR112 2
G01-054 2LpSSRH01A02[42]5
G01-073 4M15-185[43]2
G01-075 4M4-136 2
G01-079 1M4-213 1
G01-080 5NFFA030[44]2
G01-090 7NFFA036 6
G01-095 5NFFA064 5
G02-004 7NFFA087 5
G02-047 1NFFA099 3
G02-049 2NFFA114 7
G02-057 5OSW[45]7
G03-003 6pps0007[25]3
G03-079 6pps0040 4
G04-034 4pps0049 7
G05-014 4pps0080 2
G05-028 3pps0164 3
G05-033 5pps0299 6
G05-046 6pps0397 5
G05-050 2rv0244[26]1
G05-065 5rv0641 6
G05-073 7rv0706 2
G05-090 3rv0757 5
G05-092 5rye012[23]4
G05-129 2S7F7-3-4F[39]6
G07-034 3syn20738-1[22]5
LmgSSR01-02H[39]6uni001[23]3

Share and Cite

MDPI and ACS Style

Pauly, L.; Flajoulot, S.; Garon, J.; Julier, B.; Béguier, V.; Barre, P. Quantitative Trait Loci (QTL) Identification in the Progeny of a Polycross. Agronomy 2016, 6, 51. https://doi.org/10.3390/agronomy6040051

AMA Style

Pauly L, Flajoulot S, Garon J, Julier B, Béguier V, Barre P. Quantitative Trait Loci (QTL) Identification in the Progeny of a Polycross. Agronomy. 2016; 6(4):51. https://doi.org/10.3390/agronomy6040051

Chicago/Turabian Style

Pauly, Laurence, Sandrine Flajoulot, Jérôme Garon, Bernadette Julier, Vincent Béguier, and Philippe Barre. 2016. "Quantitative Trait Loci (QTL) Identification in the Progeny of a Polycross" Agronomy 6, no. 4: 51. https://doi.org/10.3390/agronomy6040051

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