[好文分享：www.ii77.com]

摘要

生命早期对免疫系统的发育很主要，对成年后的健康也有影响。固然已经确定出生后很快就获得了肠道细菌微生物组，但在生命的最初几年，我们对由噬菌体、真核RNA病毒和真核DNA病毒构成的病毒组（virome）却知之甚少。在这篇文章里，我们描述了健康双胞胎婴儿纵向队列中肠道病毒组和细菌微生物组的特征。比拟无关婴儿，同卵双胞胎之间的病毒组和细菌微生物组更相似。从出生到2岁，真核病毒组和细菌微生物组扩大，但却陪伴着噬菌体病毒组构成的收缩（contraction）和改变（shift）。 噬菌体-细菌的关系始于高捕食者-低猎物的动态（high predator–low prey dynamic），与Lotka-Volterra prey模型一致。是以，与在成年人中视察到的不乱的微生物组比拟，婴儿微生物组是高度动态的，而且与生命早期的细菌、病毒和噬菌体跟着岁数的转变相关系。 [原文来自：www.ii77.com]

配景介绍

肠道微生物组包罗细菌、真核病毒、细菌病毒(噬菌体)、真菌和古菌。已经确定，这些微生物的一部门会与免疫系统互相感化并影响宿主的健康。肠道细菌微生物组的改变与多种人类疾病有关，包罗肝硬化、糖尿病和炎症性肠病。大多数针对微生物组的治疗策略，如益生菌、益生元和粪菌移植，都旨在调节细菌微生物群落。细菌微生物组在出生后很快就竖立起来，在接下来的几年里，它的构成向一个尺度（stereotypical）的“类成人（adult-like）”的细菌群落构造改变。
这一过程会受到多种互相感化身分的影响，包罗营养、出生模式、抗生素使用和地舆情况。对双胞胎的研究表明，同卵双胞胎比无关个别共享更相似的细菌微生物组。

我们对病毒组(virome)知之甚少，这是一个由真核RNA病毒、真核DNA病毒和噬菌体构成的多样化群落。新的证据表明病毒组对人类健康有影响。 Anellovirus(真核DNA病毒)的载量（burden）与宿主免疫按捺的水平和器官移植的究竟直接相关，也是儿科发烧性疾病和AIDS的指标。致病性猴免疫缺陷病毒（Pathogenic simian immunodeficiency virus）与肠道病毒的扩张（expansion）有关，个中包罗很多真核RNA病毒。此外，慢性病毒传染可增加寄主对病原性挑战的抗击力，表明病毒或者对寄主发生有益的影响。肠道微生物还含有多种噬菌体，在健康成年人体内，这些噬菌体首要由Caudovirales 目和Microviridae 科构成。这些噬菌体平日在一段时间内连结不乱的群落。肠道噬菌体群落构成的转变与克罗恩病（Crohn’s disease）和溃疡性结肠炎（ulcerative colitis）有关。然而，与情况生态系统的噬菌体-细菌互相感化的种群动态转变遵循Lotka-Volterra “捕食者-猎物（predator-prey）”模型分歧，噬菌体和细菌之间的捕食者-猎物（predator-prey）关系尚未在人类肠道微生物组中视察到。健康婴儿肠道病毒的宏基因组学研究仅限于一项对单个婴儿的研究，个中使用适度深度的桑格测序在单个时间点剖析了DNA病毒。有针对性的PCR和RT-PCR研究已经判定出了一些真核病毒，包罗经常在健康婴儿粪便中发现的picornaviruses和anelloviruses。尽管已经报道过一些对腹泻（diarrhea）和急性弛缓性麻木（acute flaccid paralysis）等疾病患儿的肠道病毒的宏基因组研究，但迄今为止，还没有对一组健康婴儿的病毒组进行纵向剖析的研究。

鉴于细菌微生物组是在婴儿的生命早期竖立的，并或者影响历久健康，我们研究了陪伴人体发育的真核病毒和噬菌体的转变。为了说明病毒个别间和个别内变异的水平，我们对一对健康的同卵双生双胞胎和三对健康的异卵双胞胎的粪便样品进行测序。在这项研究中，我们将“健康婴儿”界说为没有显着潜在遗传或慢性疾病的婴儿，尽管这些婴儿也会发生急性病。为了确定病毒组的构成及其陪伴岁数增进的演变，我们对照了从出生到两岁的六个时间点收集的粪便样本中的肠道病毒组。此外，我们对统一粪便中的细菌16S rRNA基因进行了测序，以生成一个关于人体肠道病毒组和细菌微生物组发育关系的综合考查。我们的究竟供应重建了一个追随时间转变的深入的婴儿肠道病毒组的动态转变过程，并揭示了在健康婴儿发育过程中天然发生的捕食者-猎物，即噬菌体-细菌（predator-prey bacteriophage-bacteria）的动态转变纪律。

研究究竟

婴儿早期发育的病毒组

Virome of infants during early development

我们对栖身在美国密苏里州圣路易斯市多半会区的八名健康婴儿(四对双胞胎)的粪便样本进行了宏基因组测序(图1a)。 本研究平分析的样本从出生第1-4天(界说为第0个月)以及3、6、12、18和24个月时进行纵向收集，以确定婴儿早期发育时代的肠道微生物组。为了周全检测DNA和RNA病毒，从粪便标本中提取总核酸，并采用两种互补的扩增方式：多重置换扩增法（multiple displacement amplification，MDA）和序列自力的DNA和RNA扩增法（sequence-independent DNA and RNA amplification，SIA）。MDA法因为phi29聚合酶有高效持续的合成能力而被普遍使用。然则MDA法不克检测到RNA病毒，也会导致对小型环状DNA分子的倾向性扩增。为了对这种方式进行增补，我们采用了SIA法，检测到DNA的同时也能检测到RNA，然则它的DNA探测能力一样不如MDA法敏感(图1b)。我们归并了文库，在Illumina MiSeq..长进行测序。每个MDA样品平均获得了549,301 ± 207,521 (mean ± s.d.) 条reads，每个SIA样品平均获得了551,592 ± 229,210(mean ± s.d.) 条reads。接着，去除了这些reads的接头，进行质量过滤，分类学分派（assign），每个样品随机重抽样（randomly subsampled）获得200,000条reads。为了界说全局（global）的病毒组构成，我们生成了一个由MDA方式和SIA方式判定的病毒科水平的文件，而且生成了一个存在-缺失（presence-absence）的热图(图1c)。这和以前基于PCR的研究究竟一致，我们在婴儿粪便的样品中判定到了真核RNA病毒，例如caliciviruses, astroviruses和picornaviruses。此外，比起真核RNA病毒和真核DNA病毒，在科水平我们检测到更多的噬菌体。此外，只判定出一个古菌病毒科(Lipothrixviridae) 。

为了对照个别之间的病毒组多样性，我们在病毒属的水平用unweighted Bray-Curtis距离测定了β-多样性。病毒组的PCA剖析表明，双胞胎的相关性（co-twin relatedness）和岁数对病毒多样性的转变有进献。凝聚的分层聚类（Agglomerative hierarchical clustering）也支撑双胞胎的病毒组构成比无关个别更相似(图2)。

图1. 婴儿肠道病毒组的实验设计和宏基因组剖析

Figure 1 Study design and metagenomic analysis of the infant gut virome.

(a) 8个健康婴儿(4队双胞胎)微生物组研究的测序方式

(b) 热图展示分歧测序方式下分派至各病毒科的reads数量。可见受方式影响很大。只蒜了7个样本的方式对照，如A2-24为A2号个别24个月取样。

(c) 有无热图展示8个样本6个时间点上病毒组的状况。个别为A1-D2，时间点为0, 3, 6, 12, 18和24个月。

图2. 病毒组Beta多样性剖析

Figure 2 Analysis of virome beta diversity.

凝聚条理聚类剖析病毒组(基于真核病毒和噬菌体属水平)分歧时间点的Bray-Curtis距离。

在出生后获得真核病毒组

The eukaryotic virome is acquired after birth

接下来，我们存眷真核病毒组，真核病毒群落在最早期的样品中品貌很低，随后升高(图3a)，表明真核病毒组经由情况露出而竖立。在判定出的真核RNA病毒中，最常见的病毒属是肠病毒(enterovirus)、副肠孤病毒(parechovirus)、烟草花叶病毒(tombamovirus) 和札幌病毒(sapovirus)（图3b）。真核病毒相对稀少，这故障了对常见生态参数(包罗多样性和稀少性)的正确评估。为了确定同卵双胞胎是否在沟通的时间点携带沟通的病毒株（virus strains），我们剖析了副肠孤病毒的序列副肠孤病毒是我们的研究中检测到的最遍及的RNA病毒之一，经由组装读段(reads)并将所得的重叠群(contigs)映射（map）到副肠孤病毒基因组上。系统发育剖析表明，双胞胎中的婴儿共享几乎沟通的副肠孤病毒株（strains），然则分歧的双胞胎具有分歧的株（strains）。与这一发现相一致的是，在同卵双胞胎中也检测到了肠病毒，即另一种高度风行的真核RNA病毒。固然这些发现表明双胞胎中的婴儿经常传染沟通的病毒，但我们也视察到只有双胞胎中的一个婴儿检测到，而另一个婴儿没有检测到的病毒。例如，在婴儿 A1 3个月大的样品中检测到副肠孤病毒的序列，但在其同卵双胞胎A2中却未检测到（图3c）。为了确定视察到的纷歧致是因为测序活络度的局限照样病毒载量（viral load）的差别造成的，我们用定量逆转录-聚合酶链回响(quantitative RT-PCR，qRT-PCR)筛选了所有样品，并测量了人体parechoviruses的拷贝数。所有五个测序阳性样品都经由了qRT-PCR检测，检测呈阳性。此外的三个测序阴性样品在qRT-PCR检测中呈阳性，不外病毒载量（viral load）非常低。最后，qRT-PCR和测序都证实“纷歧致（discordant）”的婴儿A2携带与其同卵双胞胎A1沟通的人类副肠孤病毒星散物。

是以，qRT-PCR究竟自力验证了经由深度测序检测到的病毒的存在，并进一步证实了同卵双胞胎共享相似的病毒组。

指环病毒(Anelloviruses)是检测到的最遍及的真核DNA病毒科（图3d）。 几乎所有的指环病毒都是未知的，而且与以前描述的指环病毒差别很大（图3e）。因为以前的报道表明，指环病毒载量与宿主免疫状况的转变有关，所以我们假设婴儿免疫的转变(例如，母体抗体的削弱)或者反映在指环病毒风行率（prevalence）或品貌的转变上。因为检测到大量新指环病毒的序列，我们将一组非冗余的指环病毒的重叠群作为“参考基因组（要求小于95%的核酸序列相似度）”。然后，我们将每个样本的读段比对到指环病毒参考重叠群上，以确定每个重叠群的存在遍及性。这种方式与用于检测三种特定指环病毒的重叠群的PCR检测究竟的一致率为98.6%。指环病毒在3个月大之前很少被检测到，但很快显著增加，在6-12个月大时达到巅峰（图3g）。

值得注重的是，一名婴儿在12个月大时至少携带了47种指环病毒。此外，同卵双胞胎比无关婴儿共享更高比例的指环病毒（图3h）。 此外，在一些情形下，能够从统一婴儿的粪便中检测到沟通的指环病毒，这些粪便被收集的时间相隔为12个月（图3f），这表明持续存在或许是频频传染的不乱起原。

图3. 随岁数原核RNA和DNA病毒的替代以及双胞胎间共享病毒组

Figure 3 Alterations in the eukaryotic RNA and DNA viruses with age and evidence of shared viromes between co-twins

(a) 在分歧的岁数，视察到原核DNA和RNA的分类群的品貌。使用的方式是线性回来，R²和95%置信度如图所示

(b) 检测出的原核RNA病毒属的样本数量

(c) 副肠孤病毒(parechovirus)序列基于最大似然方式的系统发育剖析。Bootstrap值如图所示

(d) 包含了展示出的原核DNA病毒属的样本数量

(e) 61 个指环病毒(Anelloviruses)重叠群和12个参考株序列中的ORF1氨基酸序列比对究竟，并采用最大似然方式建树，图中展示了属名

(f) 有无热图展示了指环病毒重叠群，配色与E的系统发育树的颜色一致

(g) 在分歧的时间点指环病毒的品貌。用Wilcoxon磨练进行统计显著性剖析。

噬菌体群落的构成跟着岁数收缩和改变

Bacteriophage community contracts and shifts in composition with age

尽管在所有样品中都检测到了DNA噬菌体，然则采用SIA法发生的数据并没有发生任何RNA噬菌体。在没有RNA噬菌体的情形下，我们选择将随后的噬菌体剖析集中在采用MDA法发生的数据上，如许我们的发现能够与其他采用MDA法的噬菌体病毒研究比拟较。与真核病毒比拟，噬菌体雄厚度在0个月时最大，并跟着岁数的增进而降低（图4a）。雄厚度稀少曲线表明，细菌种类的储蓄速度跟着岁数的增进而降低，这表明雄厚度的降低不太或者是由抽样误差造成的（图4b）。同样，噬菌体多样性随岁数增进而降低（图4c）。同卵双胞胎之间噬菌体病毒组的人际变异（interpersonal variation）比无关婴儿之间的要低（图4d），即病毒组更相似。与其他研究一致，最雄厚的噬菌体来自有尾噬菌体目(Caudovirales)(包罗长尾噬菌体科Siphoviridae、丝状噬菌体科Inoviridae、肌病毒科Myoviridae 和短尾噬菌体科Podoviridae) 和微小噬菌体科(Microviridae) 。然而，在24个月大时，群落构成发生了显着的转变，微小噬菌体科 科的相对品貌增加（图4e）。

这种微小噬菌体科品貌的转变在采用SIA方式发生的数据中也发现了，表明它不是由方式误差造成的。还视察到微小噬菌体科物种雄厚度的增加，表明这种扩张（expansion）不是由特定物种驱动的。有尾噬菌体目(Caudovirales)的相对品貌与微小噬菌体科的相对品貌成反比（图4f）。*crAssphage仅存在于一个样本中，表明crAssphage不是在生命早期获得的。是以，婴儿早期发育的特点是噬菌体群落雄厚性和多样性的收缩，陪伴着向以微小噬菌体科为优势组分的改变。

图4. 岁数增进噬菌体品貌下降且构成发生改变

Figure 4 Decrease in bacteriophage richness and diversity with age coincides with a shift in bacteriophage composition

(a) 在分歧的岁数月份下噬菌体物种的品貌。使用了线性回来，R²和95%置信度如图所示

(b) 稀少曲线，展示了被识别噬菌体物种的雄厚度

(c) 在分歧时间点的噬菌体的α-多样性，使用了线性回来，标注了R²和95%的置信区间

(d) 两对(白色)噬菌体病毒体在属水平上的Bray-Curtis距离，按双胞胎内(Winthin twins，白色)与双胞胎间(Between twins，灰色)分组比拟较。统计显著性经由学生t磨练进行评估； P = 0.01–0.05，* P < 0.01。

(e) 噬菌体科水平的相对品貌。按科分类，未定名的按目分类。

(f) 该图显露了微小噬菌体科家眷品貌与有尾噬菌体目品貌的关系(n = 48个采样时间点)。线性回来和95%置信区间，并显露Spearman相关系数(r)和统计显著性

婴儿的细菌微生物组的转变

Bacterial microbiome changes in infants

跟着细菌种群雄厚度和多样性的络续增加，婴儿早期健康肠道细菌微生物组的生态特征在2-3岁时发育成熟，成为更不乱的“类成人（adult-like）”种群。 鉴于我们在噬菌体病毒组中视察到的伟大转变，我们试图认识这些转变是否与细菌微生物组的转变相关。是以，我们进行了细菌16S rRNA基因测序，平均每个样本发生67,569条序列。经由质量过滤后的序列以97%的统一性阈值被聚类到一个把持分类单元（OTU）中。未加权的UniFrac距离矩阵（unweighted UniFrac distance matrices）的PCoA剖析表明 细菌群落的转变与岁数有关。细菌雄厚度（图5a）和多样性（图5b）也跟着岁数增进而增加的现象与其他研究一致。总的来说，我们判定出越来越雄厚的梭菌属Clostridia (厚壁菌门Firmicutes)（图5c）。在此之前，芽孢杆菌纲Bacilli (Firmicutes)的OTUs在0个月占优势，Gammaproteobacteria (Proteobacteria)和Actinobacteria  (Actinobacteria)的品貌在3个月和6个月增加，Bacteroidia (Bacteroidetes)品貌在12个月、18个月和24个月增加。同卵双胞胎之间的个别间差别小于无关婴儿之间的差别。（图5d） 是以，本研究中婴儿的细菌微生物组与之前视察到的预期转变轨迹一致。

图5. 细菌群落随岁数扩张

Figure 5 Bacterial community expansion with age

(a) 在分歧的岁数下细菌物种的雄厚度。使用了线性回来，R²和95%置信度如图所示

(b) 分歧时间点的细菌的α-多样性，使用了线性回来，标注了R²和95%的置信区间

(c) 基于16S rRNA基因测序，细菌科水平的相对品貌

(d) 双胞胎之间和无关个别之间细菌群落的Unifrac距离，统计显著性经由学生的t磨练进行评估，** P < 0.01

噬菌体与细菌雷同捕食者-猎物的关系

Predator-prey–like bacteriophage and bacteria relationships

尽管海洋中噬菌体-细菌的关系显露出捕食者-猎物（predator-prey）的动态，但噬菌体和细菌种群在成年人肠道中相对不乱。然而，还没有研究过在婴儿中是否存在这种捕食者-猎物（predator-prey）的关系。在我们的研究队列中，噬菌体多样性与细菌多样性成反比（图6a）。经由搜检随时间的趋势，我们发现微生物群落在0个月时从高噬菌体-低细菌的多样性群落改变为24个月时的低噬菌体-高细菌的多样性群落。

噬菌体和细菌雄厚度由岁数决意的体式成反比也与这一发现一致（图6b）。进一步反映这些生态趋势，噬菌体和细菌的特定属之间的相关性，使用夹杂线性模型较量，首要是负相关。是以，婴儿病毒组和细菌微生物组在生命早期以动态轨迹进化。

图6. 噬菌体和细菌之间的负相关关系

Figure 6 Inverse ralationships between bacteriophages and bacteria

(a) 噬菌体多样性和细菌多样性之间的相关性(n = 48个采样时间点)。线条透露线性回来，并显露Spearman相关系数。颜色透露岁数从0到24个月不等。

(b) 噬菌体品貌和细菌品貌之间的相关图(n = 48个采样时间点)。线条透露线性回来，Spearman相关系数(r)透露。颜色透露岁数从0到24个月不等

商议

肠道微生物之间的互相感化影响宿主的生理、发育和免疫。婴儿肠道细菌微生物组的组装（assembly）或者对成年时的微生物组有决意性的感化，这会对宿主表型的某些方面发生历久影响，包罗肥胖、炎症肠病和食物过敏。然而，我们对病毒在婴儿体内的发育体式、对细菌微生物组的影响或在人类健康中的感化都知之甚少。在这个研究中，我们纵向描述了8个婴儿(4对双胞胎)的完整病毒组和细菌微生物组，并揭示了健康婴儿病毒组发育过程中的微生物里程碑。双胞胎研究设计使我们可以确定同卵双胞胎之间的婴儿微生物组(病毒组和细菌群落)比无关婴儿之间更相似。
这一究竟与一项成年同卵双胞胎的研究形成了对比，在这项研究中，DNA病毒对每个个别都是奇特的，双胞胎之间的相似水平并不比无关个别之间的相似水平高。一种或者的注释是，尽管“双胞胎（twin-ness）”这一特点在婴儿期很主要，但这或者反映了如许一个事实，即情况露出是病毒成分的首要驱动身分，因为婴儿双胞胎平日仍然共享一个配合的情况。
当在统一时间点取样时，在同卵双胞胎中检测到非常一致（near-identical）的真核病毒株就证实了这一点。

指环病毒(Anelloviruses)被提议作为功能性免疫能力（functional immunocompetence）的生物标记物，因为血清中指环病毒载量的转变与移植受体的免疫按捺水平有关。此外，还经常经由PCR在婴儿血清中检测到指环病毒。我们视察到在6个月和12个月大的时候肠道中指环病毒雄厚度的增加；这些指环病毒中的大多数与已知的指环病毒有很大分歧，强调了使用无偏的深度测序方式来系统地描述病毒组的价格，因为它们或者难以用PCR检测到。我们推想这种扩张（expansion）或者是免疫状况降低的究竟，因为跟着母体抗体的削弱，它与人类IgG的最低点相一致。

到今朝为止，没有证据表明“杀死赢家（kill-the-winner）”的动态在人类肠道微生物组中发生，因为细菌(猎物)种群的峰值并没有先于噬菌体(捕食者)的增加，随后细菌(猎物)种群削减。固然这是经典Lotka-Volterra“捕食者-猎物”模型中最遍及承认的一个方面，但该模型也描述了有限的猎物多样性掌握捕食者品貌(即捕食者峰先于猎物，也称为反向捕食者-猎物轮回)的互相关系。我们的研究表明，早期婴儿发育中噬菌体-细菌的互相感化始于Lotka-Volterra模型的后一种动态。我们假设噬菌体的多样性在出生时很高，然则噬菌体群体是弗成持续的，因为细菌假寓密度低。这导致噬菌体病毒组的收缩，从而减轻对细菌群落的攫取性压力，使其可以在肠道中形成而且定植。反过来，这鞭策了新竖立的细菌群落中选择的噬菌体构成的改变(包罗微小噬菌体科Microviridae 品貌的增加)。在一项为期2.5年的对单个成体的纵向研究中，微小噬菌体科是首要的噬菌体分类单元，这提出了一个问题，即婴儿期的动态微生物组状况若何以及何时改变为不乱的群落，固然这已经在成年人的研究中报道过。此外，尽管我们的研究无法解决早期噬菌体多样性的起原问题，但我们的数据提高了垂直或产前流传的或者性(首次采样的中央天数为生命的第2.6天)。尽管如斯，我们的数据将作为健康婴儿微生物群削发展的参考。尽管出生时与健康发育相关的噬菌体收缩或者是一种遍及的表型，但特定的肇端噬菌体构成或者因外部身分(如地舆和饮食)而分歧。

众所周知，饮食、抗生素和出生模式会影响细菌微生物组。尽管我们的队列包罗在接合性（zygosity）、母乳喂养状况和临蓐体式方面分歧的个别，但研究队列规模小清扫了这些身分在微生物组中施展感化的剖析。无论若何，我们今朝的研究供应了肠道病毒和细菌在早期发育过程中动态互相感化的具体信息。

方式

取样

Sample

该研究获得了圣路易斯华盛顿大学医学院人类研究珍爱办公室的核准。我们获得了双胞胎婴儿母亲的赞成，每月收集孩子的粪便标本，持续到三岁。粪便样品装在含有冷冻包的绝缘信封中快递给实验室，而且保留在-80℃直至剖析。收集的数据 包罗临蓐体式，婴儿给药的情形，喂养的食物和疾病发生（发烧/吐逆/腹泻）的情形，按期采访怙恃或检测双胞胎大夫的医疗记录。本研究中的4对双胞胎（8名婴儿）被选为健康婴儿的代表，这些婴儿在接合性（zygosity），母乳喂养，性别（6名男性和2名女性）和出生模式方面各不沟通。在这项研究中，我们将这些婴儿界说为“健康”，因为他们没有显着的潜在遗传或慢性疾病。没有其他清扫尺度。婴儿也有急性疾病发生。用于在研究中标记标本的定名方式是，短线前透露婴儿双胞胎，短线后是岁数（例如，A2-24是指来自第A对双胞胎的第二个小孩儿，在24个月龄时收集的样品）。为了手稿的清楚起见，生命岁数界说为0个月（平均2.6天，s.d.±1.1天），3个月（平均98.0天，s.d.± 2.7天），6个月，18个月（平均545.0天，s.d.± 17.4天）和24个月（平均718.5天，s.d.± 11.6天）。

病毒组测序

Virome sequencing

将粪便样品（约200mg）以1：6的比例稀释在磷酸盐缓冲盐水（phosphate-buffered saline,PBS）中，并经由0.45μm孔径的膜过滤。凭据制造商的介绍，在COBAS Ampliprep仪器（Roche）上从滤液中提取总核酸。使用随机数发生器随机化样品（以界说每次测序样品的分组），然后进行以下扩增方式。对总核酸进行序列无关的DNA和RNA扩增（SIA），引物由随机15-mer（15 Ns）上游的碱基均衡的16个核苷酸（nt）特异性序列构成，用于随机引物，如前所述，用于Nextera DNA试剂盒构建文库（Illumina）。

对于多重置换扩增（MDA），凭据制造商的解说用phi29聚合酶（GenomiPhi V2试剂盒，GE Healthcare）扩增总核酸，并用于Nextera DNA试剂盒构建文库（Illumina）。使用Agencourt Ampure XP珠（Beckman-Coulter）纯化和筛选文库，然后使用2100生物剖析仪（Agilent Technologies）定量。归并多重SIA文库并与多重MDA文库分隔测序。一个SIA样本（C2-0）库构建失败，未进行测序。此外，评估文库构建后或者发生的标本交叉污染水平，我们在每次测序时到场了由来自线虫Orsay病毒RNA1区段建构的奇特索引的cDNA文库。

是以，24个文库和1个Orsay病毒对照文库以等摩尔每次归并测序，并在华盛顿大学基因组科学与系统生物学中心的Illumina MiSeq..上测序（总共4次MiSeq测序运行，2×250双端读取，MiSeq v2试剂盒）。

病毒组序列处理

Virome sequence processing

在处理病毒组序列直至分类学分派（assignment）时代，研究人员对分组信息（即岁数，双胞胎，时间点）不知情。对测序读数进行样本拆分（demultiplexed），并去除了接头序列。使用ea-utils包中的fastq-join归并双端序列。以Q30 Phred质量得分的阈值下对序列进行质控，修剪和丢弃低质量核苷酸。

经由查询定制的病毒数据库来识别候选病毒读段(reads)。定制病毒数据库由具有来自NCBI NT和NR数据库（2013年11月7日下载）的“病毒”超等域分类（superkingdom classification）的所有序列构成。CD-HIT用于最小化序列冗余，发生449,469个序列的定制病毒NT数据库和621,095个序列的病毒NR数据库。使用BLASTn（e-value阈值1E-10）依次针对定制的病毒数据库比对测序这场仗，然后是BLASTx（e-value阈值1E-3）。经由使用MegaBLAST（e-value阈值1E-10），BLASTn（e-value阈值1E-10）比对到NCBI NR数据库查询候选病毒reads来过滤假阳性病毒序列，使用BLASTx（e-value值1E-3）比对到NCBI NT数据库以去除具有对应于非病毒序列（例如，人，真菌等）的最佳BLAST比对的序列。测序读段的分类学分派（assignment）由最佳BLAST究竟的分类ID确定。使用内部（in-house）perl和python剧本从BLAST的输出文件中解析科水平宁属水平的分类学分派（assignment）。使用MEGAN（version 5.8.6）中的最低配合祖先算法（lowest-common-ancestor algorithm）确定噬菌体物种分类学分派，参数如下：Min Support: 1, Min Score: 40.0, Max Expected: 0.01, Top Percent: 10.0, Min-Complexity filter: 0.44。因为读段中存在低复杂性/反复区域，剖析中省略了以下假阳性病毒家眷分类学分派：Herpesviridae (3 条reads), Mimiviridae (1 条 read) and Phycodnaviridae (2条reads)。为了评估标本交叉污染，我们评估了每个拆分（demultiplexed）数据集中是否存在映射（map）到对照Orsay病毒的读段。从对照文库自身获得的Orsay病毒测序读数的平均数为463,315条读段（SD ± 106337条读段）。48个婴儿标本中没有一个发生任何Orsay病毒测序读段。

病毒组剖析

Virome analysis

凭据测序深度，SIA法和MDA法生成的测序读段稀释（rarefied）到每个样品200000条读段（5次反复（iteration））。将针对给定病毒分类群检测的reads数量绘制为热图。所有反复（iteration）的剖析都获得了一致的究竟。是以，这解说我们获得的究竟来自于有代表性的反复（iteration）。为了界说全局（global）病毒组的构成，我们归并了经由稀少（rarefied）的SIA和MDA数据，并将其绘制为未加权（存在/不存在）（unweighted (presence/absence) ）的热图。该归并数据用于剖析全局（global）的病毒组（真核病毒和噬菌体）。在真核病毒体群落的自力剖析中，真核病毒的风行情形（prevalence）不足以正确评估最常见的生态测量指标（measurement）。为了剖析噬菌体病毒组的群落，我们使用了MDA数据，因为它更能代表DNA病毒，使我们的研究究竟能够与之前采用MDA法的其他噬菌体组的研究比拟较。使用vagen R包进行生态剖析，包罗雄厚度和多样性测量（Shannon指数，Bray-Curtis相异性），凝聚条理聚类（agglomerative hierarchical clustering）和稀少曲线剖析。使用500个置换/重排（permutations）绘制稀少曲线。使用Emperor 进行主坐标剖析（Principal-coordinate analyses，PCoA）。

为了搜检个别间病毒的株系（isolates），使用bowtie2和 PhyML将测序读段映射到参考基因组（人体parechovirus 1（FM178558），人体肠道enterovirus B（NC_001472）和 crAssphage（NC_024711.1））。用jModelTest2（version 2.16）和ProtTest（version 2.4）评估进化模型，在PHYML（version 3.00）构建最大似然的系统发育树。
经由1000个非参数bootstraps评估对ML树的支撑，而且至少进行两次剖析。

人体副肠孤病毒剖析

Human parechovirus analysis

之前描述过一个靶向基因组中5’UTR区域中的保守序列用于筛选所有样品的副肠孤病毒的泛副肠孤病毒的TaqMan RT-PCR检测。使用的引物如下：AN345F（5’-GTAACASWWGCCTCTGGGSCCAAAAG-3’）和AN344R（5’-GGCCCCWGRTCAGATCCAYAGT-3’），探针AN257（5’/ 6-FAM / CCTRYGGGTACCTYCWGGGCATCCTTC / TAMRA / -3’）。使用TaqMan Fast Virus 1-Step Master Mix（Applied Biosystems）进行qRT-PCR。20μL回响液中包罗5μL提取的样品，10 pmol每种引物和5 pmol探针。使用以下轮回前提：50℃ 5 分钟，95℃ 20秒；40 个轮回进行95℃ 3秒和58℃ 30秒。凭据制造商的方式在体外发生测定的尺度曲线，使用MEGAscript（Ambion）从一个包含方针区域的质粒中发生体外转录的RNA，生成尺度曲线。使用5×10⁶至5个拷贝的体外转录RNA的系列稀释液发生尺度曲线，并界说5个拷贝的检测下限。样品在96孔板形式顶用5个仅含水的阴性对照进行测试。

所有5个阴性对照均为副肠孤病毒阴性。我们试图从副肠孤病毒阳性样品中获得直系同源3D区域以进行系统发育对照。

使用oligo(dT)₂₀，用ThermoScript逆转录酶（Life Technologies）进行3’RACE回响，随后用对人体副肠孤病毒特异的引物（5’-CCAGGTTAACAATGAACTATGGCAG-3’）进行PCR扩增。固然在样品D1-12和D2-12经由了Taqman试验（低病毒拷贝），检测到人类副肠孤病毒，但我们无法经由3’RACE和RT-PCR从这些标本中扩增人体副肠孤病毒。

指环病毒剖析

Anellovirus analysis

很多指环病毒序列与已知的指环病毒只有非常有限的相似性（identity），这表明它们是高度分歧的。是以，我们首先整顿（curate）指环病毒的基因组序列。使用Newbler（version 2.8）从所有质量过滤的读段中从新（de novo）组装重叠群，并使用BLASTx对上述定制的病毒数据库进行查询，以判定指环病毒重叠群。将大于500 nt的指环病毒的重叠群与参考指环病毒基因组进行比对，参考基因组如下：alphatorquevirus TTV1 (AB008394), TTV-P1C1 (AF298585), TTV SIA109 (FJ426280), TTV8 genotype 22 (AB054647); betatorquevirus TTmV1 TLMV-CBD279 (AB026931), TTV-like TLMV-CLC062 (AB038625), TTV-like TTMV_LY2 (JX134045), TTmV5 TGP96 (AB041962), TTV-like LIL-y1 (EF538880), TTV-like LIL-y2 (EF538881), TTmV3 (NC_014088); gammatorquevirus TTmV1 MD1-073 (AB290918), TTmV MDJHem8-2 (AB303557), TTmV MDJN1 (AB303558)。经由采用一致序列（consensus），将具有大于95％核苷酸统一性的重叠群组合。

此外，为了确定指环病毒的系统发育关系，仅将编码ORF1的contig用于进一步剖析。从ORF1氨基酸比对构建ML系统发育树，其包罗上述anelloviruses参考基因组。经由1000个非参数自举来评估对ML树（LG + I + G + F）的支撑。剖析至少进行了两次。

为了在生物信息学方面确定指环病毒的风行率，使用bowtie 2和SAMtools，将测序读段比对到整顿好的指环病毒基因组重叠群。我们评估了采用较量机剖析和采用PCR检测风行率的一致性，选用了三种指环病毒：一种alphatorquevirus anellovirus Contig2355（正向引物 5’-GTAGCCAGAATAAGAACTATGCCC-3’，反向引物 5’-TACTGTCTAAAACCTGGAAGTTGC-3’）; betatorquevirus anellovirus Contig2737（正向引物 5’-TCCAAGAGACTTTAAACCAGGCC-3’，反向引物 5’-GGAACTCCTGGATTGTCCCATC-3’）; gammatorquevirus anellovirus Contig2393（正向引物 5’-CTGATGTAGATGATGGACATGGC-3’，反向引物 5’-CATGAGCTTTGTTGCAGAAAGTC-3’）。选择这三种指环病毒重叠群的尺度如下：（1）选择每个属的代表（α，β和γ），和（2）基于序列数据，重叠群存在于一个婴儿多个时间点的数据中或在多个婴儿的数据中都检测到了。用Taq DNA聚合酶（Life Technologies）在以下轮回前提下进行PCR：95℃ 5分钟，50个轮回95℃ 30秒，58℃ 30秒，72℃ 23秒，然后72°C，10分钟。使用2％琼脂糖凝胶电泳检测产品。

细菌16S rRNA基因测序

Bacterial 16S rRNA gene sequencing

经由珠子研磨提取粪便样品中的核酸。使用特异V4区域的Golay-barcoded引物（F515/R806）。归并等摩尔的文库并使用Illumina MiSeq测序仪（2 × 250 paired-end reads, MiSeq v2 reagent kit）在华盛顿大学的基因组科学和系统生物学中心进行测序。4个样品（A1-0，A2-0，D1-0，D2-0）发生的reads不足（<10,000个读数）。是以，我们在40个PCR轮回中对这4个样本反复16S PCR，并在随后的MiSeq测序运行中对其文库进行从新测序。

细菌16S rRNA基因剖析

Bacterial 16S rRNA gene analysis

用QIIME v1.8.0进行16S剖析。序列采用Q20进行质量掌握，然后拆分样本。基于Greengenes数据库（version 13.8）将序列分派到97％统一性阈值的有参（closed ）OTU。在处理16S rRNA基因序列直至OTU分类学分派时代，研究人员对组分派（group allocation）（即岁数，双胞胎，时间点）不知情。为了考虑样本间测序深度的变异性（variability），所有样本都稀少到每个样本10000条序列（10次反复（iterations）），跨越了之前描述的遍及接管的最小测序深度。在所有反复（iterations）中的16S剖析中获得了一致的究竟。是以，显露了从代表性反复（iterations）获得的究竟。使用QIIME较量α多样性（Faith’s phylogenetic diversity），OTU雄厚度和UniFrac距离。使用vegan R包绘制稀少曲线，使用500次置换（permutations）。使用Emperor可视化PCoA图。

相关收集

Correlation network

使用线性夹杂模型研究细菌和噬菌体随时间的相对品貌转变情形。夹杂线性模型考虑了在分歧时间点来自沟通受试者的反复测量。在该模型中，对细菌和噬菌体的相对品貌进行了对数转换，将样品收集时间指定为固定效应，并将受试者指定为随机效应。使用错误发现率（false discovery rate，q值）校正来自多个对照的P值。q值小于0.05被认为是统计学显著。使用Cytoscape绘制细菌和噬菌体相关性。剖析在R（version 3.1.2）中进行。

统计剖析

Statistics

使用SAS进行学生t磨练（双尾）。验证了组之间的正态分布和方差齐性。应用Wilcoxon磨练（成对，无参）成家的样品在0个月与24个月之间对照真核病毒雄厚度，噬菌体雄厚度，噬菌体多样性，细菌雄厚度和细菌多样性。用GraphPad Prism进行Wilcoxon磨练，线性回来和Spearman相关关系剖析。响应地显露P值和95％置信区间。

数据接见编号

Accession codes

序列数据已经留存到NCBI Sequence Read Archive中，BioProject的挂号号码为SRP058399。

译者简介

秋芒树，本科卒业于中国农业大学，硕士就读于英国帝国理工学院Computational Methods in Ecology and Evolution专业。存眷婴儿肠道微生物，肠道细菌与噬菌体的互相感化。在宏基因组公家号揭橥《Nature：TEDDY规划中幼儿肠道微生物组随时间的发育》、《Cell子刊：成年同卵双胞胎的病毒组多样性与肠道微生物组多样性相关》、《CHM：新生儿肠道微生物菌群研究》 等。迎接指摘、斧正和交流， j.wu18@imperial.ac.uk

Reference

1. Lim E S, Zhou Y, Zhao G, et al. Early life dynamics of the human gut virome and bacterial microbiome in infants[J]. Nature medicine, 2015, 21(10): 1228.

延伸阅读

1. Minot, Samuel, et al. “Rapid evolution of the human gut virome.” Proceedings of the National Academy of Sciences 110.30 (2013): 12450-12455.

2. Breitbart, Mya, et al. “Viral diversity and dynamics in an infant gut.” Research in microbiology 159.5 (2008): 367-373.

3. Bäckhed, Fredrik, et al. “Dynamics and stabilization of the human gut microbiome during the first year of life.” Cell host & microbe 17.5 (2015): 690-703.

婴儿生命早期肠道病毒组和细菌微生物组的动态

Early life dynamics of the human gut virome and bacterial microbiome in infants

翻译：秋芒树 英国帝国理工学院

责编：刘永鑫 中科院遗传发育所

原文链接：https://www.nature.com/articles/nm.3950

Nature Medine [IF = 32.621]

DOI: 10.1038/nm.3950

作者：Efrem S Lim^1，2, Yanjiao Zhou^3，4, Guoyan Zhao¹, Irma K Bauer³, Lindsay Droit^1，2, I Malick Ndao³,
Barbara B Warner³, Phillip I Tarr^1，3, David Wang^1，2 & Lori R Holtz³

1. 圣路易斯华盛顿大学，医学院，分子微生物系（Department of Molecular Microbiology, Washington University School of Medicine, St. Louis, Missouri, USA.）

2. 圣路易斯华盛顿大学，医学院，病理与免疫学系（Department of Pathology & Immunology,
   Washington University School of Medicine, St. Louis, Missouri, USA.）

3. 圣路易斯华盛顿大学，医学院，儿科系（Department of Pediatrics, Washington University School of Medicine, St. Louis, Missouri, USA.）

4. 杰克逊基因组医学实验室和布莱根妇女病院（The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut, USA, and Channing Division of Network Medicine, Department of
   Medicine, Brigham and Women’s Hospital, Boston, Massachusetts, USA.）

通信作者：
Lori R Holtz (Holtz_L@kids.wustl.edu)

David Wang (davewang@wusm.wustl.edu)

Resource： 2015-09-14

热心肠日报

链接：https://www.mr-gut.cn/papers/read/1079299880

Nature Medicine：婴儿肠道病毒和细菌组的动态转变

创作：赵弘烨 审核：蓝灿辉 2017年01月07日

判定一个纵向健康双胞胎婴儿队列的肠道病毒组和细菌组：

①比起不相关的婴儿，孪生兄弟姐妹之间的病毒组和细菌组较为相似；

②从出生到2岁，真核生物病毒组、细菌组都邑扩张，但这陪伴着噬菌体病毒组构成的萎缩和移位；

③出生后，噬菌体-细菌之间的关系即起头，并具有高捕食者-低被捕食者的状况，这相符Lotka-Volterra捕食模型；

④比起成年人相对不乱的菌群，婴儿菌群是高度转变的，并与跟着岁数转变的细菌、病毒和噬菌体构成有关。

扫一扫，存眷本公家号！

肠菌与健康 microbiota & health

带你认识肠道菌群与健康的奥秘