$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
La disponibilité de séquençage rapide, fiable et rentable en profondeur a révolutionné beaucoup d’aspects dans le domaine des sciences de la vie, permettant à grande profondeur dans les analyses axées sur le séquençage. Un défi reste réside dans la conception innovante et la création du représentant les bibliothèques de séquençage. Nous décrivons ici un protocole visant à capturer les molécules de cDNA viral naissant, spécifiquement les intermédiaires du processus de transcriptase inverse du VIH-1.
L’étape la plus critique de cette stratégie est la ligature d’un adaptateur pour l’ouvrir 3'-termini de manière quantitative et non biaisée. Efficacité des trompes entre deux terminus d’ADN simple brin, inter- et intramoléculaires, ont été étudiés et optimisés pour différentes applications11,26,27,28,29. Le choix d’utiliser un adaptateur en épingle à cheveux avec T4 DNA ligase dans les conditions décrites à l’étape 3.3 est le résultat de l’optimisation empirique dont nous avons évalué différentes ligases, adaptateurs et réactifs pour la ligature des oligonucléotides synthétiques représentant Séquences de VIH-1 (tableau 2) (données non présentées). Dans ces réactions de test in vitro , nous avons confirmé que la T4 DNA ligase médiée par ligation de l’adaptateur de l’épingle à cheveux, tel que décrit par Kwok et al. 11a un écart très faible et réalise près de ligature complète de molécules accepteur lorsque l’adaptateur est utilisé en excès. L’efficacité de la ligature n’est pas affectée par l’ajout de séquence nucléotidique pour restituer l’adaptateur compatible pour le système d’amorçage multiplex (voir Figure 4). En comparaison, nous avons constaté qu’une ligase ADN/ARN thermostable de 5' (voir « Ligase A », Table des matières pour les ligases exactes par rapport ici), qui est une ligase machiné de RNA qui a été développé en partie d’améliorer efficacité ligature avec ADN simple brin comme accepteur 27, était en effet plus efficace à la ligature des deux molécules d’ADN simple brin de ligase RNA (« Ligase B »), mais avait un biais important, avec fortes différences dans l’efficacité de la ligature, même entre les oligonucléotides avec des différences de longueur de base individuelle [tableau 2 ; HTP con mi G (a) et (b)]. En outre, nous avons trouvé seulement un biais minime dans les réactions avec « Ligase C » combiné avec un adaptateur transportant un essai randomisé 5'-termini (une stratégie pour compenser les biais de nucléotides connue de « Ligase C » ; voir par exemple Ding et al. 30). Cependant, le « Ligase C »-médiation des trompes intermoléculaires étaient incomplètes, rendant le système de ligase d’ADN T4 le meilleur choix.
Plusieurs mesures de contrôle de la qualité sur le golf du protocole et de l’inclusion des témoins positifs et négatifs permettant de détecter des problèmes potentiels avant la poursuite de l’analyse et des directives pour le dépannage des efforts. Les quantifications de qPCR en étapes 2.2.2 et 2.3.12 veiller à ce que la quantité du matériel d’entrée suffit. Nombre de copies d’ADNc typiques dans la gamme d’élution (à l’étape 2.1) 200 µL d’environ 10 000 à 300 000 par µL. L’étape de capture hybride peut entraîner une perte globale du VIH-1 cDNA quantité mais devrait se traduire par un fort enrichissement spécifiques applicables à l’ARNC VIH-1 au cours de l’ADN cellulaire, qui peut être déterminé en utilisant des amorces appropriées pour quantifier l’ADN génomique, avant et après enrichissement par qPCR ou en mesurant la concentration totale de l’ADN. Récupéré le VIH-1 cDNA après l’hybride capture étapes devrait être au moins 10 % de l’apport. Faible à partir de matériel peut expliquer autrement un témoin positif de l’oligonucléotide réussie (voir étape 3.3.2) mais seulement limité lectures atteints dans les échantillons. Faible lire nombres globale pourrait également s’expliquer par une surestimation de la concentration de la bibliothèque en raison de la présence d’espèces d’ADN non pertinents sans MiSeq adaptateurs. Cela se traduirait par une densité faible de cluster et peut être améliorée par la détermination de la concentration des séquences de VIH-1 dans la bibliothèque de qPCR en plus de la quantité d’ADN totale par dosages fluorimétrique. En raison de la nature hautement délicate de la méthode, spécial il faut pour éviter la contamination même à basse altitude, les deux autres échantillons (en particulier, depuis les stocks oligonucléotide contrôle de concentration élevée) ainsi que de matériel de laboratoire. Travaillant dans une UV de stérilisation workstation PCR est bénéfique à cet égard. L’électrophorèse automatisée de la bibliothèque finale (étape 6.1.2) est une autre mesure de contrôle de la qualité. Le calibre de l’acide nucléique généralement observé est entre 150 et 500 nt. amorces qui peuvent être détectées dans le contrôle facultatif après le PCR et avant purification (voir la note à l’étape 5.2) convient désormais absente. Un résultat représentatif, la courbe d’intensité échantillon a un pic autour de 160 à 170 nt et une nette deuxième pic autour de 320 à 350 nt. Cela reflète probablement l’abondance plu souvent vus en relativement courte transcription inverse (longueur d’insert nt 1 à 20) et forte non-stop pleine longueur (longueur de 180 à 182 nt insert) (Figure 3b).
Alors que le protocole présenté et les amorces sélectionnées sont spécifiques pour des constructions de transcriptase inverse du VIH-1 au début, la méthode est généralement applicable à toute étude visant à déterminer ouvert 3'-termini d’ADN. Les principales modifications requises dans d’autres contextes sera la méthode de capture hybride et la stratégie de conception d’amorce. Par exemple, si la cible doit être adapté aux fin transcrits du VIH-1, un plus grand nombre de différents oligonucléotides biotinylé capture recuit sur toute la longueur de l’ADNc serait souhaitable et diminuera probablement la perte lors de l’étape de capture hybride. Tel que mentionné dans l’introduction, il est important de tenir compte des limitations lors de la conception de la gamme sur laquelle est les 3'-termini à détecter pour éviter les différentes sources de biais. Tout d’abord, il peut y avoir un biais dans les réactions de PCR si les modèles avec l’adaptateur sont de longueur très variable. Deuxièmement, la plate-forme de séquençage utilisée ici (p. ex., MiSeq) a une gamme de longueur d’insertion privilégiée pour le clustering optimale et significativement plus courtes et plus de produits ne peuvent pas être séquencés avec la même efficacité. En partie, cela peut être résolu par le calcul, comme l’a fait en calculant un facteur de correction du biais de longueur linéaire (voir la Figure 4, graphique du bas). Toutefois, si la région où 3'-termini cartographie est souhaitée est longue (> 1000 nt), il est plus judicieux de diviser les réactions avec les transcriptions ligaturées et utilisent plusieurs amorces en amont pour évaluer les extrémités 3' dans les sections.
Le programme d’analyse a été écrite à l’interne dans le but spécifique de l’analyse de ces deux derniers nucléotides de la séquence de HIV-1 adjacente à la séquence de l’adaptateur fixe ainsi que la variation de base de toutes les bases de détecter des mutations. Les différentes étapes, répartis comme suit : tout d’abord, les séquences de l’adaptateur sont taillés à l’aide de l’outil fastx-0.0.13 ; Ensuite, toutes les séquences qui sont en double (c'est-à-dire des séquences identiques, y compris le code à barres) sont supprimés. Toutes les lectures uniques restants sont alors alignés sur la séquence de VIH-1 en utilisant le nœud papillon (http://bowtie-bio.sourceforge.net/index.shtml) avec le décalage maximal fixé à trois bases. La séquence de modèle est composée de la première 635 nt du VIH-1 cDNA (souche NL4.3), qui comprend la séquence - sss et le premier produit de transfert de brin jusqu'à la piste d’appariement (U5-R-U3-PPT ; voir Figure 1). Ainsi, le logiciel fourni et modèles conviennent uniquement directement si la méthode est utilisée pour la même application (détection de début transcription inverse de l' HIV-1NL4.3). Ajustements devront être faites pour les autres séquences cibles. Les positions de la 3'-termini pour chaque lecture ont été déterminées par la position dans l’alignement. Les appels de base pour chaque position sont enregistrés et des taux de mutation sont calculés à partir de la couverture totale de chaque base, qui varie, comme lectures sont de différentes longueurs et inserts de longs ne peuvent pas être entièrement couverts par le séquençage de 125-base dans Read 2.
Pour conclure, nous pensons que la méthode décrite pour être un outil précieux pour de nombreux types d’études. Des applications évidentes comprennent des enquêtes sur les mécanismes sous-jacents à l’inhibition de la transcription inverse grâce à des médicaments antirétroviraux ou facteurs de restriction cellulaires. Cependant, des ajustements relativement mineurs devraient être nécessaires d’adapter le système à 3'-termini cartographie au sein d’autres intermédiaires virales ADN simple brin, qui sont présents, par exemple, dans la réplication de parvovirus. En outre, le principe de la méthode, en particulier son étape de ligature optimisé, peut fournir une partie essentielle de la conception de préparation de bibliothèque pour la caractérisation de toute extension de 3'-ADN, y compris les allongements catalysées par cellulaire ADN à double brin polymérases.