Name: Odkrycie nowych sekwencji za pomocą genomiki subtraktywnej (Wideo) | JoVE Journal
Uploaded: 2019-01-25T13:30:00.000+00:00
Duration: 9 min 40 s
Description: Odkryj genomikę subtraktywną wykorzystującą darmowe oprogramowanie do identyfikacji nowych sekwencji genomowych trudnych do wyizolowania. Ta elastyczna metoda komputerowa znajduje zastosowanie w poszukiwaniu celów szczepionkowych bakterii i nieznanych wirusów.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

Autorzy dziękują Michelle Biederman, Alyssie Pedersen i Colinowi J. Saldanha za pomoc w projekcie genomiki zeberki na różnych etapach. Wyrażamy również uznanie dla Evgeny'ego Biska za administrowanie systemem klastrów obliczeniowych oraz grant NIH 1K22CA184297 (dla J.R.B.) i NIH NS 042767 (dla C.J.S.).

1. De novo Zmontuj Sekwencję Początkową
UWAGA: Dowolne dane sekwencji nowej generacji (NGS) mogą być użyte, o ile można z nich wyprodukować zespół. Odpowiednie dane wejściowe obejmują odczyty Illumina, PacBio lub Oxford Nanopore złożone w plik fasta. Dla ścisłości, w tej sekcji opisano zestaw transkryptomiczny oparty na Illumina, specyficzny dla przeprowadzonego przez nas badania na ziębie zeberce7; Należy jednak pamiętać, że szczegóły będą się różnić w zależności od projektu. W naszym przykładowym projekcie surowe dane pochodziły z MiSeq i uzyskano około 10 milionów sparowanych odczytów z każdej próbki.
<ol>
	<li>Użyj Trimmomatic 0.3214, aby usunąć adaptery Illumina i podstawy niskiej jakości. W wierszu polecenia wpisz: 
	java -jar trimmomatic-0.32.jar PE -phred33 forward.fq.gz reverse.fq.gz -baseout quality_and_adaptor_trimmed ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 PROWADZENIE:3 KOŃCOWE:3 SLIDINGOKNO:4:20 MINLEN:40</li>
	<li>Użyj PEAR15 v. 0.9.6, aby utworzyć wysokiej jakości scalone odczyty z sparowanych odczytów wyjściowych trimmomatic, używając domyślnych parametrów. W wierszu polecenia wpisz: 
	gruszka -f <quality_and_adaptor_trimmed_1P.fastq> -r <quality_and_adaptor_trimmed_2P.fastq></li>
	<li>Użyj Reptile v. 1.116, aby poprawić błędy odczytów generowanych przez PEAR. Postępuj zgodnie z protokołem krok po kroku opisanym w17.</li>
	<li>Użyj Trinity v. 2.4.018 w trybie domyślnym, aby złożyć poprawione sekwencje. W przypadku bibliotek specyficznych dla pasma należy użyć parametru -SS_lib_type . Wyjściem jest plik fasta (your_assembly.fasta). W wierszu polecenia wpisz: 
	Trinity --seqType fq --SS_lib_type FR –max_memory 10G –output Trinity_output --left quality_and_adaptor_trimmed_forward_paired_reads.fq –right quality_and_adaptor_trimmed_reverse_paired_reads.fq –CPU 10 
	UWAGA: Dane wyjściowe zostaną umieszczone w nowym katalogu, Trinity_output, a zestaw zostanie nazwany "Trinity.fasta", który w razie potrzeby można zmienić na Your_assembly.fasta. Więcej informacji można znaleźć na stronie internetowej Trinity: https://github.com/trinityrnaseq/trinityrnaseq/wiki/Running-Trinity.</li>
</ol>
2. WYSADZAJ zestaw w stosunku do sekwencji referencyjnej
UWAGA: Użyj tego kroku, gdy odniesienie jest zestawem lub długimi odczytami, jak Sanger; jeśli składa się z surowych odczytów Illumina, zobacz krok 3 poniżej, aby odwzorować odczyty do zapytania. Wszystkie kroki BLAST zostały wykonane w wersji 2.2.29+, chociaż polecenia powinny działać na każdej najnowszej wersji BLAST.
<ol>
	<li>Utwórz bazę danych BLAST sekwencji referencyjnej (nucleotide_reference.fasta) w wierszu poleceń. Wprowadź w wierszu poleceń następujące informacje: 
	makeblastdb -dbtype nucl -in nucleotide_reference.fasta -out nucleotide_reference.db</li>
	<li>BLAST — dopasuj zestaw zapytania (wygenerowany w kroku 1) do referencyjnej bazy danych. Aby uzyskać plik wyjściowy, użyj [-out BLAST_results.txt], a aby wygenerować tabelaryczne dane wyjściowe (wymagane dla kolejnych kroków przetwarzania za pomocą skryptów Pythona), użyj [-outfmt 6]. Opcje te można łączyć w dowolnej kolejności, więc przykładowym kompletnym poleceniem jest [blastn -query your_assembly.fasta -db nucleotide_reference.db -out BLAST_results.txt -outfmt 6]. Jeśli wymagane jest ustawienie wartości e, użyj opcji -evalue z odpowiednią liczbą, na przykład [-evalue 1e-6]. Należy jednak pamiętać, że cykl subtraktywny skutecznie odwraca ustawienie evalue zgodnie z opisem w dyskusji.</li>
	<li>Aby zwiększyć rygorystyczność, użyj sekwencji białek z zestawu jako zapytania BLAST z przetłumaczonym nukleotydem BLAST (tBLASTn), który wykonuje 6-kierunkową translację bazy danych (nukleotydów). Ta metoda jest zalecana dla większości systemów niemodelowych, unikając problemu niekompletnych adnotacji białek.
	<ol>
		<li>Upewnij się, że wybrano prawidłowy kod genetyczny dla badanego organizmu, korzystając z opcji -db_gencode. Aby uzyskać sekwencje białek dla zapytania, uruchom polecenie TransDecoder.LongOrfs (z pakietu TransDecoder v. 3.0.1) w celu zidentyfikowania najdłuższych otwartych ramek odczytu ze złożonych sekwencji zapytań. Polecenie to [TransDecoder.LongOrfs -t your_assembly.fasta]; Dane wyjściowe zostaną umieszczone w katalogu o nazwie "transcripts.transdecoder_dir" i będą zawierać plik o nazwie longest_orfs.pep zawierający najdłuższe przewidywane sekwencje białek z każdej sekwencji w your_assembly.fasta.</li>
		<li>Aby użyć tBLASTn, uruchom polecenie [tblastn -query longest_orfs.pep -db nucleotide_reference.db -out BLAST_results.txt -outfmt 6]. Jeśli dostępne jest wysokiej jakości białko referencyjne, należy użyć dopasowania białko-białko z BLASTp, a nie tBLASTn.</li>
		<li>Utwórz bazę danych BLAST z odniesieniem do białka [makeblastdb -dbtype prot -in protein_reference.fasta -out protein_reference.db], a następnie [blastp -query longest_orfs.pep -db protein_reference.db -out BLAST_results.txt -outfmt 6]. Upewnij się, że zapisujesz wyniki jako plik do dalszego przetwarzania i użyj tabelarycznego (outfmt 6), aby upewnić się, że skrypty języka Python mogą je poprawnie przeanalizować.</li>
	</ol>
	</li>
</ol>
3. Mapa jest wczytywana do zespołu
UWAGA: Ta metoda może być użyta, jeśli zestaw danych referencyjnych składa się z surowych odczytów genomowych, a nie złożonych sekwencji lub sekwencji Sangera, w którym to przypadku użyj BLAST (krok 2.1).
<ol>
	<li>Korzystając z BWA -MEM v. 0.7.1219 lub bowtie220, zamapuj pobrane nieprzetworzone odczyty (raw_reads.fastq) na zestaw zapytania. Dane wyjściowe będą w formacie .sam. Polecenia są następujące: najpierw indeksuj zestaw: [bwa index your_assembly.fasta], a następnie mapuj odczyty [bwa mem your_assembly.fasta raw_reads.fastq >mapped.sam]. (Zauważ, że symbol ">" nie jest tutaj znakiem większości niż; zamiast tego instruuje wyjście, aby trafiło do pliku mapped.sam).</li>
</ol>
4. Użyj skryptu Pythona, aby usunąć wszystkie pasujące sekwencje
UWAGA: Dostarczone skrypty działają z Pythonem 2.7.
<ol>
	<li>Wykonując krok 2, użyj subtraktywnego skryptu Pythona, używając polecenia [./Non-matching_sequences.py your_assembly.fasta BLAST_results.txt]. Przed uruchomieniem skryptu upewnij się, że plik wyjściowy BLAST jest w formacie 6 (tabelarycznym). Skrypt wypisze plik z niepasującymi sekwencjami w formacie fasta o nazwie your_assembly.fasta_non-matching_sequences_BLAST_results.txt.fasta, a także pasującymi sekwencjami dla rekordów, jako your_assembly.fasta_matching_sequences_BLAST_results.txt.fasta. Najważniejszy będzie niepasujący plik, jako źródło potencjalnych sekwencji T do badań i dalszych cykli genomiki subtraktywnej.</li>
	<li>Wykonując krok 3, uruchom skrypt języka Python removeUnmapped.py, aby przyjąć jako dane wejściowe plik sam z kroku 3.1 i zidentyfikować nazwy sekwencji zapytań bez żadnych pasujących odczytów i zapisać je w nowym pliku tekstowym. Użyj polecenia [./removeUnmapped.py mapped.sam], a dane wyjściowe zostaną mapped.sam_contigs_with_no_reads.txt. (Program wygeneruje odchudzony plik sam z usuniętymi wszystkimi niezmapowanymi odczytami; plik ten może być zignorowany dla celów tego protokołu, ale może być przydatny do innych analiz).</li>
	<li>Ponieważ danymi wyjściowymi poprzedniego kroku jest lista nazw sekwencji w pliku tekstowym o nazwie mapped.sam_contigs_with_no_reads.txt, wyodrębnij plik fasta z tymi sekwencjami: [./getContig.py your_assembly.fasta mapped.sam_contigs_with_no_reads.txt]. Wynikiem będzie plik o nazwie mapped.sam_contigs_with_no_reads.txt.fasta.</li>
</ol>
5. Zaprojektuj podstawy dla sekwencji, która pozostała
UWAGA: W tym momencie znajduje się plik fasta zawierający proponowane sekwencje T. W tej sekcji opisano qPCR w celu eksperymentalnego sprawdzenia, czy pochodzą one z T, czy z wcześniej nieznanych regionów R. Jeśli odejmowanie w kroku 4 usunęło wszystkie sekwencje, oznacza to, że albo początkowy zestaw nie zawierał T, albo odejmowanie mogło być zbyt rygorystyczne.
<ol>
	<li>Użyj Geneious21, aby ręcznie określić optymalne sekwencje starterów.
	<ol>
		<li>Wyróżnij sekwencję kandydującą 21-28 pz dla startera Forward. Unikaj runów składających się z 4 lub więcej dowolnej bazy. Spróbuj skierować się na region z dość jednolitą kombinacją wszystkich par zasad. Pojedynczy G lub C na końcu 3 stóp jest korzystny, pomagając zakotwiczyć podkład.</li>
		<li>Kliknij kartę Statystyki po prawej stronie ekranu, aby wyświetlić szacowaną temperaturę topnienia (Tm) tej sekwencji, gdy zaznaczony jest obszar kandydujący. Staraj się uzyskać temperaturę topnienia w zakresie 55-60 °C, unikając przy tym powtórzeń i długich przebiegów G/C.</li>
		<li>Wykonaj kroki 5.1.1. i 5.1.2, aby wybrać starter odwrotny, umieszczony 150-250 par zasad w odległości 3' od startera przedniego. Chociaż długości startera nie muszą się zgadzać, przewidywana Tm powinna być jak najbardziej zbliżona do Tm startera przedniego. Pamiętaj, aby odwrócić uzupełnianie sekwencji (jeśli klikniesz prawym przyciskiem myszy w Geneious, gdy sekwencja jest podświetlona, jest to opcja menu).</li>
	</ol>
	</li>
	<li>Użyj funkcji Primer Design, która znajduje się na górnym pasku narzędzi w oknie sekwencji.
	<ol>
		<li>Kliknij przycisk Projekt podkładu. Wstaw region do wzmocnienia w obszarze Region docelowy.</li>
		<li>Na karcie Charakterystyka wprowadź żądany rozmiar, temperaturę topnienia (Tm) i %GC (patrz krok 5.1.1.).</li>
		<li>Kliknij przycisk OK, aby wygenerować startery. Zamów podkłady za pośrednictwem niestandardowej usługi oligo.</li>
	</ol>
	</li>
	<li>Walidacja starterów z kontrolnym DNA (kodującym zarówno T, jak i R) w celu optymalizacji Tm i czasu wydłużenia. Użyj zwykłej elektroforezy Taq i żelu, aby zobaczyć rozmiar pasma, ale optymalizację można również przeprowadzić za pomocą qPCR zgodnie z metodami opisanymi w kroku 6.
	<ol>
		<li>Wykonać 10-krotne rozcieńczenia zarówno starterów do przodu, jak i do tyłu, tak aby startery miały stężenie 10 μM.</li>
		<li>Użyj mieszaniny PCR składającej się z 0,5 μl dNTP, 0,5 μl startera do przodu, 0,5 μl startera odwrotnego, 0,1 μl polimerazy Taq, 2 μl matrycy, 0,75 μl magnezu, 2,5 μl buforu i 18,15 μl wody, tak aby na matrycę przypadało 25 μl o stężeniu 5 ng/μl.</li>
		<li>Przetestuj startery w różnych temperaturach topnienia w programie PCR. Zazwyczaj optymalną wydajność obserwuje się w temperaturach topnienia nieco poniżej przewidywanej Tm podkładów, ale zwykle nie powyżej 60 °C. Przetestuj również optymalne czasy wydłużenia, korzystając z tego przewodnika: 1 minuta na 1000 pz (a więc zwykle 10-30 sekund w zależności od długości amplikonu).</li>
		<li>Wykonaj elektroforezę żelową w punkcie końcowym, aby potwierdzić, że startery wzmacniają oczekiwaną sekwencję. Uruchom 25 μl produktu qPCR zmieszanego z 5 μl barwnika 6X glicerolu na 2% żelu agarozowym TAE przy 200 V przez 20 minut.</li>
	</ol>
	</li>
</ol>
6. Walidacja pozostałej sekwencji metodą qPCR
UWAGA: Ten krok wymaga walidacji starterów i warunków PCR ustalonych w kroku 5.
<ol>
	<li>Uruchomić każdą matrycę w trzech egzemplarzach z następującą mieszanką: 12,5 μl mieszanki wzorcowej PowerSYBR Green, 0,5 μl startera o stężeniu 10 μM, 0,5 μl startera odwrotnego o stężeniu 10 μM, 10,5 μl wody i 1 μl matrycy DNA (w stężeniu 2 ng/μL), tak, aby każda studzienka zawierała 25 μl całkowitej objętości.</li>
	<li>Uruchom program qPCR w oparciu o zwalidowaną temperaturę i czas przedłużenia z kroku 4. Zaprojektowaliśmy i zwalidowaliśmy wszystkie podkłady tak, aby były kompatybilne z cyklem dwuetapowym, 95 °C przez 10 minut początkowego topnienia, następnie 40 cykli po 95 °C przez 30 s i 60 °C przez 1 min. Jednak trzystopniowy program (topnienie-wyżarzanie-rozciąganie) może być bardziej optymalny dla podkładów i powinien być dostosowany w razie potrzeby. Zalecamy, aby końcowe krzywe denaturacji były generowane co najmniej przy pierwszym zastosowaniu starterów w qPCR w celu walidacji amplifikacji pojedynczego produktu DNA.</li>
	<li>Zmierz zielone sygnały qPCR/SYBR w stosunku do aktyny (lub innej odpowiedniej kontroli "R") za pomocą Ct. Dla wszystkich przypadków oblicz średnią i odchylenie standardowe 2- (gen Ct - β-aktyna Ct).</li>
	<li>(Opcjonalnie) Wykonaj elektroforezę żelową w punkcie końcowym, aby potwierdzić prawidłowe wykrywanie rozmiaru produktu za pomocą qPCR. Tutaj uruchom 25 μl produktu qPCR zmieszanego z 5 μl 6x barwnika glicerolowego na 2% żelu agarozowym TAE przy 200 V przez 20 minut.</li>
</ol>
7. Powtórz z nowym odniesieniem, aby zmniejszyć dane.
UWAGA: Jeśli krok 6 potwierdził tożsamość sekwencji z T, zakończ cykl tutaj (Rysunek 2A). Jednak różne względy mogą motywować do kontynuacji cyklu, na przykład, jeśli w pliku pozostaje wiele sekwencji R lub jeśli żadna z kandydujących sekwencji T nie została zwalidowana przez qPCR w kroku 6.
<ol>
	<li>Uzyskaj nowe odniesienie. Ten krok umożliwia nową iterację cyklu i może obejmować surowe dane genomiczne, surowe dane sekwencyjne RNA lub inne zgromadzone zestawy danych. Cennymi zasobami danych referencyjnych są baza danych genomu w National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/genome), w której przechowywane są zebrane genomy dostępne za pośrednictwem FTP (ftp://ftp.ncbi.nlm.nih.gov/genomes/), oraz Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/), w którym przechowywane są surowe odczyty sekwencji nowej generacji. Projekty genomowe mogą dostarczać swoje surowe dane sekwencyjne za pośrednictwem innych stron internetowych i baz danych związanych z projektem.</li>
</ol>

<ol><li>Barh, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+Novel+Comparative+Genomics+Analysis+for+Common+Drug+and+Vaccine+Targets+in+Corynebacterium+pseudotuberculosis+and+other+CMN+Group+of+Human+Pathogens%5BTitle%5D)+AND+%22Chemical+Biology+%26+Drug+Design%22%5BJournal%5D)">A Novel Comparative Genomics Analysis for Common Drug and Vaccine Targets in Corynebacterium pseudotuberculosis and other CMN Group of Human Pathogens.</a> Chemical Biology & Drug Design. 78 (1), 73-84 (2011).</li>
<li>Sarangi, A. N., Aggarwal, R., Rahman, Q., Trivedi, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Subtractive+Genomics+Approach+for+in+Silico+Identification+and+Characterization+of+Novel+Drug+Targets+in+Neisseria+Meningitides+Serogroup+B%5BTitle%5D)+AND+%22Journal+of+Computer+Science+%26+Systems+Biology%22%5BJournal%5D)">Subtractive Genomics Approach for in Silico Identification and Characterization of Novel Drug Targets in Neisseria Meningitides Serogroup B.</a> Journal of Computer Science & Systems Biology. 2 (5), (2009).</li>
<li>Kaur, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+Druggable+Targets+for+Acinetobacter+baumannii+Via+Subtractive+Genomics+and+Plausible+Inhibitors+for+MurA+and+MurB%5BTitle%5D)+AND+%22Applied+Biochemistry+and+Biotechnology%22%5BJournal%5D)">Identification of Druggable Targets for Acinetobacter baumannii Via Subtractive Genomics and Plausible Inhibitors for MurA and MurB.</a> Applied Biochemistry and Biotechnology. 171 (2), 417-436 (2013).</li>
<li>Rathi, B., Sarangi, A. N., Trivedi, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genome+subtraction+for+novel+target+definition+in+Salmonella+typhi%5BTitle%5D)+AND+%22Bioinformation%22%5BJournal%5D)">Genome subtraction for novel target definition in Salmonella typhi.</a> Bioinformation. 4 (4), 143-150 (2009).</li>
<li>Epstein, J. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+GBV-D%2C+a+Novel+GB-like+Flavivirus+from+Old+World+Frugivorous+Bats+%28Pteropus+giganteus%29+in+Bangladesh%5BTitle%5D)+AND+%22PLoS+Pathogens%22%5BJournal%5D)">Identification of GBV-D, a Novel GB-like Flavivirus from Old World Frugivorous Bats (Pteropus giganteus) in Bangladesh.</a> PLoS Pathogens. 6 (7), (2010).</li>
<li>Kapoor, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+Rodent+Homologs+of+Hepatitis+C+Virus+and+Pegiviruses%5BTitle%5D)+AND+%22MBio%22%5BJournal%5D)">Identification of Rodent Homologs of Hepatitis C Virus and Pegiviruses.</a> MBio. 4 (2), (2013).</li>
<li>Biederman, M. K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Discovery+of+the+First+Germline-Restricted+Gene+by+Subtractive+Transcriptomic+Analysis+in+the+Zebra+Finch%2C+Taeniopygia+guttata%5BTitle%5D)+AND+%22Current+Biology%22%5BJournal%5D)">Discovery of the First Germline-Restricted Gene by Subtractive Transcriptomic Analysis in the Zebra Finch, Taeniopygia guttata.</a> Current Biology. 28 (10), 1620-1627 (2018).</li>
<li>Readhead, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multiscale+Analysis+of+Independent+Alzheimer%27s+Cohorts+Finds+Disruption+of+Molecular%2C+Genetic%2C+and+Clinical+Networks+by+Human+Herpesvirus%5BTitle%5D)+AND+%22Neuron%22%5BJournal%5D)">Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus.</a> Neuron. 99, 1-19 (2018).</li>
<li>Carroll, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+global+virome+project%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">The global virome project.</a> Science. 359 (6378), 872-874 (2016).</li>
<li>Pigozzi, M. I., Solari, A. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Germ+cell+restriction+and+regular+transmission+of+an+accessory+chromosome+that+mimics+a+sex+body+in+the+zebra+finch.+Taeniopygia+guttata%5BTitle%5D)+AND+%22Chromosome+Research%22%5BJournal%5D)">Germ cell restriction and regular transmission of an accessory chromosome that mimics a sex body in the zebra finch. Taeniopygia guttata.</a> Chromosome Research. 6, 105-113 (1998).</li>
<li>Itoh, Y., Kampf, K., Pigozzi, M. I., Arnold, A. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Molecular+cloning+and+characterization+of+the+germline-restricted+chromosome+sequence+in+the+zebra+finch%5BTitle%5D)+AND+%22Chromosoma%22%5BJournal%5D)">Molecular cloning and characterization of the germline-restricted chromosome sequence in the zebra finch.</a> Chromosoma. 118, 527-536 (2009).</li>
<li>Warren, W. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+genome+of+a+songbird%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">The genome of a songbird.</a> Nature. 464, 757-762 (2010).</li>
<li>Balakrishnan, C. N., Lin, Y. C., London, S. E., Clayton, D. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNAseq+transcriptome+analysis+of+male+and+female+zebra+finch+cell+lines%5BTitle%5D)+AND+%22Genomics%22%5BJournal%5D)">RNAseq transcriptome analysis of male and female zebra finch cell lines.</a> Genomics. 100, 363-369 (2012).</li>
<li>Bolger, A. M., Lohse, M., Usadel, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Trimmomatic%3A+a+flexible+trimmer+for+Illumina+sequence+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Trimmomatic: a flexible trimmer for Illumina sequence data.</a> Bioinformatics. 30 (15), 2114-2120 (2014).</li>
<li>Zhang, J., Kobert, K., Flouri, T., Stamatakis, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PEAR%3A+a+fast+and+accurate+Illumina+Paired-End+reAd+mergeR%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">PEAR: a fast and accurate Illumina Paired-End reAd mergeR.</a> Bioinformatics. 30, 614-620 (2014).</li>
<li>Yang, X., Dorman, K. S., Aluru, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reptile%3A+representative+tiling+for+short+read+error+correction%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Reptile: representative tiling for short read error correction.</a> Bioinformatics. 26, 2526-2533 (2010).</li>
<li>MacManes, M. D., Eisen, M. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Improving+transcriptome+assembly+through+error+correction+of+high-throughput+sequence+reads%5BTitle%5D)+AND+%22PeerJ%22%5BJournal%5D)">Improving transcriptome assembly through error correction of high-throughput sequence reads.</a> PeerJ. 1 (113), (2013).</li>
<li>Grabherr, M. G., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Full-length+transcriptome+assembly+from+RNA-seq+data+without+a+reference+genome%5BTitle%5D)+AND+%22Nature+Biotechnology%22%5BJournal%5D)">Full-length transcriptome assembly from RNA-seq data without a reference genome.</a> Nature Biotechnology. 29, 644-652 (2011).</li>
<li>Li, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Aligning+sequence+reads%2C+clone+sequences+and+assembly+contigs+with+BWA-MEM%5BTitle%5D)+AND+%22arXiv%22%5BJournal%5D)">Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM.</a> arXiv. , (2013).</li>
<li>Langmead, B., Salzberg, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fast+gapped-read+alignment+with+Bowtie+2%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Fast gapped-read alignment with Bowtie 2.</a> Nature Methods. 9, 357-359 (2012).</li>
<li>Kearse, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Geneious+Basic%3A+an+integrated+and+extendable+desktop+software+platform+for+the+organization+and+analysis+of+sequence+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data.</a> Bioinformatics. 28 (12), 1647-1649 (2012).</li>
<li>Peirson, S. N., Butler, J. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Quantitative+polymerase+chain+reaction%5BTitle%5D)+AND+%22Methods+in+Molecular+Biology%22%5BJournal%5D)">Quantitative polymerase chain reaction.</a> Methods in Molecular Biology. 362, 349-362 (2007).</li>
<li>Hunt, M., Kikuchi, T., Sanders, M., Newbold, C., Berriman, M., Otto, T. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((REAPR+citation%3A+REAPR%3A+a+universal+tool+for+genome+assembly+evaluation%5BTitle%5D)+AND+%22Genome+Biology%22%5BJournal%5D)">REAPR citation: REAPR: a universal tool for genome assembly evaluation.</a> Genome Biology. 14 (5), (2013).</li>
<li>Meyer, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+mitochondrial+genome+sequence+of+a+hominin+from+Sima+de+los+Huesos%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">A mitochondrial genome sequence of a hominin from Sima de los Huesos.</a> Nature. 505 (7483), 403-406 (2013).</li>
<li>Gunnarsdóttir, E. D., Li, M., Bauchet, M., Finstermeier, K., Stoneking, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-throughput+sequencing+of+complete+human+mtDNA+genomes+from+the+Philippines%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">High-throughput sequencing of complete human mtDNA genomes from the Philippines.</a> Genome Research. 21 (1), 1-11 (2010).</li>
<li>King, J. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-quality+and+high-throughput+massively+parallel+sequencing+of+the+human+mitochondrial+genome+using+the+Illumina+MiSeq%5BTitle%5D)+AND+%22Forensic+Science+International%3A+Genetics%22%5BJournal%5D)">High-quality and high-throughput massively parallel sequencing of the human mitochondrial genome using the Illumina MiSeq.</a> Forensic Science International: Genetics. 12, 128-135 (2014).</li>
<li>Yao, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+First+Complete+Chloroplast+Genome+Sequences+in+Actinidiaceae%3A+Genome+Structure+and+Comparative+Analysis%5BTitle%5D)+AND+%22Plos+One%22%5BJournal%5D)">The First Complete Chloroplast Genome Sequences in Actinidiaceae: Genome Structure and Comparative Analysis.</a> Plos One. 10 (6), (2015).</li>
<li>Zhang, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Complete+Chloroplast+Genome+Sequences+of+Five+Epimedium+Species%3A+Lights+into+Phylogenetic+and+Taxonomic+Analyses%5BTitle%5D)+AND+%22Frontiers+in+Plant+Science%22%5BJournal%5D)">The Complete Chloroplast Genome Sequences of Five Epimedium Species: Lights into Phylogenetic and Taxonomic Analyses.</a> Frontiers in Plant Science. 7, (2016).</li>
<li>Swart, E. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Oxytricha+trifallax+Mitochondrial+Genome%5BTitle%5D)+AND+%22Genome+Biologyogy+and+Evolution%22%5BJournal%5D)">The Oxytricha trifallax Mitochondrial Genome.</a> Genome Biologyogy and Evolution. 4 (2), 136-154 (2011).</li>
<li>Barth, D., Berendonk, T. U. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+mitochondrial+genome+sequence+of+the+ciliate+Paramecium+caudatum+reveals+a+shift+in+nucleotide+composition+and+codon+usage+within+the+genus+Paramecium%5BTitle%5D)+AND+%22BMC+Genomics%22%5BJournal%5D)">The mitochondrial genome sequence of the ciliate Paramecium caudatum reveals a shift in nucleotide composition and codon usage within the genus Paramecium.</a> BMC Genomics. 12 (1), (2011).</li>
<li>Coombe, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Assembly+of+the+Complete+Sitka+Spruce+Chloroplast+Genome+Using+10X+Genomics%27+GemCode+Sequencing+Data%5BTitle%5D)+AND+%22Plos+One%22%5BJournal%5D)">Assembly of the Complete Sitka Spruce Chloroplast Genome Using 10X Genomics' GemCode Sequencing Data.</a> Plos One. 11 (9), (2016).</li>
<li>Herschleb, J., Ananiev, G., Schwartz, D. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pulsed-field+gel+electrophoresis%5BTitle%5D)+AND+%22Nature+Protocols%22%5BJournal%5D)">Pulsed-field gel electrophoresis.</a> Nature Protocols. 2 (3), 677-684 (2007).</li>
</ol>

Autorzy nie mają nic do ujawnienia.

Chociaż genomika subtraktywna jest potężna, nie jest to podejście sztampowe, wymagające dostosowania w kilku kluczowych krokach oraz starannego doboru sekwencji referencyjnych i próbek testowych. Jeśli zestaw zapytania jest niskiej jakości, kroki filtrowania mogą izolować tylko artefakty zestawu. Dlatego ważne jest, aby dokładnie zwalidować montaż de novo przy użyciu odpowiedniego protokołu walidacji do konkretnego projektu. W przypadku sekwencjonowania RNA wytyczne znajdują się na stronie internetowej Trinity<s...

Celem tej metody jest zidentyfikowanie nowej docelowej sekwencji genomowej (T), DNA lub RNA, na podstawie kontekstu genomowego, lub odniesienia (R) (<silna klasa="xfig">Rysunek 1). Metoda ta jest najbardziej przydatna, jeśli celu nie można fizycznie oddzielić lub byłoby to kosztowne. Tylko nieliczne organizmy mają idealnie ukończone genomy do odejmowania, więc kluczową innowacją naszej metody jest połączenie metod obliczeniowych i laboratoryjnych w cykl umożliwiający naukowcom wyizolowanie sekwencji docelowych, gdy odniesienie jest niedoskonałe, lub szkic genomu z organizmu niemodelowego. Pod koniec cyklu stosuje się test qPCR w celu określenia, czy potrzebne jest większe odjęcie. Zwalidowana kandydująca sekwencja T wykaże statystycznie większą wykrywalność w znanych próbkach T-dodatnich metodą qPCR.
Wcielenia metody zostały zaimplementowane w odkrywaniu nowych celów leków bakteryjnych, które nie mają homologów gospodarza1,2,3,4 oraz identyfikacja nowych wirusów z zainfekowanych hostów5,6. Oprócz identyfikacji T, metoda może poprawić R: niedawno zastosowaliśmy tę metodę do zidentyfikowania 936 brakujących genów w genomie referencyjnym zeberki oraz nowego genu z chromosomu (T) tylko w linii zarodkowej (T)7. Genomika subtraktywna jest szczególnie cenna, gdy istnieje prawdopodobieństwo, że T jest skrajnie rozbieżne ze znanymi sekwencjami lub gdy tożsamość T jest zasadniczo niezdefiniowana, jak w przypadku chromosomu ograniczonego do linii zarodkowej zięby zebry7.
Nie wymagając wcześniej pozytywnej identyfikacji T, kluczową zaletą genomiki subtraktywnej jest to, że jest bezstronna. W niedawnym badaniu Readhead i in. zbadali związek między chorobą Alzheimera a obfitością wirusów w czterech regionach mózgu. W celu identyfikacji wirusów Readhead i in. stworzyli bazę danych 515 wirusów8, poważnie ograniczając czynniki wirusowe, które ich badanie mogło zidentyfikować. Genomika subtraktywna mogła zostać wykorzystana do porównania genomu osób zdrowych i genomu osoby z chorobą Alzheimera w celu wyizolowania ewentualnych nowych wirusów związanych z chorobą, niezależnie od ich podobieństwa do znanych czynników zakaźnych. Chociaż znane są 263 wirusy atakujące ludzi, szacuje się, że istnieje około 1,67 miliona nieodkrytych gatunków wirusów, z których 631 000-827 000 może potencjalnie infekować ludzi9.
Izolacja nowych wirusów to obszar, w którym subtraktywna genomika jest szczególnie skuteczna, ale niektóre badania mogą nie wymagać tak rygorystycznej metody. Na przykład w badaniach identyfikujących nowe wirusy zastosowano bezstronne sekwencjonowanie o wysokiej przepustowości, a następnie odwrotną transkrypcję i BLASTx dla sekwencji wirusowych5 lub wzbogacanie wirusowych kwasów nukleinowych w celu ekstrakcji i odwrotnej transkrypcji sekwencji wirusowych6. Podczas gdy w badaniach tych stosowano sekwencjonowanie i składanie de novo, odejmowanie nie było stosowane, ponieważ sekwencje docelowe zostały pozytywnie zidentyfikowane przez BLAST. Gdyby wirusy były całkowicie nowe i nie były spokrewnione (lub daleko spokrewnione) z innymi wirusami, genomika subtraktywna byłaby użyteczną techniką. Zaletą genomiki subtraktywnej jest to, że można uzyskać sekwencje, które są zupełnie nowe. Jeśli genom organizmu jest znany, można go odjąć, aby pozostawić wszelkie sekwencje wirusa. Na przykład w naszym opublikowanym badaniu wyizolowaliśmy nową sekwencję wirusa od zeberki za pomocą genomiki subtraktywnej, chociaż nie było to naszym pierwotnym zamiarem7.
Genomika subtraktywna okazała się również przydatna w identyfikacji celów szczepionek bakteryjnych, motywowanych dramatycznym wzrostem oporności na antybiotyki1,2,3,4. Aby zminimalizować ryzyko wystąpienia reakcji autoimmunologicznej, naukowcy zawęzili potencjalne cele szczepionek, odejmując wszelkie białka, które mają homologi w ludzkim gospodarzu. W jednym szczególnym badaniu, dotyczącym Corynebacterium pseudotuberculosis, przeprowadzono odjęcie genomów gospodarza kręgowców od kilku genomów bakteryjnych, aby upewnić się, że ewentualne cele leków nie wpłyną na białka w gospodarzach, prowadząc do skutków ubocznych1. Podstawowym przepływem pracy w tych badaniach jest pobranie proteomu bakteryjnego, określenie ważnych białek, usunięcie zbędnych białek, użycie BLASTp do wyizolowania niezbędnych białek i BLASTp przeciwko proteomowi gospodarza w celu usunięcia wszelkich białek z homologami gospodarza1,2,3,4. W tym przypadku genomika subtraktywna zapewnia, że opracowane szczepionki nie będą miały żadnych skutków niepożądanych w host1,2,3,4.
Użyliśmy genomiki subtraktywnej, aby zidentyfikować pierwszy gen kodujący białko na chromosomie ograniczonym do linii zarodkowej (GRC) (w tym przypadku T), który znajduje się w liniach zarodkowych, ale nie w tkance somatycznej obu płci10. Przed tym badaniem jedyną informacją genomiczną, która była znana na temat GRC, był powtarzający się region11. Składanie de novo przeprowadzono na RNA zsekwencjonowanym z tkanek jajnika i jąder (R+T) od dorosłych zeberek. Obliczeniową eliminację sekwencji przeprowadzono przy użyciu opublikowanej sekwencji genomu somatycznego (mięśniowego) (R1)12, jej surowych (Sanger) odczytanych danych (R2) i transkryptomu somatycznego (mózgowego) (R3)13. Sekwencyjne użycie trzech odniesień było napędzane przez test qPCR w kroku 5 każdego cyklu (<silna klasa="xfig">Rysunek 2A), pokazując, że wymagane było dodatkowe filtrowanie. Odkryty gen α-SNAP został potwierdzony za pomocą qPCR z DNA i RNA oraz klonowania i sekwencjonowania. Pokazujemy w naszym przykładzie, że ta metoda jest elastyczna: nie jest zależna od dopasowania kwasów nukleinowych (DNA vs RNA) i że odejmowanie można przeprowadzić za pomocą odniesień (R), które składają się z zestawów lub surowych odczytów.

<table><tbody><tr><td>Accustart II Taq DNA Polymerase</td><td>Quanta Bio</td><td>95141</td><td></td></tr><tr><td>Blasic Local Alignment Search Tool (BLAST)</td><td></td><td></td><td>https://github.com/trinityrnaseq/trinityrnaseq/wiki/Transcriptome-Assembly-Quality-Assessment</td></tr><tr><td>Bowtie 2</td><td></td><td></td><td>https://www.python.org/download/releases/2.7/</td></tr> <tr><td>BWA-MEM v. 0.7.12</td><td></td><td></td><td>https://github.com/BenLangmead/bowtie2</td></tr><tr><td>Geneious</td><td></td><td></td><td>https://blast.ncbi.nlm.nih.gov/Blast.cgi</td></tr><tr><td>PEAR v. 0.9.6</td><td></td><td></td><td>http://www.mybiosoftware.com/reptile-1-1-short-read-error-correction.html</td></tr><tr><td>Komputer osobisty</td><td>Biomatters</td><td></td><td>http://www.geneious.com/</td></tr> <tr><td>PowerSYBR qPCR mix</td><td>ThermoFisher</td><td>4367659</td><td></td></tr><tr><td>Python v. 2.7</td><td></td><td></td><td>https://sco.h-its.org/exelixis/web/software/pear/</td></tr> <tr><td>Reptile v.1.1</td><td></td><td></td><td>https://alurulab.cc.gatech.edu/reptile</td></tr><tr><td>Stratagene Mx3005P</td><td>Agilent Technologies</td><td>401456</td><td></td></tr><tr><td>TransDecoder v. 3.0.1</td><td></td><td></td><td>https://sourceforge.net/projects/bio-bwa/files/</td></tr> <tr><td>Trinity v. 2.4.0</td><td></td><td></td><td>https://github.com/TransDecoder/TransDecoder/wiki</td></tr></tbody></table>

odkrycie nowych sekwencji za pomocą genomiki subtraktywnej

Novel Sequence Discovery by Subtractive Genomics

Subtractive genomika może być wykorzystana w każdym badaniu, którego celem jest zidentyfikowanie sekwencji genu, białka lub ogólnego regionu, który jest osadzony w szerszym kontekście genomowym. Genomika subtraktywna umożliwia badaczowi wyizolowanie docelowej sekwencji będącej przedmiotem zainteresowania (T) poprzez kompleksowe sekwencjonowanie i odejmowanie znanych elementów genetycznych (odniesienie, R). Metoda może być stosowana do identyfikacji nowych sekwencji, takich jak mitochondria, chloroplasty, wirusy lub chromosomy z ograniczeniami linii zarodkowej, i jest szczególnie przydatna, gdy T nie można łatwo wyizolować z R. Począwszy od kompleksowych danych genomowych (R + T), metoda wykorzystuje narzędzie Basic Local Alignment Search Tool (BLAST) względem sekwencji referencyjnej lub sekwencji, w celu usunięcia pasujących znanych sekwencji (R), pozostawiając za sobą cel (T). Aby odejmowanie działało najlepiej, R powinno być względnie kompletnym szkicem, w którym brakuje T. Ponieważ sekwencje pozostałe po odjęciu są testowane za pomocą ilościowej reakcji łańcuchowej polimerazy (qPCR), R nie musi być kompletne, aby metoda zadziałała. W tym miejscu łączymy kroki obliczeniowe z krokami eksperymentalnymi w cykl, który można iterować w razie potrzeby, sekwencyjnie usuwając wiele sekwencji referencyjnych i udoskonalając wyszukiwanie T. Zaletą genomiki subtraktywnej jest to, że całkowicie nową sekwencję docelową można zidentyfikować nawet w przypadkach, w których fizyczne oczyszczanie jest trudne, niemożliwe lub kosztowne. Wadą metody jest znalezienie odpowiedniego odniesienia do odejmowania i otrzymywanie próbek T-dodatnich i ujemnych do badania qPCR. Opisujemy nasze zastosowanie metody w identyfikacji pierwszego genu z chromosomu ograniczonego do linii zarodkowej zeberki. W tym przypadku filtrowanie obliczeniowe obejmowało trzy odniesienia (R), sekwencyjnie usuwane w ciągu trzech cykli: niekompletny zestaw genomowy, surowe dane genomowe i dane transkryptomiczne.

Po uruchomieniu BLAST, plik wynikowy będzie zawierał listę sekwencji z zapytania, które pasują do bazy danych. Po odjęciu sekwencji za pomocą Pythona uzyskanych zostanie szereg sekwencji niedopasowanych, które zostaną przetestowane metodą qPCR. Wyniki tego etapu oraz kolejne kroki omówiono poniżej.
Wynik negatywny. Po przeprowadzeniu analizy BLAST względem sekwencji referencyjnej mogą wystąpić dwa możliwe wyniki...

Celem tego protokołu jest wykorzystanie kombinacji badań obliczeniowych i laboratoryjnych do znalezienia nowych sekwencji, które nie mogą być łatwo oddzielone od sekwencji współoczyszczającej, która może być tylko częściowo znana.

Odkrycie nowych sekwencji za pomocą genomiki subtraktywnej

Subtractive genomics can be used in any research where the goal is to identify the sequence of a gene, protein, or general region that is embedded in a ...

Nasz protokół jest istotny, ponieważ może być wykorzystany do identyfikacji sekwencji genomowych, których nie można wyizolować z sekwencji współoczyszczających, które same w sobie mogą być tylko częściowo znane. Główną zaletą tej techniki jest to, że jest niedroga, wykorzystuje głównie darmowe oprogramowanie, które można pobrać, a także jest elastyczna. Można go zastosować do wielu zagadnień biologicznych.Potencjalne zastosowania obejmują identyfikację celów szczepionek bakteryjnych i identyfikację wirusów, których sekwencje znacznie różnią się od znanych drobnoustrojów. Metoda ta może być zastosowana do każdego systemu, w którym nieznane nie może być eksperymentalnie oddzielone, o ile dostępna jest sekwencja referencyjna bez celu genomowego. Ta technika wymaga trochę prób i błędów, więc ważna jest cierpliwość.Może być konieczne rozwiązywanie problemów z programami. Możesz używać programów innych niż te, które tutaj opisujemy. W miarę możliwości korzystaj z instrukcji obsługi.Wizualna demonstracja tej metody jest pomocna, ponieważ praca obliczeniowa zależy od podstawowej wiedzy na temat struktury programowania wiersza poleceń. Na początek użyj Trimmomatic 0.32 do usunięcia adapterów illumina i podstaw niskiej jakości. Użyj Pear w wersji 0.9.11, aby utworzyć wysokiej jakości scalone odczyty ze sparowanych odczytów wyjściowych trimmomatic przy użyciu parametrów domyślnych.Następnie użyj Reptile w wersji 1.1, aby poprawić błędy odczytów generowanych przez Pear. Na koniec użyj Trinity w wersji 2.4.0 w trybie domyślnym, aby złożyć poprawione sekwencje. W przypadku bibliotek specyficznych dla pasma należy użyć parametru SS_lib_type .Wynikiem jest plik FASA, który zostanie umieszczony w nowym katalogu o nazwie trinity_output. Utwórz bazę danych BLAST sekwencji referencyjnej nucleotide_reference. fasta w wierszu poleceń.BLAST dopasował zapytanie złożone do referencyjnej bazy danych. Aby uzyskać plik wyjściowy, użyj BLAST_results. txt Aby wygenerować tabelaryczne dane wyjściowe wymagane do etapów przetwarzania podsekwencji za pomocą skryptów Pythona, użyj outfmt 6 Aby zwiększyć rygorystyczność, użyj sekwencji białek z zestawu jako zapytania BLAST z przetłumaczonym nukleotydem BLAST, który wykonuje sześciokierunkową translację bazy danych nukleotydów.Aby uzyskać sekwencje białek dla zapytania, uruchom TransDecoder. Long0rfs do identyfikowania najdłuższych otwartych ramek odczytu ze złożonych sekwencji zapytań. Teraz uruchom tblastn.W razie potrzeby upewnij się, że wybrano prawidłowy kod genetyczny dla badanego organizmu za pomocą opcji db_gencode z odpowiednim kodem. Jeśli dostępne jest białko referencyjne wysokiej jakości, należy użyć dopasowania białko-białko za pomocą blastp, a nie tblastn Utwórz bazę danych BLAST odniesienia białka. Upewnij się, że wynik został zapisany jako plik do dalszego przetwarzania i użyj danych wyjściowych tabelarycznych, aby upewnić się, że skrypty języka Python mogą je poprawnie przeanalizować.Teraz użyj subtraktywnego skryptu Pythona, aby usunąć wszelkie pasujące sekwencje. Aby zamapować odczyty na zestaw, użyj BWA-MEM w wersji 0.7.12 lub bowtie 2 w celu zamapowania pobranych nieprzetworzonych odczytów na zestaw zapytań. Najpierw zindeksuj zestaw, a następnie zamapuj odczyty.Dane wyjściowe będą w formacie SAM. Uruchom skrypt języka Python removeUnmapped. py używając pliku SAM jako danych wejściowych.Identyfikuje nazwy sekwencji zapytań bez żadnych pasujących odczytów i zapisuje je w nowym pliku tekstowym. Danymi wyjściowymi poprzedniego kroku jest lista nazw sekwencji w pliku txt. Wyodrębnij plik FASTA z tymi sekwencjami.Wyjściem będzie plik fasta. Użyj programu Genius, aby ręcznie określić optymalne sekwencje starterów. Podkreśl sekwencję kandydującą od 21 do 28 par zasad dla startera do przodu, unikając przebiegów czterech lub więcej dowolnej zasady.Spróbuj celować w region z dość jednolitą kombinacją wszystkich par zasad. Pojedyncze G lub C na trzech głównych końcach jest korzystne, pomagając zakotwiczyć starter. Kliknij kartę Statystyki po prawej stronie ekranu, aby wyświetlić szacowaną temperaturę topnienia tej sekwencji, gdy zaznaczony jest region kandydujący.Dąż do temperatury topnienia między 55 a 60 stopni Celsjusza, unikając powtórzeń i długich serii G C. Wybierz odwrócony starter w ten sam sposób, umieszczony od 150 do 250 par zasad trzy prime przedniego startera. Chociaż długości podkładu nie muszą się zgadzać, przewidywana temperatura topnienia powinna być jak najbardziej zbliżona do temperatury startera przedniego. Pamiętaj, aby odwrócić uzupełnianie sekwencji, klikając prawym przyciskiem myszy w Genius, gdy sekwencja jest podświetlona w opcji menu.Alternatywną metodą jest użycie funkcji Primer Design, która znajduje się na górnym pasku narzędzi w oknie Sekwencja. Wstaw region do wzmocnienia w obszarze Region docelowy. Na karcie charakterystyki wprowadź żądany rozmiar, temperaturę topnienia i procent G C, a następnie kliknij przycisk OK, aby wygenerować startery.Aby przeprowadzić ilościową walidację pozostałej sekwencji metodą PCR, należy najpierw przygotować mieszaninę reakcyjną dla każdej matrycy w trzech egzemplarzach, z naszą mieszanką SYBR Green Master Mix, starterami do przodu i do tyłu oraz wodą, o łącznej objętości 25 mikrolitrów. Uruchom program qPCR w oparciu o wcześniej zweryfikowaną temperaturę i czas przedłużenia. Końcowe krzywe denaturacji powinny być generowane co najmniej przy pierwszym zastosowaniu starterów w qPCR w celu walidacji amplifikacji pojedynczego produktu DNA.Zmierz zielone sygnały qPCR SYBR w stosunku do aktyny przez Ct We wszystkich przypadkach oblicz średnią i odchylenie standardowe dwa do Ct względem aktyny. Wykonaj elektroforezę żelu w punkcie końcowym, aby potwierdzić prawidłowe wykrycie rozmiaru produktu za pomocą qPCR. Tutaj uruchom 25 mikrolitrów produktu qPCR zmieszanego z pięcioma mikrolitrami barwnika 6X Glitterol na żelu agarozowym 2% TAE pod napięciem 200 woltów przez 20 minut.Jeśli qPCR wykaże, że nie zidentyfikowałeś sekwencji docelowych, powtórz cały cykl z nowym odniesieniem, które można uzyskać z internetowej bazy danych. Projekt subtraktywny w tym przypadku rozpoczął się od sekwencjonowania RNA z tkanki wyściełanej zarodkami samca i samicy dorosłego zeberki. Ostatecznie zidentyfikowano 935 genów somatycznych, które wcześniej nie były uwzględnione w adnotacji całego genomu.Po filtrowaniu obliczeniowym ilościowy PCR może dać ujemny wynik, w którym nie było różnicy w wykrywalności w tkankach ptaków. I odwrotnie, pozytywny wynik reprezentujący identyfikację prawdziwej sekwencji docelowej jest potwierdzony, gdy qPCR genomowego DNA wykazuje statystycznie większą wykrywalność w tkance będącej przedmiotem zainteresowania w stosunku do wartości referencyjnej. W tym przypadku potwierdzono, że gen alfa snap jest ograniczony do linii zarodkowej, ponieważ był zubożony w tkance somatycznej w stosunku do DNA jądra, w którym występował na poziomie równoważnym aktynie.Ważne jest, aby pamiętać o użyciu właściwych danych wejściowych podczas każdego z kroków obliczeniowych. Być może będziesz musiał wielokrotnie przejść przez odejmowanie cyklu, aby uzyskać docelową sekwencję lub sekwencje. Odkryte geny można przeprowadzić różnorodną analizę filonetyczną, strukturalną i funkcjonalną.Te dodatkowe metody dają wgląd w ewolucyjną i funkcjonalną rolę genów. Zidentyfikowaliśmy pierwszy gen na chromosomie ptaka śpiewającego z ograniczeniami linii zarodkowej, co poszerzyło zainteresowanie tym zaskakującym elementem genomu. Późniejsze prace wykazały podobne chromosomy u wielu gatunków ptaków śpiewających zawierających wiele dodatkowych genów.

novel-sequence-discovery-by-subtractive-genomics

Research

JoVE Journal

Genetics

Genetyka

Opracowaliśmy kompleksowe, nieobciążone systematycznie badania przesiewowe całego genomu, aby zrozumieć interakcje gen-lek oraz gen-środowisko. Przedstawiono metody przesiewania tych kolekcji mutantów.

Konkurencyjne przesiewy genomiczne z zastosowaniem bibliotek drożdży z kodami kreskowymi

Sekwencjonowanie puli DNA jest szybką i kosztowo efektywną strategią wykrywania rzadkich wariantów związanych ze złożonymi fenotypami w dużych kohortach. W niniejszym materiale opisujemy analizę obliczeniową sekwencjonowania nowej generacji pul 32 genów związanych z nowotworami przy użyciu pakietu oprogramowania SPLINTER. Metoda ta jest skalowalna i może być zastosowana do dowolnego analizowanego fenotypu.

Wykrywanie rzadkich wariantów genomicznych z sekwencjonowania pulującego przy użyciu SPLINTER

Scharakteryzowanie społeczności mikrobiologicznej było od dawna celem w mikrobiologii środowiskowej. Metody sekwencjonowania nowej generacji pozwalają obecnie na scharakteryzowanie społeczności drobnoustrojów na niespotykanej dotąd głębokości przy minimalnych kosztach i nakładzie pracy. Szczegółowo opisujemy tutaj nasze podejście do sekwencjonowania bakteryjnych genów rybosomalnego RNA 16S za pomocą sekwenatora laboratoryjnego.

Sekwencjonowanie nowej generacji amplikonów genu rybosomalnego RNA 16S

Sekwencjonowanie amplikonów bisulfitowych (BSAS) to metoda ilościowego określania metylacji cytozyny w docelowych regionach genomu. Metoda ta wykorzystuje konwersję wodorosiarczynów w połączeniu z amplifikacją PCR regionów docelowych przed sekwencjonowaniem nowej generacji w celu uzyskania bezwzględnej oceny ilościowej metylacji DNA na poziomie specyficznym dla zasady.

Ukierunkowana analiza metylacji DNA za pomocą sekwencjonowania nowej generacji

Opisujemy technikę łączenia cytometrii przepływowej i sekwencjonowania o wysokiej przepustowości w celu identyfikacji późno replikujących się regionów genomu.

G2-seq: Wysokoprzepustowa technika oparta na sekwencjonowaniu do identyfikacji późno replikujących się regionów genomu

Sekwencjonowanie nowej generacji (NGS) to potężne narzędzie do charakterystyki genomu, które jest ograniczone przez wysoki poziom błędów platformy (~0,5–2,0%). Opisujemy nasze metody sekwencjonowania z korekcją błędów, które pozwalają nam wyeliminować wskaźnik błędów NGS i wykryć mutacje przy frakcjach alleli wariantów tak rzadkich, jak 0,0001.

Wykrywanie rzadkich zdarzeń za pomocą sekwencjonowania DNA i RNA z korekcją błędów

Identyfikacja fizycznych interakcji między genami a elementami regulatorowymi jest trudna, ale została ułatwiona dzięki metodom wychwytywania konformacji chromosomów. Ta modyfikacja protokołu 4C-seq łagodzi stronniczość PCR, minimalizując nadmierną amplifikację matryc PCR i maksymalizując mapowalność odczytów poprzez włączenie etapu trawienia enzymu ograniczającego dodawanie.

Wysokoprzepustowa identyfikacja sekwencji regulatorowych genów przy użyciu sekwencjonowania nowej generacji wychwytywania konformacji chromosomów kołowych (4C-seq)

Protokół eksploracji genomu sterowany spektrometrią mas jest ustanowiony i opisany tutaj. Opiera się na informacjach o sekwencjach genomu i analizie LC-MS/MS i ma na celu ułatwienie identyfikacji cząsteczek ze złożonych ekstraktów mikrobiologicznych i roślinnych.

Eksploracja genomu pod kontrolą spektrometrii mas jako narzędzie do odkrywania nowych produktów naturalnych

Ten protokół opisuje proces od hodowli komórkowych ex vivo lub in vitro do wstępnego przetwarzania danych transkryptomicznych w celu opłacalnego badania przesiewowego leków opartego na transkryptomie.

Ekonomiczne badania przesiewowe leków oparte na transkryptomii

W tym manuskrypcie opisujemy protokół do funkcjonalnego badania transkrypcji i hamującej aktywności środków przeciwbakteryjnych ukierunkowanych na bakteryjną transkrypcję.

In vitro Testy transkrypcji i ich zastosowanie w odkrywaniu leków

Odkryj genomikę subtraktywną wykorzystującą darmowe oprogramowanie do identyfikacji nowych sekwencji genomowych trudnych do wyizolowania. Ta elastyczna metoda komputerowa znajduje zastosowanie w poszukiwaniu celów szczepionkowych bakterii i nieznanych wirusów.

Odkrycie nowych sekwencji za pomocą genomiki subtraktywnej

Przeglądaj więcej filmów