Skip to main content

Cytogenetic and molecular genotyping in the allotetraploid Festuca pratensis × Lolium perenne hybrids



Species of the Festuca and Lolium genera, as well as intergeneric Festuca × Lolium (Festulolium) hybrids, are valuable fodder and turf grasses for agricultural and amenity purposes worldwide. Festulolium hybrids can merge in their genomes agronomically important characteristics. However, in polyploid plants, especially in allopolyploids, the hybridization of divergent genomes could contribute to various abnormalities, such as variability in chromosome number, structural rearrangements, and/or disorders in inheritance patterns. Here we studied these issues in allotetraploid Festuca pratensis × Lolium perenne hybrids.


Cytogenetic procedures, including fluorescent in situ hybridization, genomic in situ hybridization, and molecular markers – inter-simple sequence repeats (ISSR) were exploited. This cytogenetic approach indicated the dynamics in the number and distribution of ribosomal RNA genes and structural rearrangements for both parental genomes (Festuca and Lolium) in hybrid karyotypes. The separate analysis of F. pratensis and L. perenne chromosomes in hybrid plants (F2-F3 generations of F. pratensis × L. perenne) revealed the asymmetrical level of rearrangements. Recognized structural changes were mainly located in the distal part of chromosome arms, and in chromosomes bearing ribosomal DNA, they were more frequently mapped in arms without this sequence. Based on the ISSR markers distribution, we found that the tetrasomic type of inheritance was characteristic for the majority of ISSR loci, but the disomic type was also observed. Nonetheless, no preference in the transmission of either Festuca or Lolium alleles to the following generations of allotetraploid F. pratensis × L. perenne hybrid was observed.


Our study reports cytogenetic and molecular genotyping of the F. pratensis × L. perenne hybrid and its following F2-F3 progenies. The analysis of 137 allotetraploid F. pratensis × L. perenne hybrids revealed the higher level of recombination in chromosomes derived from F. pratensis genome. The results of ISSR markers indicated a mixed model of inheritance, which may be characteristic for these hybrids.


Polyploids frequently occur within the plant kingdom, especially in domesticated and economically important plants in the world [1]. Within polyploids, two basic types auto- and allopolyploids are usually. Genomes of allopolyploids contain chromosome sets from two or more divergent genomes of donor plants. In the course of evolution, plant genomes underwent various processes that had an impact on numerical and structural changes [2, 3]. The comparison of karyotypes within the grass family and various genera revealed vast variability in the chromosome number and in their morphology. There are three different levels of changes that may occur in the hybrid genomes: (i) DNA sequence (e.g., allelic differences caused by nucleotide substitutions or indels), (ii) structure and gene order (e.g., chromosomal rearrangements), and (iii) total number of chromosomes in the complement. In allopolyploids, especially in the newly formed, smaller (at DNA sequence level) and larger (at chromosome level) genetic changes could result from recombination between homoeologous chromosomes. Pairing affinity and chiasmata formation between homoeologous chromosome sets have an enormous impact on inheritance patterns [4,5,6]. Importantly, several gametophytic and zygotic barriers have been reported, which caused deviation of allele frequencies from Mendelian ratios, for example in the interspecific crosses in rice [7].

Variation in a number of chromosomes and/or in their structure are pivotal issues in the crossing of plants derived from different genera. Plants with the expected chromosome number, as well as aneuploids may be observed within the hybrid’s progeny [8,9,10]. Recombinant chromosomes can be characterized using molecular and cytogenetic tools, such as genetic markers (e.g. [11, 12]) and genomic in situ hybridization (GISH) technique (e.g. [9, 13,14,15,16,17]). The application of GISH, usually combined with fluorescent in situ hybridization (FISH), enables the characterization of the genome composition of hybrid plants providing valuable information for breeding new cultivars on the early stages of material selection.

Lolium species (L. perenne and L. multiflorum), the major source of forage in temperate regions of the world, depict the considerable variation of quality traits, crucial implication for grass breeding programs. This group of plants is characterized by high digestibility, palatability and intensive spring growth [18]. Traits which are of great importance for adaptation to abiotic as well as biotic stresses, e.g. drought tolerance, winter hardiness and resistance to fungal diseases, are provided by Festuca species (F. pratensis and F. arundinacea) [18]. Desirable traits from Lolium species can be merged with features of Festuca species through crosses and production of intergeneric Festuca × Lolium hybrids (Festulolium). Chromosomes of these hybrids can recombine at high frequency [9, 11, 19].

The main aims of this study were the examination of the genomic composition of F. pratensis × L. perenne hybrid and its following generations (F2-F3) and the determination of inheritance pattern in this plant material. To achieve these goals, we applied the cytogenetic procedures, including FISH with ribosomal DNA (rDNA) sequences to recognize selected chromosomes in the complements, GISH to analyze the genomic architecture, and molecular markers – inter-simple sequence repeats (ISSR) to determine the types of ISSR loci inheritance in three successive generations of F. pratensis × L. perenne.


Karyotype analysis of F. pratensis × L. perenne hybrids

In total, 134 plants of F2 (combinations no. 1 and 2) and F3 (combinations no. 3 and 4) generations were analyzed (see Additional file 1: Figure S1) and the majority of them were tetraploid with 2n = 4x = 28 (83.3% in F2 and 68.2% in F3). Among aneuploid plants, both hypoploids (2n = 26, 27) and hyperploids (2n = 29, 30, 31) were recorded. Hypoploids constituted 11.1% in F2 and 25% in F3, while hyperploids – 5.6% and 6.8%, respectively (Figs. 1, 2).

Fig. 1

The distribution of the total number of chromosomes in karyotypes and rDNA-bearing chromosomes in 90 plants of F2 generation of F. pratensis × L. perenne hybrid – in combination no. 1 (a) and in combination no. 2 (b). At the top of both graphs were shown results for parents – two parental plants: P (Fp 11/59) (♀), P (Lp 08/33) (♂); four F1 plants: two plants – F1 6-2A (♀), F1 6-7A (♂) for combination no. 1 of F2, and two plants – F1 6-2A (♀), F1 6-3B (♂) for combination no. 2 of F2. On the y-axis obtained rDNA patterns (Roman numerals) and their frequency in the brackets (Arabic numerals) were shown. On x-axis were presented: a total number of chromosomes divided into two groups – Festuca and Lolium chromosomes, as well as the number of particular rDNA-bearing chromosomes

Fig. 2

The distribution of the total number of chromosomes in karyotypes and rDNA-bearing chromosomes in 44 plants of F3 generation of F. pratensis × L. perenne hybrid – in combination no. 3 (a), and in combination no. 4 (b). At the top of both graphs were shown results for parents – two F1 plants: F1 6-2A (♀), F1 6-7A (♂) and four plants of F2: two plants – F2 2A-15 (♀), F2 2A-5 (♂) for combination no. 3 of F3, and two plants – F2 2A-27 (♀), F2 2A-12 (♂) for combination no. 4 of F3. On the y-axis obtained rDNA patterns (Roman numerals) and their frequency in the brackets (Arabic numerals) were shown. On x-axis were presented: a total number of chromosomes, divided into two groups - Festuca and Lolium chromosomes, as well as the number of particular rDNA-bearing chromosomes

Allotetraploid plants with 14 chromosomes of F. pratensis and 14 chromosomes of L. perenne were distinguished in 56% of F2 and 40% of F3 generation. Within the remaining plants with 28 chromosomes, the genotypes with the higher number of L. perenne chromosomes were predominated (29.3% in F2 and 53.3% in F3). Plants with lower (2n < 28) and higher (2n > 28) total number of chromosomes, that constituted 16.7% in F2 and 31.8% in F3, were also characterized by the higher number of L. perenne chromosomes (average 43.3% for both generations) (see Additional file 2: Table S1).

rDNA-FISH mapping

rDNA-FISH mapping revealed that F. pratensis and L. perenne species, used to obtain F1 hybrids, had four loci of 5S rDNA. In the case of 35S rDNA, four and 15 loci were recorded, respectively. Among F1 plants, 5S rDNA sequence was mapped to four chromosomes in the genomes of maternal plant and both pollinators. Whereas 35S rDNA was hybridized to eight sites in the karyotype of pollinator 6-7A and nine sites in the female plant and pollinator 6-3B (Fig. 1). Four F2 plants used to produce two combinations of F3 (no. 3 and 4) were characterized by four loci of 5S rDNA and eight loci of 35S rDNA (combination no. 3), as well as four loci of 5S rDNA and nine loci of 35S rDNA (combination no. 4) (Fig. 2). Among F2 progeny, the number of 5S rDNA loci was equal to four, except for one plant with five loci of this sequence (Fig. 1b). This sequence was located in an interstitial part of chromosome 3F and near the centromeric region of chromosome 3 L. In the vast majority of F2 individuals, there were two pairs of each chromosome types, in minority three chromosomes 3F and only one chromosome 3L (see Additional file 2: Table S1). The similar trend was observed among F3 plants – the number of 5S rDNA equaled four, with the exception of one plant with three loci of 5S rDNA (Fig. 2b). In all F3 analyzed hybrids, two chromosomes 3F and two 3L were mostly observed. Only for several plants, three 3F and one 3L chromosomes as well as conversely one 3F and three 3L chromosomes were recorded (Fig. 2). In the karyotype of plants with three loci of 5S rDNA, one chromosome 3F and two chromosomes 3L were observed, while for the karyotype with five loci of this sequence three 3F and two 3L chromosomes were noticed.

The number of 35S rDNA sequence was more variable and ranged from seven to eleven loci in F2 and F3 plants (Figs. 1, 2). This sequence was located at the secondary constriction regions of 2F and 3L chromosomes, as well as on the chromosomes 1, 2 and 7L. In both generations, the variation in the number of particular chromosomes was as follows: 3-7 signals for chromosomes 1, 2, 7L, 1-3 for chromosomes 3L, and 0-3 for chromosomes 2F (Figs. 1, 2).

With regard to the type and the number of recognized rDNA-bearing chromosome pairs, the analysis showed 22 and 25 various rDNA patterns in combination no. 1 and 2 of F2 generation, respectively (Fig. 1). Whereas in F3 generation for combinations no. 3 and no. 4, subsequently 21 and 14 rDNA patterns were detected (Fig. 2). The most frequent patterns (48.9%) in F2 combination no. 1 reflected rDNA arrangements occurring in parental plants. In the remaining F2 and F3 plants, different rDNA profiles were identified. However, the most repeated patterns were connected with one more or one less 35S rDNA locus. The analysis of the total number of rDNA-bearing chromosomes for all F2 and F3 analyzed plants revealed that even number of chromosomes prevailed, what was expected, however, with the exception of the number of chromosomes 2F. In F3 combination no. 3, almost 70% of plants had in karyotype only one chromosome 2F (Fig. 2).

GISH mapping

To reveal rearrangements in the karyotypes of F. pratensis × L. perenne hybrids, the analysis with total genomic DNA (GISH) was performed and this approach did not reveal any structural changes in F1 plants. The F2 and F3 progenies of the F. pratensis × L. perenne hybrid were characterized by karyotype rearrangements. Their number increased from generation to generation and was as follows: 179 for F2 plants and 217 for F3 plants (see Additional file 2: Table S1). These numbers include the terminal, interstitial and pericentromeric types of rearrangements in total. The rearranged chromosomes were divided into two groups consisting of rearranged rDNA-bearing chromosomes and rearranged non-rDNA-bearing ones. The level of rearrangements in rDNA-bearing chromosomes equaled 38.8% in F2 and 47% in F3 generation, while in non-rDNA-bearing chromosomes 61.2% and 53%, respectively. Among identified rearrangements, terminal (92.2%), interstitial (7.1%) and pericentromeric (0.7%) types were recognized. Non-rDNA-bearing chromosomes were more frequent (61.4% in F2 and 63.2% in F3) in the group of non-recombined chromosomes (see Additional file 2: Table S1).

The separate analysis of rearranged chromosomes derived from F. pratensis and L. perenne genomes revealed that the number of recombined rDNA-bearing chromosomes was on a similar level for F2 plants, 51.7% and 48.3% for F. pratensis and L. perenne genome, respectively. More differences were found in F3 generation, 40.4% for F. pratensis and 59.6% for L. perenne. When both generations were analyzed together, slightly predominance of the rearranged rDNA-bearing chromosomes of L. perenne (55.5%) was observed. On the other hand, in non-rDNA-bearing chromosomes, the higher level of recombination was detected for F. pratensis chromosomes (67.2% in F2 and 76.7% in F3). Additionally, among analyzed plants with detected structural rearrangements, several genotypes with the uncommon location of changes were noticed. In these plants, the rearrangements colocalized with the position of rDNA locus (5S and 35S rDNA) in various type of rDNA-bearing chromosomes – chromosome 3L/non-rDNA-bearing Festuca chromosome (Fig. 3a and b), chromosome 2F (Fig. 3c and d), and chromosomes 1, 2, 7L (Fig. 3e and f). The similar position of rearrangement was also detected in chromosome 3F. Furthermore, an F2 plant without two 2F chromosomes was observed (Fig. 3g and h).

Fig. 3

Rearrangements in rDNA-bearing chromosomes among F. pratensis × L. perenne hybrids. (a, b) F2 2A-33 with recombined 3L/non-rDNA-bearing Festuca chromosome; (c, d) F2 2B-21 with recombined 2F chromosome; (e, f) F2 2B-28 with recombined one of 1, 2, 7L group of chromosomes, and (g, h) F2 2A-5 without 2F chromosomes. (a, c, e, g) FISH with 5S rDNA (red) and 35S rDNA (green); chromosomes were counterstained with DAPI (blue). (b, d, f, h) GISH with total genomic DNA of L. perenne as a probe (green) and with blocking genomic DNA of F. pratensis (orange) on the same metaphase plate; chromosomes were counterstained with propidium iodide (orange). Red dotted lines indicate the place of rearrangement. Chromosomes with rearrangements corresponding with rDNA position were exemplified in the insets. Scale bar 5 μm

The detailed analysis of recombined rDNA-bearing chromosomes in F2 and F3 plants unveiled the types and the frequency of observed changes. Rearranged regions were mostly located in the arms without rDNA-marker for both parental genomes and both generations and they always were located in the distal parts of chromosomes (Fig. 4). In F2 plants, the most recombined types were: chromosomes 2F and the group of chromosomes 1, 2, 7L, while in F3 – the group of chromosomes 1, 2, 7L and chromosomes 3F (see Additional file 2: Table S1). However, the types of rearrangements and their frequencies in the group of chromosomes 1, 2, 7L were not shown in Fig. 4, because the location of 35S rDNA and similar size of chromosomes do not allow to distinguish them, therefore we could not determine precisely the type and the frequency of changes in each pair of chromosomes separately. In the F2 generation, the highest level of recombination was recorded for 35S rDNA-bearing arm of chromosomes 3L. In the F3 generation, however, more rearrangements were found in the opposite location – in arm bearing 5S rDNA sequence. For chromosomes 2F and 3F, the majority of changes in the arms without rDNA landmarks were occurred (Fig. 4).

Fig. 4

The distribution and the frequency of rearrangements in rDNA-bearing chromosomes (3L, 2F, and 3F) among F2 and F3 generations of F. pratensis × L. perenne hybrid. Green dots on the chromosomes illustrate 35S rDNA loci; pink dots illustrate 5S rDNA; orange blocks illustrate the position and size of rearrangements

Statistical analysis of cytogenetic traits

Statistical analysis of selected cytogenetic traits in F2 generation revealed weakly significant differences between distributions and strongly significant differences between mean values for the number of 35S rDNA in L. perenne genome what is related to the number of chromosomes 1, 2, 7L (Table 1). The analysis of F2 plants and the comparison between F2 and F3 generations revealed the statistically significant differences for the number of non-rDNA-bearing chromosomes (Tables 1 and 2). In comparison between generations, it was also revealed that the statistically significant difference existed for the number of 35S rDNA in F. pratensis genome what is related to the number of chromosomes 2F (Table 2).

Table 1 Comparison of nine selected cytogenetic trait distributions in F2 generation of F. pratensis × L. perenne hybrid
Table 2 Comparison of nine selected cytogenetic trait distributions in F2 and F3 generations of F. pratensis × L. perenne hybrid

We also studied separately the chromosomes derived from F. pratensis and L. perenne genomes, divided into groups: rearranged and non-rearranged, as well as rDNA-bearing and non-rDNA-bearing chromosomes. Rearranged chromosomes with rDNA sequence, were subdivided into two classes based on the location of rearrangements – in the arm with rDNA and without rDNA landmark (Table 3). Statistical differences were observed among F2 and F3 generations for non-recombined rDNA-bearing and non-rDNA-bearing chromosomes for both F. pratensis and L. perenne, as well as for recombined non-rDNA-bearing chromosomes of both genomes, but only for F3 generation. In this analysis, we also considered the comparison between analyzed generations. Statistically significant differences were identified for non-recombined rDNA-bearing chromosomes of both genomes and non-recombined non-rDNA-bearing chromosomes of F. pratensis genome, as well as for the number of recombined rDNA- and non-rDNA-bearing chromosomes for both genomes (Table 3). Additionally, the analysis of rearrangement locations revealed statistically significant differences in arms without rDNA markers (Table 3).

Table 3 Comparison of distributions of selected cytogenetic traits (recombined and non-recombined chromosomes belonging to F. pratensis and L. perenne genomes) within and between F2 and F3 generations of F. pratensis × L. perenne hybrid

ISSR loci polymorphism between parental plants

The ISSR analysis resulted in a reliable amplification of 71 loci in parental plants (F. pratensis and L. perenne). The total number of loci amplified in subsequent generations was lower due to the occurrence of recessive genotypes in both F1 and F2 parental individuals. Each ISSR primer produced markers with similar informativeness, as calculated by Polymorphism Information Content (PIC) for each of the parental, F1 and F2 individuals, as well as for each of F2 and F3 cross combinations (see Additional file 3: Table S2). The highest level of polymorphism, accounted for nearly 72% of loci, was observed between F. pratensis and L. perenne parental forms. The frequency of polymorphic loci in F1 generation was lower than in parental forms, although the individuals used for combination no. 2 were more polymorphic compared with the parental individuals of combination no. 1. Similarly, the polymorphism between individuals from F2 generation was further reduced, but F2 parents used for combination no. 4 showed more diversity than those used for combination no. 3 (Table 4). This finding was also reflected by the mean values of PIC, which showed that ISSR markers detected in the parents of cross combination 2 and 4 were more informative, and as a result higher mean value of PIC was also found for the progeny of these cross combinations (see Additional file 3: Table S2).

Table 4 ISSR polymorphism between parental species – F. pratensis and L. perenne and between pairs of the F1 and F2 plants used as parents for subsequent generations of F. pratensis × L. perenne hybrid

Segregation of ISSR loci in F2 and F3 generations

More than one type of the segregation ratio passed the χ2 test in the majority of ISSR loci (43 loci) in the combination no. 1 of F2 generation (see Additional file 4: Table S3). For 25 loci only one, non-distorted segregation ratio was noticed and in the case of two loci a segregation distortion was found. Additionally, one locus passed the χ213:15 test, indicating its tetrasomic inheritance and the possibility of quadrivalents formation and random chromatid segregation. In the F2 combination no. 2, different proportion of each of those categories was recorded: for 33 loci only one non-distorted segregation type was found, in case of 27 loci more possibilities were detected and 11 ISSR loci showed distorted segregation (see Additional file 5: Table S4). Regardless of the above differences, similar proportions of particular inheritance types of ISSR loci were found in both F2 combinations. The tetrasomic inheritance was significantly prevailing over the disomic inheritance. Furthermore, the χ2 test results suggested that three loci (L7_300, L8_300, and L8_1400) have been inherited according to a different model in both F2 combinations (see Additional file 6: Table S5). Importantly, it was not possible to reliably assign the inheritance type to the majority of loci what may somehow bias these calculations (Table 5, see Additional file 6: Table S5). In these loci where distorted segregation was detected, the comparable proportion was skewed toward a dominant or recessive allele. In half of these cases, it was possible to point to the prevalence of the allele from specific parental species, and again similar proportion of distorted loci were skewed toward Festuca or Lolium allele, respectively (see Additional file 6: Table S5).

Table 5 Number of ISSR loci that showed particular inheritance type in F2 and F3 generations of F. pratensis × L. perenne hybrid

The similar analysis of F3 generation gave much less reliable results. This was caused by a significantly lower ISSR polymorphism between F2 plants taken as parents for the crosses. In the majority of loci, both lines were characterized by the presence of ISSR band making the genotype assignment more difficult. In addition, in 12 and 11 loci no bands amplification was observed in combination no. 3 and 4 of F3, respectively. The great majority of loci showed more than one possible segregation ratio in both F3 combinations (49 and 45 loci, for combinations no. 3 and 4, respectively). Three loci had distorted segregation in the combination no. 3 – one locus was skewed toward the recessive allele and two toward the dominant one. Five loci were distorted in the combination no. 4 – here four were skewed toward recessive allele and only one toward the dominant allele. In case of two loci, it was possible to point to the prevalence of the allele from specific parental species and in both of them, the distortion was toward Lolium allele (see Additional file 7: Table S6 and Addititional file 8: Table S7). When the occurrence of distorted loci was compared between F2 and F3 generation no specific pattern was noticed. Only one locus with distorted segregation in F2 generation was further distorted in the F3. In the remaining loci, no co-occurrence of such events in both generations was found. Moreover, the ISSR analysis of F3 generation did not allow to discover the inheritance type for the vast majority of loci (Table 5, see Additional file 9: Table S8).

In addition to the above analysis, the distribution of the frequencies of the most probable genotypes that was inferred based on ISSR loci segregation was provided (Fig. 5). It shows the gradual prevalence of dominant allele frequency in subsequent generations. Such observation may further support the tetrasomic type of ISSR inheritance, as in this type the dominant class is usually much numerous than the recessive class for the same cross combination of parental genotypes. It may also reflect the process of more frequent transmission of a dominant allele from generation to generation.

Fig. 5

The distribution of the most probable genotypes that was inferred based on the segregation of ISSR loci for parental genotypes, F1 parents of both cross combinations and the F2 parents of the cross combination no. 1. In case of ISSR loci segregation in F2 generation from cross combination no. 2 and both F3 generations, only the frequency of precedence (dominant) or absence (recessive, nulliplex) of ISSR band was taken into account


The formation of new synthetic allopolyploid plants via hybridization of different species can include: chromosomal rearrangements [17, 20], transposon activations [21,22,23], changes in rDNA loci [24, 25] and/or loss of DNA [26,27,28]. In allopolyploids, like allotetraploid F. pratensis × L. perenne hybrids, in which two genomes/species from different genera are merged, the deciphering of those processes is crucial to produce more persistent high-quality grass cultivars. For that purpose, the comprehensive analysis, including the examination of genome composition and structural reshufflings, the segregation of ISSR loci and their inheritance patterns in F2-F3 generations of allotetraploid F. pratensis × L. perenne hybrid was performed.

Genomic structure of hybrids

In rDNA-FISH experiments, the same number and distribution of 5S rDNA were observed in the vast majority of studied F2 and F3 generations (in 131 out of 134 plants). Only in three plants, the lower or higher number of 5S rDNA loci was found. Nonetheless, this phenomenon was related to the changed total number of chromosomes in the complement. The more considerable variation was noticed for the number of 35S rDNA that ranged from seven to eleven loci in F2 and F3 plants, and the statistical analysis of this trait revealed significant differences between both generations, especially for the number of chromosomes 1, 2, 7L and 2F. It has been observed for many different plant species and hybrids that the number of 5S rDNA is stable, whereas 35S rDNA sequence seems to be more polymorphic [29, 30], what reflects showed rDNA patterns among studied F. pratensis × L. perenne hybrids. Nevertheless, Książczyk et al. [17] described the unstable number (from four to six hybridization signals) of 5S rDNA in F2-F4 generations of F. pratensis × L. perenne hybrid. In this plant material, the additional smaller hybridization sites were located in different regions of F. pratensis chromosomes and this was not related to the changed total number of chromosomes. Pedrosa-Harand et al. [31] suggested that variability in chromosomal positions, numbers and signal sizes of rDNA can be caused by at least three different mechanisms: (i) structural rearrangements that affect the position of rDNA loci; (ii) amplification of satellite DNA with no impact on the order of other sequences in chromosomes; and (iii) dispersion of rDNA repeats, amplification of the new minor loci and deletion of original major loci. The first two mechanisms are connected with changes in the chromosomal position of rDNA sequences, and the third one is responsible for changing the number and size of rDNA loci. Chromosomal re-patterning can be also connected with mobile elements activity [32], such as – Enhancer/Suppressor-mutator (En/Spm) transposons in Aegilops species [33] or Wukong terminal repeat retrotransposons in miniature (TRIM) in maize germplasm [34]. In interspecific hybrids, especially in newly synthesized hybrids, the targets of structural reconstructions can be rDNA sequences [35], what was also found in our plant material. In these plants, the rearrangements between F. pratensis and L. perenne genomes in/or near rDNA position took place. Huang et al. [36] reported that fragile sites (FSs) in chromosomes of L. perenne are linked with 35S rDNA cluster. The FSs are determined as regions that are prone to forming gaps or breaks on chromosomes. Human fragile sites are known as hotspots for chromosome aberrations, especially in many tumor cells [37, 38]. However, in Festuca × Lolium plants, the analysis of the relation between FSs associated with 35S rDNA regions and the location of structural rearrangements revealed that intergenomic changes did not correspond to FSs [39].

In theory, chromosomal rearrangements should not occur in allopolyploids due to preferential homologous chromosome pairing and sophisticated molecular processes preventing homoeologous recombination, such as Ph1 locus in wheat [40] or pairing regulators in Festuca-Lolium complex [41,42,43,44]. It is worth to highlight the differences between pairing control genes in both genera – in Festuca species genes are haplo-insufficient, while in Lolium genus this genetic system is effective both at the hetero- and homozygous stage [45]. Nevertheless, closely related plant species usually display the synteny of genomes and the collinearity of gene order, what can contribute to mispairing [46]. In the analyzed F2-F3 generations of F. pratensis × L. perenne, the recombinations between both genomes took place. Rearrangements that were recognized in F2 progeny are in contrary to the idea of Ramsey and Schemske [47]. According to their theory, homoeologous pairing is suspected to be avoided during meiosis in the F1 generation following interspecific hybridization, because it can lead to infertility of progeny. Additionally, the increasing number of structural rearrangements in the next generations of F. pratensis × L. perenne hybrid was observed, what is consistent with our previous findings [17, 48, 49]. Beyond the frequent and increasing occurrence of homoeologous rearrangements among F2-F3 plants, we identified also the examples of non-homologous recombination. This type of recombination between non-rDNA bearing chromosome of F. pratensis and region of chromosome 3L (chromosomal segment bearing 5S rDNA locus) took place (Fig. 3). Non-homologous chromosome pairing and recombination could be provoked by cross hybridization and rapid genome changes. It is thought that both types of recombination probably play a significant role in genome reshuffling and stabilization. Genomes of many polyploids, such as Triticum [50], Aegilops [51], Kengyilia species [52] or × Triticosecale [53] underwent genome reshufflings during the course of evolution (especially translocation, deletions and inversions). Chromosomal aberrations are mainly located in heterochromatic regions composed of repetitive DNA. The analysis of rearrangement distribution in F. pratensis × L. perenne amphiploid forms revealed terminal location as the most frequent (see Additional file 2: Table S1). The recombination rates along individual Festuca and Lolium chromosomes evaluated by Kopecký et al. [54] and King et al. [55] revealed that recombination sites are concentrated in intercalary regions of chromosome arms and reduced in their proximal and distal parts.

rDNA-FISH combined with GISH approach allowed us to unequivocally discriminate selected rDNA-bearing chromosome pairs in the complement of analyzed plant material. For these pairs of chromosomes, the analysis was more detailed and revealed, which pair was the most recombined. Knowledge about the types of rearrangements and their frequencies is essential to better understand the evolutionary significance of genome reshufflings [56]. In both studied generations, the high frequency of recombination was calculated for Festuca rDNA-bearing chromosomes, in F2 plants – for chromosome 2F, while in F3 plants – for chromosome 3F. The high value of rearrangements in chromosome 2F was also depicted by Majka et al. [39] in F2-F9 generations of F. pratensis × L. perenne hybrid. Furthermore, the recognition of a particular chromosome pair revealed that in F3 prevailed plants with only one chromosome 2F. This odd number of the chromosome was inherited from parental forms – maternal plant had two chromosome 2F, while the pollinator did not have this pair of chromosomes. This unequivocal number of chromosomes 2F was mainly related to a proper (2n = 28) number of chromosomes. The loss of chromosome 2F could be also associated with deletion of 35S rDNA sequence or translocation of 35S rDNA into another chromosome. Jumping of rDNA sequences, especially 35S rDNA, is not the rare case in the plant kingdom [30, 33, 57].

The evolution of plant genomes, especially hybrids (including newly created by human and ancient natural hybrids) is related to genome and chromosome repatterning. In hybrid plants, the extensive interactions between genomes are observed. For hybrids derived from Festuca-Lolium complex (eight successive generations of F. pratensis × L. perenne hybrid), Zwierzykowski et al. [48, 49] observed preferential transmission of Lolium chromosomes, probably due to material selection for agronomic traits. However, the comparison of selected and unselected plants of F2-F4 generations of F. pratensis × L. perenne hybrid confirmed the dominance of Lolium over Festuca chromatin [58]. Our material was composed of plants without any selection and we also observed a higher number of Lolium chromosomes and drift in favor of Lolium in successive generations. Statistical analysis confirmed our hypothesis that asymmetrical variation between F. pratensis and L. perenne parental genomes exists. Nonetheless, Akiyama et al. [59] revealed no preferential transmission of the Lolium genome in consecutive generations of L. muliflorum (2x) × F. arundinacea (6x) hybrids.

The inheritance of ISSR loci

During meiosis in autopolyploids, chromosomes may form multivalents leading to multisomic inheritance, where all the possible allelic combinations are formed in equal frequencies. In allopolyploids chromosomes usually pair with their homologs and form bivalents resulting in disomic inheritance. It has been observed, however, that both auto- and allopolyploids may shift between those two modes of inheritance and show intermediate or mixed inheritance [47], in which not all chromosomes have the same pairing behavior [6]. In addition, in the tetraploids two types of meiotic segregation may occur: random chromosome or random chromatid segregation, depending on the formation of bivalents or quadrivalents, respectively.

Our results suggest that in the F1 hybrids most of ISSR loci segregated to gametes in a tetrasomic way, although the relatively small number of F2 individuals did not allow us to determine which type of assortment: chromosome or chromatid might take place during meiosis of the F1 hybrid. There was just one locus, which segregation suggested the possibility of quadrivalents formation and random chromatid segregation. Only in a few loci, disomic segregation could be accepted based on the χ2 test result. Their occurrence, however, proved the possibility of a mixed model of inheritance in Festulolium hybrid plants. From ISSR data only, it is not possible to decide whether there is any specific preference of chromosome pairing in particular chromosome homologue or homoeologue group. Nevertheless, from the nature of ISSR marker system, it may be assumed that the ISSR loci represent random positions in the genome and will be located in different chromosomes, suggesting that most of the chromosomes in F1 hybrids may form intergenomic bivalents or multivalents during meiosis. As indicated by Zwierzykowski et al. [60] who studied the meiosis of F. pratensis × L. perenne hybrids, the F. pratensis/L. perenne intergenomic bivalents were observed in about 33% of analyzed cells. A significant number of quadrivalents and some trivalents or univalent were also found. Such observations support the possibility of tetrasomic way of inheritance of ISSR loci, although, further analysis, for example with the use of ISSR amplification products as probes for in situ hybridization, may give more inside into this problem and shed more light on the chromosome vs. chromatid model of assortment of genetic material into the gametes.

An interesting observation was the occurrence of three ISSR loci, which seem to be inherited in different ways in both F2 combinations. It may suggest that the choice between bivalent or multivalent pairing, at least for some chromosomes in the F1 meiosis may be equally possible and take place at random.

Multivalent formation followed by a crossing-over located between the centromere and particular locus may in some cases lead to the migration of sister alleles to the same gamete during meiosis – a phenomenon known as a double reduction [61]. It may have a role in the purging of deleterious mutations through gametophytic selection [62], and could affect the distribution of gene frequency in polyploid populations [63]. In our study, a decline of ISSR polymorphism was found between subsequent generations. We understand that the dominant nature of ISSR markers does not allow to precisely measure the diversity, especially in tetraploid species, but it cannot be excluded that double reduction together with other mechanisms may be one of the reasons for this observation.

Our analysis of F3 plants shows that the mixed model of the inheritance probably persists in the subsequent generation, although this conclusion may only be speculative as almost all ISSR loci in F3 generation were not assigned to any specific model of inheritance. This high proportion of ISSR loci with an unknown type of inheritance, is another important factor to consider, as it may in some way bias our analysis. We cannot exclude, however, that irregular cytogenetic behaviors, such as the formation of univalents and trivalents, as evidenced in Zwierzykowski et al. [60] may contribute to the events with unknown inheritance mode. Here, further extending the sample size or the application of co-dominant markers and ISSR-FISH analysis may allow to specify the type of inheritance in F. pratensis × L. perenne hybrids.

Segregation distortion of ISSR loci

Our ISSR analysis showed that only the minority of ISSR loci in F2 and F3 generations were characterized by segregation distortion (SD). They accounted for about 3.0% and 15.0% in two F2 cross combinations and for 5.0% and 8.4% in the two F3 cross combinations. We are aware, however, that we might underestimate the frequency of this phenomenon, as the dominant type of ISSR markers does not allow to fully assign the genotypes at particular loci and, consequently, to test only these segregation ratios that are unambiguous. The possibility of underestimation of this result may be suggested by other studies of Festuca and/or Lolium hybrids. For example, the genetic mapping studies in L. perenne showed that depending on the mapping population used, the proportion of distorted loci ranged from 12% [64] to 60% [65]. In the mapping study of F. arundinacea, segregation distortion was observed from 23% [66] to 45% of loci [67]. The mechanisms underlying SD are only partially understood and may involve various processes, such as differential inclusion of parental alleles into gametes (meiotic drive), differential fertilization success of specific gametes or differences in survival of haploid gametes or zygotes with specific allelic composition [68]. Some evidence suggests that SD is higher in interspecific than in intraspecific crosses [69, 70] thus being a reflection of more dynamic changes in the interspecific hybrids genomes.

The use of ISSR markers in our study suggest that there were no preferences in the transmission of either Festuca or Lolium ISSR alleles to the following generations of F. pratensis × L. perenne hybrid, nor there were any preferences either for transmission of the dominant or recessive allele. Some preferences toward the transmission of the recessive allele were noticed in the F3 combination no. 4, although it may be related to the relatively low number of loci under SD detected, and happened by a chance only. However, one concurrence may be noticed – more distorted loci were found in these populations, which parents were more polymorphic. Successful amplification of ISSR locus relies on two factors: (i) the presence of two microsatellite sequences in the reverse complement orientation in DNA strand that are complementary to the oligonucleotide used for the PCR reaction, and (ii) the location of such sequences in a distance short enough to be amplified by routinely used polymerases and PCR assays. If the second mechanism will be considered, the lack of ISSR band will be a result of the too long distance between priming sites. This opens a possibility to think about the relationship between the length of Festuca and Lolium chromosomes and their possibility of successful pairing during meiosis. Distorted segregation may occur in these loci where an imbalance of a homoeologous chromosome segment length is present. Such hypothesis assumes that these chromosomes pair as multivalents and further suggests the intermediate or mixed model of inheritance in Festuca × Lolium hybrids.


The cytogenetic analysis of 134 plants of F2-F3 generations of F. pratensis × L. perenne hybrid, including rDNA-FISH and GISH together with molecular genotyping of all plants by ISSR markers allowed us to draw conclusions from data reported herein. Cytogenetic genotyping of all analyzed hybrid plants revealed discrepancies in parental genomes (Festuca and Lolium) behavior. Statistical analysis of that data implies a genome drift over successive generations in allotetraploid F. pratensis × L. perenne hybrid. We observed the increasing number of Lolium chromosomes in the next generations of hybrids and significant differences between F2-F3 generations for recombined rDNA-bearing L. perenne chromosomes. Despite the decreasing total number of F. pratensis chromosomes, their chromosomes were more rearranged and statistical differences between analyzed generations for recombined rDNA-bearing and non-rDNA-bearing ones were recognized. Nevertheless, genotyping by ISSR markers analysis evidenced no preference in the transmission of either Festuca or Lolium alleles to the following generations of F. pratensis × L. perenne hybrid. It is important, however, to bear in mind the different resolutions of both cytogenetic and molecular markers analysis. Furthermore, the analysis of ISSR loci inheritance in these hybrid plants revealed that the tetrasomic type of inheritance was characteristic for the majority of ISSR loci, but in some cases, the disomic type was also observed. This indicates that the mixed or intermediate model of alleles transmission may occur in the hybrids. Such possibility is in the agreement with the detection of non-homologous recombination events, as the multivalent formation may promote rearrangements between chromosomes from different genomes. Further examinations of the possible rearrangements combining chromosome and DNA level of resolution may give deeper insight into the nature of how the genome balance is achieved in F. pratensis × L. perenne hybrid.


Plant material

Allotetraploid F1 plants (2n = 4x = 28) were produced by crossing the autotetraploid (2n = 4x = 28) F. pratensis Huds. (cv. ‘Westa’; Fp 11/59) (♀) with the autotetraploid (2n = 4x = 28) L. perenne L. (cv.‘Solen’; Lp 08/33) (♂). The F1 maternal plant (F1 6-2A) and two different F1 pollinators were used to produce two combinations of the F2 generation: no. 1 – with F1 pollinator 6-7A and no. 2 – with F1 pollinator 6-3B. The F3 generation (combinations no. 3 and 4) was obtained by crossing four randomly selected F2 plants derived from combination no. 1. The detailed scheme of crossing experiments was presented in Additional file 1: Fig. S1. All plant material derived from the collection of the Institute of Plant Genetics, Polish Academy of Sciences in Poznan, Poland.

Preparation of root meristems

Roots from plants growing in a greenhouse were excised and pretreated in ice-cold water for 24 h, fixed in 3:1 (v/v) ethanol:glacial acetic acid and then stored at − 20 °C until required. Mitotic chromosome preparations were made using the excised roots, according to the procedure described in Majka et al. [71]. The high-quality slides were frozen and stored at 4 °C until used.

Fluorescent in situ hybridization

The FISH procedure was carried out according to the protocol published in Majka et al. [71]. The hybridization mixture consisted of 50% deionized formamide, 10% dextran sulfate, 2 × SSC, 0.5% SDS, 10 μg salmon sperm blocking DNA and 100-120 ng/slide both rDNA probes. In FISH experiments the following probes were used: (i) the 5S rDNA sequence, generated from the wheat clone pTa794 [72] and (ii) 35S rDNA sequence, generated from a 2.3-kb ClaI sub-clone of the 25S rDNA coding region of Arabidopsis thaliana [73]. 5S rDNA was labeled by PCR technique with universal M13 primers by using tetramethyl-rhodamine-5-dUTP (Roche, Mannheim, Germany), while 35S rDNA probe was prepared by nick translation with digoxigenin-11-dUTP (Roche, Mannheim, Germany). Chromosome preparations were denatured at 80 °C for 2 min in presence of hybridization mixture and hybridized overnight in a humid chamber at 37 °C. Post-hybridization washes in 0.1 × SSC at 42 °C provided an equivalent of 79% stringency. Digoxigenin-labeled probes were detected by fluorescein isothiocyanate-conjugated (FITC) anti-digoxigenin antibody (Roche, Mannheim, Germany). Chromosomes were counterstained with 2.5 mg/ml 4′,6-diamidino-2-phenylindole (DAPI, Sigma, St. Louis, Missouri) in antifade Vectashield solution (Vector Laboratories, Burlingame, CA, USA).

The chromosome nomenclature used for rDNA-bearing chromosomes was developed by Thomas et al. [74]. According to this, there are two rDNA-bearing chromosomes in F. pratensis karyotype: chromosomes 2F (bearing 35S rDNA loci) and 3F (bearing 5S rDNA). Whereas in the L. perenne karyotype, chromosomes 3L (bearing 5S and 35S rDNA in opposite arms) and a group of chromosomes 1, 2 and 7L (bearing 35S rDNA) can be recognized.

Genomic in situ hybridization

For the GISH experiments all the chromosome spreads, which had been used in FISH, were applied. Therefore, before GISH, the sets of probes from previous experiments were rinsed off. The GISH procedure was carried out according to the protocol of Majka et al. [39]. Genomic DNA of L. perenne, which was mechanically sheared to fragments of 5-10 kb by boiling for 30-45 min and labeled by nick translation with digoxigenin-11-dUTP (Roche, Mannheim, Germany), was used as a probe in GISH experiments. While the nuclear DNA of F. pratensis was used as a blocking DNA (sheared to the shorter fragments by boiling for 10-15 min). The probe was detected using a FITC conjugated anti-digoxigenin antibody. The chromosome preparations were counterstained in Vectashield (Vector Laboratories, Burlingame, CA, USA) containing 1.5 g/ml− 1 of propidium iodide.

Image capturing and data processing

Images were captured using an Olympus XM10 CCD camera attached to an Olympus BX 61 automatic epifluorescence microscope. Final image adjustments were done with Olympus Cell-F imaging software (ver. 3.1; Olympus Soft Imaging Solutions GmbH, Germany) and Micrographx Picture Publisher software (ver. 10; Corel Corporation, Canada).

Cytogenetic data were statistically processed by the χ2 test and t-test within GenStat® (ver. 16; VSN Int.). Analysis was performed to compare the distributions of selected cytogenetic traits within and between F2 and F3 generations, as well as to evaluate the difference between both parental genomes in F2-F3 generations of F. pratensis × L. perenne hybrid. The comparison between F2 and F3 generations was performed for the F2 cross combination no. 1 only, because of crossing scheme – combinations no. 3 and 4 (F3 generation) were obtained by crossing plants that derived from combination no. 1 of F2 generation.

DNA isolation and ISSR-PCR analysis

DNA samples were extracted from young leaves using CTAB method [75]. ISSR-PCRs were carried out in a volume of 25 μl containing 70 ng of genomic DNA, 2 × PCR buffer, 2.5 mM MgCl2, 2.5 mM dNTP, 10 mM of each primer and 0.25 units of Taq DNA Polymerase (Thermo Scientific, Waltham, Massachusetts). The thermal profile for ISSR-PCR was as follows: 94 °C for 2 min – initial denaturation, then 35 cycles of 95 °C for 30 s, the temperature of annealing was depended on primer (48–57 °C) for 1 min and elongation at 72 °C for 1 min. A final extension step – 5 min at 72 °C. The temperature of annealing and sequences of analyzed oligonucleotides were shown in Additional file 10: Table S9. ISSR amplification products were separated in 2.0% agarose gels (EURx, Gdańsk, Poland), stained with ethidium bromide (Sigma, St. Louis, Missouri) and visualized under UV light and photographed (Syngen UV visualizer). GeneRuler Ladder (Novazym, Poznań, Poland) was used to determine the size of the DNA fragments.

The analysis of ISSR loci polymorphism

In order to estimate the level of polymorphism between parental plants: F. pratensis (Fp 11/59) and L. perenne (Lp 08/33) and between pairs of the F1 and F2 individuals used as parents for subsequent generations, all the bands produced by each ISSR primer were counted. The entire set of bands was composed of monomorphic products present in both analyzed forms and polymorphic products characteristic for individual genotypes. The level of polymorphism was defined as the percentage of polymorphic bands for all the primer combinations (see Additional file 11: Table S10 and Additional file 12: Table S11). The Polymorphism Information Content (PIC) was calculated using iMEC program [76].

The analysis of ISSR loci segregation

Based on the presence or absence of ISSR band, the possible genotypes were assigned to parental F. pratensis and L. perenne plants and to the F1 and F2 genotypes which were used as parents for subsequent cross combinations. As both parental plants (Fp 11/59 and Lp 08/33) are autotetraploids, we assumed that their chromosomes will segregate to gametes in a tetrasomic way and we have taken into account the possibility of random chromosome or random chromatic segregation. The individuals of F1 and F2 generation were allotetraploids, but our observations suggest that both tetrasomic and disomic inheritance of chromosomes is possible, although it was found that bivalent formation is more frequent in the hybrid [77]. Taking into account this observation as well as the dominant type of ISSR markers that were used in our study and the relatively low sample size, which does not allow to distinguish between chromosome vs. chromatid assortment, we made most of the later calculations based on the simplified assumption of the random chromosome segregation model. This allowed us to divide the types of inheritance into tetrasomic and disomic, based on the different possible rates of band segregation, depending on the genotypes in each ISSR locus in parental individuals (see Additional file 13: Table S12). All possible segregation ratios for random chromosome segregation model were tested using the χ2 test. We assumed that the lack of segregation distortion is more probable than its occurrence. In all cases, where only one segregation ratio passed the χ2 test (P < 0.05), the one was selected as describing the most probable inheritance type in particular ISSR locus. When more segregation rates passed the χ2 tests and some of them were characteristic to different inheritance types, we assumed such result as ambiguous with unknown inheritance type. When none of the segregation rates passed the χ2 test, such locus was treated as segregation distortion and a possible signature of genomic changes between parental plants and the subsequent F2 and/or F3 generations. To confirm the probability of segregation distortion of such locus it was additionally analyzed for those segregation ratios that are characteristic for random chromatid segregation model, which assumes quadrivalent formation (see Additional file 13: Table S12).





Fluorescent in situ hybridization


Fragile Sites


Genomic in situ hybridization


Inter-Simple Sequence Repeats


Polymorphism Information Content


Segregation Distortion


  1. 1.

    Leitch AR, Leitch IJ. Genomic plasticity and the diversity of polyploid plants. Science. 2008;320(5875):481–3.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Murat F, Xu JH, Tannier E, Abrouk M, Guilhot N, Pont C, et al. Ancestral grass karyotype reconstruction unravels new mechanisms of genome shuffling as a source of plant evolution. Genome Res. 2010;20:1545–57.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Salse J. Deciphering the evolutionary interplay between subgenomes following polyploidy: a paleogenomics approach in grasses. Am J Bot. 2016;103(7):1167–74.

    Article  PubMed  Google Scholar 

  4. 4.

    Wendel JF. Genome evolution in polyploids. Plant Mol Biol. 2000;42:225–49.

    CAS  Article  Google Scholar 

  5. 5.

    Catalan P, Gabriel SMJ, Palop-Esteban M, Moreno C, Gonzalez-Candelas F. A Bayesian approach for discriminating among alternative inheritance hypotheses in plant polyploids: the allotetraploid origin of genus Borderea (Dioscoreaceae). Genetics. 2006;172:1939–53.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Jeridi M, Perrier X, Rodier-Goud M, Ferchichi A, D'Hont A, Bakry F. Cytogenetic evidence of mixed disomic and polysomic inheritance in an allotetraploid (AABB) Musa genotype. Ann Bot. 2012;110(8):1593–606.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Reflinur KB, Jang SM, Chu S. Analysis of segregation distortion and its relationship to hybrid barriers in rice. Rice. 2014;7:3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Thomas HM, Thomas BJ. Meiosis of triploid Lolium. I. Synaptonemal complex formation and chromosome configurations at metaphase I in aneuploid autotriploid L. multiflorum. Genome. 1994;37:181–9.

    CAS  Article  Google Scholar 

  9. 9.

    Zwierzykowski Z, Tayyar R, Brunell M, Lukaszewski AJ. Genome recombination in intergeneric hybrids between tetraploid Festuca pratensis and Lolium multiflorum. J Hered. 1998;89(4):324–8.

    Article  Google Scholar 

  10. 10.

    Zhang H, Bian Y, Gou X, Zhu B, Xu C, et al. Persistent whole-chromosome aneuploidy is generally associated with nascent allohexaploid wheat. Proc Natl Acad Sci U S A. 2013;110:3447–52.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Kopecký D, Bartos J, Lukaszewski AJ, Baird JH, Černoch V, Kölliker R, et al. Development and mapping of DArT markers within the Festuca-Lolium complex. BMC Genomics. 2009;10:473.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Jang TS, Emadzade K, Parker J, Temsch EM, Leitch AR, Speta F, Weiss-Schneeweiss H. Chromosomal diversification and karyotype evolution of diploids in the cytologically diverse genus Prospero (Hyacinthaceae). BMC Evol Biol. 2013;13:136.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Thomas HM, Morgan WG, Meredith MR, Humphreys MW, Thomas H, Leggett JM. Identification of parental and recombined chromosomes in hybrid derivatives of Lolium multiflorum x Festuca pratensis by genomic in situ hybridization. Theor Appl Genet. 1994;88:909–13.

    CAS  Article  Google Scholar 

  14. 14.

    King IP, Morgan WG, Armstead IP, Harper JA, Hayward MD, Bollard A, et al. Introgression mapping in the grasses I. introgression of Festuca pratensis chromosomes and chromosome segments into Lolium perenne. Heredity. 1998;81:462–7.

    CAS  Article  Google Scholar 

  15. 15.

    Kopecký D, Loureiro J, Zwierzykowski Z, Ghesquière M, Doležel J. Genome constitution and evolution in Lolium × Festuca hybrid cultivars (Festulolium). Theor Appl Genet. 2006;113:731–42.

    CAS  Article  Google Scholar 

  16. 16.

    Kosmala A, Zwierzykowski Z, Gasior D, Rapacz M, Zwierzykowska E, Humphreys MW. GISH/FISH mapping of genes for freezing tolerance transferred from Festuca pratensis to Lolium multiflorum. Heredity. 2006;96:243–51.

    CAS  Article  Google Scholar 

  17. 17.

    Książczyk T, Zwierzykowska E, Molik K, Taciak M, Krajewski P, Zwierzykowski Z. Genome-dependent chromosome dynamics in three successive generations of the allotetraploid Festuca pratensis × Lolium perenne hybrid. Protoplasma. 2015;252:985–96.

    Article  PubMed  Google Scholar 

  18. 18.

    Thomas H, Humphreys MO. Progress and potential of inter-specific hybrids of Lolium and Festuca. Agric Sci. 1991;117:1–8.

    Article  Google Scholar 

  19. 19.

    Humphreys MW, Zare AG, Panakinskiene I, Thomas H, Rogers WJ, Collin HA. Interspecific genomic rearrangements in androgenic plants derived from a Lolium multiforum × Festuca arundinacea (2n = 5x = 35) hybrid. Heredity. 1998;80:78–82.

    CAS  Article  Google Scholar 

  20. 20.

    Udall JA, Quijada PA, Osborn TC. Detection of chromosomal rearrangements derived from homeologous recombination in four mapping populations of Brassica napus L. Genetics. 2005;169:967–79.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Madlung A, Tyagi AP, Watson B, Jiang HM, Kagochi T, Doerge RW, et al. Genomic changes in synthetic Arabidopsis polyploids. Plant J. 2005;41:221–30.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Vitte C, Panaud O. LTR retrotransposons and flowering plant genome size: emergence of the increase/decrease model. Cytogen Genom Res. 2005;110:91–107.

    CAS  Article  Google Scholar 

  23. 23.

    Vitte C, Peterson T, Bennetzen JL. Transposable elements and the evolving structure of the maize genome, in: Advances in maize. In: Prioul JL, Thevenot C, Molnar T, editors. the Society for Experimental Biology, London, UK: Taylor & Francis, vol. 3, 2-36. p. 2011.

  24. 24.

    Lim KY, Kovarik A, Matyasek R, Bezdeek M, Lichtenstein CP, Leitch AR. Gene conversion of ribosomal DNA in Nicotiana tabacum is associated with undermethylated, decondensed and probably active gene units. Chromosoma. 2000;109:161–72.

    CAS  Article  Google Scholar 

  25. 25.

    Pontes O, Neves N, Silva M, Lewis MS, Madlung A, Comai L, et al. Chromosomal locus rearrangements are a rapid response to formation of the allotetraploid Arabidopsis suecica genome. Proc Natl Acad Sci U S A. 2004;101:18240–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Ozkan H, Levy AA, Feldman M. Allopolyploidy-induced rapid genome evolution in the wheat (Aegilops-Triticum) group. Plant Cell. 2001;13:1735–47.

    CAS  Article  Google Scholar 

  27. 27.

    Shaked H, Kashkush K, Ozkan H, Feldman M, Levy AA. Sequence elimination and cytosine methylation are rapid and reproducible responses of the genome to wide hybridization and allopolyploidy in wheat. Plant Cell. 2001;13:1749–59.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Renny-Byfield S, Kovařík A, Chester M, Nichols RA, Macas J, Novák P, Leitch AR. Independent, rapid and targeted loss of highly repetitive DNA in natural and synthetic allopolyploids of Nicotiana tabacum. PLoS One. 2012:e36963.

    CAS  Article  Google Scholar 

  29. 29.

    Linares C, Gonzalez J, Ferrer E, Fominaya A. The use of double fluorescence in situ hybridization to physically map the positions of 5S rDNA genes in relation to the chromosomal location of 18S-5.8S-26S rDNA and a C genome specific DNA sequence in the genus Avena. Genome. 1996;39:535–42.

    CAS  Article  Google Scholar 

  30. 30.

    Lan T, Albert VA. Dynamic distribution patterns of ribosomal DNA and chromosomal evolution in Paphiopedilum, a lady’s slipper orchid. BMC Plant Biol. 2011;11:126.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Pedrosa-Harand A, Souza de Almeida CC, Mosiolek M, Blair MW, Schweizer D, Guerra M. Extensive ribosomal DNA amplification during Andean common bean (Phaseolus vulgaris L.) evolution. Theor Appl Genet. 2006;112:924–33.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Zhang J, Yu C, Krishnaswamy L, Peterson T. Transposable elements as catalysts for chromosome rearrangements. Methods Mol Biol. 2010;701:315–26.

    CAS  Article  Google Scholar 

  33. 33.

    Raskina O, Belyayev A, Nevo E. Quantum speciation in Aegilops: molecular cytogenetic evidence from rDNA cluster variability in natural populations. Proc Natl Acad Sci U S A. 2004;101:14818–23.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Liu Z, Li X, Messing J, Xu JH. The Wukong terminal-repeat Retrotransposon in miniature (TRIM) elements in diverse maize Germplasm. G3. 2015;5:1585–92.

    Article  PubMed  Google Scholar 

  35. 35.

    Volkov RA, Komarova NY, Hemleben V. Ribosomal DNA in plant hybrids: inheritance, rearrangement, expression. Syst Biodivers. 2007;5(3):261–76.

    Article  Google Scholar 

  36. 36.

    Huang J, Ma L, Yang F, S-z F, Li L. 45S rDNA regions are chromosome fragile sites expressed as gaps in vitro on metaphase chromosomes of root-tip meristematic cells in Lolium spp. PLoS One. 2008;3:e2167.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Debatisse M, Le Tallec B, Letessier A, Dutrillaux B, Brison O. Common fragile sites: mechanisms of instability revisited. Trends Genet. 2012;28(1):22–32.

    CAS  Article  Google Scholar 

  38. 38.

    Glover TW. Fragile sites in cancer: more than meets the eye. Nat Rev Cancer. 2017;17(8):489–501.

    CAS  Article  Google Scholar 

  39. 39.

    Majka J, Zwierzykowski Z, Majka M, Kosmala A. Karyotype reshufflings of Festuca pratensis × Lolium perenne hybrids. Protoplasma. 2018;255:451–8.

    Article  PubMed  Google Scholar 

  40. 40.

    Riley R, Chapman V. Genetic control of cytologically diploid behaviour of hexaploid wheat. Nature. 1958;182:713–5.

    Article  Google Scholar 

  41. 41.

    Jauhar PP. Genetic control of diploid-like meiosis in hexaploid tall fescue. Nature. 1975;254:595–7.

    CAS  Article  Google Scholar 

  42. 42.

    Jauhar PP. Genetic regulation of diploid-like chromosome pairing in the hexaploid species, Festuca arundinacea Schreb. And F. rubra L. (Gramineae). Chromosoma. 1975;52:363–82.

    Article  Google Scholar 

  43. 43.

    Evans GM, Davies EW. The genetics of meiotic chromosome pairing in Lolium temulentum × Lolium perenne tetraploids. Theor Appl Genet. 1985;71:185–92.

    CAS  Article  Google Scholar 

  44. 44.

    Armstead IP, Bollard A, King IP, Forster JW, Hayward MD, Evans GM, Thomas HM. Chromosome pairing in Lolium perenne × L. temulentum diploid hybrids: genetic and cytogenetic evaluation. Heredity. 1999;83:298–303.

    Article  Google Scholar 

  45. 45.

    Jenczewski E, Alix K. From diploids to allopolyploids: the emergence of efficient pairing control genes in plants. Crticial Rev Plant Sci. 2004;23(1):21–45.

    CAS  Article  Google Scholar 

  46. 46.

    Taylor SJ, Rojas LD, Ho SW, Martin NH. Genomic collinearity and the genetic architecture of floral differences between the homoploid hybrid species Iris nelsonii and one of its progenitors, Iris hexagona. Heredity. 2013;110:63–70.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Ramsey J, Schemske DW. Neopolyploidy in flowering plants. Annu Rev Ecol Syst. 2002;33:589–639.

    Article  Google Scholar 

  48. 48.

    Zwierzykowski Z, Kosmala A, Zwierzykowska E, Jones N, Jokś W, Bocianowski J. Genome balance in six successive generations of the allotetraploid Festuca pratensis × Lolium perenne. Theor Appl Genet. 2006;113:539–47.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Zwierzykowski Z, Zwierzykowska E, Taciak M, Kosmala A, Jones RN, Zwierzykowski W, et al. Genomic structure and fertility in advanced breeding populations derived from an allotetraploid Festuca pratensis × Lolium perenne cross. Plant Breed. 2011;130:476–80.

    Article  Google Scholar 

  50. 50.

    Ma J, Stiller J, Berkman PJ, Wei Y, Rogers J, Feuillet C, et al. Sequence-based analysis of translocations and inversions in bread wheat (Triticum aestivum L.). PLoS One. 2013;8(11):e79329.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Badaeva ED, Amosova AV, Samatadze TE, Zoshchuk SA, Shostak NG, Chikida NN, et al. Genome differentation in Aegilops. 4. Evolution of the U-genome cluster. Plant Syst Evol. 2004;246:45–76.

    CAS  Article  Google Scholar 

  52. 52.

    Wang Q, Liu H, Gao A, Yang X, Liu W, Li X, Li L. Intergenomic rearrangements after polyploidization of Kengyilia thoroldiana (Poaceae: Triticeae) affected by environmental factors. PLoS One. 2012;7:e31033.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Bento M, Pereira HS, Rocheta M, Gustafson P, Viegas W, Silva M. Polyploidization as a retraction force in plant genome evolution: sequence rearrangements in triticale. PLoS One. 2008;3(1):e1402.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Kopecký D, Havránková M, Loureiro J, Castro S, Lukaszewski AJ, Bartoš J, et al. Physical distribution of homoeologous recombination in individual chromosomes of Festuca pratensis in Lolium multiflorum. Cytogen Gen Res. 2010;129:162–72.

    Article  Google Scholar 

  55. 55.

    King J, Armstead I, Harper J, Ramsey L, Snape J, Waugh R, et al. Exploitation of interspecific diversity for monocot crop improvement. Heredity. 2013;110:475–83.

    CAS  Article  Google Scholar 

  56. 56.

    Badaeva ED, Dedkova OS, Gay G, Pukhalskyi VA, Zelenin AV, Bernard S, Bernard M. Chromosomal rearrangements in wheat: their types and distribution. Genome. 2007;50:907–26.

    CAS  Article  Google Scholar 

  57. 57.

    Schubert I, Wobus U. In situ hybridization confirms jumping nucleolus organizing regions in Allium. Chromosoma. 1985;92:143–8.

    Article  Google Scholar 

  58. 58.

    Zwierzykowski Z, Książczyk T, Taciak M, Zwierzykowska E, Jones N, Kosmala A. Genome constitution in selected and unselected plants of F2-F4 generations derived from an allotetraploid Festuca pratensis × Lolium perenne hybrid. In: Barth S, Milbourne D, editors. Breeding Strategies for Sustainable Forage and Turf Grass Improvement: Springer Science+Business Media, Dodrecht; 2012. p. 75–9.

    Google Scholar 

  59. 59.

    Akiyama Y, Ueyama Y, Hamada S, Kubota A, Kato D, Yamada-Akiyama H, et al. Utilization of flow cytometry for Festulolium breeding (Lolium multiflorum (2x) × Festuca arundinacea (6x)). Breed Sci. 2016;66:234–43.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Zwierzykowski Z, Zwierzykowska E, Taciak M, Jones N, Kosmala A, Krajewski P. Chromosome pairing in allotetraploid hybrids of Festuca pratensis × Lolium perenne revealed by genomic in situ hybridization (GISH). Chromosom Res. 2008;16(4):575–85.

    CAS  Article  Google Scholar 

  61. 61.

    Gupta PK. Cytogenetics: An Advanced Study. Meerut New Delhi: Rastogi Publications; 2015.

    Google Scholar 

  62. 62.

    Butruille DV, Boiteux LS. Selection-mutation balance in polysomic tetraploids: impact of double reduction and gametophytic selection on the frequency and subchromosomal localization of deleterious mutations. Proc Natl Acad Sci U S A. 2000;97:6608–13.

    CAS  Article  Google Scholar 

  63. 63.

    Wright S. The distribution of gene frequencies in populations of polyploids. Proc Natl Acad Sci U S A. 1938;24:372–7.

    CAS  Article  Google Scholar 

  64. 64.

    Bert PF, Charmet G, Sourdille P, Hayward MD, Balfourier F. A high-density molecular map for ryegrass (Lolium perenne) using AFLP markers. Theor Appl Genet. 1999;99:445–52.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Jensen LB, Andersen JR, Frei U, Xing Y, Taylor C, Holm PB, Lübberstedt T. QTL mapping of vernalisation response in perennial ryegrass (Lolium perenne L.) reveals co-location with an orthologue of wheat VRN1. Theor Appl Genet. 2005;110:527–36.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Saha MC, Mian R, Zwonitzer JC, Chekhovskiy K, Hopkins AA. An SSR- and AFLP-based genetic linkage map of tall fescue (Festuca arundinacea Schreb.). Theor Appl Genet. 2005;110:323–36.

    CAS  Article  Google Scholar 

  67. 67.

    Dierking R, Azhaguvel P, Kallenbach R, Saha M, Bouton J, Chekhovskiy K, et al. Linkage maps of a mediterranean × continental tall fescue population and their comparative analysis with other Poaceae species. Plant Genome. 2015, 2015;8.

    CAS  Article  Google Scholar 

  68. 68.

    Fishman L, Willis JH. A novel meiotic drive locus almost completely distorts segregation in mimulus (monkeyflower) hybrids. Genetics. 2005;169(1):347–53.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Causse MA, Fulton TM, Cho YG, Ahn SN, Chunwongse J, Wu K, et al. Saturated molecular map of the rice genome based on an interspecific backcross population. Genetics. 1994;138(4):1251–74.

    CAS  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Jenczewski E, Gherardi M, Bonnin I, Prosperi JM, Olivieri I, Huguet T. Insight on segregation distortions in two intraspecific crosses between annual species of Medicago (Leguminosae). Theor Appl Genet. 1997;94:682–91.

    Article  Google Scholar 

  71. 71.

    Majka J, Książczyk T, Kiełbowicz-Matuk A, Kopecký D, Kosmala A. Exploiting repetitive sequences and BAC clones in Festuca pratensis karyotyping. PLoS One. 2017;12(6):e0179043.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Gerlach WL, Dyer TA. Sequence organisation of the repeating units in the nucleus of wheat which contain 5S rRNA genes. Nucl Acids Res. 1980;8:4851–65.

    CAS  Article  Google Scholar 

  73. 73.

    Unfried I, Gruendler P. Nucleotide sequence of the 5.8S and 25S rRNA genes and the internal transcribed spacers from Arabidopsis thaliana. Nucl Acids Res. 1990;1:4011.

    Article  Google Scholar 

  74. 74.

    Thomas HM, Harper JA, Meredith MR, Morgan WG, King IP. Physical mapping of ribosomal DNA sites in Festuca arundinacea and related species by in situ hybridization. Genome. 1997;40:406–10.

    CAS  Article  Google Scholar 

  75. 75.

    Doyle JJ, Doyle JL. A rapid total DNA preparation procedure for fresh plant tissue. Focus. 1990;12:13–5.

    Google Scholar 

  76. 76.

    Amiryousefi A, Hyvönen J, Poczai P. iMEC: online marker efficiency calculator. Appl Plant Sci. 2018;6:e01159.

    Article  Google Scholar 

  77. 77.

    Canter PH, Pasakinskiene I, Jones RN, Humphreys MW. Chromosome substitutions and recombination in the amphiploid Lolium perenne × Festuca pratensis cv prior (2n=4x=28). Theor Appl Genet. 1999;98:809–14.

    Article  Google Scholar 

Download references


We thank Professor Arkadiusz Kosmala (Institute of Plant Genetics, Polish Academy of Sciences, Poznań, Poland) for valuable comments on the manuscript.


JM acknowledges scholarship from the National Science Center (Poland); Etiuda project (no. 2018/28/T/NZ9/00074).

Availability of data and materials

All data generated and analyzed during this study are included in this published article (and its supplementary information files).

Author information




JM and TK conceived and designed the experiments. JM performed the experiments. JM, KB, AJ, PK, HCK and TK analyzed the data. JM, KB and AJ drafted the manuscript. JM, PK and ZZ revised the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Joanna Majka.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. The scheme of crossing experiments. (TIF 1581 kb)

Additional file 2:

Table S1. Cytogenetic genotyping in F. pratensis, L. perenne, F. pratensis × L. perenne hybrids (F1) and in the following F2-F3 generations. (XLSX 85 kb)

Additional file 3:

Table S2. The Polymorphism Information Content (PIC) calculated for individual ISSR primer and the PIC mean value. PIC values are provided for individuals that were used as parents for each of the cross combinations, as well as for the F2 and F3 populations. (XLSX 11 kb)

Additional file 4:

Table S3. The results of χ2 test of ISSR bands segregation in combination no. 1 of F2 generation of F. pratensis × L. perenne hybrid. (XLSX 38 kb)

Additional file 5:

Table S4. The results of χ2 test of ISSR bands segregation in combination no. 2 of F2 generation of F. pratensis × L. perenne hybrid. (XLSX 39 kb)

Additional file 6:

Table S5. List of the most probable genotypes of P1, P2 and F1 parental plants for each of SSR locus. (XLSX 17 kb)

Additional file 7:

Table S6. The results of χ2 test of ISSR bands segregation in combination no. 3 of F3 generation of F. pratensis × L. perenne hybrid. (XLSX 45 kb)

Additional file 8:

Table S7. The results of χ2 test of ISSR bands segregation in combination no. 4 of F3 generation of F. pratensis × L. perenne hybrid. (XLSX 47 kb)

Additional file 9:

Table S8. List of the most probable genotypes of P1, P2, F1, and F2 parental plants for each of ISSR locus. (XLSX 18 kb)

Additional file 10:

Table S9. Specification of analyzed oligonucleotides. (XLSX 10 kb)

Additional file 11:

Table S10. ISSR band pattern in combinations no. 1 and 2 of F2 generation of F. pratensis × L. perenne hybrid. (XLSX 40 kb)

Additional file 12:

Table S11. ISSR band pattern in combinations no. 3 and 4 of F3 generation of F. pratensis × L. perenne hybrid. (XLSX 30 kb)

Additional file 13:

Table S12. Possible segregation ratios for the disomic and tetrasomic inheritance of ISSR loci depending on the possible genotypes and band pattern in parental species – F. pratensis and L. perenne and in the F2-F3 generations of F. pratensis × L. perenne hybrid. For simplicity, the F1 genotypes were described using one letter code (A or a), regardless of the origin of particular allele from Festuca or Lolium parent. (XLSX 10 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Majka, J., Bzdęga, K., Janiak, A. et al. Cytogenetic and molecular genotyping in the allotetraploid Festuca pratensis × Lolium perenne hybrids. BMC Genomics 20, 367 (2019).

Download citation


  • Festuca
  • Festulolium hybrids
  • Lolium
  • Genomic composition
  • Inheritance pattern
  • Rearrangements
  • Recombination