Preparação da amostra e análise de dados da expressão do Gene RNASeq-baseado de Zebrafish

Summary

Este protocolo apresenta uma abordagem para a análise do transcriptoma toda de zebrafish embriões, larvas, ou células de classificados. Nós incluímos o isolamento do RNA, análise de caminho de dados de RNASeq e validação de qRT-PCR-baseados de mudanças de expressão do gene.

Abstract

A análise das alterações de expressão do gene global é uma ferramenta valiosa para identificar novos caminhos subjacentes fenótipos observados. O peixe-zebra é um excelente modelo para avaliação rápida do transcriptoma toda de populações de células individuais ou todo animal devido à facilidade de isolamento do RNA de um grande número de animais. Aqui apresenta-se um protocolo para análise de expressão de gene global em embriões de peixe-zebra usando a sequenciação do ARN (RNASeq). Descrevemos a preparação do RNA de embriões de todo ou de populações de células obtidas usando célula de triagem em animais transgénicos. Também descrevemos uma abordagem para análise de dados RNASeq para identificar caminhos enriquecidos e termos de Gene Ontology (GO) em conjuntos de dados de expressão de gene global. Finalmente, nós fornecemos um protocolo para a validação das alterações de expressão do gene usando transcriptase reversa quantitativo PCR (qRT-PCR). Esses protocolos podem ser usados para análise comparativa de controle e experimental de moda de zebrafish para identificar alterações de expressão do gene de romance e fornecem a introspecção molecular de fenótipos de interesse.

Introduction

Análise comparativa da expressão gênica global é uma ferramenta valiosa para identificar novos genes, contribuindo para fenótipos observados. Tais análises normalmente contam com avaliação quantitativa da abundância de transcrição em comparação entre experimental e amostras de controlo. Abordagens específicas, tais como qRT-PCR são relativamente rápida e exata para investigação de alterações de expressão de gene único. A sequenciação do ARN (RNASeq) oferece uma abordagem ampla, livre de hipótese para identificar alterações significativas na expressão gênica entre as amostras, tornando-se agora o padrão para tais investigações através de sistemas experimentais.

Zebrafish têm emergido como um modelo de destaque em várias áreas da doença. Originalmente desenvolvido para sua utilidade em estudos de biologia do desenvolvimento, devido a sua alta fecundidade e relativamente baixo custo de manutenção, utilização experimental de zebrafish evoluiu para incluir uma ampla gama de fenótipos de embrionárias para estágios adultos, bem como um ampla gama de molecular ensaios^1,2,3. Com efeito, estas vantagens fazem estudos moleculares mecanicistas rápida e econômica, devido à facilidade de adquirir grandes quantidades de material combinado com a facilidade de manipulação genética e ambiental em todas as fases da vida. Além disso, a natureza transparente de zebrafish embriões e larvas tornam ideal para a geração de linhas de transgénicos repórter células e tecidos específicos permitindo visualização na vivo de populações de células discretas⁴. Exploração de tais linhas permite a análise da expressão de gene global em tipos específicos de célula isolada baseados na expressão do gene repórter.

Aqui nós apresentamos um protocolo completo para análise de expressão do gene global usando o RNASeq após a cultura de embriões de peixe-zebra. Manipulações genéticas experimentais, incluindo Morpholinos (MO)-knockdown com base transitória do gene ou genoma mediada por CRISPR edição, foram apresentadas em outro lugar^5,6,7. Nós, portanto, concentrar-se em um protocolo detalhado para isolamento de RNA de embriões todo ou classificados células transgénicas repórter-expressando seguidas por simples análise computacional de RNASeq resultados usando ferramentas de caminho e gene ontology (GO) termos. Finalmente, nós incluímos uma estratégia para a validação das alterações de expressão do gene por transcriptase reversa quantitativo PCR (qRT-PCR). Estes protocolos são aplicáveis aos embriões de zebrafish, sujeitados a uma ampla gama de condições experimentais, incluindo comparação de mutantes genéticos ou condições ambientais.

Protocol

todos os protocolos animais descritos abaixo estão em conformidade com e aprovado pela Universidade de Maryland institucional Animal Care e Comissão de utilização (IACUC).

1. preparação do embrião

1. gerar embriões através de acasalamento natural
   1. embriões de cultura para 3 meses de idade, reprodutiva maturidade ^{5 , 8} .
   2. Separe peixes adultos masculinos e femininos da estirpe desejada em tanques de acasalamento divididos na noite antes da coleta de embrião e adicionar 2 machos e 3 fêmeas para cada tanque.
      Nota: Uso da cepa transgênica insulin2a:mCherry fluorescente repórter permitiu a análise das células β-pancreáticas.
   3. Transferência de peixes para o acasalamento do tanque com água fresca do sistema e remover o divisor imediatamente após as luzes acendem-se na manhã seguinte.
   4. Permitir que os peixes acasalar, naturalmente, até que os embriões são observados no reservatório inferior. Colete embriões em intervalos de 30 min até a quantidade desejada é coletada. Coletados de loja cada ponto do tempo em pratos de Petri separado em embriões em 28,5 ° C.
   5. Executar microinjection de material genético ou colocação em meios de cultura experimental ⁶, se desejado e cultura de embriões frescos Hank ' s embrião médio ⁸ em pratos de Petri de 10 cm a 28,5 ° C.
      Nota: para a análise de expressão do gene em Morpholinos (MO) injetado ou animais mutantes, observe que cada manipulação pode ter impactos acidentais na expressão gênica. Mutantes podem apresentar compensação genética ao nível da transcrição que não foi observada pelo gene baseados em MO direcionamento ⁹.
2. Estágio
   1. cultura de embriões em grupos de 50 – 75 embriões por prato de Petri de 10 cm para promover o sincronismo do desenvolvimento consistente de todos os embriões.
   2. Monitorar a morfologia do desenvolvimento usando microscópio dissecação em blastômero, epiboly e somite fases para garantir a progressão do desenvolvimento ¹⁰.
      Nota: Remova qualquer moribundos ou malformados embriões para evitar o atraso de desenvolvimento no prato.
   3. Separar os embriões baseados na idade do desenvolvimento. Medir a idade do embrião usando somite número após a segmentação (fertilização post pós-gastrulação, 10,33 h (hpf)) até aproximadamente 24 hpf. Encenar os embriões e larvas usando o comprimento total do corpo após 24 hpf.
      Nota: As somitas são tecidos em forma de chevron mesodérmica presentes na parte dorsal do embrião.
   4. Colocar os embriões classificados em uma incubadora de 28,5 ° C e permitir o desenvolvimento de progredir para a idade desejada.

2. Célula única dissociação: Todo embrião e classificada de populações de células

1. dissociação todo embrião
   1. eutanásia os embriões de zebrafish na fase desejada, colocando o prato de Petri no gelo por 5-15 min, até que nenhum movimento é observado .
   2. Transferir uma piscina de 20 embriões para um tubo de microcentrifugadora rotulado de 1,5 mL. Remover qualquer excesso de embriões do tubo do microcentrifuge.
   3. Lise de adicionar 200 µ l de reagente (ver Tabela de materiais) para o tubo contendo embriões. Homogeneizar os embriões no reagente de Lise mecanicamente, utilizando um pilão. Adicione o reagente de Lise 800 µ l para trazer o volume total de 1 mL. Loja a-80 ° C até a extração do RNA.
2. De isolamento da célula de fluxo-classificados ¹¹
   1. transferir os embriões do prato de Petri em um tubo de microcentrifugadora de 1,5 mL e remover tanto médio de embrião possível.
      Nota: Dependendo da idade do embrião, dissociação mecânica pode igualmente ser necessária neste passo. Com embriões > 48 hpf, recomendamos o uso de um pilão para romper a integridade do embrião antes de prosseguir.
   2. Incubar os embriões com 1 mL de tampão de dissociação 1 (ver Tabela de materiais) no tubo de 1,5 mL microcentrifuge. Incube durante 15-30 min à temperatura ambiente. Triture a solução suavemente pipetando acima e para baixo usando uma dica P1000 para misturar a cada 3-4 min durante a etapa de incubação. FAZER não VORTEX.
   3. Recolher os embriões por centrifugação por 3 min a 300 x g e remover o sobrenadante. Re-suspenda os embriões em 1 mL de tampão de dissociação 2 (ver Tabela de materiais). Incubar durante 15-30 min à temperatura ambiente com trituração de pipeta regular.
   4. Avaliar a digestão diluindo-se 1-2 µ l da suspensão em fluorescência-ativado da pilha classificação Reserva (FACS) sob microscópio.
      Nota: Digestão completa ocorreu quando a alíquota examinada mostra células únicas com clusters de célula mínima e pedaços de tecido embrionário.
   5. Coletar as células por centrifugação a 300 x g, durante 5 min. remover sobrenadante. Ressuspender as células em 1 mL FACS buffer e filtro de gravidade através de um filtro de célula de 40 µm para remover o tecido não digerido da amostra.
   6. Contar as células usando um hemocytometer e diluir com tampão de FACS para aproximadamente 1 x 10 ⁶ células/mL em um tubo de FACS.
      Nota: Para os tipos de células com um número relativamente pequeno por embrião (menos de 5% do total), tais como células β-pancreáticas, usar um mínimo de 1.000 embriões estágio correspondente.
   7. Fornecer os FACS classificação facilidade com o FACS tubos contendo 1 mL de amostras de suspensão de células, bem como o tubo de FACS contendo apenas os FACS Buffer ou reagente de Lise. Portão do fluxo FACS-classificador para coletar apenas única células que expressam o repórter fluorescente.
   8. Executar o fluxo FACS-classificação até o número de células desejado ter sido colhido.
      Nota: Para executar RNASeq, há um mínimo necessário de RNA total. Encontramos que 3.000-5.000 células são suficientes para isolar a alta qualidade, alta concentração de RNA.
   9. Armazenar as células coletadas no gelo e proceder de imediato à preparação do RNA.

3. Preparação do RNA

1. Extrair o RNA
   1. incubar a amostra coletada (da etapa 2) com o reagente de Lise por 5 min à temperatura ambiente.
   2. Adicionar 0,2 mL de clorofórmio para cada 1,0 mL de reagente de Lise usado e inverter os tubos à mão por 15 s. Incubar por 2-3 min à temperatura ambiente. Centrifugar as amostras a 12.000 x g por 15 min a 4 ° C.
   3. Transferir a fase aquosa, separada para um novo tubo e adicionar 0,5 mL de isopropanol para cada 1,0 mL de reagente de Lise usado. Incube durante 10 minutos à temperatura ambiente. Centrifugar as amostras a 12.000 x g por 10 min a 4 ° C.
   4. Lavagem com 75% de etanol e centrifugar as amostras a 7.500 x g por 5 min a 4 ° C. Retire o sobrenadante completamente e deixe secar à temperatura ambiente. Ressuspender o RNA em 15-30 µ l pyrocarbonate de dietilo (DEPC)-água tratada.
2. Purify alta qualidade RNA
   1. combinar a suspensão de RNA com acetato de sódio 3M volume 1/10 (pH 5,5) e 1 volume de isopropanol. Incubar durante 20 min à temperatura ambiente e centrifugar as amostras a 12.500 x g por 10 minutos a 4 ° C.
   2. Lavar a pelota com etanol a 70% gelada na água tratada DEPC e centrifugar as amostras a 10.000 x g por 5 min a 4 ° C. repetir a etapa de lavagem uma vez.
   3. Remover o sobrenadante e deixar para secar à temperatura ambiente. Re-suspender em água tratada DEPC.
   4. Do ensaio o RNA extraído para a concentração (ng / µ l) e pureza (proporção de 260/230) usando um espectrômetro de absorção. Certifique-se de que o valor de 260/230 é ~ 2.0 antes de prosseguir com a RNASeq. Armazenar a-80 ° C até nose (até 6 meses).
   5. Enviar as amostras de RNA para um fornecedor ou núcleo para RNASeq e expressão gênica alterar análise com base na quantificação das leituras de sequenciamento. Os resultados retornados do fornecedor normalmente são fornecidos como uma lista de genes, exibindo uma dobra-mudança significativa na transcrição lê entre amostras.
      Nota: Uma coluna pode ser usada se o método acima produz RNA de má qualidade. A quantidade, concentração e número de integridade de RNA (RIN) para amostras do RNA devem ser confirmadas com o fornecedor ou o núcleo. O fornecedor ou o núcleo avaliará também tipicamente RIN.

4. Percurso e análise do termo GO

Nota: consulte a Figura 1 para saída representativa da análise de expressão de gene após RNASeq fornecido pelo fornecedor ou núcleo.

1. Sorting diferencialmente expressos genes
   1. comparando uma única condição experimental versus controle
      1. Abrir dados (resultados) em um software de gestão de folha de cálculo (ver tabela de materiais ).
      2. Marque a soltar Seta para baixo ao lado do ' tipo ' e seleccione " Custom tipo " dentro da planilha que contém os genes diferencialmente expressos no experimental versus condição controle.
      3. Marque a caixa em ' coluna ' na janela que aparece em cima e escolha classificar pelo ' LFC ' coluna. Certifique-se de ' valores ' é selecionado no ' tipo na ' coluna. Selecione esta opção para classificar do ' maior para o menor ' na ' ordem ' coluna e clique " Okey ".
         Nota: Isto irá classificar os genes diferencialmente expressos, para que aqueles com expressão aumentada (LFC positivo) irão aparecer no topo da lista, enquanto aqueles com diminuição da expressão (LFC negativo) aparecerá na parte inferior da lista.
   2. Comparar as várias condições experimentais para um único controle da seguinte maneira.
      1. Selecione a primeira célula em uma coluna vazia no 1 Experimental contra a planilha de controle para determinar diferencialmente expressos genes encontrados em duas condições experimentais. Digite a seguinte equação para a célula:
         
         Nota: esta equação é uma combinação de IF, ÉERROS e funções do jogo. (i) a função CORRESP, MATCH (A1, Experimental2_vs_Control! A:A, 0), irá procurar o valor em A1 (a ID de ENSEMBL) na coluna A (A:A) da planilha chamada Experimental2_vs_Control (Experimental2_vs_Control!). O " 0 " indica para procurar uma correspondência exata, e se for encontrada uma correspondência exata, a função retornará um valor de " 1 ", enquanto se nenhuma correspondência for encontrada, a função retornará um erro " N/A ". (ii) o resultado da função MATCH, em seguida, é entrado a função ÉERROS, ÉERROS (MATCH (A1, Experimental2_vs_Control! A:A, 0)), onde, se a entrada for " 1 " indicando uma correspondência foi encontrada, ele irá retornar " falso ", e se a entrada for " N/A " indicando uma correspondência não foi encontrada, ele irá retornar " verdade ". (iii) o " TRUE " ou " falso " resultado então é entrada para a função se, se (ISERROR (MATCH (A1, Experimental2_vs_Control! A:A, 0)), " ", " duplicar "), onde, se a entrada for " TRUE " indicando uma correspondência não foi encontrada, ele irá retornar o valor entre o primeiro conjunto de aspas (" ") que está vazio e, portanto, a célula será vazia. Se a entrada for " falso " indicando uma correspondência foi encontrada, ele irá retornar o valor entre o segundo conjunto de aspas (" duplicar "), e a célula vai dizer " duplicar ".
      2. Imprensa " Enter " ou " voltar " para executar a equação. Selecione a célula que contém a equação e clique na caixa no canto inferior direito da célula. Manter o mouse clicado e arraste a área selecionada toda a maneira para baixo a coluna até a última ' recurso ID ', para copiar a equação para cada célula na coluna.
      3. Select " Custom tipo " novamente, adicionar um segundo nível de classificação clicando o " + " ícone no canto inferior esquerdo. No primeiro nível, " classificar por ", em Selecionar coluna " duplicar ", em selecione Classificar na ' valores ' e em selecionar ordem ' Z a '. No segundo nível, " por ", em Selecionar coluna ' LFC ', em selecione Classificar na ' valores ' e em selecionar ordem ' maior para o menor ' e selecione " Okey ".
         Nota: Isto irá classificar a lista de genes em 4 grupos, com base na direção da mudança de expressão. É importante verificar os genes encontrados em ambas as condições experimentais para garantir que a mudança na expressão é na mesma direção em ambos os casos, ou em direções opostas.
      4. Remover qualquer parênteses da coluna de símbolos do gene, antes de prosseguir ou os resultados não serão encontrados nos bancos de dados. Selecione a coluna inteira que contém símbolos de gene. Selecione " editar " então " substituir " no menu drop-down ' arquivo ' menu. Tipo " (*) " no " encontrar o que: " barra do substituir janela e deixar o " substituir com: " bar vazio. Selecione ' substituir todos os ' para remover todas as instâncias de parênteses.
         Nota: O * representa qualquer número de caracteres, colocando este caractere entre dois parênteses, o programa é solicitado a localizar qualquer instância nas células realçadas e qualquer instância de parênteses será substituído por nada, assim, removê-los.
2. Determinantes enriquecido vias
   1. copiar os símbolos de gene do conjunto para o qual pretende determinar os caminhos enriquecidos para a área de transferência. Ir para ConsensusPathDB ¹². Selecione " Gene definir análise " na barra lateral esquerda da página, seguido de " análise de sobre-representação ".
   2. Cole a lista de genes no " colar uma lista de identificadores de gene/proteína " caixa. Símbolo do gene selecione a " tipo de identificador do Gene/proteína " caixa e clique em " continuar ". Marque a caixa ao lado de " caminhos definidos pelos bancos de dados via " sob a base de caminho define seção.
      Nota: Um número de opções para a análise aparecerá incluindo a lista de bancos de dados disponíveis para pesquisar. Recomenda-se de selecionar todos os bancos de dados, exceto os bancos de dados Kegg e Reactome, como eles são os mais estabelecida e abrangente. A sobreposição mínima com entrada configurações lista e p-valor de corte também pode ser ajustada com base nas necessidades. Recomendamos deixar essas configurações na sobreposição mínima de 2 gene padrão e um limite de valor 0,01 p.
   3. Select " encontrar enriquecido de moda " para obter a lista de caminhos que contêm genes na lista de entrada.
      Nota: A saída será em formato de tabela, incluindo o nome de cada percurso seguido por " definir tamanho " (indicando o número de genes totais na via), " candidatos contidos " (indicando o número de genes na lista de entrada que estão no caminho que ), bem como o valor-p e q-valor (FDR) e o banco de dados do qual foi identificado o pathway.
3. Geração de redes via
   1. determinar caminhos enriquecidos, como descrito acima.
   2. Todos os caminhos enriquecidos de Visualizar em uma rede de caminho selecionando cada caixa ao lado dos nomes de caminho para ser visualizado ou clicando em selecionar " todos os " sob " selecione " no cabeçalho da coluna acima as caixas de seleção, selecione " Visualizar conjuntos selecionados ".
      Nota: Recomendamos a visualização de todos os percursos se houver menos de 30; Se houver maior que 30 vias recomendamos selecionar as top 30 vias mais enriquecidas (aqueles que contêm o maior número de genes da lista de entrada).
   3. Ajustar o " sobreposição relativa " e " compartilhado candidatos " filtros, selecionando ou caixa no centro superior da página e introduzir desirpor cento de Ed relativo se sobrepõem ou número de candidatos compartilhados, para ser mais rigorosas a fim de tornar o mais pertinente sobrepõe dentro das redes de caminho mais claras e selecionadas " aplicar ".
      Nota: Recomendamos uma sobreposição mínima relativa de 0,2, indicando uma sobreposição de gene de 20% entre 2 caminhos para conectá-los e um mínimo de 2 candidatos compartilhados. A legenda do gráfico pode ser vista clicando-se na lenda de gráfico no canto superior esquerdo da página.
4. Determinação dos termos GO enriquecidos
   1. copiar os símbolos do gene do grupo para o qual deseja determinar as ontologias gene enriquecido para a área de transferência. Vá para o ' ferramenta de análise de enriquecimento GO ' no Gene Ontology consórcio ¹³. Cole a lista de símbolos de gene na caixa do lado esquerdo da página sob " seu gene IDs aqui … "
   2. selecionar qual conjunto de termos de ir usar, listados por baixo da caixa de IDs de gene: processo biológico, função molecular ou componente celular. Escolha (recomendado) ' processo biológico '. Selecione ' Danio rerio ' debaixo do GO termos caixa e clique em " enviar ".
      Nota: Recomendamos a utilização do corte de p-valor padrão de 0,05.

5. Verificação por qRT-PCR

Nota: genes individuais identificados com as mudanças de expressão de gene significativo no RNASeq devem ser verificados por qRT-PCR alvo em replicar experiências.

1. síntese do cDNA
   1. extrair o RNA de embriões utilizando o método descrito na secção 3.1. Extração do RNA.
   2. Converter 1 µ g de RNA de cDNA usando cDNA conversão kit (veja a Tabela de materiais). Combine o RNA, dNTPs e enzima transcriptase reversa. Realizar a reação de thermocycler de acordo com o fabricante ' especificações do s.
   3. Diluir cDNA 1:3 em água tratada DEPC. Loja a 4 ° C por até 1-2 meses.
2. verificação de qRT-PCR
   1. selecionar genes para verificar por qRT-PCR de resultados do sequenciamento de RNA. Identifica genes com alta prega mudanças na expressão. Determinar qual destes genes também contêm altas contagens de leitura, já que indica expressão abundante.
   2. Selecionar 10-12 alvos da mudança alta dobra, lista de alta contagem de leitura de genes. Desenho de primers para amplificar o gene-alvo que se estendem pelo menos 1 limite de intron-exon.
   3. Entra o gene ou região de destino do gene em qPCR cartilha projeto local ¹⁴. Selecione o " 2 primeiras demão qPCR de design " e solicitar o site para criar cartilha conjunto.
      Nota: Certifique-se de incluir iniciadores de controle interno, tais como β-actina, para normalizar as comparações entre as amostras. Os primers de β-actina são 5 f: ' - TCGAGCTGTCTTCCCATCCA - 3 ' e 5 r: ' - TCACCAACGTAGCTGTCTTTCTG - 3 '.
   4. Combinar o cDNA, encaminhar cartilha, cartilha reversa e polimerase. Executar a amplificação de qRT-PCR para determinar a expressão relativa dos genes ¹⁵.
   5. Analisar a expressão de mRNA relativo dos genes de interesse pela quantificação relativa determinação ¹⁵.
   6. Comparar o qRT-PCR para RNASeq resultados para garantir que a direcionalidade de expressão diferencial é consistente.

Results

Classificação de diferencialmente expressos Genes:

Para identificar genes diferencialmente expressos na fase larval de zebrafish modelos de síndrome de Alström e síndrome de Bardet-Biedl (BBS), escolhemos a qualquer alms1 ou bbs1 transcrições injetando validados anteriormente da tala-MOs em bloqueio zebrafish do selvagem-tipo embriões^16,17. Em 5 dias pós f...

Discussion

A abordagem descrita no presente protocolo oferece uma estratégia relativamente rápida e econômica para a análise do transcriptoma-nível de animais inteiros ou populações de célula específica classificada. O peixe-zebra fornece um modelo de vantajoso para este tipo de estudo devido a facilidade e a rapidez na geração de grandes quantidades de começar o material, a facilidade de aplicação genética ou ambientais condições experimentais e a disponibilidade de um grande espectro de linhas transgénicas repórter permitindo iso...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Reagentes comerciais
TriZol	Reagente de lise	de 15596026	Thermo Scientific
TrypLE	Gibco	12604013	tampão de dissociação 1
FACSMax	Genlantis	T200100	tampão de dissociação 2
Água tratada com DEPC	Sigma	95284
Kit de conversão de cDNA FirstStrand	Thermo Scientific	K1621	cDNA
2X SYBR Green Master Mix	Roche	4707516001	qRT-PCR Master Mix
FACS buffer	Fisher Scientific	50-105-9042
clorofórmio	Sigma Aldrich	288306
acetato de sódio	Sigma Aldrich	S2889
Nome	Empresa	Número de Catálogo	Comentários
Linhagens de Peixe-Zebra
Tuebingen	ZIRC	ZL57
ins2a:mCherry	ZIRC	ZL1483
Nome	< strong>Empresa	Número de catálogo	Comentários
Equipamento
40 filtro de células de mícron	Sigma	CLS431750
tubo FACS	BD Falcon	352063
hemocitômetro	Sigma	Z359629
Microscópio de dissecação	Zeiss
Microscópio invertido	Zeiss
Nanodrop	Thermo Scientific
Illumina HiSeq	Illumina
LightCycler 480	Roche
Tanques de acasalamento 1.0L Conjunto de tanques de cruzamento	Aquaneering	ZHCT100
tubo FACS 5 mL tubo de polipropileno	BD Falcon	352063
Nome	Empresa	Número de catálogo<	strong>Comentários
Software
Excel	Microsoft
Consensus Path DB	http://cpdb.molgen.mpg.de/
GO Análise de Enriquecimento	http://geneontology.org/page/go-enrichment-analysis