The first division of meiosis is unique in its capacity to halve the ploidy of future gametes. To this end, one key innovation compared with mitosis is the monopolar orientation of the pairs of sister kinetochores required for the proper separation of homologs at meiosis I. How monopolar orientation is imposed is unclear and seems to vary in eukaryotes. Here, we performed a forward genetic screen in Arabidopsis thaliana, specifically designed to identify the molecular components imposing monopolar orientation, based on mutants’ ability to restore fertility in spo11 osd1 haploid plants. We show that monopolar orientation involves all four cohesin subunits (REC8, SCC3, SMC1, and SMC3), the cohesion establishment factors CTF18 and DCC1, the cohesin protectors SGO1/2 and PATRONUS (PANS1), the inner kinetochore protein CENP-C, and the deSUMOylase SUMO PROTEASE RELATED TO FERTILITY 2 (SPF2). The mutants show bipolar orientation of achiasmatic chromosomes; however, monopolar orientation is maintained in the presence of crossovers, despite most of them displaying the splitting of sister kinetochores and reduced levels of cohesin at metaphase I. Taken together, the findings demonstrate that cohesion establishment and maintenance, kinetochore function, and deSUMOylation, together with crossovers, promote monopolar orientation in plants and support a cohesion-driven model of kinetochore orientation at meiosis I that is conserved across kingdoms.

Background Meiotic recombination creates genetic diversity through reciprocal exchange of haplotypes between homologous chromosomes. Scalable and robust methods for mapping recombination breakpoints are essential for understanding meiosis and for genetic mapping. Single cell sequencing of gametes offers a direct approach to recombination mapping, yet the effect of technical differences between single-cell sequencing methods for crossover detection remains unclear. Results We benchmark single cell methods for droplet-based chromatin accessibility and RNA sequencing and plate-based whole-genome amplification for mapping meiotic recombination in Arabidopsis thaliana . For this purpose we introduce two novel open-source tools coelsch\_mapping\_pipeline and coelsch for haplotype-aware alignment and per-cell crossover detection, using them to recover known recombination frequencies and quantify the effects of coverage sparsity. We subsequently apply our approach to a panel of 40 recombinant F₁ hybrids derived from crosses of 22 diverse natural accessions, successfully recovering genetic maps for 34 F1s in a single dataset. This analysis reveals substantial variation in recombination rate and identifies a ∼10 Mb pericentric inversion in the accession Zin-9, the largest natural inversion reported in A. thaliana to date. Conclusions These results demonstrate the applicability and scalability of single-cell gamete sequencing for high-throughput mapping of meiotic recombination, and highlight the strengths and limitations of different single-cell modalities. The accompanying open-source tools provide a framework for haplotyping and crossover detection analysis using sparse single-cell sequencing data. Our methodology enables parallel analysis of large numbers of hybrids in a single dataset, removing a major technical barrier to large-scale studies of natural variation in recombination rate. ### Competing Interest Statement The authors declare the following competing interests: Takara provided Takara Shasta Whole Genome Amplification kits free of charge. The company had no role in study design, analysis, or publication decisions. A patent was granted to Institut National de la Recherche Agronomique (INRA) on the use of RECQ4 to manipulate meiotic recombination in crops, with RM listed among the inventors (US10,920,237/EP3149027). A patent was filed by the Max Planck Society on the combined use of RECQ4 and ZYP1 to manipulate recombination in crops, with RM listed among the inventors (EP23179262. 14.06.2023).

In eukaryotes, chromosome segregation relies on attachment to the spindle, ensured by the kinetochore. The outer kinetochore attaches to the microtubules and is named after three sub-complexes KNL1C, MIS12C, and NDC80C (KMN). While the KMN complex comprises ten proteins in humans S. cerevisiae, its conservation in more distant eukaryotes is unclear. Here, we aimed to define the KMN complex in the plant Arabidopsis using affinity purification and identified thirteen KMN proteins. Seven were previously known to have a conserved function (atMIS12, atNNF1, atNDC80, atSPC24, atSPC25, atNUF2, and atKNL1) and six were uncharacterized. These six proteins show remote similarity to yeast/human KMN-associated proteins, whose homologs have not yet been characterized in plants. We named them atDSN1, atCSM1, atNSL1.1/.2, and atZWINT1.1/.2. We confirmed kinetochore localization for atDSN1, atCSM1, atNSL1.1, and atZWINT1.1 in planta. In addition, atDSN1, atCSM1, and ZWINT1.1/.2 are essential, further supporting their kinetochore function. AlphaFold3 predicts an alike3D organization of the KMN complex in plants and mammals. We conclude that the KMN complex is globally conserved with a matching composition and similar organization in distant eukaryotes, with some local variations, suggesting its presence in the common ancestors of all living eukaryotes.

The first division of meiosis is unique in its capacity to halve the ploidy of the future gametes. To this end, one key innovation compared to mitosis is the monopolar orientation of the pairs of sister kinetochores required for the proper separation of homologs at meiosis I. How monopolar orientation is imposed is unclear and seems to vary in eukaryotes. Here we performed a forward genetic screen in Arabidopsis thaliana , specifically designed to identify the molecular components imposing monopolar orientation, based on mutants’ ability to restore fertility in spo11 osd1 haploid plants. We show that monopolar orientation involves all four cohesin subunits (REC8, SCC3, SMC1, SMC3), the cohesion establishment factors CTF18 and DCC1, and the cohesin protectors SGO1/2 and PANS1, the inner kinetochore protein CENP-C, and the desumoylase SPF2. The mutants show bipolar orientation of achiasmatic chromosomes, but monopolar orientation is maintained in the presence of crossovers despite splitting of sister kinetochores at metaphase I and reduced levels of cohesin. Taken together, the findings demonstrate that cohesion establishment and protection, kinetochore function, and deSUMOylation, together with crossovers, enforce monopolar orientation in plants and support a cohesion-driven model of kinetochore orientation at meiosis I that is conserved across kingdoms. ### Competing Interest Statement The authors have declared no competing interest.

Potatoes were first brought to Europe in the sixteenth century1,2. Two hundred years later, one of the species had become one of the most important food sources across the entire continent and, later, even the entire world3. However, its highly heterozygous, autotetraploid genome has complicated its improvement since then4-7. Here we present the pan-genome of European potatoes generated from phased genome assemblies of ten historical potato cultivars, which includes approximately 85% of all haplotypes segregating in Europe. Sequence diversity between the haplotypes was extremely high (for example, 20× higher than in humans), owing to numerous introgressions from wild potato species. By contrast, haplotype diversity was very low, in agreement with the population bottlenecks caused by domestication and transition to Europe. To illustrate a practical application of the pan-genome, we converted it into a haplotype graph and used it to generate phased, megabase-scale pseudo-genome assemblies of commercial potatoes (including the famous French fries potato 'Russet Burbank') using cost-efficient short reads only. In summary, we present a nearly complete pan-genome of autotetraploid European potato, we describe extraordinarily high sequence diversity in a domesticated crop, and we outline how this resource might be used to accelerate genomics-assisted breeding and research.

Background Plant meristems are structured organs consisting of distinct layers of stem cells, which differentiate into new plant tissue. Mutations in meristematic layers can propagate into large sectors of the plant. However, the characteristics of meristematic mutations remain unclear, limiting our understanding of the genetic basis of somaclonal phenotypic variation. Results Here, we analyse the frequency and distribution of somatic mutations in an apricot tree. We separately sequence the epidermis (developing from meristem layer 1) and the flesh (developing from meristem layer 2) of several fruits sampled across the entire tree. We find that most somatic mutations (> 90%) are specific to individual layers. Interestingly, layer 1 shows a higher mutation load than layer 2, implying different mutational dynamics between the layers. The distribution of somatic mutations follows the branching of the tree. This suggests that somatic mutations are propagated to developing branches through axillary meristems. In turn, this leads us to the unexpected observation that the genomes of layer 1 of distant branches are more similar to each other than to the genomes of layer 2 of the same branches. Finally, using single-cell RNA sequencing, we demonstrate that layer-specific mutations were only transcribed in the cells of the respective layers and can form the genetic basis of somaclonal phenotypic variation. Conclusions Here, we analyse the frequency and distribution of somatic mutations with meristematic origin. Our observations on the layer specificity of somatic mutations outline how they are distributed, how they propagate, and how they can impact clonally propagated crops.

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Crossovers (CO) shuffle genetic information and physically connect homologous chromosomal pairs, ensuring their balanced segregation during meiosis. COs arising from the major class I pathway require the activity of the well-conserved group of ZMM proteins, which, in conjunction with MLH1, facilitate the maturation of DNA recombination intermediates specifically into COs. The HEI10 Interacting Protein 1 (HEIP1) was identified in rice and proposed to be a new, plant-specific member of the ZMM group. Here, we establish and decipher the function of the Arabidopsis thaliana HEIP1 homolog in meiotic crossover formation and report its wide conservation in eukaryotes. We show that the loss of Arabidopsis HEIP1 elicits a marked reduction in meiotic COs and their redistribution toward chromosome ends. Epistasis analysis showed that AtHEIP1 acts specifically in the class I CO pathway. Further, we show that HEIP1 acts both prior to crossover designation, as the number of MLH1 foci is reduced in heip1 , and at the maturation step of MLH1-marked sites into COs. Despite the HEIP1 protein being predicted to be primarily unstructured and very divergent at the sequence level, we identified homologs of HEIP1 in an extensive range of eukaryotes, including mammals.

Abstract Although originally primarily a model for functional biology, the thale cress Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse abiotic and biotic environments, developed into a powerful model in population genomics. Here, we present chromosome-level genome assemblies of 72 diverse accessions located across the global species range. We found that genomic co-linearity is remarkably conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare, typically present only in a single accession. This indicates that the highly re-arranged karyotype of Arabidopsis thaliana (as compared to other Brassicaceae species), is quasi-fixed and that rearrangements in chromosome arms are counter-selected. In contrast, centromeric regions show much higher structural dynamics, and divergences in core centromeres account for most of the variation in genome size. Pan-genome analyses of the gene space uncovered 40,382 distinct genes, 52% being present in all accessions and 48% appearing to be dispensable, including 23% private to a single accession, indicating a vast unexplored genic diversity. The availability of 72 new Arabidopsis thaliana genome assemblies will empower future genetic research.

At meiosis, programmed meiotic DNA double-strand breaks are repaired via homologous recombination, resulting in crossovers (COs). From a large excess of DNA double-strand breaks that are formed, only a small proportion gets converted into COs because of active mechanisms that restrict CO formation. The Fanconi anemia (FA) complex proteins AtFANCM, MHF1, and MHF2 were previously identified in a genetic screen as anti-CO factors that function during meiosis in Arabidopsis thaliana . Here, pursuing the same screen, we identify FANCC as a new anti-CO gene. FANCC was previously only identified in mammals because of low primary sequence conservation. We show that FANCC, and its physical interaction with FANCE-FANCF, is conserved from vertebrates to plants. Further, we show that FANCC, together with its subcomplex partners FANCE and FANCF, regulates meiotic recombination. Mutations of any of these three genes partially rescues CO-defective mutants, which is particularly marked in female meiosis. Functional loss of FANCC, FANCE, or FANCF results in synthetic meiotic catastrophe with the pro-CO factor MUS81. This work reveals that FANCC is conserved outside mammals and has an anti-CO role during meiosis together with FANCE and FANCF.

Meiotic recombination frequency varies along chromosomes and strongly correlates with sequence divergence. However, the causal relationship between recombination landscapes and polymorphisms is unclear. Here, we characterize the genome-wide recombination landscape in the quasi-absence of polymorphisms, using Arabidopsis thaliana homozygous inbred lines in which a few hundred genetic markers were introduced through mutagenesis. We find that megabase-scale recombination landscapes in inbred lines are strikingly similar to the recombination landscapes in hybrids, with the notable exception of heterozygous large rearrangements where recombination is prevented locally. In addition, the megabase-scale recombination landscape can be largely explained by chromatin features. Our results show that polymorphisms are not a major determinant of the shape of the megabase-scale recombination landscape but rather favour alternative models in which recombination and chromatin shape sequence divergence across the genome.

The oomycete Phytophthora infestans causes late blight of potato, which can completely destroy the crop. Therefore, for the past 160 years, late blight has been the most important potato disease worldwide. The identification of cultivars with high and durable field resistance to P. infestans is an objective of most potato breeding programs. This type of resistance is polygenic and therefore quantitative. Its evaluation requires multi-year and location trials. Furthermore, quantitative resistance to late blight correlates with late plant maturity, a negative agricultural trait. Knowledge of the molecular genetic basis of quantitative resistance to late blight not compromised by late maturity is very limited. It is however essential for developing diagnostic DNA markers that facilitate the efficient combination of superior resistance alleles in improved cultivars. We used association genetics in a population of 184 tetraploid potato cultivars in order to identify single nucleotide polymorphisms (SNPs) that are associated with maturity corrected resistance (MCR) to late blight. The population was genotyped for almost 9000 SNPs from three different sources. The first source was candidate genes specifically selected for their function in the jasmonate pathway. The second source was novel candidate genes selected based on comparative transcript profiling (RNA-Seq) of groups of genotypes with contrasting levels of quantitative resistance to P. infestans. The third source was the first generation 8.3k SolCAP SNP genotyping array available in potato for genome wide association studies (GWAS). Twenty seven SNPs from all three sources showed robust association with MCR. Some of those were located in genes that are strong candidates for directly controlling quantitative resistance, based on functional annotation. Most important were: a lipoxygenase (jasmonate pathway), a 3-hydroxy-3-methylglutaryl coenzyme A reductase (mevalonate pathway), a P450 protein (terpene biosynthesis), a transcription factor and a homolog of a major gene for resistance to P. infestans from the wild potato species Solanum venturii. The candidate gene approach and GWAS complemented each other as they identified different genes. The results of this study provide new insight in the molecular genetic basis of quantitative resistance in potato and a toolbox of diagnostic SNP markers for breeding applications.

BACKGROUND:The soil borne, obligate biotrophic fungus Synchytrium endobioticum causes tumor-like tissue proliferation (wart) in potato tubers and thereby considerable crop damage. Chemical control is not effective and unfriendly to the environment. S. endobioticum is therefore a quarantined pathogen. The emergence of new pathotypes of the fungus aggravate this agricultural problem. The best control of wart disease is the cultivation of resistant varieties. Phenotypic screening for resistant cultivars is however time, labor and material intensive. Breeding for resistance would therefore greatly benefit from diagnostic DNA markers that can be applied early in the breeding cycle. The prerequisite for the development of diagnostic DNA markers is the genetic dissection of the factors that control resistance to S. endobioticum in various genetic backgrounds of potato.RESULTS:Progeny of a cross between a wart resistant and a susceptible tetraploid breeding clone was evaluated for resistance to S. endobioticum pathotypes 1, 2, 6 and 18 most relevant in Europe. The same progeny was genotyped with 195 microsatellite and 8303 single nucleotide polymorphism (SNP) markers. Linkage analysis identified the multi-allelic locus Sen1/RSe-XIa on potato chromosome XI as major factor for resistance to all four S. endobioticum pathotypes. Six additional, independent modifier loci had smaller effects on wart resistance. Combinations of markers linked to Sen1/RSe-XIa resistance alleles with one to two additional markers were sufficient for obtaining high levels of resistance to S. endobioticum pathotypes 1, 2, 6 and 18 in the analyzed genetic background.CONCLUSIONS:Potato resistance to S. endobioticum is oligogenic with one major and several minor resistance loci. It is composed of multiple alleles for resistance and susceptibility that originate from multiple sources. The genetics of resistance to S. endobioticum varies therefore between different genetic backgrounds. The DNA markers described in this paper are the starting point for pedigree based selection of cultivars with high levels of resistance to S. endobioticum pathotypes 1, 2, 6 and 18.

Starch accumulation and breakdown are vital processes in plant storage organs such as seeds, roots, and tubers. In tubers of potato (Solanum tuberosum L.) a small fraction of starch is converted into the reducing sugars glucose and fructose. Reducing sugars accumulate in response to cold temperatures. Even small quantities of reducing sugars affect negatively the quality of processed products such as chips and French fries. Tuber starch and sugar content are inversely correlated complex traits that are controlled by multiple genetic and environmental factors. Based on in silico annotation of the potato genome sequence, 123 loci are involved in starch-sugar interconversion, approximately half of which have been previously cloned and characterized. By means of candidate gene association mapping, we identified single-nucleotide polymorphisms (SNPs) in eight genes known to have key functions in starch-sugar interconversion, which were diagnostic for increased tuber starch and/or decreased sugar content and vice versa. Most positive or negative effects of SNPs on tuber-reducing sugar content were reproducible in two different collections of potato cultivars. The diagnostic SNP markers are useful for breeding applications. An allele of the plastidic starch phosphorylase PHO1a associated with increased tuber starch content was cloned as full-length cDNA and characterized. The PHO1a-H-A allele has several amino acid changes, one of which is unique among all known starch/glycogen phosphorylases. This mutation might cause reduced enzyme activity due to impaired formation of the active dimers, thereby limiting starch breakdown.

Resistance to pathogens is essential for survival of wild and cultivated plants. Pathogen susceptibility causes major losses of crop yield and quality. Durable field resistance combined with high yield and other superior agronomic characters are therefore important objectives in every crop breeding program. Precision and efficacy of resistance breeding can be enhanced by molecular diagnostic tools, which result from knowledge of the molecular basis of resistance and susceptibility. Breeding uses resistance conferred by single R genes and polygenic quantitative resistance. The latter is partial but considered more durable. Molecular mechanisms of plant pathogen interactions are elucidated mainly in experimental systems involving single R genes, whereas most genes important for quantitative resistance in crops like potato are unknown. Quantitative resistance of potato to Phytophthora infestans causing late blight is often compromised by late plant maturity, a negative agronomic character. Our objective was to identify candidate genes for quantitative resistance to late blight not compromised by late plant maturity. We used diagnostic DNA-markers to select plants with different field levels of maturity corrected resistance (MCR) to late blight and compared their leaf transcriptomes before and after infection with P. infestans using SuperSAGE (serial analysis of gene expression) technology and next generation sequencing. We identified 2034 transcripts up or down regulated upon infection, including a homolog of the kiwi fruit allergen kiwellin. 806 transcripts showed differential expression between groups of genotypes with contrasting MCR levels. The observed expression patterns suggest that MCR is in part controlled by differential transcript levels in uninfected plants. Functional annotation suggests that, besides biotic and abiotic stress responses, general cellular processes such as photosynthesis, protein biosynthesis and degradation play a role in MCR.

Tuber yield, starch content, starch yield and chip color are complex traits that are important for industrial uses and food processing of potato. Chip color depends on the quantity of reducing sugars glucose and fructose in the tubers, which are generated by starch degradation. Reducing sugars accumulate when tubers are stored at low temperatures. Early and efficient selection of cultivars with superior yield, starch yield and chip color is hampered by the fact that reliable phenotypic selection requires multiple year and location trials. Application of DNA-based markers early in the breeding cycle, which are diagnostic for superior alleles of genes that control natural variation of tuber quality, will reduce the number of clones to be evaluated in field trials. Association mapping using genes functional in carbohydrate metabolism as markers has discovered alleles of invertases and starch phosphorylases that are associated with tuber quality traits. Here, we report on new DNA variants at loci encoding ADP-glucose pyrophosphorylase and the invertase Pain-1, which are associated with positive or negative effect with chip color, tuber starch content and starch yield. Marker-assisted selection (MAS) and marker validation were performed in tetraploid breeding populations, using various combinations of 11 allele-specific markers associated with tuber quality traits. To facilitate MAS, user-friendly PCR assays were developed for specific candidate gene alleles. In a multi-parental population of advanced breeding clones, genotypes were selected for having different combinations of five positive and the corresponding negative marker alleles. Genotypes combining five positive marker alleles performed on average better than genotypes with four negative alleles and one positive allele. When tested individually, seven of eight markers showed an effect on at least one quality trait. The direction of effect was as expected. Combinations of two to three marker alleles were identified that significantly improved average chip quality after cold storage and tuber starch content. In F1 progeny of a single-cross combination, MAS with six markers did not give the expected result. Reasons and implications for MAS in potato are discussed.

The oomycete Phytophthora infestans causes late blight, the most relevant disease of potato (Solanum tuberosum) worldwide. Field resistance to late blight is a complex trait. When potatoes are cultivated under long day conditions in temperate climates, this resistance is correlated with late plant maturity, an undesirable characteristic. Identification of natural gene variation underlying late blight resistance not compromised by late maturity will facilitate the selection of resistant cultivars and give new insight in the mechanisms controlling quantitative pathogen resistance. We tested 24 candidate loci for association with field resistance to late blight and plant maturity in a population of 184 tetraploid potato individuals. The individuals were genotyped for 230 single nucleotide polymorphisms (SNPs) and 166 microsatellite alleles. For association analysis we used a mixed model, taking into account population structure, kinship, allele substitution and interaction effects of the marker alleles at a locus with four allele doses. Nine SNPs were associated with maturity corrected resistance (P < 0.001), which collectively explained 50% of the genetic variance of this trait. A major association was found at the StAOS2 locus encoding allene oxide synthase 2, a key enzyme in the biosynthesis of jasmonates, plant hormones that function in defense signaling. This finding supports StAOS2 as being one of the factors controlling natural variation of pathogen resistance.

Genetic diversity of crop plants resulting from breeding and selection is preserved in gene banks. Utilization of such materials for further crop improvement depends on knowledge of agronomic performance and useful traits, which is usually obtained by phenotypic evaluation. Associations between DNA markers and agronomic characters in collections of crop plants would (i) allow assessment of the genetic potential of specific genotypes prior to phenotypic evaluation, (ii) identify superior trait alleles in germplasm collections, (iii) facilitate high resolution QTL mapping and (iv) validate candidate genes responsible for quantitative agronomic characters. The feasibility of association mapping was tested in a gene bank collection of 600 potato cultivars bred between 1850 and 1990 in different countries. The cultivars were genotyped with five DNA markers linked to previously mapped QTL for resistance to late blight and plant maturity. Specific DNA fragments were tested for association with these quantitative characters based on passport evaluation data. Highly significant association with QTL for resistance to late blight and plant maturity was detected with PCR markers specific for R1, a major gene for resistance to late blight, and anonymous PCR markers flanking the R1 locus at 0.2 Centimorgan genetic distance. The marker alleles associated with increased resistance and later plant maturity were traced to an introgression from the wild species S. demissum. These DNA markers are the first marker that are diagnostic for quantitative agronomic characters in a large collection of cultivars.

A candidate gene approach has been used as a first step to identify the molecular basis of quantitative trait variation in potato. Sugar content of tubers upon cold storage was the model trait chosen because the metabolic pathways involved in starch and sugar metabolism are well known and many of the genes have been cloned. Tubers of two F1 populations of diploid potato grown in six environments were evaluated for sugar content after cold storage. The populations were genotyped with RFLP, AFLP, and candidate gene markers. QTL analysis revealed that QTL for glucose, fructose, and sucrose content were located on all potato chromosomes. Most QTL for glucose content mapped to the same positions as QTL for fructose content. QTL explaining >10% of the variability for reducing sugars were located on linkage groups I, III, VII, VIII, IX, and XI. QTL consistent across populations and/or environments were identified. QTL were linked to genes encoding invertase, sucrose synthase 3, sucrose phosphate synthase, ADP-glucose pyrophospholase, sucrose transporter 1, and a putative sucrose sensor. The results suggest that allelic variants of enzymes operating in carbohydrate metabolic pathways contribute to the genetic variation in cold sweetening.