Subskrypcja JoVE jest wymagana do oglądania tego materiału. Zaloguj się lub rozpocznij bezpłatny okres próbny.

Artykuł metodologiczny

Metodologia dokładnego wykrywania metylacji mitochondrialnego DNA

13.3K wyświetleń

DOI:

10.3791/57772

20 maja 2018

W tym artykule

Podsumowanie

Tutaj prezentujemy protokół umożliwiający dokładne ilościowe określenie metylacji mitochondrialnego DNA (mtDNA). W tym protokole opisujemy enzymatyczne trawienie DNA za pomocą BamHI w połączeniu z procesem analizy bioinformatycznej, który można wykorzystać do uniknięcia przeszacowania poziomów metylacji mtDNA spowodowanych drugorzędową strukturą mtDNA.

Streszczenie

Kwantyfikacja metylacji DNA może być osiągnięta za pomocą sekwencjonowania wodorosiarczynowego, które wykorzystuje właściwość wodorosiarczynu sodu do przekształcania niemetylowanej cytozyny w uracyl, w kontekście jednoniciowego DNA. Sekwencjonowanie wodorosiarczynów może być ukierunkowane (przy użyciu PCR) lub przeprowadzone na całym genomie i zapewnia bezwzględną kwantyfikację metylacji cytozyny w rozdzielczości pojedynczej zasady. Biorąc pod uwagę odrębny charakter DNA jądrowego i mitochondrialnego, zwłaszcza w strukturze drugorzędowej, należy dostosować metody sekwencjonowania wodorosiarczynów do badania metylacji cytozyny w mtDNA. Drugorzędowa i trzeciorzędowa struktura mtDNA może rzeczywiście prowadzić do artefaktów sekwencjonowania wodorosiarczynów, co prowadzi do wyników fałszywie dodatnich z powodu niepełnej denaturacji, słabego dostępu wodorosiarczynu do jednoniciowego DNA. W tym miejscu opisujemy protokół wykorzystujący enzymatyczne trawienie DNA za pomocą BamHI w połączeniu z procesem analizy bioinformatycznej, aby umożliwić dokładne ilościowe określenie poziomów metylacji cytozyny w mtDNA. Ponadto podajemy wytyczne dotyczące projektowania starterów sekwencjonowania wodorosiarczynów specyficznych dla mtDNA, aby uniknąć celowania w niepożądane segmenty miTochondrial NUclear (NUMT) wprowadzone do genomu jądrowego.

Wprowadzenie

Genom mitochondrialny to okrągła, dwuniciowa struktura o długości około 16,5 kilo podstawy (kb), składająca się z ciężkiej i lekkiej nici. Genom mitochondrialny jest obecny w wielu kopiach w każdej komórce, dziedziczony przez matkę i koduje niezbędne składniki kompleksów łańcucha oddechowego1. Podobnie jak genomy bakteryjne i w przeciwieństwie do genomu jądrowego, genom mitochondrialny jest zorganizowany w liczne struktury drugorzędowe i trzeciorzędowe, takie jak struktury zwinięte i superskręcone2, co może utrudnić dostęp podczas eksperymentów sekwencjonowania3.

W jądrze, metylacja DNA jest szeroko badanym znakiem epigenetycznym, który odgrywa rolę w wielu procesach, zwłaszcza w regulacji ekspresji genów. W genomach ssaków metylacja DNA zachodzi głównie w pozycji 5 pierścienia pirymidynowego deoksycytydyn, głównie na dinukleotydach CG (lub CpG). Metylacja cytozyny występuje w 70% wszystkich CpG w genomie komórek somatycznych i stanowi ~1% wszystkich zasad DNA4. Metylacje DNA zostały również opisane w kontekstach innych niż CpG, takich jak CpA, CpT i CpC i występują w różnych ilościach w jądrowym DNA, z wartościami do 25% wszystkich metylowanych cytozyn w embrionalnych komórkach macierzystych5,6,7.

Chociaż metylacja genomu jądrowego cytozyny jest powszechnie akceptowana, istnienie metylacji mitochondrialnego DNA (mtDNA) jest nadal kontrowersyjne. Pierwsze badanie dotyczące metylacji mtDNA przeprowadzono w hodowanych komórkach, w których metylacja mtDNA była łatwo wykrywalna, chociaż na niższych poziomach w porównaniu z jądrowym DNA8. Zarówno w komórkach ludzkich, jak i mysich metylację mtDNA wykryto również na niskich poziomach (2-5%). Korzystając z testów opierających się na wychwytywaniu 5 metylocytozyny, takich jak immunoprecypitacja metylowanego DNA (MeDIP), a następnie ilościowy PCR, metylację mtDNA wykryto również w różnych liniach myszy i ludzi oraz komórkach9,10,11,12. Używając przeciwciał przeciwko 5-metylocytozynie w teście ELISA lub spektrometrii mas, wykryto znaczne poziomy metylacji DNA z oczyszczonych frakcji mitochondrialnych13,14,15,16. Jednak większość testów w wyżej wymienionych badaniach wykorzystywała techniki, które nie zostały zaprojektowane w celu zapewnienia bezwzględnej kwantyfikacji metylacji DNA w rozdzielczości pojedynczej zasady.

Ilościowa i rezolutna analiza metylacji DNA może być osiągnięta za pomocą techniki zwanej "sekwencjonowaniem dwusiarczynowym", która wykorzystuje właściwość wodorosiarczynu sodu do przekształcania niemetylowanej cytozyny w uracyl w jednoniciowym DNA Context17. Korzystając z sekwencjonowania wodorosiarczynów, konstelacja badań wykryła obecność metylacji cytozyny na różnych poziomach. Metylacja mtDNA w regionie pętli D, regionie 12S lub 16S została łatwo wykryta w tkankach i komórkach human18,19,20,21,22,23 i mouse24 tkankach i komórkach, jednak z intrygującą zmiennością, wynoszącą 1-20% całkowitych cytozyn we wszystkich badaniach.

W porównaniu z tymi licznymi badaniami, tylko kilka badań, w tym z naszej grupy, kwestionowało obecność metylacji mtDNA3,25,26,27 lub kwestionowały biologiczne znaczenie bardzo niskich poziomów poziomów mtDNA (poniżej 2%)28. Niedawno informowaliśmy o obserwacji potencjalnego artefaktu sekwencjonowania wodorosiarczynów w całym sekwencjonowaniu wodorosiarczynów w mitochondriach3. Dostarczyliśmy dowodów na to, że drugorzędowa struktura mitochondrialnego DNA może prowadzić do wyników fałszywie dodatnich w sekwencjonowaniu wodorosiarczynów, przeszacowując w ten sposób poziomy metylacji. Przedstawiamy tutaj protokół zapobiegający artefaktowi konwersji mtDNA przez wodorosiarczyn. Protokół ten wykorzystuje proste enzymatyczne trawienie DNA w celu zakłócenia drugorzędowych struktur mtDNA i umożliwienia pełnego dostępu do wodorosiarczynu po protokole sekwencjonowania wodorosiarczynu. Ponadto zapewniamy towarzyszący proces bioinformatyczny do analizy sekwencjonowania wodorosiarczynów.

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Protokół

1. Traktowanie enzymami restrykcyjnymi

  1. Zlinearyzuj mtDNA, traktując całkowite ludzkie DNA enzymem restrykcyjnym BamHI, który tnie w pozycji 14258 w ludzkim mitochondrialnym DNA.
    UWAGA: W przypadku DNA mysiego należy użyć enzymu restrykcyjnego BglII w identycznych warunkach, jak opisano poniżej.
  2. Dla każdej próbki, w której ma zostać oceniona metylacja mtDNA, przygotuj jedną probówkę reakcyjną o pojemności 0,2 mL i dodaj następującą mieszaninę: 3 μg genomic DNA (ilość określona metodą fluorometrii), 15 μL Buffer 3, 3 μL BamHI oraz wodę do uzyskania objętości 150 μL.
  3. Umieść probówki w termocyklerze na 4 h w temperaturze 37 °C.

2. Konwersja pirosiarczynowa

  1. Przeprowadzić konwersję DNA traktowanego enzymem BamHI za pomocą pirosiarczynu sodu zgodnie z opisem zamieszczonym w innym opracowaniu29.
  2. Do każdej reakcji konwersji użyć 200-500 ng DNA traktowanego enzymem BamHI w całkowitej objętości 20 μL. Odzysk DNA wynosi 50-70%.
  3. Eluować DNA po konwersji pirosiarczynem w 10 μL buforu do elucji.
  4. Ilościowo oznaczyć jednoniciowe DNA za pomocą fluorometrii.
  5. Opcjonalnie: DNA po konwersji pirosiarczynem przechowywać w temperaturze -20 °C przed przystąpieniem do pozostałych etapów protokołu. Do długotrwałego przechowywania użyć temperatury -80 °C

3. Projektowanie starterów do sekwencjonowania z pirosiarczynem

  1. Zaprojektuj startery do sekwencjonowania po konwersji pirosiarczynowej, korzystając z narzędzi online Methprimer30 oraz BiSearch31,32.
    UWAGA: Methprimer i BiSearch umożliwiają wyszukiwanie sekwencji wiązania starterów w sekwencjach genomicznych, w których cytozyny zostały przekształcone w tyminy, analogicznie do procesu po konwersji pirosiarczynowej.
  2. Aby uzyskać optymalną i nieobciążoną amplifikację, unikaj dinukleotydów CpG w starterach i zaprojektuj amplikon o rozmiarze od 100 do 300 bp.
  3. Upewnij się, że zaprojektowane startery są specyficzne dla mtDNA, korzystając z funkcji Primer Search w programie BiSearch. Wprowadź wcześniej zaprojektowane startery do DNA po konwersji pirosiarczynowej, zaznacz pole bisulfite oraz podaj genom referencyjny i parametry PCR.
    UWAGA: Wyświetlona zostanie lista możliwych produktów PCR.
  4. Skopiuj 1–5 najlepszych sekwencji starterów i wklej je do formularza zamówienia syntezy oligonukleotydów.
    UWAGA: Lista zwalidowanych w laboratorium sekwencji starterów do mtDNA człowieka i myszy znajduje się w Tabeli 1.
  5. Przetestuj startery, wykonując PCR z DNA po konwersji pirosiarczynowej, a następnie elektroforezę w żelu agarozowym zgodnie z opisem w kroku 4. Przetestuj i zoptymalizuj wszystkie pary starterów oddzielnie i upewnij się, że na żelu agarozowym pojawia się amplikon o prawidłowym rozmiarze.
    UWAGA: Jeśli wymagane są startery do PCR multipleksowego, dodaj kilka starterów do jednej reakcji PCR i zwizualizuj produkty PCR za pomocą elektroforezy w żelu agarozowym lub metodą o wyższej rozdzielczości, np. elektroforezą w żelu poliakrylamidowym, aby rozróżnić produkty o podobnej wielkości.

4. PCR po konwersji pirosiarczynowej i ekstrakcja z żelu

  1. Aby przeprowadzić amplifikację obszarów zainteresowania, wykonaj reakcję PCR z użyciem polimerazy Taq typu hot start. W skrócie: wymieszaj 100 ng przekształconego DNA, 500 µM starterów sensu i anty-sensu, 1 µL mieszaniny dNTP (po 10 µM każdego dNTP) oraz polimerazę w ilości 7,5 U/reakcję w całkowitej objętości 50 µL. Przeprowadź PCR w następujących warunkach: 5 min w 95 °C (60 s 94 °C; 60 s 55 °C*; 60 s 72 °C) Schemat równowagi statycznej, ΣFx=0, MA=0; siły, równowaga momentów sił; ilustracja edukacyjna z fizyki. 35 cykli; 10 min 72 °C. Użyj starterów osobno lub w układzie multipleksowym.
  2. Opcjonalnie: Podczas stosowania starterów w układzie multipleksowym użyj 200 ng przekształconego DNA i następujących warunków cyklicznych: 5 min 95 °C (60 s 94 °C; 90 s 55 °C*; 90 s 72 °C) Schemat równowagi statycznej, ΣFx=0, MA=0; siły, równowaga momentów sił; ilustracja edukacyjna z fizyki. 35 cykli; 10 min 72 °C.
    UWAGA: Temperatura może się różnić w zależności od temperatury topnienia starterów.
  3. Produkty PCR poddaj elektroforezie w 2% żelu agarozowym przy napięciu 100 V.
  4. Pod słabym światłem UV wytnij produkty PCR skalpelem i przenieś je do mikroprotek. Nie wystawiaj produktów PCR na nadmierne działanie światła UV, aby uniknąć uszkodzenia DNA. Unikaj czasu ekspozycji przekraczającego 30 s.
  5. Oczyść produkty PCR metodą oczyszczania z żelu, zgodnie z opisem w punkcie3.
  6. Oznacz ilość oczyszczonego DNA z kroku 4.5 za pomocą fluorometrii. Przejdź do kroku 5 lub przechowuj próbki w -20°C.

5. Przygotowanie biblioteki do sekwencjonowania bisulfitowego

  1. Przeprowadź przygotowanie biblioteki z produktów PCR po konwersji bisulfitowej.
    Opcjonalnie: produkty multiplex PCR na tym etapie.
  2. Przeprowadzić naprawę końców, wykorzystując do 100 ng produktu PCR. Dodać 3 µL enzym do przygotowania końców i 6.5 µL Bufor do reakcji naprawy końców (10x) do 55.5 µL produkt PCR w probówce 0,2 ml. Umieścić probówkę w termocyklerze i inkubować przez 30 min w temperaturze 20 °C a następnie 30 min w temperaturze 65 °C.
  3. Ligować adaptery sekwencyjne, mieszając DNA po naprawie końców z 15 µL mieszanina ligazy Blunt/TA, 2.5 µL Adapter i 1 µL Wzmacniacz ligacji. W przypadku stosowania <100 ng produktu PCR jako materiał wyjściowy, rozcieńczyć adaptor 10-krotnie.
  4. Umieść mieszaninę (krok 5.3) w termocyklerze i inkubuj przez 15 min w temperaturze 20 °CWstrzymaj pracę termocyklera i dodaj 3 µl Dodać enzym do probówki. Wymieszać i umieścić probówkę z powrotem w termocyklerze, a następnie inkubować przez 15 min w temperaturze 37 °C.
  5. Przeprowadź selekcję wielkości DNA z przyligowanymi adapterami, wykorzystując kule do odwracalnej immobilizacji w fazie stałej (SPRI), stosując różne proporcje kul do DNA w zależności od wielkości produktu PCR.
  6. Dodaj 13.5 µL H2Do tego dodaj DNA z ligowanymi adaptorami z etapu 5.4 i 55 µL resuspender kulki SPRI. Inkubować w temperaturze pokojowej przez 5 min, a następnie umieścić probówkę na statywie magnetycznym. Po całkowitym oczyszczeniu roztworu przenieść nadsącz do nowej probówki i odrzucić probówkę z kulkami.
    UWAGA: Nadfiltr zawiera DNA z przyłączonymi adaptorami, a duże, niepożądane fragmenty są związane z odrzuconymi kuleczkami.
  7. Dodaj 25 µL dodać resuspenderowane kulki SPRI do nadsączu z kroku 5.6, wymieszać i inkubować przez 5 min w temperaturze pokojowej. Umieścić probówkę na statywie magnetycznym i odpiplikować nadsącz, gdy roztwór stanie się klarowny.
    UWAGA: Odrzucony nadsącz zawiera niepożądane DNA, natomiast DNA z przyłączonymi adapterami jest związane z kuleczkami.
  8. Gdy probówka znajduje się na statywie magnetycznym, należy dodać 200 µL 80% etanol (świeżo przygotowany) do przemycia kulek. Inkubować przez 30 s, a następnie usunąć i odrzucić nadsącz. Powtórzyć ten krok, aby w sumie wykonać 2 przemycia.
  9. Suszyć kuleczki na powietrzu przez 5 min, pozostawiając probówkę z otwartą zakrętką na statywie magnetycznym.
  10. Wyjmij probówkę z magnesu, a następnie eluuj DNA do 23 µL bufor do elucji i wymieszać. Inkubować w temperaturze pokojowej przez 2 min.
  11. Umieść probówkę na statywie magnetycznym, a gdy roztwór stanie się klarowny, przenieś 2 µL do płytki PCR z 96 dołkami do wstępnej amplifikacji PCR (krok 5.14).
  12. Przenieść pozostałą część w przybliżeniu 21 µL do nowej probówki o pojemności 0,2 ml do amplifikacji PCR (krok 5.17).
    UWAGA: Należy upewnić się, aby nie przenieść żadnych kulek, ponieważ mogą one hamować reakcje enzymatyczne w kolejnych etapach.
  13. Przejdź do kroku 5.14 lub przechowuj próbki w temperaturze -20°C.
  14. Przeprowadź amplifikację przed-PCR metodą RT-qPCR w celu oszacowania liczby cykli niezbędnych do amplifikacji DNA z ligowanymi adaptorami. Do każdej reakcji w 96-dołkowej płytce PCR dodaj następujące składniki: 2 µL DNA (krok 5.11): 0.4 µL Starter indeksowy, 0.4 µL uniwersalny primer PCR, 7.2 µL H2O, 10 µL Master mix do RT-qPCR.
  15. Umieść płytkę w urządzeniu do PCR w czasie rzeczywistym i poddaj ją następującym warunkom: 30 s w temperaturze 98 °C (10 s 98 °C; 75 s 65 °C) Schemat równowagi statycznej, ΣFx=0, MA=0; siły, równowaga momentów sił; ilustracja dydaktyczna z fizyki. 20 cykli; 5 min 65 °C.
  16. Oszacuj liczbę cykli amplifikacji niezbędną dla każdej próbki, odejmując od wartości Ct uzyskanej w preamplifikacyjnej reakcji RT-qPCR wartość Ct odpowiadającą końcowi fazy liniowej reakcji PCR.
  17. Przeprowadź amplifikację DNA z przyligowanymi adapterami z kroku 5.12, mieszając: 21 µL DNA, 25 µL gotowa mieszanina do PCR (master mix), 1 µL starter indeksowy, 1 µL Uniwersalna starter do PCR i 2 µL H2O. W termocyklerze poddać każdą próbkę następującym warunkom: 30 s w temperaturze 98 °C (10 s 98 °C; 75 s 65 °C) Schemat równowagi statycznej, ΣFx=0, MA=0; siły, równowaga momentów sił; ilustracja edukacyjna z fizyki. [obliczona wartość Ct z kroku 5.16] cykli; 5 min 65 °C.
    UWAGA: Należy pamiętać, aby do każdej próbki użyć tylko jednego startera indeksującego, co umożliwi multipleksowanie. Indeks jest wprowadzany podczas reakcji PCR.
  18. Oczyść namnożone DNA za pomocą kulek SPRI i eluuj do 22 µL bufor elucyjny
  19. Kontroluj jakość biblioteki za pomocą elektroforezy żelowej lub alternatywnej metody. Sprawdź, czy rozmiary pików biblioteki są prawidłowe (rozmiar produktu/produktów PCR plus adapter).
  20. Oszacuj średnią wielkość par zasad produktu/produktów biblioteki za pomocą analizy rozmycia (smear analysis): W zaawansowanych ustawieniach globalnych kliknij dwukrotnie w „Tablepodczas analizy rozmycia (smear analysis). Zdefiniuj zakres od 100 do 1000 bp i kliknij Proszę dostarczyć tekst źródłowy do tłumaczenia.Pod Tabela regionówpojawia się zdefiniowany obszar i obliczana jest średnia wielkość biblioteki (bp).
  21. Oznacz ilość bibliotek metodą fluorometrii. Przejdź do kroku 6 lub przechowuj biblioteki w temperaturze -20 °C.

6. Sekwencjonowanie nowej generacji

  1. Oblicz stężenie bibliotek w nM, korzystając ze wzoru:
    Wzór obliczania stężenia DNA; konwersja masy cząsteczkowej na jednostki nanomolarne; równanie.
  2. Rozcieńcz biblioteki do tej samej koncentracji w nM, np. 2 nM (wybierz 4 nM, 2 nM, 1 nM lub 0,5 nM w zależności od stężeń bibliotek). Połącz wszystkie biblioteki, aby uzyskać pulę bibliotek o stężeniu 2 nM.
  3. Zdenaturować i rozcieńczyć biblioteki zgodnie z instrukcją obsługi sekwenconera. W skrócie: zdenaturować pulę 2 nM poprzez dodanie 0,2 N NaOH i inkubować w temperaturze pokojowej przez 5 min. Rozcieńczyć zdenaturowaną pulę do końcowego stężenia 10 pM. Zdenaturować i rozcieńczyć dostępną komercyjnie bibliotekę kontrolną (Control library) do końcowego stężenia 12,5 pM. W nowej probówce zmieszać pulę 10 pM z 20% biblioteki kontrolnej 12,5 pM.
    UWAGA: Konwersja pirosiarczynowa wprowadza bardzo niską złożoność. Dodatek 20% biblioteki kontrolnej zapewni lepszą jakość sekwencjonowania.
  4. Przeprowadzić sekwencjonowanie paired-end 150 bp przy użyciu urządzenia do głębokiego sekwencjonowania.

7. Analiza obliczeniowa – szacowanie poziomów metylacji

  1. Wygeneruj pliki .fastq (tutaj nazwane sample1_R1.fastq.gz oraz sample1_R2.fastq.gz), wygeneruj indeks Bismark dla całego analizowanego genomu (tutaj w folderze o nazwie bismarkIndex) i upewnij się, że sekwencja genomowa chrM jest dostępna w formacie fasta (tutaj w folderze o nazwie chrM).
  2. Przeprowadź wstępne przetwarzanie odczytów w celu usunięcia adapterów i błędów (bias) wprowadzonych przez primery: użyj narzędzia Trim_galore33. Usuń dwa nukleotydy z końca 5' zarówno z odczytów forward, jak i reverse.
    trim_galore --clip_R1 2 --clip_R2 2 -o ./ --trim1 --paired sample1_R1.fastq.gz sample1_R2.fastq.gz
  3. Zmapuj wstępnie przetworzone odczyty do całego genomu przy użyciu Bismark34
    bismark --dovetail --bam -o ./ ./bismarkIndex / -1 ./sample1_R1_val_1.fq.gz -2 ./sample1_R2_val_2.fq.gz
    UWAGA: Jeśli używany jest inny program do mapowania, należy upewnić się, że odczyty, których nie można jednoznacznie zmapować, zostaną odrzucone.
  4. Wyodrębnij odczyty mapujące się do chromosomu M, tutaj przy użyciu Samtools35
    samtools sort -l 0 -O BAM -o sample1_sorted.bam
    sample1_R1_val_1_bismark_bt2_pe.bam
    samtools index sample1_sorted.bam
    samtools view -u sample1_sorted.bam chrM | samtools sort -n -O BAM -o sample1_chrM.bam
  5. Wyekstrahuj informacje o metylacji
    bismark_methylation_extractor -p --CX --cytosine_report --comprehensive --genome_folder ./chrM ./sample1_chrM.bam
  6. Wykorzystaj pliki .bedGraph, .cov lub CX_report.txt do dalszej analizy. Użyj pliku sample1_chrM.CX_report.txt do wizualizacji i testów.
  7. Jeśli analizowane jest wiele regionów i na wykresach M-bias wygenerowanych podczas mapowania widoczny jest błąd wprowadzony przez primery, przeprowadź dodatkowe przycinanie (clipping) na etapie wstępnego przetwarzania. Więcej informacji znajduje się w dokumentacji Bismark.

8. Analiza obliczeniowa - Test różnic

  1. Upewnij się, że plik CX_report.txt zawiera 7 kolumn: chromosom (tutaj chrM), pozycję, nić, liczbę odczytów zmetylowanych, liczbę odczytów niezmetylowanych, kontekst C oraz sekwencję otaczającą dla każdej cytyzyny (C) na chromosomie M.
  2. Ponieważ raport CX_report obejmuje cały chrM, wydziel podzbiór odpowiadający regionowi (regionom), które zostały poddane amplifikacji.
  3. Do wykazania różnic między próbkami zastosuj testy nieparametryczne, takie jak dokładny test Fishera dla poszczególnych C lub test znakowy dla całego regionu.
    UWAGA: Poziom metylacji często będzie bliski 0%, dlatego założenie o rozkładzie normalnym jest nieodpowiednie.
  4. Podobnie w przypadku wizualizacji, przedstaw kwantyle lub oblicz przedziały ufności dla proporcji dwumianowej.

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Wyniki

Dwa kroki w tym protokole są kluczowe podczas badania metylacji mtDNA. 1) Rozplecenie struktury drugorzędowej oraz 2) zaprojektowanie starterów specyficznych dla mitochondrialnego DNA.

Poprzez trawienie ludzkiego genomowego DNA za pomocą enzymu restrykcyjnego BamHI (Rysunek 1), struktura mitochondrialnego DNA zostanie przecięta na pozycji nukleotydowej 14 258, a struktura drugorzędowa zostanie otwar...

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Dyskusja

W tym miejscu udostępniamy protokół sekwencjonowania wodorosiarczynów, który jest specjalnie zaprojektowany do badania metylacji mtDNA. Różnice w stosunku do protokołów sekwencjonowania wodorosiarczynów stosowanych do genomowego DNA polegają na wykorzystaniu wcześniejszego etapu trawienia enzymem restrykcyjnym i analizie bioinformatycznej wykluczającej fałszywie dodatnie wyniki wynikające z sekwencji NUMT.

Zapewniamy protokół zapobiegający artefaktom sekwencjonowania wodorosiarczynów podczas b...

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Oświadczenia

Autorzy nie mają żadnych konkurencyjnych interesów finansowych, które mogliby zadeklarować.

Podziękowania

Centrum Podstawowych Badań Metabolicznych Fundacji Novo Nordisk jest niezależnym ośrodkiem badawczym na Uniwersytecie Kopenhaskim, częściowo finansowanym z nieograniczonej darowizny od Fundacji Novo Nordisk.

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Materiały

Lista materiałów użytych w tym artykule
NazwaFirmaNumer katalogowyKomentarze
BamHINew England BioLabs# R0136
EZ Zestaw do metylacji DNA i błyskawicyZymo Research# D5030 
Test Qubit ssDNAThermo Fisher Scientific# Q10212
Probówki do testów kubitowychThermo Fisher Scientific# Q32856
HotStarTaq plus zestaw polimerazy DNAQiagen# 203603 
Zestaw do ekstrakcji żelu QIAquickQiagen# 28704
NEBNext Ultra Zestaw biblioteki DNA dla IlluminaNew England BioLabs# E7370S
NEBNext Multiplex Oligoes dla IlluminaNew England BioLabs# E7335S
Koraliki AMPure XPBeckman Coulter# A63881
Chip DNA o wysokiej czułościAgilent# 5067-4626
Test dsDNA o wysokiej czułości kubituThermo Fisher Scientific# Q33230
Zestaw odczynników MiSeq v2 300 cykliIllumina# MS-102-2003
PhiX control v3Illumina# FC-110-3001
Wodorotlenek sodu Sigma# S5881
Termocykler C1000Biorad# 1851148
CFX96 System wykrywania PCR w czasie rzeczywistymBiorad #1855195
Fluorometr kubitowyThermo Fisher Scientific# Q33226
Bioanalizator 2100Agilent# G2939BA
MiSeq instrumentIllumina# SY-410-1003

Bibliografia

  1. Falkenberg, M., Larsson, N. -G., Gustafsson, C. M. DNA replication and transcription in mammalian mitochondria. Annu Rev Biochem. 76, 679-699 (2007).
  2. Kolesar, J. E., Wang, C. Y., Taguchi, Y. V., Chou, S. -H., Kaufman, B. A. Two-dimensional intact mitochondrial DNA agarose electrophoresis reveals the structural complexity of the mammalian mitochondrial genome. Nucleic Acids Res. 41 (4), e58(2013).
  3. Mechta, M., Ingerslev, L. R., Fabre, O., Picard, M., Barres, R. Evidence Suggesting Absence of Mitochondrial DNA Methylation. Front Genet. 8, 166(2017).
  4. Ehrlich, M., Gama-Sosa, M. A., et al. Amount and distribution of 5-methylcytosine in human DNA from different types of tissues of cells. Nucleic Acids Res. 10 (8), 2709-2721 (1982).
  5. Lister, R., Pelizzola, M., et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 462 (7271), 315-322 (2009).
  6. Yan, J., Zierath, J. R., Barres, R. Evidence for non-CpG methylation in mammals. Exp Cell Res. 317 (18), 2555-2561 (2011).
  7. Patil, V., Ward, R. L., Hesson, L. B. The evidence for functional non-CpG methylation in mammalian cells. Epigenetics. 9 (6), 823-828 (2014).
  8. Nass, M. K. Differential methylation of mitochondrial and nuclear DNA in cultured mouse, hamser and virus trasnforemd hamser cells in vivo and in vitro methylation. J Mol Biol. 80 (1), 155-175 (1973).
  9. Shock, L. S., Thakkar, P. V., Peterson, E. J., Moran, R. G., Taylor, S. M. DNA methyltransferase 1, cytosine methylation, and cytosine hydroxymethylation in mammalian mitochondria. PNAS. 108 (9), 3630-3635 (2011).
  10. Bellizzi, D., D'Aquila, P., et al. The Control Region of Mitochondrial DNA Shows an Unusual CpG and Non-CpG Methylation Pattern. DNA Res. 20 (6), 537-547 (2013).
  11. Jia, Y., Li, R., et al. Maternal Low-Protein Diet Affects Epigenetic Regulation of Hepatic Mitochondrial DNA Transcription in a Sex-Specific Manner in Newborn Piglets Associated with GR Binding to Its Promoter. PLoS ONE. 8 (5), e63855(2013).
  12. Jia, Y., Song, H., Gao, G., Cai, D., Yang, X., Zhao, R. Maternal Betaine Supplementation during Gestation Enhances Expression of mtDNA-Encoded Genes through D-Loop DNA Hypomethylation in the Skeletal Muscle of Newborn Piglets. J Agric Food Chem. 63 (46), 10152-10160 (2015).
  13. Infantino, V., Castegna, A., Iacobazzi, F., Spera, I. Impairment of methyl cycle affects mitochondrial methyl availability and glutathione level in Down's syndrome. Mol Genet Metab. 102 (3), 378-382 (2011).
  14. Chen, H., Dzitoyeva, S., Manev, H. Effect of valproic acid on mitochondrial epigenetics. Eur J Pharmacol. 690 (1-3), 51-59 (2012).
  15. Dzitoyeva, S., Chen, H., Manev, H. Effect of aging on 5-hydroxymethylcytosine in brain mitochondria. Neurobiol Aging. 33 (12), 2881-2891 (2012).
  16. Menga, A., Palmieri, E. M., et al. SLC25A26 overexpression impairs cell function via mtDNA hypermethylation and rewiring of methyl metabolism. FEBS J. 284 (6), 967-984 (2017).
  17. Frommer, M., McDonald, L. E., et al. A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. PNAS. 89 (5), 1827-1831 (1992).
  18. Byun, H. -M., Panni, T., et al. Effects of airborne pollutants on mitochondrial DNA Methylation. Part Fibre Toxicol. 10 (1), 1(2013).
  19. Janssen, B. G., Byun, H. -M., Gyselaers, W., Lefebvre, W., Baccarelli, A. A., Nawrot, T. S. Placental mitochondrial methylation and exposure to airborne particulate matter in the early life environment: An ENVIRONAGE birth cohort study. Epigenetics. 10 (6), 536-544 (2015).
  20. Zheng, L. D., Linarelli, L. E., et al. Insulin resistance is associated with epigenetic and genetic regulation of mitochondrial DNA in obese humans. Clin Epigenetics. 7 (1), 1739(2015).
  21. Byun, H. -M., Colicino, E., Trevisi, L., Fan, T., Christiani, D. C., Baccarelli, A. A. Effects of Air Pollution and Blood Mitochondrial DNA Methylation on Markers of Heart Rate Variability. J Am Heart Assoc. 5 (4), (2016).
  22. Bianchessi, V., Vinci, M. C., et al. Mitochondrion. MITOCH. 27 (C), 40-47 (2016).
  23. Wijst, M. G. P., van Tilburg, A. Y., Ruiters, M. H. J., Rots, M. G. Experimental mitochondria-targeted DNA methylation identifies GpC methylation, not CpG methylation, as potential regulator of mitochondrial gene expression. Sci Rep. 7 (1), 177(2017).
  24. Wong, M., Gertz, B., Chestnut, B. A., Martin, L. J. Mitochondrial DNMT3A and DNA methylation in skeletal muscle and CNS of transgenic mouse models of ALS. Front cellr Neurosci. 7, 279(2013).
  25. Dawid, I. B. 5-methylcytidylic acid: absence from mitochondrial DNA of frogs and HeLa cells. Science. 184 (4132), 80-81 (1974).
  26. Gama-Sosa, M. A. Levels and distribution of 5-methylcytosine in Chordate DNA. , 1-201 (1985).
  27. Hong, E. E., Okitsu, C. Y., Smith, A. D., Hsieh, C. -L. Regionally specific and genome-wide analyses conclusively demonstrate the absence of CpG methylation in human mitochondrial DNA. Mol Cell Biol. 33 (14), 2683-2690 (2013).
  28. Liu, B., Du, Q., et al. CpG methylation patterns ofhuman mitochondrial DNA. Nature Publishing Group. , 1-10 (2016).
  29. Donkin, I., Versteyhe, S., et al. Obesity and Bariatric Surgery Drive Epigenetic Variation of Spermatozoa in Humans. Cell Metab. 23 (2), 369-378 (2016).
  30. Li, L. C., Dahiya, R. MethPrimer: designing primers for methylation PCRs. Bioinformatics. 18 (11), 1427-1431 (2002).
  31. Tusnády, G. E., Simon, I., Váradi, A., Arányi, T. BiSearch: primer-design and search tool for PCR on bisulfite-treated genomes. Nucleic Acids Res. 33 (1), e9(2005).
  32. Arányi, T., Váradi, A., Simon, I., Tusnády, G. E. The BiSearch web server. BMC bioinformatics. 7, 431(2006).
  33. Krueger, F. Trim Galore!. , Available from: http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ (2018).
  34. Krueger, F., Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 27 (11), 1571-1572 (2011).
  35. Li, H., Handsaker, B., et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 25 (16), 2078-2079 (2009).
  36. Ramos, A., Barbena, E., et al. Nuclear insertions of mitochondrial origin: Database updating and usefulness in cancer studies. Mitochondrion. 11 (6), 946-953 (2011).
  37. Warnecke, P. M., Stirzaker, C., Song, J., Grunau, C., Melki, J. R., Clark, S. J. Identification and resolution of artifacts in bisulfite sequencing. Methods (San Diego, Calif). 27 (2), 101-107 (2002).

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Przedruki i uprawnienia

Tagi

Sekwencjonowanie z u yciem pirosiarczyn wtrawienie enzymem BamHIunikanie NUMTprojektowanie starter welektroforeza elowakuleczki SPRIprzygotowanie bibliotekianaliza Bioanalyzeremkwantyfikacja metylacji