Name: Um Pipeline bioinforático para investigar a evolução molecular e a expressão genética usando RNA-seq (Vídeo) | JoVE Journal
Uploaded: 2021-05-28T13:00:00.000+00:00
Duration: 7 min 9 s
Description: Aprenda um pipeline bioinformático para analisar evolução molecular e expressão gênica usando RNA-seq. Este método aplica-se a qualquer organismo e família gênica, auxiliando pesquisas em biologia evolutiva.

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

Agradecemos a Adriana Briscoe, Gil Smith, Rabi Murad e Aline G. Rangel por conselhos e orientações na incorporação de algumas dessas etapas em nosso fluxo de trabalho. Também somos gratos a Katherine Williams, Elisabeth Rebboah e Natasha Picciani por comentários sobre o manuscrito. Este trabalho foi apoiado em parte por uma bolsa de pesquisa médica da Fundação George E. Hewitt para a A.M.M.

Este protocolo segue as diretrizes de cuidados com animais da UC Irvine.
1. Preparação da biblioteca RNA-seq
<ol><li>Isole o RNA usando os seguintes métodos.<ol>
<li>Coletar amostras. Se o RNA for extraído posteriormente, congele a amostra ou coloque na solução de armazenamentoRNA 15 (Tabela de Materiais).</li>
		<li>Eutanize e disseque o organismo para separar tecidos de interesse.</li>
		<li>Extrair RNA total usando um kit de extração e purificar o RNA usando um kit de purificação de RNA(Tabela de Materiais) NOTA: Existem protocolos e kits que podem funcionar melhor para diferentes espécies e tipos de tecidos16,17. Extraímos RNA de diferentes tecidos corporais de uma borboleta18 e uma Hidragelatinosa 19 (ver discussão).</li>
		<li>Meça a concentração e a qualidade do RNA de cada amostra(Tabela de Materiais). Use amostras com números de integridade de RNA (RIN) superiores a 8, idealmente mais próximo de 920 para construir bibliotecas cDNA.</li>
	</ol>
</li>
	<li>Construa biblioteca cDNA e sequência da seguinte forma.<ol>
<li>Construa bibliotecas cDNA de acordo com o manual de instruções de preparação da biblioteca (ver discussão).</li>
		<li>Determinar concentração e qualidade do CDNA(Tabela de Materiais).</li>
		<li>Multiplex as bibliotecas e sequenciá-las.</li>
	</ol>
</li>
</ol>2. Acesse um cluster de computador
NOTA: A análise do RNA-seq requer manipulação de arquivos grandes e é melhor feita em um cluster de computador(Tabela de Materiais).
<ol><li>Faça login na conta de cluster do computador usando o username@clusterlocation de comando em uma janela de aplicativo terminal (Mac) ou PuTTY (Windows).</li>
</ol>3. Obter leituras de RNA-seq
<ol><li>Obtenha leituras de RNA-seq da instalação de sequenciamento ou, para dados gerados em uma publicação, do repositório de dados onde foi depositado (3.2 ou 3.3).</li>
	<li>Para baixar dados de repositórios como o ArrayExpress faça o seguinte:<ol>
<li>Pesquise no site usando o número de adesão.</li>
		<li>Encontre o link para baixar os dados e clique à esquerda e selecione Copiar link.</li>
		<li>Na janela do terminal, digite wget e selecione Colar link para copiar os dados no diretório para análise.</li>
	</ol>
</li>
	<li>Para baixar os dados do NCBI Short Read Archive (SRA) siga estas etapas alternativas:<ol>
<li>No terminal baixe SRA Toolkit v. 2.8.1 usando wget. NOTA: Baixar e instalar programas no cluster do computador pode exigir acesso raiz, entrar em contato com o administrador do cluster do computador se a instalação falhar.</li>
		<li>Termine de instalar o programa digitando tar -xvf $TARGZFILE.</li>
		<li>Pesquise ncbi para obter o número de adesão sra para as amostras que você deseja baixar, ele deve ter o formato SRRXXXXXX.</li>
		<li>Obtenha os dados do RNA-seq digitando [local de srtoolkit]/bin/prefetch SRRXXXXXX NA janela do terminal.</li>
		<li>Para o tipo de arquivos emparelhados [sedetoolkit localização]/bin/fastq-dump -split-files SRRXXXXXX para obter dois arquivos de fastq (SRRXXXXXX_1.FASTQ e SRRXXXXXX_2.FASTQ). NOTA: Para fazer um conjunto Trinity de novo use o comando [sedetoolkit localização]/bin/fastq-dump --defline-seq '@$sn[_$rn]/$ri' -arquivos divididos SRRXXXXXX
</li>
	</ol>
</li>
</ol>4. Aparar adaptadores e leituras de baixa qualidade (opcional)
<ol><li>Instale ou carregue Trimmomatic21 v. 0.35 no cluster de computação.</li>
	<li>No diretório onde os arquivos de dados RNA-seq estão localizados, digite um comando que inclua a localização do arquivo de frasco trimmomatic, os arquivos FASTQ de entrada, arquivos FASTQ de saída e parâmetros opcionais, como comprimento de leitura e qualidade. NOTA: O comando variará pela qualidade e comprimento brutos e desejados das leituras. Para leituras de Illumina 43 bp com primers Nextera, usamos: java -jar /data/apps/trimmomatic/0.35/trimmomatic-0.35.jar PE $READ 1. FASTQ $READ 2. paired_READ1 FASTQ. unpaired_READ1 FASTQ. paired_READ2 FASTQ. unpaired_READ2 FASTQ. FASTQ ILLUMINACLIP:adapters.fa:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:4:17 MINLEN:30.
</li>
</ol>5. Obter montagem de referência
<ol><li>Pesquise no Google, EnsemblGenomes e NCBI Genomas e Nucleotídeos TSA (Transcriptome Shotgun Assembly) para obter um genoma de referência ou transcriptome montado para as espécies de interesse(Figura 1). NOTA: Se um genoma de referência ou transcriptome não estiver disponível ou de baixa qualidade, proceda ao PASSO 6 para gerar um conjunto de novo.</li>
	<li>Se existir um genoma de referência ou transcriptome montado, baixe-o como um arquivo fasta para onde a análise será realizada seguindo as etapas abaixo.<ol>
<li>Encontre o link para baixar o genoma, clique à esquerda e Copy Link.</li>
		<li>Na janela do terminal digite wget e cole o endereço de link. Se disponível, copie também o arquivo GTF e o arquivo FASTA de proteína para o genoma de referência.</li>
	</ol>
</li>
</ol>6. Gerar um conjunto de novo (Alternativa ao Passo 5)
<ol><li>Combine os arquivos RNA-seq READ1 e READ2 fastq para todas as amostras digitando gato *READ1. FASTQ > $all_READ1. FASTQ e gato *READ2. > all_READ2 FASTQ. FASTQ na janela do terminal.</li>
	<li>Instale ou carregue o Trinity22 v.2.8.5 no cluster de computação.</li>
	<li>Gerar e montagem digitando no terminal: Trinity --seqType fq --max_memory 20G --esquerda $all_READ1. FASTQ --direita $all_READ2. O FASTQ.
</li>
</ol>7. Mapa lê para o genoma (7.1) ou de novo transcriptome (7.2)
<ol><li>O mapa lê para o genoma de referência usando STAR23 v. 2.6.0c e RSEM24 v. 1.3.0.<ol>
<li>Instale ou carregue STAR v. 2.6.0c. e RSEM v. 1.3.0 para o cluster de computação.</li>
		<li>Indexe o genoma digitando rsem-prepare-referência --gtf $GENOME. GTF -estrela -p 16 $GENOME. $OUTPUT FASTA.
</li>
		<li>Mapa lê e calcula a expressão para cada amostra digitando rsem-calculate-expression -p 16 --star --emparelhado-end $READ 1. FASTQ $READ 2. $INDEX $OUTPUT FASTQ.
</li>
		<li>Renomeie o arquivo de resultados para algo descritivo usando mv RSEM.genes.results $sample.genes.results.
</li>
		<li>Gerar uma matriz de todas as contagens digitando rsem-generate-data-matrix *[genes/isoforms.results] > $OUTPUT.</li>
	</ol>
</li>
	<li>Mapeie RNA-seq para o conjunto Trinity de novo usando RSEM e bowtie.<ol>
<li>Instale ou carregue Trinity22 v.2.8.5, Bowtie25 v. 1.0.0 e RSEM v. 1.3.0.</li>
		<li>Mapa lê e calcula a expressão para cada amostra digitando [trinity_location]/align_and_estimate_abundance.pl --prep-reference --transcrições $TRINITY. FASTA --seqType fq --esquerda $READ 1. FASTQ - direita $READ 2. FASTQ --est_method RSEM --aln_method bowtie --trinity_mode --output_dir $OUTPUT.
</li>
		<li>Renomeie o arquivo de resultados para algo descritivo usando mv RSEM.genes.results $sample.genes.results.
</li>
		<li>Gerar uma matriz de todas as contagens digitando [trinity_location]/abundance_estimates_to_matrix.pl --est_method RSEM *[genes/isoforms].resultados
</li>
	</ol>
</li>
</ol>8. Identificar genes de interesse
NOTA: As seguintes etapas podem ser feitas com arquivos nucleotídeos ou proteínas FASTA, mas funcionam melhor e são mais simples com sequências proteicas. Pesquisas de BLAST usando proteína para proteína são mais propensas a dar resultados na busca entre diferentes espécies.
<ol><li>Para um genoma de referência, use o arquivo FASTA de proteína do STEP 5.2.2 ou consulte Materiais Suplementares para gerar um recurso genético personalizado GTF.</li>
	<li>Para um de novo transcriptome, gere uma proteína FASTA usando TransDecoder.<ol>
<li>Instale ou carregue TransDecoder v. 5.5.0 no cluser do computador.</li>
		<li>Encontre o quadro de leitura aberto mais longo e a sequência de peptídeos prevista digitando [transdecoder location]/TransDecoder.LongOrfs -t $TRINITY. O FASTA.
</li>
	</ol></li>
	<li>Procure no NCBI Genbank por homólogos em espécies intimamente relacionadas.<ol>
<li>Abra uma janela do navegador de internet e vá para https://www.ncbi.nlm.nih.gov/genbank/.</li>
		<li>Na barra de pesquisa digitam o nome do gene de interesse e o nome de espécies intimamente relacionadas que foram sequenciadas ou gêneros ou filo. À esquerda da barra de pesquisa selecione proteína e clique em pesquisar.</li>
		<li>Extrair sequências clicando em Enviar e, em seguida, selecionar Arquivo. Em Formato, selecione FASTA e clique em Criar Arquivo.
</li>
		<li>Mova o arquivo FASTA de homólogos para o cluster do computador digitando scp $FASTA username@clusterlocation:/$DIR em uma janela de terminal local ou use FileZilla para transferir arquivos de e para computador e cluster.</li>
	</ol></li>
	<li>Procure genes de candidatos usando BLAST+26.<ol>
<li>Instale ou carregue BLAST+ v. 2.8.1 no cluster do computador.</li>
		<li>No cluster do computador, faça um banco de dados BLAST a partir do genoma ou transcriptome proteína traduzida FASTA digitando [localização BLAST+]/makeblastdb -em $PEP. FASTA -dbtype prot -out $OUTPUT
</li>
		<li>BLAST as sequências genéticas homólogos do NCBI para o banco de dados da espécie de interesse digitando [localização BLAST+]/blastp -db $DATABASE -query $FASTA -evalue 1e-10 -outfmt 6 -max_target_seqs 1 -out $OUTPUT.</li>
		<li>Exibir o arquivo de saída usando maiso comando . Copie iDs genéticos exclusivos da espécie de interesse para um novo arquivo de texto.</li>
		<li>Extrair as sequências de genes candidatos digitando perl -ne 'se(/^>(\S+)/){$c=$i{$1}}$c?print:chomp;$i{$_}=1 se @ARGV' $gene_id.txt $PEP. > $OUTPUT FASTA.
</li>
	</ol></li>
	<li>Confirme a anotação genética usando BLAST recíproco.<ol>
<li>No navegador da internet vá para https://blast.ncbi.nlm.nih.gov/Blast.cgi.</li>
		<li>Selecione tblastn, em seguida, cole as sequências do candidato, selecione o banco de dados de sequência de proteínas não redundante e clique em BLAST.</li>
	</ol></li>
	<li>Identifique genes adicionais anotando todos os genes do genoma ou transcriptome com termos de ontologia genética (GO) (ver discussão).<ol>
<li>Transfira a proteína FASTA para o computador local.</li>
		<li>Baixe e instale o Blast2GO27,28,29 v. 5.2 para o computador local.</li>
		<li>Abra o Blast2GO, clique em Arquivo,vá para Carga,vá para Sequências de carga,clique em Load Fasta File (fasta). Selecione o arquivo FASTA e clique em Carregar.</li>
		<li>Clique em Blast, escolha NCBI Blaste clique em Next. Editar parâmetros ou clicar em Next, editar parâmetros e clicar em Executar para encontrar a descrição genética mais semelhante.</li>
		<li>Clique em mapear e clique em Executar para pesquisar anotações de Gene Ontology para proteínas semelhantes.</li>
		<li>Em seguida, clique em interpro, selecione EMBL-EBI InterProe clique em Next. Editar parâmetros ou clicar em Seguire clique em Executar para procurar assinaturas de famílias e domínios de genes conhecidos.</li>
		<li>Exporte as anotações clicando em Arquivo,selecione Exportar,clique em Tabela de Exportação. Clique em Procurar, nomeie o arquivo, clique em Salvar, clique em Exportar.</li>
		<li>Pesquise na tabela de anotação termos de interesse do GO para identificar genes adicionais de candidatos. Extrair as sequências do arquivo FASTA (STEP 8.4.5)</li>
	</ol></li>
</ol>9. Árvores filogenéticas
<ol><li>Baixe e instale MEGA30 v. 7.0.26 no seu computador local.</li>
	<li>Abra MEGA, clique em Alinhar,clique em Editar/Construir Alinhamento,selecione Criar um novo alinhamento clique em OK, selecione Proteína.</li>
	<li>Quando a janela de alinhamento for aberta, clique em Editar,clique em Inserir sequências do arquivo e selecione o FASTA com sequências proteicas de genes candidatos e prováveis homólogos.</li>
	<li>Selecione todas as sequências. Encontre o símbolo do braço e passe o mouse sobre ele. Deve-se dizer Alinhar sequências usando algoritmo MUSCLE31. Clique no símbolo do braço e clique em Alinhar proteína para alinhar as sequências. Editar parâmetros ou clicar em OK para alinhar usando parâmetros padrão.</li>
	<li>Inspecione visualmente e faça quaisquer alterações manuais e, em seguida, salve e feche a janela de alinhamento.</li>
	<li>Na janela PRINCIPAL MEGA, clique em Modelos,clique em Encontrar os melhores modelos de DNA/proteína (ML), selecione o arquivo de alinhamento e selecione parâmetros correspondentes como: Análise: Seleção de Modelo (ML), Árvore para usar: Automática (árvore de junção vizinha), Método Estatístico: Probabilidade Máxima, Tipo de Substituição: Aminoácido, Tratamento de dados Gap/ausente: Use todos os sites, Filtro do site do ramo: Nenhum.
</li>
	<li>Uma vez determinado o melhor modelo para os dados, vá para a janela MEGA principal. Clique em Phylogeny e clique em Árvore de Máxima Probabilidade e selecione o alinhamento, se necessário. Selecione os parâmetros apropriados para a árvore: Método estatístico: Máxima Probabilidade, Teste de Filogenia: Método Bootstrap com 100 réplicas, tipo de substituição: aminoácido, modelo: LG com Freqs. (+F), taxas entre os locais: gama distribuída (G) com 5 categorias gama discretas, tratamento de dados gap/missing: use todos os locais, método heurístico ML: Near-Neighbor-Interchange (NNI).
</li>
</ol>10. Visualize a expressão genética usando TPM
<ol><li>Para Trinity, no cluster do computador vá para o diretório onde abundance_estimates_to_matrix.pl foi executado e uma das saídas deve ser matriz. TPM.not_cross_norm. Transfira este arquivo para o seu computador local. NOTA: Consulte Materiais Suplementares para normalização da amostra cruzada.</li>
	<li>Para TPMs de uma análise de genoma siga os passos abaixo.<ol>
<li>No cluster do computador, vá para o local de instalação do RSEM. Copie a matriz rsem-generate-data,digitando scp rsem-generate-data-matrix rsem-generate-TPM-matrix. Use nano para editar o novo arquivo e alterar "meu $offsite = 4" de 4 para 5 para TPM, ele agora deve ler "meu $offsite = 5".</li>
	</ol>
</li>
	<li>Vá para o diretório onde os arquivos de saída RSEM .genes.results estão e agora use rsem-generate-TPM-matrix *[genes/isoforms.results] > $OUTPUT para gerar uma matriz TPM. Transferir resultados para um computador local.</li>
	<li>Visualize os resultados em ggplot2.<ol>
<li>Baixe R v. 4.0.0 e RStudio v. 1.2.1335 para um computador local.</li>
		<li>Abra rstudio à direita da tela vá para a guia Pacotes e clique em Instalar. Digite ggplot2 e clique em instalar.</li>
		<li>Na janela do script R lida na tabela TPM digitando dados<-read.table ("$tpm.txt", cabeçalho = T)
</li>
		<li>Para gráficos de barra semelhantes à Figura 4 digitam algo semelhante a: p<-ggplot() + geom_bar(aes(y=TPM, x=Symbol, fill=Tissue), data=data, stat="identity") preenchimento<-c("#d7191c", "#fdae61", "#ffffbf", "#abd9e9", "#2c7bb6") p<-p+scale_fill_manual(valores=preenchimento) p + tema (axis.text.x = element_text(ângulo = 90))
</li>
	</ol>
</li>
</ol>

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

O objetivo deste protocolo é fornecer um esboço das etapas para caracterizar uma família genética usando dados RNA-seq. Estes métodos têm sido comprovados para funcionar para uma variedade de espécies e conjuntos de dados4,34,35. O gasoduto aqui estabelecido foi simplificado e deve ser fácil o suficiente para ser seguido por um novato em bioinformática. O significado do protocolo é que ele delineia todas as etapas e programas ne...

Os avanços nas tecnologias de sequenciamento facilitaram o sequenciamento de genomas e transcriptomes de organismos não-modelos. Além da maior viabilidade do sequenciamento do DNA e do RNA de muitos organismos, uma abundância de dados está disponível publicamente para estudar genes de interesse. O objetivo deste protocolo é fornecer passos bioinforáticos para investigar a evolução molecular e a expressão de genes que possam desempenhar um papel importante no organismo de interesse.
Investigar a evolução de um gene ou gene familiar pode fornecer uma visão da evolução dos sistemas biológicos. Membros de uma família genética são tipicamente determinados pela identificação de motivos conservados ou sequências genéticas homólogos. A evolução da família genética foi previamente investigada usando genomas de organismos modelos distantes1. Uma limitação para essa abordagem é que não está claro como essas famílias genéticas evoluem em espécies intimamente relacionadas e o papel de diferentes pressões seletivas ambientais. Neste protocolo, incluímos uma busca por homólogos em espécies intimamente relacionadas. Ao gerar uma filogenia a um nível de filogênio, podemos notar tendências na evolução da família genética, como a de genes conservados ou duplicações específicas de linhagem. Neste nível, também podemos investigar se genes são ortologs ou paralogs. Embora muitos homólogos provavelmente funcionem de forma semelhante entre si, isso não é necessariamente o caso2. A incorporação de árvores filogenéticas nesses estudos é importante para resolver se esses genes homólogos são ortologs ou não. Nos eucariotes, muitos ortologos mantêm funções semelhantes dentro da célula, como evidenciado pela capacidade das proteínas mamíferas de restaurar a função dos ortologs de levedura3. No entanto, há casos em que um gene não ortologos realiza uma função caracterizada4.
As árvores filogenéticas começam a delinear relações entre genes e espécies, mas a função não pode ser atribuída apenas com base nas relações genéticas. Estudos de expressão genética combinados com anotações funcionais e análise de enriquecimento fornecem forte suporte para a função genética. Casos em que a expressão genética pode ser quantificada e comparada entre indivíduos ou tipos de tecidos podem ser mais reveladores da função potencial. O protocolo a seguir segue métodos usados na investigação de genes de opsina em Hydra vulgaris7, mas eles podem ser aplicados a qualquer espécie e qualquer família genética. Os resultados desses estudos fornecem uma base para uma investigação mais aprofundada sobre a função genética e redes genéticas em organismos não-modelo. Como exemplo, a investigação da filogenia das opsinas, que são proteínas que iniciam a cascata de fototransdução, dá contexto à evolução dos olhos e da detecção de luz8,9,10,11. Neste caso, organismos não-modelos, especialmente espécies de animais basais, como cnidários ou ctenoforos, podem elucidar a conservação ou alterações na cascata de fototransdução e na visão através de claes12,13,14. Da mesma forma, determinar a filogenia, expressão e redes de outras famílias genéticas nos informará sobre os mecanismos moleculares subjacentes às adaptações.

<table><tbody><tr><td>Kit Bioanalyzer-DNA</td><td>Agilent</td><td>5067-4626</td><td>materiais de laboratório úmido</td></tr><tr><td>Kit Bioanalyzer-RNA</td><td>Agilent</td><td>5067-1513</td><td>materiais de laboratório úmido</td></tr><tr><td>BLAST+ v. 2.8.1</td><td></td><td></td><td>No cluster de computadores* https://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/</td></tr> <tr><td>Blast2GO (no seu PC)</td><td></td><td></td><td>No computador local https://www.blast2go.com/b2g-register-basic</td></tr><tr><td>boost v. 1.57.0</td><td></td><td></td><td>No cluster de computadores</td></tr><tr><td>Bowtie v. 1.0.0</td><td></td><td></td><td>No cluster de computadores https://sourceforge.net/projects/bowtie-bio/files/bowtie/1.3.0/</td></tr> <tr><td>Cluster de computação (altamente recomendado)</td><td></td><td></td><td>NOTA: As análises de dados genômicos são melhor feitas em um cluster de computação de alto desempenho porque os arquivos são muito grandes.</td></tr><tr><td>Abotoaduras v. 2.2.1</td><td></td><td></td><td>No cluster de computadores</td></tr><tr><td>edgeR v. 3.26.8 (in R)</td><td></td><td></td><td>No Rstudio https://bioconductor.org/packages/release/bioc/html/edgeR.html</td></tr><tr><td>gcc v. 6.4.0</td><td></td><td></td><td>No cluster de computadores</td></tr><tr><td>Java v. 11.0.2</td><td></td><td></td><td>No cluster de computadores</td></tr><tr><td>MEGA7 (no seu PC)</td><td></td><td></td><td>No computador local https://www.megasoftware.net</td></tr><tr><td>MEGAX v. 0.1</td><td></td><td></td><td>No computador local https://www.megasoftware.net</td></tr><tr><td>Kit NucleoSpin RNA II</td><td>Macherey-Nagel</td><td>740955.5</td><td>materiais de laboratório úmido</td></tr><tr><td>perl 5.30.3</td><td></td><td></td><td>No cluster de computador</td></tr><tr><td>python</td><td></td><td></td><td>No cluster de computador</td></tr><tr><td>Qubit 2.0 Fluorômetro</td><td>ThermoFisher</td><td>Q32866</td><td>materiais de laboratório úmido</td></tr><tr><td>R v.4.0.0</td><td></td><td></td><td>Em cluster de computador https://cran.r-project.org/src/base/R-4/</td></tr> <tr><td>RNAlater</td><td>ThermoFisher</td><td>AM7021</td><td>materiais de laboratório úmido</td></tr><tr><td>RNeasy kit</td><td>Qiagen</td><td>74104</td><td>materiais de laboratório úmido</td></tr><tr><td>RSEM v. 1.3.0</td><td></td><td></td><td>Software de computador https://deweylab.github.io/RSEM/</td></tr><tr><td>RStudio v. 1.2.1335</td><td></td><td></td><td>No computador local https://rstudio.com/products/rstudio/download/#download</td></tr><tr><td>Samtools v. 1.3</td><td></td><td></td><td>Software de computador</td></tr><tr><td>SRA Toolkit v. 2.8.1</td><td></td><td></td><td>No cluster de computadores https://github.com/ncbi/sra-tools/wiki/01.-Downloading-SRA-Toolkit</td></tr><tr><td>STAR v. 2.6.0c</td><td></td><td></td><td>No cluster de computadores https://github.com/alexdobin/STAR</td></tr><tr><td>StringTie v. 1.3.4d</td><td></td><td></td><td>Em cluster de computadores https://ccb.jhu.edu/software/stringtie/</td></tr> <tr><td>Transdecoder v. 5.5.0</td><td></td><td></td><td>Em cluster de computadores https://github.com/TransDecoder/TransDecoder/releases</td></tr><tr><td>Trimmomatic v. 0.35</td><td></td><td></td><td>Em cluster de computadores http://www.usadellab.org/cms/?page=trimmomatic</td></tr><tr><td>Trinity v.2.8.5</td><td></td><td></td><td>Em cluster de computador https://github.com/trinityrnaseq/trinityrnaseq/releases</td></tr><tr><td>TRIzol</td><td>ThermoFisher</td><td>15596018</td><td>materiais de laboratório úmido</td></tr><tr><td>TruSeq RNA Library Prep Kit v2</td><td>Illumina</td><td>RS-122-2001</td><td>materiais de laboratório úmido</td></tr><tr><td>Kit TURBO sem DNA</td><td>ThermoFisher</td><td>AM1907</td><td>materiais</td></tr> <tr><td></td><td>de laboratório úmido</td><td></td><td></td></tr><tr><td>*Downloads e instalação no cluster de computadores podem exigir acesso root. Entre em contato com o administrador da rede.</td><td></td><td></td><td></td></tr></tbody></table>

um pipeline bioinforático para investigar a evolução molecular e a expressão genética usando rna-seq

A Bioinformatics Pipeline for Investigating Molecular Evolution and Gene Expression using RNA-seq

Destilar e relatar grandes conjuntos de dados, como dados de genoma inteiro ou transcriptome, é muitas vezes uma tarefa assustadora. Uma maneira de quebrar resultados é focar em uma ou mais famílias genéticas que são significativas para o organismo e estudar. Neste protocolo, delineamos etapas bioinformáticas para gerar uma filogenia e quantificar a expressão de genes de interesse. As árvores filogenéticas podem dar uma visão de como os genes estão evoluindo dentro e entre espécies, bem como revelar ortologia. Esses resultados podem ser aprimorados usando dados de RNA-seq para comparar a expressão desses genes em diferentes indivíduos ou tecidos. Estudos de evolução molecular e expressão podem revelar modos de evolução e conservação da função genética entre espécies. A caracterização de uma família genética pode servir de trampolim para estudos futuros e pode destacar uma importante família genética em um novo genoma ou papel transcriptome.

Os métodos acima são resumidos na Figura 1 e foram aplicados a um conjunto de dados de tecidos hydra vulgaris. H. vulgaris é um invertebrado de água doce que pertence ao filo Cnidaria que também inclui corais, águas-vivas e anêmonas do mar. H. vulgaris pode se reproduzir assexualmente brotando e eles podem regenerar a cabeça e o pé quando bissecto. Neste estudo, buscou-se investigar a evolução e expressão dos genes opsin na Hydra...

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

Distilling and reporting large datasets, such as whole genome or transcriptome data, is often a daunting task. One way to break down results is to focus ...

Este protocolo descreve etapas bioinformáticas para investigar a evolução molecular e a expressão dos genes candidatos. Aqui fornecemos instruções minuciosas para que qualquer pessoa com experiência bioinformática mínima possa passar por este protocolo. Este oleoduto pode ser aplicado a qualquer organismo e qualquer família genética.Uma questão comum ao fazer bioinformática é que os scripts shell falham. Ao tentar este protocolo, certifique-se de ter o software mais atualizado, leia os arquivos de erro e verifique o manual cuidadosamente. Para começar, faça login na conta de cluster do computador em uma janela de aplicativo terminal ou PuTTY.No terminal, baixe a versão 2.8.1 do SR Toolkit usando o Wget e, em seguida, termine de instalar o programa. Pesquise NCBI para obter o número de adesão sra para as amostras desejadas e, em seguida, obtenha os dados de sequência de RNA na janela do terminal. Obtenha dois arquivos FASTQ para o tipo de arquivos de extremidade emparelhado.Encontre o genoma de referência on-line se existir. Para obter um conjunto de referência, digite wget na janela do terminal e cole o endereço de link. Se disponível, copie também o arquivo GTF e o arquivo FASTA de proteína para o genoma de referência.Indexe o genoma, então mapeia leituras e calcula a expressão para cada amostra. Renomeie o arquivo de resultados para algo descritivo e gere uma matriz de todas as contagens. Abra uma janela do navegador de internet e vá para o NCBI GenBank.Na barra de pesquisa, digite o nome do gene de interesse e o nome de espécies intimamente relacionadas que foram sequenciadas. À esquerda da barra de pesquisa, selecione Proteína e clique em Pesquisar. Extrair as sequências clicando em Enviar e, em seguida, selecione Arquivo.Em Formato, selecione FASTA e clique em Criar arquivo. Mova o arquivo FASTA de homologues para o cluster do computador usando uma janela de terminal local ou FileZilla. Em seguida, procure genes candidatos usando BLAST+No cluster do computador, faça um banco de dados BLAST a partir do genoma ou proteína traduzida por transcriptome, FASTA.BLAST as sequências genéticas homólogos do NCBI para o banco de dados da espécie de interesse, em seguida, visualize o arquivo de saída usando mais o comando. Copie iDs genéticos exclusivos da espécie de interesse para um novo arquivo de texto. Extrair as sequências de genes candidatos.Para confirmar a anotação genética usando o BLAST recíproco, vá para a Ferramenta de Pesquisa de Alinhamento Local BLAST, selecione BLASTP e cole as sequências de candidatos, selecione o banco de dados de sequência de proteínas não redundante e clique em BLAST. Abra MEGA, clique em Alinhar, em seguida, Editar Alinhamento de compilação, selecionar Criar um novo alinhamento e clique em OK. Selecione Proteína. Quando a janela Alinhamento for aberta, clique em Editar.Clique em Inserir sequências do arquivo e selecione o FASTA com sequências proteicas de genes candidatos e possíveis homólogos. Selecione Todas as sequências. Encontre o símbolo do braço e passe o mouse sobre ele.Deve dizer alinhar sequências usando algoritmo muscular. Clique no símbolo do braço e clique em Alinhar Proteína para alinhar as sequências Editar parâmetros ou clicar em OK para usar parâmetros padrão. Este protocolo foi aplicado aos tecidos da Hydra vulgaris, que é um invertebrado de água doce que pertence ao filo Cnidaria.Genes de opsina foram investigados para obter uma visão sobre a evolução dos olhos e a detecção de luz em animais. Sequências para genes relacionados à opsina de H.vulgaris e outras espécies foram extraídas em um arquivo FASTA do NCBI GenBank. Os genes de opsina foram alinhados no MEGA, possibilitando identificar as opsinas hydra que estavam faltando um aminoácido conservado de lisina necessário para ligar uma molécula sensível à luz.Uma árvore de maior probabilidade foi gerada usando sequências de opsina da Hydra vulgaris e outras espécies. A filogenia sugere que os genes da opsina estão evoluindo por duplicações específicas de linhagem em cnidários, e potencialmente por duplicação em H.vulgaris. Em seguida, foi realizada uma análise de expressão diferencial no edgeR para investigar a expressão absoluta dos genes opsina.Para determinar se uma ou mais opsinas são reguladas no hipostome, ou cabeça, foram realizadas comparações par-wise de hipostome versus coluna corporal, zona de brotamento, pé e tentáculos foram realizadas. Verificou-se que 1.774 transcrições foram expressas diferencialmente entre o hipostome e a coluna corporal. Os genes que foram regulados em várias comparações foram determinados, e um enriquecimento funcional no Blast2GO foi realizado.Finalmente, a expressão absoluta dos genes opsin foi investigada em diferentes tecidos durante diferentes estágios de brotação, e durante diferentes pontos de tempo de regeneração. A inspeção visual do alinhamento e da árvore confirmará se os genes candidatos pertencem à família de interesse. Genes que são muito diferentes em sequência ou um grupo fora de tudo o resto, provavelmente fazem parte de uma família genética diferente.Os resultados deste protocolo podem ser considerados geradores de hipóteses. Este pipeline pode destacar genes candidatos para estudar funcionalmente em estudos futuros. Depois de explorar a expressão de hidra opsina, agora estamos usando técnicas semelhantes para investigar genes relacionados entre espécies, a fim de identificar semelhanças e diferenças de função.

a-bioinformatics-pipeline-for-investigating-molecular-evolution-gene

Research

JoVE Journal

Biology

Biologia

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

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, 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

Apresentado aqui é um procedimento stepwise para diferenciação in vitro de queratinócitos primários humanos por inibição de contato seguido de caracterização a nível molecular pela análise de RNA-seq.

Caracterização da diferenciação in vitro dos queratinócitos primários humanos pela análise do RNA-Seq

Este método descreve as etapas para melhorar a qualidade e a quantidade de dados sequenciais que podem ser obtidos a partir de amostras de RNA incorporadas à parafina fixa-formalina (FFPE). Descrevemos a metodologia para avaliar com mais precisão a qualidade das amostras de FFPE-RNA, preparar bibliotecas de sequenciamento e analisar os dados das amostras de FFPE-RNA.

Otimização para Sequenciamento e Análise de Amostras de RNA FFPE-RNA degradadas

Foi fornecido um protocolo detalhado de métodos de análise de expressão diferencial para sequenciamento de RNA: limma, EdgeR, DESeq2.

Três métodos diferenciais de análise de expressão para sequenciamento de RNA: limma, EdgeR, DESeq2

Este protocolo permite o controle de qualidade inicial para experimentos de RNA-seq para biólogos de laboratório úmido com experiência limitada em bioinformática.

Rup (RNA-seq Usability Assessment Pipeline) - Controle de qualidade para experimentos de RNA-seq em massa em eucariotos

Aqui nós relatamos a geração de Tre recombinase através dirigido evolução, molecular. Tre recombinase reconhece uma seqüência alvo pré-definido dentro das sequências LTR do pró-vírus HIV-1, resultando na excisão e erradicação do provírus de células humanas infectadas. Enquanto ainda em sua infância, dirigido a evolução molecular irá permitir a criação de enzimas personalizadas que servirão como ferramentas de cirurgia e medicina molecular molecular.

Evolução Molecular do Recombinase Tre

Este procedimento descreve como estabelecer culturas primárias de células hematopoiéticas da medula óssea de camundongos e é seguido por transfecção usando o Gene Pulser sistema de eletroporação MXCell.

Elaboração de culturas de células primárias hematopoéticas de medula óssea de camundongos para Eletroporação

A demonstração do isolamento do cabo do mouse neonatal espinhal para estudos eletrofisiológicos.

Spinal Cord Eletrofisiologia

Aprenda um pipeline bioinformático para analisar evolução molecular e expressão gênica usando RNA-seq. Este método aplica-se a qualquer organismo e família gênica, auxiliando pesquisas em biologia evolutiva.

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

Explorar mais vídeos