La théorie du séquençage de l'ADN désigne l'ensemble des travaux qui tentent de jeter les bases analytiques pour déterminer l'ordre des nucléotides spécifiques dans une séquence d'ADN donnée, autrement connu sous le nom de séquençage de l'ADN. Les aspects pratiques de cette théorie tournent autour de la conception et de l'optimisation de projets de séquençage (connus sous le nom de "génomique stratégique"), de la prévision des performances du projet, de la résolution des problèmes rencontrés dans les résultats expérimentaux, de la caractérisation de facteurs tels que le biais séquentiel et les effets des algorithmes des logiciels de traitement, et la comparaison de diverses méthodes de séquençage entre elles. En ce sens, elle pourrait être considérée comme une branche des domaines de l'ingénierie des systèmes ou de la recherche opérationnelle. Les archives permanentes de ces travaux sont principalement mathématiques, bien que des calculs numériques soient également souvent effectués pour traiter des problèmes particuliers. La théorie du séquençage de l'ADN traite des processus physiques liés au séquençage de l'ADN et ne doit pas être confondue avec les théories de l'analyse des séquences d'ADN résultantes, comme par exemple l'alignement des séquences. Les publications scientifiques. ne font parfois pas de distinction minutieuse entre les deux, mais la deuxième catégorie est principalement concernée par les problèmes algorithmiques. La théorie du séquençage est basée sur des éléments issus des mathématiques, de la biologie et de l'ingénierie des systèmes, ce qui en fait donc un domaine hautement interdisciplinaire. Le sujet peut être étudié dans le cadre de la biologie numérique .
Théorie et stratégies de séquençage
Le séquençage comme problème de couverture Toutes les méthodes traditionnelles de séquençage d'ADN reposent sur la lecture de petits fragments d'ADN puis la reconstruction ultérieure de ces données pour déduire la cible d'ADN d'origine, soit par assemblage, soit par alignement sur une séquence de référence. L'abstraction commune à ces méthodes est celle d'un problème de couverture mathématique.. Par exemple, on peut imaginer un segment de ligne représentant la cible et un processus ultérieur où des segments plus petits sont "déposés" sur des emplacements aléatoires de la cible. La cible est considérée comme "séquencée" lorsqu'une couverture (covering en anglais) adéquate y est appliquée (par exemple, lorsqu'il n'y a plus de lacunes ou gaps en anglais).
Les propriétés abstraites de cette couverture sont étudiées par les mathématiciens depuis plus d'un siècle.. Cependant, l'application directe de ces résultats n'est généralement pas possible. Les solutions mathématiques de forme fermée, en particulier pour les distributions de probabilité, ne peuvent souvent pas être facilement évaluées, c'est-à-dire qu'elles impliquent des quantités excessivement longues de temps de processing informatique pour les paramètres caractéristiques du séquençage d'ADN. La configuration de Stevens en est un exemple.. Les résultats obtenus du point de vue des mathématiques pures ne tiennent pas compte non plus des facteurs qui sont réellement importants dans le séquençage, comme par exemple le chevauchement détectable dans les fragments de séquençage, le double brin, les effets de bord et la multiplicité des cibles. Par conséquent, le développement de la théorie du séquençage a progressé davantage sous l'angle de la philosophie des mathématiques appliquées. En particulier, il a été axé sur les problèmes rencontrés et fait un usage judicieux des approximations, des simulations, etc.
Premières utilisations dérivées de la théorie des probabilités élémentaires Le premier résultat peut être trouvé directement à partir de la théorie des probabilités élémentaires. Supposons que nous modélisons le processus ci-dessus en prenant et comme longueur de fragment et longueur cible, respectivement. La probabilité de "recouvrir" un emplacement donné sur la cible avec un fragment particulier est alors . (Cela suppose , qui est souvent valable, mais pas pour tous les cas issus du monde réel.) La probabilité qu'un seul fragment ne couvre pas un emplacement donné sur la cible est donc , et pour fragments. La probabilité de couvrir un emplacement donné sur la cible avec au moins un fragment est donc
Cette équation a d'abord été utilisée pour caractériser les banques de plasmides. mais elle peut apparaître sous une forme modifiée. Dans la plupart des projets de séquençage on a , de sorte que, à un bon degré d'approximation, on obtient
où est appelé la redondance. À noter l'importance de la redondance comme représentant le nombre moyen de fois qu'une position est couverte par des fragments. À noter également qu'en considérant le processus de couverture sur toutes les positions dans la cible, cette probabilité est identique à la valeur attendue de la variable aléatoire , qui est la fraction de la couverture cible. Le résultat final,
reste largement utilisé comme estimateur de "l'arrière de l'enveloppe" et prédit que la couverture pour tous les projets évoluera le long d'une courbe universelle qui n'est fonction que de la redondance.
Théorie de Lander-Waterman En 1988, Eric Lander et Michael Waterman ont publié un article important. examinant le problème de couverture du point de vue des lacunes (gaps). Bien qu'ils se soient concentrés sur le dénommé problème de cartographie génétique, l'abstraction du séquençage est à peu près la même. Ils ont fourni un certain nombre de résultats utiles qui ont été adoptés comme théorie standard dès les premiers pas du séquençage du génome à "grande échelle".. Leur modèle a également été utilisé dans la conception du projet du génome humain et continue de jouer un rôle important dans le séquençage de l'ADN.
En fin de compte, l'objectif principal d'un projet de séquençage est de combler toutes les lacunes, de sorte que la "perspective des lacunes" représentait une base logique pour développer un modèle de séquençage. Un des résultats les plus fréquemment utilisés de ce modèle est le nombre attendu de contigs, étant donné le nombre de fragments séquencés. Si l'on néglige la quantité de séquences qui est essentiellement "gaspillée" au cours de la détection des chevauchements, leur théorie donne
En 1995, Roach publié des améliorations à cette théorie, lui permettant d'être appliquée à des projets de séquençage dont le but était de séquencer complètement un génome cible. Michael Wendl et Bob Waterston. confirmé, sur la base de la méthode de Stevens, que les deux modèles produisaient des résultats similaires lorsque le nombre de contigs était important, comme dans les projets de cartographie à faible couverture ou de séquençage. À mesure que les projets de séquençage montaient en puissance dans les années 1990 et que les projets étaient en passe d'être achevés, les approximations de faible couverture devenaient inadéquates et l'utilisation du modèle exact de Roach était devenue nécessaire. Cependant, comme le coût du séquençage a chuté, les paramètres des projets de séquençage sont devenus plus faciles à tester directement de manière empirique, et l'intérêt de même que le financement pour la génomique stratégique ont diminué.
Les idées de base de la théorie de Lander-Waterman ont conduit à l'obtention d'un certain nombre de résultats supplémentaires concernant des variations particulières des techniques de cartographie... Cependant, les progrès technologiques ont rendu les théories de la cartographie largement obsolètes, sauf dans les organismes autres que les organismes modèles très étudiés (par exemple, les levures, les mouches, les souris et les humains).
Stratégie de stationnement La stratégie de stationnement dans le cadre du séquençage ressemble au processus de stationnement des voitures le long d'un trottoir. Chaque voiture correspond à un clone séquencé et le trottoir correspond à la cible génomique.. Chaque clone séquencé est criblé pour garantir que les clones séquencés ultérieurement ne chevauchent aucun clone précédemment séquencé. Aucun effort de séquençage n'est redondant dans cette stratégie. Cependant, tout comme les espaces entre les voitures garées, des gaps non séquencés inférieurs à la longueur d'un clone s'accumulent entre les clones séquencés, et il peut s'avérer très coûteux de les combler.
Séquençage d'extrémités par paires En 1995, Roach et al.. ont proposé et démontré, à travers des simulations, une généralisation d'un ensemble de stratégies explorées précédemment par Edwards et Caskey.. Cette méthode de séquençage du génome entier est devenue extrêmement populaire car elle a été défendue par Celera et utilisée pour séquencer plusieurs organismes modèles avant que Celera ne l'applique au génome humain. Aujourd'hui, la plupart des projets de séquençage utilisent cette stratégie, souvent appelée séquençage d'extrémité par paires.
Avancées depuis le projet sur le génome humain Les processus physiques et les protocoles de séquençage de l'ADN ont continué d'évoluer, en grande partie grâce aux progrès des méthodes biochimiques, de l'instrumentation et des techniques d'automatisation. Le séquençage de l'ADN a permis de résoudre un grand nombre de problèmes, y compris la métagénomique et le séquençage médical (cancer). Il y a des facteurs importants dans ces scénarios que la théorie classique ne prend pas en compte. Des travaux récents ont commencé à se concentrer sur la résolution des effets causés par certains de ces problèmes. Le niveau des mathématiques devient proportionnellement plus sophistiqué.
Divers artefacts de séquençage à grand insert Les biologistes ont mis au point des méthodes pour filtrer des régions de génomes hautement répétitives et essentiellement non séquençables. Ces procédures sont importantes lorsqu'on travaille sur des organismes dont le génome se compose principalement d'un tel ADN, par exemple le maïs. Ils produisent une multitude de petits "îlots" de produits d'ADN séquençables. Wendl et Barbazuk. ont proposé une extension de la théorie de Lander-Waterman pour tenir compte des gaps dans la cible dues au filtrage et à ce que l'on appelle "l'effet de bord" (edge-effect en anglais). Ce dernier est un biais d'échantillonnage spécifique à une position ; par exemple, la position de la base terminale (dernière base d'un fragment séquencé) n'a que chance d'être couverte, alors que pour les positions plus intérieures à la séquence, la chance qu'elles soient couvertes est de . Pour , la théorie classique de Lander-Waterman donne toujours de bonnes prévisions, mais la dynamique change pour des redondances plus élevées.
Les méthodes de séquençage modernes séquencent généralement les deux extrémités d'un plus grand fragment, ce qui fournit des informations de liaison pour l'assemblage de novo et de meilleures probabilités d'alignement sur la séquence de référence. Les chercheurs pensent généralement que des longueurs de données (longueurs de lecture) plus longues améliorent les performances lorsqu'on séquence de très grandes cibles d'ADN, une idée cohérente avec les prévisions des modèles de distribution.. Cependant, Wendl. a montré que des fragments plus petits offrent une meilleure couverture sur de petites cibles linéaires car ils réduisent l'effet de bord dans les molécules linéaires. Ces résultats ont des implications pour le séquençage des produits des procédures de filtrage de l'ADN. La lecture par paires (ou non) et la taille des fragments ont évidemment une influence négligeable lorsqu'on séquence des grandes cibles comme celles de la classe du génome entier.
Séquençage individuel et de population