Tous les échantillons utilisés pour déterminer les cadres de lecture ouverts (ORF) sélectionnés dans cette étude ont été recueillies à partir de l'île de Starbuck, site 7 (STAR7) et Caroline Atoll, Site 9 (CAR9) des îles de la Ligne du sud. On estime à 100 L d'eau de mer à partir de ces sites ont été recueillis ci-dessous les coraux couche limite à l'aide de pompes de cale, comme décrit précédemment 5. Sommaire des pompes étaient soumis à fragmentation par les grands filtres à pores pour éliminer les petites eucaryotes et ensuite concentrées en utilisant 100 kDa filtres de flux tangentiel ne laissant que les microbes et les virus comme des particules (VLP). Pour séparer les VLP, restant l'eau de mer a été adoptée de 0,45 um filtres résultant dans le virome. Le chloroforme a été introduit à cette fraction virale pour interrompre la croissance de toutes les cellules restantes et stocké à 4 ° C.
VLP ont été purifiées en utilisant le procédé au chlorure de césium, dans laquelle des gradients de densité séparés par centrifugation et permettent la récupération de vir~ ions à 1,35 g / ml à 1,5 g / ml 3. L'ADN viral a été extrait en utilisant un CTAB / phénol: protocole de chloroforme et amplifié par une amplification de déplacement multiples utilisant des réactifs phi29. Le séquençage de l'virome a été réalisée avec la technologie disponible dans le commerce pyroséquençage.
Bioinformatique utilisés dans le traitement et la sélection des ORF viraux pour cette étude sont les suivantes. Trois étapes de pré-traitement ont été utilisés sur les métagénomes virales CAR9 et STAR7. Tout d'abord, le logiciel du public a été utilisée pour éliminer les séquences d'étiquette résultant de l'amplification de l'ADN viral 27 avant le séquençage. Deuxièmement, les artefacts de séquençage communs tels que des doublons de séquence et faible nombre de copies ont été filtrées hors de l'ensemble de données par l'intermédiaire d'un programme de bioinformatique 28 supplémentaire. Enfin, l'élimination de la contamination de séquence étrangère 29 a été réalisée pour les séquences qui avaient ≥ 90% de couverture et ≥ 94% d'identité avec des séquences dans les bases de données suivantes: RefSeq vgénomes irus; Human - Référence GRCh37; Human - Celera Genomics; Human - Craig Venter (HuRef); Human - Seong-Jin Kim (Corée); Human - Chromosome 7 version 2 (GACT); et humain - James Watson, Yanhuang (YH; Asie), Yoruba (NA18507; africaine) séquences de référence 21. Suite à ces processus, des séquences à partir des échantillons CAR9 totalisé 591 600 et 939 311 séquences ont totalisé STAR7. Ces séquences ont été téléchargés sur MGRAST et assemblés par le logiciel assembleur en utilisant les paramètres par défaut. Contigs ont été traduites en 6 cadres de lecture et putatifs cadres de lecture ouverts (pORFs) ont été identifiés en utilisant des scripts, comme décrit précédemment 21.
Pour identifier ORF inconnus un certain nombre de recherches de similarité base ont été effectuées pour retirer ORF de fonction connue. En bref, les recherches suivantes ont été réalisées avec leurs critères de Recherche en 21 correspondants:
- Similarité significative de ≥ 95% d'identité sur ≥40 paires de base (pb) par MGRAST BLAT dans la base M5NR.
- Similarité significative (valeur e ≤ 0,001) par TBLASTN contre toutes métagénomes publics dans Mon ressources de base de données de Metagénome.
- Similarité significative (valeur e ≤ 0,001) par BLASTP et TBLASTN contre NR base de données.
- Similarité significative (valeur e ≤ 0,001) par RPS-BLAST contre la base de données Conserves de domaine.
- traduction de protéines à partir d'une fraction de sélection de chaque ensemble de données comparés à ceux de structures de protéines résolues dans la Protein Data Bank.
- Calcul des fréquences dinucléotidiques utilisant package Dinucleotide Signatures.
Les pORFs résultants ont été conçus pour l'expression dans E. coli en utilisant la disposition du public un logiciel de conception de gène. Retour-traduction des séquences amino-acide utilisé une utilisation Universal Table Codon conçu pour accueillir l'expression dans E. coli avec un seuil minimum d'utilisation de 2%. séquence de reconnaissance d'enzyme de restrictions pour BamHI et HindIII ont été exclus des séquences pour faciliter le clonage. Une entreprise de l'extérieur synthétisé le gène modifié séquences 30 et puis ORF ont été clonées dans un vecteur de promoteur numéro de pBAD moyen copie, pEMB11, via le clonage d'enzyme de restriction standard. Tous les clones ont été transformés en E. coli souche K-12 BW 27784 23.
Multi-phénotype Assay Plaques (MAP)
Un pipeline à haut débit et un logiciel robuste a été mobilisé pour l'analyse des cartes, l'PMAnalyzer 24. Le pipeline a été développé dans un environnement de serveur Linux et effectue diverses étapes, notamment: analyse de fichiers de densité optique, le formatage des données dans des fichiers textes, de pré-traitement des courbes de croissance pour l'assurance qualité (AQ), et de la scène de techniques de modélisation mathématique pour analyser les courbes de croissance . Les scripts de modélisation primaires ont été développés dans la version Python 2.7.5 de faire usage du module PyLab.
Cartes reproductibilité a été évaluée en utilisant l'erreur standard (SE) pour les données répliquées (figure 2A). Matières premières courbes de croissance ont été comparées à des courbes de croissance de logistique afin de déterminer si le programme PMAnalyzer croissance de clones avec précision paramétrée au cours de l'expérimentation et modélisées (données non présentées). Pour de plus amples informations sur l'exactitude et la validité des cartes et PMAnalyzer voir Cuevas et al. 24
Après la validation de la méthode, les données ont été analysées MAP en utilisant les plusieurs paramètres, tels que la vitesse maximale de croissance (μ max) et le niveau de croissance (GL), fourni par le pipeline de traitement. Visualisation comparative des courbes de croissance est souvent utilisé pour l'interprétation des données de la croissance; toutefois, le nombre de courbes qui peut être visualisée à la fois pour la comparaison a ses limites. Pour analyser de nombreuses courbes de croissance simultanément, carte de chaleur parcelles dérivées ont été imploré de comparer des dizaines de clones cultivés sur une seule substratee à l'encontre réponse moyen que les conditions de (figure 2B). L'influence placé par la surexpression d'un roman protéine de phage est observé à travers les changements dans les paramètres de la courbe, spécifiquement: phase de latence, la phase exponentielle et le rendement de la biomasse maximum (asymptote). A titre d'exemple, la montée raide de la phase de latence en phase exponentielle modélisé dans la courbe de croissance pour la protéine de capside (figure 2A) est reproduit par un changement rapide de l'intensité de la couleur du noir au blanc pour le même clone dans le complot dynamique de la figure 2B .
Pour obtenir une image globale de la distribution de clone travers substrats, classifications phénotypiques provenant de la GL ont été utilisés (Figure 3). Ici, les quatre phénotypes sont séparés en quatre tableaux où la hauteur de chaque barre représente le nombre de clones présentant ce phénotype pour un substrat spécifique. Les valeurs aberrantes dans les données sont reconnues comme des clones de tomber dans le «gain de function "ou" perte de la catégorie de fonction ". Les valeurs aberrantes peuvent ensuite être recherchées individuellement et étudiés de plus près expérimentalement. En outre, l'analyse globale reconnaît biais de substrat dans le dosage. Des substrats tels que la phenylalanine, l'acide malique, la glycine et ont donné lieu à une "absence de croissance" classification. Substrats qui tombent toujours dans le pas de classification de la croissance, dans tous les clones, ne sont pas fortement pondérés dans la caractérisation fonctionnelle en aval.
Métabolomique
Les produits cataboliques de clones exprimant des gènes de phages inconnus ont été identifiés en utilisant la métabolomique. En bref, les clones ont été cultivés sous une culture en continu, soit un passage en série ou en lot milieu de culture avant d'être envoyé à l'analyse GC-TOFMS à un centre de noyau de la métabolomique. Pour plus de détails sur le traitement des échantillons, l'analyse et la normalisation pour GC-TOFMS mises en œuvre par l'installation de base choisie voir Fiehn et al. 31 BrieFly, 1 ml de solvant d'extraction à froid est ajouté à chaque échantillon, après quoi les échantillons sont vortexés et soniquées dans un bain froid pendant 5 min. Les échantillons sont centrifugés et enfin la moitié de l'échantillon est décanté et séché vers le bas pour analyse. Extraits sont purifiés et dopés avec des marqueurs d'index de rétention interne avant d'être chargé sur le chromatographe en phase gazeuse et ensuite transférée vers le spectromètre de masse. Les données de chaque échantillon sont analysées de telle sorte que les intensités de signal de tous les signaux détectés dans le chromatogramme sont rapportés. Pour la normalisation, l'abondance des pics pour chaque échantillon est résumée et le total des abondances de pointe sont en moyenne entre tous les échantillons dans l'ensemble. Abondances métabolites par échantillon sont divisés par le pic de l'abondance de l'échantillon, puis multipliées par l'abondance de pointe moyenne de l'ensemble de l'échantillon. Les données résultantes sont utilisées pour l'analyse du métabolome de la recherche examinée.
Validation de la métabolomique reproductibilité a été nécessaire pourdéterminer la taille appropriée de l'échantillon pour chaque méthode de culture. Pour détecter la précision vu dans les échantillons et la variation vu à travers des tailles d'échantillon de l'erreur standard de la moyenne (S M) à la fois pour n = 3 et n = 6 ensembles de données ont été examinés (figure 5B). Quelle que soit la culture en continu (CC) la taille de l'échantillon, moins de 1% des données avait un M S ≤ 1,5. Median S M s étaient 221 et 300, et les valeurs variaient de 0 à 7,55 x 10 5 et 3,74 x 10 5, pour n = 3 et n = 6, respectivement. S M s ont également été calculées pour chaque ensemble de l'échantillon répétitions dans la méthode de la culture de série (SC). Encore une fois, moins de 1% des données a eu un S M ≤ 1,5, une médiane S M de 137, et une gamme de 0 à 3,51 x 10 5. Pour comparer les distributions S M entre chaque ensemble de l'échantillon (CC n = 3 vs . CC n = 6, CC n = 3 vs SC n = 3, et CC n = 6 vs SC n = 3) un test de permutation a été effectuée. La distribution soit culture continue S M valeurs données n'a pas été significativement différente de la distribution de la culture série S M valeurs (valeur p = 0,0). Cependant, la distribution des valeurs S M pour la culture en continu des données n = 3 était significativement différente de celle des valeurs S M pour la culture en continu des données n = 6 (p = 1,908804 x 10 -49). Enfin, le coefficient de variation par métabolite a été comparée avant et après la mise en œuvre d'une étape d'assurance qualité (QA) (Figure 5C). Au total, 210 métabolites ont été éliminés après application de la canalisation de QA (40% des données). Moins de 1% de l'effacement des données avait une abondance de métabolite de zéro, ~ 2% était données standard internes, ~ 5% a été données de métabolites jamais observée dans E. coli, et les métabolites restants (> 30%) avaient un coefficient de variation supérieur à 1.
Comme pour l'analyse des cartes, observatio mondialens fourni une compréhension initiale de la profondeur de la métabolomique d'information offres. Pour obtenir une image globale, les clones ont été regroupés hiérarchiquement en fonction de leurs abondances relatives de métabolites fournissant des informations sur les profils de clones-métabolite, clones potentiels avec des fonctions connexes, et les valeurs aberrantes de clones-métabolite (figure 6). Pour mettre en évidence les fonctions des protéines, métabolites sont séparées et regroupées en fonction des voies métaboliques communes. En utilisant cette analyse avec les résultats préliminaires, il était évident que la métabolomique est capable de séparer les gènes de différentes classes (figure 6, mis en évidence clones). En outre, l'identification des valeurs aberrantes avec les données de métabolomique a été déterminée par le calcul des scores standard (z scores) pour chaque paire de clone-métabolite. Pour assurer la signification statistique, les valeurs aberrantes ont été définis comme une paire de clone-métabolite avec une valeur de Z score de 2, ce qui représente seulement 5 pour cent des données (données non présentées).
ether.within-page = "always"> 
Figure 1. Définitions des classifications de phénotype. (A) Rapport entre le niveau de croissance (GL) et le taux de croissance maximal. Les points de données entourées en rouge représentent les courbes de croissance montrant peu ou pas de l'utilisation du substrat. (B) Représentation Boxplot définissant le seuil de croissance basée sur la distribution des courbes de croissance avec un taux de croissance minime (<0,15 OD / h). (C) La variance et l'écart type est calculé pour GL substrat D-galactose. Pointillés courts représentent deux écarts types de la moyenne. S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.
highres.jpg "/>
Figure 2. Cartes validation par la précision et la différenciation. (A) Les courbes de croissance pour les clones structurelles (capside) et métaboliques (annotés thiorédoxine), deux nouveaux clones métaboliques (EDT2440, EDT2441), et la réponse moyenne de clones cultivés sur le saccharose, D- galactose et D-mannose dans les cartes. Les lignes bleues indiquent l'erreur-type observée entre données répliquées (n = 3). (B) Les courbes de croissance pour les 47 clones différents sont représentés sur les localisations de chaleur pour le saccharose, le D-galactose et D-mannose. La structure (cercle vert) annoté et clones métaboliques (cercle orange), deux clones métaboliques nouvelles (cercles bleus sombres et claires), et la réponse moyenne (cercle rouge) sont mis en évidence. S'il vous plaît, cliquez ici pour voir une version plus grande de cette figure.
pload / 52854 / 52854fig3highres.jpg "width =" 700 "/>
Figure 3. Répartition de Clone pour chaque phénotype sur plusieurs substrats. Le phénotype-clone comptent pour 47 clones dans 72 des conditions de croissance spécifiques carbone. Tableau fournit les chiffres directs pour chaque phénotype. S'il vous plaît cliquez ici pour voir une version plus grande de cette figure.

Figure 4. Schéma décrivant la construction de l'appareil de culture continue. (A) Les étapes utilisées pour construire les ports α-γ du réacteur de culture continue, (B) les étapes utilisées pour construire l'orifice d'écoulement du réacteur de culture en continu, et (C ) les mesures pour construire des ports δ et ε de la bouteille culture de l'alimentation continue. = "Https://www.jove.com/files/ftp_upload/52854/52854fig4highres.jpg" target = "_ blank"> S'il vous plaît cliquez ici pour voir une version plus grande de cette figure.

Figure 5. Comparaison des méthodes Phenomic présentés. (A) Flux de travail pour la préparation des multi-phénotype Assay Plaques (MAP), des cultures et des cultures continues de série. (B) Le pourcentage d'erreur standard de la moyenne (S M) compte à la fois pour n = 3 et n = 6 tailles d'échantillon pour la culture en continu (CC) et la culture de série (SC) les méthodes de préparation pour la métabolomique. L'axe des ordonnées est sur une échelle logarithmique. (C) Les distributions des coefficients de variation (CV) par un métabolite, avant et après la mise en œuvre de la canalisation d'assurance qualité pour le procédé de culture en continu (CC).g5highres.jpg "target =" _ blank "> S'il vous plaît cliquez ici pour voir une version plus grande de cette figure.

Figure 6. profils métabolomique de clones cultivés en culture continue. Les abondances métabolites médian pour un ensemble de métabolites sont tracées pour 84 clones cultivés dans des cultures continues. Profils métabolites pour les clones annotés structurelles (capside) et métaboliques (thiorédoxine), deux nouveaux clones métaboliques (EDT2440, EDT2441), et la réponse métabolique moyenne sont surlignés en rouge. S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.
| Composé | Carbone | Azote | Soufre | Phosphore |
| Glycérol | - | 0,40% | 0,40% | 0,40% |
| Le chlorure d'ammonium | 9,5 mM | - | 9,5 mM | 9,5 mM |
| Le sulfate de sodium | 0,250 mM | 0,250 mM | - | 0,250 mM |
| Le sulfate de magnésium | 1,0 mM | 1,0 mM | - | 1,0 mM |
| phosphate de potassium | 1,32 mM | 1,32 mM | 1,32 mM | - |
| Le chlorure de magnésium | - | - | * | - |
| Le chlorure de potassium | 10 mM | 10 mM | 10 mM | 10 mM |
| Le chlorure de calcium | 0,5 pM | 0,5 pM | 0,5 pM | 0,5 pM |
| Le chlorure de sodium | 5 mM | 5 mM | 5 mM | 5 mM |
| Chlorure ferrique | 6 pM | 6 pM | 6 pM | 6 pM |
| L- arabinose | 0,10% | 0,10% | 0,10% | 0,10% |
| MOPS pH 7,4 | 1x | 1x | 1x | 1x |
Tableau 1. Les composés et les concentrations des différents milieux de base utilisés dans les cartes. * 1,0 mM de chlorure de magnésium est substitué. MOPS 1x = 40 mM de MOPS, 4 mM de Tricine.
| des substrats de carbone | substrats d'azote | substrats de soufre | Psubstrats hosphorus |
| 2 désoxy-D-ribose | 2-désoxy-D-ribose | L'acide butane-1-sulfonique | l'adénosine-5-monophoshate |
| Acide acétique hydroxy-4 phényl | acétamide | acétylcystéine | bêta-glycérophosphate |
| acide acétique | adénine | D-cystéine | creatinephosphate |
| l'adénosine-monophosphate 5- | adénosine | D-méthionine | D-glucose-6-phosphate |
| adonitol | allantoïne | diéthyl-dithiophosphate | diéthyl-dithiophosphate |
| alpha-D-glucose | la bêta-phényléthylamine | DL-éthionine | DL-alpha-glycérophosphate |
| alpha-D-lactose | biuret | glutathion | phosphate de potassium |
| alpha-D-melebiose | cytidine | l'acide iséthionique | le pyrophosphate de sodium |
| acide citrique | cytosine | l'acide L-cystéique | thiophosphate de sodium |
| D-alanine | D-alanine | L-cystéine | |
| D-arabinose | D-asparagine | l'acide L-djenkolique |
| D-arabitol | D-aspartate | L-méthionine |
| D-asparagine | D-cystéine | du sulfate de magnésium |
| D-aspartate | D-glucosamine | l'acide méthane-sulfonique |
| D-cellobiose | l'acide D-glutamique | N-acétyl-DL-méthionine |
| D-cystéine | Acide DL-alpha-amino-n-butyrique | N-acétyl-L-cystéine |
| D-fructose | D-méthionine | potassium-tétra-thionate |
| D-galactose | D-sérine | du thiosulfate de sodium |
| D-glucosamine | D-valine | l'acide sulfanic |
| D-glucose | l'acide gamma-amino-N-butyrique | taurine |
| D-glucose-6-phosphate | glycine | l'acide taurocholique |
| D-glutamate | guanidine | thiourée |
| D-mannose | histamine | |
| D-raffinose | inosine |
| D-ribose | L-alanine |
| D-salicine | L-arginine |
| D-sérine | L-asparagine |
| D-tréhalose | L-citrulline |
| D-xylose | L-cystéine |
| dulcitol | l'acide L-glutamique |
| glycérol | L-glutamine |
| glycine | L-glutathion |
| i-érythritol | L-histidine |
| inosine | L-isoleucine |
| L-alanine | L-leucine |
| L-arabinose | L-lysine |
| L-arabitol | L-méthionine |
| L-asparagine | L-ornithine |
| L-aspartate | L-phényl-alanine |
| l'acide L-cystéique | L-proline |
| L-cystéine | L'acide L-pyro-glutamique |
| L-fucose | L-sérine; L-thréonine |
| Acide L-glutamique | L-tryptophane |
| L-glutamine | L-valine |
| L-isoleucine | N-acétyl-D-glucosamine |
| L-leucine | putrescine |
| L-lysine | thiourée |
| L-méthionine | thymidine |
| L-phénylalanine | thymine |
| L'acide L-pyro-glutamique | tyramine |
| L-rhamnose | tyrosine |
| L-sérine | uridine |
| L-sorbose | |
| L-thréonine |
| L-tryptophane |
| L-valine |
| L-xylose |
| lactate |
| lactulose |
| malate |
| myo-inositol |
| l'acide oxalique |
| sorbate de potassium |
| l'acide propionique |
| putrescine |
| l'acide quinique |
| pyruvate de sodium |
| le succinate de sodium |
| saccharose |
| thymidine |
| xylitol |
Tableau 2. Liste des substrats utilisés dans les expériences du PAM.