Todas as amostras utilizadas para determinar as grelhas de leitura aberta (ORF) selecionados neste estudo foram coletados de Starbuck Island, local 7 (STAR7) e Caroline Atoll, Sítio 9 (CAR9) da linha ilhas do sul. Estima-se 100 L de água do mar a partir destes locais foi recolhido abaixo as bombas de esgoto camada limite usando coral, como descrito anteriormente 5. Conteúdo das bombas estavam sujeitos a fragmentação através de grandes filtros de poros para remover pequenos eucariontes e posteriormente concentrada utilizando 100 kDa filtros de fluxo tangencial deixando apenas micróbios e viral como partículas (VLP). Para separar as VLP, a água do mar remanescente foi passado através de filtros de 0,45 um, resultando na virome. O clorofórmio foi introduzido para esta fracção viral para parar o crescimento de quaisquer células remanescentes e armazenado a 4 ° C.
As VLP foram purificadas utilizando o método de cloreto de césio, em que os gradientes de densidade de separar por centrifugação e permitem a recuperação de Vir~ iões a 1,35 g / ml a 1,5 g / ml 3. O ADN virai foi extraído utilizando um CTAB / fenol: clorofórmio protocolo e amplificado através de amplificação por deslocamento de múltiplos utilizando reagentes phi29. A sequenciação do virome foi realizado com a tecnologia disponível no mercado pyrosequencing.
Bioinformatics utilizados no processamento e selecção de ORFs virais para este estudo são as seguintes. Três etapas de pré-processamento foram utilizados nas metagenomes virais CAR9 e STAR7. Em primeiro lugar, o software público foi usado para remover sequências marcadoras que resultaram da amplificação do ADN viral antes da sequenciação 27. Em segundo lugar, artefatos de sequenciamento comuns, tais como duplicatas de seqüência e baixo número de cópias foram filtrados para fora do conjunto de dados através de um programa de bioinformática adicional 28. Por fim, a remoção da contaminação sequência externa 29 foi realizada para as sequências que tinham ≥ 90% de cobertura e ≥ 94% de identidade com sequências em bases de dados: v RefSeqgenomas irus; Humano - Referência GRCh37; Humano - Celera Genomics; Humano - Craig Venter (HuRef); Humano - Seong-Jin Kim (coreano); Humano - Chromosome 7 versão 2 (TCAG); e Humano - James Watson, Yanhuang (YH; asiática), Yoruba (NA18507; Africano) sequências de referência 21. Na sequência destes processos, sequências de amostras CAR9 totalizou 591.600 e sequências STAR7 totalizaram 939.311. Estas sequências foram enviados para MGRAST e montados via software assembler usando as configurações padrão. Contigs foram traduzidos em seis quadros de leitura e frames putativos de leitura aberta (pORFs) foram identificados usando scripts, como descrito anteriormente 21.
Para identificar ORF desconhecidos foram realizadas uma série de pesquisas de similaridade baseado remover ORFs de função conhecida. Resumidamente, os seguintes foram realizadas pesquisas com os seus critérios de pesquisa correspondentes 21:
- Similaridade significativa de ≥ 95% de identidade ao longo ≥40 pares de bases (pb) por MGRAST BLAT na base de dados M5NR.
- Similaridade significativa (e-valor ≤ 0,001) por TBLASTN contra todas metagenomes públicas em My Resource metagenoma banco de dados.
- Similaridade significativa (e-valor ≤ 0,001) por BLASTP e TBLASTN contra banco de dados NR.
- Similaridade significativa (e-valor ≤ 0,001) por RPS-BLAST contra o banco de dados Conserved Domain.
- Traduções de proteína a partir de uma fracção seleccionada de cada conjunto de dados de estruturas de proteínas em comparação contra resolvidos no Protein Data Bank.
- Cálculo das freqüências dinucleotidicos usando pacote Dinucleotide Assinaturas.
Os pORFs resultantes foram concebidos para expressão em E. coli usando software de design gene disponível publicamente. Back-tradução das sequências de amino-ácidos empregou uma tabela de uso Universal Codon projetado para acomodar a expressão em E. coli com um limiar mínimo a utilização de 2%. Sequência de reconhecimento da enzima de restriçãos para BamHI e HindIII foram excluídas das sequências para facilitar a clonagem. Uma empresa fora sintetizado o gene obtido por engenharia sequências de 30 e, em seguida, as ORFs foram clonados num vector de cópia de médio número promotor pBAD, pEMB11, através de clonagem de enzima de restrição padrão. Todos os clones foram transformados em E. coli K-12 estirpe BW 27784 23.
Placas Multi-fenótipo ensaio (mapas)
Um gasoduto de alto rendimento e software robusto foi aproveitado para análise dos mapas, o PMAnalyzer 24. O gasoduto foi desenvolvido em um ambiente de servidor Linux e executa várias etapas, incluindo: análise de arquivos de densidade óptica, a formatação de dados em arquivos de texto legível, pré-processamento de curvas de crescimento para a garantia de qualidade (QA), e realizando técnicas de modelagem matemática para analisar as curvas de crescimento . Os scripts modelagem primários foram desenvolvidos em Python versão 2.7.5 para fazer uso do módulo Pylab.
Mapas reprodutibilidade foi avaliada utilizando o erro padrão (SE) para dados replicados (Figura 2A). As curvas de crescimento foram comparadas com matérias de curvas de crescimento de logística para determinar se o programa PMAnalyzer crescimento parametrizado com precisão e modelado clone decurso da experimentação (dados não apresentados). Para mais informações sobre a exactidão ea validade dos mapas e PMAnalyzer ver Cuevas 24 et al.
Após a validação do método, mapas de dados foi analisada utilizando os múltiplos parâmetros, como taxa máxima de crescimento (μ max) e nível de crescimento (GL), fornecidos pelo pipeline de processamento. Visualização comparativa de curvas de crescimento é muitas vezes usado para interpretação de dados de crescimento; no entanto, o número de curvas que podem ser visualizados por vez para a comparação tem limitações. Para analisar numerosas curvas de crescimento simultaneamente, mapa de calor parcelas derivados foram implorou para comparar dezenas de clones cultivados em um único substrate contra a resposta média que 'condições (Figura 2B). A influência colocada por sobre-expressão de uma nova proteína de fago é observado através de mudanças nos parâmetros da curva, especificamente: lag fase, fase exponencial ea produção de biomassa máxima (asymptote). Como um exemplo, a subida íngreme da fase de latência em fase exponencial modelado na curva de crescimento para a proteína da cápside (Figura 2A) é reproduzida por uma rápida mudança na intensidade de cor de preto para branco para o mesmo clone na trama dinâmica da Figura 2B .
Para obter uma imagem global de distribuição clone através substratos, classificações fenotípicas derivados do GL foram usadas (Figura 3). Aqui, as quatro fenótipos são separados em quatro gráficos onde a altura de cada barra representa o número de clones que exibem o fenótipo de um substrato específico. Outliers nos dados são reconhecidos como clones que caem no "ganho de function "ou" perda de função categoria ". Outliers pode, então, ser procurados individualmente e investigado mais de perto experimentalmente. Além disso, reconhece polarizações análise global de substrato no ensaio. Os substratos tais como fenilalanina, ácido málico, glicina e resultou num "nenhum crescimento" classificação. Substratos que constantemente caem na nenhuma classificação de crescimento, em todos os clones, não são fortemente ponderada em caracterização funcional jusante.
Metabolomics
Os produtos catabólicos de clones que expressam genes de fago desconhecidos foram identificados utilizando metabolômica. Resumidamente, os clones foram cultivadas sob quer uma cultura contínua ou passagem em série em meio de cultura de lote antes de ser enviado para análise GC-TOFMS em uma instalação do núcleo metabolômica. Para mais detalhes sobre o processamento da amostra, análise e normalização para GC-TOFMS implementadas pela instalação de núcleo escolhida ver Fiehn et al. 31 Briefly, 1 ml de solvente de extracção frio são adicionados a cada amostra, após o que as amostras são agitadas e sonicada num banho frio durante 5 minutos. As amostras são centrifugadas e, finalmente, a metade da amostra é decantado e secou-se para baixo para análise. Os extractos são purificados e com adição de marcadores de índice de retenção interna antes de ser carregado para o cromatógrafo de fase gasosa e, em seguida, subsequentemente transferido para o espectrómetro de massa. Os dados de cada amostra é analisada de modo a que as intensidades de sinal de todos os sinais detectados no cromatograma são relatados. Para a normalização, a abundância de picos para cada uma das amostras é resumida e as abundâncias totais dos picos são em média entre todas as amostras no conjunto. Abundâncias de metabolitos por amostra são divididos por pico de abundância da amostra e, em seguida, multiplicado pelo pico de abundância média do conjunto da amostra. Os dados resultantes são utilizados para a análise do metaboloma na pesquisa discutidos.
Validação de metabolômica reprodutibilidade foi obrigado adeterminar o tamanho da amostra adequado para cada método de cultivo. Para detectar a precisão visto dentro das amostras e a variação observada entre tamanhos de amostra o erro padrão da média (M S), tanto para n = 3 e n = 6 conjuntos de dados foram analisados (Figura 5B). Independentemente do tamanho da amostra de cultura contínua (CC), menos de 1% dos dados teve um S M ≤ 1,5. Mediana S M S foram 221 e 300, e os valores variaram de 0-7,55 x 10 5 3,74 x 10 e 5, para n = 3 e n = 6, respectivamente. S S M foram também calculados para cada conjunto de amostra se replica no método de cultura de série (SC). Mais uma vez, menos de 1% dos dados teve um S M ≤ 1,5, um número médio de S M 137, e uma gama de 0-3,51 x 10 5. Para comparação das distribuições S m entre cada conjunto de amostras (n = 3 CC vs . CC n = 6, CC N = 3 x SC n = 3, e CC n = 6 versus SC n = 3), foi realizado um teste de permutação. O distribuição de cultura contínua ou S dos valores m conjunto de dados não era significativamente diferente da distribuição dos valores de cultura de série S M (valor de p = 0,0). No entanto, a distribuição de valores S M para cultura contínua de dados n = 3 era significativamente diferente da dos valores de S M para a cultura contínua n = 6 dados (valor de p = 1,908804 x 10 -49). Finalmente, o coeficiente de variação por metabolito foi comparado antes e depois da execução de uma etapa de garantia de qualidade (QA) (Figura 5C). No total, 210 metabolitos foram removidos após a implementação do gasoduto QA (40% dos dados). Menos de 1% dos dados removidas tinham uma abundância metabolito de zero, ~ 2% foi dados do padrão interno, ~ 5% foi a partir de dados de metabolitos nunca antes observados em E. coli, e os metabolitos restantes (> 30%) teve um coeficiente de variação maior do que 1.
Tal como acontece com a análise mapas, observatio mundialns fornecida uma compreensão inicial da profundidade da informação metabolômica ofertas. Para obter uma imagem global, os clones foram hierarquicamente agrupadas com base em suas abundâncias de metabolitos relativos fornecendo informações sobre perfis de clone-metabolito, clones potenciais com funções relacionadas, e de outliers clone-metabolito (Figura 6). Para destacar as funções das proteínas, metabólitos são separadas e agrupadas com base em vias metabólicas comuns. Utilizando esta análise com resultados preliminares, foi evidente que metaboloma é capaz de separar os genes de classes diferentes (Figura 6, em destaque clones). Além disso, com a identificação de aberrações dados metaboloma foi determinada calculando pontuações normais (pontuações z) para cada par clone de metabolito. Para garantir significado estatístico, valores aberrantes foram definidos como um par clone-metabolito, com um valor de pontuação Z 2, correspondendo a apenas 5 por cento dos dados (dados não mostrados).
ether.within-page = "always"> 
Figura 1. Definições de classificações fenótipo. (A) Relação entre o nível de crescimento (GL) e taxa máxima de crescimento. Os pontos de dados circulado em vermelho representam as curvas de crescimento que mostram pouca ou nenhuma utilização do substrato. (B) representação Boxplot definição de limiar de crescimento com base na distribuição de curvas de crescimento com uma taxa de crescimento mínima (<0,15 OD / hr). (C) A variância eo desvio padrão de GL é calculado para substrato de D-galactose. Linhas tracejadas curtas representam dois desvios padrão de distância a partir da média. Por favor clique aqui para ver uma versão maior desta figura.
highres.jpg "/>
Figura 2. Os mapas de validação por meio de precisão e diferenciação. (A) curvas de crescimento de clones estruturais (cápside) e metabólicos (anotados tioredoxina), dois novos clones metabólicas (EDT2440, EDT2441), e a resposta média de clones cultivados em sacarose, D- galactose e D-manose nos mapas. As linhas azuis indicam o erro padrão visto entre dados replicados (n = 3). (B) As curvas de crescimento para 47 clones diferentes são representadas como mapas de calor para sacarose, D-galactose e D-manose. A estrutural (círculo verde) anotada e clones metabólicas (círculo laranja), dois clones metabólicas novos (círculos azuis claros e escuros), ea resposta média (círculo vermelho) são destacados. Por favor clique aqui para ver uma versão maior desta figura.
pload / 52.854 / 52854fig3highres.jpg "width =" 700 "/>
Figura 3. Distribuição Clone para cada fenótipo em vários substratos. A contagem de fenótipo-clone para 47 clones em 72 condições de crescimento específicas de carbono. Tabela fornece contagens directas para cada fenótipo. Por favor clique aqui para ver uma versão maior desta figura.

Figura 4. Diagrama que detalha a construção do aparelho de cultura contínua. (A) passos utilizados para construir portas α-γ do reactor cultura contínua, (B) etapas utilizadas para construir a porta de fluxo para fora do reactor cultura contínua, e (C ) Os passos para construir portos δ e ε da mamadeira cultura contínua. = "Https://www.jove.com/files/ftp_upload/52854/52854fig4highres.jpg" target = "_ blank"> Clique aqui para ver uma versão maior desta figura.

Figura 5. Comparação dos métodos apresentados Phenomic. (A) do fluxo de trabalho para a preparação das placas multi-fenótipo de ensaio (MAPs), culturas contínuas e de culturas em série. (B) O percentual de erro padrão da média (S M) conta tanto para n = 3 e n = 6 tamanhos de amostra para a cultura contínua (CC) e da cultura de série (SC) métodos de preparação para metabolômica. O eixo y está numa escala logarítmica. (C) As distribuições de coeficientes de variação (CV) por metabolito, antes e depois da implantação do gasoduto GQ para o método de cultura contínua (CC)."target =" _ blank g5highres.jpg "> Clique aqui para ver uma versão maior desta figura.

Figura 6. Perfis metabolômico de clones crescidos em cultura contínua. As abundâncias metabolito mediana de um conjunto de metabolitos são representados por 84 clones crescidos em culturas contínuas. Perfis de metabólitos de clones anotados estruturais (capsídeo) e metabólicos (thioredoxin), dois novos clones metabólicas (EDT2440, EDT2441), ea resposta metabólica média são destacadas em vermelho. Por favor clique aqui para ver uma versão maior desta figura.
| Composto | Carbono | Azoto | Enxofre | Fósforo |
| Glicerina | - | 0,40% | 0,40% | 0,40% |
| Cloreto de amônio | 9,5 mM | - | 9,5 mM | 9,5 mM |
| O sulfato de sódio | 0,250 mM | 0,250 mM | - | 0,250 mM |
| O sulfato de magnésio | 1,0 mM | 1,0 mM | - | 1,0 mM |
| Fosfato de potássio | 1,32 mM | 1,32 mM | 1,32 mM | - |
| O cloreto de magnésio | - | - | * | - |
| Cloreto de potássio | 10 mM | 10 mM | 10 mM | 10 mM |
| Cloreto de cálcio | 0,5 uM | 0,5 uM | 0,5 uM | 0,5 uM |
| Cloreto de sódio | 5 mM | 5 mM | 5 mM | 5 mM |
| Cloreto férrico | 6? M | 6? M | 6? M | 6? M |
| L- arabinose | 0,10% | 0,10% | 0,10% | 0,10% |
| MOPS pH 7,4 | 1x | 1x | 1x | 1x |
Tabela 1. Os compostos e concentrações dos diferentes meios de comunicação basais utilizados nos mapas. * 1,0 mM de cloreto de magnésio é substituído. 1x MOPS = 40 mM de MOPS, Tricina 4 mm.
| Substratos de Carbono | Substratos de nitrogênio | Substratos de enxofre | Psubstratos hosphorus |
| 2 desoxi-D-ribose | 2-desoxi-D-ribose | Ácido 1-butano-sulfónico | adenosina-5-monophoshate |
| Ácido acético 4-hidroxi-fenil | acetamida | acetilcisteína | beta-glicerofosfato |
| ácido acético | adenina | D-cisteína | creatinephosphate |
| adenosina-5-monofosfato | adenosina | D-metionina | D-glicose-6-fosfato |
| adonitol | alantoína | dietil-ditiofosfato | dietil-ditiofosfato |
| alfa-D-glicose | beta-feniletilamina | DL-etionina | DL-alfa-glicerofosfato |
| alfa-D-lactose | biureto | glutationa | fosfato de potássio |
| alfa-D-melebiose | cytidine | ácido isetiónico | pirofosfato de sódio |
| ácido cítrico | citosina | Ácido L-cisteico | tiofosfato de sódio |
| D-alanina | D-alanina | L-cisteína | |
| D-arabinose | D-asparagina | Ácido L-djenkolic |
| D-arabitol | D-aspartato | L-metionina |
| D-asparagina | D-cisteína | sulfato de magnésio |
| D-aspartato | D-glucosamina | ácido sulfónico metano |
| D-cellubiose | Ácido D-glutâmico | N-acetil-DL-metionina |
| D-cisteína | Ácido DL-alf a-amino-n-butírico | N-acetil-L-cisteína |
| D-frutose | D-metionina | de potássio-tetra-thionate |
| D-galactose | D-serina | tiossulfato de sódio |
| D-glucosamina | D-valina | ácido sulfanic |
| D-glicose | ácido gama-amino-N-butírico | taurina |
| D-glicose-6-fosfato | glicina | ácido taurocólico |
| D-glutamato | guanidina | tioureia |
| D-manose | histamina | |
| D-rafinose | inosine |
| D-ribose | L-alanina |
| D-salicin | L-arginina |
| D-serina | L-asparagina |
| D-trealose | L-citrulina |
| D-xilose | L-cisteína |
| dulcitol | Ácido L-glutâmico |
| glicerina | L-glutamina |
| glicina | L-glutationa |
| i-eritritol | L-histidina |
| inosine | L-isoleucina |
| L-alanina | L-leucina |
| L-arabinose | L-lisina |
| L-arabitol | L-metionina |
| L-asparagina | L-ornitina |
| L-aspartato | L-fenil-alanina |
| Ácido L-cisteico | L-prolina |
| L-cisteína | Ácido L-piro-glutico |
| L-fucose | L-serina; L-treonina |
| Ácido L-glutâmico | L-triptofano |
| L-glutamina | L-valina |
| L-isoleucina | N-acetil-D-glucosamina |
| L-leucina | putrescina |
| L-lisina | tioureia |
| L-metionina | timidina |
| L-fenilalanina | timina |
| Ácido L-piro-glutico | tiramina |
| L-ramnose | tirosina |
| L-serina | uridine |
| L-sorbose | |
| L-treonina |
| L-triptofano |
| L-valina |
| L-xilose |
| lactato |
| lactulose |
| malato |
| mio-inositol |
| ácido oxálico |
| sorbato de potássio |
| ácido propiónico |
| putrescina |
| ácido quínico |
| piruvato de sódio |
| succinato de sódio |
| sacarose |
| timidina |
| xilitol |
Tabela 2. Lista de substratos utilizados nos experimentos MAP.