Kaiju
Classificação TaxonômicaKaiju é um classificador metagenômico em nível de proteína que traduz leituras de DNA em sequências de aminoácidos e as compara com um banco de dados de proteínas de referência usando a transformada de Burrows-Wheeler (BWT). [1]
Como Obter o Arquivo de Saída do Modelo
Abaixo está um breve fluxo de trabalho que a equipe executou para obter os exemplos de arquivo de saída que apresentamos na página de ferramentas.
Entrada
Arquivos FASTQ/FASTA (leituras single-end ou paired-end)
Saída
Classificação separada por tabulação com IDs de táxons, além de tabelas de resumo em cada rank taxonômico
conda install -c bioconda kaiju
Imagem Docker: quay.io/biocontainers/kaiju:1.10.1--h43eeafb_0
Amostra 1 leituras de microbioma intestinal (subconjunto de 200K) classificadas com Kaiju contra um Banco de dados curado de microbioma intestinal (644 genomas representativos/referência do RefSeq do NCBI).
- 1
Baixar leituras do NCBI SRA
prefetch SRR14092310 && fasterq-dump SRR14092310 -O /data --split-files && gzip /data/SRR14092310_*.fastq
- 2
Subamostrar para 200K pares de leituras
zcat SRR14092310_1.fastq.gz | head -800000 > subset_R1.fastq && zcat SRR14092310_2.fastq.gz | head -800000 > subset_R2.fastq
O subconjunto de 200K mantém a saída abaixo do limite de upload de 50 MB do IntMeta.
- 3
Executar classificação Kaiju
kaiju -t nodes.dmp -f kaiju_db_gut.fmi -i subset_R1.fastq -j subset_R2.fastq -o kaiju_raw.tsv -z 12
- 4
Adicionar nomes de táxons
kaiju-addTaxonNames -t nodes.dmp -n names.dmp -i kaiju_raw.tsv -o kaiju_output.tsv -r superkingdom,phylum,class,order,family,genus,species -u
Materiais Utilizados
Arquivos de Saída de Exemplo
Baixe os arquivos de saída usados nas demonstrações da página de ferramentas. Você pode enviá-los diretamente para IntMeta para explorar as visualizações.
Única & Comparação
Análise de Grupo
Referência de Gráficos
Descrições detalhadas para todas as 30 visualizações geradas por Kaiju em IntMeta.
distribution
Gráfico de barras dos principais organismos classificados por contagem de leituras da classificação em nível de proteína de Kaiju. Como Kaiju traduz leituras para aminoácidos antes de comparar, pode detectar organismos divergentes que métodos de k-mers de nucleotídeo perdem.
composition
Gráfico de pizza/rosca mostrando os táxons principais como percentual da sua contagem combinada de leituras no rank selecionado, derivado das atribuições em nível de proteína BWT de Kaiju.
richness
Conta cada organismo distinto com pelo menos uma leitura atribuída em cada rank principal (Domain até Species). A classificação em nível de proteína tende a detectar táxons adicionais divergentes ou novos em comparação com ferramentas de k-mers de DNA.
diversity
Índices de diversidade alfa calculados a partir de proporções de contagem de leituras por táxon: Shannon H = −Σ(pᵢ · ln pᵢ), Simpson D = 1 − Σ(nᵢ(nᵢ−1))/(N(N−1)), e uniformidade de Pielou J = H / ln(riqueza). Mede riqueza e uniformidade da comunidade a partir de dados de abundância de Kaiju.
multilevel-composition
Gráfico de barras empilhadas com os principais táxons (por contagem de leituras) em cada rank principal; táxons restantes agrupados como 'Outro'. Revela mudanças composicionais do nível de Filo ao nível de Espécie.
dependency-wheel
Diagrama de cordas vinculando táxons pai a táxons filho através dos ranks da saída de linhagem de Kaiju. A espessura da conexão é proporcional à contagem compartilhada de leituras; arestas abaixo do limite mínimo de cobertura são filtradas.
sankey-flow
Diagrama de fluxo Sankey mostrando como as leituras se distribuem de ranks superiores a ranks inferiores através da hierarquia de classificação de Kaiju. A largura da banda é igual à contagem de leituras; arestas de baixa cobertura são removidas.
comp-grouped-abundance
Gráfico de barras agrupadas colocando os táxons principais de cada amostra lado a lado no rank selecionado. Cada grupo contém uma barra por amostra, colorida pela identidade da amostra, permitindo comparação visual direta das contagens absolutas de leituras para o mesmo organismo entre amostras classificadas por Kaiju.
comp-relative-abundance
Gráfico de barras empilhadas 100% onde cada barra representa uma amostra e os segmentos mostram a contribuição proporcional de cada táxon. Útil para comparar composição da comunidade quando as amostras têm profundidades de sequenciamento muito diferentes, pois todas as barras são normalizadas para a mesma altura.
comp-abundance-heatmap
Mapa de calor de matriz de cores com táxons em um eixo e amostras no outro. A intensidade da cor da célula é proporcional à abundância de leituras no rank selecionado. O agrupamento hierárquico em ambos os eixos agrupa amostras similares e táxons co-ocorrentes.
comp-diversity-indices
Gráfico multi-painel exibindo entropiaShannon, diversidade Simpson, Riqueza Observada e Uniformidade de Pielou para cada amostra em sua própria escala de eixo y. Permite comparação rápida entre amostras de diversidade alfa sem distorção de escala.
comp-shared-taxa
Diagrama de Venn mostrando a contagem de táxons compartilhados entre amostras versus táxons exclusivos para cada amostra individual. Calculado no rank taxonômico selecionado usando presença/ausência de leituras ≥ 1.
group-alpha-diversity
Boxplots de entropiaShannon, diversidade Simpson, Riqueza Observada e Uniformidade de Pielou por grupo. Cada caixa resume a variação intra-grupo. Valores p de Kruskal-Wallis testam diferenças significativas entre grupos.
group-pcoa
Análise de Coordenadas Principais (PCoA) na matriz de dissimilaridade Bray-Curtis. Pontos são coloridos por atribuição de grupo com elipses de confiança de 95%. Rótulos de eixo mostram o percentual de variância explicada por cada coordenada.
group-nmds
Escala Multidimensional Não-Métrica (NMDS) com valor de stress exibido. Pontos coloridos por grupo com elipses de confiança. Stress menor (<0,2) indica uma boa representação das distâncias originais em 2D.
group-distance-boxplots
Boxplots de distâncias Bray-Curtis intra-grupo versus inter-grupo. PERMANOVA R² e p-valor quantificam a fração da variância explicada pelo agrupamento. Estatística R de ANOSIM mede o grau de separação dos grupos.
group-relative-abundance
Abundância relativa média do grupo como barras empilhadas 100% no rank taxonômico selecionado. Cada barra mostra a composição proporcional média em todas as amostras do grupo, permitindo comparação direta em nível de grupo.
group-differential-abundance
Gráfico de barras de táxons com diferenças significativas de abundância entre grupos (teste de Kruskal-Wallis, correção BH FDR). Cor indica qual grupo tem maior abundância. Apenas táxons que passam o limite de significância são mostrados.
group-lefse
LEfSe (Análise de Discriminante Linear Tamanho de Efeito) gráfico de biomarcador. Barras horizontais mostram escores LDA para táxons que discriminam significativamente entre grupos. Escores LDA mais altos indicam associação mais forte com o respectivo grupo.
group-shared-taxa
Diagrama de Venn mostrando táxons compartilhados e únicos através de grupos experimentais. O conjunto de táxons de cada grupo é a união de todos os táxons detectados em qualquer amostra pertencente a esse grupo. Clique nas regiões para visualizar a lista específica de táxons com contagens de leituras.
group-classification
Gráfico de colunas empilhadas mostrando leituras médias classificadas versus não classificadas por grupo experimental. O tooltip mostra porcentagem e contagem de amostras por grupo.
group-distribution
Gráfico de colunas agrupadas mostrando contagens médias de leituras para os N principais táxons no nível taxonômico selecionado, agrupados por grupo experimental. Permite comparação de padrões de abundância taxonômica entre grupos.
group-heatmap
Mapa de calor de contagens médias de leituras transformadas em log₁₀ por grupo experimental (colunas) e táxons (linhas). A intensidade da cor reflete abundância, permitindo comparação visual rápida de perfis taxonômicos entre grupos.
rarefaction
Curva de rarefação plotando táxons observados versus profundidade de leitura subamostrada no rank selecionado. Calculada analiticamente usando expectativa hipergeométrica. Uma curva que se estabiliza indica profundidade de sequenciamento suficiente; uma curva ainda em ascensão sugere subamostragem.
rank-abundance
Curva de classificação-abundância (Gráfico de Whittaker) com táxons classificados por abundância relativa decrescente em eixo y logarítmico. O comprimento da curva ao longo do eixo x reflete riqueza, enquanto a inclinação indica uniformidade — uma queda abrupta significa que alguns poucos táxons dominam.
comp-rarefaction
Curvas de rarefação sobrepostas para todas as amostras carregadas. Permite comparação visual do esforço de sequenciamento e saturação taxonômica entre amostras.
comp-rank-abundance
Curvas de classificação-abundância sobrepostas para todas as amostras carregadas. Curvas mais íngremes indicam comunidades mais desiguais dominadas por menos táxons.
group-taxonomic-sunburst
Sunburst interativo exibindo a hierarquia taxonômica completa como anéis concêntricos. Nós de biomarcador LEfSe são coloridos pelo seu grupo enriquecido; nós não significativos são cinzentos. Clique em qualquer segmento de anel para fazer zoom.
group-volcano
Gráfico vulcão: log₂ alteração de dobro (x) versus −log₁₀ p-valor ajustado (y). Pontos acima da linha horizontal passam significância (BH FDR). Pontos coloridos estão significativamente enriquecidos no respectivo grupo.
comp-classification
Gráfico de barras empilhadas mostrando leituras totais classificadas versus não classificadas por amostra. Permite comparação rápida das taxas de classificação entre amostras — frações não classificadas grandes podem indicar organismos novos ou limitações do banco de dados.