FISH landmarks reflecting meiotic recombination and structural alterations of chromosomes in wheat (Triticum aestivum L.)
BMC Plant Biology volume 21, Article number: 167 (2021)
DNA sequence composition affects meiotic recombination. However, the correlation between tandem repeat composition and meiotic recombination in common wheat (Triticum aestivum L.) is unclear.
Non-denaturing fluorescent in situ hybridization (ND-FISH) with oligonucleotide (oligo) probes derived from tandem repeats and single-copy FISH were used to investigate recombination in three kinds of the long arm of wheat 5A chromosome (5AL). 5AL535–18/275 arm carries the tandem repeats pTa-535, Oligo-18, and pTa-275, 5AL119.2–18/275 arm carries the tandem repeats pSc119.2, Oligo-18 and pTa-275, and 5AL119.2 arm carries the tandem repeats pSc119.2. In the progeny of 5AL535–18/275 × 5AL119.2, double recombination occurred between pSc119.2 and pTa-535 clusters (119–535 interval), and between pTa-535 and Oligo-18/pTa-275 clusters (535–18 interval). The recombination rate in the 119–535 interval in the progeny of 5AL535–18/275 × 5AL119.2–18/275 was higher than that in the progeny of 5AL535–18/275 × 5AL119.2. Recombination in the 119–535 interval produced 5AL119 + 535 segments with pTa-535 and pSc119.2 tandem repeats and 5ALNo segments without these repeats. The 5AL119 + 535 and 5ALNo segments were localized between the signal sites of the single-copy probes SC5A-479 and SC5A-527. The segment between SC5A-479 and SC5A-527 in the metaphase 5ALNo was significantly longer than that in the metaphase 5AL119 + 535.
The structural variations caused by tandem repeats might be one of the factors affecting meiotic recombination in wheat. Meiotic recombination aggregated two kinds of tandemly repeated clusters into the same chromosome, making the metaphase chromosome more condensed. To conclude, our study provides a robust tool to measure meiotic recombination and select parents for wheat breeding programs.
Increasing the chromosomal recombination rate can accelerate new cultivars’ development in common wheat (Triticum aestivum L.). Researchers have extensively studied the molecular mechanisms of meiotic recombination in plants and reported many genes involved in the process , which provide an opportunity to manipulate the mechanisms for crop improvement . Mutation in the anticrossover gene FANCM (Fanconi anemia complementation group M) resulted in a two-fold increase in meiotic recombination in hybrid rice and pea . However, the fancm mutation has almost no effect on recombination in hybrid Arabidopsis but resulted in a three-fold increase in inbreds [3,4,5]. Combining the anticrossover mutants recq4 (RecQ like helicase 4) and figl1 (FIDGETIN-like-1) resulted in a 7.8-fold increase in crossover frequency, while the fancm, recq4 and figl1 triple mutant displayed less recombination . These findings indicated that factors other than genes help control meiotic recombination. Both genes and chromatin structure control the recombination rate. In humans and animals, the recombination rate is related to cytogenetic structures that were displayed by staining intensity of G bands and sequence compositions including repetitive elements, GC content, CpG density and poly(A)/poly(T) stretches [6, 7]. In Arabidopsis thaliana, A-rich, CCN-repeat and CTT-repeat motifs are enriched in the crossover regions . Reduced DNA methylation at the CG sites increases recombination rate in the euchromatic, but not in the pericentric heterochromatin regions . Heterochromatin plays a role in meiotic recombination. It represses centromeric meiotic recombination in fission yeast . The role of pericentric heterochromatin in suppressing meiotic recombination has been widely studied in eukaryotes . The disruption of H3K9me2 and non-CG DNA methylation pathways via gene mutations increased pericentromeric crossovers in hybrid and inbred Arabidopsis . Meanwhile, in mice, a tandem array of mo-2 minisatellite conferred higher-order structures crucial for recombination in the pseudoautosomal region of sex chromosomes . Therefore, the effect of tandem repeats on chromosomal recombination should be investigated. Common wheat can be used as a reference model to study the role of tandem repeats in chromosomal recombination. Previous studies have indicated a correlation between the wheat chromosomal crossover and gene-rich regions [14,15,16,17,18]. However, different aspects of recombination should be considered for common wheat because of allopolyploidy and repetitive sequences .
Some new tandem repeats were discovered from common wheat [20, 21]. In our previous work, oligonucleotide (oligo) probes derived from these tandem repeats displayed large structural variations in the 5AL arms of common wheat [22, 23]. However, the correlation between the composition of tandem repeats and the rate of recombination in wheat is still unclear. This study investigates the effect of genomic structural variations reflected by tandem repeats on recombination frequency in the wheat chromosome 5A.
FISH karyotypes of 5A chromosomes
In wheat, the 5A chromosomes can be distinguished from the other chromosomes based on the signal patterns of Oligo-713, Oligo-pSc119.2–1, and Oligo-pTa535–1 [22, 24]. The Oligo-pSc119.2–1 signals on the telomeric region of 5BS (the short arm of 5B chromosome) are stronger than that on 5AS (the short arm of 5A chromosome), and the signal patterns of Oligo-pTa535–1 on 5D chromosomes are different from that on 5A chromosomes [22, 24]. However, when Oligo-pSc119.2–1 signals occur on both 5AS and 5AL, it is challenging to distinguish 5A and 2B chromosomes. Then, the 5A chromosomes are identified based on the submetacentric feature. The pericentromeric region of 7AS (the short arm of the 7A chromosome) and 7DS (the short arm of the 7D chromosome) contains Oligo-713 signals, however, the 7A and 7D chromosomes can be distinguished from 5A chromosomes based on the signal patterns of Oligo-pSc119.2–1 and Oligo-pTa535–1 .
The non-denaturing fluorescence in situ hybridization (ND-FISH) of this study indicated that the signals of Oligo-713 occurred on the pericentromeric region of 5AS arms of all the nine wheat cultivars/lines (Fig. 1). The probe Oligo-pSc119.2–1 produced signals on 5AS arms of eight wheat cultivars/lines, except for CM90, and on 5AL arms of 14 T141–2, XKM8, M1403, and MM43 (Fig. 1). The signals of Oligo-pTa535–1 occurred on 5AL arms of 14 T105–1, CM36, and CM39 (Fig. 1). Both Oligo-pSc119.2–1 and Oligo-pTa535–1 probes produced signals on 5AL arms of KCM2 (Fig. 1). The signals of Oligo-275.1 and Oligo-18 occurred on both 5AS and 5AL arms of 14 T105–1, 14 T141–2, CM90, and KCM2 (Fig. 1), and 5AL arms of CM39, XKM8 and CM36 (Fig. 1). Although Oligo-275.1 and Oligo-18 signals appeared on 5AS arms of M1403 and MM43, they did not occur on 5AL arms of the two cultivars (Fig. 1). Therefore, we named the 5AL arms as follows: 14 T105–1, CM39, and CM36 as 5AL535–18/275; 14 T141–2 and XKM8 as 5AL119.2–18/275; M1403 and MM43 as 5AL119.2; CM90 as 5AL18/275; and KCM2 as 5AL119.2–535-18/275 (Fig. 1).
FISH karyotypes and recombination rate of the F2 generation
A total of 420 F2 generation seeds (840 5A chromosomes) derived from the CM39 × M1403, CM39 × MM43, CM39 × XKM8, and 14 T105–1 × 14 T141–2 hybrid combinations were analyzed to determine the recombination rate of 5AL arms. These included 106 seeds (212 5A chromosomes) derived from CM39 × M1403 (5AL535–18/275 × 5AL119.2), 110 (220 5A chromosomes) from CM39 × MM43 (5AL535–18/275 × 5AL119.2), 100 (200 5A chromosomes) from CM39 × XKM8 (5AL535–18/275 × 5AL119.2–18/275) and 104 (208 5A chromosomes) from 14 T105–1 × 14 T141–2 (5AL535–18/275 × 5AL119.2–18/275). For the progeny derived from CM39 × M1403 and CM39 × MM43, two recombination intervals were detected (Fig. 2 and Additional files 1, 2). One recombination occurred between the signal sites of Oligo-pSc119.2–1 and Oligo-pTa535–1, and this region was named 119–535 interval. The other recombination occurred between the signal sites of Oligo-pTa535–1 and Oligo-18/Oligo-275.1, and this region was named 535–18 interval (Fig. 2 and Additional files 1, 2). For the progeny of CM39 × XKM8 and 14 T105–1 × 14 T141–2, only the recombination occurred in the 119–535 interval was observed because the Oligo-275.1 and Oligo-18 signals were detected on both 5AL535–18/275 and 5AL119.2–18/275 arms of (Fig. 3 and Additional file 3). The 106 seeds of CM39 × M1403 had 22 combination types of 5A chromosomes (Fig. 2), and the recombination occurred in 119–535 and 535–18 intervals on 20 and 88 5A chromosomes, respectively, including the seven chromosomes with double recombination (Fig. 2). The 110 seeds of CM39 × MM43 had 25 combination types of 5A chromosomes (Fig. 2). Of the 220 5A chromosomes, recombination occurred in 119–535 and 535–18 intervals on 31 and 96 5A chromosomes, respectively, including the 15 chromosomes with double recombination (Fig. 2). Only eight combination types of 5A chromosomes were observed among the 100 seeds (200 chromosomes) derived from CM39 × XKM8, and recombination in the 119–535 interval was observed on 37 5A chromosomes (Fig. 3). The signals of Oligo-713, Oligo-275.1, and Oligo-18 were not shown on these 5A chromosomes because they could not reflect recombination in the 535–18 interval. Similarly, only 10 combination types of 5A chromosomes were observed among the 104 seeds (208 5A chromosomes) derived from 14 T105–1 × 14 T141–2, and recombination in the 119–535 interval occurred on 42 5A chromosomes (Fig. 3).
Sequential single-copy FISH and ND-FISH assays indicated that the signal sites of Oligo-pSc119.2–1 and SC5A-479 were close each other (Fig. 4a, c and Additional file 4c, f) and that of Oligo-pTa535–1 and SC5A-527 were close each other (Fig. 4b, c and Additional file 4c, f). The signal sites of SC5A-479 and SC5A-527 on 5AL119.2–535-18/275 arms of KCM2 confirmed this observation (Fig. 4e and Additional file 4i). Additionally, the signals of Oligo-18 and SC5A-586 were close (Fig. 4f and Additional file 4 l) (The details of obtaining the three single copy probes are described in Methods section.). The 119–535 and 535–18 intervals were determined as 48 Mb and 59 Mb long, respectively, based on the physical distance of the single-copy probes SC5A-479, SC5A-527, and SC5A-586. The recombination rate of the 5AL arms in the four kinds of F2 populations is shown in Table 1. The recombination rate in 119–535 interval in the progeny derived from CM39 × XKM8 (5AL535–18/275 × 5AL119.2–18/275) and 14 T105–1 × 14 T141–2 (5AL535–18/275 × 5AL119.2–18/275) was higher than that in the progeny derived from the other two combinations (5AL535–18/275 × 5AL119.2) (Table 1), and a significant difference was observed between the progeny of 5AL535–18/275 × 5AL119.2–18/275 and 5AL535–18/275 × 5AL119.2 (Fig. 5a). However, no significant difference was observed for the recombination rate in 535–18 interval between the progeny of CM39 × M1403 (5AL535–18/275 × 5AL119.2) and CM39 × MM43 (5AL535–18/275 × 5AL119.2) (Fig. 5b). Our analysis indicated that crossover interference existed between 119 and 535 and 535–18 intervals (Table 1), however, no significant difference was observed for the double recombination rate between the progeny of CM39 × M1403 (5AL535–18/275 × 5AL119.2) and CM39 × MM43 (5AL535–18/275 × 5AL119.2) (Fig. 5c) (The details of the calculating of recombination and interference are described in Methods section.).
Relative length of the metaphase chromosome segment between SC5A-479 and SC5A-527
Figures 2 and 3 show that the recombination in the 119–535 interval formed the recombination types Rec3, Rec4 and DRec2. The signals of both Oligo-pSc119.2–1 and Oligo-pTa535–1 occurred on type Rec4 5AL arms that disappeared from type Rec3 5AL arms (Figs. 2, 3 and Additional files 1, 2, 3). The 119–535 interval corresponds to the segment between the signal sites of SC5A-479 and SC5A-527 on metaphase 5AL arms (Fig. 4 and Additional file 4). We named the segment between SC5A-479 and SC5A-527 on 5AL119.2, 5AL535–18/275, Rec3, and Rec4 5AL arms as 5AL119, 5AL535, 5ALNo, and 5AL119 + 535, respectively. The relative metaphase length (RML) of the four kinds of metaphase 5AL segments (5AL119, 5AL535, 5AL119 + 535, and 5ALNo) was measured (Fig. 6a). From the four types of F2 generations, a total of 44, 36, 43 and 40 segments of 5AL119, 5AL535, 5AL119 + 535 and 5ALNo, respectively, were successfully measured (Fig. 6b). The 5ALNo segment was significantly longer than 5AL119 + 535 (Fig. 6b). Further, the 5AL119, 5AL535, 5AL119 + 535, and 5ALNo segments (Fig. 1) in the wheat cultivars XKM8, CM36, KCM2 and CM90 were also measured to confirm this (Fig. 6c). A total of 80, 84, 82, and 86 5AL119, 5AL535, 5AL119 + 535, and 5ALNo segments were obtained and measured from the four wheat cultivars (Fig. 6c). The 5AL119, 5AL535, and 5ALNo segments were significant longer than 5AL119 + 535 (Fig. 6c).
Correlation between tandem repeats and meiotic recombination
The four kinds of tandem repeats pSc119.2, pTa-535, pTa-275 and Oligo-18 displayed apparent structural variations in the 5AL arms of wheat. In this study, the cross combinations CM39 × XKM8, CM39 × M1403 and CM39 × MM43 have the same parent CM39. However, the recombination rate in the 119–535 interval was lower in the progeny derived from CM39 × M1403 (5AL535–18/275 × 5AL119.2) and CM39 × MM43 (5AL535–18/275 × 5AL119.2) than that in the progeny derived from CM39 × XKM8 (5AL535–18/275 × 5AL119.2–18/275). Besides, the 5AL arms of M1403 and MM43 showed ND-FISH signal patterns different from that of XKM8. It was speculated that recombination rate variation in the 119–535 interval is consistent with the ND-FISH signal patterns on 5AL arms, supported by the progeny derived from 14 T105–1 × 14 T141–2 (5AL535–18/275 × 5AL119.2–18/275). Therefore, a correlation between the composition of tandem repeats in 5AL arms and meiotic recombination might existed.
Researchers have observed a higher rate of motifs of DNA transposons in the recombination intervals. In potato, the Stowaway family of miniature inverted-repeat transposable elements spanned the crossover regions . Darrier et al. observed a higher frequency of a DNA motif specific to the TIR-Mariner DNA transposon in common wheat recombinant intervals . These studies found differences in DNA transposon composition between recombination and non-recombination hotspots along the same chromosome, a phenomenon observed in diverse populations [16, 25]. Putative sequence divergence or insertions/deletions in wheat led to significant differences in crossover frequency along 3B chromosomes between two different F2 segregating populations . Meanwhile, a nested association mapping (NAM) population used to map the QTL affecting the crossover distribution and frequency indicated a lower recombination rate in the regions with more single nucleotide polymorphisms (SNPs) than with fewer SNPs in wheat . Besides, similar crossover patterns were detected in populations derived from closely related parents . Even in the mutants of anticrossover genes, recombination was prevented in regions with the highest sequence divergence, displayed by SNP . In Arabidopsis, the msh2 (MutS-related heterodimers) mutant displayed significantly reduced SNP enrichment around crossovers compared with the wild type . These studies indicated that SNP polymorphisms and indels are important factors affecting meiotic recombination. Although a significant difference was observed in the recombination rate in the 119–535 interval between the progeny derived from 5AL535–18/275 × 5AL119.2 and 5AL535–18/275 × 5AL119.2–18/275 (Fig. 5a), other sequence differences among the 5AL arms, especially in 119–535 and 535–18 intervals, that result in different genetic distances cannot be ignored. In addition to the tandem repeats, these sequence components may also affect the recombination rate between the different 5AL arms. Meanwhile, the recombination rate in the 119–535 interval in the progeny derived from 5AL535–18/275 × 5AL119.2–18/275 was higher than that in the progeny derived from 5AL535–18/275 × 5AL119.2. Therefore, the effects of tandem repeats on meiotic recombination should be further investigated.
Effects of juxtaposed heterozygous and homozygous regions on recombination
Although it cannot be concluded that the composition of tandem repeats alone resulted in a comparatively low recombination rate in the 119–535 interval in the progeny of CM39 × M1403 (5AL535–18/275 × 5AL119.2) and CM39 × MM43 (5AL535–18/275 × 5AL119.2), the variations in the recombination rate in this interval are consistent with the structural differences in the 5AL arms. Different ND-FISH patterns occurred at the pSc119.2 and pTa-535 sites between 5AL535–18/275 and 5AL119.2–18/275. However, the different ND-FISH patterns were observed between 5AL535–18/275 and 5AL119.2 at the pSc119.2, pTa-535, and Oligo-18/pTa-275 sites. In Arabidopsis, heterozygous regions increase crossover in juxtaposed megabase homozygous and heterozygous regions, and heterozygosity can increase crossover interference . Thus, our findings, together with these previous reports, indicate the need to investigate the effects of juxtaposition of heterozygous and homozygous intervals on meiotic recombination.
Effects of tandem repeats on metaphase chromosome condensation
In this study, recombination formed the 5AL119 + 535 and 5ALNo metaphase segments. A significant difference in length was observed only between 5AL119 + 535 and 5ALNo segments derived from the F2 populations; 5AL119 + 535 of KCM2 was significantly shorter than that of XKM8 (5AL119), CM36 (5AL535), and CM90 (5ALNo) (Fig. 6b, c). These observations indicate other sequence differences in the 119–535 interval of the 5AL arms, affecting the length of 5AL119, 5AL535, 5AL119 + 535 and 5ALNo segments. However, for both F2 populations and wheat cultivars, the 5ALNo segments were significantly longer than the 5AL119 + 535 segments. These findings indicate that the metaphase chromosome segments containing two tandemly repeated clusters close to each other in the same chromosome will be more condensed. Conversely, the metaphase chromosome segments with fewer tandem repeats will be more relaxed. Moreover, the cellSens Dimension software-based method used in this study to measure the length is reliable. Tandemly repeated satellite DNA sequences can drive population and species divergence by inducing alterations in heterochromatin and/or centromere . Satellite DNA sequences are the major heterochromatin components that play an essential role in heterochromatin formation and regulation . In mouse, the heterochromatin alters the loop size of the chromatin of mitotic chromosome . Therefore, the 5AL119 + 535 segment of this study might be highly heterochromatinized by the aggregating pSc119.2 and pTa-535 making the segment more condensed.
Application of ND-FISH in wheat breeding programs
Understanding the factors affecting the variations in recombination rate and recombination point is important in wheat breeding programs. Generally, crossover frequency decreases from telomere to centromere in eukaryotes , including wheat [15, 16, 32]. Mutations in the anticrossover genes do not increase the recombination in regions close to centromere . However, increasing the interstitial recombination rate reduces deleterious genetic load, which will benefit crop improvement . Therefore, we should understand the mechanisms and methods to increase the recombination rate in the proximal regions . The differences in recombination rate detected in this study reflect differences in sequence composition in 119–535 and 535–18 intervals derived from different crosses. ND-FISH assay displays the differences in tandem repeat composition in wheat chromosomes and can be used to predict the recombination rate.
Earlier, Derrier et al. used chromosome 3B pseudomolecule and high-throughput SNP detection to map 252 crossover events at intervals of < 26 kb, a high-resolution crossover location in common wheat . Although high-throughput sequencing and DNA markers have great advantages in studying chromosomal recombination [15,16,17,18, 32], large number of individual recombinants need to be sequenced. Cytological markers have advantages of visualization and intuition in studying chromosomal recombination [33,34,35]. ND-FISH technology based on oligo probes can be used to establish abundant FISH karyotypes of wheat chromosomes conveniently [21, 36], and oligo probes derived from tandem repeats can display structural variations in wheat chromosomes [22, 23, 36, 37]. Besides, chromosomal recombination indicated by FISH karyotypes is visual comprehend and can be used conveniently in wheat breeding programs [34, 35]. Therefore, high-throughput sequencing and DNA markers should be combined with FISH karyotypes to study the wheat chromosomal recombination rules.
The present study confirms the effects of tandem repeats on meiotic recombination in wheat, however, studies should be carried out to identify other factors determining the recombination rate in 5AL arms. Our study provides a robust visual tool based on ND-FISH to measure meiotic recombination and crossover interference in wheat.
Four wheat cultivars Chuanmai 39 (CM39), Xikemai 8 (XKM8), Mian 1403 (M1403) and Mianmai 43 (MM43), and two wheat lines 14 T105–1 and 14 T141–2 were used as parents for hybridization. Four cross combinations, CM39 × M1403, CM39 × MM43, CM39 × XKM8, and 14 T105–1 × 14 T141–2, were carried out, and four F2 generations were obtained by selfing. The wheat cultivars Chuanmai 36 (CM36), Chuanmai 90 (CM90), and Kechengmai 2 (KCM2) were used to measure the distance between the FISH signals of the single-copy probes SC5A-479 and SC5A-527 on the 5AL arms at metaphase.
Non-denaturing FISH (ND-FISH)
Oligo-pSc119.2–1, Oligo-pTa535–1 , Oligo-713, Oligo-275.1 , and Oligo-18  were used as oligo probes for ND-FISH. The information on these oligo probes is listed in Additional file 5. These oligo probes can replace their original sequences’ roles in identifying wheat chromosomes [22, 24]. The sequence of the Oligo-18 probe is identical to the original one, both are 18 bp long . The metaphase chromosomes were prepared from the root tips according to the method described by Han et al. . ND-FISH was performed according to the method described by Tang et al. . An epifluorescence microscope (BX51, Olympus Corporation, Tokyo, Japan) with cellSens Dimension software (Olympus Corporation, Tokyo, Japan) was used to capture images.
The Tandem Repeat Finder (TRF, Version 4.09)  program was used to filter the tandem repeats in the 5A chromosome (IWGSC RefSeq Version 2.0). An in-house R package was then used to filter other repeated sequences, including retrotransposon and transposon elements. The remaining sequences were used as the query to align with the full-length sequence of chromosome 5A of the bread wheat variety Chinese Spring (IWGSC RefSeq Version 2.0) using the in-house R package. The query sequences aligned to themselves at a single site were kept as preparatory probes. The sequences of Oligo-pSc119.2–1, Oligo-pTa535–1 , and Oligo-18  were also used as the query to align with the full-length sequence of 5A chromosome using BLAST in the B2DSC web server (http://mcgb.uestc.edu.cn/b2dsc) . The sequences of Oligo-pTa535–1 and Oligo-pSc119.2–1 hit their targets with high copy numbers at 431–432 Mbp and 507–508 Mbp sites on the 5AL arm, respectively. These were used to narrow down the options for single copy probes. Finally, two single-copy sequences, SC5A-479 (479272790–479,275,822 bp) and SC5A-527 (527288389–527,289,202 bp), were selected as probes. The sequence of Oligo-18 hit the target with a high copy number at one site (584–585 Mbp) on the 5AL arm. Then a single copy sequence SC5A-586 (586379936–586,380,903 bp) was selected as the probe. The primer pairs of SC5A-479 (5’TCGTTGACTAGAAAGACGTG TGT3’, 5’ACGCCTGTGTTAAGTTAAGTGAC3’), SC5A-527 (5’TGCGTACATAGGGTGAGTGTATG3’, 5’GGCCTCTGGAAGAACGTTTTAT3’), and SC5A-586 (5’TTGCTCGTGTCCACCATTGA3’, 5’TGTGGAATACTTACCGCGCA3’) were used to amplify the three single-copy probe sequences. The target sequences were cloned into the TSINGKE pClone007 vector (TSINGKE, Chengdu, China) and labeled with Texas-Red-5-dUTP (PerkinElmer, USA) according to the method described by Han et al. . The root tip metaphase chromosomes were prepared, and the hybridization was performed as described by Han et al.  with slight modifications. The probe mixture contained 5–6 ng/μL of each probe and 1 × ENZO buffer (ENZO Life Science Inc. USA). Slides were washed in 2 × SSC buffer containing 0.1% NP-40 detergent (Solarbio Life Sciences Ltd. China) for 3–5 min at 45–50 °C. An epifluorescence microscope (BX51, Olympus Corporation, Tokyo, Japan) with cellSens Dimension software (Olympus Corporation, Tokyo, Japan) was used to capture the images.
Calculation of recombination rate and crossover interference
Recombination was determined based on the FISH patterns of the probes Oligo-pSc119.2–1, Oligo-pTa535–1, Oligo-275.1, and Oligo-18. The recombination rate was represented in cM/Mb. In this study, recombination was observed in two intervals: the 119–535 interval between Oligo-pSc119.2–1 and Oligo-pTa535–1 signal sites and the 535–18 interval between Oligo-pTa535–1 and Oligo-18/Oligo-275.1 signal sites. The recombination rate and crossover interference were calculated according to the methods described by Ziolkowski et al.  using the following formulas: 119–535 cM = number of 5AL arms with recombination in 119–535 interval/total number of 5AL arms × 100, 535–18 cM = number of 5AL arms with recombination in 535–18 interval/total number of 5AL arms × 100, Recombination rate in 119–535 interval = 119–535 cM/48 Mb × 100, Recombination rate in 535–18 interval = 535–18 cM/59 Mb × 100, Observed double crossover (Obs DCO) = number of 5AL arms with recombination in both 119–535 and 535–18 intervals/total number of 5AL arms × 100, Expected double crossover (Exp DCO) = (119–535 cM/100) × (535–18 cM/100) × total number of 5AL arms, Coefficient of Coincidence (CoC) = Obs DCO/Exp DCO, and Crossover interference = 1 − CoC. Further, t-test and Chi-square test were carried out to determine significant differences in recombination rate and crossover interference, respectively, as described by Girard et al. . An interference value close to zero indicates interference, and value close to 1, indicates no interference.
Measurement of length between SC5A-479 and SC5A-527 probes on metaphase chromosomes
The length of metaphase 5AL segment (MSL) between the signals of the two single-copy probes SC5A-479 and SC5A-527 was measured. Four kinds of metaphase 5AL segments were measured: segments with Oligo-pSc119.2–1 signal (named 5AL119), segments with Oligo-pTa535–1 signal (named 5AL535), segments with signals of the two probes (named 5AL119 + 535), and segments without the signals of the two probes (named 5ALNo). These four kinds of segments were derived from the four F2 generations, wheat cultivars CM36, CM90, XKM8, and KCM2. The length of the metaphase 5AL segment between SC5A-479 and SC5A-527 (MSL) was determined as relative metaphase length (RML) to avoid errors caused by chromosome condensation. RML = MSL/WL where WL represents the whole length of the 5AL arm from the central point of the centromere to the distal end. The MSL and WL were measured using cellSens Dimension software (Olympus Corporation, Tokyo, Japan). A one-way ANOVA was carried out using GraphPad Prism software (Version 5) for pairwise comparisons among the groups. Graphs were also plotted using GraphPad Prism software.
Availability of data and materials
The materials used and/or analyzed in the current study are available from the corresponding author on reasonable request.
Non-denaturing fluorescence in situ hybridization
- Obs DCO:
Observed double crossover
- Exp DCO:
Expected double crossover
Coefficient of coincidence
Length of metaphase 5AL segment
Relative metaphase length
Whole length of 5AL arm
Wang Y, Copenhaver GP. Meiotic recombination: mixing it up in plants. Annu Rev Plant Biol. 2018;69(1):577–609. https://doi.org/10.1146/annurev-arplant-042817-040431.
Blary A, Jenczewski. Manipulation of crossover frequency and distribution for plant breeding. Theor Appl Genet. 2019;132(3):575–92. https://doi.org/10.1007/s00122-018-3240-1.
Mieulet D, Aubert G, Bres C, Klein A, Droc G, Vieille E, et al. Unleashing meiotic crossovers in crops. Nat Plants. 2018;4(12):1010–6. https://doi.org/10.1038/s41477-018-0311-x.
Girard C, Chelysheva L, Choinard S, Froger N, Macaisne N, Lemhemdi A, et al. AAA-ATPase FIDGETIN-LIKE 1 and helicase FANCM antagonize meiotic crossovers by distinct mechanisms. PLoS Genet. 2015;11(7):e1005369. https://doi.org/10.1371/journal.pgen.1005369.
Fernandes JB, Séguéla-Arnaud M, Larchevêque C, Lloyd AH, Mercier R. Unleashing meiotic crossovers in hybrid plants. Proc Natl Acad Sci U S A. 2018;115(10):2431–6. https://doi.org/10.1073/pnas.1713078114.
Kong A, Gudbjartsson DF, Sainz J, Jonsdottir GM, Gudjonsson SA, Richardsson B, et al. A high-resolution recombination map of the human genome. Nat Genet. 2002;31(3):241–7. https://doi.org/10.1038/ng917.
Jensen-Seaman MI, Furey TS, Payseur BA, Lu Y, Roskin KM, Chen CF, et al. Comparative recombination rates in the rat, mouse, and human genomes. Genome Res. 2004;14(4):528–38. https://doi.org/10.1101/gr.1970304.
Shilo S, Melamed-Bessudo C, Dorone Y, Barkai N, Levy AA. DNA crossover motifs associated with epigenetic modifications delineate open chromatin regions in Arabidopsis. Plant Cell. 2015;27(9):2427–36. https://doi.org/10.1105/tpc.15.00391.
Melamed-Bessudo C, Levy AA. Deficiency in DNA methylation increases meiotic crossover rates in euchromatic but not in heterochromatic regions in Arabidopsis. Proc Natl Acad Sci U S A. 2012;109(16):E981–8. https://doi.org/10.1073/pnas.1120742109.
Ellermeier C, Higuchi EC, Phadnis N, Holm L, Geelhood JL, Thon G, et al. RNAi and heterochromatin repress centromeric meiotic recombination. Proc Natl Acad Sci U S A. 2010;107(19):8701–5. https://doi.org/10.1073/pnas.0914160107.
Nambiar M, Smith GR. Repression of harmful meiotic recombination in centromeric regions. Semin Cell Dev Biol. 2016;54:188–97. https://doi.org/10.1016/j.semcdb.2016.01.042.
Underwood CJ, Choi K, Lambing C, Zhao X, Serra H, Borges F, et al. Epigenetic activation of meiotic recombination near Arabidopsis thaliana centromeres via loss of H3K9me2 and non-CG DNA methylation. Genome Res. 2018;28(4):519–31. https://doi.org/10.1101/gr.227116.117.
Acquaviva L, Boekhout M, Karasu ME, Brick K, Pratto F, Li T, et al. Ensuring meiotic DNA break formation in the mouse pseudoautosomal region. Nature. 2020;582(7812):426–31. https://doi.org/10.1038/s41586-020-2327-4.
Erayman M, Sandhu D, Sidhu D, Dilbirligi M, Baenziger PS, Gill KS. Demarcating the gene-rich regions of the wheat genome. Nucleic Acids Res. 2004;32(12):3546–65. https://doi.org/10.1093/nar/gkh639.
Saintenac C, Faure S, Remay A, Choulet F, Ravel C, Paux E, et al. Variation in crossover rates across a 3-Mb contig of bread wheat (Triticum aestivum) reveals the presence of a meiotic recombination hotspot. Chromosoma. 2011;120(2):185–98. https://doi.org/10.1007/s00412-010-0302-9.
Darrier B, Rimbert H, Balfourier F, Pingault L, Josselin AA, Servin B, et al. High-resolution mapping of crossover events in the hexaploid wheat genome suggests a universal recombination mechanism. Genetics. 2017;206(3):1373–88. https://doi.org/10.1534/genetics.116.196014.
Jordan KW, Wang S, He F, Chao S, Lun Y, Paux E, et al. The genetic architecture of genome-wide recombination rate variation in allopolyploid wheat revealed by nested association mapping. Plant J. 2018;95(6):1039–54. https://doi.org/10.1111/tpj.14009.
Gardiner LJ, Wingen LU, Bailey P, Joynson R, Brabbs T, Wright J, et al. Analysis of the recombination landscape of hexaploid bread wheat reveals genes controlling recombination and gene conversion frequency. Genome Biol. 2019;20(1):69. https://doi.org/10.1186/s13059-019-1675-6.
Lambing C, Franklin CH, Wang CJR. Understanding and manipulating meiotic recombination in plants. Plant Physiol. 2017;173(3):1530–42. https://doi.org/10.1104/pp.16.01530.
Komuro S, Endo R, Shikata K, Kato A. Genomic and chromosomal distribution patterns of various repeated DNA sequences in wheat revealed by a fluorescence in situ hybridization procedure. Genome. 2013;56(3):131–7. https://doi.org/10.1139/gen-2013-0003.
Tang S, Tang Z, Qiu L, Yang Z, Li G, Lang T, et al. Developing new oligo probes to distinguish specific chromosomal segments and the a, B, D genomes of wheat (Triticum aestivum L.) using ND-FISH. Front. Plant Sci. 2018;9:1104.
Tang S, Qiu L, Xiao Z, Fu S, Tang Z. New oligonucleotide probes for ND-FISH analysis to identify barley chromosomes and to investigate polymorphisms of wheat chromosomes. Genes (Basel). 2016;7(12):118. https://doi.org/10.3390/genes7120118.
Jiang M, Xiao ZQ, Fu SL, Tang ZX. FISH karyotype of 85 common wheat cultivars/lines displayed by ND-FISH using oligonucleotide probes. Cereal Res Commun. 2017;45(4):549–63.
Tang ZX, Yang ZJ, Fu SL. Oligonucleotides replacing the roles of repetitive sequences pAs1, pSc119.2, pTa-535, pTa71, CCS1, and pAWRC.1 for FISH analysis. J Appl Genet. 2014;55(3):313–8. https://doi.org/10.1007/s13353-014-0215-z.
Marand AP, Jansky SH, Zhao H, Leisner CP, Zhu X, Zeng Z, et al. Meiotic crossovers are associated with open chromatin and enriched with Stowaway transposons in potato. Genome Biol. 2017;18(1):203. https://doi.org/10.1186/s13059-017-1326-8.
Blackwell AR, Dluzewska J, Szymanska-Lejman M, Desjardins S, Tock AJ, Kbiri N, et al. MSH2 shapes the meiotic crossover landscape in relation to interhomolog polymorphism in Arabidopsis. EMBO J. 2020;39(21):e104858. https://doi.org/10.15252/embj.2020104858.
Ziolkowski PA, Berchowitz LE, Lambing C, Yelina NE, Zhao X, Kelly KA, et al. Juxtaposition of heterozygous and homozygous regions causes reciprocal crossover remodelling via interference during Arabidopsis meiosis. eLife. 2015;4:e03708. https://doi.org/10.7554/eLife.03708.
Feliciello I, Akrap I, Brajković J, Zlatar I, Ugarković Đ. Satellite DNA as a driver of population divergence in the red flour beetle Tribolium castaneum. Genome Biol Evol. 2014;71(1):228–39.
Pezer Ž, Brajković J, Feliciello I, Ugarkovć Đ. Satellite DNA-mediated effects on genome regulation. Genome Dyn. 2012;7:153–69. https://doi.org/10.1159/000337116.
Fitz-James MH, Tong P, Pidoux AL, Ozadam H, Yang L, White SA, et al. Large domains of heterochromatin direct the formation of short mitotic chromosome loops. Elife. 2020;9:e57212. https://doi.org/10.7554/eLife.57212.
Haenel Q, Laurentino TG, Roesti M, Berner D. Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics. Mol Ecol. 2018;27(11):2477–97. https://doi.org/10.1111/mec.14699.
Saintenac C, Falque M, Martin OC, Paux E, Feuillet C, Sourdille P. Detailed recombination studies along chromosome 3B provide new insights on crossover distribution in wheat (Triticum aestivum L.). Genetics. 2009;181(2):393–403. https://doi.org/10.1534/genetics.108.097469.
Lukaszewski AJ, Curtis CA. Physical distribution of recombination in B-genome chromosomes of tetraploid wheat. Theor Appl Genet. 1993;86(1):121–7. https://doi.org/10.1007/BF00223816.
Du H, Li S, He Y, Tang Z, Fu S, Ren Y. Structural variation of 6VS/6AL translocation chromosomes in Mianmai 37 and its derivatives. J Triticeae Crops. 2019;39(6):659–65.
Xi W, Jiang J, Tang Z, Wang S. Structural variation of chromosomes of Nanmai wheat varieties. J Triticeae Crops. 2019;39(5):567–73.
Huang X, Zhu M, Zhuang L, Zhang S, Wang J, Chen X, et al. Structural chromosome rearrangements and polymorphisms identified in Chinese wheat cultivars by high resolution multiplex oligonucleotide FISH. Theor Appl Genet. 2018;131(9):1967–86. https://doi.org/10.1007/s00122-018-3126-2.
Guo J, Gao D, Gong W, Li H, Li J, Li G, et al. Genetic diversity in common wheat lines revealed by fluorescence in situ hybridization. Plant Syst Evol. 2019;305(3):247–54. https://doi.org/10.1007/s00606-019-1567-y.
Han FP, Lamb JC, Birchler A. High frequency of centromere inactivation resulting in stable dicentric chromosomes of maize. Proc Natl Acad Sci U S A. 2006;103(9):3238–43. https://doi.org/10.1073/pnas.0509650103.
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80. https://doi.org/10.1093/nar/27.2.573.
Lang T, Li G, Wang H, Yu Z, Chen Q, Yang E, et al. Physical location of tandem repeats in the wheat genome and application for chromosome identification. Planta. 2019;249(3):663–75. https://doi.org/10.1007/s00425-018-3033-4.
We are thankful to Ennian Yang (Crop Research Institute, Sichuan Academy of Agricultural Sciences, China) and Yong Ren (Mianyang Branch of National Wheat Improvement Center, Mianyang Institute of Agricultural Sciences, China) for their kindly providing seeds of CM39, CM36, CM90, KCM2, M1403, MM43 and XKM8.
This manuscript is provided by the Science and Technology Project of Sichuan, China (No. 2020YJ0128). The funding body has no role in study design, data collection, analysis, and interpretation, and manuscript writing.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The recombination types of 5AL arms in the F2 progeny derived from CM39 × M1403. ‘5A-Rec1’ and ‘5A-Rec2’ indicate 5AL arms formed by recombination in the 535–18 interval. ‘5A-Rec3’ and ‘5A-Rec4’ indicate the 5AL arms formed by recombination in the 119–535 interval. ‘5A-DRec1’ and ‘5A-DRec2’ indicate the 5AL arms formed by recombination in both 535–18 and 119–535 intervals. Chromosomes were counterstained with DAPI (blue). Scale bar: 10 μm.
The recombination types of 5AL arms in the F2 progeny derived from CM39 × MM43. ‘5A-Rec1’ and ‘5A-Rec2’ indicate 5AL arms formed by recombination in the 535–18 interval. ‘5A-Rec3’ and ‘5A-Rec4’ indicate the 5AL arms formed by recombination in the 119–535 interval. ‘5A-DRec1’ and ‘5A-DRec2’ indicate the 5AL arms formed by recombination in both 535–18 and 119–535 intervals. ‘5A’ indicates the 5A chromosome derived from CM39. Chromosomes were counterstained with DAPI (blue). Scale bar: 10 μm.
The recombination types of 5AL arms in the F2 progeny derived from CM39 × XKM8 and 14 T105–1 × 14 T141–2. ‘5A-Rec3’ and ‘5A-Rec4’ indicate the 5AL arms formed by recombination in the 119–535 interval. ‘5A’ in (a) and (d) indicates the 5A chromosomes derived from XKM8 and 14 T141–2. ‘5A’ in (b) and (c) indicates the 5A chromosomes derived from CM39 and 14 T105–1. Chromosomes were counterstained with DAPI (blue). Scale bar: 10 μm.
Sequential single-copy FISH and ND-FISH assays to determine the position of the signals of SC5A-479, SC5A-527, SC5A-586, Oligo-pSc119.2–1, Oligo-pTa535–1 and Oligo-18 probes. (a, b, c) Cells derived from CM39 × MM43 showing the signal sites of SC5A479 and SC5A527 close to those of Oligo-pSc119.2–1 and Oligo-pTa535–1, respectively, and 5AL119 and 5AL535 segments corresponding to the 119–535 interval between SC5A479 and SC5A527 signal sites. (d, e, f) Cells derived from CM39 × MM1403 showing 5AL119 + 535 and 5ALNo segments corresponding to the 119–535 interval between SC5A479 and SC5A527 signal sites. (g, h, i) Cells of KCM2 showing 5AL119 + 535 segment between SC5A479 and SC5A527 sites. (j, k, l) Cells derived from CM39 × MM43 showing the signal sites of Oligo-18 and SC5A-586 close to each other. Chromosomes were counterstained with DAPI (blue). Scale bar: 10 μm.
Oligonucleotide probes used for ND-FISH assays in this study.
About this article
Cite this article
Zou, Y., Wan, L., Luo, J. et al. FISH landmarks reflecting meiotic recombination and structural alterations of chromosomes in wheat (Triticum aestivum L.). BMC Plant Biol 21, 167 (2021). https://doi.org/10.1186/s12870-021-02947-1