Name: Sequenciamento de nova geração de RNA e um fluxo de análise de bioinformática para identificar expressão de LINE-1s a nível locus-específico (Vídeo) | JoVE Journal
Uploaded: 2019-05-19T15:00:00.000+00:00
Duration: 11 min 4 s
Description: Aprenda um fluxo de bioinformática para identificar expressão locus-específica de LINE-1s usando RNA-Seq e ferramentas como Bowtie1, SAMtools e IGV. Este método rigoroso elimina ruído de fundo e permite quantificação confiável em células tumorais.

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

Gostaríamos de agradecer ao Dr. Yan Dong pelas células tumorais da próstata DU145. Gostaríamos de agradecer ao Dr. Nathan Ungerleider por sua orientação e aconselhamento na criação de scripts de supercomputador. Parte deste trabalho foi financiado pela NIH Grants R01 GM121812 para PD, R01 AG057597 para VPB, e 5TL1TR001418 para TK. Também gostaríamos de reconhecer o apoio do cancer Crusaders e do centro de bioinformática do Tulane Cancer Center.

1. extração de RNA citoplasmático
<ol><li>Obter células através dos seguintes métodos.<ol>
<li>Colete células ao vivo de 2,75% – 100% confluentes, frascos T-75.<ol>
<li>Lave o balão 2 vezes em 5 mL de PBS frio e, na última lavagem, raspe as células e transfira para um tubo cônico de 15 mL. Centrifugue durante 2 min a 1.000 x g e 4 ° c e retire e elimine cuidadosamente o sobrenadante (tabela de materiais).</li>
		</ol>
</li>
		<li>Colete células de espécimes de tecido.<ol>
<li>Prepare o tecido para extração de RNA citoplasmático dentro de uma hora de ser dissecado e sempre manter no gelo. Para o armazenamento a longo prazo, use soluções do inibidor do RNA para armazenar o tecido por até 72 horas após a dissecção que segue o protocolo do fabricante (tabela dos materiais).</li>
			<li>Dados uma amostra de 10 μm3 e homogeneizar a amostra fresca com 5 ml de PBS frio em um homogeneizador de com estéril, transferência para um tubo cônico de 15 ml, centrífuga por 2 min a 1.000 x g a 4 ° c, e remover e descartar cuidadosamente o sobrenadante (tabela de materiais < /C8 >).</li>
		</ol>
</li>
	</ol>
</li>
	<li>Adicione 2 mL de tampão de Lise ao pellet-Mix celular e incubar no gelo por 5 min.<ol>
<li>Prepare a reserva de Lise fresca com 150 mM de NaCl, 50 mM de HEPES (pH 7,4) e 25 μg/mL de digitonina (tabela de materiais).</li>
		<li>Como a concentração mínima de digitonina no tampão de Lise necessária para penetrar na membrana plasmática pode variar de acordo com o tipo de célula, confirme microscopicamente que as células tratadas com tampão de Lise perdem a membrana plasmática e mantêm a membrana nuclear intacta.</li>
		<li>Pouco antes de usar, adicionar 1.000 U/mL inibidor de RNase (tabela de materiais).</li>
	</ol>
</li>
	<li>Centrifugue por 1 minuto em 1.000 x g e 4 ° c, e colete o sobrenadante.</li>
	<li>Adicionar sobrenadante a pré-refrigerados 7,5 mL de Trizol e 1,5 mL de clorofórmio. Todos os passos que necessitam de clorofórmio devem ser feitos dentro de uma capa química limpa (tabela de materiais).</li>
	<li>Centrifugador para 35 min em 3.220 x g e 4 ° c.</li>
	<li>Transfira a porção aquosa (camada superior) para um tubo de 15 mL pré-refrigerado fresco.</li>
	<li>Adicionar 4,5 mL de clorofórmio e Vortex.</li>
	<li>Centrifugador por 10 min a 3.220 x g e 4 ° c.</li>
	<li>Transfira a porção aquosa para o tubo pré-refrigerado fresco.</li>
	<li>Adicionar 4,5 mL de isopropanol, agitar bem e incubar a-80 ° c durante a noite (tabela de materiais).</li>
	<li>Centrifugador em 3.220 x g e 4 ° c por 45 minutos.</li>
	<li>Remover isopropanol, adicione 15 mL de 100% de etanol (tabela de materiais).</li>
	<li>Centrifugador a 3.220 x g durante 10 min.</li>
	<li>Remova o etanol, escorra e seque por aproximadamente 1 h.<ol>
<li>Use um cotonete de algodão estéril para apagar qualquer etanol restante (tabela de materiais).</li>
	</ol>
</li>
	<li>Resuspenda a amostra em 100 a 200 μL de água livre de RNase dependendo do tamanhoda pelota (tabela de materiais).</li>
	<li>Fractionate amostras usando a tecnologia da electroforese para determinar a qualidade e a concentração das amostras de acordo com intructions do fabricante23 (tabela dos materiais).<ol>
<li>As amostras se qualificam para análise de RNA-Seq se RIN > 824.</li>
	</ol>
</li>
</ol>2. sequenciamento da próxima geração
<ol><li>Envie amostras de RNA citoplasmáticos a serem sequenciadas usando a plataforma de sequenciamento de próxima geração destinada a gerar pelo menos 50 milhões leituras de 100 BP de extremidade emparelhada.</li>
	<li>Selecione para RNAs poli-adeniladas e sequenciamento específico da vertente.</li>
</ol>3. criar anotações (opcional se um tiver uma anotação existente)
<ol><li>Crie uma anotação L1 de comprimento total ou faça o download da anotação L1 de comprimento total (arquivo suplementar 1a-b).<ol>
<li>Transfira anotações do Masker da repetição para os elementos LINE-1 do navegador do genoma de UCSC com a ferramenta do navegador da tabela (https://genome.ucsc.edu/cgi-bin/hgTables). Especifique o clado de mamíferos, o genoma humano, o conjunto hg19 (ou hg38 para um genoma mais atualizado) e filtre para "LINE1" em nome da classe. Baixe como um arquivo. GTF e Rotule como FL-L1-BLAST. GTF.</li>
		<li>Execute uma pesquisa BLAST local da primeira 300 BP do elemento L1 de comprimento total L 1.3 abrangendo a região promotora no genoma humano e adicione 6.000 BP downstream para criar um fim das coordenadas L1 para o arquivo de anotação. Salve em um arquivo GTF e Rotule como FL-L1-RM. GTF.</li>
		<li>Intercepta a anotação RepeatMasker e a anotação L1 baseada em promotor usando bedtools e Rotule como FL-L1-BLAST_RM. txt (pacotes de software).<ol>
<li>Use este comando no terminal Linux: bedtools Intersect-a fl-L1-Blast. GTF-b FL-L1-RM. gtf > FL-L1-BLAST_RM. txt.
</li>
		</ol>
</li>
		<li>Separe a anotação FL-L1 intersectada pela vertente superior e inferior.<ol>
<li>Copie sobre o FL-L1-BLAST_RM. txt em software de planilha e classificar pela "menos" e "mais" vertente e, em seguida, ordenar por localização do cromossomo.</li>
			<li>Crie dois novos documentos de planilha, um com as coordenadas intersectadas para o comprimento total L1s na vertente menos e uma na vertente inferior e salve como FL-L1-BLAST_RM_minus. xls e FL-L1-BLAST_RM_plus. xls.</li>
			<li>Guarde os dois novos documentos como ficheiros. txt.</li>
		</ol>
</li>
		<li>Use o programa mac2unix para converter os arquivos. txt para os arquivos de anotação corretos (pacotes de software).<ol>
<li>Use este comando no terminal: Mac2unix.sh FL-L1-BLAST_RM_minus. GFF.
</li>
			<li>Use este comando no terminal: Mac2unix.sh FL-L1-BLAST_RM_plus. GFF.
</li>
			<li>Salve novos arquivos com a extensão. GFF.</li>
		</ol>
</li>
		<li>Alternativamente, use AWK para filtrar linhas associadas com a vertente + e –.<ol>
<li>Use o seguinte comando para obter o + vertente: awk '/+/' FL-L1_BLAST_RM. gtf > FL-L1_BLAST_RM_plus. GTF.</li>
			<li>Use a seguinte linha de comando para obter o-Strand: awk '/-/' FL-L1_BLAST_RM. gtf > FL-L1_BLAST_RM_minus. GTF.</li>
		</ol>
</li>
	</ol>
</li>
</ol>4. Leia o pipeline de alinhamento para identificar as L1s expressas
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always"><tbody>
<tr>
<td>Opção</td>
			<td>Descrição</td>
		</tr>
<tr>
<td>– p</td>
			<td>Isso detalha o número de threads que o computador deve usar executando o alinhamento. Maior memória do computador permitirá mais threads e deve ser empiricamente d.</td>
		</tr>
<tr>
<td>– m 1</td>
			<td>Isso informa ao programa para aceitar apenas leituras que têm uma correspondência no genoma que é melhor do que qualquer outra correspondência de genoma.</td>
		</tr>
<tr>
<td>– y</td>
			<td>Este é o interruptor de tryhard que faz a busca do mapeamento para todas as correspondências possíveis e para não permitir que feche depois que um número fixo de fósforos é alcangado.</td>
		</tr>
<tr>
<td>– v 3</td>
			<td>Isso só permite que o programa utilize memória para leituras mapeadas com 3 ou menos correspondências para o genoma.</td>
		</tr>
<tr>
<td>– X 600</td>
			<td>Isso permite somente leituras emparelhadas que mapeiam dentro de 600 bases umas das outras. Isto certifica-se que os pares lidos são co-lineares no genoma e seleciona de encontro a s que envolvem moléculas processadas do RNA.</td>
		</tr>
<tr>
<td>– chunkmbs 8184</td>
			<td>Esse comando atribui memória extra para manipular a grande quantidade de alinhamentos possíveis para cada leitura relacionada a L1.</td>
		</tr>
</tbody></table>Tabela 1: opções de linha de comando para o bowtie.
<ol><li>Execute os arquivos fastq de sequenciamento de fim emparelhado de alinhamento com a amostra de RNA-Seq de interesse usando bowtie. Nota: Bowtie1 deve ser usado e não Bowtie2 porque os parâmetros necessários para o alinhamento exclusivo são especificamente encontrados apenas nesta versão do bowtie (pacotes de software). Bowtie é usado sobre alinhadores Splice-Aware como STAR em ordem avaliar concordante, contíguo lê mais relevante para a biologia L1 e expressão.<ol>
<li>Use esta linha de comando no terminal Linux: bowtie-p 10-m 1-S-y-v 3-X 600--chunkmbs 8184 hg_X_Y_M_index-1 hg_sample_1. FQ-2 hg_sample_2. FQ | samtools View-hbuS-| samtools Sort – hg_sample_sorted. bam. Consulte a tabela 1 para obter uma descrição das opções de linha de comando para o bowtie.</li>
	</ol>
</li>
	<li>Vertente separar o arquivo de saída Bam usando samtools (pacotes de software) e os seguintes comandos do Linux. Observe que os valores de sinalizador real podem variar se um não estiver usando protocolos de seqüenciamento de próxima geração padrão.<ol>
<li>Use esta linha de comando para selecionar para a vertente superior: samtools View-h hg_sample_sorted. bam | awk ' substr ($ 0, 1, 1) = = "@" | | $2 = = 83 | | $2 = = 163 {Print} ' | samtools View-bS-> hg_sample_sorted_topstrand. bam.</li>
		<li>Use esta linha de comando para selecionar para a vertente inferior: samtools View-h hg_sample_sorted. bam | awk ' substr ($ 0, 1, 1) = = "@" | | $2 = = 99 | | $2 = = 147 {Print} ' | samtools View-bS-> hg_sample_sorted_bottomstrand. bam.</li>
	</ol>
</li>
	<li>Gere contagens de leitura contra anotações para L1 loci usando bedtools (pacotes de software).<ol>
<li>Use esta linha de comando para gerar contagens de leitura para L1s na direção do sentido na vertente superior: bedtools Coverage-Abam FL-L1-BLAST_RM_plus. GFF-b hg_sample_sorted_topstrand. bam > hg_sample_sorted_bowtie_tryhard_plus_top. txt.</li>
		<li>Use esta linha de comando para gerar contagens de leitura para L1s na direção do sentido na vertente inferior: bedtools Coverage-Abam FL-L1-BLAST_RM_minus. GFF-b hg_sample_sorted_bottomstrand. bam > hg_sample_sorted_bowtie_tryhard_minus_bottom. txt.</li>
	</ol>
</li>
	<li>Indexar o arquivo Bam da etapa 5.1.1 para torná-lo visível no Visualizador de genômica Integrativa (IGV)25 (pacotes de software).<ol>
<li>Use esta linha de comando: samtools index hg_sample_sorted. bam
</li>
	</ol>
</li>
	<li>Para usar um modo de lote para aumentar o número de amostras de RNA-Seq canalizada por um tempo, use um script de supercomputador para concluir a etapa 4,1 chamada human_bowtie. sh, um script para concluir as etapas 4.2-4.3 foi criado chamado human_L1_pipeline. sh e um script para concluir etapa 4,4 foi criado chamado bam_index. sh. Esses scripts podem ser encontrados no arquivo suplementar 2 com comandos supercomputador associados para executar os scripts.</li>
</ol>5. Curação manual
<ol><li>Crie uma planilha para leituras mapeadas para cada Locus L1 anotado.<ol>
<li>Cópia sobre hg_sample_sorted_bowtie_tryhard_minus_bottom. txt criado na etapa 4.3.2 e página de rótulo como "menos-inferior."<ol>
<li>Classifique todas as colunas com base no número mais alto para o menor de leituras encontradas na coluna J.</li>
		</ol>
</li>
		<li>Copie hg_sample_sorted_bowtie_tryhard_plus_top. txt criado na etapa 4.3.1 e Rotule como "Top-Plus" em outra planilha.<ol>
<li>Classifique todas as colunas com base no número mais alto para o menor de leituras encontradas na coluna J.</li>
		</ol>
</li>
		<li>Crie uma terceira página rotulada como "combinada" e adicione todos os loci com dez ou mais leituras das páginas "menos-inferior" e "Plus-Top".<ol>
<li>Classifique todas as colunas com base no número mais alto para o menor de leituras encontradas na coluna J.</li>
		</ol>
</li>
		<li>Carregue os seguintes arquivos no IGV25 (pacotes de software): 1) genoma de referência de interesse para visualizar genes anotados, 2) FL-L1-BLAST_RM. GFF para visualizar a anotação L1, 3) hg_sample_sorted. bam para visualizar transcrições mapeadas de amostra de interesse, e 4) hg_genomicDNA_sorted. bam para avaliar a mappability de regiões genomic.</li>
		<li>Remova a cobertura e as linhas de junção associadas a cada arquivo Bam.</li>
		<li>Comprimir hg_sample_sorted. Bam e hg_genomicDNA_sorted. bam para que todas as faixas IGV caber em uma tela.</li>
	</ol>
</li>
	<li>Manualmente curate.<ol>
<li>Usando as coordenadas de loci listado na planilha "combinado" página, ver chamado loci em IGV25 (pacotes de software).</li>
		<li>Curate um locus a ser expressado autenticamente fora de seu próprio se não há nenhuma leitura ascendente na direção L1 até 5 KB.<ol>
<li>Rotule a linha verde na cor e anote porque é um L1 autenticamente expressado. Nota: uma excepção a esta regra existe se a região a montante do L1 não é mappable. Se for esse o caso, rotule a linha de cor vermelha e observe que a expressão da região a montante do promotor L1 não pode ser avaliada e, portanto, a expressão L1's não é capaz de ser determinada com confiança.</li>
		</ol>
</li>
		<li>Curate um locus para não ser autenticamente expressa fora de seu próprio promotor, se houver leituras upstream até 5 KB.<ol>
<li>Rotule a linha de cor vermelha e observe por que não é um L1 autenticamente expresso.</li>
			<li>Curate um locus como o falso se é expressado dentro de um intron de um gene expresso no mesmo sentido com leituras a montante do L1, se é a jusante de um gene expresso no mesmo sentido com leituras a montante do L1, ou para os testes padrões un-anotated da expressão com re anúncios a montante do L1. Observação: uma exceção a essa regra se aplica quando há leituras mínimas diretamente sobrepondo o local de início do promotor L1, mas ligeiramente a montante do L1. Se não há nenhum outro lê a montante de um caso L1 como este, considere este L1 para ser expressado autenticamente. Rotule a cor verde da linha e observe por que ele é um L1 autenticamente expresso.</li>
		</ol>
</li>
		<li>Curate um locus L1 como provável ser falso se o teste padrão de leituras traçadas ao Locus não se correlacionar com as regiões L1's específicas do mappability. Nota: por exemplo, se um L1 é altamente mapeável mas tem somente uma pilha acima das leituras em uma região condensada dentro do L1, é menos provável estar relacionada à expressão L1 fora de seu próprio promotor e mais provável ser das fontes un-anotated como exons ou ltrs. Em casos como este, cura o loci como a laranja e anote porque o locus é suspeito. Verifique as fontes de pilha suspeitos, verificando a localização L1 no navegador do genoma UCSC.</li>
		<li>Curate um locus para não ser autenticamente expressado se está dentro de um ambiente genomic de regiões un-anotated esporadicamente expressas Nota: por exemplo, leituras podem ser expressas 10 KB upstream do L1, mas cada 10 KB ou assim há leituras mapeadas e algumas dessas leituras alinhar com o L1. Estes L1s são menos propensos a ser expressa fora de seu próprio promotor, e mais propensos a ter mapeado leituras devido a un-anotado padrões de expressão genômica. Em casos como este, cura o loci como a laranja e anote porque o locus é suspeito.</li>
	</ol>
</li>
</ol>6. Leia a estratégia do alinhamento para avaliar o mappability no genoma da referência (opcional se um tem um conjunto de dados genomic alinhado existente do ADN)
<ol><li>Transfira arquivos inteiros da seqüência do ADN do genoma e converta-os aos arquivos. FQ<ol>
<li>Vá para o site NCBI encontrado aqui: https://www.ncbi.nlm.nih.gov/sra</li>
		<li>Digite WGS HeLa final emparelhado.</li>
		<li>Selecione para Homo sapiens resultados por táxon.</li>
		<li>Selecione um exemplo que está emparelhado final e tem leituras com 100 ou mais BP como o exemplo a seguir: https://www.ncbi.nlm.nih.gov/sra/ERX457838 [ACCN]</li>
		<li>Confirme o comprimento de leitura selecionando executar e, em seguida, metadados como mostrado aqui: https://Trace.NCBI.nlm.nih.gov/TRACES/Sra/?Run=ERR492384</li>
		<li>Para baixar todos os dados da sequência de DNA do genoma, digite este comando no terminal Linux: sratoolkit. 2.9.2-mac64/bin/prefetch-X 100g ERR492384 Nota: a função de pré-busca do Toolkit SRA transfere o número de adesão "ERR492384" encontrado no site NCBI (pacotes de software). O "100G" limita a quantidade de dados baixados para 100 gigabytes.</li>
		<li>Digite este comando no terminal Linux: fastq-dump--Split-files ERR492384 Observação: isso divide o conjunto de dados de DNA genômica baixado em dois arquivos fastq.</li>
	</ol>
</li>
	<li>Execute o alinhamento usando bowtie.<ol>
<li>Use este comando no Linux para alinhamento: bowtie-p 10-m 1-S-y-v 3-X 600--chunkmbs 8184 hg_X_Y_M_index-1 hg_genomicDNA_1. FQ-2 hg_genomicDNA_2. FQ | samtools View-hbuS-| samtools Sort – hg_genomicDNA_sorted. bam.<ol>
<li>Consulte a etapa 4,1 para entender os parâmetros usados no alinhamento do bowtie (pacotes de software).</li>
			<li>Transfira o arquivo Bam genômico alinhado para avaliar o mappability disponível em cima do pedido do autor.</li>
		</ol>
</li>
	</ol>
</li>
	<li>Indexar o arquivo Bam da etapa 4.2.1 usando samtools para torná-lo visível no IGV25 (pacotes de software) para informar ainda mais a Curação manual.<ol>
<li>Use esta linha de comando no Linux: samtools index hg_genomicDNA_sorted. bam
</li>
	</ol>
</li>
	<li>Avaliar a mappability de cada L1 loci<ol>
<li>Determine o número de leituras mapeadas exclusivamente para L1 loci usando o programa bedtools, a anotação FL-L1 e os dados de sequência genômica alinhados (pacotes de software).<ol>
<li>Use esta linha de comando no Linux: bedtools Coverage-Abam FL-L1-BLAST_RM. GTF – b hg_genomicDNA_sorted. bam &Gt; L1_Mappability_hg_genomicDNA. txt.</li>
		</ol>
</li>
		<li>Designe um locus L1 para ter o mappability cheio da cobertura quando 400 leituras originais são alinhadas a ela.</li>
		<li>Determine o fator exigido para escalar acima ou para baixo ADN genomic alinhado lê a 400 para cada L1 individual.</li>
		<li>Para ter uma medida dimensionada de expressão de acordo com a mappability individual do locus L1, multiplique o fator determinado na etapa 6.4.3 ao número de leituras do transcrito do RNA que alinham ao L1s autenticamente expressado determinado nas seções 4 – 5.</li>
	</ol>
</li>
</ol>

<ol><li>International Human Genome Sequencing. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Initial+sequencing+and+analysis+of+the+human+genome%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Initial sequencing and analysis of the human genome.</a> Nature. 409, 860(2001).</li>
<li>Brouha, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Hot+L1s+account+for+the+bulk+of+retrotransposition+in+the+human+population%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Hot L1s account for the bulk of retrotransposition in the human population.</a> Proceedings of the National Academy of Sciences of the United States of America. 100 (9), 5280-5285 (2003).</li>
<li>Dombroski, B. A., Mathias, S. L., Nanthakumar, E., Scott, A. F., Kazazian, H. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Isolation+of+an+active+human+transposable+element%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">Isolation of an active human transposable element.</a> Science. 254 (5039), 1805(1991).</li>
<li>Swergold, G. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification%2C+characterization%2C+and+cell+specificity+of+a+human+LINE-1+promoter%5BTitle%5D)+AND+%22Molecular+and+Cellular+Biology%22%5BJournal%5D)">Identification, characterization, and cell specificity of a human LINE-1 promoter.</a> Molecular and Cellular Biology. 10 (12), 6718-6729 (1990).</li>
<li>Speek, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Antisense+promoter+of+human+L1+retrotransposon+drives+transcription+of+adjacent+cellular+genes%5BTitle%5D)+AND+%22Molecular+and+Cellular+Biology%22%5BJournal%5D)">Antisense promoter of human L1 retrotransposon drives transcription of adjacent cellular genes.</a> Molecular and Cellular Biology. 21 (6), 1973-1985 (2001).</li>
<li>Deininger, L., Batzer, M. A., Hutchison, C. A., Edgell, M. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Master+genes+in+mammalian+repetitive+DNA+amplification%5BTitle%5D)+AND+%22Trends+in+Genetics%22%5BJournal%5D)">Master genes in mammalian repetitive DNA amplification.</a> Trends in Genetics. 8 (9), 307-311 (1992).</li>
<li>Boissinot, S., Chevret, P., Furano, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((L1+%28LINE-1%29+Retrotransposon+Evolution+and+Amplification+in+Recent+Human+History%5BTitle%5D)+AND+%22Molecular+Biology+and+Evolution%22%5BJournal%5D)">L1 (LINE-1) Retrotransposon Evolution and Amplification in Recent Human History.</a> Molecular Biology and Evolution. 17 (6), 915-918 (2000).</li>
<li>Khazina, E., Weichenrieder, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Non-LTR+retrotransposons+encode+noncanonical+RRM+domains+in+their+first+open+reading+frame%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Non-LTR retrotransposons encode noncanonical RRM domains in their first open reading frame.</a> Proceedings of the National Academy of Sciences of the United States of America. 106 (3), 731-736 (2009).</li>
<li>Martin, S. L., Bushman, F. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Nucleic+acid+chaperone+activity+of+the+ORF1+protein+from+the+mouse+LINE-1+retrotransposon%5BTitle%5D)+AND+%22Molecular+and+Cellular+Biology%22%5BJournal%5D)">Nucleic acid chaperone activity of the ORF1 protein from the mouse LINE-1 retrotransposon.</a> Molecular and Cellular Biology. 21 (2), 467-475 (2001).</li>
<li>Feng, Q., Moran, M. H., Kazazian, H. H., Boeke, J. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Human+L1+Retrotransposon+Encodes+a+Conserved+Endonuclease+Required+for+Retrotransposition%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Human L1 Retrotransposon Encodes a Conserved Endonuclease Required for Retrotransposition.</a> Cell. 87 (5), 905-916 (1996).</li>
<li>Mathias, S. L., Scott, A. F., Kazazian, H. H., Boeke, J. D., Gabriel, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reverse+transcriptase+encoded+by+a+human+transposable+element%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">Reverse transcriptase encoded by a human transposable element.</a> Science. 254 (5039), 1808(1991).</li>
<li>Luan, D. D., Korman, M. H., Jakubczak, J. L., Eickbush, T. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reverse+transcription+of+R2Bm+RNA+is+primed+by+a+nick+at+the+chromosomal+target+site%3A+A+mechanism+for+non-LTR+retrotransposition%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Reverse transcription of R2Bm RNA is primed by a nick at the chromosomal target site: A mechanism for non-LTR retrotransposition.</a> Cell. 72 (4), 595-605 (1993).</li>
<li>van den Hurk, J. A. J. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Novel+types+of+mutation+in+the+choroideremia+%28CHM%29+gene%3A+a+full-length+L1+insertion+and+an+intronic+mutation+activating+a+cryptic+exon%5BTitle%5D)+AND+%22Human+Genetics%22%5BJournal%5D)">Novel types of mutation in the choroideremia (CHM) gene: a full-length L1 insertion and an intronic mutation activating a cryptic exon.</a> Human Genetics. 113 (3), 268-275 (2003).</li>
<li>Miné, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+large+genomic+deletion+in+the+PDHX+gene+caused+by+the+retrotranspositional+insertion+of+a+full-length+LINE-1+element%5BTitle%5D)+AND+%22Human+Mutation%22%5BJournal%5D)">A large genomic deletion in the PDHX gene caused by the retrotranspositional insertion of a full-length LINE-1 element.</a> Human Mutation. 28 (2), 137-142 (2007).</li>
<li>Solyom, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pathogenic+orphan+transduction+created+by+a+nonreference+LINE-1+retrotransposon%5BTitle%5D)+AND+%22Human+Mutation%22%5BJournal%5D)">Pathogenic orphan transduction created by a nonreference LINE-1 retrotransposon.</a> Human Mutation. 33 (2), 369-371 (2012).</li>
<li>Hancks, D. C., Kazazian, H. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Roles+for+retrotransposon+insertions+in+human+disease.+Mobile+DNA%5BTitle%5D)+AND+%22Mobile+DNA%22%5BJournal%5D)">Roles for retrotransposon insertions in human disease. Mobile DNA.</a> Mobile DNA. 7, 9-9 (2016).</li>
<li>Tubio, J. M. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mobile+DNA+in+cancer.+Extensive+transduction+of+nonrepetitive+DNA+mediated+by+L1+retrotransposition+in+cancer+genomes%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">Mobile DNA in cancer. Extensive transduction of nonrepetitive DNA mediated by L1 retrotransposition in cancer genomes.</a> Science. 345 (6196), 1251343-1251343 (2014).</li>
<li>Ewing, A. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Widespread+somatic+L1+retrotransposition+occurs+early+during+gastrointestinal+cancer+evolution%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">Widespread somatic L1 retrotransposition occurs early during gastrointestinal cancer evolution.</a> Genome Research. 25 (10), 1536-1545 (2015).</li>
<li>Beck, C. R., Garcia-Perez, J. L., Badge, R. M., Moran, J. V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((LINE-1+elements+in+structural+variation+and+disease%5BTitle%5D)+AND+%22Annual+Review+of+Genomics+and+Human+Genetics%22%5BJournal%5D)">LINE-1 elements in structural variation and disease.</a> Annual Review of Genomics and Human Genetics. 12, 187-215 (2011).</li>
<li>Philippe, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Activation+of+individual+L1+retrotransposon+instances+is+restricted+to+cell-type+dependent+permissive+loci%5BTitle%5D)+AND+%22eLife%22%5BJournal%5D)">Activation of individual L1 retrotransposon instances is restricted to cell-type dependent permissive loci.</a> eLife. 5, e13926(2016).</li>
<li>Deininger, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+comprehensive+approach+to+expression+of+L1+loci%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">A comprehensive approach to expression of L1 loci.</a> Nucleic Acids Research. 45 (5), e31-e31 (2017).</li>
<li>Jin, Y., Tam, O. H., Paniagua, E., Hammell, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((TEtranscripts%3A+a+package+for+including+transposable+elements+in+differential+expression+analysis+of+RNA-seq+datasets%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">TEtranscripts: a package for including transposable elements in differential expression analysis of RNA-seq datasets.</a> Bioinformatics. 31 (22), 3593-3599 (2015).</li>
<li><a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=%22Agilent+RNA+6000+Nano+Kit+Guide%22">Agilent RNA 6000 Nano Kit Guide</a>. , Agilent. (2017).</li>
<li>Mueller, O. L., Schroeder, A. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aMueller%2C+OL+%22RNA+Integrity+Number+%28RIN%29+%E2%80%93Standardization+of+RNA+Quality+Control%22">RNA Integrity Number (RIN) –Standardization of RNA Quality Control</a>. , Agilent Technologies. (2016).</li>
<li>Robinson, J. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrative+genomics+viewer%5BTitle%5D)+AND+%22Nature+Biotechnology%22%5BJournal%5D)">Integrative genomics viewer.</a> Nature Biotechnology. 29, 24(2011).</li>
<li>Speek, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Antisense+promoter+of+human+L1+retrotransposon+drives+transcription+of+adjacent+cellular+genes%5BTitle%5D)+AND+%22Molecular+Cellular+Biology%22%5BJournal%5D)">Antisense promoter of human L1 retrotransposon drives transcription of adjacent cellular genes.</a> Molecular Cellular Biology. 21 (6), 1973-1985 (2001).</li>
<li>Belancio, V. P., Deininger, L., Roy-Engel, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((LINE+dancing+in+the+human+genome%3A+transposable+elements+and+disease%5BTitle%5D)+AND+%22Genome+Medicine%22%5BJournal%5D)">LINE dancing in the human genome: transposable elements and disease.</a> Genome Medicine. 1 (10), 97-97 (2009).</li>
<li>Iskow, R. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Natural+Mutagenesis+of+Human+Genomes+by+Endogenous+Retrotransposons%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Natural Mutagenesis of Human Genomes by Endogenous Retrotransposons.</a> Cell. 141 (7), 1253-1261 (2010).</li>
<li>Scott, E. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+hot+L1+retrotransposon+evades+somatic+repression+and+initiates+human+colorectal+cancer%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">A hot L1 retrotransposon evades somatic repression and initiates human colorectal cancer.</a> Genome Research. 26 (6), 745-755 (2016).</li>
<li>Kines, K. J., Sokolowski, M., deHaro, D. L., Christian, C. M., Belancio, V. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Potential+for+genomic+instability+associated+with+retrotranspositionally-incompetent+L1+loci%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">Potential for genomic instability associated with retrotranspositionally-incompetent L1 loci.</a> Nucleic Acids Research. 42 (16), 10488-10502 (2014).</li>
</ol>

Os autores não têm nada a revelar.

A atividade L1 tem demonstrado causar dano genético e instabilidade contribuindo para a doença27,28,29. Das aproximadamente 5.000 cópias L1 full-length, somente algumas dúzias evolutionariamente jovens L1s respondem para a maioria da atividade do retrotransposição2. Entretanto, há uma evidência que mesmo alguns mais velhos, retrotransposicionalmente-incompentent L1s são ainda capazes de produzir proteínas prejudiciais do ADN30. Para apreciar inteiramente o papel de L1s na instabilidade e na doença genomic, a expressão L1 a nível locus-específico deve ser compreendida. Entretanto, o fundo elevado de seqüências L1-relacionadas incorporados em outros RNAs não relacionados à retrotransposição L1 levanta um desafio significativo em interpretar a expressão L1 autêntica. Outro desafio na identificação e, portanto, compreensão de padrões de expressão de loci L1 individuais ocorre por causa de sua natureza repetitiva que não permite que muitas seqüências de leitura curtas para mapear para um único Locus único. Para superar esses desafios, desenvolvemos a abordagem descrita acima na identificação da expressão de Locos L1 individuais usando dados de RNA-Seq.
Nossa abordagem filtra o alto nível (mais de 99%) de ruído transcricional gerado a partir de sequências L1 que não estão relacionadas com a retrotransposição L1, tomando uma série de passos. O primeiro passo envolve a preparação do RNA citoplasmático. Selecionando para o RNA cytoplasmic, as leituras L1-relacionadas encontradas dentro do mRNA intronic expressado no núcleo são esgotadas significativamente. Na preparação da biblioteca de sequenciamento, outra etapa tomada para reduzir o ruído transcricional não relacionada ao L1s incluem a seleção de transcrições poliadeniladas. Isso remove o ruído de transcrição relacionado ao L1 encontrado em espécies não mRNA. Uma outra etapa inclui arranjar em seqüência vertente-específica a fim identificar e eliminar transcritos L1-relacionados do antisense. O uso de uma anotação para o L1s completo com regiões de promotor funcionais ao identificar o número de transcritos RNA-Seq que mapeiam para L1s também elimina o ruído de fundo que de outra forma se origina de L1s truncado. Finalmente, a última etapa crítica na eliminação do ruído transcricional das sequências L1 não relacionadas à retrotransposição L1 é a Curação manual de L1s de comprimento total identificada para ter mapeado transcrições de RNA-Seq. A Curação manual envolve a visualização de cada lócus L1 bioinformaticamente identificado-a-ser-expresso no contexto de seu ambiente genómico circunvizinho para confirmar que a expressão se origina do promotor L1. Esta aproximação foi aplicada a DU145, uma linha celular do tumor da próstata. Mesmo com todas as etapas relacionadas à preparação adotadas para reduzir o ruído de fundo, aproximadamente 50% dos Locos L1 identificados bioinformaticamente em DU145 foram rejeitados como ruído de fundo L1 proveniente de outras fontes transcricional (Figura 4), enfatizando o rigor exigido para produzir resultados confiáveis. Essa abordagem usando a Curação manual é trabalhosa, mas necessária no desenvolvimento deste pipeline para avaliar e compreender o ambiente genómico em torno de um L1 de comprimento total. Os próximos passos incluem a redução da quantidade de curadoria manual necessária automatizando algumas das regras de curadoria, embora devido à natureza ainda não completamente conhecida da expressão genômica, fontes de expressão não anotadas no genoma de referência, regiões de baixa mappability, e mesmo complicando fatores envolvidos com a construção de um genoma da referência não é ser possível automatizar inteiramente a Curação L1 neste tempo.
O segundo desafio na identificação da expressão de Locos L1 individuais com sequenciamento relaciona-se com o mapeamento de transcrições L1 repetitivas. Nesta estratégia de alinhamento, é necessário que uma transcrição deve alinhar de forma única e colinearmente ao genoma de referência para ser mapeado. Ao selecionar para sequências de extremidade emparelhada que mapeiam concoremente, a quantidade de transcrições que se alinham exclusivamente a Locos L1 encontrados no genoma de referência aumenta. Essa estratégia de mapeamento exclusivo fornece confiança no chamado de mapeamento de leituras especificamente para um único Locus L1, embora potencialmente subestimar a quantidade de expressão de cada L1 repetitivo identificado para ser autenticamente expresso. Para aproximadamente corrigir essa subestimação, um escore de "mappability" para cada Locus L1 baseado em sua mappability foi desenvolvido e aplicado ao número de leituras de transcrição exclusivamente mapeadas (Figura 6). É de notar que, idealmente, a mappability deve ser pontuada para a cobertura completa lê através do comprimento total L1 de acordo com a amostra WGS correspondente. Aqui, nós usamos WGS de pilhas de HeLa para determinar contagens do mappability de cada Locus L1 a fim inflar ou esvaziar leituras que mapeiam ao loci L1 em linhas de pilha do tumor da próstata DU145. Este cálculo do mappability é uma contagem bruta da correção, mas o ' mappability completo da cobertura escolhida de 400 leituras foi determinado com a natureza dinâmica de linhas de pilha do tumor na mente. Pode ser observado na Figura 1 suplementar, que há alguns Locos L1 com HELA WGS com número extremamente alto de leituras mapeadas. Estes provavelmente vêm das seqüências duplicadas do cromossoma dentro de HeLa que não estão dentro do genoma da referência, que é porque aqueles loci não foram escolhidos para ser representativos da cobertura completa do mappability. Em vez disso, foi determinado que a média de 100% de cobertura de leitura ocorre em torno de 400 leituras de acordo com a Figura 1 suplementar e, em seguida, foi assumido que esta média se aplica à linha de células da próstata tumor DU145 também.
Esta estratégia do alinhamento com 100-200 BP lê da tecnologia do RNA-Seq seleciona preferencialmente também para evolutionarmente mais velho L1s dentro do genoma da referência como L1s mais velhos acumularam sobre mutações originais do tempo que os fazem mais mappable. Esta abordagem, portanto, tem sensibilidade limitada quando se trata de identificar o mais jovem de L1s, bem como não-referência, L1s polimóica. Para identificar o mais novo de L1s, sugerimos o uso de 5 ' RACE seleção de transcritos L1 e tecnologia de sequenciamento como PacBio que fazem uso de leituras mais21. Isto permite o mapeamento mais original e conseqüentemente a identificação confiável do expressado, jovem L1s. o uso de RNA-Seq e PacBio abordagens em conjunto pode levar a uma lista mais abrangente de L1s autenticamente expressa. Para identificar os L1s polimoróricos autenticamente expressos, os primeiros passos seguintes incluem a construção e inserção de sequências popoliméricos no genoma de referência.
Os desafios biológicos e técnicos no estudo de sequências de repetição são grandes, embora com o procedimento rigoroso acima para remover o ruído transcricional das sequências L1 não relacionadas ao retrotransposição usando a tecnologia de sequenciamento de RNA, começamos a peneirar os grandes níveis de ruído de fundo transcricional e ser para identificar com confiança e de forma estrita os padrões de expressão L1 e a quantidade no nível do locus individual.

Retrotransposons são elementos repetitivos do ADN que podem "saltar" no genoma em um mecanismo da copiar-e-pasta através dos intermediários do RNA. Um subconjunto de retrotransposons é sabido como elementos INterspersed Long-1 (LINEs/L1s) e compõe um sexto do genoma humano com sobre 500, 0000 cópias1. Apesar de sua abundância, a maioria dessas cópias são defeituosas e truncados com apenas um estimado 80-120 L1 elementos pensado para ser ativo2. Um L1 de comprimento total é de cerca de 6 KB de comprimento com 5 ' e 3 ' regiões não traduzidas, um promotor interno e promotor anti-Sense associado, dois quadros de leitura aberta não sobrepostas (ORFS), e um sinal e cauda Polya3,4,5 . Nos seres humanos, L1s são compo das subfamílias distinguidas pela idade evolutiva com as famílias mais velhas que acumularam umas mutações mais originais da seqüência sobre o tempo comparada à subfamília a mais nova, L1HS6,7. L1s são os únicos retrotransposons autônomos, humanos e seus ORFS codificam uma transcriptase reversa, um endonuclease, e um rnps com as atividades do RNA-ligação e do acompanhante exigidas retrotranspose e inserção no genoma em um processo referido como alvo-aprontado transcrição reversa8,9,10,11,12.
O retrotransposição de L1s foi relatado para causar doenças humanas do germline por uma variedade de mecanismos que incluem o mutagenesis do insertional, as eliminações do local de destino, e os rearranjos13,14,15, 16. recentemente, foi supor que L1s pode desempenhar um papel na oncogênese e/ou progressão tumoral como o aumento da expressão e eventos de inserção deste elemento mutagênico têm sido observados em uma variedade de cânceres epiteliais17,18 . Estima-se que haja uma nova inserção de L1 em cada 200 nascimentos19. Portanto, é imperativo compreender melhor a biologia do expressar ativamente L1s. A natureza repetitiva e abundância de cópias defeituosas encontradas dentro de transcrições de outros genes tornaram este nível de análise desafiador.
Felizmente, com o advento de tecnologias de sequenciamento de alta taxa de transferência, foram feitos avanços para analisar e identificar autenticamente expressando L1s no nível específico do locus. Há umas filosofias de deferimento em como identificar melhor expressado L1s usando o sequenciamento da próxima geração do RNA. Houve apenas duas abordagens razoáveis sugeridas para mapear transcrições L1 no nível Locus-specific. Um centra-se somente na transcrição potencial que lê através do sinal da poliadenilação L1 e em seqüências flanqueando20. Nossa aproximação aproveita-se das diferenças pequenas da seqüência entre elementos L1 e mapeia somente aqueles RNA-Seq lê esse mapa excepcionalmente a um locus21. Ambos os métodos têm limitações em termos de quantificação dos níveis de transcrição. A quantitation pode ser melhorada potencialmente adicionando uma correção para a "mappability original" de cada lócus L121, ou usando algoritmos mais complexos que redistribuir as leituras de vários mapeados que não puderam ser mapeadas exclusivamente para um locus específico22. Aqui, iremos detalhar de forma passo a passo a extração de RNA e o protocolo de sequenciamento e bioinformática da próxima geração para identificar elementos L1 expressos no nível locus-específico. Nossa aproximação toma a vantagem máxima de nosso conhecimento da biologia de elementos funcionais de L1. Isso inclui saber que os elementos funcionais L1 devem ser gerados a partir do promotor L1, iniciado no início do elemento L1, deve ser traduzido no citoplasma e que suas transcrições devem ser colineares com o genoma. Momentaneamente, nós coletamos o RNA fresco, cytoplasmic, selecionamos para transcritos poliadenilado, e utilizamos análises Strand-specific do RNA-Seq para mapear excepcionalmente leituras ao loci L1 no genoma humano da referência. Estas leituras alinhadas então ainda exigem o curadoria manual extensivo para determinar se as leituras do Transcript originam do promotor L1 antes de designar um locus como um L1 autenticamente expressado. Nós aplicamos esta aproximação na amostra da linha celular do tumor da próstata DU145 para demonstrar como identifica relativamente poucos membros L1 ativamente transcritos da massa de cópias inativas.

<table><tbody><tr><td>1 M HEPES</td><td>Affymetrix</td><td>AAJ16924AE</td><td></td></tr><tr><td>5 M NaCl</td><td>Invitrogen</td><td>AM9760G</td><td></td></tr><tr><td>Agilent bioanalyzer 2100</td><td>Agilent technologies</td><td></td><td></td></tr><tr><td>Agilent RNA 6000 Nano Kit</td><td>Agilent technologies</td><td>5067-1511</td><td></td></tr><tr><td>bedtools.26.0</td><td></td><td></td><td>https://bedtools.readthedocs.io/en/latest/content/installation.html</td></tr><tr><td>gravata borboleta-0.12.8</td><td></td><td></td><td>https://sourceforge.net/projects/bowtie-bio/files/bowtie/0.12.8/</td></tr> <tr><td>raspador de células</td><td>Olympus plastics</td><td>25-270</td><td></td></tr><tr><td>Clorofórmio</td><td>Fisher</td><td>C298-500</td><td></td></tr><tr><td>Digitonin</td><td>Research Products International Corp</td><td>50-488-644</td><td></td></tr><tr><td>Etanol</td><td>Fisher</td><td>A4094</td><td></td></tr><tr><td>Gibco (solução salina tamponada com fosfato)</td> <td>Invitrogen</td><td>10-010-049</td><td></td></tr><tr><td>Homogeneizador</td><td>Thomas Scientific</td><td>BBI-8541906</td><td></td></tr><tr><td>IGV 2.4</td><td></td><td></td><td>https://software.broadinstitute.org/software/igv/download</td></tr><tr><td>Isopropanol</td><td>Fisher</td><td>A416-500</td><td></td></tr><tr><td>mac2unix</td><td></td><td></td><td>https://sourceforge.net/projects/cs-cmdtools/files/mac2unix/</td></tr> <tr><td>Q-tips</td><td>Fisher</td><td>23-400-122</td><td></td></tr><tr><td>RNAse solução posterior</td><td>Invitrogen</td><td>AM7022</td><td></td></tr><tr><td>RNaseZap Solução de descontaminação RNase</td><td>Invitrogen</td><td>AM9780</td><td></td></tr><tr><td>samtools-1.3</td><td></td><td></td><td>https://sourceforge.net/projects/samtools/files/</td></tr> <tr><td>sratoolkit.2.9.2</td><td></td><td></td><td>https://github.com/ncbi/sra-tools/wiki/Downloads</td></tr><tr><td>SUPERase&middot; No inibidor de RNase</td><td>Invitrogen</td><td>AM2694</td><td></td></tr><tr><td>Trizol</td><td>Invitrogen</td><td>15-596-018</td><td></td></tr><tr><td>Água (DNASE, livre de RNASE)</td> <td>Fisher</td><td>BP2484100</td><td></td></tr></tbody></table>

sequenciamento de nova geração de rna e um fluxo de análise de bioinformática para identificar expressão de line-1s a nível locus-específico

RNA Next-Generation Sequencing and a Bioinformatics Pipeline to Identify Expressed LINE-1s at the Locus-Specific Level

Os elementos INterspersed longos-1 (LINEs/L1s) são elementos repetitivos que podem copiar e inserir aleatòria no genoma tendo por resultado a instabilidade e a mutagenese genomic. Compreender os padrões de expressão de Locos L1 no nível individual irá emprestar para a compreensão da biologia deste elemento mutagênico. Este elemento autônomo compõe uma parcela significativa do genoma humano com mais de 500.000 cópias, embora 99% são truncados e defeituosos. Entretanto, sua abundância e número dominante de cópias defeituosas fazem-lhe o desafio identificar autenticamente expressado L1s das seqüências L1-relacionadas expressadas como parte de outros genes. Também é desafiador identificar qual Locus L1 específico é expresso devido à natureza repetitiva dos elementos. Superando esses desafios, apresentamos uma abordagem Bioinformatica de RNA-Seq para identificar a expressão de L1 no nível específico do locus. Em resumo, nós coletamos o RNA cytoplasmic, selecionamos para transcritos poliadenilado, e utilizamos análises Strand-specific do RNA-Seq para mapear excepcionalmente leituras ao loci L1 no genoma humano da referência. Nós visualmente cura cada Locus L1 com leituras excepcionalmente mapeadas para confirmar a transcrição de seu próprio promotor e ajustar as leituras de transcrição mapeadas para dar conta da mappability de cada Locus L1 individual. Esta aproximação foi aplicada a uma linha celular do tumor da próstata, DU145, para demonstrar a habilidade deste protocolo de detectar a expressão de um número pequeno dos elementos L1 full-length.

As etapas descritas acima e descritas graficamente na Figura 1 foram aplicadas a uma linha celular de tumor prostático humano DU145. A amostra do RNA era citomicamente preparado e era Next-Gen seqüenciada em um poli-a selecionou, vertente-específico, protocolo do par-fim. Usando o bowtie, os arquivos de sequenciamento de fim pareado foram alinhados, permitindo apenas correspondências exclusivas nas quais a leitura de extremidade emparelhada correspondeu melhor a uma localização genômica em comparação a qualquer outro local genômico. Os arquivos de seqüência DU145 foram alinhados ao genoma de referência humana criando um arquivo Bam, que está disponível mediante solicitação do autor. Usando bedtools, os dados foram extraídos dos arquivos Bam separados por vertente DU145 no número de leituras que mapeadas para o comprimento total L1s. Essas leituras foram classificadas em uma planilha de maior a menor e manualmente curadoria examinando o ambiente genómico em torno de cada lócus L1 em IGV para confirmar sua autenticidade (tabela suplementar 1). Se uma amostra foi curada para ser autenticamente expressada, era verde codificado por cores com uma explicação para sua aceitação na coluna mais à direita. Exemplos de loci L1 aceitos para serem expressos autenticamente seguindo as diretrizes descritas na seção de métodos são mostrados na Figura 2a-b. Se uma amostra foi rejeitada para ser autenticamente expressa, ele foi codificado por cores como vermelho com o motivo da rejeição na coluna mais à direita. Exemplos de loci L1 rejeitados por causa da expressão de um promotor diferente de suas próprias diretrizes descritas na seção de métodos são detalhados na Figura 2c-e.
Aqui, somente o L1s completo com uma região intacta do promotor foi estudado. Se esta distinção não é feita, uma grande fonte do ruído transcricional que origina de L1s truncado é introduzida. Exemplos de L1s truncados em DU145 são mostrados na Figura 3a-b onde foram identificados como tendo leituras de RNA-Seq exclusivamente mapeadas. Em IGV, entretanto, é aparente que aquelas transcrições não estiveram iniciadas do L1 truncado, mas da inclusão da seqüência L1 em um gene ou a jusante de um gene expressado.
Globalmente, em DU145, a percentagem de loci L1 de comprimento total e leituras que são rejeitadas como autenticamente expressas L1s após a Curação manual é de aproximadamente 50% (tabela suplementar 2) demonstrando o alto nível de leituras de transcrição mapeada L1 que seria caso contrário, ser gravado como falsos positivos sem Curação manual. Especificamente, em DU145 havia 114 total Full-Length L1 loci para ter exclusivamente mapeado leituras na direção sentido com um total de 3.152 leituras, mas havia apenas 60 loci identificado para ser expressa fora de seu próprio promotor após a Curação manual com 1.879 leituras ( Tabela suplementar 1). Este é o caso, mesmo quando as medidas foram tomadas para reduzir a expressão irrelevante para a biologia L1, selecionando para mRNA citoplasmática. Note-se que o locus com o mais alto nível de transcrições mapeadas em DU145 foi rejeitado por não ser um L1 autenticamente expresso (Figura 4). Em geral, o número de transcrições mapeadas para intervalos específicos de L1 localizam similarmente entre os Locos L1 aceitos e rejeitados como autenticamente expressos após a Curação manual (Figura 4).
Após a Curação manual, o número de leituras que mapeiam exclusivamente para o loci L1 específico autenticamente expresso em DU145 variam de 175 leituras a um corte mínimo arbitrariamente escolhido de 10 leituras (Figura 5). Essa abordagem de identificar as leituras de transcrição mapeadas exclusivamente para L1s limita a capacidade de quantificar a expressão com precisão. Para isso, foi criado um fator de correção para cada Locus baseado em sua mappability. Para criar esse fator de correção, os primeiros bedtools foram usados para extrair o número de leituras mapeadas exclusivamente do arquivo de Bam genômica HeLa que alinhou a todos os Locos L1 de comprimento total e grafou aqueles loci das leituras de transcrição mapeadas mais altas para as mais baixas (suplementar Figura 1). Foi designado arbitrariamente que L1s com 400 leituras tiveram o mappability cheio da cobertura. O número de leituras capazes de mapear para um locus L1 na amostra de sequenciamento genômica de HeLa foi dimensionado em relação a 400 leituras e que o número em escala foi então multiplicada para o número de leituras que mapearam para cada loci L1 autenticamente expresso em DU145 (tabela suplementar 2) . Como esperado, os elementos L1 que apresentaram maiores escores de correção para a mappability vieram de subfamílias mais jovens como L1PA2 (tabela suplementar 2). Uma vez que as leituras foram ajustadas para escores de mappability em cada locus, a quantitação para a expressão para a maioria dos Locos aumentou (Figura 6). O número de leituras que mapearam exclusivamente para o loci L1 específico autenticamente expresso com correções de mappability em DU145 variou de 612 a 4 leituras e houve uma re-ordenação de loci expressando o mais alto para o mais baixo (Figura 6).
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/59771/59771fig01v2.jpg"> Figura 1: esquema de fluxo de trabalho. Graficamente descritos são as etapas para identificar L1s expressas em uma amostra humana. Observe que as etapas 1 e 2 não precisam ser repetidas se os arquivos apropriados já estão disponíveis. Estes arquivos apropriados podem ser baixados de suplemento arquivo 1a-b e suplemento arquivo 2. As caixas em vermelho indicam as etapas em que o programa de cobertura de bedtools é usado para contar o número de leituras mapeando para L1s na mesma direção de sentido. Estes loci com sentido orientado mapeamento leituras são o L1s que deve ser manualmente curadoria. <a href="https://www.jove.com/files/ftp_upload/59771/59771fig01large.jpg" target="_blank">Por favor clique aqui para ver uma versão maior desta figura.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/59771/59771fig02v2.jpg"> Figura 2: exemplos de Locos L1 com curadoria em DU145. Carregado em IGV são o genoma de referência, o arquivo de anotação de comprimento total L1 GFF correspondente a versão do genoma de referência (suplemento arquivo 1), o arquivo DU145 Bam e, por último, o arquivo genômica Hela Bam para avaliar a mappability, que estão todos disponíveis em cima do autor Solicitação. As setas foram adicionadas para auxiliar na visualização da direção do L1 anotado. Setas e leituras em vermelho são orientadas em seqüência da direita para a esquerda. Setas e leituras em azul são orientadas em seqüência da esquerda para a direita. a) em IGV, este Locus L1 parece ser expressa fora de seu próprio promotor, pois não há leituras a montante do L1 na orientação sentido para mais de 5 KB. Este L1 tem baixa mappability, não está em um gene, e tem a evidência da atividade esperada do promotor do antisense26. b) em IGV, este Locus L1 parece ser expressa fora de seu próprio promotor como não há leituras a montante do L1 na orientação sentido para mais de 5 KB. Este L1 tem baixa mappability e está dentro de um gene de direção oposta. c) em IGV, este Locus L1 foi rejeitado como um L1 expresso porque há leituras ascendentes na mesma orientação dentro de 5 KB. Este L1 está dentro de um gene do mesmo sentido assim que o Transcript lê é originando o mais provável do promotor do gene expressado. d) em IGV, este Locus L1 foi rejeitado como um L1 expresso porque há leituras ascendentes na mesma orientação dentro de 5 KB. Este L1 é a jusante de um gene altamente expressado no mesmo sentido assim que o Transcript lê é originando o mais provável do promotor daquele gene expresso e estendendo além do terminador normal do gene. e) em IGV, este Locus L1 foi rejeitado como um L1 expresso porque há leituras ascendentes na mesma orientação dentro de 5 KB. Este L1 não está dentro ou perto de um gene anotado no gene de referência assim que a origem destas transcrições dentro e a montante do elemento L1 sugerem um promotor un-anotated. <a href="https://www.jove.com/files/ftp_upload/59771/59771fig02large.jpg" target="_blank">Por favor clique aqui para ver uma versão maior desta figura.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/59771/59771fig03v2.jpg"> Figura 3: o ruído de fundo origina do L1s truncado também. Nossa anotação L1 não inclui truncado L1s como eles são uma fonte principal de ruído de fundo. As setas foram adicionadas para auxiliar na visualização da direção do L1 anotado. Setas e leituras em azul são orientadas em seqüência da esquerda para a direita. a) demonstrado é um exemplo de um L1 truncado na SUFAMÍLIA L1MB5 que é 2706 bps. Em IGV é aparente que as leituras originam da extensão a jusante de um gene expressado. b) mostrado é outro exemplo de um L1 truncado. Este L1 é um L1PA11 que é 4767 bps de comprimento. Em IGV é evidente que as leituras que mapeam excepcionalmente ao L1 originam do exon expressado, que o L1 está dentro. <a href="https://www.jove.com/files/ftp_upload/59771/59771fig03large.jpg" target="_blank">Por favor clique aqui para ver uma versão maior desta figura.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/59771/59771fig04v2.jpg"> Figura 4: o transcrito lê esse mapa excepcionalmente a todo o L1s intacto completo no genoma humano expressado na linha de pilha do tumor da próstata DU145. No preto estão os Locus específicos a ser identificados como autenticamente expressos depois que o curadoria manual e no vermelho é o loci específico a ser rejeitado como leituras autenticamente expressas após o curadoria manual. Em cinza são loci com menos de dez leituras de mapeamento para cada um. Como esses loci representam uma pequena fração de leituras de transcrição, eles não foram manualmente curate. As marcas de escala do eixo x denotam cada 100 de comprimento total, L1s intacto. aproximadamente 4.500 loci não são mostrados graficamente como eles tinham zero mapeado leituras. <a href="https://www.jove.com/files/ftp_upload/59771/59771fig04large.jpg" target="_blank">Por favor clique aqui para ver uma versão maior desta figura.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/59771/59771fig05v2.jpg"> Figura 5: transcrição lê que o mapa exclusivamente para autenticamente expressos full-length L1s intacta na linha celular tumor de próstata DU145. São mostrados os números de leituras de transcrição que mapeiam para Locos específicos em células DU145 após a Curação manual. <a href="https://www.jove.com/files/ftp_upload/59771/59771fig05large.jpg" target="_blank">Por favor clique aqui para ver uma versão maior desta figura.</a>
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/59771/59771fig06v2.jpg"> Figura 6: lê o mapeamento para o L1 autenticamente expresso quando ajustado por mappability. São mostrados os números de leituras de transcrição ajustadas por pontuações de mappability específicas do loci que mapeiam para Locos L1 manualmente curados em células DU145. <a href="https://www.jove.com/files/ftp_upload/59771/59771fig06large.jpg" target="_blank">Por favor clique aqui para ver uma versão maior desta figura.</a>
Arquivo suplementar 1: anotações para o corpo inteiro, L1s humano intacto de acordo com a orientação. a) FL-L1-BLAST_RM_minus. GFF. b) FL-L1-BLAST_RM_plus. GFF. <a href="https://www.jove.com/files/ftp_upload/59771/SupplementalFile1.zip" target="_blank">Por favor, clique aqui para baixar este arquivo.</a>
Arquivo suplementar 2: scripts de supercomputador usados para automatizar o pipeline de Bioinformática detalhado na seção 4. <a href="https://www.jove.com/files/ftp_upload/59771/SupplementalFile_2.docx" target="_blank">Por favor, clique aqui para baixar este arquivo.</a>
Figura suplementar 1: amostra de DNA genômica usada para determinar a mappability L1. São mostrados o número de leituras de transcrição genômica da amostra de linha de célula HeLa que mapeiam exclusivamente para todos os 5.000 locos de comprimento total L1 no genoma. Foi designado que um L1 tem o mappability cheio da cobertura quando 400 lê o mapa ao L1. <a href="https://www.jove.com/files/ftp_upload/59771/SupplementalFigure_1.pdf" target="_blank">Por favor clique aqui para baixar esta figura.</a>
Tabela suplementar 1: Curação manual de L1s em DU145. <a href="https://www.jove.com/files/ftp_upload/59771/Supplemental_Table_1.xlsx" target="_blank">Por favor clique aqui para baixar esta tabela.</a>
Tabela suplementar 2: L1s curadoria em DU145 com ajuste de mappability. <a href="https://www.jove.com/files/ftp_upload/59771/Supplemental_Table_2.xlsx" target="_blank">Por favor clique aqui para baixar esta tabela.</a>

Aqui, apresentamos uma abordagem Bioinformatica e análises para identificar a expressão da linha 1 no nível específico do locus.

Sequenciamento de nova geração de RNA e um fluxo de análise de bioinformática para identificar expressão de LINE-1s a nível locus-específico

Long INterspersed Elements-1 (LINEs/L1s) are repetitive elements that can copy and randomly insert in the genome resulting in genomic instability and ...

Elementos móveis são uma das principais fontes de instabilidade genética humana. Entender sua expressão em diferentes tecidos e condições é fundamental para entender seu impacto no genoma. A vasta das transcrições L1 são o resultado da inclusão passiva de sequências relacionadas ao L1 em outras transcrições que não têm papel no ciclo de vida L1.Nossa abordagem elimina esse fundo irrelevante. Este protocolo pode ser adaptado a estudos de qualquer elemento móvel, ou mesmo vírus em qualquer genoma sequencial. É preciso que haja pelo menos alguma variação de sequência para permitir a discriminação entre loci.A demonstração visual deste método é fundamental para ilustrar a stringency e os cuidados necessários para identificar com confiança elementos repetitivos L1 expressos no nível específico do lócus. Inicie este procedimento com extração de RNA citoplasmático e sequenciamento de próxima geração, conforme descrito no protocolo de texto. Ao selecionar para RNA citoplasmático, as leituras relacionadas ao L1 encontradas no mRNA intronic expresso no núcleo são significativamente esgotadas.Na preparação da biblioteca de sequenciamento, outro passo dado para reduzir o ruído transcricional não relacionado aos L1s inclui a seleção de transcrições poliadenylated. Isso remove o ruído da transcrição relacionada ao L1 encontrado em espécies não-mRNA. Execute o paradigma de alinhamento sequenciando arquivos FASTQ com a amostra de interesse RNA seq usando bowtie1 digitando a linha de comando no terminal Linux.Esta estratégia de alinhamento exige que as transcrições sejam únicas e estejam quase alinhadas com uma busca genômica exaustiva. Esta estratégia fornece confiança no mapeamento de leituras especificamente para um único lócus L1. Strand separe os arquivos BAM de saída Usando samtools e comandos Linux para selecionar para o fio superior e o fio inferior.Observe que os valores reais da bandeira podem variar se não estiver usando protocolos de sequenciamento padrão de próxima geração. Esta etapa de separação da vertente funciona para filtrar o ruído transcricional gerado dentro das sequências L1 que não estão relacionadas à retrotransposição L1, eliminando potenciais leituras mapeadas relacionadas ao Antísceno L1. Gerar contagem de leituras contra anotações para L1 loci usando ferramentas de cama.Primeiro digite a linha de comando para gerar leitura conta para L1s na direção do sentido no fio superior e, em seguida, digite a linha de comando para gerar contagem de leitura para L1s na direção do sentido no fio inferior. As anotações usadas para identificar L1s denotam L1s de comprimento total com regiões de promotores funcionais que trabalham para eliminar ruídos de fundo que de outra forma se originam de L1s truncados. Crie uma planilha para leituras mapeadas para cada lócus L1 anotado.Copiar sobre o arquivo de texto gerado conta leitura que foi criado para a cadeia inferior e rotular a página como minus_bottom. Classifique todas as colunas com base no maior número de leituras encontradas na coluna J.Copy sobre o arquivo de texto gerado que foi criado para a cadeia superior. Classifique todas as colunas com base no maior número de leituras encontradas na coluna J.E rotule a página como top_plus.Crie uma terceira página rotulada como combinada e adicione todos os loci com 10 ou mais leituras de minus_bottom e páginas plus_top. Classificar todas as colunas com base no maior número de leituras encontradas na coluna J.To auxiliar a maplicabilidade de regiões genômicas, especificamente dentro ou perto de L1 loci, arquivos de genoma inteiro emparelhados e sequenciamento das espécies de interesse foram baixados do NCBI e convertidos em arquivos FASTQ, conforme descrito no protocolo de texto. Agora, indexe os arquivos BAM para torná-los visualizados no Visualizador de Genômica Integrativa, abreviado IGV, antes de carregar os arquivos.No IGV carregue o genoma de referência de interesse para visualizar genes anotados. Também carregue o arquivo de anotação para elementos L1 de comprimento total para visualizar a anotação L1, o arquivo BAM para expressão RNA humana, para visualizar transcrições mapeadas a partir da amostra de interesse e o arquivo BAM para mappability genoma humano para avaliar a mappability de regiões genômicas. Remova as linhas de cobertura e junção associadas a cada arquivo BAM.Comprimir os arquivos BAM para a expressão RNA humana e para a mlicabilidade do genoma humano para que todas as faixas de IGV se encaixem em uma tela. O último passo crítico para eliminar o ruído transcricional de sequências L1 não relacionadas à retrotransposição L1 é a criação manual de L1s de comprimento completo identificados como mapeados RNA buscar transcrições. A curadoria manual envolve a visualização de cada lócus L1 expresso no contexto de seu ambiente genômico circundante para confirmar que a expressão se origina do promotor L1.Usando coordenadas de L1 loci listadas na página combinada da planilha, faça a curadoria manual de cada lócus L1 com transcrições exclusivamente mapeadas examinando seu ambiente genômico circundante no IGV. Faça a curadoria de um lócus para ser autenticamente expresso por conta própria se não houver leituras rio acima na direção L1 até cinco quilobases. Rotule a linha de cor verde e observe por que é um L1 autenticamente expresso. Existe uma exceção a esta regra se a região a montante do L1 não for mappable.Se esse for o caso, rotule a linha de cor vermelha e note que a expressão da região a montante do promotor L1 não pode ser avaliada e, portanto, a expressão do L1 não é capaz de ser determinada com confiança. Faça a curadoria de um lócus para não ser expresso autenticamente de seu próprio promotor se houver leituras rio acima de cinco quilobases. Rotule a linha de cor vermelha e observe por que não é um L1 autenticamente expresso. Faça a curadoria de um lócus como falso se for expresso dentro de um intron de um gene expresso na mesma direção, com leituras a montante do L1, se for a jusante de um gene expresso na mesma direção com leituras a montante do L1, ou para padrões de expressão não anotados com leituras a montante da L1. Uma exceção a esta regra se aplica quando há leituras mínimas sobrepondo diretamente o site de início do promotor L1, mas ligeiramente a montante do L1. Se não houver outras leituras a montante de um caso L1 como este, considere este L1 como expresso autenticamente.Rotule a linha verde e note por que é um L1 autenticamente expresso. Cura de um lócus L1 como provável que seja falso se o padrão de leituras mapeadas para o lócus não se correlacionar com as regiões específicas de Mappability L1. Se um L1 é altamente mappable, mas só tem um acúmulo de leituras em uma região condensada dentro do L1, é menos provável que esteja relacionado com a expressão L1 de seu próprio promotor e mais provável de ser de fontes não anotadas como exons ou LTRs. Em casos como este, faça a curadoria do loci como laranja e note por que o lócus é suspeito.Verifique as fontes de acúmulos suspeitos verificando a localização L1 no Navegador de Genomas UCSC. Cura um lócus para não ser expresso autenticamente se estiver dentro de um ambiente genômico de regiões esporadicamente expressas não notadas. As leituras podem ser expressas a 10 quilobases a montante do L1. Mas a cada 10 quilobases ou mais, há leituras mapeadas e algumas dessas leituras se alinham com o L1. Estes L1s provavelmente têm leituras mapeadas devido a padrões não anotados de expressão genômica.Em casos como este, faça a curadoria do loci como vermelho e note por que o lócus é suspeito. Para ajudar a mappability de cada loci L1 determinar o número de leituras exclusivamente mapeadas para L1 loci usando o programa bedtools, a anotação FL-L1 e os dados de sequência genômica alinhada. Designe um lócus L1 para ter cobertura completa quando 400 leituras exclusivas estiverem alinhadas a ele.Determine o fator necessário para escalar para cima ou para baixo leituras genômicas alinhadas a 400 para cada L1 individual. Para ter uma medida de expressão dimensionada de acordo com a mapibilidade de locus L1 individual, multiplique o fator pelo número de leituras de transcrição de RNA que se alinham aos L1s expressos autenticamente. Cada passo é usado para destacar diferenças entre elementos L1 expressos em seu próprio promotor, e todas as maneiras que os elementos L1 podem ser incluídos em outras transcrições que não estão relacionadas com o ciclo de vida L1. Aqui estão as leituras de transcrição que mapeiam exclusivamente todos os L1s intactos de comprimento total no genoma humano expressos na linha celular tumoral da próstata DU145.Em preto estão os loci específicos identificados como autenticamente expressos após a curadoria manual. E em vermelho estão os loci específicos rejeitados como leituras autenticamente expressas após a curadoria manual. Em cinza estão loci com menos de 10 leituras mapeando para cada um.Como esses loci representam uma pequena fração de leituras de transcrição, eles não foram curados manualmente. Aproximadamente 4500 loci não são mostrados graficamente, pois não tinham leituras mapeadas. Após a curadoria manual, o número de leituras que mapeiam exclusivamente para loci L1 específico expresso autenticamente em DU145 variam de 175 leituras a um corte mínimo arbitrariamente escolhido de 10 leituras.Uma vez que as leituras foram ajustadas para escores de mappability em cada lócus, a quantitação para expressão para a maioria dos loci aumentou. O número de leituras que mapearam exclusivamente para loci L1 específico expresso autenticamente com correções de mappability em DU145 variou de 612 a quatro leituras e houve uma reordenação de loci de maior a menor expressagem. Cada passo desempenha um papel crucial na redução do alto nível de ruído de fundo transcricional.No entanto, o passo mais crítico é a curadoria manual de cada lócus L1 para confirmar a transcrição de seu próprio promotor. Aproximadamente 50% de L1 loci identificado bioinformáticamente em células DU145 foram rejeitados como ruído de fundo L1 originário de outras fontes transcricionais, enfatizando o rigor necessário para produzir resultados confiáveis. Para identificar o mais novo dos L1s, sugerimos usar a seleção RACE cinco-prime de transcrições L1 e tecnologia de sequenciamento como o PacBio que fazem uso de leituras mais longas e permite mapeamento mais exclusivo.Com esta abordagem, podemos identificar e quantificar com rigor e confiança os padrões de expressão L1. Isso abre caminho para uma melhor compreensão da regulação do loci L1 individual e do impacto potencial.

rna-next-generation-sequencing-bioinformatics-pipeline-to-identify

Research

JoVE Journal

Genetics

Genética

RNA-de sequenciamento e bioinformática análises foram utilizados para identificar fatores de transcrição significativamente e diferencialmente expressos em subpopulações Lin-CD34 + e Lin-CD34 @ de EMLcells do mouse. Estes fatores de transcrição pode desempenhar um papel importante na determinação da chave entre as células Lin-CD34 @ auto-renovação Lin-CD34 + e parcialmente diferenciados.

Identificação dos fatores chave que regula a auto-renovação e diferenciação em células precursoras hematopoéticas EML por Análise RNA-seqüenciamento

O sequenciamento de amplicon de bissulfito (BSAS) é um método para quantificar a metilação da citosina em regiões genômicas alvo de interesse. Este método usa conversão de bissulfito emparelhada com amplificação de PCR de regiões-alvo antes do sequenciamento de próxima geração para produzir quantificação absoluta da metilação do DNA em um nível específico de base.

Alvejado DNA análise de metilação por Sequenciação de próxima geração

Aqui nós utilizamos a metodologia FISH para rastrear LINE-1 retrotransposição no nível núcleos único dos spreads de cromossomos de linhas celulares HepG2 expressam de forma estável LINE-1 sintético.

Análise de LINE-1 retrotransposição no Núcleo de nível único

Descrevemos um método baseado em sequenciamento de RNA direcionado que inclui a preparação de bibliotecas de cDNA indexadas, hibridização e captura com sondas personalizadas e análise de dados para interrogar transcritos selecionados quanto à expressão gênica, mutações e fusões gênicas. O RNAseq direcionado permite uma avaliação rápida e econômica de transcritos selecionados em um sequenciador de desktop.

Targeted RNA Sequenciação Ensaio para Caracterizar expressão gênica e genômica Alterações

Aqui nós descrevemos um protocolo teve como objetivo investigar o impacto do splicing aberrante na resistência aos medicamentos em tumores sólidos e doenças malignas hematológicas. Para este objetivo, foram analisados os perfis transcriptomic de modelos parentais e resistentes in vitro através de RNA-seq e estabeleceu um método baseado qRT-PCR para validar genes candidatos.

Usando RNA-sequenciação para detectar Novel Splice variantes relacionadas à resistência aos medicamentos em<em&gt; In Vitro</em&gt; modelos de cancro

Nós descrevemos uma técnica de combinação de citometria de fluxo e sequenciamento de alto throughput para identificar regiões de replicação finais do genoma.

G2-seq: Uma alta taxa de transferência baseada no sequenciamento técnica para identificar a tarde replicando regiões do genoma

Alvo da próxima geração sequenciamento é uma abordagem de tempo e custo-eficiente que está se tornando cada vez mais popular na pesquisa da doença e diagnóstico clínico. O protocolo descrito aqui apresenta o fluxo de trabalho complexo necessário para sequenciamento e bioinformática processo usado para identificar variantes genéticas que contribuem para a doença.

Sequenciamento de nova geração direcionado e fluxo de análises de bioinformática para avaliar determinantes genéticos de doenças constitucionais

Aqui, descrevemos um em situ ensaio de hibridização que permite detecção sensível e específica de sequências tão curtas quanto 50 nucleótidos com single-nucleotide resolução a nível de célula única. O ensaio, que pode ser executado manualmente ou automaticamente, pode habilitar a visualização de variantes da tala, sequências curtas e mutações dentro do contexto de tecido.

Visualização de variantes genéticas, metas de curtas e mutações pontuais no contexto morfológico do tecido com um ensaio de RNA In Situ da hibridação

O propósito do presente protocolo é usar uma combinação de computacional e pesquisa do banco para encontrar romance sequências que não podem ser separadas facilmente uma sequência co purificante, que pode ser apenas parcialmente conhecida.

Sequência de romance descoberta pela genômica subtrativa

O objetivo deste protocolo é investigar a evolução e expressão dos genes candidatos usando dados de sequenciamento de RNA.

Um Pipeline bioinforático para investigar a evolução molecular e a expressão genética usando RNA-seq

Aprenda um fluxo de bioinformática para identificar expressão locus-específica de LINE-1s usando RNA-Seq e ferramentas como Bowtie1, SAMtools e IGV. Este método rigoroso elimina ruído de fundo e permite quantificação confiável em células tumorais.

Sequenciamento de nova geração de RNA e um fluxo de análise de bioinformática para identificar expressão de LINE-1s a nível locus-específico (Vídeo) | JoVE Journal

Sequenciamento de nova geração de RNA e um fluxo de análise de bioinformática para identificar expressão de LINE-1s a nível locus-específico

Explorar mais vídeos