Name: Aproveitando os Recursos CyVerse para De Novo Transcriptômica Comparativa de Organismos Não Servidos (não Modelos) (Vídeo) | JoVE Journal
Uploaded: 2017-05-09T14:00:00.000+00:00
Duration: 10 min 41 s
Description: Aprenda a montar e anotar transcriptomas de novo usando recursos CyVerse para organismos não modelos. Este fluxo de trabalho interativo permite analisar dados de RNA-Seq e comparar expressão gênica diferencial.

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

Os autores gostariam de agradecer o financiamento da concessão USDA-NIFA 2013-00984, NSF grant IOS - 1339156, IOS - 1444490 e CyVerse (NSF: DBI - 1265383).

NOTA: O protocolo global foi numerado de acordo com pastas que serão criadas e nomeadas no passo 1.2 ( Figura 1 e 2 ). Este protocolo representa um padrão comparativo de novo transcriptome análise, e cada passo detalhado aqui pode não ser necessário para todos os pesquisadores. Este fluxo de trabalho é documentado completamente em um tutorial complementar wiki, que também contém todos os arquivos adicionais e links para documentos de interesse 3 º partido desenvolvedores para cada pacote de análise ( Tabela 1 ). Links para este material serão incluídos em todo este protocolo para facilitar o acesso a esta informação. Melhores práticas são notas fornecidas aos usuários como sugestões para a melhor maneira de realizar tarefas ou para os usuários a considerar e serão comunicadas através de notas no protocolo. Uma pasta de entrada de dados de exemplo e saída analítica está publicamente disponível para os usuários, e é organizada como sugerido no protocolo ( de novo </Em&gt; transcriptoma montagem e análise. 1. Configure o projeto, faça o upload de seqüências de leitura em sequência e avalie leituras usando FastQC <ol><li> Obtenha acesso à atmosfera e ao ambiente de descoberta. <ol><li> Solicite uma conta CyVerse gratuita navegando para a página de registro ( por exemplo, person@institution.edu). </li><li> Preencha as informações necessárias e envie. </li><li> Navegue até a página principal (http://www.cyverse.org/) e selecione &quot;Iniciar sessão&quot; na barra de ferramentas superior. Selecione &quot;Cyverse Login&quot; e faça login usando suas credenciais do CyVerse. </li><li> Navegue até a guia Aplicativos e Serviços e solicite acesso ao Atmosphere. O acesso ao Ambiente de Descoberta é automaticamente concedido. </li></ol></li><li> Configurar o projeto e mover dados para o armazenamento de dados. <ol><li> Efetue login no ambiente de descoberta (https://de.iplantcollaborative.org/de). Selecione a guia &quot;Dados&quot; para abrir um menu contendo todas as pastas no Data Store. </Li <li> Crie uma pasta de projeto principal que abrigará todos os dados associados ao projeto. Localize a barra de ferramentas na parte superior da janela de dados e selecione Arquivo | Nova pasta. Não use espaços ou caracteres especiais nos nomes das pastas ou quaisquer nomes de arquivos de entrada / saída, por exemplo &quot;! @ # () [] {}:; $% ^ &amp; *.&quot; Em vez disso, use sublinhados ou traços, ou seja , &quot;_&quot; ou &quot;-&quot;, quando apropriado. </li><li> Crie cinco pastas dentro da pasta principal do projeto para organizar as análises ( Figura 1 ) Nomeie as pastas da seguinte maneira sem vírgulas ou aspas: &quot;1_Raw_Sequence&quot;, &quot;2_High_Quality_Sequence&quot;, &quot;3_Assembly&quot;, &quot;4_Differential_Expression&quot;, &quot;5_Annotated_Assembly&quot;. As subpastas serão colocadas em cada uma das pastas principais do projeto ( Figura 2 ). </li></ol></li></ol><img alt="figura 1" src="/files/ftp_upload/55009/55009fig1.jpg" /> Figura1: Uma Visão Geral da Organização da Pasta do Projeto e o Workflow de Montagem e Análise do Novo Transcriptoma. Os usuários enviarão as leituras de seqüência em bruto para a pasta principal do projeto no Data Store e, em seguida, colocarão os resultados de cada etapa em pastas separadas. <a href="http://ecsource.jove.com/files/ftp_upload/55009/55009fig1large.jpg" target="_blank">Clique aqui para ver uma versão maior desta figura.</a> <img alt="Figura 2" src="/files/ftp_upload/55009/55009fig2.jpg" /> Figura 2: Uma visão geral detalhada do fluxo de trabalho de montagem e análise do transcriptoma De Novo que ocorre na Cyberinfrastructure CyVerse. Todo o fluxo de trabalho de montagem e análise será concluído em cinco etapas, cada uma com sua própria pasta (ícones em negrito e numerados). Cada uma das cinco pastas de etapa de fluxo de trabalho numeradas tem subpastas contendo dados de saída de análises bioinformáticas (pastaÍcones). Entradas para análise vêm de uma subpasta e, em seguida, mover para outra pasta através da saída de um programa de análise (caixas de retângulo). Os dados finais dos três primeiros passos são comparados e preparados para publicação. Em última análise, este esquema produz uma pasta de projeto principal que tem análise passo a passo para colaboradores e / ou revisores manuscritos podem rapidamente entender o fluxo de trabalho e repeti-lo usando cada arquivo, se necessário. <a href="http://ecsource.jove.com/files/ftp_upload/55009/55009fig2large.jpg" target="_blank">Clique aqui para ver uma versão maior desta figura.</a> <ol start="3"><li> Suba arquivos de seqüência FASTQ raw para a pasta &quot;1_Raw_Sequence&quot; em uma subpasta intitulada &quot;A_Raw_Reads&quot; usando um dos três métodos a seguir. <ol><li> Utilize o recurso de carregamento simples do Armazenamento de dados para navegar até a barra de ferramentas da janela Dados, clicando no botão de dados na área de trabalho principal do DE e selecione Upload | Upload Simples do Desktop. Selecione o botão ProcurarPara navegar para os arquivos de seqüenciamento FASTQ em bruto no computador local. Este método é adequado apenas para arquivos com menos de 2 GB. </li><li> Selecione o botão Carregar na parte inferior da tela para enviar o upload. Uma notificação será registrada no canto superior direito do DE no ícone de sino que o upload foi enviado. Outra notificação será registrada quando o upload estiver concluído. </li><li> Como alternativa, use o Cyberduck para transferir arquivos maiores (https://wiki.cyverse.org/wiki/x/pYcVAQ). Instale o Cyberduck e execute-o como um programa no ambiente de trabalho do computador local. </li><li> Por fim, faça o download do iCommands e instale no computador local de acordo com as instruções (https://wiki.cyverse.org/wiki/display/DS/Using+iCommands). </li></ol></li><li> Avaliar as leituras de seqüência de dados carregadas, usando o aplicativo FastQC no DE. <ol><li> Selecione o botão &quot;Aplicativos&quot; na área de trabalho principal do DE para abrir uma janela contendo todos os aplicativos de análise disponíveis no DE. </li><li> Pesquisar e abrir a vitóriaDow para a ferramenta FastQC na barra de ferramentas de pesquisa na parte superior da janela. Abra a versão de vários arquivos se houver mais de um arquivo FASTQ. Selecione Arquivo | Nova Pasta para criar uma pasta chamada &quot;B_FastQC_Raw_Reads&quot; e selecione esta pasta como a pasta de saída. </li><li> Carregue os arquivos de leitura FASTQ na janela de ferramentas chamada &quot;Selecionar dados de entrada&quot; e selecione &quot;Análise de Inicialização&quot;. </li><li> Abra o arquivo .html ou .pdf para exibir os resultados quando a análise estiver concluída. O FastQC executa várias análises que testam diferentes aspectos dos arquivos de leitura ( Figura 3 ). </li></ol></li></ol> 2. Trim e filtro de qualidade Raw lê para produzir seqüência de alta qualidade Nota: Use o aplicativo Trimmomatic ou o aplicativo Sickle. <ol><li> Procure o aplicativo Trimmomatic programável no DE e abra-o como antes. <ol><li> Faça o upload da pasta de arquivos de leitura FASTQ em bruto na seção &quot;Configurações&quot;. </li><li> Selecione se o seOs arquivos de quence são single-ou emparelhados-end. </li><li> Use o arquivo de controle padrão fornecido selecionando o botão Procurar e colando / iplant / home / shared / Trinity_transdecoder_trinotate_databases na caixa &quot;Viewing:&quot;. Selecione o arquivo denominado Trimmomaticv0.33_control_file e inicie a análise. O arquivo pode ser baixado, as configurações editadas e, em seguida, carregado na segunda pasta de projeto para criar um script de corte personalizado. </li><li> Opcional: Se a análise FastQC identificou sequências de adaptador, use a configuração ILLUMINACLIP para recortar adaptadores Illumina. Selecione o arquivo de adaptador apropriado na pasta / iplant / home / shared / Trinity_transdecoder_trinotate_databases como acima. </li></ol></li><li> Seqüência de corte de qualidade lê usando Sickle. <ol><li> Pesquise e abra o aplicativo Sickle no DE. Selecione as leituras FASTQ aparadas como leituras de entrada e renomeie os arquivos de saída. Inclua configurações de qualidade nas opções. Configurações típicas são Formato de qualidade: illumina, sanger, solexa; Qualidade tLimiar: 20; Comprimento mínimo: 50. </li><li> Mova toda a saída para a pasta aparada e filtrada (2_High_Quality_Sequence). </li></ol></li><li> Avalie as leituras finais usando FastQC e compare com relatórios FastQC anteriores. Selecione o arquivo .html para exibir uma página da Web de todos os resultados. Selecione a pasta de arquivos de imagem (.png) que são fornecidos na saída, se isso não puder ser exibido. </li></ol> 3. De Novo Transcriptome Assembly Usando a Trindade na Atmosfera <ol><li> Abra a versão mais atual da instância do Atmosphere navegando para a página wiki (https://wiki.cyverse.org/wiki/x/dgGtAQ). Selecione o link para a versão mais recente da imagem Trinity e Trinotate. Como alternativa, procure &quot;Trinotate&quot; na ferramenta de pesquisa de imagens Atmosphere (https://atmo.iplantcollaborative.org/application/images) para exibir todas as versões das imagens Trinity e Trinotate. <ol><li> Selecione o botão &quot;Log in to launch&quot; e nomeie o Atmosphere iNstance. </li><li> Selecione um tamanho de instância de &quot;medium3&quot; (CPU: 4, Mem: 32GB) ou &quot;grande3&quot; (CPU: 8, Mem: 64 GB). Inicie a instância e aguarde que ela seja compilada. Em alguns casos raros CyVerse sofre manutenção para atualizar plataformas. As instâncias existentes estão disponíveis durante essas atualizações, mas pode não ser possível criar novas instâncias. Visite a página CyVerse Status para ver o estado atual de qualquer plataforma (http://status.cyverse.org/). </li></ol></li><li> Abra a instância assim que estiver pronta clicando no nome e selecionando &quot;Remote Desktop&quot; na parte inferior do menu à direita. Permitir Java e VNC Viewer se solicitado. Selecione o botão &quot;Conectar&quot; na janela do Visualizador do VNC e selecione &quot;Continuar&quot;. <ol><li> Faça login para abrir uma janela separada que será a nova instância de computação em nuvem. </li><li> Mova os arquivos de leitura FASTQ cortados e / ou filtrados para a instância usando um dos três métodos descritos nas etapas 1.3.1 a 1.3.4. NosE o navegador da Internet para acessar o DE e baixar arquivos exatamente como antes no computador local. Ou use iCommands instalados nessas imagens para transferir rapidamente grandes conjuntos de dados. </li></ol></li><li> Running Trinity para montar leituras de alta qualidade. <ol><li> Configure a pasta de análise na instância Atmosphere. Use o script disponível no DE (/ iplant / home / shared / Trinity_transdecoder_trinotate_databases) ou copie e cole os comandos da página wiki (https://wiki.cyverse.org/wiki/x/dgGtAQ). A explicação de todos os comandos pode ser encontrada na página wiki. </li><li> Uma vez estabelecida a pasta de análise e os bancos de dados do Trinotate, execute o assemblage Trinity usando os comandos acima. Existem vários arquivos de saída, mas o mais importante é o arquivo de montagem final intitulado &quot;Trinity.fasta&quot;. Renomeie este arquivo FASTA para ser exclusivo para o organismo e tratamento das leituras montadas antes de movê-lo para o Data Store (pasta 3_Assembly) para minimizar a confusão potencial. NOTA: A saída conta as tabelas para análise de expressão gênica diferencial em uma pasta (4_Differential_Expression). </li></ol></li><li> Avaliar a montagem usando rnaQUAST ( Figura 4 ). <ol><li> Mova os arquivos de saída da Trinity para a pasta &quot;3_Assembly&quot; no DE e marque a pasta &quot;A_Trinity_de_novo_assembly&quot;. Dê a cada transcriptome que foi montado uma subpasta dentro da pasta &quot;A_Trinity_de_novo_assembly&quot; com nomes exclusivos, incluindo o nome científico dos organismos e tratamentos associados a cada transcriptome. Crie outra subpasta chamada &quot;B_rnaQUAST_Output&quot; na &quot;pasta 3_Assembly&quot;. </li><li> Abra o aplicativo intitulado &quot;rnaQUAST 1.2.0 (baseado em denovo)&quot; e nomeie a análise e selecione &quot;B_rnaQUAST_Output&quot; como a pasta de saída. <ol><li> Adicione o novo arquivo de montagem FASTA arquivos para a seção &quot;Entrada de dados&quot;. Na seção &quot;Saída de dados&quot;, digite um nome exclusivo para o </eM&gt; assembly. Isso criará uma pasta de arquivos de saída rnaQUAST dentro da pasta &quot;B_rnaQUAST_Output&quot;. </li></ol></li><li> Selecione opções adicionais nas seções &quot;GenemarkS-T Gene Prediction&quot;, &quot;BUSCO&quot; e &quot;Parameters&quot;. <ol><li> Selecione prokaryote na seção &quot;GenemarkS-T Gene Prediction&quot; se o organismo não é eucariótico. </li><li> Execute o BUSCO para selecionar o botão de navegação e copie o caminho iplant / home / shared / iplantcollaborative / example_data / BUSCO.sample.data na caixa &quot;Viewing:&quot; e pressione enter. Selecione a pasta BUSCO mais específica que está disponível para o organismo. NOTA: BUSCO irá avaliar a montagem para a linhagem de genes específicos do núcleo, e saída que a percentagem de genes núcleo é encontrado. Existem pastas gerais, por exemplo , eucariotas e linhagens mais específicas, por exemplo , artrópodes. </li></ol></li></ol></li><li> Procure por &quot;Transcript decoder&quot; e execute Transdecoder no de novO Trinity assembly arquivo FASTA de saída no ambiente de descoberta. </li><li> Mova o arquivo .pep de saída para a pasta de novo assembly (3_Assembly) para uso na etapa 5 Annotation. </li></ol> 4. Expressão Diferencial Pairwise Usando DESeq2 no DE <ol><li> Abra o aplicativo DESeq2 no DE, como descrito anteriormente. Nomeie a análise e selecione a pasta de saída como 4_Differential_Expression. </li><li> Na seção &quot;Entradas&quot;, selecione o arquivo de tabela de contagens da execução de montagem Trinity ea coluna que os nomes contig podem ser encontrados na tabela de contagem. </li><li> Insira os cabeçalhos de coluna do arquivo de tabela de dados de contagens para determinar quais colunas são comparadas. Inclua as vírgulas entre cada uma das condições. Não inclua o cabeçalho da primeira coluna que contém os nomes contig. </li><li> Para repetições, repita o mesmo nome ( por exemplo , Tratamento1rep1, Tratamento1rep2, Tratamento1rep3 se tornaria Tratamento1, Tratamento1, Tratamento1). EmE segunda linha, fornecer os nomes das duas condições a serem comparados ( por exemplo , Tratamento1, Tratamento2). Faça corresponder os nomes de cabeçalho de coluna fornecidos na primeira linha. Observação: esses cabeçalhos de coluna devem ser alfanuméricos e não podem conter caracteres especiais. </li></ol> 5. Anotações Usando Trinotate <ol><li> Execute cada parte do Trinotate na instância de computação em nuvem do Atmosphere. Nota: Os comandos Bash são fornecidos em um arquivo txt para serem copiados, colados e modificados antes de serem executados no DE (/ iplant / home / shared / Trinity_transdecoder_trinotate_databases) ou na página wiki (https://wiki.cyverse.org/ Wiki / x / dgGtAQ). Se anotar vários assemblies, anotar cada assembly um de cada vez e, em seguida, transferir arquivos de anotações concluídas de volta para a pasta &quot;5_Annotation&quot; cada um com uma pasta única correspondente com o nome do assembly. <ol><li> Execute o comando bash para pesquisar transcrições da Trindade. Alterar o número de threads para coincidir com quantos CPUs sãoNa instância, ou seja , médio tem 4 CPUs e grande tem 8 CPUs. Consulte a etapa 3.1.2 para obter mais detalhes. Altere o comando Trinity.fasta para coincidir com o nome do arquivo de montagem FASTA. OBSERVAÇÃO: as pesquisas BLAST + exigirão mais tempo. Pode ser dias antes de ser concluída. A atividade do computador em nuvem pode ser verificada em Atmosfera sem ter que abrir o Visualizador do VNC. </li><li> Execute o comando bash para pesquisar as proteínas preditas pelo Transdecoder. Como antes, altere o número de threads eo nome do arquivo para que correspondam às condições em 5.2.1. </li><li> Execute o comando bash para HMMER e altere o número de threads como acima. </li><li> Execute o comando bash para signalP e tmHMM, se necessário. SignalP irá prever peptídeos de sinal e tmHMM prevê motivos de proteína transmembrana. </li></ol></li><li> Carregando resultados no banco de dados SQLite <ol><li> Uma vez que todas as análises acima são concluídas, execute o comando bash para carregar os arquivos de saída em um banco de dados de anotação SQLite final. Remover quaisquer comandosPara análises que não foram executadas. </li><li> Exporte o banco de dados SQLite para um arquivo .xls para visualização em visualizadores de tabelas populares. </li></ol></li></ol>

<ol><li>Hasselmann, M., Ferretti, L., Zayed, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Beyond+fruit-flies%3A+population+genomic+advances+in+non-Drosophila+arthropods%5BTitle%5D)+AND+%22Brief.+Funct.+Genomics%22%5BJournal%5D)">Beyond fruit-flies: population genomic advances in non-Drosophila arthropods.</a> Brief. Funct. Genomics. 14 (6), 424-431 (2015).</li>
<li>Scholz, M. B., Lo, C. -C., Chain, P. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Next+generation+sequencing+and+bioinformatic+bottlenecks%3A+the+current+state+of+metagenomic+data+analysis%5BTitle%5D)+AND+%22Anal.+Biotech%22%5BJournal%5D)">Next generation sequencing and bioinformatic bottlenecks: the current state of metagenomic data analysis.</a> Anal. Biotech. 23 (1), 9-15 (2012).</li>
<li>Merchant, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+iPlant+Collaborative%3A+Cyberinfrastructure+for+Enabling+Data+to+Discovery+for+the+Life+Sciences%5BTitle%5D)+AND+%22PLoS+Biol%22%5BJournal%5D)">The iPlant Collaborative: Cyberinfrastructure for Enabling Data to Discovery for the Life Sciences.</a> PLoS Biol. 14 (1), e1002342(2016).</li>
<li>Oliver, S. L., Lenards, A. J., Barthelson, R. A., Merchant, N., McKay, S. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Using+the+iPlant+collaborative+discovery+environment%5BTitle%5D)+AND+%22Cur.+Protoc.+Bioinformatics%22%5BJournal%5D)">Using the iPlant collaborative discovery environment.</a> Cur. Protoc. Bioinformatics. , 1-22 (2013).</li>
<li>Skidmore, E., Kim, S., Kuchimanchi, S., Singaram, S., Merchant, N., Stanzione, D. <a target="_blank" href="https://www.researchgate.net/publication/254003391_iPlant_atmosphere_A_gateway_to_cloud_infrastructure_for_the_plant_sciences">iPlant atmosphere: a gateway to cloud infrastructure for the plant sciences.</a> Proc. 2011 ACM. , 59-64 (2011).</li>
<li>Altschul, S. F., Gish, W., Miller, W., Myers, E. W., Lipman, D. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Basic+local+alignment+search+tool%5BTitle%5D)+AND+%22J.+Mol.+Bio%22%5BJournal%5D)">Basic local alignment search tool.</a> J. Mol. Bio. 215 (3), 403-410 (1990).</li>
<li>Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., Zdobnov, E. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((BUSCO%3A+assessing+genome+assembly+and+annotation+completeness+with+single-copy+orthologs%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs.</a> Bioinformatics. , (2015).</li>
<li>Camacho, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((BLAST%2B%3A+architecture+and+applications%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">BLAST+: architecture and applications.</a> BMC Bioinformatics. 10, 421(2009).</li>
<li>Eddy, S. R. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/9918945">Profile hidden Markov models.</a> Bioinformatics. 14 (9), 755-763 (1998).</li>
<li>Krogh, A., Larsson, B., von Heijne, G., Sonnhammer, E. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Predicting+transmembrane+protein+topology+with+a+hidden+markov+model%3A+application+to+complete+genomes%5BTitle%5D)+AND+%22J.+Mol.+Biol%22%5BJournal%5D)">Predicting transmembrane protein topology with a hidden markov model: application to complete genomes.</a> J. Mol. Biol. 305 (3), 567-580 (2001).</li>
<li>Finn, R. D., Coggill, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Pfam+protein+families+database%3A+towards+a+more+sustainable+future%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">The Pfam protein families database: towards a more sustainable future.</a> Nucleic Acids Res. 44 (D1), D279-D285 (2016).</li>
<li>Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., Tanabe, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((KEGG+as+a+reference+resource+for+gene+and+protein+annotation%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">KEGG as a reference resource for gene and protein annotation.</a> Nucleic Acids Res. 44 (D1), D457-D462 (2016).</li>
<li>Kanehisa, M., Goto, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((KEGG%3A+Kyoto+Encyclopedia+of+Genes+and+Genomes%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">KEGG: Kyoto Encyclopedia of Genes and Genomes.</a> Nucleic Acids Res. 28 (1), 27-30 (2000).</li>
<li>Wolfien, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((TRAPLINE%3A+a+standardized+and+automated+pipeline+for+RNA+sequencing+data+analysis%2C+evaluation+and+annotation%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">TRAPLINE: a standardized and automated pipeline for RNA sequencing data analysis, evaluation and annotation.</a> BMC Bioinformatics. 17, 21(2016).</li>
<li>Bolger, A. M., Lohse, M., Usadel, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Trimmomatic%3A+a+flexible+trimmer+for+Illumina+sequence+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Trimmomatic: a flexible trimmer for Illumina sequence data.</a> Bioinformatics. 30 (15), 2114-2120 (2014).</li>
<li>Robinson, M. D., McCarthy, D. J., Smyth, G. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((edgeR%3A+a+Bioconductor+package+for+differential+expression+analysis+of+digital+gene+expression+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.</a> Bioinformatics. 26 (1), 139-140 (2010).</li>
<li>Anders, S. <a target="_blank" href="https://pdfs.semanticscholar.org/09a4/33297b5c0476a77ec4a34bba145f9149c327.pdf">Analysing RNA-Seq data with the DESeq package.</a> Mol. Biol. 43 (4), 1-17 (2010).</li>
<li>Love, M. I., Huber, W., Anders, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Moderated+estimation+of+fold+change+and+dispersion+for+RNA-seq+data+with+DESeq2%5BTitle%5D)+AND+%22Genome+Bio%22%5BJournal%5D)">Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.</a> Genome Bio. 15 (12), 1-21 (2014).</li>
</ol>

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

Existem cinco passos críticos no protocolo que cada um criará sua própria pasta separada dentro da pasta principal do projeto ( Figuras 1 e 2 ). Todos os dados de seqüenciamento primário são sacrossanct: deve ser carregado e mantido na primeira pasta rotulada &quot;1_Raw_Sequence&quot; e não alterado de forma alguma. Os dados podem ser carregados de uma das três maneiras. A interface DE pode ser usada para carregar arquivos diretamente. Esta é a maneira mais fácil de fazer o upload de dados, mas também levará mais tempo para transferir. Cyberduck tem uma interface gráfica e permite aos usuários arrastar e soltar arquivos para transferir para o DE. O iCommands é uma ferramenta de linha de comando que pode ser usada para transferir dados de e para o Data Store, fazer diretórios e gerenciar conjuntos de dados, e é provavelmente a maneira mais rápida de transferir arquivos de dados. Todos os dados no Data Store podem ser compartilhados com outros usuários CyVerse (https://wiki.cyverse.org/wiki/display/DEmanual/Sharing+Data+Files+and+Folders+Via+the+DiscoveRy + Environment), tornados públicos através de um URL gerado (https://wiki.cyverse.org/wiki/display/DEmanual/Sharing+Data+Files+Via+Public+Links), ou podem ser hospedados publicamente e anonimamente ( Nenhum nome de usuário requerido) dados disponíveis da comunidade (http://data.iplantcollaborative.org; http://mirrors.cyverse.org). Dentro dessa pasta, as leituras de sequência em bruto são analisadas com FastQC (http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/) para avaliar como aparar e filtrar as leituras para gerar leituras de alta qualidade. Depois de aparar e filtrar a qualidade, é útil comparar as saídas FastQC para determinar se a qualidade de leitura foi alterada para determinar se melhorou sem perder informações ( Figura 3 ). Observe que o eixo x de FastQC não é linear, mas é binned para muitos gráficos de saída, o que pode levar a interpretação errada dos resultados. As leituras aparadas e filtradas são então usadas para montar transcriptomes de novo usando uma instância de computação em nuvem Atmosphere. esteO computador da nuvem usa a tela, o teclado eo mouse do computador local, mas tem seu próprio software (Trinity e Trinotate) e hardware instalado. A execução de programas na instância do computador em nuvem não afetará o computador local de qualquer forma. A montagem de novo ea anotação a jusante provavelmente serão as duas etapas mais longas neste fluxo de trabalho. Portanto, eles são concluídos na atmosfera para evitar comum laboratório compartilhado problemas do computador que iria interromper a análise como interrupções de energia, reinicia após atualizações automáticas tarde da noite, ou falhas causadas por outros usuários. A anotação de trinotate usa BLAST + 8 , HMMER 9 , tmHMM 10 e PFAM 11 . A saída final da anotação é um banco de dados SQLite e um arquivo .xls. As saídas podem ser usadas fora CyVerse em plataformas de análise downstream, como KEGG 12 , 13 . Este fluxo de trabalhoEstá pronto para usar no DE e Atmosfera. Isso elimina a necessidade de gastar tempo instalando, configurando e solucionando problemas de cada pacote de análise e todas as dependências que cada ferramenta requer. Isso agiliza as análises dos pesquisadores, minimiza o esforço desperdiçado e reduz a barreira de entrada para muitos cientistas. Este fluxo de trabalho especificamente monta leituras single-ou pares-end da plataforma de sequenciamento Illumina, mas muitas ferramentas existem no DE e Atmosfera para lidar com outros tipos de tecnologias de seqüenciamento. Ferramentas neste fluxo de trabalho podem ser facilmente substituídas por uma ferramenta alternativa correspondente para lidar com qualquer tipo de tecnologia de sequenciamento de entrada. Isso também é verdade para novas versões de ferramentas de análise ou ferramentas completamente novas. Este fluxo de trabalho é projetado especificamente para montar, comparar e anotar apenas alguns transcriptomes de cada vez. Portanto, os usuários podem achar que é demorado para montar múltiplos transcriptomes para genética populacional comparativa. AnálisePipelines estará disponível para os usuários de genética populacional em um futuro próximo eo link para o pipeline pode ser encontrado na página wiki (https://wiki.cyverse.org/wiki/x/dgGtAQ). O passo de análise de expressão genética diferencial pode lidar com repetições, mas é uma comparação em pares e não avaliará com precisão múltiplos fatores ( por exemplo , condições que variam ao longo do tempo, mais de dois tratamentos). Existem fluxos de trabalho automatizados para organismos com genomas de referência ( por exemplo , TRAPLINE 14 ). Enquanto fluxos de trabalho automatizados são os mais fáceis de usar para iniciantes, conjuntos de novo exigem avaliação e consideração para cada etapa delineada aqui. Além disso, os usuários são obrigados a usar pipelines automatizados como eles são construídos e, portanto, são inerentemente não flexíveis para atender às demandas de mudança dos usuários. Como a maior parte deste protocolo é realizado através da Internet, os usuários podem ter problemas com as configurações do navegador. Primeiramente,Os bloqueadores de pop-up podem impedir que as janelas se abram, ou podem impedir que as janelas sejam abertas até que seja dada permissão ao CyVerse no navegador. O Atmosphere usa o VNC para acessar desktops remotos, mas outros softwares podem ser usados. Todo esse protocolo foi conduzido no Firefox versão 45.0.2 e deve funcionar com todos os navegadores de Internet populares, mas algumas inconsistências podem aparecer. O fluxo de trabalho será atualizado quando a Trinity lançar novas versões (https://github.com/trinityrnaseq/trinityrnaseq/wiki). As versões mais recentes e as informações atualizadas sobre o fluxo de trabalho podem ser encontradas na página do tutorial wiki ( Tabela 1 , https://wiki.cyverse.org/wiki/x/dgGtAQ). Os usuários podem contatar o suporte diretamente ou postar perguntas no Ask CyVerse (ask.cyverse.org/) para solucionar problemas com o fluxo de trabalho. No DE existem várias aplicações para realizar cada passo deste protocolo. Por exemplo, os usuários podem querer executar Scythe (https://github.com/najoshi/sickle) em vez de Trimmomatic15 para ler recortar ou executar EdgeR 16 em vez de DESeq 17 , 18 . Embora fora do escopo deste manuscrito, os aplicativos DE podem ser copiados, editados e liberados pelos usuários (https://wiki.cyverse.org/wiki/display/DEmanual/Creating,+Copying,+and+Editing+DE+ Apps) ou novos aplicativos podem ser adicionados pelos usuários (https://wiki.cyverse.org/wiki/display/DEmanual/Dockerizing+Your+Tools+for+the+CyVerse+Discovery+Environment). As imagens do Atmosphere também podem ser modificadas e reimaginadas para criar fluxos de trabalho novos ou modificados que atendam às necessidades dos usuários mais especificamente (https://wiki.cyverse.org/wiki/x/TwHX). Este trabalho serve como uma introdução à utilização da linha de comando para mover dados e executar análises. Os usuários podem considerar a utilização de recursos de linha de comando mais avançados, como interfaces de programação de aplicativos CyVerse (APIs) (http://www.cyverse.org/science-apis), ou projetar seus próprios aplicativos DE, que exigem conhecimentoSobre como a ferramenta de análise é executada na linha de comando (https://wiki.cyverse.org/wiki/display/DEmanual/Creating+a+New+App+Interface).

O Homo sapiens e várias espécies-chave de animais como Drosophila melanogaster , Mus musculus e Danio rerio representam a maioria do trabalho genômico funcional atual e passado. No entanto, o custo rapidamente decrescente da tecnologia de sequenciamento de alto rendimento está proporcionando oportunidades para a genômica funcional em espécies animais não-modelo ( também conhecidas como &quot;negligenciadas&quot; ou &quot;subatendidas&quot;) 1 . Esta é uma transição importante na genômica, uma vez que os organismos não-modelo freqüentemente representam espécies economicamente relevantes ( por exemplo , ostras, camarão, caranguejo) e oferecem oportunidades para investigar novos fenótipos e sistemas biológicos fora do escopo daqueles encontrados em espécies-modelo. Embora os organismos subatendidos apresentam uma oportunidade atraente para investigar sistemas biológicos únicos, vários desafios enfrentam os pesquisadores particularmente durante a análise bioinformática. Alguns dosOs desafios são inatos ao processamento de grandes conjuntos de dados, enquanto outros resultam da falta de recursos genéticos disponíveis para os pesquisadores que trabalham em organismos subatendidos, como um genoma de referência, ontologias específicas de organismos, etc. Os desafios do isolamento e seqüenciamento de ácidos nucleicos são freqüentemente rotineiros. Comparação com os da análise de dados, e como tais análises bioinformáticas geralmente revela-se o custo mais subestimado de projetos de seqüenciamento 2 . Por exemplo, uma análise bioinformática de sequenciação básica de próxima geração pode consistir nas seguintes etapas: filtragem de qualidade e corte de leituras em seqüência em bruto, montagem de leituras curtas em peças contíguas maiores e anotação e / ou comparações com outros sistemas para obter entendimento biológico. Embora aparentemente simples, este fluxo de trabalho de exemplo requer conhecimento especializado e recursos computacionais além do escopo de um computador de laboratório, colocando-o fora do alcance de muitos cientistas estudando não-Organismos-modelo. Os desafios inerentes podem ser baseados na infra-estrutura ou no conhecimento. Um desafio clássico da infra-estrutura é o acesso a recursos computacionais apropriados. Por exemplo, a montagem ea anotação dependem de algoritmos computacionalmente intensivos que exigem computadores poderosos ou clusters de computadores, com grande quantidade de RAM (256 GB-1 TB) e vários processadores / núcleos para serem executados. Infelizmente, muitos pesquisadores não têm acesso a tais recursos de computação ou não têm o conhecimento necessário para interagir com esses sistemas. Outros pesquisadores podem ter acesso a clusters de computação de alto desempenho por meio de suas universidades ou instituições, mas o acesso a esses recursos pode ser limitado e às vezes resulta em encargos por hora computacional, ou seja , o número de processadores de CPU multiplicado pelo número de clock Horas &quot;que esses processadores estão em execução. Aproveitando um sistema de infra-estrutura cibernética financiado pela US National Science Foundation sComo o CyVerse 3, que oferece acesso gratuito aos recursos de computação para pesquisadores, nos Estados Unidos e em todo o mundo, pode ajudar a aliviar os desafios de infra-estrutura, como será demonstrado aqui. Um exemplo de um desafio típico baseado no conhecimento é entender o software necessário para análises completas. Para efetivamente realizar um projeto baseado em seqüenciamento, os pesquisadores precisam estar familiarizados com a miríade de ferramentas de software que foram desenvolvidas para análises bioinformáticas. Aprender cada pacote é difícil por si só, mas é exacerbado pelo fato de que os pacotes estão constantemente sendo atualizados, rereleased, reunidos em novos fluxos de trabalho, e às vezes ficam restritos para uso sob novas licenças. Além disso, ao vincular as entradas e saídas dessas ferramentas, às vezes, é necessário transformar os tipos de dados para torná-los compatíveis, adicionando outra ferramenta ao fluxo de trabalho. Finalmente, também é difícil saber qual é o pacote deE melhor &quot;para uma análise, e freqüentemente identificar o melhor software para determinadas condições experimentais é uma questão de sutis diferenças. Em alguns casos, as revisões úteis do software estão disponíveis, mas devido à liberação contínua de atualizações novas e de opções do software, estas rapidamente saem da data. Para os pesquisadores que investigam organismos subatendidos, esses desafios inatos vêm além dos desafios associados à análise de dados em um novo organismo. Esses desafios específicos do organismo mal atendidos são melhor ilustrados durante a anotação do gene. Por exemplo, os organismos subatendidos frequentemente não têm um organismo modelo estreitamente relacionado que possa ser razoavelmente utilizado para identificar a ortografia e a função dos genes ( por exemplo, invertebrados marinhos e Drosophila ). Muitas ferramentas bioinformáticas também exigem &quot;treinamento&quot; para identificar motivos estruturais, que podem ser usados ​​para identificar a função do gene. No entanto, os dados de treinamento normalmente só estão disponíveis para modOs organismos e formação de modelos de Markov ocultos (HMMs) está fora do âmbito de biólogos, e até mesmo muitos bioinformáticos. Por fim, mesmo se as anotações podem ser realizadas usando dados de organismos modelo, algumas ontologias de genes associadas a organismos modelo não fazem sentido quando a biologia ea história natural do organismo subatendido é considerada ( por exemplo , transferindo informações de Drosophila para camarão). À luz desses desafios, os recursos bioinformáticos precisam ser desenvolvidos com pesquisadores que realizam análises de novo sobre os organismos carentes especificamente em mente. Os próximos anos de projetos de sequenciação de genômica funcional ajudarão a reduzir o hiato entre os organismos modelo e subatendidos ( <a href="https://genome10k.soe.ucsc.edu/">https://genome10k.soe.ucsc.edu/</a> ), mas há muitas ferramentas que precisarão ser desenvolvidas para enfrentar os desafios Consideradas acima. CyVerse é dedicado à criação de ecossistemas de iA interoperabilidade ligando a infra-estrutura existente e as aplicações de terceiros para fornecer gerenciamento de dados, ferramentas de análise bioinformática e visualizações de dados para cientistas de vida. A interoperabilidade ajuda a suavizar as transições entre aplicações e plataformas bioinformáticas fornecendo recursos computacionais escaláveis ​​e limitando as conversões de formatos de arquivos ea quantidade de dados transferidos entre plataformas. CyVerse oferece várias plataformas, incluindo o Discovery Environment (DE 4 , Atmosphere 5 e Data Store 3. O DE é baseado na web e tem muitas ferramentas analíticas bioinformáticas comuns convertidas em formatos user-friendly point-and-click &quot;), E é a interface gráfica do usuário (GUI) para o armazenamento de dados onde os grandes conjuntos de dados ( ou seja , leituras de seqüenciamento em bruto, genomas montados) são armazenados e gerenciados.A atmosfera é um serviço de computação em nuvem que oferece aos pesquisadores maior flexibilidade paraUsando recursos computacionais da Máquina Virtual, que possuem uma ampla gama de ferramentas bioinformáticas pré-instaladas. Ambas as plataformas estão vinculadas ao Data Store e podem ser usadas em conjunto para criar fluxos de trabalho como o descrito aqui. Este relatório centra-se em um novo transcriptome assembly e análise diferencial expressão genética workflows, e ainda aborda algumas das melhores práticas associadas ao desenvolvimento e realização de bioinformática análises. Uma explicação da missão mais ampla do CyVerse ( <a href="http://www.cyverse.org/about">http://www.cyverse.org/about</a> ) e descrições detalhadas da plataforma ( <a href="http://www.cyverse.org/learning-center">http://www.cyverse.org/learning-center</a> ) estão publicamente disponíveis. Todas as análises aqui descritas utilizam o Ambiente de Descoberta 4 (DE) e Atmosfera 5 , e são apresentadas de forma a torná-las acessíveis a investigadores de todos os níveis computacionais. DE workflows e AtmosphEre imagens podem ser referenciadas diretamente usando URLs para garantir proveniência a longo prazo, reutilização e reprodutibilidade.

<table><tbody><tr><td>Trimmomatic v0.33</td><td>USADELLAB.org</td><td>https://github.com/timflutre/trimmomatic</td><td>https://de.iplantcollaborative.org/de/?type=apps&app-id=9c2a30dc-028d-11e6-a915-ab4311791e69</td></tr><tr><td>Sickle</td><td>Joshi e Fass</td><td>https://github.com/najoshi/sickle</td><td>https://de.iplantcollaborative.org/de/?type=apps&app-id=68b278f8-d4d6-414d-9a64-b685a7714f7c</td></tr><tr><td>Trinity</td><td>Instituto Broad e Universidade Hebraica de Jersalem</td><td>https://github.com/trinityrnaseq/trinityrnaseq/wiki</td><td>https://atmo.iplantcollaborative.org/application/images/1261</td></tr><tr><td>rnaQUAST v1.2.0</td><td>Laboratório de Biologia Algorítmica, Universidade Acadêmica de São Petersburgo da Academia Russa de Ciências</td><td>http://spades.bioinf.spbau.ru/rnaquast/release1.2.0/manual.html</td><td>https://de.iplantcollaborative.org /de/?type=apps& app- id=980dd11a-1666-11e6-9122- 930ba8f23352</td></tr><tr><td>Transdecoder</td>, <td>Broad Institute e Commonwealth Scientific and Industrial Research Organization</td><td>, https://transdecoder.github.io</td><td></td></tr><tr><td>https://de.iplantcollaborative.org/de/?type=apps&app-id=5a0ba87e-b0fa-4994-92a2-0d48ee881179 EdgeR</td><td>, Robinson et al., 2010.</td><td>https://bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeR.pdf</td><td>https://de.iplantcollaborative.org/de/?type=apps&app-id=5aa9e294-6f95-42f9-98e9-c9c96b44f499</td></tr><tr><td>Trinotate</td><td>Broad Institute e a Universidade Hebraica de Jersalem</td><td>https://trinotate.github.io/</td> <td>https://atmo.iplantcollaborative.org/application/images/1261</td></tr></tbody></table>

aproveitando os recursos cyverse para<em&gt; de novo</em&gt; transcriptômica comparativa de organismos não servidos (não modelos)

Leveraging CyVerse Resources for De Novo Comparative Transcriptomics of Underserved (Non-model) Organisms

Este fluxo de trabalho permite que pesquisadores novatos aproveitem recursos computacionais avançados, como computação em nuvem, para realizar transcriptômicos comparativos em pares. Ele também serve como um primário para biólogos para desenvolver dados cientista competências computacionais, por exemplo , executar comandos bash, visualização e gestão de grandes conjuntos de dados. Todos os códigos de linha de comando e explicações adicionais de cada comando ou etapa podem ser <a href="https://wiki.cyverse.org/wiki/x/dgGtAQ">encontrados</a> no wiki ( <a href="https://wiki.cyverse.org/wiki/x/dgGtAQ">https://wiki.cyverse.org/wiki/x/dgGtAQ</a> ). As plataformas Discovery Environment e Atmosphere são conectadas entre si através do CyVerse Data Store. Como tal, uma vez que os dados iniciais de sequenciação em bruto foram carregados, não há mais necessidade de transferir arquivos de dados grandes por meio de uma conexão com a Internet, minimizando o tempo necessário para realizar análises. Este protocolo foi concebido para analisar apenas dois tratamentos ou condições experimentais. A análise diferencial da expressão genética éOugh pairwise comparações, e não será adequado para testar vários fatores. Este fluxo de trabalho também é projetado para ser manual, em vez de automatizado. Cada passo deve ser executado e investigado pelo usuário, resultando em um melhor entendimento dos dados e resultados analíticos e, portanto, melhores resultados para o usuário. Uma vez completado, este protocolo produzirá transcriptome (s) montado (s) de novo para organismos não servidos (sem modelo) sem a necessidade de mapear para genomas de referência previamente montados (que normalmente não estão disponíveis em organismos subatendidos). Estes transcriptomas de novo são ainda utilizados na análise de expressão de genes diferenciais em pares para investigar genes que diferem entre duas condições experimentais. Os genes expressos diferencialmente são então funcionalmente anotados para compreender os organismos de resposta genética têm condições experimentais. No total, os dados derivados deste protocolo são utilizados para testar hipóteses sobre respostas biológicas de organismos mal atendidos.

Uma vez que os arquivos de organização do projeto foram criados ( Figura 1 e 2 ), a primeira tarefa neste fluxo de trabalho é avaliar os arquivos de seqüência em bruto e, em seguida, limpá-los por recorte e filtragem de qualidade. O FastQC irá gerar estatísticas de resumo legível sobre as pontuações de qualidade eo comprimento das sequências do formato de arquivo FASTQ. As figuras do FastQC são então comparadas antes e depois do corte para avaliar se as leituras finais são de alta qualidade e, portanto, adequadas para montagem. &quot;Por qualidade de sequência de base&quot; mostra a qualidade média de leituras através de cada par de bases de sequenciação. É melhor ter um índice de qualidade de phred acima de 20-28 indicado pelas cores nos números de FastQC. &quot;Por pontuação de qualidade de seqüência&quot; determina se a filtragem de qualidade de leituras pode ser necessária. Se muitas leituras tiverem uma pontuação média abaixo de 20-25, então pode ser necessário filtrar com base na qualidade de leitura média. &quot;Por conteúdo de seqüência de base&quot; deve mostrar uma distribuição uniforme em todas as quatro bases nucleotídicas. Se houver viés no conteúdo de nucleótidos é mostrado, então as extremidades de corte podem ser necessárias. &quot;Por conteúdo de GC de base também deve ser uniforme em todas as posições.Se houver um oscilação as leituras podem precisar ser aparadas como em 1.4.4.3.&quot; Por conteúdo de GC seqüência &quot;deve ser uma distribuição normal.Adaptador ou reação em cadeia da polimerase (PCR ) Pode contaminação na biblioteca de seqüenciamento e enviesar a distribuição normal.Em este caso, ajuste de adaptador pode ser necessário.Resolução de comprimento de seqüência &quot;dá a média de comprimentos de todas as leituras .Reads menor do que 35-45 pares de bases são geralmente filtradas. Os &quot;níveis de duplicação de sequências&quot; mostram quantas vezes a seqüência de uma determinada leitura é vista dentro da biblioteca.Altamente duplicada ler seqüência e contagem são fornecidos na seção &quot;Overrepresented seqüências&quot; FastQC também tenta identificar se as leituras duplicadasSï¿½ sequï¿½cias adaptadoras ou outras sequï¿½cias conhecidas associadas a plataformas de sequenciaï¿½o. Um rótulo de &quot;No Hit&quot; significa que a sequência deve ser investigada utilizando NCBI BLAST 6 para determinar se é uma sequência biologicamente relevante, ou se deve ser removida. O DE também tem várias versões do BLAST disponíveis. O aplicativo DE BLASTn está disponível em: <a href="https://de.iplantcollaborative.org/de/?type=apps&amp;amp;app-id=6f94cc92-6d28-45c6-aef1-036be697671d">https://de.iplantcollaborative.org/de/?type=apps&amp;amp;app-id=6f94cc92-6d28-45c6-aef1-036be697671d</a> . Depois de seqüenciamento em bruto foram selecionados para produzir leituras de alta qualidade, as leituras precisam ser montadas para criar sequências contíguas (contigs). Em suma, as montagens são criadas alinhando todas as leituras de seqüência curta para encontrar seqüências semelhantes. Áreas de seqüência similar maiores do que um certo comprimento são consideradas as saMe porque a probabilidade de uma seqüência similar de ocorrência aleatória de um certo comprimento é quase zero. Trinity irá produzir arquivos de log, arquivos fasta para cada etapa no processo de montagem. No entanto, a saída mais importante é o arquivo de montagem final que contém os contigs, que é rotulado como &quot;Trinity.fasta&quot; e encontrado na pasta principal. Este arquivo contém todos os contigs montados, e em si mesmo não é praticamente &quot;legível por humanos&quot;. Portanto, a ferramenta rnaQUAST pode ser usada para entender a montagem com mais profundidade. A ferramenta rnaQUAST produzirá figuras que permitirão aos usuários comparar montagens para determinar quais são as mais completas ( Figura 4 ). Informações adicionais sobre cada figura do rnaQUAST podem ser encontradas no wiki ( <a href="https://wiki.cyverse.org/wiki/x/fwuEAQ">https://wiki.cyverse.org/wiki/x/fwuEAQ</a> ). Se BUSCO 7 foi executado, de particular interesse é o arquivo specificity.txt que mostra o número de completo e pARC BUSCO e o número de GeneMarkS-T predições do gene em uma montagem. Os genes BUSCO são conjuntos curados de genes comuns a um grupo de organismos. Eles podem ser usados ​​para avaliar quão bem uma montagem está capturando conjuntos de genes que se espera que estejam presentes em qualquer tipo de organismo, que é baseado em clados filogenéticos. Um aplicativo BUSCO autônomo também está disponível no DE ( <a href="https://de.iplantcollaborative.org/de/?type=apps&amp;amp;app-id=112b8a52-efd8-11e5-a15c-277125fcb1b1">https://de.iplantcollaborative.org/de/?type=apps&amp;amp;app-id=112b8a52-efd8-11e5-a15c-277125fcb1b1</a> ). A análise diferencial da expressão genética identifica os transcritos que têm diferentes padrões de expressão entre tratamentos a partir de contagens simples por tabelas de transcrições reunidas. DESeq2 utiliza um modelo linear generalizado (GLM) para determinar a variação de uma média normalizada. As experiências com repetições são preferidas de modo que variações técnicas frOm sequenciamento pode ser normalizado pelo algoritmo DESeq2. A análise DESeq2 DEG produz números e um arquivo de relatório .html que contém todos os números de saída e uma descrição. Como alternativa, EdgeR pode ser usado em vez de DESeq2, e o mesmo relatório .html será gerado com visualizações EdgeR em vez disso. Os pesquisadores podem desejar executar ambos DESeq2 e EdgeR para encontrar genes expressos diferencialmente identificados por ambos os algoritmos para qualquer experiência. O Trinotate criará um arquivo .xls de saída que pode ser aberto em qualquer programa de software de planilha. Os arquivos .txt DEG e o arquivo .xls de anotação podem ser analisados ​​e visualizados em vários aplicativos downstream que existem fora da plataforma CyVerse. <img alt="Figura 3" src="/files/ftp_upload/55009/55009fig3.jpg" /> Figura 3: Relatórios FastQC de Leituras de Seqüência Bruta, Leituras Recortadas e Leis Finais Recortadas e Filtradas. Comparação sistemática de seqüenciamento lidaS após cada passo de pré-processamento. Leituras de alta qualidade são necessárias para montar transcriptomas de novo . FastQC pode ajudar os pesquisadores a entender a qualidade inicial de seus dados de seqüenciamento e controlar a eficiência com que as leituras foram pré-processadas. Os resultados do FastQC dependerão dos organismos e amostras que estão sendo sequenciados, mas a uniformidade em todas as amostras que serão comparadas a jusante é o objetivo primário do pré-processamento. Um vídeo tutorial e documentação estão disponíveis a partir dos autores e desenvolvedores do FastQC. <a href="http://ecsource.jove.com/files/ftp_upload/55009/55009fig3large.jpg" target="_blank">Clique aqui para ver uma versão maior desta figura.</a> <img alt="Figura 4" src="/files/ftp_upload/55009/55009fig4.jpg" /> Figura 4: rnaQUAST Relatórios de Três Assembléias Separadas. RnaQUAST pode ser usado para comparar vários assemblies de leitura usando o mesmo assembler, ou múltiplo a<html> Usando as mesmas leituras iniciais. RnaQUAST aproveita BUSCO para gerar estatísticas sumárias sobre montagens baseadas em genes de núcleo conhecidos presentes em clados taxonômicos. O número de desajustes por transcrição e quantos transcritos correspondem aos genes canônicos, fração correspondente, fornecem insight sobre a precisão dos montadores. As quatro últimas subparcelas apresentadas aqui fornecem estatísticas resumidas do comprimento do contig e da isoforma ea cobertura das isoformas esperadas. NAx representa a porcentagem (x) de contigs com um comprimento maior que o comprimento (bp) no eixo y. A fração montada é o maior transcrito montado individual dividido pelo seu comprimento. A fracção coberta é a percentagem de transcritos / isoformas completos montados conforme esperado pelos genes procarióticos ou eucarióticos centrais da BUSCO. Uma descrição de todos os gráficos gerados pelo rnaQUAST está disponível ( <a href="https://wiki.cyverse.org/wiki/x/fwuEAQ">https://wiki.cyverse.org/wiki/x/fwuEAQ</a> ).09 / 55009fig4large.jpg &quot;target =&quot; _ blank &quot;&gt; Clique aqui para ver uma versão ampliada desta figura. <table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always"><tbody><tr><td> Nome do aplicativo </td><td> Plataforma CyVerse </td><td> Documentação de terceiros </td><td> Documentação CyVerse </td><td> Tempo de execução estimado para conjunto de dados de amostra </td><td> Link para o aplicativo </td></tr><tr><td> FastQC </td><td> DE </td><td> Http: //www.bioinformatics. Babraham.ac.uk/projects/fastqc/ https://www.youtube.com/watch?v=bz93ReOv87Y </td><td> Https://wiki.cyverse.org/wiki/pages/viewpage.action?pageId=9316768 </td><td> 15 min </td><td> Https: //de.iplantcollaborative. Org / de /? Type = apps &amp; app-id = 112b9aa8-c4a7-11e5-8209- 5f3310948295 </td></tr><tr><td> Trimmomatic v0.33 </td><td> DE </td><td> Https://github.com/timflutre/trimmomatic </td><td> Https://wiki.cyverse.org/wiki/display/DEapps/Trimmomatic-prOgrammable-0,33 </td><td> 30 minutos </td><td> Https: //de.iplantcollaborative. Org / de /? Type = apps &amp; app-id = 9c2a30dc-028d- 11e6-a915-ab4311791e69 </td></tr><tr><td> Foice </td><td> DE </td><td> Https://github.com/najoshi/sickle </td><td> Https://wiki.cyverse.org/wiki/display/DEapps/Sickle-quality-based-trimming </td><td> 30 minutos </td><td> Https: //de.iplantcollaborative. Org / de /? Type = apps &amp; app-id = 68b278f8-d4d6-414d-9a64-b685a7714f7c </td></tr><tr><td> Trindade </td><td> Atmosfera </td><td> Https://github.com/trinityrnaseq/trinityrnaseq/wiki </td><td> Https: //pods.iplantcollaborative Org / wiki / display / atmman / Trinity + - + Trinotate + Atmosfera + Imagem </td><td> 1 semana </td><td> Https: //atmo.iplantcollaborative. Org / application / images / 1261 </td></tr><tr><td></td><td> DE </td><td></td><td> Https://wiki.cyverse.org/wiki/display/DEapps/Trinity-64GB-2.1.1 </td><td> 2-5 dias </td><td> Https: // wikI.cyverse.org/wiki/display/DEapps/Trinity-64GB-2.1.1 </td></tr><tr><td> RnaQUAST v1.2.0 </td><td> DE, Atmosfera </td><td> Http://spades.bioinf.spbau.ru/rnaquast/release1.2.0/manual.html </td><td> Https: //pods.iplantcollaborative Org / wiki / display / TUT / rnaQUAST + 1.2.0 +% 28denovo + baseado% 29 + usando + DE </td><td> 30 minutos </td><td> Https: //de.iplantcollaborative. Org / de /? Type = apps &amp; app-id = 980dd11a-1666- 11e6-9122-930 Ba8f23352 </td></tr><tr><td> Transdecoder </td><td> DE </td><td> Https://transdecoder.github.io </td><td> Https://wiki.cyumper.org/wiki/display/DEapps/Transcript+decoder+2.0 </td><td> 2-3 horas </td><td> Https: //de.iplantcollaborative. Org / de /? Type = apps &amp; app-id = 5a0ba87e-b0fa-4994-92a2- 0d48ee881179 </td></tr><tr><td> DESeq2 </td><td> DE </td><td> Https://bioconductor.org/packages/release/bioc/html/DESeq2.html </td><td> Https: //pods.iplantcollaborative Org / wiki / páginas /Viewpage.action? PageId = 28115142 </td><td> 2-3 horas </td><td> Https: //de.iplantcollaborative. Org / de /? Type = apps &amp; app-id = 9574e87c-4f90- 11e6-a594-008 Cfa5ae621 </td></tr><tr><td> EdgeR </td><td> DE </td><td> Https://bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeR.pdf </td><td> Https://wiki.cyverse.org/wiki/pages/viewpage.action?pageId=28115144 </td><td> 2-3 horas </td><td> Https: //de.iplantcollaborative. Org / de /? Type = apps &amp; app-id = 4a08ceda-54fe- 11e6-862f-008 Cfa5ae621 </td></tr><tr><td> Trinotate </td><td> Atmosfera </td><td> Https://trinotate.github.io/ </td><td> Https: //pods.iplantcollaborative Org / wiki / display / atmman / Trinity + - + Trinotate + Atmosfera + Imagem </td><td> 1 semana </td><td> Https: //atmo.iplantcollaborative. Org / application / images / 1261 </td></tr></tbody></table> Tabela 1: Programas de Análise, Plataformas em que estão disponíveis, umNd Recursos adicionais disponíveis para os fluxos de trabalho em ordem pela primeira aparência. Todas as versões do pacote são atuais a partir de abril de 2016.

Este protocolo descreve um fluxo de trabalho comparativo de montagem e anotação de transcriptoma de novo para bioinformáticos iniciantes. O fluxo de trabalho está disponível gratuitamente inteiramente através do CyVerse e conectado pelo Data Store. A linha de comando e as interfaces gráficas do usuário são usadas, mas todo o código necessário está disponível para copiar e colar.

Aproveitando os Recursos CyVerse para<em&gt; De Novo</em&gt; Transcriptômica Comparativa de Organismos Não Servidos (não Modelos)

This workflow allows novice researchers to leverage advanced computational resources such as cloud computing to carry out pairwise comparative ...

O objetivo geral deste procedimento é avaliar, montar, anotar e comparar a expressão gênica diferencial por meio da transcriptômica De Novo, partindo de arquivos FASTQ brutos. Este método pode ajudar a responder questões em biologia comparativa e molecular, incluindo quais transcritos estão presentes em um organismo, o que esses transcritos estão fazendo dentro desses organismos e quais são as diferenças entre condições experimentais. A principal vantagem desta técnica é que ela fornece um ambiente interativo.Ele fornece recursos computacionais sob demanda e permite que pesquisadores comecem imediatamente a analisar seus dados de RNA-Seq. Este método é particularmente útil para pesquisadores que comparam experimentos dentro de um único organismo, envolvendo múltiplos tecidos, condições ou pontos temporais, para compreender como os sistemas biológicos mudam. Este método concentra-se em organismos não-modelo sem genomas, mas também pode ser aplicado a organismos com montagens genômicas disponíveis, mesmo aqueles com dezenas ou centenas de milhares de scaffolds em sua montagem.Para começar, obtenha acesso ao Atmosphere no Ambiente de Descoberta. Solicite uma conta gratuita no CyVerse acessando a página de registro. Utilize um e-mail institucional para se registrar na conta.Em seguida, acesse a guia de aplicativos e serviços e solicite acesso ao Atmosphere. O acesso ao Ambiente de Descoberta é concedido automaticamente. Faça login no Ambiente de Descoberta, abreviado como DE. Depois, selecione a guia Dados para exibir um menu contendo todas as pastas no repositório de dados.Crie uma pasta principal do projeto que abrigará todos os dados associados ao projeto. Localize a barra de ferramentas na parte superior da janela de dados e selecione Arquivo, Nova Pasta. Não use espaços nem caracteres especiais nos nomes das pastas ou em quaisquer nomes de arquivos de entrada e saída.Em vez disso, use sublinhados ou traços onde apropriado. Faça o upload dos arquivos brutos de sequência FASTQ e da pasta 1_Raw_Sequence para uma subpasta denominada Pasta A_Raw_Reads. Para arquivos com menos de dois gigabits, utilize o recurso de upload simples do armazenamento de dados para acessar a barra de ferramentas da janela de dados clicando no botão Dados no desktop principal do DE.Selecione Upload, Upload Simples da área de trabalho. Em seguida, selecione o botão Procurar para navegar até os arquivos de sequenciamento FASTQ brutos no computador local. Avalie as leituras de sequenciamento brutas carregadas utilizando o aplicativo FastQC no DE. Selecione o botão Aplicativos na área de trabalho principal do DE para abrir uma janela contendo todos os aplicativos de análise disponíveis no DE. Pesquise na janela a ferramenta FastQC na barra de pesquisa na parte superior da janela.Abra a versão de múltiplos arquivos se houver mais de um arquivo FASTQ. Selecione Arquivo e crie uma nova pasta, então escolha esta pasta como pasta de saída. Carregue os arquivos de leitura FASTQ na janela da ferramenta chamada Selecionar Dados de Entrada e selecione Iniciar Análise.Procure o aplicativo programável Trimmomatic no DE e abra-o. Faça o upload da pasta com os arquivos brutos de leituras FASTQ na seção de configurações. Selecione se os arquivos de sequenciamento são de extremidade simples ou pareados.Utilize o arquivo de controle padrão fornecido, selecionando o botão Procurar e colando o caminho do arquivo na caixa de visualização. Selecione o arquivo de controle do Trimmomatic e inicie a análise. Para o corte de qualidade das leituras de sequência, pesquise e abra o aplicativo Sickle no DE. Selecione as leituras FASTQ recortadas como leituras de entrada e renomeie os arquivos de saída.Inclua as configurações de qualidade nas opções. Abra a versão mais recente da instância do Atmosphere navegando até a página da wiki. Selecione o link para a versão mais recente da imagem do Trinity e do Trinotate.Selecione o botão Fazer login para iniciar e, em seguida, nomeie a instância Atmosphere. Selecione um tamanho de instância entre medium3 ou large3. Inicie a instância e aguarde a sua criação.Se uma imagem do Atmosphere não iniciar corretamente, você pode tentar solicitar uma instância menor ou solicitar ao Jetstream uma alocação maior. Todos os detalhes estão no wiki complementar. Mova os arquivos de saída do Trinity para a pasta 3_Assembly no DE e rotule a pasta como A_Trinity_de_novo_assembly.Executar o Trinity requer conhecimento de linha de comando e vários dias, ou possivelmente semanas, para concluir análises extensas. Existem recursos gratuitos disponíveis, vinculados à wiki, para ajudar a compreender a linha de comando. Dê a cada transcriptoma que foi montado uma subpasta dentro da pasta A_Trinity_de_novo_assembly.Utilize nomes exclusivos incluindo os nomes científicos dos organismos e os tratamentos associados a cada transcriptoma, em seguida crie outra subpasta chamada Pasta B_rnaQUAT_Output na pasta 3_Assembly. Abra o aplicativo intitulado De Novo rnaQUAST. Nomeie a análise e selecione a Pasta B_rnaQUAST_Output como pasta de saída.Procure por transcript decoder e execute o transdecoder no arquivo fasta de saída do Assembly De Novo Trinity no Ambiente de Descoberta. Abra o aplicativo deseq2 no DE. Nomeie a análise e selecione a pasta de saída como 4_Differential_Expresssion. Na seção de Entrada, selecione o arquivo da tabela de contagens do resultado do Assembly Trinity.Além disso, selecione a coluna onde os nomes dos contigs podem ser encontrados. Insira os cabeçalhos das colunas do arquivo da tabela de contagens para determinar quais colunas serão comparadas. Inclua vírgulas entre cada uma das condições.Não inclua o cabeçalho da primeira coluna que contém os nomes dos contigs. Para réplicas, repita o mesmo nome. Na segunda linha, forneça os nomes das duas condições a serem comparadas.Corresponda aos nomes dos cabeçalhos das colunas fornecidos na primeira linha. Aqui é mostrada uma comparação sistemática das leituras de sequenciamento após cada etapa de pré-processamento. Após o corte, a leitura deve apresentar conteúdo de GC menos distorcido e conteúdo de sequência mais equilibrado, além de uma proporção maior de leituras com pontuação de qualidade elevada.Leituras de alta qualidade são necessárias para montar transcriptomas De Novo. Os resultados da análise rápida de qualidade dependem dos organismos e amostras que estão sendo sequenciados. A uniformidade entre todas as amostras que serão comparadas posteriormente é o objetivo principal do pré-processamento das leituras.rnaQUAST utiliza o código boost para gerar estatísticas resumidas sobre montagens com base em genes centrais conhecidos em clados taxonômicos. A precisão dos montadores é revelada pelo número de incompatibilidades por transcritos e por quantos transcritos correspondem aos genes canônicos. Os últimos quatro subgráficos apresentados aqui fornecem estatísticas resumidas do comprimento de contigs e isoformas, bem como da cobertura das isoformas esperadas.NAx representa a porcentagem de contigs com comprimento maior que o comprimento do eixo y. A fração montada é o maior transcritos montado individualmente dividido pelo seu comprimento. Já a fração coberta é a porcentagem de isoformas de transcritos completos montados esperadas pelos genes procariontes ou eucariontes centrais do BUSCO.Após assistir a este vídeo, você deverá ter uma boa compreensão de como montar e analisar transcriptomas. Além disso, este protocolo permitirá que você detecte expressão gênica diferencial entre duas condições. Geralmente, as pessoas têm dificuldade com pacotes bioinformáticos porque existem muitos deles, há uma grande quantidade de configurações e variáveis associadas e, normalmente, é necessário ter conhecimento da linha de comando para executá-los.É importante rotular e organizar suas entradas de dados e saídas de análise para que outros pesquisadores possam entender o que foi feito. Você deve incluir a ordem em que os passos foram concluídos, as versões dos programas e as informações das amostras. Além disso, omita quaisquer espaços nos nomes de pastas ou arquivos.Ferramentas novas e novas versões das ferramentas estão sendo integradas constantemente, mas as versões antigas das ferramentas também estão sendo mantidas. Todas as alterações serão registradas na wiki complementar. Seguindo este procedimento, outros métodos bioinformáticos, como análise de redes, enriquecimento de GO e identificação de vias metabólicas, podem ser realizados para ajudar a responder perguntas como variação fenotípica, condições que alteram perfis de expressão e identificação de genes de interesse para genômica funcional.

leveraging-cyverse-resources-for-de-novo-comparative-transcriptomics

Research

JoVE Journal

Genetics

Genética

RNA-Seq as análises estão se tornando cada vez mais importantes para identificar os fundamentos moleculares de características adaptativas em organismos não modelo. Aqui, um protocolo para identificar genes diferencialmente expressos entre mosquitos Aedes albopictus diapausados e não diapausados é descrito, desde a criação de mosquitos até o sequenciamento de RNA e análises de bioinformática de dados de RNA-Seq.

Uma Experimental e Bioinformática Protocolo de Análises de RNA-seq de fotoperiódicas diapausa no Tigre Asiático Mosquito, Aedes albopictus

Galaxy e David surgiram como ferramentas populares que permitem que os investigadores sem formação bioinformática para analisar e interpretar dados de RNA-Seq. Descreve-se um protocolo para C. elegans investigadores para executar ARN-Seq experiências, o acesso e processar o conjunto de dados usando Galaxy e obter informação biológica significativa a partir das listas de genes utilizando DAVID.

Análise transcriptomic de<em&gt; C</em&gt;.<em&gt; elegans</em&gt; RNA Sequenciação de dados através da Suíte Tuxedo no Projeto Galaxy

Nós fornecemos um protocolo padronizado para a utilização da análise de enriquecimento conjunto do gene do transcriptomic de dados para identificar um modelo de mouse ideal para investigação de translação.Este protocolo pode ser usado com DNA microarray e dados de sequenciamento de RNA e mais pode ser estendido a outros dados omics se dados estão disponíveis.

Um protocolo para o uso de Gene definir análise de enriquecimento para identificar o modelo Animal adequado para investigação de translação

Aqui, apresentamos um protocolo para acessar e analisar muitas bases de dados de organismos humanos e modelo de forma eficiente. Este protocolo demonstra o uso de MARRVEL para analisar variantes causadoras de doenças candidatas identificadas a partir de esforços de sequenciamento de próxima geração.

Navegando no MARRVEL, uma ferramenta Web que integra informações da genômica humana e genética de organismos modelo

O objetivo deste protocolo é delinear o desenho e o desempenho de experimentos in vivo em Drosophila melanogaster para avaliar as conseqüências funcionais de variantes genéticas raras associadas a doenças humanas.

Estudo funcional in vivo de variantes humanas raras associadas à doença usando Drosophila

Este trabalho descreve um protocolo semi-alta de supressão que permite imagens simultâneas de lapso de tempo 3D de embriaênese embrionária em embriões de 80-100 C. elegans em uma única corrida noturna. Além disso, as ferramentas de processamento e visualização de imagens são incluídas para simplificar a análise de dados. A combinação destes métodos com tensões personalizadas do repórter permite a monitoração detalhada do embryogenesis.

Um método de imagem semi-alta-fonte e um kit de ferramentas de visualização de dados para analisar o desenvolvimento embrionário de C. elegans

O protocolo aqui apresentado descreve um pipeline completo para analisar dados de transcriptome de sequenciamento de RNA desde leituras brutas até análises funcionais, incluindo etapas de controle de qualidade e pré-processamento para abordagens analíticas estatísticas avançadas.

Análise de transcriptome de alto rendimento para investigar interações host-pathogen

Este protocolo descreve um fluxo de trabalho de culturas de células ex vivo ou in vitro para pré-processamento de dados transcriptômicos para triagem de drogas baseada em transcriptoma custo-efetivo.

Triagem Transcriptomic-Based Drug Screening Custo-Eficiente

O presente trabalho descreve um protocolo para executar o algoritmo Pathway2Targets, um script R que prevê e prioriza alvos terapêuticos com base no perfil das vias de sinalização intracelular geradas pela comparação de amostras de caso versus controle de um experimento de sequenciamento de RNA em massa.

Usando listas de genes diferencialmente expressos em humanos para realizar análise de enriquecimento de vias a jusante e priorização de alvos

Este artigo descreve a aplicação de metabolômica irrelevantes, transcriptomics e análise estatística multivariada para transcrições baga da uva e metabolitos, a fim de ter uma visão sobre o conceito terroir, ou seja, o impacto do meio ambiente sobre características de qualidade de baga.

O Conceito Terroir Interpretado através da baga da uva Metabolomics e transcriptómica

Aprenda a montar e anotar transcriptomas de novo usando recursos CyVerse para organismos não modelos. Este fluxo de trabalho interativo permite analisar dados de RNA-Seq e comparar expressão gênica diferencial.

Aproveitando os Recursos CyVerse para

Explorar mais vídeos