# Location and Identification on Chromosome 3B of Bread Wheat of Genes Affecting Chiasma Number

Benoit Darrier, Isabelle Colas, Hélène Rimbert, Frédéric Choulet, Jeanne

Bazile, Aurélien Sortais, Eric Jenczewski, Pierre Sourdille

► To cite this version:

Benoit Darrier, Isabelle Colas, Hélène Rimbert, Frédéric Choulet, Jeanne Bazile, et al. Location and Identification on Chromosome 3B of Bread Wheat of Genes Affecting Chiasma Number. Plants, 2022, 11, 10.3390/plants11172281. [https://doi.org/10.3390/plants11172281](https://doi.org/10.3390/plants11172281)

HAL Id: hal-03766369

https://hal.inrae.fr/hal-03766369

Submitted on 1 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Location and Identification on Chromosome 3B of Bread Wheat of Genes Affecting Chiasma Number

Benoit Darrier <fnref n="1" />,<fnref n="2" />,†, Isabelle Colas <fnref n="3" />,†, Hélène Rimbert <fnref n="1" />, Frédéric Choulet <fnref n="1" />, Jeanne Bazile <fnref n="1" />, Aurélien Sortais <fnref n="1" />, Eric Jenczewski <fnref n="4" /> and Pierre Sourdille <fnref n="1" />,*

**Abstract:** Understanding meiotic crossover (CO) variation in crops like bread wheat (*Triticum aestivum* L.) is necessary as COs are essential to create new, original and powerful combinations of genes for traits of agronomical interest. We cytogenetically characterized a set of wheat aneuploid lines missing part or all of chromosome 3B to identify the most influential regions for chiasma formation located on this chromosome. We showed that deletion of the short arm did not change the total number of chiasma genome-wide, whereas this latter was reduced by ~35% while deleting the long arm. Contrary to what was hypothesized in a previous study, deletion of the long arm does not disturb the initiation of the synaptonemal complex (SC) in early meiotic stages. However, progression of the SC is abnormal, and we never observed its completion when the long arm is deleted. By studying six different deletion lines (missing different parts of the long arm), we revealed that at least two genes located in both the proximal (C-3BL2-0.22) and distal (3BL7-0.63-1.00) deletion bins are involved in the control of chiasmata, each deletion reducing the number of chiasmata by ~15%. We combined sequence analyses of deletion bins with RNA-Seq data derived from meiotic tissues and identified a set of genes for which at least the homoeologous copy on chromosome 3B is expressed and which are involved in DNA processing. Among these genes, eight (CAP-E1/E2, DUO1, MLH1, MPK4, MUS81, RTEL1, SYN4, ZIP4) are known to be involved in the recombination pathway.

**Keywords:** recombination; synapsis; wheat; chiasmata; 3D-SIM; meiosis; cytogenetics; deletion bin

## 1. Introduction

Mastering genes controlling meiosis and recombination is essential for wheat breeders to better exploit the wealth of diversity that exists in wild-grass resources [1–3]. Meiosis is a specialized cell division that halves the chromosome complement, or ploidy, during the production of gametes in sexually reproducing eukaryotes. To ensure that each gamete has a full complement of genetic material, homologous chromosomes must pair and then separate in a coordinated manner during meiosis. With a few exceptions (e.g., *Bombix mori* in females [4], Drosophila in males [5] and *Rhynchospora tenuis* [6]), this is mediated by the formation of genetic crossovers (COs; for a review see [7]) that reshuffles genetic material between chromosomes. Chiasmata are the cytogenetic manifestation of the COs and are essential to ensure homologue alignment on the metaphase I plate. At least one CO is mandatory per homologous pair (obligate CO [7]), and the absence of COs results in the random segregation of homologous chromosomes at anaphase I. This absence of COs leads to the formation of aneuploid gametes at the end of meiosis because of

**Citation:** Darrier, B.; Colas, I.; Rimbert, H.; Choulet, F.; Bazile, J.; Sortais, A.; Jenczewski, E.; Sourdille, P. Location and Identification on Chromosome 3B of Bread Wheat of Genes Affecting Chiasma Number. *Plants* **2022**, *11*, 2281.
https://doi.org/10.3390/plants11172281

Academic Editor: Kazuo N. Watanabe

Received: 25 May 2022

Accepted: 8 August 2022

Published: 31 August 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ UMR 1095 Genetics, Diversity and Ecophysiology of Cereals, 5, INRAE-Université Clermont-Auvergne, Chemin de Beaulieu, 63000 Clermont-Ferrand, France</footnote>

<footnote n="2">² Present address: Syngenta, Toulouse Innovation Centre 12 Chemin de l'Hobit, 31790 Saint-Sauveur, France</footnote>

<footnote n="3">³ Cell and Molecular Sciences, The James Hutton Institute, Invergowrie, Dundee, Scotland DD2 5DA, UK</footnote>

<footnote n="4">⁴ Institut Jean-Pierre Bourgin, INRAE, AgroParisTech, CNRS, Université Paris-Saclay, 78000 Versailles, France</footnote>

<footnote>† These authors contributed equally to the work.</footnote>

<footnote>* Correspondence: pierre.sourdille@inrae.fr; Tel.: +33-443-761-517</footnote>
the occurrence of univalent chromosomes resulting from homologous-pairing failure dur-
ing the prophase I stage.

Hexaploid bread wheat (*Triticum aestivum* L.; AABBD; 2n = 6X = 42) derives from two successive interspecific crosses [8–10] involving the three diploid species *T. urartu* (AA), a yet unknown species related to *Aegilops speltoides* (SS, close to BB [11–13]) and *Ae. tauschii* (DD). Bread wheat has, therefore, three sets (A, B, D) of highly similar (called ho-moeologues) chromosomes. Due to this polyploid structure, bread wheat tolerates aneu-ploidy quite well, probably because of the redundancy of the genetic information carried by the homoeologous chromosomes. Aneuploid lines missing a pair of one of the 21 chro-mosomes (nullisomic lines), a pair of chromosome arms (ditelosomic lines) or part of the chromosome (deletion lines) were produced [14–17]. They were further used for gene identification [18] as well as for molecular marker assignment [19,20]. Interestingly, the line missing chromosome 3B (nullisomic 3B) shows an aberrant meiotic behaviour [15] with numerous univalents, and it generates many aneuploid descents (monosomics or trisomics). Similarly, Naranjo [21] showed that deletions of the short arm of chromosome 3B lead to a reduction in chiasmata within this arm but do not affect genome-wide chias-mata numbers, while different deletions of the long arm of chromosome 3B result in a genome-wide reduction in chiasmata. This suggests that chromosome 3B carries essential genes required to ensure a faithful and complete meiotic recombination and that the other two chromosome-3A and -3D copies cannot compensate for the loss of 3B-homoeologous copies of these genes. However, neither of these two studies considered the complete syn-apsis of the aneuploids nor evidenced any candidate genes [15,21].

Chromosome 3B from bread wheat has been the most studied for more than 15 years.
Due to its huge physical size, this chromosome was the first to be isolated through chro-
mosome sorting and the first for which a BAC library was constructed [22]. Consequently,
chromosome 3B was also the first for which an anchored physical map was elaborated
[23] as well as a completely assembled and annotated pseudomolecule [24]. This unique
and well-documented resource was further used to study the global structure, functioning
and evolution of chromosome 3B [24–27]. This also provides an opportunity to identify
the genes located on this chromosome that affect meiotic behaviour [15,21].

Wheat chromosomes 3A, 3B and 3D are stated as playing a role in recombination. The most well-known locus is *Ph2* (for Pairing *homoeologous* 2), a gene located on the short arm of chromosome 3D (3DS; [28–30]). Evidence suggests that *Ph2* acts on homologous recombination through synaptic progression [31,32], but without affecting the recombination rate in wheat [33–35]. Actually, *Ph2* rather affects homoeologous recombination in wheat haploids (ABD), pentaploids (AABBD) and wheat-rye allo-haploid hybrids (ABDR) [28,29]. It was recently shown that *TaMsh7*, a gene involved in mismatch repair [36,37], corresponds to *Ph2* [38]. Another gene preventing homoeologous recombination is also mentioned on 3AS, but it has a weaker effect than *Ph2* [30,39]. Its position suggests that it could be a homoeologous copy of *Ph2*. Deficiency for both 3AS and 3DS results in a level of homoeologous recombination almost as high as that expected for plants missing *Ph1* [40]. *Ph1* is a cloned locus [41] involved in homoeologous recombination and recently characterized as being a wheat-specific duplicated copy of *Zip4* (*TaZip4-B2*; [42,43]). Another unknown gene preventing homoeologous recombination and which does not interact with *Ph1* was located on the long arm of chromosome 3B [44–46] as well as probable homoeo-alleles on the long arms of chromosomes 3A and 3D [29,30,39].

The genomic resources available for chromosome 3B can boost cloning of the genes involved in recombination and located on chromosome 3B, as well as on chromosomes 3A and 3D, when homoeologous copies exist [18,29,30,39]. The objectives of this paper were to study the variation in chiasma number in the set of deletion lines from chromosome 3B [17] in order to locate more precisely on this chromosome the genes that are involved in this trait. Our results showed a reduction in chiasma number in the deletion lines from the long arm. The difference between the lines suggests the presence of at least
two genes involved in this trait. In a second step, we used the reference sequence of chromosome 3B [47] to map in silico genes that are known to be involved in meiosis in model species in order to identify putative candidate genes implicated in recombination. We identified 10 genes affecting genome stability or cell cycle, located in the two deletion bins associated with a reduction in chiasma number. In addition, we present original results on the meiotic behaviour in the deletion line 3BL7-0.63-1.00 (hereafter named 3BL7) and in the ditelosomic line missing the long arm of chromosome 3B (Dt3BS), showing that deletion of the long arm does not prevent the initiation of the synaptonemal complex (SC) but hampers its completion.

## 2. Results

### 2.1. Chromosome 3B is Essential for Chiasma Number but Cannot be Compensated by Chromosomes 3A or 3D

Numbers of bivalents (rods and rings), univalents, multivalents and chiasmata were assessed on a set of 50 meiotic (metaphase I) nuclei in all genotypes (Table 1; raw data in Table S1; statistics in Table S2, S3, S4). Wild type euploid Chinese Spring variety (CS; 2n = 6x = 42 chromosomes; 14 A, 14 B and 14 D chromosomes) exhibited ~21 bivalents (~20 rings and ~1 rod; Table 1; Figure 1a) and a mean of ~41 chiasmata.

<table><caption>**Table 1.** Number of chiasmata in the different aneuploid lines. Average numbers of chiasmata, univalent, bivalent (rods and rings) and multivalent chromosomes for the wild-type Chinese Spring line (CS) and for the aneuploid lines (nullisomic-3B (N3B); nullisomic-tetrasomic-(NxxTxx); ditelosomic Dt3BL and Dt3BS; hybrids CS x ditelosomic (F1CSx3BS and F1CSx3BL); deletion lines (3BSxx or 3BLxx). % arm: percentage of the chromosome 3B or chromosome arm present (L: Long arm; S: Short arm).</caption><thead><tr><th>Line</th><th>% arm</th><th>2n</th><th>Chiasmata</th><th>Univalents</th><th>Bivalents</th><th>Rods</th><th>Rings</th><th>Multivalents</th></tr></thead><tbody><tr><td>CS</td><td>100</td><td>42</td><td>40.98 ± 1.12<br>A,B</td><td>0.04 ± 0.28<br>E</td><td>20.98 ± 0.14<br>B</td><td>1.02 ± 1.12<br>NA</td><td>19.98 ± 1.12<br>A,B</td><td>0<br>C</td></tr><tr><td>N3B</td><td>0</td><td>40</td><td>35.22 ± 1.99<br>H,IJ</td><td>0.44 ± 0.84<br>D,E</td><td>19.78 ± 0.42<br>J,K,L</td><td>4.34 ± 1.76<br>B,C</td><td>15.44 ± 1.83<br>F,G,H</td><td>0<br>C</td></tr><tr><td>N3AT3B</td><td>100</td><td>42</td><td>38.44 ± 2.35<br>D,E,F,G</td><td>0.04 ± 0.28<br>E</td><td>20.42 ± 0.91<br>B,C,D,E,F</td><td>2.34 ± 1.49<br>D</td><td>18.08 ± 1.69<br>C,D,E</td><td>0.28 ± 0.45<br>B</td></tr><tr><td>N3AT3D</td><td>100</td><td>42</td><td>38.82 ± 2.09<br>C,D,E,F,G</td><td>0.28 ± 0.70<br>D,E</td><td>20.22 ± 0.95<br>D,E,F,G,H,IJ</td><td>2.90 ± 1.93<br>C,D</td><td>17.32 ± 1.96<br>D,E,F</td><td>0.32 ± 0.47<br>A,B</td></tr><tr><td>N3BT3A</td><td>0</td><td>42</td><td>36.74 ± 1.98<br>G,H,IJ</td><td>0.88 ± 1.14<br>C,D</td><td>19.46 ± 1.16<br>IJ,K,L</td><td>4.34 ± 1.49<br>B,C</td><td>15.12 ± 1.75<br>G,H</td><td>0.56 ± 0.50<br>A</td></tr><tr><td>N3BT3D</td><td>0</td><td>42</td><td>37.32 ± 1.91<br>F,G,H,I</td><td>0.52 ± 0.89<br>D,E</td><td>19.70 ± 0.91<br>G,H,IJ,K,L</td><td>4.08 ± 1.91<br>B,C,D</td><td>15.62 ± 2.18<br>F,G,H</td><td>0.52 ± 0.50<br>A,B</td></tr><tr><td>N3DT3A</td><td>100</td><td>42</td><td>39.94 ± 1.56<br>B,C,D</td><td>0.08 ± 0.40<br>E</td><td>20.28 ± 1.07<br>B,C,D,E,F,G,<br>H</td><td>1.92 ± 1.56<br>NA</td><td>18.36 ± 1.75<br>C,D,E</td><td>0.34 ± 0.52<br>A,B</td></tr><tr><td>N3DT3B</td><td>100</td><td>42</td><td>39.82 ± 1.48<br>B,C,D</td><td>0.08 ± 0.40<br>E</td><td>20.24 ± 0.96<br>C,D,E,F,G,H,<br>I</td><td>1.98 ± 1.38<br>NA</td><td>18.26 ± 1.58<br>C,D,E</td><td>0.36 ± 0.48<br>A,B</td></tr><tr><td>Dt3BL</td><td>L</td><td>40+2t</td><td>39.44 ± 1.32<br>B,C,D</td><td>0.08 ± 0.40<br>E</td><td>20.96 ± 0.20<br>B</td><td>2.48 ± 1.23<br>D</td><td>18.48 ± 1.27<br>B,C,D,E</td><td>0<br>C</td></tr><tr><td>Dt3BS</td><td>S</td><td>40+2t</td><td>25.84 ± 2.94<br>K</td><td>6.16 ± 3.61<br>A</td><td>17.92 ± 1.81<br>L</td><td>10.00 ± 2.11<br>A</td><td>7.92 ± 1.82<br>I</td><td>0<br>C</td></tr><tr><td>F1CSx3BS</td><td>100</td><td>41+t</td><td>36.44 ± 1.89<br>G,H,IJ</td><td>2.00 ± 1.47<br>A,B</td><td>19.84 ± 0.99<br>F,G,H,IJ,K,L</td><td>4.16 ± 1.77<br>B,C,D</td><td>15.68 ± 1.89<br>F,G,H</td><td>0.44 ± 0.51<br>A,B</td></tr><tr><td>F1CSx3BL</td><td>100</td><td>41+t</td><td>39.30 ± 1.09<br>A,B</td><td>0.16 ± 0.55<br>A,B</td><td>20.92 ± 0.27<br>F,G,H,IJ,K,L</td><td>2.54 ± 1.11<br>B,C,D</td><td>18.38 ± 1.07<br>F,G,H</td><td>0<br>C</td></tr></tbody></table>
<table><thead><tr><th></th><th></th><th></th><th>C,D,E,F</th><th>D,E</th><th>B,C,D</th><th>D</th><th>C,D,E</th><th>C</th></tr></thead><tbody><tr><td>3BS3</td><td>87</td><td>42</td><td>39.64 ± 1.32<br>B,C,D</td><td>0.04 ± 0.28<br>E</td><td>20.98 ± 0.14<br>B</td><td>2.32 ± 1.33<br>D</td><td>18.66 ± 1.32<br>B,C,D</td><td>0<br>C</td></tr><tr><td>3BS8</td><td>78</td><td>42</td><td>39.92 ± 1.63<br>B,C,D</td><td>0.12 ± 0.48<br>E</td><td>20.94 ± 0.24<br>B,C</td><td>1.96 ± 1.58<br>NA</td><td>18.98 ± 1.58<br>B,C,D</td><td>0<br>C</td></tr><tr><td>3BS7</td><td>75</td><td>42</td><td>40.40 ± 1.21<br>A,B,C</td><td>0.12 ± 0.48<br>E</td><td>20.94 ± 0.24<br>B,C</td><td>1.48 ± 1.16<br>NA</td><td>19.46 ± 1.16<br>A,B,C</td><td>0<br>C</td></tr><tr><td>3BS2</td><td>57</td><td>42</td><td>40.32 ± 1.20<br>A,B,C,D</td><td>0.04 ± 0.28<br>E</td><td>20.98 ± 0.14<br>B</td><td>1.64 ± 1.16<br>NA</td><td>19.34 ± 1.71<br>A,B,C</td><td>0<br>C</td></tr><tr><td>3BS4</td><td>55</td><td>42</td><td>39.78 ± 1.53<br>B,C,D</td><td>0.04 ± 0.2<br>E8</td><td>20.98 ± 0.14<br>B</td><td>2.18 ± 1.48<br>D</td><td>18.80 ± 1.50<br>B,C,D</td><td>0<br>C</td></tr><tr><td>3BS1</td><td>33</td><td>42</td><td>40.20 ± 0.97<br>A,B,C,D</td><td>0.04 ± 0.28<br>E</td><td>20.98 ± 0.14<br>B</td><td>1.76 ± 0.96<br>NA</td><td>19.22 ± 0.95<br>A,B,C</td><td>0<br>C</td></tr><tr><td>3BS5</td><td>7</td><td>42</td><td>39.64 ± 1.50<br>B,C,D</td><td>0.08 ± 0.57<br>E</td><td>20.96 ± 0.28<br>B</td><td>2.28 ± 1.34<br>D</td><td>18.68 ± 1.39<br>B,C,D</td><td>0<br>C</td></tr><tr><td>3BL2</td><td>22</td><td>42</td><td>34.00 ± 2.48<br>I,J,K</td><td>2.20 ± 1.91<br>B</td><td>19.90 ± 0.95<br>F,G,H,I,J,K</td><td>5.80 ± 1.84<br>A,B</td><td>14.10 ± 1.92<br>H,I</td><td>0<br>C</td></tr><tr><td>3BL8</td><td>28</td><td>42</td><td>37.54 ± 1.84<br>E,F,G,H</td><td>0.36 ± 0.78<br>D,E</td><td>20.82 ± 0.39<br>B,C,D,E</td><td>4.10 ± 1.69<br>B,C,D</td><td>16.72 ± 1.73<br>E,F,G</td><td>0<br>C</td></tr><tr><td>3BL1</td><td>31</td><td>42</td><td>36.80 ± 2.16<br>G,H,I,J</td><td>1.22 ± 1.11<br>B,C</td><td>20.20 ± 0.83<br>E,F,G,H,I,J</td><td>4.52 ± 1.95<br>B,C</td><td>15.68 ± 2.03<br>F,G,H</td><td>0.46 ± 0.50<br>A,B</td></tr><tr><td>3BL9</td><td>38</td><td>41</td><td>33.10 ± 2.38<br>J,K</td><td>2.18 ± 1.51<br>A,B</td><td>19.38 ± 0.83<br>K,L</td><td>5.70 ± 1.82<br>A,B</td><td>13.68 ± 2.00<br>H,I</td><td>0.02 ± 0.14<br>C</td></tr><tr><td>3BL10</td><td>50</td><td>42</td><td>35.90 ± 2.50<br>H,I,J</td><td>1.04 ± 1.47<br>C,D</td><td>20.44 ± 0.76<br>B,C,D,E,F,G</td><td>5.04 ± 2.02<br>B</td><td>15.40 ± 2.15<br>F,G,H</td><td>0.02 ± 0.14<br>C</td></tr><tr><td>3BL7</td><td>63</td><td>41</td><td>33.46 ± 2.57<br>J,K</td><td>2.12 ± 1.22<br>A,B</td><td>19.44 ± 0.61<br>K,L</td><td>5.42 ± 2.13<br>A,B</td><td>14.02 ± 2.28<br>G,H,I</td><td>0<br>C</td></tr></tbody></table>

A–L: Multiple comparisons by pairs following the Dunn procedure/Bilateral test (α = 0.05; see Table S3 for details). Each comparison test is considered by columns only.

NA: Grouping could not be done properly as significance of differences was not transitive in this case.

As expected, almost no univalent was observed (2 in 50 PMCs), which confirmed that euploid CS exhibits a normal meiotic behaviour. On the contrary, the line missing chromosome 3B (N3B; Table 1; Figure 1b) had more univalents (mean 0.44 ± 0.84) and a significantly increased number of rod bivalents (4.34 ± 1.76; Table 1), leading to a significantly reduced number of chiasmata compared to CS (35.22 ± 1.99; Mann–Whitney *p* = 7.9 10⁻¹⁸; Table S4). This reduction was higher than expected (41–2 = ~39 chiasmata expected), showing that the effect of the loss of chromosome 3B was genome-wide and not limited to chromosome 3B itself. This confirmed the previous results from the literature [12,18], indicating that chromosome 3B of bread wheat is carrying essential genes involved in the control of chiasma number.

Chiasma number was also evaluated in Nulli-Tetrasomic (NT) lines (lines for which one pair of homologous chromosomes is replaced by one pair of their homoeologues) missing either chromosomes 3A, 3B or 3D. If the genes located on homoeologous group 3 chromosomes fully compensate for each other, then the same number of chiasmata is expected in all NT lines as in the euploid line (~41). Contrary to this prediction, all of the NT lines had a significantly reduced number of chiasmata compared to the euploid line (Table S3, S4). The aneuploid lines missing chromosome 3B (N3BT3A and N3BT3D) showed the highest difference (36.74 ± 1.98 and 37.32 ± 1.91, respectively) followed by the aneuploid lines missing chromosomes 3A (N3AT3B and N3AT3D; 38.44 ± 2.35 and 38.82 ± 2.09) and
3D (N3DT3A and N3DT3B; 39.94 ± 1.56 and 39.82 ± 1.48). This confirmed that homoeologous group 3 chromosomes bear important genes for chiasma number. This also suggests that none of the homoeologous pairs of chromosomes can fully compensate for the one missing. This finally indicates that the genes (or copies of genes) involved in chiasma number on chromosome 3B have stronger effects than those from chromosomes 3A and 3D. We thus focused our analysis on this specific chromosome.

<figure><img src="image_4.png" /></figure>

**Figure 1.** Meiotic behaviour of pollen mother cells (PMC) of wild-type, Nulli 3B and ditelosomic Dt3BS Chinese Spring lines. Wild-type cells (a) exclusively show ring bivalents (ri) at metaphase I. Nullisomic 3B cells (b) exhibit a range of ring or rod (ro) bivalents and rare univalents (u), while ditelosomic Dt3BS cells (c) contain ring-bivalents and many rod bivalents (ro) and univalents (u), suggesting loss of obligate crossovers. Scale bar = 10 µm.
## 2.2. *The Long Arm of Chromosome 3B is Essential for Chiasma Number in Wheat*

In order to narrow down the location of this (these) gene (s), similar analyses were made on double ditelosomic lines having either two copies of the short (S; Dt3BS) or the long (L; Dt3BL) arms of chromosome 3B. Since, for these two lines, one arm is missing, we expected to observe the same number of bivalents but with an additional rod and with one ring and one chiasma missing compared to CS, i.e., ~21 bivalents, ~2 rods, ~19 rings and ~40 chiasmata. For Dt3BL, we observed ~21 bivalents (2.5 rods and 18.5 rings) and 39.44 ± 1.32 chiasmata (Table 1), which is not statistically different from expectation (Wilcoxon–Mann–Whitney; α*altered* = 0.00185; respective *p*-values = 0.57; 0.04; 0.03; 0.03). This indicates that genes on this arm affecting chiasma number have only a limited effect or that their absence is fully compensated for by homoeologous genes on chromosomes 3A and 3D. On the contrary, when Dt3BS (Figure 1c) was analysed compared to normal CS, we observed only 17.92 ± 1.81 bivalents (10.00 ± 2.11 rods and 7.92 ± 1.82 rings) on average and 25.84 ± 2.94 chiasmata (p = 3.7 10<sup>-18</sup>). Especially, the number of rod bivalents was largely increased (10 times) at the expense of ring bivalents compared to normal CS. At the same time, the mean number of univalents was 6.16 ± 3.61 in Dt3BS while none were observed in euploid CS and only 0.44 ± 0.84 for N3B. These results suggest that at least one gene with a strong effect on chiasma number is present on the long arm of chromosome 3B. It is also worth noting that loss of 3BL (Dt3BS) has a stronger negative effect on chiasma frequency than loss of the entire chromosome 3B (N3B; Table 1). Since deletion of the entire short arm does not result in a decrease in chiasma number and, since deletion of the long arm only has a stronger effect than suppression of the entire chromosome, this suggests that the two arms have an opposite effect on chiasma frequencies and that the short arm of chromosome 3B may carry anti-crossover (anti-CO) genes.

We then analysed meiosis in F1 hybrids derived from the cross between CS and either Dt3BL or Dt3BS. In the first case, two doses of 3BL for one of 3BS resulted in a mean number of chiasmata (39.30 ± 1.09) that was similar (p = 0.4) to that observed in Dt3BL (39.44 ± 1.32). This is consistent with our hypothesis that 3BS does not carry any gene with a significant effect on chiasma number. On the contrary, one dose of 3BL for two doses of 3BS resulted in a mean frequency of chiasmata (36.44 ± 1.89), intermediate between those of Dt3BS (25.84 ± 2.94) and CS (40.98 ± 1.12). This suggests that the genes affecting chiasma number on the long arm show haplo-insufficiency (a situation in which the product of a single allele, although active, is synthesized in an insufficient quantity to allow the normal functioning of the cell).

From these results, we can conclude that (i) at least one gene located on the long arm of chromosome 3B is required to ensure WT crossover numbers; (ii) this gene shows haplo-insufficiency (i.e., is dosage-dependent) and (iii) at least another gene present on the short arm of chromosome 3B encodes an anti-crossover protein. We decided to locate this (these) gene(s) more precisely by evaluating bivalent/chiasma formation using a wider range of 3B-deletion lines from CS.

## 2.3. *The Long Arm of Chromosome 3B Bears at Least Two Genes Involved in Chiasma Number*

Seven deletion lines missing from 13% (3BS3) to 93% (3BS5) of the short arm of chromosome 3B were analysed. In all these lines, we observed almost no univalent (0 to 0.4/cell; Table 1; Figure S1), an average number of bivalents of ~21 and a slight variation in chiasma number between the lines, from 39.6 to 40.4 (Figure 2). These differences were not statistically significant compared to the wild type of CS (Table S3), except the number of chiasmata. However, the small drop in chiasmata and the way in which it is accentuated with larger deletions suggest an effect of the deletion in cis, i.e., deletions of the terminal part of the chromosome prevent chiasma formation in this region that is particularly prone to chiasmata in the WT. This was confirmed when the number of chiasmata of CS was reduced by one and the difference was no more significant (Table S3; 0.320 < p < 0.951). As we previously suggested from the Dt3BL result, this confirmed that if genes
involved in chiasma formation exist on 3BS, they must have a limited effect, or their loss
is fully compensated for by the presence of their homoeologues.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Average numbers of chiasmata (red crosses) for the various Chinese Spring lines: wild-type (CS), ditelosomic (Dt3BS and Dt3BL), deletion lines (3BS-xx or 3BL-xx for short and long arms respectively), hybrids between CS and Dt (F1xxx) and nulli-tetrasomics. **\*\*\***: *p* < 0.001.</figcaption></figure>

Similarly, six deletion lines missing from 37% (3BL7) to 78% (3BL2) of the long arm of chromosome 3B were analysed. They all highly significantly differed from CS in the number of chiasmata, that ranged from 33.10 ± 2.38 to 36.80 ± 2.16 (Table 1; Figure 2; *p* < 0.0001, Table S3). They also all significantly differed from Dt3BS (*p* < 0.05), except 3BL-9, that was just above the threshold (α = 0.05; *p* = 0.052). Interestingly, they did not differ significantly from N3B (0.101 < *p* < 0.543, Table S3), except 3BL-1, that was just below the threshold (α = 0.05; *p* = 0.039). This difference in chiasma number was explained by (1) the occurrence of more univalents that were observed in all the 3BL deletion lines ranging from 0.36 ± 0.78 to 2.20 ± 1.91/cell (Table 1) and (2) the increased number of rod bivalents (4.10 ± 1.69-5.80 ± 1.84/cell) at the expense of ring bivalents (13.68 ± 2.00–16.72 ± 1.73/cell). Given that the increase in univalents was paralleled by a decrease in chiasmata on bivalents, our results suggest that the deletions have an effect on genome-wide chiasma formation (in trans) in addition to an effect (in cis) at the chromosome 3B level.

The decrease in chiasma number in 3BL deletion lines and in the Dt3BS line suggests that at least two 3BL deletion bins contain genes involved in the control of this trait. Figure 2 represents the effect of these deletions on chiasma number and reveals two main drops: from ~41 to ~33 when deleting the distal bin 3BL7-0.63 and from ~33 to ~26 when deleting the proximal bin (between the centromere and 3BL2-0.22). These results suggest that the first region is located distally, in deletion bin 3BL7-0.63-1.00; this region is lost in all deletion lines, which explains why they all show less chiasmata than CS. The second region would be located proximally in deletion bin C-3BL2-0.22; this region is present in all deletion lines and only lost in Dt3BS, which explains why this line showed less chiasmata than any deletion lines. From the results of the analysis of the chromosome 3B deletion lines, we can conclude that: (1) no gene with a significant effect regarding chiasma reduction is located on the short arm; (2) the long arm of chromosome 3B carries at least two genes involved in chiasma number.

## 2.4. Comparative Analysis of Meiotic Behaviour in CS, 3bl7 and dt3bs Lines

We analysed the meiotic behaviour of lines Chinese Spring (CS; control), Dt3BS and 3BL-7. Using antibodies raised against HvZYP1 [48] and the axial element associated protein TaASY1 [49], combined with Structured Illumination Microscopy (SIM), we com-
pared synapsis in CS, Dt3BS and 3BL7 (Figure 3). Axis formation, tracked by TaASY1 (magententa), appeared normal in all three lines, as shown by the regular ASY1 signal. In CS, at the onset of the zygotene stage, the synaptonemal complex (visualized with HvZYP1 antitibody) starts to polymerize at one side of the nucleus in CS (Figure 3a), presumably at the telomere region [50,51], and brings the chromosomes in a zipper-like manner (Figure 3b,c). At the pachytene stage, the chromosomes aligned along each other via the synaptonemal complex (Figure 3d). The line 3BL-7 displays normal synapsis from zygotene to pachytene stages (Figure 3e-h), which is not significantly different from the wild type. In the Dt3BS line, synapsis appears similar to CS at the leptotene/zygotene stage (Figure 3i) and at the mid-zygotene stage (Figure 3j). However, by the late zygotene stage (Figure 3k), Dt3BS exhibits ZYP1 aggregates, which resemble a polycomplex, suggesting a defect in the maturation of the synaptonemal complex (ZYP1-polycomplex means that the SC is not normal, but that the ZYP1 protein is still in the nucleus and forms a polymer; [52]). At the pachytene stage (Figure 3l), although the chromosomes are fully aligned, ZYP1 is not fully polymerized, revealing an abnormal synapsis.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Comparison of synapsis in CS, 3BL7 and Dt3BS by 3D super resolution imaging. Chromosomes axes are labelled with ASY1 (magenta), and ZYP1 (green) is used to follow synapsis. In CS (a–d), synapsis starts at one side of the nucleus (white star; a) and progressively brings the homologues together at zygotene (b,c). Chromosomes are fully synapsed at the pachytene/diplotene stage (d). In 3BL7 (e–h), synapsis resembled that of the wild type of CS. In Dt3BS, the initiation of synapsis (i) and the progression at zygotene (j) are similar to CS. However, during zygotene, ZYP1 forms polycomplexes (white arrow; k) during synapsis, suggesting a difficulty to polymerize the central element of the synaptonemal complex, although chromosomes seem aligned with a dotty synapsis at pachytene (l). Scale bar = 5 µm.</figcaption></figure>

Thus, contrary to what was hypothesized in a previous study, Dt3BS is not fully asyn-
aptic, since we observed the formation of the synaptonemal complex (SC) in the earlier
stages. However, progression of the SC is abnormal in this line, and we never found a complete synapsis.

We also checked for possible non-homologous chromosome pairing using fluorescent DNA in situ hybridization of the rDNA probe 45S [53,54] in the same three lines, CS, 3BL-7 and Dt3BS (Figure 4). We detected the four largest 45S rDNA sites, two per homologous bivalents, on every metaphase. In CS (Figure 4a), we detected two ring bivalents carrying the two 45S signal at the telomere region. In 3BL-7 line (Figure 4b), we also detected bivalents (ring or rod) carrying two signals for 45S, suggesting that, in the 3BL7 line, homologous chromosome pairing is not affected. Similarly, the Dt3BS line (Figure 4c) seems to have normal homologous chromosome associations. Chromosome segregation occurred normally in CS (Figure 4d) and 3BL7 (Figure 4e), as shown by the same number of 45S-rDNA (red) signals equally segregating. Occasional anaphase bridges were visible, even in CS (wild type). We were not able to find the corresponding anaphase for Dt3BS, as meiotic staging is difficult in this line. However, in Dt3BS, the chromosomes appeared "fuzzier", suggesting difficulties in chromatin compaction/condensation (Figure 4f). Altogether, our results suggest that chromosome recognition is not affected in Dt3BS and 3BL-7 and that the reduction in chiasmata (crossover) is likely due to a defect in the crossover pathway during meiosis.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Homologous chromosome pairing and segregation at Metaphase I and Anaphase I in lines CS (a,d), 3BL7 (b,e) and Dt3BS (c,f). Homologous pairing was monitored using the 45S-rDNA probe. All lines exhibit correct homologous chromosome pairing as shown by the localisation of two 45s signal (red, yellow arrows) per bivalent (DAPI blue). Homologous chromosomes halve at anaphase I in both CS (d) and 3BL7 (e), with occasional chromosome bridges (white arrows) in 3BL7. We could not find a similar stage for Dt3BS. Moreover, Dt3BS chromosomes were "fuzzier" and stickier at Metaphase I (f) than in the other lines, making chromosome spread for 3D DNA in situ hybridization in this line more challenging. Scale bar 10 µm.</figcaption></figure>
## 2.5. In Silico Mapping of Genes Expressed During Wheat Meiosis or Known as Involved in Meiosis in Arabidopsis

We used the pseudomolecule of chromosome 3B [47] and our RNA-Seq data from wheat developing anthers [55] to identify genes expressed during meiosis and located within the different deletion bins (Table 2). Among the 6060 high-confidence (HC) or manually annotated (MA) genes located in our deletion bins on chromosome 3B (Table S5), 3805 and 2255 were located on the long arm and the short arm, respectively. RNA-Seq data revealed 2099 genes (mean 34.64%; 33.08% and 35.56% for short and long arm, respectively) expressed during meiosis (FPKM > 5) and physically assigned to deletion bins (Table 2). The number of expressed genes in the different deletion bins varied greatly and ranged from 18 (3BS7-0.75-0.78) to 583 (3BL7-0.63-1.00), reflecting the difference in bin size ranging from 4.5 Mb (3BS3-0.87-1.00) to 205 Mb (3BL7-0.63-1.00).

<table><caption>**Table 2.** Results of expression for the genes assigned to chromosome 3B. Number of genes expressed during meiosis assigned to the different deletion bins of chromosome 3B. Estimated size (Mb): estimated size of each deletion bin; Number of genes mapped in bins: Total number of high confidence (HC) genes mapped per bin; Number of genes expressed: Number of genes expressed (FPKM > 5) per bin; Percentage of genes expressed: associated percentage of expressed genes; Number of genes on 3B overexpressed: number of genes for which the 2^log2^ (Fold Change FPKM) is >2 in both group 3 B-A and B-D comparisons. Genes without homoeologous copies were included.</caption><thead><tr><th>Arm</th><th>Deletion Bin</th><th>Estimated Size (Mb)</th><th>Number of Genes Mapped in Bins</th><th>Number of Genes Expressed</th><th>Percentage of Genes Expressed</th><th>Number of Genes on 3B Overexpressed</th></tr></thead><tbody><tr><td rowspan="8">Short</td><td>3BS3-0.87-1.00</td><td>4.46</td><td>93</td><td>28</td><td>30.11</td><td>7</td></tr><tr><td>3BS8-0.78-0.87</td><td>31.38</td><td>519</td><td>100</td><td>19.27</td><td>23</td></tr><tr><td>3BS7-0.75-0.78</td><td>5.93</td><td>58</td><td>18</td><td>31.03</td><td>3</td></tr><tr><td>3BS2-0.57-0.75</td><td>75.56</td><td>644</td><td>189</td><td>29.35</td><td>38</td></tr><tr><td>3BS4-0.55-0.57</td><td>20.31</td><td>117</td><td>47</td><td>40.17</td><td>8</td></tr><tr><td>3BS1-0.33-0.55</td><td>105.52</td><td>608</td><td>241</td><td>39.64</td><td>43</td></tr><tr><td>3BS5-0.07-0.33</td><td>53.08</td><td>149</td><td>85</td><td>57.05</td><td>12</td></tr><tr><td>C-3BS5-0.07</td><td>48.04</td><td>67</td><td>38</td><td>56.72</td><td>11</td></tr><tr><td rowspan="6">Long</td><td>C-3BL2-0.22</td><td>90.60</td><td>416</td><td>223</td><td>53.61</td><td>40</td></tr><tr><td>3BL2-0.22-0.28</td><td>28.46</td><td>189</td><td>67</td><td>35.45</td><td>19</td></tr><tr><td>3BL8-0.28-0.32</td><td>17.85</td><td>106</td><td>48</td><td>45.28</td><td>11</td></tr><tr><td>3BL1-0.32-0.38</td><td>43.64</td><td>260</td><td>124</td><td>47.69</td><td>16</td></tr><tr><td>3BL9-0.38-0.50</td><td>44.96</td><td>359</td><td>164</td><td>45.68</td><td>25</td></tr><tr><td>3BL10-0.50-0.63</td><td>48.48</td><td>314</td><td>144</td><td>45.86</td><td>26</td></tr><tr><td></td><td>3BL7-0.63-1.00</td><td>204.93</td><td>2161</td><td>583</td><td>26.98</td><td>115</td></tr><tr><td colspan="2">Total Short Assigned</td><td>344.28</td><td>2255</td><td>746</td><td>33.08</td><td>145</td></tr><tr><td colspan="2">Total Long Assigned</td><td>478.92</td><td>3805</td><td>1353</td><td>35.56</td><td>252</td></tr><tr><td colspan="2">Total Assigned</td><td>823.20</td><td>6060</td><td>2099</td><td>34.64</td><td>397</td></tr></tbody></table>

Based on gene ontology (http://www.informatics.jax.org/vocab/gene_ontology/ (accessed on 23 August 2022); 5578 genes, 36.7%) or on their putative function (Table S6), at least 70 HC genes (Table S5; 32 on 3BS and 38 on 3BL) could play a role in the meiosis pathway. For example, two of them are good candidates that may explain the meiotic behaviour we observed when the short arm is absent. TraesCS3B02G048300 (3BS8-0.78-0.87) is annotated as a member of the Fanconi anemia group M protein, related to FANCM, that acts as an anti-crossover protein [56]. TraesCS3B02G136600 (3BS4-0.55-0.57)
is annotated as a DNA mismatch repair from the MutS family, that could also restrain crossover formation between non-identical chromosomes [38,57,58].

To confirm this, we selected a set of 104 genes from Arabidopsis, known to be involved in meiosis or playing a role in DNA repair (supplementary Table S7), to identify orthologous genes on 3B. Five mapped on the short arm of chromosome 3B (Table 3; BRCA2, MSH7, DUO1-3S, CYCA1;2/TAM1 and MUS81-3S), but none of them were known as being involved in the anti-CO pathway. Eleven of them mapped on the long arm (Table 3), among which five were located in bin C-3BL2-0.22 (RTEL1, DUO1-3L, and three copies of MPK4) and five in bin 3BL7-0.63-1.00 (CAP-E1/E2, SYN4, ZIP4, MUS81-3L and MLH1).

<table><caption>**Table 3.** Known meiotic genes assigned to chromosome 3B. Meiotic gene detected using model species protein (Arabidopsis) and exonerate alignment (score >500 identity >40% similarity >60%). Gene name determined according to Arabidopsis. Gene ID: ID of 3B-homoeologous copy gene in *Triticum aestivum* annotation. Deletion bin: Deletion bin where the gene is assigned. FPKM Values (average of four meiotic stages and two replicates [55]) are indicated for the three homoeologous copies when available; the most expressed copy is in bold.</caption><thead><tr><th>Related</th><th>Deletion Bin</th><th>Gene ID</th><th>FPKM3A</th><th>FPKM3B</th><th>FPKM3D</th></tr></thead><tbody><tr><td>BRCA2</td><td>3BS2-0.57-0.75</td><td>TraesCS3B02G11<br>5500</td><td><strong>5.26 ± 2.72</strong></td><td>1.75 ± 0.99</td><td>0.93 ± 0.51</td></tr><tr><td>MSH7</td><td>3BS4-0.55-0.57</td><td>TraesCS3B02G13<br>6600</td><td><strong>4.13 ± 2.54</strong></td><td>3.08 ± 1.76</td><td>3.37 ± 1.61</td></tr><tr><td>DUO1-3S</td><td>3BS1-0.33-0.55</td><td>TraesCS3B02G17<br>8200</td><td>0.45 ± 0.19</td><td>0.74 ± 0.39</td><td><strong>7.37 ± 4.01</strong></td></tr><tr><td>CYCA1;2/TA<br>M1</td><td>3BS1-0.33-0.55</td><td>TraesCS3B02G18<br>3400</td><td></td><td><strong>24.66 ± 7.16</strong></td><td></td></tr><tr><td>MUS81-3S</td><td>3BS5-0.07-0.33</td><td>TraesCS3B02G21<br>8300</td><td>6.29 ± 3.17</td><td>6.41 ± 3.93</td><td><strong>21.67 ± 7.85</strong></td></tr><tr><td>RTEL1</td><td>C-3BL2-0.22</td><td>TraesCS3B02G24<br>2700</td><td><strong>5.87 ± 2.53</strong></td><td>2.92 ± 1.55</td><td>2.29 ± 1.26</td></tr><tr><td>DUO1-3L</td><td>C-3BL2-0.22</td><td>TraesCS3B02G25<br>4800</td><td><strong>4.77 ± 2.45</strong></td><td>4.07 ± 2.52</td><td>3.30 ± 1.77</td></tr><tr><td>MPK4-1</td><td>C-3BL2-0.22</td><td>TraesCS3B02G25<br>6700</td><td><strong>21.08 ± 7.62</strong></td><td>9.83 ± 3.93</td><td>7.61 ± 4.64</td></tr><tr><td>MPK4-2</td><td>C-3BL2-0.22</td><td>TraesCS3B02G26<br>0900</td><td>4.07 ± 2.39</td><td>4.47 ± 2.64</td><td><strong>7.44 ± 3.32</strong></td></tr><tr><td>MPK4-3</td><td>C-3BL2-0.22</td><td>TraesCS3B02G27<br>0200</td><td>13.01 ± 7.36</td><td><strong>17.93 ± 8.58</strong></td><td>14.06 ± 7.12</td></tr><tr><td>ASK1/ASK2</td><td>3BL1-0.32-0.38</td><td>TraesCS3B02G30<br>8600</td><td></td><td></td><td></td></tr><tr><td>CAP-E1/E2</td><td>3BL7-0.63-1.00</td><td>TraesCS3B02G42<br>3800</td><td>11.63 ± 4.45</td><td><strong>18.08 ± 9.11</strong></td><td>11.23 ± 6.71</td></tr><tr><td>SYN4</td><td>3BL7-0.63-1.00</td><td>TraesCS3B02G42<br>9700</td><td>20.58 ± 9.22</td><td>16.23 ± 9.71</td><td><strong>31.22 ± 8.59</strong></td></tr><tr><td>ZIP4</td><td>3BL7-0.63-1.00</td><td>TraesCS3B02G43<br>4600</td><td>5.79 ± 2.94</td><td>5.91 ± 3.28</td><td><strong>6.07 ± 2.13</strong></td></tr><tr><td>MUS81-3L</td><td>3BL7-0.63-1.00</td><td>TraesCS3B02G53<br>5000</td><td>9.96 ± 6.53</td><td>19.29 ± 3.30</td><td><strong>36.59 ± 11.18</strong></td></tr><tr><td>MLH1*</td><td>3BL7-0.63-1.00</td><td>TraesCS3B02G56<br>4100</td><td>3.90 ± 1.89</td><td><strong>13.29 ± 4.19</strong></td><td>12.06 ± 4.84</td></tr></tbody></table>

<footnote>* There are two copies of *Mlh1* gene on chromosome 3D, and the expression value is a mean of the two.</footnote>
Interestingly, DUO1 and MUS81 have duplicates on both arms of this chromosome and MPK4 has triplicates in the same deletion bin. We compared the expression of these candidate genes with that of their homoeologous copies using a set of already available RNA-Seq data from the literature as transcript reference [55]. Copies from chromosome 3B were the most highly expressed for CAP-E1/E2, MLH1 and MPK4-3 (Figure 5). For ZIP4, the expression level of the three copies was balanced. These genes can be considered as good candidates responsible for the control of chiasma formation.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Level of expression (FPKM values; average of the four stages and two replicates) in meiotic tissues [55] of the different homoeologous copies (A: orange; B: yellow; D: green) of genes known as being involved in recombination in model species and located on chromosome 3B.</figcaption></figure>

We evaluated the potential functionality of the three copies for all these genes by comparing the protein sequence of eleven wheat varieties with those from the diploid species (*T. urartu* and *Ae. tauschii*) that should be functional to ensure a good fertility of these two species (Table S8). Based on these analyses, we estimate that all A and D copies of the 16 genes should be functional except BRCA2-A, MPK4-D1, MPK4-A3 and ZIP4-A. We then compared the A, B and D copies, and results suggest that MSH7-3B, RTEL-3B and MUS81-3BL may not be functional. However, this must be confirmed using a genetic approach with appropriate mutants.

### 3. Discussion

#### 3.1. *Variation in Chiasma Number in the Aneuploid Stock From Chromosome 3B*

Chromosome 3B has been found to carry genes playing a role in pairing in many studies [18,21,44–46,59]. In this study, we used the aneuploid stocks of wheat, namely, the nullisomic, nulli-tetrasomic, ditelosomic and deletion lines related to chromosome 3B, to locate the genes involved in recombination and chiasma number on this chromosome. Due to the polyploid nature of the wheat genome, homoeologous copies can compensate for the lack of one copy. Thus, only genes that are present in a single copy, those that are sensitive to a dosage effect or those with the 3B copy having the major effect, could be detected in our
study. Moreover, our observations rely on chiasma counting of Pollen Mother Cells (PMC) that underwent a complete meiosis until metaphase I. However, the genes concerned by the deletions could have affected the earlier phases of meiocyte development and, perhaps, a large number of cells could have died before reaching metaphase. We observed a fair amount of synapsis even in Dt3BS, although we could never find a complete synapsis in this line. The timing of meiosis may be disturbed, and we saw cells that could have been late at synapsis. A time course analysis would allow us to answer this question. We did not observe loss in fertility in any of the aneuploid lines we used. This would have been important information, since mutated genes could reach an earlier control for pairing and recombination with regard to the observation of meiotic behaviour.

We observed a significant but limited effect of the partial or total loss of the short arm of the chromosome 3B on chiasma number. This suggests that this chromosome arm carries genes involved in this trait, but with only a tiny effect. Interestingly, deletion of the long arm only (Dt3BS) has a stronger effect than suppression of the entire chromosome (N3B) with~26 and ~35 chiasmata, respectively. This suggests that the two arms have an opposite effect on chiasma frequencies and that anti-crossover genes locate on the short arm of chromosome 3B. Anti-crossover protein-encoding genes have been found in Arabidopsis, such as FANCM, FIGL1, RECQ4 and TOP3α [56,60–65]. None of the true wheat orthologues of these genes are located on the short arm of homoeologous group 3 chromosomes. However, we found high-confidence genes annotated as members of the FANCM and TOP3α proteins. There are also several other DNA helicase-encoding genes that have been predicted. This suggests that if anti-crossover protein-encoding genes are present on the short arm of chromosome 3B, they are different from those that have already been isolated. This may also reflect the fact that some genes from the short arm of chromosome 3B act as a complex with those from the long arm and/or regulate their activity. Once the short arm only is present, the genes it carries may more strongly regulate the homoeologous copies from chromosome arms 3AL and 3DL.

On the contrary, partial or total deletion of the long arm of chromosome 3B revealed more contrasted results. Deletion of the terminal bin (205 Mb) strongly disturbs the number of chiasmata. In addition, complete removal of the long arm increases this perturbation of pairing and decreases recombination. This suggests the presence of at least two genes, one located in deletion bin 3BL7-0.63-1.00 and the second in C-3BL2-0.22. Naranjo [21] also observed a reduced number of chiasmata in deletion lines 3BL6-0.54 and 3BL3-0.41 (23.3 and 29.6, respectively), while the deletion line 3BL11-0.81 exhibited normal behaviour (40.2 chiasmata). This allows for refining the position of the gene located distally in the interval determined by the breakpoints from the deletion lines 3BL7-0.63 and 3BL11-0.81, which halves the size of the segment from 205 Mb to 100 Mb.

The presence of a pairing promoter gene located on the long arm of chromosome 3B is necessary for normal synapsis and chiasma formation [18,59]. This gene is probably the one we identified in bin 3BL7-0.63-1.00. Evaluation of wheat radiation hybrids indicated the presence of a gene (*TaDES2*) which, when deleted, leads to a reduced fertility that could be due to a problem in pairing [66], which would confirm our results. In addition, Dt3BS was previously described as asynaptic, which implies a defect in pairing [67]. Implication of the long arm of chromosome 3B in pairing is in accordance with previous results which mentioned the effect of this region, either directly in wheat [18,35,46] or with crosses with related species such as *Aegilops* [68,69]. Moreover, the homoeologous chromosome arms, 3AL and 3DL, carry promoters of pairing that could be homoeo-alleles of the gene on 3BL [28,29,39].

### 3.2. Identification of Candidate Genes Located in the Two Deletion Bins Concerned

Combining in silico assignment of genes known to be involved in meiosis and RNA-Seq data led to the identification of several candidate genes in deletion bins C-3BL2-0.22 (RTEL1, DUO1 and MPK4) and 3BL7-0.63-1.00 (MLH1, ZIP4, SYN4, CAP-E1/E2 and
MUS81), among which the 3B copies of MPK4, MLH1 and CAP-E1/E2 were probably func-
tional and more expressed than their homoeologues.

MPK4 (Mitogen-activated protein kinase) was initially described as playing a role in regulating cortical microtubule bundling [70], but, recently, it was evidenced that *mpk4* mutants develop normal microspore mother cells that further fail to complete meiotic cytokinesis, leading to aborted pollen grains [71]. Here, we find a related MPK4 protein expressed during meiosis, suggesting an implication in regulating meiosis division. The meiotic phenotype we observed does not suggest that cytokinesis are affected in our deletion lines. Thus, MPK4 would not be the best candidate, even if we could not exclude that it can be involved in the regulation of another gene that, if less or differentially expressed, will lead to the appropriate phenotype.

MLH1 belongs to the DNA mismatch repair system (MMR) and is the homologue to the *MutL* gene originally described in *Escherichia coli*. Three homologues of *MutL* exist in Arabidopsis (*AtMLH1*, *AtMLH3* and *AtPMS1*; [72,73]). It was shown that *AtMLH1* and *AtMLH3* act as dimers during the pachytene stage of meiosis [74]. The mutation of *AtMLH1* evidenced that it is required for homologous recombination and that it also acts to limit recombination between divergent sequences [75]. In bread wheat, it was shown that the *Ph1* locus [40] prevents MLH1 sites from becoming crossovers on paired homoeologues during meiosis [76]. In our study, MLH1 could be a good candidate to explain the reduction in chiasma number since it contributes to resolving the crossovers at the end of meiosis. If MLH1 is inactivated, we should thus observe less chiasmata compared to the wild type.

CAP-E1 is a functional orthologue of SMC2 (Structural Maintenance of Chromosomes 2), a gene that belongs to the condensin complex involved in mitotic chromosome condensation and dosage compensation in yeast [77]. Two copies (CAP-E1 and CAP-E2) were identified in Arabidopsis, and mutant analysis revealed that the genes are expressed during meiosis and that double heterozygous mutants exhibited a reduction in condensation at metaphase and anaphase I and a mild cut-type phenotype at anaphase I [78]. Since chiasma resolution requires that sister chromatids decatenate from each other along the chromosome length, when CAP-E1 and/or CAP-E2 are mutated, this decatenation fails, and the failure results in a lower number of chiasmata and a higher number of univalents. Moreover, chiasma frequency is significantly reduced from 9.05 chiasma per cell in the wild type to 8.15–8.44 chiasma per cell in AtSMC4 RNAi mutants [79], suggesting a role of the condensin complex during meiotic prophase I. CAP-E1 could thus be a good candidate as well to explain our reduced number of chiasmata in our deletion lines.

Despite the fact that the homoeologous copy from chromosome 3B of these three
genes is the most expressed among the three copies, and that they thus constitute prefer-
ential candidates, we cannot exclude that, when they are deleted, the expression of their
3A and/or 3D homoeologous copies could be upregulated to compensate for their absence.
This would lead to an absence of any meiotic phenotype. Looking at the expression of
these genes in different deletion stocks such as N3B, Dt3BL or Dt3BS could enable us to
dissect whether their expression level may indeed play a role in the observed phenotype.

We cannot exclude that the other genes that we located in the two deletion bins also play
a role in the phenotype observed. For example, ZIP4 is involved in the formation of the syn-
aptonemal complex in rice [80] but not in Arabidopsis [81] and SYN4 in centromere cohesion
[82]. MUS81 is involved in class II crossovers and resolves aberrant intermediates induced
during replication as well as for efficient synthesis-dependent strand annealing (SDSA) and,
to a smaller extent, for single-strand annealing (SSA) [83]. It could thus affect homologous
recombination frequency, and its removal would explain the reduced number of chiasmata
that we observed in the lines missing part or all of the long arm of chromosome 3B.

Moreover, we had more than 800 genes that were expressed during meiosis in the 2 deletion bins C-3BL2-0.22 and 3BL7-0.63-1.00. We cannot exclude that at least one of these genes that is currently unknown as playing a role in meiosis could be responsible for the observed phenotype. Then, mutants for these genes will have to be analysed individually, and the resources developed for mutant identification (http://www.wheat-tilling.com/,
addressed on 20 April 2022) will be of main importance to accelerate this process. Finally, all the genes that are located in the two bins may also act in complexes or may regulate each other, which will largely complicate the identification of those having the largest impact on the phenotype we observed.

### 3.3. Conclusion

Our study showed that at least two genes impact chiasma number and locate distally (3BL7-0.63-1.00) and proximally (C-3BL2-0.22) on the long arm of chromosome 3B of bread wheat. The formation of bivalents at meiosis results from three major processes: chromosome pairing itself, which is the interaction between homologous chromosomes resulting in their alignment; synapsis, which is the formation of the synaptonemal complex between homologues; and crossovers, which occur at the end of prophase I. Meiotic behaviour of our mutants needs to be better deciphered to understand which part of the process is impacted and to better focus on the appropriate candidates. The recent genomic resources (whole-genome sequence [46], expression data [84]) that are now available for wheat provided us with the opportunity to find putative candidate genes. They now need to be studied in more detail by searching for mutants in TILLING populations (http://www.wheat-tilling.com/, 20 April 2022) or by developing transgenic plants (knock-out or over-expression) and evaluating the chiasma number and recombination behaviour of the mutants.

## 4. Materials and Methods

### 4.1. Plant Material

To estimate the location of genes located on chromosome 3B involved in chiasma number at meiosis, we used the reference line Chinese Spring as a control and the aneuploid stock of chromosome 3B described earlier [16–18,85]. The aneuploid lines consisted of the nullisomic 3B (N3B), the nullisomic–tetrasomic (NT) lines, two ditelocentric lines (Dt3BL and Dt3BS) and thirteen chromosome-3B deletion lines (Table 1), seven for the short arm (3BS3-0.87, 3BS8-0.78, 3BS7-0.75, 3BS2-0.57, 3BS4-0.55, 3BS1-0.33, 3BS5-0.07) and six for the long arm (3BL2-0.22, 3BL8-0.28, 3BL1-0.31, 3BL9-0.38, 3BL10-0.50, 3BL7-0.63). Each deletion line is defined by the fraction length of the arm that remains present on the deleted chromosome (for example, line 3BS4-0.55 conserves proximal 55% and lacks the distal 45% of the short arm of chromosome 3B). These deletion lines together with Dt3BL and Dt3BS define 15 deletion bins: 3BS3-0.87-1.00, 3BS8-0.78-0.87, 3BS7-0.75-0.78, 3BS2-0.57-0.75, 3BS4-0.55-0.57, 3BS1-0.33-0.55, 3BS5-0.07-0.33, C-3BS5-0.07, C-3BL2-0.22, 3BL2-0.22-0.28, 3BL8-0.28-0.31, 3BL1-0.31-0.38, 3BL9-0.38-0.50, 3BL10-0.50-0.63 and 3BL7-0.63-1.00. Each deletion bin is defined by the fraction length of two successive deletion lines (for example, the deletion bin 3BS4-0.55-0.57 is the chromosome fragment covering the region between 55% and 57% of the short arm of chromosome 3B and delimited by deletion lines 3BS4-0.55 and 3BS2-0.57). Plants grew in growth chambers until the first two leaves emerged. They were then potted and grown for two weeks in cool conditions (17 °C day and 12 °C night, 12 h of day) and then for two months at 24 °C during the day and 18 °C during the night with 16h of day until the spikes arrive at the appropriate stage to observe chiasmata (metaphase I). Two plants for each studied line were sown to check for the homogeneity of the deletion and one of the two plants was randomly selected for the cytogenetic analyses.

### 4.2. Chiasma Observation and Statistical Analysis

Young spikes were collected and the metaphase I stage was determined by evaluating one anther per flower through a rapid squash under carmine–acetic coloration. The two other anthers were fixed in a solution of alcohol/acetic acid (3/1) for 48h and were then dissected in aceto-carmine added with FeCl3 using an Axioskop (Zeiss) dissecting microscope. Only pollen mother cells (PMC) were conserved for further analyses. The
slides were heated until separation of the chromosomes and were washed with acetic wa-
ter (45%). They were then slightly pressed to spread out the chromosomes and they could
be stored in the fridge for rapid observation or in alcohol at -20°C after removal of the
slide. Cytoplasmic proteins were removed with pepsin followed with a paraformalde-
hyde treatment and the slides were washed twice (2X SSC) and fixed in successive alcohol
solutions (70%, 95% and 100%). Squashed anther preparations were then viewed using
the Axioskop. To estimate chiasma frequencies, 50 PMCs in metaphase I were scored for
each plant and numbers of univalents, rod and ring bivalents and multivalents were rec-
orded (Table S1). To estimate the mean chiasma frequencies, we assumed that univalents,
rod bivalents, ring bivalents, trivalents and quadrivalents contribute respectively with
zero, one, two, three and four chiasmata.

Normality of the data was assessed using a Shapiro-Wilk test (Table S2). Significant
differences between mutant and corresponding wild-type control chiasma frequencies
were computed using either non-parametric Kruskal-Wallis (Table S3) or Mann-Whitney
(Table S4) tests adjusted (Bonferroni correction) for multiple comparisons.

4.3. *Rna-Seq Data Production and in silico Gene Mapping*

To see if the genes that we located in the relevant deletion bins were indeed expressed, and if we could find any difference between gene expression of each of the homoeologous copies, we used already available RNA-Seq data from the literature as a transcript reference [55]. These latter were derived from four different stages of developing anthers (Latent/Leptotene; Zygotene/Pachytene; Diplotene/Diakinesis; Metaphase I; data publicly available: http://wheat-urgi.versailles.inra.fr/Seq-Repository/Expression; accessed on 20 December 2021). Briefly, total RNAs were extracted from 30 to 50 anthers (depending on the stage) of the hexaploid wheat cv Chinese Spring and from 10 mg of anthers using the Macherey-Nagel Nucleospin RNA-XS kit, according to the manufacturer's instructions. RNA-Seq non-oriented libraries were constructed in duplicates using the TruSeq kit (Illumina). The eight libraries were sequenced (GATC, Konstanz, Germany) on two lanes (four samples per lane) of HiSeq2000 (Illumina) with paired-end sequence (500 bp) in 2 × 100 bp, which generated 40 to 50 million pairs of reads per sample. For alignment of the reads, we used the sequences available through the International Wheat Genome Sequencing Consortium (http://www.wheatgenome.org/, 20 April 2022). All the reads from the RNA-Seq libraries were mapped on the scaffolds representing the gene models produced from the assembly of the reads [47]. We used TopHat2 v2.0.8 (http://tophat.cbcb.umd.edu/; [86]) with the default parameters, except that we tolerated no mismatches or splice-mismatches. PCR duplicates that are generated during library construction were removed from the analysis using SAMTOOLS (rmdup option). The transcripts reconstruction and expression determination (FPKM; [87]) were analysed with Cufflinks v2.0.2 (http://cufflinks.cbcb.umd.edu/). The genes were previously annotated, and mapping of the reads was performed on exons only, which allowed discrimination of expression of the homoeologous copies.

For candidate-gene identification, we selected the genes expressed during meiosis
and a set of 104 genes (see supplementary Table S7) known to be involved in meiosis in
Arabidopsis, rice, maize and *Aegilops*. We mapped all these genes onto the pseudomole-
cule of chromosome 3B [47,88] using their protein sequence determined through an exon-
erate (V2.2.0) analysis [89]. We confirmed their position using RNA-Seq data. Only those
showing the best hits for alignment (scores > 500, identity > 40%, similarity > 60%) were
considered as relevant. Only predictions with an FPKM value > 5 in at least one condition
were considered for expression analysis. In addition, a comparative analysis of expression
between homoeologous copies was conducted. Expression on homoeologous chromos-
omes 3A and 3D of genes identified on chromosome 3B was searched for using the
IWGSC survey protein database and blastp. Homoeologous gene copies were selected ac-
cording to the closest similarity with the copy from chromosome 3B.
For comparison of the functionality of the different copies (Table S8), we used available sequences from the diploid species *T. urartu* (AA), *Ae. tauschii* (DD), the tetraploid species *T. durum* (cv. Svevo; AABB), *T. diccocoides* (cv. Zavitan; AABB) and 11 hexaploid wheat varieties (AABBDD; CS, Arina-Lr-Forno, Jagger, Julius, Lancer, Landmark, Mace, Norin-61, Spelt, Stanley, SY-Mattis [90]). Various copies were retrieved through a blasp-p analysis [88] using the sequence from the corresponding proteins in Arabidopsis (Table S7). The different copies were then aligned to identify differences. We hypothesize that the copies are functional in the diploid progenitors. First, the A and D copies of *T. urartu* and *Ae. tauschii* were compared to identify neutral mutations that do not affect functionality. These latter were removed from further analyses. Then, all A and D copies from the tetraploid and hexaploid varieties were compared separately to the copies of *T. urartu* and *Ae. tauschii*, respectively, to identify potential recurrent deleterious mutations in the A or D copies that may suggest loss of functionality in either of the two. Finally, the B copies of the tetraploid and hexaploid varieties were compared to their homoeologous A and D copies to check for their putative functionality.

### 4.4. Immunostaining

Samples were prepared according to Colas et al. [91] with slight modifications for wheat [48]. Four rounds (sown at two weeks intervals in the same growing condition) of three plants per pot for each genotype were used. Three to four tillers of each plant (excluding the main tiller) were examined, and each floret was staged under acetocarmine before proceeding to the cytological experiment. Since it is not possible to find all the different stages in only one plant, the experiments were conducted over several weeks. In wheat, each floret corresponds to a different meiotic stage along the spike, but all three anthers per floret were synchronized [92]. For each line, we selected anthers of the same size, and we observed them to evaluate the meiotic stage before proceeding to the next steps. Each experiment was conducted in both wild type and deletion lines as described in [48]. For each line, between 10 and 30 cells were observed.

Staged anthers were collected in 1X-PBS containing 0.5% Triton® X100 in an embryo dish and subsequently fixed in freshly made 4% formaldehyde solution for 30 min at room temperature. Fixed anthers were washed twice in 1X-PBS/0.5% Triton® X100 and tapped to release the meiocytes in the embryo dish. Meiocytes (30 µL) were transferred onto a Polysine® slide (Poly-L-Lysine coated slides) and left to air-dry gently. Slides were first incubated for 1h at room temperature in 1X-PBS, 1% Triton® X100 for 30-45 min to permeabilize the cells in a blocking solution consisting of 3% BSA in 1X-PBS, 0.1% TritonX100. Slides were incubated in the primary antibody solution in a wet chamber for 1h at room temperature followed by 24 to 36 hours in the fridge (4-6 °C). The primary antibody solution consisted of one anti-TaASY1 (rabbit, 1:5000) [49,93] and anti-HvZYP1 (rat, 1:500) [48]. After incubation, slides were let for an additional 1h at room temperature before washing them for 15min in 1X-PBS and incubated for 90-120 minutes at room temperature in the secondary antibody solution consisting of a mixture of anti-rabbit Alexa Fluor (488 or 568) and/or anti-rat Alexa Fluor (568 or 488) diluted in 1X-PBS (1:300). After a 15min wash in 1X-PBS, the slides were counterstained with 1µg/mL Hoescht 33342 for 10-15 min, mounted in Vectashield (H-1000) and sealed with nail varnish.

### 4.5. Fluorescence DNA in situ Hybridization

The protocol is adapted from Colas et al. [91]. Briefly, anthers were fixed in Ethanol:Acetic Acid (3:1) for at least 24h and then transferred into 70% ethanol for storage. Anthers were spread on polysine slides, according to [51], and slides were air-dried. The 45S probe was prepared from the 45S rDNA clone from wheat (pTa71 [94]) and labelled with biotin, according to [95] and [96]. The hybridization mix (50% deionised formamide, 20% dextran sulfate, 1X PIPES/EDTA buffer (100:10), 0.3 M NaCl, 500 ng of salmon sperm blocking DNA, and 50 ng of each probe) was denatured for 8 min at 100°C and immediately placed on ice to cool down for 5 min. Slides were successively dehydrated in 30, 50,
70 and 100% ethanol. Ten µL of the probe mixture was applied onto the tissue, covered with a plastic coverslip and denatured on a hotplate for 5 min at 75 °C. Slides were moved into a wet chamber and incubated at 37 °C overnight. The next day, slides were washed at 42 °C in 20% formamide, 0.1X SSC for 10 min, then in 2X SSC for 10 min. This was followed by a wash at room temperature in 2X SSC for 10 min, and in 4X SCC, 0.2% Tween 20 for 10 min. A blocking solution, consisting of 5% BSA in 4X SSC, 0.2% Tween 20 was applied for 5 min in a humidity chamber at room temperature. Biotin-labelled probes were detected with extravidin-cy3. Slides were counter stained in 1µg/mL Hoescht 33342.

## 4.6. *Microscopy and Imaging*

The 3D-SIM acquisition of individual meiocytes was sequentially obtained using an OMX Blaze structured illumination microscope (GE Applied Precision) for three channels (405, 488 and 568 nm); SIM image reconstruction and processing were performed with the commercial software Softworx (GE Applied Precision) [97]. The 3D Confocal stack images were acquired with LSM-Zeiss 710 for the three channels sequentially, with an averaging of at least four, and de-convolved with the Tikhonov-Miller algorithm from the DeconvolutionLab package [98] for Fiji/ImageJ [99]. Projections of 3D pictures and light brightness/contrast adjustment were performed with the public domain program Fiji and/or Imaris 7.1 (Bitplane).

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants11172281/s1; Figure S1: Scatterplots of the number of chiasmata, univalents, bivalents, rods, rings and multivalents for all the deletion and aneuploid lines; Table S1: Raw data of the number of chiasmata, univalents, bivalents, rods, rings and multivalents for all the deletion and aneuploid lines; Table S2: Normality tests (Shapiro-Wilk) of the number of chiasmata, univalents, bivalents, rods, rings and multivalents for all the deletion and aneuploid lines; Table S3: Kruskal-Wallis tests of the number of chiasmata, univalents, bivalents, rods, rings and multivalents for all the deletion and aneuploid lines; Table S4: Mann-Whitney tests of the number of chiasmata, univalents, bivalents, rods, rings and multivalents for all the deletion and aneuploid lines. Table S5: List of genes present on chromosome 3B; analysis of those that are over-expressed; expression of those known to be involved in recombination; Table S6: Gene ontology terms; Table S7: ID and protein sequences (Arabidopsis) of the genes known as involved in recombination pathway; Table S8: Comparative analyses of the different homoeologous copies (proteins) of the putative candidate genes to assess for their functionality;

**Author Contributions:** P.S., E.J. and F.C. conceptualized the experiments; B.D. and I.C. performed the cytogenetic experiments and analysed the data; H.R. and F.C. conducted the bioinformatics analyses; J.B. and A.S. calculated the statistics and analysed the functionality of the proteins; P.S. and E.J. wrote the original draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received funding for the OMX microscope through the MRC Next Generation Optical Microscopy Award (ref: MR/K015869/1) at the CAST facility, University of Dundee. IC was funded by the European Community's Seventh Framework Programme FP7/2007–2013 (no. 222883), ERC Shuffle (Project ID: 669182) and BBSRC BB/T008636/1.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All raw data are provided in Table S1.

**Acknowledgments:** We are grateful to B. Gill and J. Raupp (Kansas State University, USA) for providing the deletion lines. We would like to thank our colleagues from greenhouses (VégéPôle) and GENTYANE for helping in preparing the material and for the support in the expression analysis. BD was funded by Meta-Programme SELGEN (project BO-DeLiRe) and Region Auvergne. Georges Gay is greatly acknowledged for collecting metaphases and counting chiasmata. Heidi Serra is also acknowledged for her help in statistics. We are also grateful to Michel Bernard and Catherine Feuillet for having supported the project.

**Conflicts of Interest:** The authors declare no conflict of interest.
## References

1. Balfourier, F.; Bouchet, S.; Robert, S.; De Oliveira, R.; Rimbert, H.; Kitt, J.; Choulet, F.; Paux, E.; International Wheat Genome Sequencing Consortium; BreedWheat Consortium. Worldwide phylogeography and history of wheat genetic diversity. *Sci. Adv.* **2019**, *5*, eaav0536. https://doi.org/10.1126/sciadv.aav0536.

2. He, F.; Pasam, R.; Shi, F.; Kant, S.; Keeble-Gagnere, G.; Kay, P.; Forrest, K.; Fritz, A.; Hucl, P.; Wiebe, K.; et al. Exome sequencing highlights the role of wild-relative introgression in shaping the adaptive landscape of the wheat genome. *Nat. Genet.* **2019**, *51*, 896–904. https://doi.org/10.1038/s41588-019-0382-2.

3. Pont, C.; Leroy, T.; Seidel, M.; Tondelli, A.; Duchemin, W.; Armisen, D.; Lang, D.; Bustos-Korts, D.; Goué, N.; Balfourier, F.; et al. Tracing the ancestry of modern bread wheats. *Nat. Genet.* **2019**, *51*, 905–911. https://doi.org/10.1038/s41588-019-0393-z.

4. Rasmussen, S.W. Meiosis in *Bombyx mori* Females. *Phil. Trans. R. Soc. Lond.* **1977**, *277*, 343–350.

5. Tsai, J.H.; McKee, B.D. Homologous pairing and the role of pairing centers in meiosis. *J. Cell Sci.* **2011**, *124*, 1955–1963. https://doi.org/10.1242/jcs.006387.

6. Cabral, G.; Marques, A.; Schubert, V.; Pedrosa-Harand, A.; Schlögelhofer, P. Chiasmatic and achiasmatic inverted meiosis of plants with holocentric chromosomes. *Nat. Comm.* **2014**, *5*, 5070. https://doi.org/10.1038/ncomms6070.

7. Mercier, R.; Mézard, C.; Jenczewski, E.; Macaisne, N.; Grelon, M. The Molecular Biology of Meiosis in Plants. *Annu. Rev. Plant Biol.* **2015**, *66*, 297–327. https://doi.org/10.1146/annurev-arplant-050213-035923.

8. Blake, N.K.; Lehfeldt, B.R.; Lavin, M.; Talbert, L.E. Phylogenetic reconstruction based on low copy DNA sequence data in an allopolyploid: The B genome of wheat. *Genome* **1999**, *42*, 351–360. https://doi.org/10.1139/G98-136.

9. Feldman, M.; Lupton, F.; Miller, T. Wheats. In *Evolution of Crop Plants*; Smartt, J., Simmonds, N.W., Eds.; Longman Group: London, UK, **1995**; pp. 184–192.

10. Huang, S.; Sirikhachornkit, A.; Faris, J.D.; Su, X.; Gill, B.S.; Haselkorn, R.; Gornicki, P. Phylogenetic analysis of the acetyl-CoA carboxylase and 3-phosphoglycerate kinase loci in wheat and other grasses. *Plant Mol. Biol.* **2002**, *48*, 805–820. https://doi.org/10.1023/A:1014868320552.

11. Glémin, S.; Scornavacca, C.; Dainat, J.; Burgarella, C.; Viader, V.; Ardisson, M.; Sarah, G.; Santoni, S.; David, J.; Ranwez, V. Pervasive hybridizations in the history of wheat relatives. *Sci. Adv.* **2019**, *5*, eaav9188. https://doi.org/10.1126/sciadv.aav9188.

12. Li, L.-F.; Zhang, Z.-B.; Wang, Z.-H.; Li, N.; Sha, Y.; Wang, X.-F.; Ding, N.; Li, Y.; Zhao, J.; Wu, Y.; et al. Genome sequences of five Sitopsis species of Aegilops and the origin of polyploid wheat B subgenome. *Mol. Plant* **2022**, *15*, 488–503. https://doi.org/10.1016/j.molp.2021.12.019.

13. Avni, R.; Lux, T.; Minz-Dub, A.; Millet, E.; Sela, H.; Distelfeld, A.; Deek, J.; Yu, G.; Steuernagel, B.; Pozniak, C.; et al. Genome sequences of three Aegilops species of the section Sitopsis reveal phylogenetic relationships and provide resources for wheat improvement. *Plant J.* **2022**, *110*, 179–192. https://doi.org/10.1111/tpj.15664.

14. Sears, E.R. Cytogenetic Studies with Polyploid Species of Wheat. I. Chromosomal Aberrations in the Progeny of a Haploid of *Triticum vulgare*. *Genetics* **1939**, *24*, 509–523. https://doi.org/10.1093/genetics/24.4.509.

15. Sears, E.R. Cytogenetic studies with polyploid species of wheat. II. Additional chromosomal aberrations in *Triticum vulgare*. *Genetics* **1944**, *29*, 232–246. https://doi.org/10.1093/genetics/29.3.232.

16. Sears, E.R.; Sears, L. The telocentric chromosomes of common wheat. In Proceedings of the 5th Wheat Genetics Symposium, New Delhi, India, 23–28 February 1978.

17. Endo, T.R.; Gill, B.S. The Deletion Stocks of Common Wheat. *J. Hered.* **1996**, *87*, 295–307. https://doi.org/10.1093/oxfordjour-nals.jhered.a023003.

18. Sears, E.R. *The Aneuploids of Common Wheat*; University of Missouri College of Agriculture, Agricultural Experiment Station: Columbia, MO, USA, 1954; pp. 1–58.

19. Sharp, P.J.; Chao, S.; Desai, S.; Gale, M.D. The isolation, characterization and application in the *Triticaceae* of a set of wheat RFLP probes identifying each homoeologous chromosome arm. *Theor. Appl. Genet.* **1989**, *78*, 342–348. https://doi.org/10.1007/BF00265294.

20. Anderson, J.A.; Churchill, G.A.; Autrique, J.E.; Tanksley, S.D.; Sorrells, M.E. Optimizing parental selection for genetic linkage maps. *Genome* **1993**, *36*, 181–186. https://doi.org/10.1139/g93-024.

21. Najanjo, T. Forcing the shift of the crossover site to proximal regions in wheat chromosomes. *Theor. Appl. Genet.* **2015**, *128*, 1855–1863. https://doi.org/10.1007/s00122-015-2552-7.

22. Šafár, J.; Bartoš, J.; Janda, J.; Bellec, A.; Kubaláková, M.; Valárik, M.; Pateyron, S.; Weiserová, J.; Tušková, R.; Cíhalíková, J.; et al. Dissecting large and complex genomes: Flow-sorting and BAC cloning of individual chromosomes from bread wheat. *Plant J.* **2004**, *39*, 960–968. https://doi.org/10.1111/j.1365-313X.2004.02179.x.

23. Paux, E.; Sourdille, P.; Salse, J.; Saintenac, C.; Choulet, F.; Leroy, P.; Korol, A.; Michalak, M.; Kianian, S.; Spielmeyer, W.; et al. A physical map of the 1-gigabase bread wheat chromosome 3B. *Science* **2008**, *322*, 101–104. https://doi.org/10.1126/science.1161847.

24. Choulet, F.; Alberti, A.; Theil, S.; Glover, N.M.; Barbe, V.; Daron, J.; Pingault, L.; Sourdille, P.; Couloux, A.; Paux, E.; et al. Structural and Functional Partitioning of Bread Wheat Chromosome 3B. *Science* **2014**, *345*, 1249721. https://doi.org/10.1126/science.1249721.

25. Daron, J.; Glover, N.M.; Pingault, L.; Theil, S.; Jamilloux, V.; Paux, E.; Barbe, V.; Mangenot, S.; Alberti, A.; Wincker, P.; et al. Organization and evolution of transposable elements along the bread wheat chromosome 3B. *Genome Biol.* **2014**, *15*, 546. https://doi.org/10.1186/s13059-014-0546-4.
26. Glover, N.M.; Daron, J.; Pingault, L.; Vandepoele, K.; Paux, E.; Feuillet, C.; Choulet, F. Small-scale gene duplications played a major role in the recent evolution of wheat chromosome 3B. *Genome Biol.* **2015**, *16*, 188. https://doi.org/10.1186/s13059-015-0754-6.

27. Pingault, L.; Choulet, F.; Alberti, A.; Glover, N.M.; Wincker, P.; Feuillet, C.; Paux, E. Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome. *Genome Biol.* **2015**, *16*, 29. https://doi.org/10.1186/s13059-015-0601-9.

28. Mello-Sampayo, T. Genetic Regulation of Meiotic Chromosome Pairing by Chromosome 3D of *Triticum aestivum*. *Nat. New Biol.* **1971**, *230*, 22–23. https://doi.org/10.1038/newbio230022a0.

29. Mello-Sampayo, T.; Canas, P. Suppressors of meiotic chromosome pairing in common wheat. In Proceedings of the 4th International Wheat Genetics Symposium, Columbia, MO, USA, 6–11 August 1973; pp. 709–713.

30. Mello-Sampayo, T.; Lorente, R. The role of chromosome 3D in the regulation of meiotic pairing in hexaploid wheat. *EWAC Newsl.* **1969**, *2*, 19–24.

31. Martinez, M.; Cuñado, N.; Carcelén, N.; Romero, C. The *Ph1* and *Ph2* loci play different roles in the synaptic behaviour of hexaploid wheat *Triticum aestivum*. *Theor. Appl. Genet.* **2001**, *103*, 398–405. https://doi.org/10.1007/s00122-001-0543-3.

32. Prieto, P.; Moore, G.; Reader, S. Control of conformation changes associated with homologue recognition during meiosis. *Theor. Appl. Genet.* **2005**, *111*, 505–510. https://doi.org/10.1007/s00122-005-2040-6.

33. Sears, E.R. Genetic Control of Chromosome Pairing in Wheat. *Annu. Rev. Genet.* **1976**, *10*, 31–51.

34. Sears, E.R. An induced mutant with homoeologous pairing in common wheat. *Can. J. Genet. Cytol.* **1977**, *19*, 585–593. https://doi.org/10.1139/g77-063.

35. Sears, E.R. A wheat mutation conditioning an intermediate level of homoeologous chromosome pairing. *Can. J. Genet. Cytol.* **1982**, *24*, 715–719. https://doi.org/10.1139/g82-076.

36. Dong, C.; Whitford, R.; Langridge, P. A DNA mismatch repair gene links to the *Ph2* locus in wheat. *Genome* **2002**, *45*, 116–124. https://doi.org/10.1139/G01-126.

37. Sutton, T.; Whitford, R.; Baumann, U.; Dong, C.; Able, J.A.; Langridge, P. The *Ph2* pairing homoeologous locus of wheat (*Triticum aestivum*): Identification of candidate meiotic genes using a comparative genetics approach. *Plant J.* **2003**, *36*, 443–456. https://doi.org/10.1046/j.1365-313x.2003.01891.x.

38. Serra, H.; Svačina, R.; Baumann, U.; Whitford, R.; Sutton, T.; Bartoš, J.; Sourdille, P. *Ph2* encodes the mismatch repair protein MSH7-3D that inhibits wheat homoeologous recombination. *Nat. Comm.* **2021**, *12*, 803. https://doi.org/10.1038/s41467-021-21127-1.

39. Driscoll, C.J. Genetic suppression of homoeologous chromosome pairing in hexaploid wheat. *Can. J. Genet. Cytol.* **1972**, *14*, 39–42. https://doi.org/10.1139/g72-004.

40. Riley, R.; Chapman, V. Genetic Control of the Cytologically Diploid Behaviour of Hexaploid Wheat. *Nature* **1958**, *182*, 713–715. https://doi.org/10.1038/182713a0.

41. Griffiths, S.; Sharp, R.; Foote, T.N.; Bertin, I.; Wanous, M.; Reader, S.; Colas, I.; Moore, G. Molecular characterization of *Ph1* as a major chromosome pairing locus in polyploid wheat. *Nature* **2006**, *439*, 749–752. https://doi.org/10.1038/nature04434.

42. Martín, A.C.; Rey, M.D.; Shaw, P.; Moore, G. Dual effect of the wheat *Ph1* locus on chromosome synapsis and crossover. *Chromosoma* **2017**, *126*, 669–680. https://doi.org/10.1007/s00412-017-0630-0.

43. Rey, M.D.; Martín, A.C.; Higgins, J.; Swarbreck, D.; Uauy, C.; Shaw, P.; Moore, G. Exploiting the ZIP4 homologue within the wheat *Ph1* locus has identified two lines exhibiting homoeologous crossover in wheat-wild relative hybrids. *Mol Breed.* **2017**, *37*, 95. https://doi.org/10.1007/s11032-017-0700-2.

44. Miller, T.E.; Chapman, V. Aneuhaploids in bread wheat. *Genet. Res.* **1976**, *28*, 37–45. https://doi.org/10.1017/S0016672300016700.

45. Viegas, W.S. The effect of B-chromosomes of rye on the chromosome association in F1 hybrids *Triticum aestivum* × *Secale cereale* in the absence of chromosomes 5B or 5D. *Theor. Appl. Genet.* **1980**, *56*, 193–198. https://doi.org/10.1007/BF00295448.

46. Denis, J.-B.; Bernard, M.; Arnoux, J.; Cauderon, Y. Analyse statistique des configurations méiotiques lors de la création d'une série monosomique du blé tendre " Courtot ". *Agronomie* **1982**, *2*, 701–708. https://doi.org/10.1051/agro:19820803.

47. The International Wheat Genome Sequence Consortium Shifting the limits in wheat research and breeding using a fully annotated reference genome. *Science* **2018**, *361*, eaar7191. https://doi.org/10.1126/science.aar7191.

48. Colas, I.; Darrier, B.; Arrieta, M.; Mittmann, S.U.; Ramsay, L.; Sourdille, P.; Waugh, R. Observation of Extensive Chromosome Axis Remodeling during the "Diffuse-Phase" of Meiosis in Large Genome Cereals. *Front. Plant Sci.* **2017**, *8*, 1235. https://doi.org/10.3389/fpls.2017.01235.

49. Boden, S.A.; Langridge, P.; Spangenberg, G.; Able, J.A. TaASY1 promotes homologous chromosome interactions and is affected by deletion of *Ph1*. *Plant J.* **2009**, *57*, 487–497. https://doi.org/10.1111/j.1365-313X.2008.03701.x.

50. Colas, I.; Shaw, P.; Prieto, P.; Wanous, M.; Spielmeyer, W.; Mago, R.; Moore, G. Effective chromosome pairing requires chromatin remodeling at the onset of meiosis. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 6075–6080. https://doi.org/10.1073/pnas.0801521105.

51. Barakate, A.; Higgins, J.D.; Vivera, S.; Stephens, J.; Perry, R.M.; Ramsay, L.; Colas, I.; Oakey, H.; Waugh, R.; Franklin, F.C.H.; et al. The Synaptonemal Complex Protein ZYP1 Is Required for Imposition of Meiotic Crossovers in Barley. *Plant Cell* **2014**, *26*, 729–740. https://doi.org/10.1105/tpc.113.121269.

52. Voelkel-Meiman, K.; Taylor, L.F.; Mukherjee, P.; Humphreys, N.; Tsubouchi, H.; MacQueen, A.J. SUMO Localizes to the Central Element of Synaptonemal Complex and Is Required for the Full Synapsis of Meiotic Chromosomes in Budding Yeast. *PLoS Genet.* **2013**, *9*, e1003837. https://doi.org/10.1371/journal.pgen.1003837.

53. Heslop-Harrison, J.S. Comparative Genome Organization in Plants: From Sequence and Markers to Chromatin and Chromosomes. *Plant Cell* **2000**, *12*, 617–636. https://doi.org/10.2307/3870990.
54. El-Twab, M.H.A. Physical mapping of the 45S rDNA on the chromosomes of *Triticum turgidum* and *T. aestivum* using fluorescence in situ hybridization for chromosome ancestors. *Arab. J. Biotech.* **2007**, *10*, 69-80.

55. Lloyd, A.; Ranoux, M.; Vautrin, S.; Glover, N.M.; Fourment, J.; Charif, D.; Choulet, F.; Lassalle, G.; Marande, W.; Tran Ngoc Danh, J.; et al. Meiotic gene evolution: Can you teach a new dog new tricks? *Mol. Biol. Evol.* **2014**, *31*, 1724-1727.
https://doi.org/10.1093/molbev/msu119.

56. Crismani, W.; Girard, C.; Froger, N.; Pradillo, M.; Santos, J.L.; Chelysheva, L.; Copenhaver, G.P.; Horlow, C.; Mercier, R. FANCM limits meiotic crossovers. *Science* **2012**, *336*, 1588-1590.

57. Emmanuel, E.; Yehuda, E.; Melamed-Bessudo, C.; Avivi-Ragolsky, N.; Levy, A.A. The role of AtMSH2 in homologous recombination in *Arabidopsis thaliana*. *EMBO Rep.* **2006**, *7*, 100-105. https://doi.org/10.1038/sj.embor.7400577.

58. Tam, S.M.; Hays, J.B.; Chetelat, R.T. Effects of suppressing the DNA mismatch repair system on homeologous recombination in tomato. *Theor. Appl. Genet.* **2011**, *123*, 1445-1458. https://doi.org/10.1007/s00122-011-1679-4.

59. Kempanna, C.; Riley, R. Relationships between the genetic effects of deficiencies for chromosome III and V on meiotic pairing in *Triticum aestivum*. *Nature* **1962**, *195*, 1270-1273. https://doi.org/10.1038/1951270a0.

60. Girard, C.; Crismani, W.; Froger, N.; Mazel, J.; Lemhemdi, A.; Horlow, C.; Mercier, R. FANCM-associated proteins MHF1 and MHF2, but not the other Fanconi anemia factors, limit meiotic crossovers. *Nucleic Acids Res.* **2014**, *42*, 9087-9095.

61. Girard, C.; Chelysheva, L.; Choinard, S.; Froger, N.; Macaisne, N.; Lemhemdi, A.; Mazel, J.; Crismani, W.; Mercier, R. AAA-ATPase FIDGETIN-LIKE 1 and helicase FANCM antagonize meiotic crossovers by distinct mechanisms. *PLoS Genet.* **2015**, *11*, e1005369. https://doi.org/10.1371/journal.pgen.1005369.

62. Séguéla-Arnaud, M.; Crismani, W.; Larchevêque, C.; Mazel, J.; Froger, N.; Choinard, S.; Lemhemdi, A.; Macaisne, N.; van Leene, J.; Gevaert, K.; et al. Multiple mechanisms limit meiotic crossovers: TOP3α and two BLM homologs antagonize crossovers in parallel to FANCM. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, 4713-4718. https://doi.org/10.1073/pnas.1423107112.

63. Séguéla-Arnaud, M.; Choinard, S.; Larchevêque, C.; Girard, C.; Froger, N.; Crismani, W.; Mercier, R. RMI1 and TOP3α limit meiotic CO formation through their C-terminal domains. *Nucleic Acids Res.* **2017**, *45*, 1860-1871.
https://doi.org/10.1093/nar/gkw1210.

64. Fernandes, J.B.; Séguéla-Arnaud, M.; Larchevêque, C.; Lloyd, A.H.; Mercier, R. Unleashing meiotic crossovers in hybrid plants. *Proc. Natl. Acad. Sci. USA* **2018**, *115*, 2431-2436. https://doi.org/10.1073/pnas.1713078114.

65. Mieulet, D.; Aubert, G.; Bres, C.; Klein, A.; Droc, G.; Vieille, E.; Rond-Coissieux, C.; Snachez, M.; Dalmais, M.; Mauxion, J.-P.; et al. Unleashing meiotic crossovers in crops. *Nat. Plant* **2018**, *4*, 1010-1016 https://doi.org/10.1038/s41477-018-0311-x.

66. Bassi, F.M.; Kumar, A.; Zhang, Q.; Paux, E.; Huttner, E.; Kilian, A.; Dizon, R.; Feuillet, C.; Xu, S.S.; Kianian, S.F. Radiation hybrid QTL mapping of *Tdes2* involved in the first meiotic division of wheat. *Theor. Appl. Genet.* **2013**, *126*, 1977-1990.
https://doi.org/10.1007/s00122-013-2111-z.

67. Devos, K.; Sorrells, M.; Anderson, J.; Miller, T.E.; Reader, S.M.; Lukaszewski, A.J.; Dubcovsky, J.; Sharp, P.J.; Faris, J.; Gale, M.D. Chromosome aberrations in wheat nullisomic-tetrasomic and ditelosomic lines. *Cereal Res. Commun.* **1999**, *27*, 231-239.

68. Dvořák, J. Genetic variability in *Aegilops speltoides* affecting homoeologous pairing in wheat. *Can. J. Genet. Cytol.* **1972**, *14*, 371-380. https://doi.org/10.1139/g72-046.

69. Dvorak, J.; Deal, K.R.; Luo, M.-C. Discovery and Mapping of Wheat *Ph1* Suppressors. *Genetics* **2006**, *174*, 17-27.
https://doi.org/10.1534/genetics.106.058115.

70. Beck, M.; Komis, G.; Muller, J.; Menzel, D.; Samaj, J. Arabidopsis Homologs of Nucleus- and Phragmoplast-Localized Kinase 2 and 3 and Mitogen-Activated Protein Kinase 4 Are Essential for Microtubule Organization. *Plant Cell* **2010**, *22*, 755-771.
https://doi.org/10.1105/tpc.109.071746.

71. Zeng, Q.; Chen, J.-G.; Ellis, B.E. AtMPK4 is required for male-specific meiotic cytokinesis in Arabidopsis. *Plant J.* **2011**, *67*, 895-906. https://doi.org/10.1111/j.1365-313X.2011.04642.x.

72. Alou, A.H.; Jean, M.; Domingue, O.; Belzile, F.J. Structure and expression of AtPMS1, the Arabidopsis ortholog of the yeast DNA repair gene PMS1. *Plant Sci.* **2004**, *167*, 447-456. https://doi.org/10.1016/j.plantsci.2004.04.012.

73. Franklin, F.C.H.; Higgins, J.D.; Sanchez-Moran, E.; Armstrong, S.J.; Osman, K.E.; Jackson, N.; Jones, G.H. Control of meiotic recombination in Arabidopsis: Role of the *MutL* and *MutS* homologues. *Biochem. Soc. Trans.* **2006**, *34*, 542-544.
https://doi.org/10.1042/BST0340542.

74. Rogacheva, M.V.; Manhart, C.M.; Chen, C.; Guarne, A.; Surtees, J.; Alani, E. Mlh1-Mlh3, a Meiotic Crossover and DNA Mismatch Repair Factor, Is a Msh2-Msh3-stimulated Endonuclease. *J. Biol. Chem.* **2014**, *289*, 5664-5673.
https://doi.org/10.1074/jbc.M113.534644.

75. Dion, É.; Li, L.; Jean, M.; Belzile, F. An Arabidopsis MLH1 mutant exhibits reproductive defects and reveals a dual role for this gene in mitotic recombination. *Plant J.* **2007**, *51*, 431-440.

76. Martín, A.C.; Shaw, P.; Phillips, D.; Reader, S.; Moore, G. Licensing MLH1 sites for crossover during meiosis. *Nat. Commun.* **2014**, *5*, 4580.

77. Sutani, T.; Yuasa, T.; Tomonaga, T.; Dohmae, N.; Takio, K.; Yanagida, M. Fission yeast condensin complex: Essential roles of non-SMC subunits for condensation and Cdc2 phosphorylation of Cut3/SMC4. *Genes Dev.* **1999**, *13*, 2271-2283.

78. Siddiqui, N.U.; Stronghill, P.E.; Dengler, R.E.; Hasenkampf, C.A.; Riggs, C.D. Mutations in Arabidopsis condensin genes disrupt embryogenesis, meristem organization and segregation of homologous chromosomes during meiosis. *Development* **2003**, *130*, 3283-3295. https://doi.org/10.1242/dev.00542.
79. Smith, S.J.; Osman, K.; Franklin, F.C.H. The condensin complexes play distinct roles to ensure normal chromosome morphogenesis during meiotic division in Arabidopsis. *Plant J.* **2014**, *80*, 255–268.

80. Shen, Y.; Tang, D.; Wang, K.; Wang, M.; Huang, J.; Luo, W.; Luo, Q.; Hong, L.; Li, M.; Cheng, Z. ZIP4 in homologous chromosome synapsis and crossover formation in rice meiosis. *J. Cell Sci.* **2012**, *125*, 2581–2591. https://doi.org/10.1242/jcs.090993

81. Chelysheva, L.; Gendrot, G.; Vezon, D.; Doutriaux, M.P.; Mercier, R.; Grelon, M. Zip4/Spo22 is required for class I CO formation but not for synapsis completion in *Arabidopsis thaliana*. *PLoS Genet.* **2007**, *3*, 802–813. https://doi.org/10.1371/journal.pgen.0030083.

82. Yuan, L.; Yang, X.; Ellis, J.L.; Fisher, N.M.; Makaroff, C.A. The Arabidopsis SYN3 cohesin protein is important for early meiotic events. *Plant J.* **2012**, *71*, 147–160. https://doi.org/10.1111/j.1365-313X.2012.04979.x.

83. Hartung, F.; Suer, S.; Bergmann, T.; Puchta, H. The role of AtMUS81 in DNA repair and its genetic interaction with the helicase AtRecQ4A. *Nucleic Acids Res.* **2006**, *34*, 4438–4448. https://doi.org/10.1093/nar/gkl576.

84. Ramirez-Gonzalez, R.H.; Borrill, P.; Lang, D.; Harrington, S.A.; Brinton, J.; Venturini, L.; Davey, M.; Jacobs, J.; van Ex, F.; Pasha, A.; et al. The transcriptional landscape of polyploid wheat. *Science* **2018**, *361*, eaar6089.

85. Sears, E.R. Nullisomic-Tetrasomic Combinations in Hexaploid Wheat. In *Chromosome Manipulations and Plant Genetics*; Riley, R., Lewis, K.R., Eds; Springer US: Boston, MA, USA, **1966**; pp. 29–45. https://doi.org/10.1007/978-1-4899-6561-5_4.

86. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. *Nat. Biotechnol.* **2010**, *28*, 511–515. https://doi.org/10.1038/nbt.1621.

87. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. *Nat. Methods* **2008**, *5*, 621–628.

88. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. *Nucleic Acids Res.* **1997**, *25*, 3389–3402. https://doi.org/10.1093/nar/25.17.3389.

89. Slater, G.; Birney, E. Automated generation of heuristics for biological sequence comparison. *BMC Bioinform.* **2005**, *6*, 31. https://doi.org/10.1186/1471-2105-6-31.

90. Walkowiak, S.; Gao, L.; Monat, C.; Haberer, G.; Kassa, M.T.; Brinton, J.; Ramirez-Gonzalez, R.H.; Kolodziej, M.C.; Delorean, E.; Thambugala, D.; et al. Multiple wheat genomes reveal global variation in modern breeding. *Nature* **2020**, *588*, 277–283.

91. Colas, I.; Macaulay, M.; Higgins, J.D.; Phillips, D.; Barakate, A.; Posch, M.; Armstrong, S.J.; Franklin, F.C.H.; Halpin, C.; Waugh, R.; et al. A spontaneous mutation in MutL-Homolog 3 (HvMLH3) affects synapsis and crossover resolution in the barley desynaptic mutant *des10*. *New Phytol.* **2016**, *212*, 693–707. https://doi.org/10.1111/nph.14061.

92. Darrier, B.; Arrieta, M.; Mittmann, S.U.; Sourdille, P.; Ramsay, L.; Waugh, R.; Colas, I. Following the Formation of Synaptonemal Complex Formation in Wheat and Barley by High-Resolution Microscopy. *Methods Mol. Biol.* **2020**, *2061*, 207–215. https://doi.org/10.1007/978-1-4939-9818-0_15.

93. Boden, S.A.; Shadiac, N.; Tucker, E.J.; Langridge, P.; Able, J.A. Expression and functional analysis of *TaASY1* during meiosis of bread wheat (*Triticum aestivum*). *BMC Mol. Biol.* **2007**, *8*, 65. https://doi.org/10.1186/1471-2199-8-65.

94. Gerlach, W.L.; Bedbrook, J.R. Cloning and characterization of ribosomal RNA genes from wheat and barley. *Nucleic Acids Res.* **1979**, *7*, 1869–1885. https://doi.org/10.1093/nar/7.7.1869.

95. Aragon-Alcaide, L.; Beven, A.; Moore, G.; Shaw, P. The use of vibratome sections of cereal spikelets to study anther development and meiosis. *Plant J.* **1998**, *14*, 503–508. https://doi.org/10.1046/j.1365-313X.1998.00148.x.

96. Higgins, J.D.; Perry, R.M.; Barakate, A.; Ramsay, L.; Waugh, R.; Halpin, C.; Armstrong, S.J.; Franklin, F.C.H. Spatiotemporal Asymmetry of the Meiotic Program Underlies the Predominantly Distal Distribution of Meiotic Crossovers in Barley. *Plant Cell* **2012**, *24*, 4096–4109. https://doi.org/10.1105/tpc.112.102483.

97. Das, R.M.; Storey, K.G. Apical Abscission Alters Cell Polarity and Dismantles the Primary Cilium During Neurogenesis. *Science* **2014**, *343*, 200–204.

98. Vonesch, C.; Unser, M. A Fast Thresholded Landweber Algorithm for Wavelet-Regularized Multidimensional Deconvolution. *IEEE Trans. Image Process.* **2008**, *17*, 539–549.

99. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saafeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. *Nat. Methods* **2012**, *9*, 676–682.