brmlab

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== Rendo ======


{{template>:project:infobox|
name=Rendo|
image=:rendo:img_6194.jpg?400|
founder= [[user:magn3zy]]|
interested=[[user:sumie-dh]],[[user:sachy]],[[user:joe]]|
sw=FIXME|
hw=N/A|
status=active}}
~~META:
status = in progress
&relation firstimage = :project:rendo:img_6194.jpg
~~

===== Vědecká hypotéza a cíle práce =====

==== Vědecká hypotéza ====

Restrikční analýza (RFLP) patří k nejstarším metodám využívaným k odlišení různých alel, SNV či indelů. Pomocí softwaru, který by analyzoval nejenom několik různých sekvencí, ale také další variabilitu z referenčních verzí genomů a spolehlivých relevatních sekvenačních dat z GenBank by bylo možné automatizovaně navrhovat vhodné restrikční endonukleázy k analýze konkrétních alel, SNV a indelů.

==== Cíle práce ====

  - Naprogramovat software umožňující navrhnutí vhodných restrikčních endonukleáz k analýze uživatelem zadaného úseku dle referenční verze genomu a dalších dostupných sekvencí k doplnění variability v zadaném úseku.
  - Vybrat vhodné testovací úseky genomu modelových organismů (//Arabidopsis thaliana//).
  - Ověřit funkčnost navrhnutých enzymů softwarem pro konkrétní analýzu variability.

===== Metodika =====

==== Tvorba softwaru ====

Aplikace pro návrh vhodných restrikčních endonukleáz byla programována pomocí programovacího jazyku Python a VS Code IDE, poskytujícího mnoho vhodných knihoven pro zajištění fukčnosti a možnosti stahování dat z online databází. Aplikace byla strukturována na 3 dílčí aplikace umožňující základní a pokročilou analýzu. Aplikace také obsahuje 259 restikčních endonukleáz, které jsou využívány ke všem analýzám. Tento počet byl určen na základě průzkumu komerčně nabízených enzymů firmami specializujími se na jejich výrobu jako Thermo Fisher Scientific a New England Biosciences. Mnoho restrikčních endonukleáz je prodáváno pod více komerčními názvy, proto bylo zajištěno, aby pro jednu rozpoznávací sekvenci a pozici řezu byla k dispozici pouze jedna restrikční endonukleáza pod nejznámějším názvem.

==== Výběr testovacích genů a návrh primerů ====

Dle referenční verze genomu //Arabidopsis thaliana// TAIR10 byly identifikovány vysoce konzervované oblasti genomu. Jako testovací geny byly vybrány geny ovlivňující dobu kvetení (FLC, FRI) a také geny z rodiny cytokinin dehydrogenáz/oxidáz (CKX1 - CKX7).

Primery byly navrhnuty pomocí softwaru BLAST a dále byly analyzovány softwarem Primer3. Primery byly dále vybrány i s ohledem na podobné anelační teploty během PCR ke zmešení počtu prováděných PCR reakcí, všechny až na jednu vyjímku mají optimální anelační teplotu v rozmezí 50 - 53°C. Primery jsou zobrazeny v //Tab.// níže.


^Název^Sekvence 5′‑3′^Chromozom^RefSeq pozice^PCR produkt^
|FLC5‑F|ACAGAGACAACACAGAAGAGG|5|3 174 992|132|
|FLC5‑R|ACAACCGTGCTGCTTTTGTT|5|3 175 104|--|
|FLC35‑F|GGCAGTCTCAAGGTGTTCCT|5|3 175 198|472|
|FLC35‑R|CTCTCCAGCCTGGTCAAGAT|5|3 175 650|--|
|FRI2‑F|TCAAGTATCAAGCGTGGAATGC|4|270 369|184|
|FRI2‑R|TGTGCAAGGGAAGGGTTCAT|4|270 533|--|
|CKX13‑F|GGTTAAATGGATCAGGGTACTCTAC|2|17 315 739|328|
|CKX13‑R|GCCAACGGTTGATTTGGAGC|2|17 316 047|--|
|CKX42‑F|GGAAAAGGTGAAATGATGACTTGCT|4|14 567 352|120|
|CKX42‑R|TGCATGATCCAACGCAATCC|4|14 567 452|--|
|CKX33‑F|CGGACCTCGTTGACTGTGTAT|5|23 045 707|268|
|CKX33‑R|AGGCCAAGTGGTTAAGGTTTCTA|5|23 045 952|--|
|CKX512‑F|GAGCATCTCATCACCTATCCT|1|28 317 349|327|
|CKX512‑R|CGTCCGACCGGATGAAATGT|1|28 317 656|--|
|CKX73‑F|CGGACGGCTAAGATTTAATGAAGG|5|7 228 755|217|
|CKX73‑R|TGGACCATTGACTCCGTTCT|5|7 228 952|--|


Dle referenční verze genomu \textit{Zea mays} Zm-B73-REFERENCE-NAM-5.0 byly identifikovány vysoce konzervované oblasti genomu. Jako testovací geny byly vybrány geny EF1A, TB1, ADH, TGA1, HMGR1. Primery byly navrhnuty pomocí softwaru BLAST a dále byly analyzovány softwarem Primer3. V //Tab.// jsou všechny využívané primery zobrazeny.

^Název^Sekvence 5′‑3′^Chromozom^RefSeq pozice^PCR produkt^
|EF1A‑F|TGGGTTTTGCTGTTTTACAGATGG|—|—|990|
|EF1A‑R|CGCGACAACTGCAAGGATG|—|—|--|
|EF1AS‑F|TCATCATGAACCACCCTGGC|—|—|440|
|EF1AS‑R|GTTCACGCGACAACTGCAAG|—|—|--|
|ADH1‑F|AATACCACAGCTAAGGACGC|—|—|302|
|ADH1‑R|GAGAGGGTGTGACTGACGTA|—|—|--|
|TB1‑F|GCACAGCAAGATATGCACCG|—|—|853|
|TB1‑R|AGCTCACGTTGAGGGAAAGG|—|—|--|
|HMGR1‑F|TGACGCCTGTCTTTGACTCG|—|—|701|
|HMGR1‑R|CACTGGCCTAGTATGGACGC|—|—|--|
|TGA11‑F|CATCGCACCTCACACCTCTC|—|—|732|
|TGA11‑R|ACAGGGACGAGCTACACTACA|—|—|--|



==== Materiál ====

Jednotlivé rostliny Arabidopsis thaliana a Zea mays byly získány na polních lokalitách v Středočeském a Karlovarském kraji. Z jedné lokality byly vždy vyizolovány 2 až 3 rostliny.

=== Izolace DNA ===

Genomická DNA //Arabidopsis thaliana// byla izolována za pomoci izolačního kitu Plant Genomic DNA Fast Kit (Mebep Biosicence, China) dle protokolu přiloženého ke kitu za využití mírných modifikací doporučených výrobcem:

  - Do zkumavky typu Eppendorf 1,5 ml bylo umístěno 100 mg rozmělněné tkáně.
  - Do zkumavky bylo přidáno 400 μl AP1 Bufferu a 4 μl RNase A a následně byl obsah vortexován.
  - Zkumavka byla umístěna do inkubátoru na 10 minut při teplotě 65°C a během inkubace byl obsah 3 krát promíchán.
  - Po inkubaci bylo do zkumavky přidáno 130 μl AP2 Bufferu, obsah zkumavky byl promíchán, poté byla zkumavky na 5 minut nechána na ledu a poté byla zkumavka centrifugována 5 minut při 13 000 rpm.
  - Supernatant ze zkumavky byl přenesen do nové zkumavky typu Eppendorf 1,5 ml. Supernatant byl opět centrifugován 5 minut při 13 000 rpm a přenesen opět do nové zkumavky pro zajištění vyšší čistoty DNA.
  - Do zkumavky byl přidán isopropanol pokojové teploty v objemu 70% získaného supernatantu, zkumavka byla 30 krát obracením promíchána.
  - Zkumavka byla centrifugována 2 minuty při 12 000 rpm a následně byl odstraněn supernatant.
  - Do zkumavky byl přidán 1 ml 70% ethanolu a zkumavka byla několikrát promíchána a následně centrifugována 1 minutu při 12 000 rpm. Supernatant byl odstraněn a DNA byla nechána odvětrat.
  - Do zkumavky s DNA izolátem byl přidán Buffer DD a DNA byla uložena do -20°C pro dlouhodobé skladování.

Genomická DNA kukuřice byla izolována pomocí klasické fenol-chloroformové metody pomocí předem připravených roztoků dle následujícího protokolu:

    - Příprava CTAB lyzačního pufru: do kádinky nalejte 30ml deionizovené vody a přidejte: 1g CTAB; 4,09 NaCl; 0,6055g Tris; 0,372g EDTA. Roztok doplňte na celkový objem 50ml deionizovanou vodou. Upravte pomocí koncentrované HCL na pH 8.0. Zahřejte roztok na 65°C a promíchejte - roztok by měl být čirý.
    - Příprava rehydratačního TE pufru: Do autoklávovatelné nádoby nalejte 100ml deionizované vody a přídejte 1,211g Tris a 0,372 EDTA. Roztok vyautoklávujte.
    - Připavte si roztoky RNázy A a Proteinázy K v 1ml TE pufru v následující koncentraci: RNáza A 10mg/ml a Proteináza K 20mg/ml.
    - Příprava roztoku fenol:chloroform:isoamylakohol v poměru 25:24:1 finální množství připravte dle počtu prováděných izolací - na jednu izolaci je třeba 800μl.
    - Příprava roztoku chloroform:isoamylalkohol v poměru 24:1 finální množství připravte dle počtu prováděných izolací - na jednu izolaci je třeba 600μl.
    - Do zkumavky typu Eppendorf 1,5 ml bylo umístěno 100 mg rozmělněné tkáně.
    - K rozmělněným mladým listům přidejte 800μl CTAB lyzačního roztoku a 2μl Proteinázy K. Zvortexujte a 30 minut inkubujte při 65°C.
    - Přidejte 800μl fenol:chloroform:isoamylalkoholu a 12krát invertujte. Centifugujte po dobu 10minut při 12 000rpm.
    - Do nové zkumavky 1,5 ml typu Eppendorf přeneste přibližně 600μl vodné fáze a přidejte 1μl RNázy A a inkubujte 15 minut při 37°C.
    - Přidejte 600μl chloroform:isoamylalkoholu invertujte 10krát a centrifugujte 5 minut při 12 000rpm.
    - Přesuňte 500μl vodné fáze do nové zkumavky a přidejte 300μl isopropanolu. Invertujte 10krát a nechte 10 minut při pokojové teplotě. Poté centrifugujte 15 minut při 12 000rpm.
    - Odstraňte supernatant a přidejte 1ml 70% ethanolu, krátce promíchejte a centrifugujte 5 minut při 12 000rpm.
    - Odstraňte ethanol a peletu DNA nechte 5 minut odvětrat.
    - Přidejte 50 až 100μl TE rehydratačního roztoku a nechte DNA rozpustit (při 70°C 30 minut; při 4°C přes noc) a poté skladujte při teplotě -20°C.


==== Testování amplifikace vybraných úseků ====

Pro vybrané primery pro amplifikaci vybraných úseků k restrikční analýze bylo třeba ověřit amplifikační schopnost těchto primerů. Proto vždy jeden vzorek z jedné lokality byl testován všemi primery. Ty které amplifikovali byly dále testovány u všech jedinců ze všech lokalit. Produkty PCR reakcí byly ověřovány na gelu pomocí agarózové elektroforézy, ověření amplifikace správného úseku bylo provedeno porovnáním produktů PCR s 100bp DNA Ladderem (GDSBio, China).

=== PCR ===

První kolo testování amplifikační schopnosti 8 párů primerů pro //Arabidopsis thaliana// a 6 párů primerů pro //Zea mays// bylo provedeno pro vždy jednu rostlinu z každé lokality. Všechny reakce byly provedeny v reakčním objemu 25 ul. Byla pro to použita DreamTaq polymeráza (Thermo Fisher Scientific, USA). PCR reakce byly prováděny pomocí termocykleru Thermo Hybaid PXE 0.2 (Thermo Electron Corporation, USA).

Druhá vlna testování amplifikační schopnost primerů zahrnovala dle konkrétního problému při žádné amplifikaci, či při vzniku nespecifických PCR produktů využití jiné DNA polymerázy. Využívána k tomu byla Pfu DNA Polymerase (GDSBio, China) a také Taq Plus DNA Polymerase (GDSBio, China) obsahující kombinaci vysoce specifické DNA polymerázy a také Pfu DNA polymerázy či FirePol DNA Polymerase (Solis BioDyne, Estonia). Podmínky při kterých byly tyto PCR reakce prováděny jsou zobrazeny v //Tab.//.

^Fáze           ^Teplota (°C)        ^Délka (s)  ^Počet opakování  ^
|Predenaturace  |95                  |300        |1                |
|Denaturace     |95                  |30         |33               |
|Anelace        |50, 51, 52, 53, 56  |30, 45, 60 |33               |
|Elongace       |72                  |60         |33               |
|Postelongace   |72                  |600        |1                |





Poslední vlna testování amplifikace dosud neoptimalizovaných podmínek pro jednotlivé primery zahrnovala gradientní testování anelačních teplot vždy po 2 stupních pro konkrétní primer a dále také prodloužení anelační fáze. Amplifikující primery při optimalizovaných podmínkách byly využity pro všechny rostliny a jejich amplifikace byla ověřena na gelu.

PCR amplifikace poté byla provedena pro všechny amplifikující primery a pro všechny vzorky. Tyto vzorky nebyly již kontrolovány na gelu, byly dále využity pro štěpení restrikčními endonukleázami pro samotné testování vytvořeného RENDO softwaru.

=== Gelová alektroforéza ===

První kolo testování aplifikační schopností vybraných PCR primerů bylo vždy testováno u jedné rostliny z jedné lokality. Délka výsledných PCR produktů se výrazněji odlišovala, proto byly testovány na 1,4% agarózovém gelu. Pro přípravu agarózového gelu byla používána PCR Agaróza (Top-Bio, Czech Republic) a 1x TAE pufr. Samotná gelová elektroforéza byla provedena pomocí SGE-014-02 Econosubmarine Electrophoresis Gell Kit (CBS Scientific, USA). Gely byly obarveny barvivem DSView Nucleic Acid Stain 20,000X (GDSBio, China). Zobrazení výsledku elektroforézy bylo prováděno pomocí transiluminátoru LB0100 (Invitrogen by Thermo Fisher Scientific, USA). Gely byly fotografovány v tmavé komoře pomocí fotoaparátu Canon 6D MarkII.

Druhé kolo testování primerů proběhlo pro všechny primery, které amplifikovali pro všechny rostliny. Využito bylo stejné přístrojové a chemické vybavení jako v prvním kole testování, lišilo se pouze procentuální zastoupení agarózy využité pro přípravu gelu, dle potřebné přesnosti rozlišení.

==== Restrikční analýza ====

=== Štěpení ===

Amplifikující vybrané úseky byly využity pro analýzu restrikčními endonukleázami. Dle referenční verze genomu a dostupných dat z databáze GenBank byly vybrány místa vykazující potvrzenou variabilitu i v České republice. Tyto místa a amplifikovaný úsek genomu byl pomocí softwaru RENDOVarPro analyzován a navržené enzymy tímto softwarem byly využity k štěpení produktů PCR reakcí.

=== Gelová elektroforéza ===

Výsledky štěpení PCR produktů navrhnutými endonukleázami byly analyzovány pomocí gelové elektroforézy. Zjištěná variabilita byla ověřena porovnáním vzniklých fragmentů s 100bp DNA Ladderem (GDSBio, China) umožňující ověření softwarem předpovězených výsledků restrikční analýzy.

===== Výsledky =====

==== Tvorba softwaru ====

=== RENDOSim ===

Uživatelské rozhraní této dílčí aplikace je zobrazeno na obrázku níže. Uživatel může nahrát sekvenci ve formátu fasta ze svého zařízení. Dále si z nabídky může vybrat jednu z restrikčních endonukláz. V nabídce jsou rozděleny na endonukleázy tvořící tupé nebo lepivé konce a seřazeny jsou abecedně pro rychlé vyhledání konkrétní restrikční endonukleázy.

{{ :project:rendo:rendosim.png}}

Po nahrání vlákna a zvolení endonukleázy započne analýza sekvence stlačením tlačítka „Cut“ a lineárně je vlákno prohledáno na rozpoznávací sekvenci zvolené restrikční endonukleázy. V případě nalezení sekvence je vlákno rozřezáno dle pozice řezu pro vybranou endonukleázu a výsledky jsou zobrazeny uživateli. Pokud nedojde k nalezení rozpoznávací sekvence je uživateli zobrazeno pouze původní vlákno. Před dalším vyhledáváním je doporučeno si uživatelské rozhraní zresetovat pomocí tlačítka „Clear“, ale není to nutné je možné pouze změnit vybranou restrikční endonukleázu, či změnit nahraný fasta soubor a analyzovat jiný soubor stejnou endonukleázou.

=== RENDOVar ===

Dílčí aplikace RENDOVar a její uživatelské rozhraní je zobrazeno na obrázku níže. Tato aplikace je určená pro rozlišení známých variant jednotlivých genomických úseků, klonů či geneticky modifikovaných organismů. Zde není zohledňována variabilita, kromě té kterou softwaru ve formátu fasta souborů poskytnete. Po nahrání libovolného počtu fasta souborů tlačítkem „Analyze“ započnete analýzu. V okně s výsledky se zobrazí všechny restikční endonukleázy nabízející rozlišení jednotlivých variant. Seřazeny jsou dle vzniklého poměrového počtu fragmentů od nejlépe rozlišitelných možností po ty hůře rozlišitelné. 

{{ :project:rendo:rendovar.png}}

=== RENDOVarPro ===

Nejdůležitější dílčí aplikací je aplikace RENDOVarPro. Zohledňující variabilitu z konkrétní referenční verze genomu a také variabilitu z dalších dostupných sekvenačních dat z GenBank. Uživatel zadá úsek který by pomocí PCR amplifikoval, název organismu, referenční verzi genomu, pozici indelu či SNV kterou chce analyzovat. Důležité je analyzované pozice zadávat s ohledem na referenční verzi genomu dle, které se pracuje pro stáhnutí správných a relevantních souborů. Pro přesnější analýzu je třeba zadat také číslo chromozomu. Dále musí uživatel zadat svůj email, tento email není nijak uchováván, slouží pouze jako parametr požadavku pro NCBI Entrez API, která dle emailu zajišťuje kontrolu dodržení počtu požadavků.

Po vyplnění těchto údajů začne analýza nejprve stažením z NCBI databáze pomocí Entrez API referenční sekvence, dále je z EVA databáze stažen celogenomický vcf soubor pro zadaný organismus a referenční verzi genomu. Z těchto dat je poté vytvořen vcf soubor zahrnující variabilitu z referenční verze genomu. Poté probíhá návrh vhodných restrikčních endonukleáz k analýze SNV. Uživatelské rozhraní RENDOVarPro aplikace je zobrazeno na obrázku níže. Samotná analýza probíhá různě dlouhou dobu. Způsobeno je to různou velikostí celogenomového vcf souboru a také rychlostí připojení k internetu. Pokud pro zjištění konkrétní variability neexistuje vhodný enzym je v poli pro výsledky zobrazeno hlášení „No enzymes found“.

{{ :project:rendo:rendovarpro.png}}

RENDOVarPro je však ještě stále vylepšována, nejdůležitějším úkolem je vytvoření vizualizace výsledků na gelu při štěpení konkrétním enzymem pro všechny varianty a rozšíření této dílčí aplikace o jistý manuální a polomanuální mód umožňující uživatelům více zasahovat do nastavení a parametrů hledání vhodných restrikčních endonukleáz. Po pilotních testech softwaru bude tento software celý dostupný na GitHub jako veřejný opensource nástroj.

==== Restrikční analýza ====

=== Arabidopsis thaliana ===

Softwarem RENDOVarPro byly navrhnuty vhodné restrikční enzymy pro analýzu konkrétních SNV, jednotlivé úseky genů uvedených a štepených na obrázku byly tímto enzymem štěpeny dle metodiky.

{{ :project:rendo:13_5_2025_combined_for_article.jpg?600 |}}

===== Licence =====

Copyright © Eliška Korbová 2025

Tento software a jeho dokumentace jsou chráněny autorským právem. Všechna práva vyhrazena.

Tento projekt je ve fázi vývoje a není určen k veřejnému použití. Jakékoli šíření, kopírování, úpravy, prezentace nebo použití tohoto softwaru (ani jeho částí) bez výslovného písemného souhlasu autora jsou zakázány.

Použití v akademických, výukových, komerčních i nekomerčních kontextech není momentálně dovoleno.

Autor si vyhrazuje právo změnit podmínky šíření, konkrétně zveřejněním projektu na GitHub přejde projekt pod licenci Creative Commons Attribution-NonCommercial 4.0 International (CC BY-NC 4.0).

Kontakt: eliskakorbova@gmail.com



Copyright © Eliška Korbová 2025

This software and its documentation are protected by copyright. All rights reserved.

This project is currently under development and is not intended for public use. Any distribution, copying, modification, presentation, or use of this software (or any part thereof) without the express written consent of the author is strictly prohibited.

Use in academic, educational, commercial, or non-commercial contexts is currently not permitted.

The author reserves the right to change the terms of distribution. Specifically, upon publication of the project on GitHub, the project will be released under the Creative Commons Attribution-NonCommercial 4.0 International (CC BY-NC 4.0) license.

Contact: eliskakorbova@gmail.com

== Vysvětlivky ==

Tento projekt se nachází ve vývojové fázi, a proto v současnosti není určen k distribuci ani veřejnému využití. Po úspěšném dokončení pilotních laboratorních testů bude plánováno jeho zveřejnění jako open-source nástroj určený pro nekomerční použití.
Autor tímto krokem chrání původní nápad a zároveň si ponechává možnost využití výsledků projektu v rámci akademických prací. Oddálením veřejného zpřístupnění je také zajištěna vyšší spolehlivost softwaru při jeho budoucím použití.



This project is currently in the development phase and is therefore not intended for distribution or public use. After successful completion of initial laboratory testing, it is planned to be released as an open-source tool for non-commercial use.
This approach protects the original idea and allows the author to potentially use the project results in academic work. Delaying public availability also helps ensure greater software reliability upon future deployment.


===== Poděkování =====

Ráda bych poděkovala zaměstnancům Katedry genetiky a šlechtění ČZU za možnost využívání jejich laboratoře a přístrojového vybavení, které zde není dostupné (spektrofotometr pro mikrolitrové objemy, restrikční endonukleázy).