FR 3084374 A1 20200131 FR 1857115 A 20180730 L'invention concerne un procédé de contrôle-qualité dans le domaine des semences et de la pureté variétale. La commercialisation de semences est soumise au contrôle du taux de pureté de celles-ci. Ce taux est propre à chaque espèce mais doit être de 98% du poids ou plus (directive 66/402/CEE concernant la commercialisation des semences de céréales), cette norme s'applique également aux semences qui sont commercialisées pour la production de semences de bases, pré-base, la production de semences certifiées ou la production d'hybrides. Cette pureté variétale est principalement contrôlée par inspection sur pied, dans le cas de production de semences hybrides avec un parent géniteur male stérile, le taux de pureté de ce parent doit être encore supérieur (99,9 % pour le maïs). La disponibilité d'une solution de contrôle de la qualité alternative au contrôle sur pied est d'intérêt pour les sociétés semencières, notamment par la nécessité d'avoir une évaluation rapide, sans attendre le développement des plantes nécessaire à une l'évaluation phénotypique. Par ailleurs pour ces sociétés, le contrôle de la pureté variétale n'est pas limité aux stades cités ci-dessus, chaque étape en amont de la production de semences de base est concernée par cette exigence de pureté variétale. Il est rappelé que le taux de pureté variétale se définit comme le pourcentage de plantes provenant d'un lot et qui sont conformes à la description de la variété. Ce pourcentage s'exprime en poids de graines. En production de semences hybrides, l'amélioration de la qualité de la production de semences agricoles passe par la vérification de la pureté génétique des lots de semences de base (lignées parentales utilisées pour la production d'hybrides) utilisés dans la production des semences commerciales. Cette pureté est évaluée par la détection et l'identification de grains contaminants dans un échantillon de lot de semence des géniteurs. Les contaminants sont des graines de la même espèce, mais présentant des variations génétiques à certains /oc/ de leur génome, par rapport au génotype 35 attendu pour les graines du lot considéré. Dans le processus de production de lots de semences, la présence de contaminants est réduite, de par la vigilance dans 2 les étapes de production en amont, pratiques culturales, épuration, isolement, et les contrôles réalisés tout au long du procédé. Ainsi, la presque-totalité des graines du lot présente le même génotype, les contaminants étant présents à un pourcentage généralement faible et en effet le taux toléré dans un lot pour qu'il 5 puisse être commercialisé devra être inférieur à 2%. L'identification de traits génétiques d'intérêt est également importante en commercialisation de semences, en effet certains traits assurant par exemple la tolérance à un herbicide ou à un pathogène (par exemple le Mildiou en Tournesol) apportent une valeur ajoutée certaine à un lot de semence et lorsqu'une variété 10 sera commercialisée comme porteuse de ce trait, une vérification de la présence de ce trait dans le lot de semence sera intéressante. Par trait on entend forme allélique d'un lad lié à un caractère phénotypique. Une problématique similaire porte sur la présence fortuite d'OGM ou de toute autre altération dans le génome. La commercialisation de plante non OGM 15 nécessite de faire la preuve de l'absence d'OGM ou de la présence d'un taux inférieur à un pourcentage déterminé par la réglementation. A l'opposé la réglementation dans certains pays, pour certains traits OGM, résistance à l'encontre d'insectes notamment, prévoit que les semences contenant l'OGM soient vendues avec un certain taux de semences ne possédant pas le trait OGM, 20 de façon à assurer des zones refuges pour l'insecte. Le développement massif de marqueurs de type SNP (Single Nucleotide Polymorphism) et des technologies de génotypage à haut débit a permis de favoriser l'essor de la sélection assistée par marqueurs. Le génotypage est 25 classiquement réalisé au moyen de différentes technologies, par PCR (Kasp - LGC Genomics, Taqman - Life Technologies) ou hybridation sur des puces à ADN (Axiom - Life Technologies, Infinium - Illumina). Si la technologie de PCR quantitative Taqman est aujourd'hui considérée comme la référence pour la détection de présence fortuite de plantes OGM dans 30 un mélange de plantes non OGM, celle-ci repose sur la détection d'un polymorphisme de type présence/absence d'une séquence donnée, et non sur un polymorphisme entre différentes formes alléliques d'un SNP. Ainsi, dans ce cas particulier de détection d'OGM, le polymorphisme porte sur la présence d'un trait qui pourra être amplifié (amplicon) et donc facilement identifiable. L'estimation de la pureté de lots de graines, entendue comme l'absence de trait OGM, a été travaillée par Remund (Seed Science Research (2001) 11, 101- 3 119), deux solutions ont été identifiées par ces auteurs pour limiter les ressources nécessaires à ces vérifications et notamment l'analyse en pool. Ils indiquent que cette méthode est efficace lorsque l'on recherche l'absence d'un individu particulier, par contre quand un taux de pureté est recherché il est préférable de 5 travailler graine à graine. Ces auteurs ont élaboré un outil Seedcalc, qui permet notamment une approche quantitative en jouant sur le nombre de pools et le nombre de graines par lot, cette méthode est notamment adaptée pour de la PCR temps réel (Laffont, Seed Science Research (2005) 15, 197-204). Un exemple d'utilisation de pools de graine pour vérifier de la pureté 10 variétale existe cependant. La demande WO 2015/110472 propose d'analyser des lots de graines par prélèvement manuel ou semi-automatique d'un volume d'échantillon déterminé à partir d'une ou plusieurs graines, ce volume étant déterminé pour permettre l'analyse d'au moins un constituant de la graine ou des graines. Le tissu prélèvé à partir de plusieurs graines est placé dans un puits 15 identifié et traçable, puis on effectue l'analyse dudit constituant sur le contenu du ou des puits. Cette méthode de constitution de bulk permet de faire de la pureté variétale (exemple 6) cette pureté est évaluée par la méthode Kaspar (KBioscience) à partir de bulks de 5 et 10 graines, la présence d'un contaminant dans ces bulks se caractérise par la présence d'un cluster hétérozygote, 20 cependant les auteurs indiquent que ce cluster est proche du cluster homozygote et qu'il est plus facile à identifier pour un bulk de 5 graines que pour un bulk de 10 graines. Le développement des technologies de séquençage à haut débit, ou NGS (Next Generation Sequencing) a révolutionné le monde de la génomique, 25 permettant la découverte massive de marqueurs SNP entre lignées d'une espèce donnée. Ces techniques permettent un nombre important de lectures de séquences possibles en une seule expérience. La profondeur de séquençage, permet d'identifier un allèle faiblement représenté lors de l'identification de formes alléliques pour un ensemble d'individus 30 en pool. Elle peut permettre par ailleurs d'identifier un nombre de formes alléliques supérieures à deux pour un même locus. Ainsi, le séquençage d'amplicons permet d'étudier de manière ciblée des /oc/ d'intérêt, d'identifier des SNP et de caractériser la composition allélique d'un individu ou d'un mélange d'individus. Une application en recherche est la détection de mutations rares au sein d'une population 35 mutagénéisée (TILLING, Targeting Induced Local Lesions in Genomes). Dans ces applications l'identification d'allèles rares en pool peut être combinée avec des 4 pools d'individus en 2D ou 3D permettant une diminution du nombre de pools à analyser (Tsai et al, Plant Physiol. 2011 Ju1;156(3)1257-68 ; Taheri et al, Mol Breeding (2017) 37:40 ; Gupta et al, The Plant Journal (2017) 92, 495-508) W02014134729, EP 2 200 424). Cette approche peut être également appliquée à 5 l'identification de mutations par des méthodes de Gene Editing (Kumar et al, Mol Breeding (2017) 37:14). Ces approches restent toutefois qualitatives, il n'y a pas de considération quantitative. La possibilité d'utiliser du génotypage par séquençage en pool a été 10 évaluée pour l'identification de fréquences alléliques sur des populations par Gautier (Mol Ecol. 2013 Ju1;22(14):3766-79). Toutefois cette approche est particulièrement adaptée à l'analyse de populations larges sur un grand nombre de SNP, et ne semble pas adaptée à la détection d'allèles rares (généralement inférieur à 5%). 15 Une des difficultés liées à la recherche d'allèle rare est la fiabilité du résultat, la fréquence de l'allèle rare s'approchant du taux d'erreur de séquençage. Dans le cas du contrôle qualité de lots de semences, l'objectif est de détecter la présence d'un contaminant, d'en estimer précisément le taux au sein du 20 lot de semences dont est issu l'échantillon analysé, et de préférence d'en déterminer le profil génétique pour mieux en comprendre l'origine. La détection peut être réalisée par l'analyse de /oc/ d'intérêt, choisis par l'homme du métier, en fonction de sa connaissance du matériel génétique à qualifier et du matériel génétique susceptible de le contaminer. 25 Ainsi, Chen et al (2016, PLOS ONE 11(6)) ont développé, pour le maïs, deux séries de SNP pour le contrôle qualité : un set de marqueurs pour un contrôle rapide, employant un nombre réduit de SNP (50-100) pour identifier les erreurs potentielles d'étiquetage des paquets ou des parcelles de semences, et un set de marqueurs plus large, et utilisé à des fins de caractérisation et de discrimination 30 plus fine du matériel génétique. Dans cet exemple, l'échantillonnage de 192 individus analysés individuellement permettrait d'avoir une probabilité proche de 100% de détecter une contamination de 5 % dans un lot, mais cette probabilité devient inférieure à 90% si l'on s'intéresse à une contamination de 1 %. Dans le cas du contrôle qualité de lots de semences de base, la pureté génétique 35 attendue est élevée, ainsi que la précision d'estimation recherchée, qui dépend à la fois du nombre de graines échantillonnées (testés) et du nombre de graines du 5 lot de semences de bases. Par exemple, si 200 grains sont analysés et que le taux d'impureté est de 0%, l'intervalle de confiance de cette valeur s'étend de 0% à 1.49%. L'effectif analysé est donc trop faible pour garantir un taux de pureté suffisant en analysant seulement 200 grains. En revanche, lors de l'analyse de 5 2000 grains, un taux de 0% d'impureté a un intervalle de confiance de 0% à 0.15%. Toutefois, même si les coûts de génotypage ont considérablement baissé, un tel échantillonnage, associé à du traitement plante à plante, n'est pas économiquement viable pour du contrôle qualité. La société Genia (Montevideo, Uruguay) propose une méthode de 10 détermination de la pureté génétique sur des lots de lignées, et d'identification des contaminants, par l'analyse d'un mélange unique de 10 000 graines et séquençage d'amplicons ciblant environ 350 SNP. Cette société revendique le fait de déterminer la pureté variétale avec une sensibilité de 0,8 % et un intervalle de confiance de 99%. Cette approche est similaire à celle développée par Gautier et 15 al., en ce qu'elle repose sur un modèle statistique d'estimation des fréquences alléliques sur un nombre important (350) de SNP, à partir de laquelle est réalisée une estimation de la fréquence des différents profils génétiques présents dans le mélange. Toutefois, une telle approche ne permet pas de détecter avec fiabilité un allèle rare pour un SNP donné, ce qui est nécessaire dans la recherche d'une 20 contamination pour un trait donné. Il convient donc de disposer d'une méthode économique, permettant l'analyse d'un nombre important d'individus, afin de déterminer précisément la pureté génétique d'un lot de semence donné et ce notamment pour des lots de 25 semences ayant un taux de pureté élevé. La méthode présentée ici repose sur l'estimation de la pureté d'un lot de semences à partir de l'analyse qualitative binaire (présence/absence d'un contaminant) de plusieurs sous-lots d'échantillons. L'analyse sur chaque sous-lot 30 consiste à détecter la présence d'un allèle alternatif à un ou plusieurs loci d'intérêt, par séquençage d'amplicons. Le nombre de sous-lots, ainsi que la taille de chaque sous-lot sont définis en fonction du taux de pureté attendu (estimé par l'opérateur) et de la précision recherchée, et de manière à ce que l'on ait préférentiellement une probabilité statistique de trouver au maximum un contaminant dans un sous-lot 35 donné. Cela signifie que, à partir, d'un nombre donné de graines que l'utilise pour le test, on forme au moins autant de sous-lots que le nombre de contaminants 6 estimé, préférentiellement exactement autant de sous-lots que le nombre de contaminants estimé. Par ailleurs, du fait de l'analyse de plusieurs sous-lots, la méthode permet de distinguer une contamination par un hybride (ségrégation) et une contamination par une lignée (pas de ségrégation), en comparant les profils 5 contaminants des différents sous-lots. Cependant cette méthode ne se limite pas à cette approche binaire, en effet l'utilisation du séquençage permet de ne pas limiter la méthode à l'identification de deux formes alléliques et dans ce contexte la méthode permet également d'identifier des contaminants dans des lots de semences hétérozygote 10 pour l'allèle considéré, le contaminant étant hétérologue aux formes alléliques de cet individu. L'invention se rapporte ainsi à un procédé de détermination de la quantité de contaminants à au moins un locus d'intérêt, présents dans un lot de semences 15 d'une variété d'intérêt, caractérisé en ce a) on regroupe des graines d'un lot de semences par sous-lots d'au moins 10 graines, le nombre sous-lots ainsi obtenus étant supérieur ou égal à 10 b) on effectue, pour chaque sous-lot, un séquençage ciblé d'au moins la 20 région du génome des graines contenant le locus d'intérêt, c) on détermine, pour chaque sous-lot, la présence d'un contaminant de façon qualitative en cas de détection d'un allèle alternatif à ou aux l'allèle(s) attendus (il peut y avoir plusieurs allèles attendus à un seul locus, notamment si les semences sont des semences d'une plante 25 hybride) pour chaque région génomique séquencée (présence / absence d'un allèle alternatif) d) on détermine la quantité de contaminants dans le lot global par la compilation des résultats qualitatifs obtenus pour l'ensemble des sous-lots. 30 De façon optionnelle, et pour réaliser le séquençage, on amplifie par PCR la région correspondant au locus d'intérêt entre l'étape a) et l'étape b). Cette étape d'amplification est effectuée directement sur l'ensemble des graines dans chaque sous-lot. L'estimation de l'impureté du lot est obtenue suivant la formule : 7 d dans laquelle n est le nombre de pools ; m est le nombre de grains dans un pool ; d est le nombre de pools dans lesquels un contaminant a été identifié. Cette formule est la formule proposée par Remund (2001, op. cit.), qui 5 permet notamment de prendre en compte le fait que les recherches de contaminants sont effectuées uniquement sur un échantillon du lot de semences et donc de prendre en compte les biais potentiellement induits par cet échantillonnage. Ce procédé permet donc de calculer le pourcentage de contaminants dans 10 le lot de semences (et donc la pureté du lot de semences : 1- P ). Un contaminant est une graine présentant un allèle différent de l'allèle attendu au locus d'intérêt donné dans ce lot de semences. Toutefois, lorsque l'on met en oeuvre la méthode sur plusieurs loci d'intérêt, on peut décider que l'on a contamination d'un lot que lorsque l'on observe, dans ce lot, des allèles non- 15 attendus à plus d'un locus, par exemple à 2 ou 3 loci. De façon préférée, à l'étape a), on utilise un nombre maximum de graines, calculé de telle sorte qu'au plus un seul contaminant soit présent dans chaque échantillon (sous-lot) de graines, d'un point de vue statistique. Dans les méthodes 20 de production industrielles, on observe généralement un taux de pureté supérieur à 99 %. Ainsi, avec un effectif d'environ 100 graines, par exemple entre 80 et 120, l'on peut s'attendre à détecter majoritairement une graine contaminante. Les méthodes décrites ci-dessus sont en effet mises en oeuvre pour des lots de semences homogènes, c'est-à-dire pour lesquels au moins 95 %, de préférence au 25 moins 96 %, de façon plus préférée au moins 97 % de façon encore plus préférée au moins 98 %, de façon la plus préférée au moins 99 % des semences ont le même génotype. Selon la pureté estimée du lot de semence, les sous-lots contiennent au maximum 20, ou au maximum 50, ou au maximum 80, ou au maximum 100, voire au maximum 200, ou 2000 graines. 30 L'étape b) du procédé consiste en le séquençage ciblé d'au moins une région génomique, contenant le locus d'intérêt pour lequel on cherche la présence d'un contaminant. Il est clair que cette étape de séquençage est effectuée sur de l'acide nucléique. Ainsi, on prépare l'ADN des lots, par exemple en écrasant les graines et 8 utilisant la farine ou isolant l'ADN à partir de cette farine. Ces méthodes sont connues dans l'art. Cette étape de séquençage est préférentiellement effectuée par séquençage à haut débit (NGS). Différentes technologies (Illumina®, Roche 454, 5 Ion torrent: Proton / PGM (ThermoFisher) ou SOLiD(Applied BioSystems)). En résumé, ces technologies de NGS présentent 2 étapes communes : - une étape d'amplification, par PCR une étape de séquençage, cette étape étant réalisée par des approches différentes en fonction de la technologie utilisée. 10 La technologie Illumina®, utilise l'amplification clonale et le séquençage par synthèse (SBS). On génère une banque d'ADN double brin à partir de l'échantillon à analyser par amplification PCR et ajout d'adaptateurs spécifiques aux extrémités, puis on dénature l'ADN en simple brin, et l'on fixe l'extrémité des simples brins 15 aléatoirement à la surface du « flowcell », sur laquelle on effectue une PCR « bridge » en phase solide (création de groupes (clusters) denses où les fragments sont amplifiés). Le séquençage s'effectue en ajoutant les 4 terminateurs réversibles marqués, les amorces et l'ADN polymérase, puis on lit la fluorescence émise par 20 chaque cluster, permettant de déterminer la première base. On effectue alors plusieurs cycles afin de lire l'ensemble de la séquence. Pour la mise en oeuvre de la technologie 454, on obtient une banque d'ADN simple brin matrice, des adaptateurs spécifiques étant ajoutés aux extrémités 3' et 25 5', et chaque brin d'ADN étant immobilisé sur une bille (un fragment d'ADN = une bille). Ces billes sont ensuite intégrées avec les produits d'amplification dans une émulsion eau-huile, afin de créer des « microréacteurs » (chaque goutte d'eau dans l'huile) contenant une seule bille. La PCR est réalisé dans cette émulsion l'ensemble de la banque étant amplifiée en parallèle, permettant d'obtenir plusieurs 30 millions de copies par bille. Puis on purifie les billes et charge les fragments sur des plaques telles que le diamètre des puits ne permet l'entrée que d'une seule bille à la fois. On ajoute les enzymes de séquençage et on envoie les nucléotides marqués individuels les uns après les autres. La détection de la séquence est effectuée par une caméra 35 CCD en fonction du signal luminescent. 9 Pour la technologie SOLiD, on prépare les banques, ajoute les adaptateurs et effectue une PCR dans une émulsion, comme dans la méthode 454. Puis on effectue un enrichissement des billes amplifiées, on modifie l'extrémité 3' des ADN pour permettre une fixation covalente sur une lame, et on dépose les billes sur la 5 lame. Le séquençage est effectué par ligation : des amorces s'hybrident sur les adaptateurs présents sur la matrice. Un jeu de 4 sondes de 2 bases marquées en fluorescence sont associées aux amorces. La spécificité des sondes de 2 bases s'effectue avec les lé' et 2"des bases de chaque réaction de ligation. Plusieurs cycles de ligation, détection et clivages sont effectués. Dans ce processus chaque 10 base est détectée par deux réactions de ligation indépendantes par deux différentes amorces. Le système de codage de la lecture sur deux bases permet une très grande fidélité de la lecture des résultats. Cette méthode permet de faire la différence entre les erreurs de séquençages et les variants réels (SNP, insertions et délétions). 15 Pour la technologie lonTorrent, on prépare des banques et on ajoute des adaptateurs. On effectue des PCR en émulsion. Le séquençage ne s'appuie pas sur la détection de fluorescence de nucléotides ou de leurs résidus de polymérisation par un capteur optique CCD, mais utilise un capteur CMOS qui 20 détecte les ions H+ dégagés lors de la polymérisation de l'ADN. Le capteur CMOS mesure le pH dans chacun des puits, ce qui indique la présence d'une ou plusieurs bases qui ont été intégrées dans l'ADN en cours d'analyse. On ajoute les bases les unes après les autres pour détecter laquelle est intégrée puis on rince et on recommence. 25 D'autres technologies de séquences existent telles que la technique MinION de Oxford Nanopore technologies (https://nanoporetech.com/products#minion, Mikheyev et Tin (2014). Molecular Ecology Resources. 14(6):1097-102.) ou Pac Bio de Pacific bioscience 30 (https://www.pacb.com/products-and-services/pacbio-systems/). Le procédé décrit ici permet de limiter le risque de détection d'un faux-positif (on conclue par erreur à la présence de l'allèle alternatif) ou d'un faux-négatif (on conclue par erreur à l'absence de l'allèle alternatif) que ces méthodes 35 de séquençage NGS peuvent présenter du fait du taux d'erreur de séquençage inhérent à chaque technologie. En effet, l'étape c) consiste en la détermination de 10 l'absence ou de la présence, pour un échantillon, d'une séquence non attendue dans les produits de séquençage. En cas de présence d'une telle séquence non attendue (correspondant à la présence d'un contaminant), il n'est nul besoin de quantifier la quantité de séquence non attendue par rapport à la quantité de 5 séquence attendue (correspondant à la séquence des graines correctes du lot de semences). La détection est donc uniquement qualitative (c'est-à-dire binaire : présence / absence d'une séquence d'un allèle alternatif à ou aux allèle(s) attendus. Le fait d'utiliser des sous-lots de graines permet également d'augmenter le nombre de graines étudiées pour chaque réaction de séquençage et ainsi 10 d'avoir un échantillon suffisant de graines tout en maîtrisant les coûts. La présence d'une telle séquence d'un allèle alternatif est indicatrice de la présence d'un contaminant pour cet allèle. Cette analyse est effectuée pour chaque région génomique analysée, c'est-à-dire pour chaque locus d'intérêt préalablement déterminé par l'homme du métier, 15 et permettant de caractériser le lot de semences. De fait, lorsque l'on choisit le nombre de graines dans chaque sous-lot de telle sorte qu'un seul contaminant soit présent (de façon statistique) au sein de ce sous-lot, la présence d'un allèle alternatif est suffisante pour conclure à la présence d'un unique contaminant. 20 L'étape suivante du procédé est le calcul du pourcentage effectif des contaminants dans le lot de semences. Ceci est effectué par la compilation des résultats qualitatifs obtenus pour l'ensemble des sous-lots. Le taux de pureté du lot de semence est alors estimé en considérant le 25 nombre de sous-lots contaminés, le nombre total de sous-lot analysés, et l'effectif de chacun des sous-lots. L'estimation de l'impureté du lot est obtenue suivant la formule : J=1 - (1 - - î dans laquelle n est le nombre de pools ; m est le nombre de grains dans un 30 pool d est le nombre de pools dans lesquels un contaminant a été identifié. On peut aussi déterminer l'intervalle de confiance de cette estimation par toute méthode statistique appropriée, notamment par une distribution F, tel qu'appliqué dans l'outil SeedCal utilisé dans le cadre de l'ISTA (International Seed Test Association) et tel qu'explicité dans Remund (2001). 11 Dans un mode de réalisation préféré, on effectue, à l'étape b), le séquençage ciblé de plusieurs régions du génome contenant plusieurs /oc/ 5 d'intérêt. Ceci permet de mieux garantir l'identité des semences présentes dans chaque échantillon et de détecter, de façon plus fine, la présence de contaminants. Ainsi, on peut séquencer de manière ciblée, au moins 2, de préférence, au moins 5, de préférence, au moins 10, de façon plus préférée au moins 100, 50, 40, 15 loci d'intérêt, voir au moins 20 loci d'intérêt. Même s'il n'existe pas de limite 10 supérieure au nombre de loci d'intérêt que l'on peut évaluer, on préfère limiter ceux-ci. En effet, il est possible de caractériser une variété avec un nombre de marqueurs (spécifiques de loci) limité (compris entre 15 et 20), et d'utiliser ce jeu de marqueurs pour discriminer les plantes de cette variété d'autres plantes. Une variété s'entend comme un ensemble de plantes ayant un même fond génétique, 15 la variété peut être une variété commercialisée, mais aussi une lignée non encore inscrite au catalogue, lignée de base, de pre-base ou lignée en cours de multiplication. Le nombre optimal de loci d'intérêt est défini par l'homme du métier, en fonction du matériel végétal considéré, mais également en fixant le nombre 20 minimal de loci discriminant toute paire de variétés donnée. Ainsi, le nombre minimal de loci discriminant toute paire de variétés peut être fixé à trois, limitant le risque de confondre une contamination réelle et un faux-positif expérimental. Différents algorithmes sont décrits par Rosenberg et al. (Journal of Computational Biology 12 (9), 2005, 1183-1201) pour sélectionner un jeu de marqueurs 25 discriminants. On peut améliorer ou modifier ces algorithmes pour prendre en compte d'autres critères tels que la qualité des marqueurs choisis (par qualité on entend leur aptitude à être amplifiés, identifiés sans équivoque). Des groupes ou catégories de marqueurs pourront être identifiées et définir un sous-groupes de 30 marqueurs qui contiendra préférentiellement des marqueurs issus d'un groupe donné ou issu de groupes différents. On peut ainsi définir un set de marqueurs que l'on souhaite utiliser. 12 L'algorithme peut aussi prendre en compte la qualité statistique de ces marqueurs définie comme le nombre minimum de marqueurs discriminants pour déclarer un couple d'individus comme différents. À partir de ce critère, la qualité de discrimination d'un set de marqueurs peut être évaluée par le nombre de couples 5 d'individus que ce set est capable de discriminer, idéalement la totalité des individus gérés par le producteur. Dans le contexte de la présente invention, on mettra préférentiellement en oeuvre la méthode sur des tact d'intérêt permettant à la fois de discriminer la variété d'intérêt (s'assurer de la constance et la concordance du fond génétique 10 entre les plantes) et d'identifier la présence ou l'absence d'autres /oc/ d'intérêt (notamment liés à des traits d'intérêt). Dans ce mode de réalisation, c'est-à-dire lorsque l'on effectue un séquençage de plusieurs régions du génome, on peut décider de considérer qu'il y a présence d'un contaminant dans un lot seulement si on observe la présence de 15 séquences non attendues pour plus d'un locus d'intérêt dans ce lot. En d'autres termes, on peut décider que, si l'on observe, dans un lot donné, la présence d'un unique allèle alternatif (une séquence non attendue pour une seule région du génome, alors que les séquences obtenues pour les autres régions sont celles attendues), on considère que la présence d'un contaminant n'est pas avérée. 20 La méthode décrite ici permet donc de déterminer la présence de contaminants dans un lot de semences, en particulier contrôler la pureté variétale lors d'un processus de production industriel. Cette méthode peut également être mise en oeuvre afin de vérifier le taux 25 de pureté d'un trait qui est recherché à l'état homozygote dans le lot de graines. Dans ce mode de réalisation, on évalue préférentiellement uniquement la région du génome contenant le trait particulier que l'on souhaite suivre. Plusieurs traits peuvent être suivis de façon simultanée, par utilisation de marqueurs spécifiques de chaque trait. 30 Par trait on entend forme allélique spécifique à un locus donné, dans ce contexte cette forme allélique peut être native, liée à une mutation identifiée par Tilling ou Ecotilling, mutation liée à l'empreinte d'un élément transposable, mutation obtenue par Gene Editing (édition de gène) ou par n'importe quelle autre méthode... dans ce contexte la mutation qu'elle soit une mutation ponctuelle, une 35 insertion ou une délétion implique un nombre limité de bases. Cette méthode pourra également s'appliquer sur un trait recherché à l'état hétérozygote, le 13 contaminant correspondra alors à une forme alternative aux formes alléliques attendues chez cet individu. Dans un mode de réalisation préféré, un trait (qui peut être lié à un seul allèle ou à plusieurs allèles) fournit à la plante un caractère phénotypique d'intérêt 5 (tel que résistance à la sécheresse, résistance au stress biotique, résistance au manque d'azote, augmentation du rendement...). Lorsque le trait est lié à une mutation impliquant une insertion de grande taille, telle qu'un trait OGM, un mutant obtenu par insertion d'un élément transposable ou un mutant obtenu par Gene Editing, la méthode pourra être mise 10 en oeuvre en recherchant la présence de la forme allélique ne contenant pas l'insertion ou la mutation considérée. La présence de cette forme allélique indiquant que la présence du trait lié à la mutation sous une forme homozygote dans le lot de graines n'est pas totalement garantie. Cette méthode pourra être utilisée par exemple lorsque la mutation correspond à l'introgression d'un fragment 15 d'ADN issu d'une autre espèce, ce cas particulier se rencontrera par exemple pour vérifier la pureté de lignées restauratrices de fertilité chez le colza. Cette méthode permet également de faire de la recherche de présence fortuite d'un trait, le trait dont on recherchera la présence fortuite pourra être un OGM, une mutation liée à du Gene Editing ou l'introgression d'un fragment 20 provenant d'une espèce hétérologue, cette recherche sera faite par amplification puis séquençage d'une région spécifique du T-DNA, ou de l'insertion. Par extension cette méthode peut s'appliquer à des traits liés à des mutation de faible taille si des amorces permettant s'amplifier spécifiquement la région lorsque l'on est en présence de la forme allélique mutée peuvent être définies. En adaptant le 25 protocole nombre de lots et nombre de graines par lot le protocole peut être étendu à l'identification de la présence de traits pour des fréquences allant par exemple jusqu'à 10% et dans ce contexte on pourra vérifier par exemple la présence de 10% de graines sauvages dans un lot de graines OGM (législation sur les zones refuge). Ces applications ne sont pas limitées aux OGM, le trait suivit par cette 30 méthode peut être l'introgression dans une lignée d'un fragment issue d'une autre espèce, la présence d'un locus restaurateur de fertilité issu du radis chez le Colza par exemple. De façon identique la vérification pourra permettre de vérifier que cette introgression est bien à l'état homozygote. Dans un autre mode de réalisation, la méthode peut être utilisée pour 35 détecter la présence fortuite (non désirée) d'OGM ou d'autre mutation liée à l'insertion d'un fragment de taille conséquente, dans un lot de semences. Cette 14 mutation peut être liée à la présence d'un élément transposable ou à une insertion obtenue notamment par Gene Editing. Dans ce mode de réalisation, on utilisera des amorces spécifiques d'un transgène ou de l'insertion particulier (si on suspecte une contamination particulière) ou différentes amorces génériques permettant de 5 détecter différents transgènes sans a priori. Dans le cas de la pureté variétale, on peut aussi ajouter des marqueurs liés à ces traits à la liste des marqueurs utilisés pour caractériser la variété. Ainsi, dans un mode de réalisation préféré, on effectue, les étapes b), c) et 10 d) pour plusieurs régions du génome contenant plusieurs loci d'intérêt. Dans ce mode de réalisation, on préfère quand un sous-ensemble de plusieurs loci permet de discriminer ou identifier une variété d'intérêt. Ainsi que vu plus haut, ce nombre de /oc/ est variable et ces loci peuvent être déterminés par l'homme du métier notamment selon les enseignements de Rosenberg (cité ci- 15 dessus). Dans un mode particulier de l'invention, il pourra intégrer des informations concernant le plan de production, impliquant des contrôles et des mesures particulières : distances d'isolement, zones de bordures, castration, ce qui implique que le risque de contamination va être limité et le lot de semences sera a priori non contaminé ou faiblement contaminé. Par ailleurs de part ces mesures, une 20 contamination sera très probablement issue d'un contaminant connu, notamment d'une lignée parentale, y compris les lignées parentales impliquées dans la production des semences de bases et pre base. Dans ce contexte particulier le nombre de marqueur permettant d'identifier la pureté d'une lignée peut être très réduit, il pourra notamment être de 20 ou moins. 25 Ainsi que vu plus haut, dans un mode de réalisation, on déclare un lot comme contenant un contaminant si l'on observe un allèle alternatif à l'allèle attendu pour un seul locus d'intérêt. Dans un autre mode de réalisation, on déclare un lot comme contenant un contaminant si l'on observe un allèle alternatif à l'allèle attendu pour plus d'un locus d'intérêt (notamment 2 ou 3 loci). 30 Dans un mode de réalisation, au moins ou exactement un locus d'intérêt est lié à un caractère d'intérêt (trait). Dans un autre mode de réalisation, c'est une combinaison de loci qui est est liée à un caractère d'intérêt (trait). Dans un mode de réalisation, au moins un locus d'intérêt est lié à un trait spécifique a priori non présent dans les graines du lot. Dans ce mode de 35 réalisation, on recherche la présence fortuite de ce trait. On ajoute donc des marqueurs pour vérifier l'absence du trait. Dans ce mode de réalisation, la 15 méthode est essentiellement qualitative. L'intégration de ces marqueurs dans le protocole revendiqué permet de faire dans une seule expérience des contrôles supplémentaires nécessaires par ailleurs. Dans un mode de réalisation préféré, la quantité de graines dans chaque 5 sous-lot préparé à l'étape a) est comprise entre 80 et 120. Dans un mode de réalisation préféré, on peut identifier le(s) contaminant(s) présent(s) dans le lot de semence. Pour chaque sous-échantillon, il est possible de définir un profil moléculaire 10 correspondant à la compilation des données de chaque locus d'intérêt. Le profil de chaque sous-échantillon peut alors être comparé au profil moléculaire attendu, et un profil moléculaire contaminant peut être déduit par soustraction. Ainsi, un locus d'intérêt ne présentant pas d'allèle alternatif sera considéré identique au locus entre la variété attendue et le contaminant, tandis qu'un locus présentant un allèle 15 alternatif sera défini comme potentiellement homozygote pour l'allèle alternatif, ou hétérozygote allèle attendu/allèle alternatif. Ces profils moléculaires contaminants peuvent alors être comparés à une base de données de référence afin d'identifier la nature du contaminant, et éventuellement le moment où il est entré dans le cycle de production. 20 Ainsi, on envisage un procédé d'identification du contaminant, qui met en oeuvre la méthode telle que décrite ci-dessus, et qui comprend en outre les étapes consistant à i) définir le profil moléculaire du contaminant de chaque sous-lot 25 contaminé par comparaison du profil observé dans ce sous-lot avec le profil attendu en l'absence de contaminant, et ii) comparer le profil obtenu en i) avec ceux d'une base de données de référence. 30 Alternativement, on considère un procédé de détermination du degré de pureté, tel que défini ci-dessus, caractérisé en ce que l'on procède en outre à l'identification du contaminant pour chaque sous-lot contaminé en i) déduisant le profil moléculaire du contaminant dans un sous-lot contaminé par comparaison du profil observé dans ce sous-lot avec le profil 35 attendu en l'absence de contaminant et en 16 ii) Comparant le profil obtenu en i) avec ceux d'une base de données de référence. On obtient donc un ou plusieurs profils de contaminants pour le lot de semences de départ, correspondant à la somme des contaminants de chaque 5 sous-lot contaminé. DESCRIPTION DES FIGURES Figure 1 : résultat de l'analyse en Taqman pour un SNP, comprenant deux formes alléliques détectées respectivement par les fluorochromes FAM et VIC, dans des 10 échantillons de mais homozygotes (A, B) ou hétérozygote pour le SNP (C). A : échantillon homozygote pour la forme allélique détecté en FAM. B : échantillon homozygote pour la forme allélique détecté en VIC. C : échantillon hétérozygote pour les formes alléliques détectés en FAM et VIC. Figure 2 : Fréquence relative, dans chaque sous-lot, de l'allèle alternatif pour le 15 SNP10. Les sous-lots 3, 14 et 16 montrent une fréquence de l'allèle alternatif significative. Figure 3: Profil qualitatif (présence/absence d'un allèle contaminant) Profil de présence d'un allèle alternatif pour les 17 marqueurs (ligne) (16 marqueurs discriminants et un marqueur associé à un trait) au sein des 16 sous-lots (colonne). 20 La présence d'un allèle alternatif est détectée pour au moins 3 SN P dans les sous-lots 3, 14 et 16. Ces sous-lots sont déclarés contaminés. Les 13 autres sous-lots sont déclarés non contaminés. Figure 4 : profils moléculaires obtenus sur les 17 SNP (16 marqueurs discriminants et un marqueur associé à un trait) obtenus sur les 16 sous-lots analysés. Le profil 25 de la première ligne correspond au profil majoritaire, les profils suivants aux profils contaminés observés pour les lots 3, 14 et 16 respectivement. EXEMPLES Exemple 1 : Détection de contaminants par Taqman 30 Cet exemple évalue la possibilité de détecter une graine contaminante dans un sous-lots de graines de maïs, par génotypage utilisant la technologie Taqman (Applied Biosystem). La figure 1 montre le résultat de l'analyse en Taqman pour un SNP, comprenant deux formes alléliques détectées respectivement par les fluorochromes FAM et 35 VIC, dans des échantillons de maïs homozygotes ou hétérozygote au SNP, et met en évidence la présence de signal avec la sonde FAM dans un échantillon 17 homozygote pour l'allèle VIC (B), c'est-à-dire un signal non spécifique, ne permettant pas de distinguer un signal faux positif d'un signal lié à une contamination réelle dans un échantillon. Ces résultats montrent que la méthode Taqman ne permet pas de détecter des 5 contaminants de façon fiable. Exemple 2 : Détection de contaminants par génotypage sur puce Dans cet exemple des lots de 200 graines issues d'une lignée A contenant 10%, 20%, 30%, 40%, et jusqu'à 90% de contaminants par une lignée B ont été 10 élaborés et un échantillon de 15 graines issues de ce lot ont été analysées par génotypage sur une puce Infinium (Illumina), afin d'évaluer la faisabilité de l'identification d'une contamination. On arrive à détecter les contaminations supérieures à 10 %, mais les mélanges contenant 10% de contamination ne sont pas distinguables des témoins non contaminés. A fortiori, les contaminations moins 15 importantes ne seront pas détectables. Exemple 3: mise en oeuvre de la méthode selon l'invention sur un set de marqueurs Dans cet exemple, un set de 16 marqueurs discriminants (SNP) a été 20 utilisé, permettant d'identifier sans ambiguïté la présence d'une variété autre que celle attendue. Ce set de 16 marqueurs a été défini à partir de données de génotypage de référence sur plusieurs milliers de marqueurs pour les variétés d'intérêt, et permet de différencier chacune d'elle par rapport aux autres grâce à au moins 3 marqueurs discriminants. Dans ce cas c'est le profil moléculaire global sur 25 les 16 marqueurs qui détermine l'identité de chaque variété. Chaque marqueur est spécifique d'un locus d'intérêt. Dans une expérience en condition de contamination contrôlée, 24 graines d'une lignée pure Li ont été introduites dans un lot de 2376 graines d'une lignée pure L2, le lot ainsi obtenu a un taux de pureté de 99%, les graines ont été 30 distribuées aléatoirement en vingt-quatre sous-lots de 100 grains (soit 2400 grains analysés) chaque lot de graines ainsi obtenu est broyé indépendamment et l'ADN est extrait des broyats. Ainsi, on a en moyenne 1 contaminant par lot : le nombre de sous-lots est en effet égal au nombre de contaminants présents dans le lot complet de semences. Du fait de la distribution statistique aléatoire, on sait 35 toutefois que certains sous-lots ne contiendront pas de contaminants, et que 18 d'autres sous-lots contiendront plusieurs contaminants, du fait de l'échantillonnage en formant les sous-lots Pour chacun des 16 marqueurs, un amplicon de 70 à 120bp a été défini, et les 16 marqueurs co-amplifiés par PCR multiplex. Un index unique (TAG) est 5 utilisé pour chaque échantillon d'ADN, permettant ainsi de faire un séquençage de l'ensemble des amplicons et d'attribuer les séquences obtenues à leur lot d'origine. Les amplicons ont été séquencés par la technologie. Illumina sur un séquenceur Miniseq. Des séquences appariées de 75 bases ont été générées, attribuées aux ADN d'origine par une étape de démultiplexage. Après suppression 10 des séquences des adaptateurs et des bases de mauvaise qualité (seuil Q30), chaque paire de séquences est rassemblée en une séquence unique, puis alignée sur le génome de référence du maïs (RefGenV4). Pour chaque SNP, les fréquences alléliques relatives de l'allèle majoritaire et de l'allèle alternatif ont été calculées, et correspondent au nombre de lectures contenant l'allèle d'intérêt par 15 rapport à la somme des lectures de chaque allèle. On considère qu'il y a contamination pour un marqueur SNP si, dans un sous-lot, la séquence d'une forme allélique, qui n'est pas celle de l'allèle attendu pour la variété testée, apparaît de façon supérieure au bruit de fond. Un échantillon est déclaré contaminé dès lors qu'il contient au moins 3 SNP 20 pour lesquels un allèle alternatif est détecté. Ainsi, on conclut que, parmi ces 24 sous-lots, 13 sont considérés comme contaminés et 11 comme purs. Le nombre de sous-lots contaminés permet d'estimer la pureté variétale du lot analysé, ce calcul est réalisé à l'aide du logiciel Seed Calc qui utilise les formules de Remund (2001). Dans cet exemple, la pureté estimée est de 99.22% 25 (98.64%-99.6%), pour une pureté réelle contrôlée de 99%. L'estimation de l'impureté fi du lot est obtenue suivant la formule : d dans laquelle n est le nombre de pools ; m est le nombre de grains dans un pool ; d est le nombre de pools dans lesquels un contaminant a été identifié. 30 Dans le cas ci-dessus : 141-13/24)001= 1-0.9922 = 0.0078 soit une pureté de 99,22. L'intervalle de confiance est également calculé selon les procédés décrits dans Remund 2001. 19 Exemple 4 : identification du contaminant Dans cet exemple, des lots de semences de base de maïs ont été analysés selon la même approche que celle citée dans l'exemple 3. Pour un lot, 16 sous-lots de 100 graines ont été constitués. 5 Les graines de chaque sous-lot ont été broyées et l'ADN extrait. Un set de 17 marqueurs, dont 16 SNP discriminants (permettant d'identifier sans ambiguïté la présence d'une variété autre que celle attendue) et un marqueur associé à un trait, a été identifié. Pour chaque marqueur, un amplicon de 70-120bp a été défini, et les 17 marqueurs ont été co-amplifiés par PCR multiplex. Un index (Tag) unique est 10 utilisé pour chaque échantillon d'ADN, permettant ainsi de faire un séquençage de l'ensemble des amplicons et d'attribuer les séquences obtenues à leur lot d'origine. Les amplicons ont été séquencés par la technologie Illumina sur un séquenceur Miniseq. Des séquences appariées de 75 bases ont été générées, attribuées aux ADN d'origine par une étape de démultiplexage. Après suppression 15 des séquences des adaptateurs et des bases de mauvaise qualité (seuil Q30), chaque paire de séquences est rassemblée en une séquence unique, puis alignée sur le génome de référence du maïs (RefGenV4). Pour chaque SNP, les fréquences alléliques relatives de l'allèle majoritaire et de l'allèle alternatif ont été calculées, et correspondent au nombre de lectures contenant l'allèle d'intérêt par 20 rapport à la somme des lectures de chaque allèle. La figure 2 montre, pour un SNP (SNP10), la fréquence de l'allèle alternatif dans chacun des sous-lots (c'est-à-dire la fréquence d'apparition de la séquence de l'allèle alternatif). Dans cet exemple, les sous-lots 3, 14 et 16 montrent une présence significative de l'allèle alternatif (au-dessus du bruit de fond représenté 25 par la ligne horizontale). Cette analyse est réalisée pour chaque SNP, et la figure 3 montre le profil qualitatif (présence/absence de l'allèle alternatif) obtenu pour chaque SNP dans chaque sous-lot. On confirme la présence d'un allèle alternatif pour au moins 3 SNP dans les sous-lots 3, 14 et 16. Ces 3 sous-lots sont déclarés contaminés. Les 13 autres sous-lots sont déclarés non contaminés. Le taux de 30 pureté variétal estimé avec SeedCalc est de 99.79% (intervalle de confiance à 95% : 99.39% - 99.96%). En parallèle, le même lot a été analysé sur 558 graines individuelles. Pour chaque graine, un fragment est prélevé en poinçonnant l'embryon à l'aide d'un emporte-pièce, puis l'ADN extrait et le génotypage réalisé avec la technologie 35 KASP (LGC Genomics) sur 16 marqueurs discriminants. Cette analyse permet 20 d'estimer une pureté de 99.46% (intervalle de confiance à 95% 98.42% 99.89%). Le marqueur SNP17 a été analysé séparément et permet d'estimer la pureté du trait associé. 5 La figure 3 montre que les sous-lots 3 et 16 présentent une fréquence significative de l'allèle alternatif. Ces 2 sous-lots sont déclarés contaminés, conduisant à une estimation de la pureté de trait de 99.87% (intervalle de confiance à 95% : 99.52 99.98%). Le profil moléculaire identifié sur les sous-lots non contaminés est dans un premier 10 temps utilisé pour vérifier sa conformité par rapport au profil attendu pour la variété analysée (l'étape précédente permet de vérifier la pureté variétale du lot, cette étape permet de vérifier que la variété identifiée est bien celle attendue). Ensuite, sur les sous-lots 3, 14 et 16 montrant une contamination, un profil moléculaire contaminant est déduit du profil moléculaire observé, par soustraction du profil 15 attendu. Pour chaque marqueur SNP montrant une contamination, les 2 allèles observés sont reportés (Figure 4). Le contaminant peut ainsi être homozygote pour l'allèle minoritaire, ou bien hétérozygote. Chaque profil moléculaire contaminant est ensuite comparé à une base de 20 données de référence afin de l'identifier. Si ce génotype correspond à une accession connue, celle-ci est proposée comme contaminant potentiel, sinon le génotype contaminant est déclaré non identifiable. Cette base de données de référence peut être affinée en fonction du plan de production notamment, cette base contiendra alors prioritairement l'ensemble REVENDICATIONS1. Procédé de détermination de la quantité de contaminants à au moins un locus d'intérêt, présents dans un lot de semences d'une variété d'intérêt, caractérisé en ce a) on regroupe des graines d'un lot de semences par sous-lots d'au moins 10 graines, le nombre sous-lots ainsi obtenus étant supérieur ou égal à 10 b) on effectue, pour chaque sous-lot, un séquençage ciblé d'au moins la région du génome des graines, contenant le locus d'intérêt, c) on détermine, pour chaque sous-lot, la présence d'un contaminant de façon qualitative, en cas de détection d'un allèle alternatif à ou aux l'allèle(s) attendus pour chaque région génomique séquencée (présence / absence de / des l'allèle(s) attendus) d) on détermine la quantité de contaminants dans le lot global par la compilation des résultats qualitatifs obtenus pour l'ensemble des sous-lots. 2. Procédé selon la revendication 1, caractérisé en ce que l'on effectue, les étapes b), c) et d) pour plusieurs régions du génome correspondant à plusieurs tact d'intérêt. 3. Procédé selon la revendication 2, caractérisé en ce qu'un sous ensemble de ces loci d'intérêt est suffisant pour permettre d'identifier la variété d'intérêt. 4. Procédé selon la revendication 3, caractérisé en ce que l'on déclare un lot comme contenant un contaminant si l'on observe un allèle alternatif à ou aux allèle(s) attendus pour un seul locus d'intérêt. 5. Procédé selon la revendication 3, caractérisé en ce que l'on déclare un lot comme contenant un contaminant si l'on observe un allèle alternatif à ou aux allèle(s) attendus pour plus d'un locus d'intérêt. 6. Procédé selon l'une des revendications 1 à 5, caractérisé en ce qu'au moins un locus d'intérêt est lié à un caractère d'intérêt (trait). 7. Procédé selon la revendication 2, caractérisé en ce qu'une combinaison de tact est liée à des caractères d'intérêt (trait). 22 8. Procédé selon la revendication 2, caractérisé en ce qu'une combinaison de loci est liée à un caractère d'intérêt (trait). 5 9. Procédé selon l'une des revendications 1 à 8, caractérisé en que au moins un locus d'intérêt est lié à un trait spécifique a priori non présent dans les graines du lot, afin de détecter la présence fortuite de ce trait. 10. Procédé selon l'une des revendications 1 à 9, caractérisé en ce que la quantité 10 de graines dans chaque sous-lot préparé à l'étape a) est comprise entre 80 et 120. 11 Procédé selon l'une des revendications 1 à 10, caractérisé en ce que l'on procède en outre à l'identification du contaminant pour chaque sous-lot 15 contaminé en i) déduisant le profil moléculaire du contaminant dans un sous-lot contaminé par comparaison du profil observé dans ce sous-lot avec le profil attendu en l'absence de contaminant et en ii) Comparant le profil obtenu en i) avec ceux d'une base de données 20 de référence.