CAP 3 Contig Assembly Program George Darmiton da
CAP 3 (Contig Assembly Program) George Darmiton da Cunha Cavalcanti (gdcc@cin. ufpe. br) UFPE – CIn Junho de 2001
Roteiro n Introdução n Arquitetura n O Programa CAP 3 – Entrada e Saída n Pontos Fortes e Fracos do CAP 3 n Resultados do CAP 3 n Bibliografia – CAP 3 versus PHRAP
Introdução n n É um programa utilizado para montar cadeias de DNA Desenvolvido por Xiaoqiu Huang, – Department of Computer Science, Michigan Technological University n n Versão 3 foi desenvolvida em 1999 Atualmente encontra-se na versão 4, sendo esta comercial
Arquitetura Remoção de regiões pobres Calcular sobreposição 1º Fase Remoção de falsas sobreposições Construção de Contigs 2º Fase Construção do Consenso 3º Fase
Primeira Fase n Composta de três etapas – Remoção de regiões pobres – Cálculo da sobreposição – Remoção de falsas sobreposições n Antes dessas etapas é necessário identificar sobreposições entre fragmentos.
Primeira Fase Identificação de sobreposição n Criar a seqüência combinada – Os fragmentos f 1, f 2, . . . , fn são concatenados – Caracter especial separa dois fragmentos – f 1 # f 2 #. . . # fn n Para cada fragmentos fx e o seu complemento reverso rx, encontrar o par (fx, fy) e (rx, fy) – tal que x<y e que os fragmentos tenham uma sobreposição relevante n Para cada par com sobreposição uma faixa diagonal da matriz da programação dinâmica e calculado. (Smith e Waterman 1981) – Essa matriz será usada mais adiante por motivo de eficiência
Primeira Fase Remoção de Regiões Pobres Remoção de regiões pobres Calcular sobreposição 1º Fase Remoção de falsas sobreposições Construção de Contigs 2º Fase Construção do Consenso 3º Fase
Primeira Fase Remoção de Regiões Pobres (cont. ) Posição de remoção 5’ Fragmento h Fragmento f Fragmento g Posição de remoção 3’
Primeira Fase Remoção de Regiões Pobres (cont. ) n O algoritmo de alinhamento local de Smith e Waterman foi generalizado para usar valores de qualidade de base Sendo p a probabilidade do erro estimado para a base m * min(q 1, q 2) n * min(q 1, q 2) -g * min(q 1, q 2)
Primeira Fase Remoção de Regiões Pobres (cont. ) n Os valores que indicam a qualidade da base são usados para permitir que: – Matches em bases que possuem altos valores de qualidade recebem alta pontuação positiva; – Mismatches em bases que possuem altos valores de qualidade recebem alta pontuação negativa; – Matches e mismatches em bases que possuem baixos valores de qualidade recebem pontuações baixas positivas e negativas, respectivamente
Primeira Fase Remoção de Regiões Pobres (cont. ) Se os valores de qualidade de base foram informados qualpos 5 crange (-y) qualpos 3 crange Fragmento f Maioria dos valores de qualidade são maiores que qualcut (-c) n Caso os valores de qualidade da base não sejam informados – qualpos 5 = 1 e qualpos 3 = tamanho de f n A cobertura mínima é determinada por gdepth (-z)
Primeira Fase Cálculo das Sobreposições Remoção de regiões pobres Calcular sobreposição 1º Fase Remoção de falsas sobreposições Construção de Contigs 2º Fase Construção do Consenso 3º Fase
Primeira Fase Cálculo das Sobreposições (cont. ) n O alinhamento global é utilizado para calcular a sobreposição entre fragmentos n Alinhamento global versus Alinhamento local – Utilizando o método global é possível identificar falsas sobreposições. • mostrar que algumas regiões dos fragmentos não similares, indicando que esta sobreposição é falsa. – O alinhamento local está restrito a regiões similares.
Primeira Fase Remoção de Sobreposições Falsas Remoção de regiões pobres Calcular sobreposição 1º Fase Remoção de falsas sobreposições Construção de Contigs 2º Fase Construção do Consenso 3º Fase
Primeira Fase – Remoção de Sobreposições Falsas (cont. ) n Cada sobreposição é avaliada por 5 (cinco) medidas – 1ª Medida • Informa o comprimento mínimo para a sobreposição (-o) – 2ª Medida • Determina que o percentual de identidade não deve ser menor que o valor estabelecido pela opção –p – 3ª Medida • • Determina o valor de similaridade da sobreposição (-s) m * min(q 1, q 2) n * min(q 1, q 2) -g * min(q 1, q 2)
Primeira Fase – Remoção de Sobreposições Falsas (cont. ) – 4ª Medida • Se a sobreposição contiver um número grande de diferenças entre bases de altos valores de qualidade, essa sobreposição é provavelmente falsa. (-b e -d) – 5ª Medida • Se o número de diferenças em uma sobreposição for maior que o esperado, então é provável que essa sobreposição seja falsa. (-e) n Caso uma dessas medidas falhe, a sobreposição é considerada falsa.
Segunda Fase Construção dos Contigs Remoção de regiões pobres Calcular sobreposição 1º Fase Remoção de falsas sobreposições Construção de Contigs 2º Fase Construção do Consenso 3º Fase
Segunda Fase Construção dos Contigs (cont. ) n 1ª Etapa – Um layout inicial é gerado • Método guloso n 2ª Etapa – A qualidade do layout corrente é avaliada – O número de restrições satisfeitas e não satisfeitas é calculado para cada sobreposição – Restrições não satisfeitas são particionadas em grupos • cada grupo possui restrições associadas com uma sobreposição não usada ou com um par de contigs
Segunda Fase Construção dos Contigs (cont. ) n 3ª Etapa – O grupo com o maior número de restrições não satisfeitas é selecionado – Caso 1 • grupo associado a uma sobreposição não usada – Caso 2 • grupo associado a um par de contigs – Se nenhuma correção for feita • o processo é repetido com os grupos restantes – Caso contrário • a 2ª etapa é repetida para o novo layout
Terceira Fase Construção do Consenso Remoção de regiões pobres Calcular sobreposição 1º Fase Remoção de falsas sobreposições Construção de Contigs 2º Fase Construção do Consenso 3º Fase
Terceira Fase Construção do Consenso (cont. ) n A soma ponderada dos valores de qualidade é calculada para cada base – Os valores de qualidade são divididos em dois grupos, um para cada sentido (5’ 3’ - 3’ 5’) – cada grupo é ordenado em ordem decrescente – pesos: w 1=1, wi = 0. 5, para i>1 – Ex: 20+, 40 -, 30+ e 10 • grupo 1: 30+ e 20+ grupo 2: 40 - e 10 • soma ponderada = 30(1) + 40(1) + 20(0. 5) + 10(0. 5) = 85
Terceira Fase Construção do Consenso (cont. ) n Cálculo da média do valor de qualidade
Terceira Fase Construção do Consenso (cont. )
Parâmetros de Entrada do CAP 3 n CAP 3 recebe um arquivo com as seqüências de fragmentos no formato FASTA – Uso: cap 3 arquivo_de_fragmentos [opções] n Arquivos opcionais – arquivo contendo os quality values no formato FASTA, usando extensão. qual – arquivo contendo restrições forward-reverse, usando extensão. con. • Pode ser gerado usando o programa FORMCON • Formato: Read. A Read. B Min. Distance Max. Distance
Opções Valor das penalidades Parâmetro descrição -g N gap N>0 (6) -m N match N>0 (2) -n N mismatch N<0 (-5) Faixa Diagonal Parâmetro descrição -a N N>10 (20)
Resposta do Programa n Consenso no formato ace – arquivo com extensão. ace n Consenso – arquivo com extensão. contigs n Quality values do consenso – arquivo com extensão. contigs. qual n Fragmentos não são usados na montagem – arquivo com extensão. singlets
Resposta do Programa (cont. ) n Informações adicionais sobre a montagem – arquivo com extensão. info n Satisfação das restrições – arquivo com extensão. results
Pontos Fortes do CAP 3 n Uso de forward-reverse constraints para corrigir erros de montagem – Objetivo: localizar e corrigir erros no layout da seqüência e ligar contigs separados por gaps – Dois fragmentos devem estar em direções opostas na molécula de DNA e a uma determinada distância. – O algoritmo usado no CAP 3 é tolerante a restrições erradas
Pontos Fortes do CAP 3 (cont. ) n Geração do resultado da montagem no formato ace para Consed – ferramenta gráfica para editar seqüências n CAP 3 pode ser usado no GAP 4 do pacote Staden. – GAP 4(Genome Assembly Program) é uma ferramenta gráfica do pacote de ferramentas Staden
Pontos Fortes do CAP 3 (cont. ) n Uso de base quality values – Usados no alinhamento de fragmentos e na construção do consenso. – Melhora a qualidade na geração do consenso n Remoção de regiões pobres, 5’ e 3’ – Objetivo • Utilizar apenas regiões ‘boas’ do fragmento na montagem.
Pontos Fracos do CAP 3 n A remoção de regiões ditas pobres, pode excluir áreas importantes no processo de alinhamento n Tempo de processamento n A ferramenta não possui interface gráfica (entretanto pode ser usada no pacote GAP 4)
Resultados do CAP 3 Conjunto de dados BAC (Bacterial Artificial Chromossome) Dados Gen. Bank Nº Nº de frag Tam médio frag Tam seq 203 AC 004669 1812 598 89. 779 216 AC 004638 2353 614 124. 645 322 F 16 AF 111103 4297 1011 159. 179 526 N 18 AF 123462 3221 965 180. 182 Dados Tempo (min) Nº de contigs Tam seq CAP Nº de <> 203 37 1 90. 292 0 216 154 1 132. 057 17 322 F 16 127 1 157. 982 11 526 N 18 73 2 180. 128 10
CAP 3 versus PHRAP n Conjunto de dados BAC – PHRAP, normalmente, produz cadeias mais longas de contigs – CAP 3 produz menos erros no consenso n Caso os valores de qualidade da base não estejam disponíveis – CAP 3 é uma boa escolha já que trata redundância
Bibliografia n Xiaoqiu Huang, Anup Madan. CAP 3: A DNA Sequence Assembly Program. Genome Research 9: 869 -877, 1999. n Xiaoqiu Huang. Na Improved Sequence Assembly Program. Genomics 33, 21 -31, 1996. n Site oficial na Internet. http: //genome. cs. mtu. edu/cap 3. html n Staden Package WWW site. – http: //www. mrc-lmb. cam. ac. uk/pubseq/staden_home. html – http: //www. mrc-lmb. cam. ac. uk/pubseq/contig. html
- Slides: 34