Alle monsters gebruikt om de open reading frames (ORF) is geselecteerd in dit onderzoek vast te stellen werden verzameld van Starbuck Island, Site 7 (STAR7) en Caroline Atoll, Site 9 (CAR9) van de zuidelijke Line-eilanden. Naar schatting 100 L van zeewater van deze sites werd onder het koraal grenslaag met behulp van lenspompen verzameld, zoals eerder 5 beschreven. Inhoud van de pompen werden onderworpen aan fragmentatie met grote poriën filters om kleine eukaryoten te verwijderen en vervolgens geconcentreerd onder toepassing van 100 kDa tangentiële stroming filters waardoor alleen microben en virale deeltjes (VLP's). Om VLP scheiden, werd de resterende zeewater geleid door 0,45 urn filters waardoor de virome. Chloroform werd aan deze virale fractie groei van resterende cellen arrestatie en bewaard bij 4 ° C.
VLP's werden gezuiverd met de cesiumchloridewerkwijze, waarbij dichtheidsgradiënten scheiden door middel van centrifugatie en het oog op het herstel van virionen bij ~ 1,35 g / ml tot 1,5 g / ml 3. Viraal DNA werd geëxtraheerd onder toepassing van een CTAB / fenol: chloroform protocol en geamplificeerd via meerdere verdringing amplificatie middels Phi29 reagentia. Sequencing van de virome werd bewerkstelligd met commercieel verkrijgbaar pyrosequencing technologie.
Bioinformatica die bij de verwerking en selectie van virale ORF's voor deze studie zijn de volgende. Drie stappen pre-processing werden gebruikt op de CAR9 en STAR7 virale metagenomes. Eerst werd openbaar software waarmee tag sequenties die het gevolg van amplificatie van het virale DNA voorafgaand aan sequencing 27 verwijderen. Ten tweede, gemeenschappelijke sequencing artefacten zoals opeenvolging duplicaten en lage aantal kopieën werden gefilterd uit de dataset via een extra bioinformatica-programma 28. Tot slot, het verwijderen van buitenlandse opeenvolging verontreiniging 29 werd uitgevoerd voor die sequenties die ≥ 90% dekking en ≥ 94% identiteit met sequenties in de volgende databases had: RefSeq virus genomen; Human - Reference GRCh37; Human - Celera Genomics; Human - Craig Venter (HuRef); Human - Seong-Jin Kim (Koreaans); Human - Chromosome 7 versie 2 (TCAG); en Human - James Watson, Yanhuang (YH; Aziatische), Yoruba (NA18507; Afrikaanse) referentie-sequenties 21. Naar aanleiding van deze processen, sequenties uit de CAR9 monsters bedroeg 591.600 en STAR7 sequenties bedroeg 939.311. Deze sequenties werden geüpload naar MGRAST en gemonteerd via assembler software met de standaardinstellingen. Contigs werden omgezet in 6 reading frames en vermeende open reading frames (pORFs) werden geïdentificeerd via scripts, zoals eerder beschreven 21.
Om onbekende ORF een aantal gelijkenis gebaseerd zoekopdrachten werden uitgevoerd om ORF van bekende functie te verwijderen identificeren. In het kort werden de volgende zoekopdrachten die met hun overeenkomstige zoeken criteria 21:
- Aanzienlijke gelijkenis van ≥ 95% identiteit dan ≥40 basenparen (bp) van MGRAST BLAT in M5NR database.
- Aanzienlijke gelijkenis (e-waarde ≤ 0.001) door TBLASTN tegen alle openbare metagenomes Mijn metagenoom Database Resource.
- Aanzienlijke gelijkenis (e-waarde ≤ 0.001) door BLASTP en TBLASTN tegen NR database.
- Aanzienlijke gelijkenis (e-waarde ≤ 0.001) door RPS-BLAST tegen de Behouden Domain Database.
- Eiwit vertalingen uit een select fractie van elke dataset vergeleken opgelost eiwitstructuren in de Protein Data Bank.
- Berekening van dinucleotide frequenties met Dinucleotide Handtekeningen pakket.
De resulterende pORFs zijn ontworpen voor expressie in E. coli met behulp van publiek beschikbare gen design software. Back-translatie van de aminozuursequenties gebruik gemaakt van een Universal Codon Usage Table ontworpen om tegemoet expressie in E. coli met een minimum verbruik drempel van 2%. Restrictie-enzym herkenningssequenties voor BamHI en HindIII werden uit de sequenties uitgesloten klonering. Een extern bedrijf gesynthetiseerd de gemanipuleerde gen sequenties 30 en vervolgens ORF's werden gekloond in een medium-copy number pBAD promotor vector, pEMB11, via standaard restrictie-enzym klonen. Alle klonen werden getransformeerd in E. coli K-12 stam BW 27784 23.
Multi-fenotype Assay platen (MAP)
Een hoge doorvoer en robuuste software pijplijn werd ingezet voor de analyse van de kaarten, de PMAnalyzer 24. De pijpleiding werd ontwikkeld in een Linux-server-omgeving en voert verschillende stappen, waaronder: parsing van optische dichtheid van bestanden, het formatteren van gegevens in leesbare tekst bestanden, pre-processing van groeicurven voor de kwaliteitsborging (QA), en het uitvoeren van wiskundige modellering technieken om groeicurven analyseren . De primaire modellering scripts ontwikkeld in Python versie 2.7.5 gebruik van de PyLab module maken.
MAP reproduceerbaarheid werd bepaald met behulp van de standaard fout (SE) voor het repliceren van gegevens (Figuur 2A). Raw groeikrommen werden vergeleken met logistische groeicurven te bepalen of het programma PMAnalyzer nauwkeurig geparametriseerd en gemodelleerd kloon groei tijdens proeven (gegevens niet getoond). Voor meer informatie over de juistheid en geldigheid van de kaarten en PMAnalyzer zie Cuevas et al. 24
Na validatie van de werkwijze werd MAPs gegevens geanalyseerd met de verschillende parameters, zoals de maximale groeisnelheid (μ max) en groeipeil (GL), door verwerking pijpleiding. Vergelijkende visualisatie van groeicurves wordt vaak gebruikt voor de groei van data interpretatie; echter, het aantal bochten dat kan worden gevisualiseerd op een tijdstip voor de vergelijking heeft beperkingen. Om tal van groeicurven tegelijkertijd te analyseren, werden heat map afgeleid percelen smeekte om te vergelijken tientallen klonen gegroeid op een enkele substrate tegen gemiddeld reactie dat de omstandigheden '(figuur 2B). De invloed geplaatst door overexpressie van een nieuwe faag eiwit wordt waargenomen door veranderingen in de curve parameters, in het bijzonder: vertraging fase exponentiële fase en de maximale biomassa opbrengst (asymptoot). Als voorbeeld wordt de klim van opstartfase in exponentiële fase gemodelleerd in de groeicurve voor het capside-eiwit (Figuur 2A) gereproduceerd door een snelle verandering in kleurintensiteit van zwart naar wit voor dezelfde kloon in de dynamische grafiek van figuur 2B .
Een algemeen beeld van kloon verdeling over substraten te verkrijgen, werden fenotypische classificaties afgeleid uit de GL gebruikt (figuur 3). Hier zijn de vier fenotypes gescheiden in vier grafieken waarin de hoogte van elke staaf geeft het aantal klonen tonen die fenotype voor een specifiek substraat. Uitschieters in de data worden erkend als klonen vallen in de "winst van functies "of" het verlies van de categorie functie ". Uitschieters kan dan afzonderlijk worden gezocht en onderzocht nauwer experimenteel. Bovendien globale analyse herkent substraat biases in de assay. Substraten zoals fenylalanine, appelzuur en glycine leidde tot een "geen groei" classificatie. Substraten die consequent vallen in de geen groei classificatie, in alle klonen, zijn niet zwaar gewogen in downstream functionele karakterisatie.
Metabolomics
De katabole producten van klonen die faaggenen onbekend werden geïdentificeerd met behulp van metabolomics. In het kort werden klonen gegroeid onder ofwel een continue cultuur of seriële passage in batch cultuur omgeving alvorens te worden verzonden voor GC-TOFMS analyse op een metabolomics kern faciliteit. Voor meer informatie over het monster verwerking, analyse en normalisatie voor GC-TOFMS geïmplementeerd door de gekozen kernfaciliteit zien Fiehn et al. 31 Briefly, 1 ml koude extractie-oplosmiddel wordt aan elk monster toegevoegd, waarna monsters gewerveld en gesonificeerd in een koud bad gedurende 5 minuten. De monsters worden tenslotte gecentrifugeerd en de helft van het monster wordt gedecanteerd en neer gedroogd voor analyse. Extracten zijn gezuiverd en verrijkt met interne retentie index markers voor op de gaschromatograaf worden geladen en dan vervolgens overgebracht naar de massaspectrometer. De gegevens van elk monster geanalyseerd zodanig dat signaal intensiteit voor alle gedetecteerde signalen in het chromatogram worden gerapporteerd. Voor de normalisatie, is de overvloed van pieken voor elk monster opgeteld en de totale piek abundanties zijn gemiddeld tussen alle monsters in de set. Metaboliet abundanties per monster worden gedeeld door het monster piek overvloed en vervolgens vermenigvuldigd met de gemiddelde piek overvloed aan de sample set. De resulterende gegevens worden gebruikt voor de analyse van metabolomics in het onderzoek besproken.
Validatie van metabolomics reproduceerbaarheid was nodig omhet bepalen van de juiste steekproefomvang voor elke kweken methode. Om de precisie gezien binnen monsters en de variatie gezien in het monster te detecteren maten de standaardafwijking van het gemiddelde (S M) voor zowel n = 3 en n = 6 datasets werden beoordeeld (Figuur 5B). Ongeacht de continue kweek (CC) monstergrootte, minder dan 1% van de gegevens hadden S M ≤ 1,5. Mediane S M s waren 221 en 300, en waarden varieerde van respectievelijk 0 tot 7,55 x 10 5 en 3,74 x 10 5 voor n = 3 en n = 6. S M s werden ook berekend voor elke set van het monster repliceert in de methode seriële cultuur (SC). Nogmaals, minder dan 1% van de gegevens had een S M ≤ 1,5, een mediane S M van 137, en een scala 0-3,51 x 10 5. Om de S M distributies tussen elke sample set (CC n = 3 vs vergelijken . CC n = 6, n = 3 CC vs. SC n = 3, en CC n = 6 vs. SC n = 3) een permutatie test werd uitgevoerd. De distribution hetzij continukweek S M waarden dataset was niet significant verschillend van de verdeling van de seriële kweek S M waarden (p-waarde = 0,0). De verdeling van de M waarden voor continue kweek n = 3 data significant verschilde van dat van de S M waarden voor de continue kweek n = 6 data (p-waarde = 1,908804 x 10 -49). Ten slotte, de variatiecoëfficiënt per metaboliet werd vergeleken voor en na implementatie van een kwaliteitsborging (QA) stap (figuur 5C). In totaal werden 210 metabolieten verwijderd na de uitvoering van de QA pipeline (40% van de gegevens). Minder dan 1% van de verwijderde gegevens had metaboliet overvloed aan nul, ~ 2% was interne standaard data, ~ 5% werden gegevens van metabolieten niet eerder waargenomen in E. coli, en de resterende metabolieten (> 30%) een variatiecoëfficiënt groter dan 1.
Net als bij de MAP analyse globale observations voorzien in een eerste inzicht in de diepte van informatie metabolomics aanbiedingen. Een globaal beeld te krijgen, werden klonen hiërarchisch gegroepeerd op basis van hun relatieve abundanties metaboliet voorlichting over kloon-metabolietprofielen, potentiële klonen met verwante functies en kloon-metaboliet uitschieters (figuur 6). Om eiwit functies te benadrukken, zijn metabolieten gescheiden en gegroepeerd op basis van gemeenschappelijke metabole routes. Met deze analyse voorlopige resultaten was het duidelijk dat metabolomics kan scheiden genen uit verschillende klassen (figuur 6, gemarkeerd klonen). Bovendien uitbijter identificatie met metabolomics data werd bepaald door standaard scores (z scores) voor elke kloon-metaboliet pair. Om statistische significantie te waarborgen, werden outliers gedefinieerd als een kloon-metaboliet pair met een Z score waarde van 2, goed voor slechts 5 procent van de data (data niet getoond).
ether.within-page = "always"> 
Figuur 1. Definities van fenotype classificaties. (A) Verhouding tussen de groei-niveau (GL) en de maximale groeisnelheid. Datapunten rood omcirkeld vertegenwoordigen groeicurves tonen weinig tot geen substraatgebruik. (B) Boxplot vertegenwoordiging definiëren groei drempel gebaseerd op de verdeling van de groeicurven met een minimale groei (<0,15 OD / hr). (C) de variantie en de standaardafwijking van GL wordt berekend substraat D-galactose. Korte stippellijnen vertegenwoordigen twee standaarddeviaties van het gemiddelde. Klik hier om een grotere versie van deze figuur te bekijken.
highres.jpg "/>
Figuur 2. MAP validatie via precisie en differentiatie. (A) groeicurven voor geannoteerde structurele (capside) en metabole klonen (Thioredoxine), twee nieuwe metabole klonen (EDT2440, EDT2441), en de gemiddelde respons van klonen gekweekt op sucrose, D- galactose en D-mannose in de kaarten. Blauwe lijnen geven de standaardfout waargenomen bij herhaalde gegevens (n = 3). (B) De groeicurves van 47 verschillende klonen worden voorgesteld als warmte map sucrose, D-galactose en D-mannose. De geannoteerde structurele (groene cirkel) en metabole klonen (oranje cirkel), twee nieuwe metabole klonen (donker en licht blauwe cirkels), en de gemiddelde respons (rode cirkel) worden gemarkeerd. Klik hier om een grotere versie van deze figuur te bekijken.
Pbelasting / 52.854 / 52854fig3highres.jpg "width =" 700 "/>
Figuur 3. Kloon de distributie voor elke fenotype over meerdere ondergronden. Het fenotype-kloon tellen voor 47 klonen over 72 koolstof-specifieke groeiomstandigheden. Tabel geeft direct telt voor elke fenotype. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4. Schema waarin de constructie van de continue kweek apparaat. (A) stap naar een haven α-γ van de continue kweek reactor bouwen, (B) de stappen gebruikt om de stroompoort van de continue kweek reactor en bouw (C ) de stappen om havens δ en ε van de continue cultuur zuigfles te bouwen. = "Https://www.jove.com/files/ftp_upload/52854/52854fig4highres.jpg" target = "_ blank"> Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5. Vergelijking van de phenomic werkwijzen gepresenteerd. (A) workflow voor de bereiding van de Multi-fenotype testplaten (MAP), continu kweken en seriële culturen. (B) Het percentage van de standaard fout van de gemiddelde (S M) geldt zowel voor n = 3 en n = 6 steekproefomvang voor de continue cultuur (CC) en seriële cultuur (SC) bereidingswijzen voor metabolomics. De y-as is op logaritmische schaal. (C) De verdelingen van variatiecoëfficiënten (CV) per metaboliet, voor en na de uitvoering van de QA strategie voor de continue kweek (CC) methode.g5highres.jpg "target =" _ blank "> Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6. Metabolomic profielen van klonen gegroeid in continue cultuur. De mediane metaboliet abundanties voor een set van metabolieten worden uitgezet voor 84 klonen gegroeid in continue culturen. Metabolietprofielen voor geannoteerde structurele (Capsid) en metabole klonen (Thioredoxine), twee nieuwe metabole klonen (EDT2440, EDT2441), en de gemiddelde metabole respons zijn gemarkeerd in het rood. Klik hier om een grotere versie van deze figuur te bekijken.
| Samenstelling | Koolstof | Stikstof | Zwavel | Fosfor |
| Glycerol | - | 0,40% | 0,40% | 0,40% |
| Ammoniumchloride | 9,5 mM | - | 9,5 mM | 9,5 mM |
| Natriumsulfaat | 0,250 mM | 0,250 mM | - | 0,250 mM |
| Magnesiumsulfaat | 1,0 mM | 1,0 mM | - | 1,0 mM |
| Kaliumfosfaat | 1,32 mM | 1,32 mM | 1,32 mM | - |
| Magnesium chloride | - | - | * | - |
| Kaliumchloride | 10 mM | 10 mM | 10 mM | 10 mM |
| Calciumchloride | 0,5 uM | 0,5 uM | 0,5 uM | 0,5 uM |
| Sodium chloride | 5 mM | 5 mM | 5 mM | 5 mM |
| IJzerchloride | 6 uM | 6 uM | 6 uM | 6 uM |
| L- arabinose | 0,10% | 0,10% | 0,10% | 0,10% |
| MOPS pH 7,4 | 1x | 1x | 1x | 1x |
Tabel 1. De verbindingen en concentraties van de verschillende basale media in de MAP. * 1,0 mM magnesiumchloride gesubstitueerd. 1x MOPS = 40 mM MOPS, 4 mM Tricine.
| Koolstofsubstraten | Stikstof substraten | Zwavel substraten | Phosphorus substraten |
| 2 deoxy-D-ribose | 2-deoxy-D-ribose | 1-butaan-sulfonzuur | adenosine-5-monophoshate |
| 4-hydroxy- fenylazijnzuur | aceetamide | acetylcysteïne | beta-glycerofosfaat |
| azijnzuur | adenine | D-cysteïne | creatinephosphate |
| adenosine-5-monofosfaat | adenosine | D-methionine | D-glucose-6-fosfaat |
| adonitol | allantoïne | diethyl-dithiofosfaat | diethyl-dithiofosfaat |
| alfa-D-glucose | beta-fenylethylamine | DL-ethionine | DL-alpha-glycerofosfaat |
| alfa-D-lactose | biureet | glutathion | kaliumfosfaat |
| alfa-D-melebiose | cytidine | isethionzuur | natriumpyrofosfaat |
| citroenzuur | cytosine | L-cysteïnezuur | natrium thiofosfaat |
| D-alanine | D-alanine | L-Cysteïne | |
| D-arabinose | D-asparagine | L-djenkolic zuur |
| D-arabitol | D-aspartaat | L-methionine |
| D-asparagine | D-cysteïne | magnesiumsulfaat |
| D-aspartaat | D-glucosamine | methaansulfonzuur |
| D-cellubiose | D-glutaminezuur | N-acetyl-DL-methionine |
| D-cysteïne | DL-a-amino-N-boterzuur | N-acetyl-L-cysteine |
| D-fructose | D-methionine | kalium-tetra-thionate |
| D-galactose | D-serine | natriumthiosulfaat |
| D-glucosamine | D-valine | sulfanic zuur |
| D-glucose | gamma-amino-N-boterzuur | taurine |
| D-glucose-6-fosfaat | glycine | taurocholinezuur |
| D-glutamaat | guanidine | thioureum |
| D-mannose | histamine | |
| D-raffinose | inosine |
| D-ribose | L-alanine |
| D-salicine | L-arginine |
| D-serine | L-asparagine |
| D-trehalose | L-citrulline |
| D-xylose | L-cysteine |
| dulcitol | L-glutaminezuur |
| glycerol | L-glutamine |
| glycine | L-glutathion |
| i-erythritol | L-histidine |
| inosine | L-isoleucine |
| L-alanine | L-leucine |
| L-arabinose | L-lysine |
| L-arabitol | L-methionine |
| L-asparagine | L-ornithine |
| L-aspartaat | L-fenyl- alanine |
| L-cysteïnezuur | L-proline |
| L-cysteine | L-pyro-glutaminezuur |
| L-fucose | L-serine; L-threonine |
| L-glutaminezuur | L-tryptofaan |
| L-glutamine | L-valine |
| L-isoleucine | N-acetyl-D-glucosamine |
| L-leucine | putrescine |
| L-lysine | thioureum |
| L-methionine | thymidine |
| L-fenylalanine | thymine |
| L-pyro-glutaminezuur | tyramine |
| L-rhamnose | tyrosine |
| L-serine | uridine |
| L-sorbose | |
| L-threonine |
| L-tryptofaan |
| L-valine |
| L-xylose |
| lactaat |
| lactulose |
| malaat |
| myo-inositol |
| oxaalzuur |
| kaliumsorbaat |
| propionzuur |
| putrescine |
| kinazuur |
| natriumpyruvaat |
| natriumsuccinaat |
| sucrose |
| thymidine |
| xylitol |
Tabel 2 Lijst van materialen in de MAP experimenten.