高级检索+

繁殖障碍母驴子宫细菌多样性研究——基于16S高通量测序

Study on uterine bacterial diversity of female donkey with reproductive disorders——Based on 16S high-throughput sequencing

  • 摘要: 为了研究繁殖障碍母驴子宫细菌多样性,试验将11头3~5岁成年母驴按健康状况分为健康组(HG组)、子宫炎组(MG组)和屡配不孕组(SG组),分别采集每组的子宫分泌物,提取样本总DNA后对细菌16S rRNA的V3~V4可变区进行PCR扩增,利用Illumina MiSeq PE300平台进行高通量测序。对测序数据进行统计和OTU聚类并绘制稀释曲线,计算Alpha多样性指数(Chao、Ace、Shannon、Simpson、Coverage指数),并基于OTU聚类结果进行物种注释和统计、Beta多样性和菌群相对丰度差异分析。结果表明:共得到有效序列436 487条,其中HG组152 882条,MG组161 781条,SG组121 824条;有效序列长度主要分布在421~460 bp之间,平均序列长度为437 bp,其中HG组约为437 bp, MG组约为438 bp, SG组约为436 bp;根据序列97%相似水平划分得到的OTU数目为511个。SG组和MG组的Ace、Chao指数均高于HG组;Shannon指数以SG组为最高,MG组为最低;Simpson指数以MG组为最高,SG组为最低;三组Coverage指数均为0.99。三组细菌相对丰度高于1%的共10个门,占菌群相对丰度的95.74%~97.98%,分别为拟杆菌门(Bacteroidetes)、硬壁菌门(Firmicutes)、变形菌门(Proteobacteria)、放线菌门(Actinobacteria)、梭杆菌门(Fusobacteria)、Synergistetes菌门、Saccharibacteria菌门、螺旋菌门(Spirochaetae)、SR1_Absconditabacteria菌门和软皮菌门(Tenericutes),优势菌门为拟杆菌门、硬壁菌门、变形菌门和放线菌门;三组细菌相对丰度高于1%的共37个属,占菌群相对丰度的66.60%~78.67%,分别为紫单胞菌属(Porphyromonas)、丝状杆菌属(Filobacterium)、弯曲杆菌属(Campylobacter)、梭杆菌属(Fusobacterium)、Oceanivirga菌属、棒状杆菌属(Corynebacterium)、嗜胨菌属(Peptoniphilus)、Ezakiella菌属、活动弯曲杆菌属(Mobiluncus)、蛋白胨链球菌属(Peptostreptococcus)等;优势菌属为紫单胞菌属、丝状杆菌属、弯曲杆菌属和梭杆菌属。经PCA法分析HG组、MG组和SG组在属水平PC1轴上的菌群组成差异解释度为48.54%,在PC2轴上的菌群组成差异解释度为15.02%;经PCoA法分析HG组、SG组和MG组在属水平PC1轴上的菌群组成差异解释度为37.76%,在PC2轴上的菌群组成差异解释度为23.44%;三组总体上菌群组成相似,但MG组个别样本与其他样本菌群组成差异较大。HG组的纤毛菌属细菌相对丰度在三组中最高,与其他两组差异显著(P<0.05);SG组的棒状杆菌属和克里斯滕森菌属细菌相对丰度在三组中均最高,与其他两组差异显著(P<0.05);MG组的拟杆菌属细菌相对丰度在三组中最高,与其他两组差异显著(P<0.05)。说明不同健康状况下母驴子宫细菌多样性具有明显差异。

     

    Abstract: In order to study the bacterial diversity of female donkeys with reproductive disorders, eleven adult female donkeys aged 3-5 years were divided into healthy group(HG group), uterine inflammation group(MG group) and infertility group(SG group) according to their health conditions. The uterine secretions of each group were collected, and the V3-V4 variable region of bacterial 16 S rRNA was amplified by PCR after total DNA was extracted. High-throughput sequencing was performed using Illumina MiSeq PE300 platform. Statistical analysis and OTU clustering were performed on sequencing data, and dilution curves were drawn. Alpha diversity index(Chao, Ace, Shannon, Simpson and Coverage index) was calculated. Species annotation and statistics were performed based on OTU clustering results, and differences in Beta diversity and relative abundance of flora were analyzed. The results showed that a total of 436 487 valid sequences were obtained, including 152 882 sequences in HG group, 161 781 sequences in MG group and 121 824 sequences in SG group. The effective sequence length mainly ranged from 421 to 460 bp, with an average of 437 bp, including 437 bp in HG group, 438 bp in MG group and 436 bp in SG group. The number of OTUs was 511 according to 97% similarity level of sequences. Ace and Chao indexes in SG and MG groups were higher than those in HG group. The Shannon index was highest in SG group and lowest in MG group. Simpson index was the highest in MG group and the lowest in SG group. The Coverage index of the three groups was 0.99. There were 10 phyla with relative abundance higher than 1% in the three groups, accounting for 95.74%-97.98% of the relative abundance of bacteria. They were Bacteroidetes, Firmicutes, Proteobacteria, Actinobacteria, Fusobacteria, Synergistetes, Saccharibacteria, Spirochaetae, Absconditabacteria_SR1 and Tenericutes, and the dominant bacteria were Bacteroidetes, Firmicutes, Proteobacteria and Actinomycetes. There were 37 genera with relative abundance higher than 1% in the three groups, accounting for 66.60%-78.67% of the relative abundance of bacteria. They were Porphyromonas, Filobacterium, Campylobacter, Fusobacterium, Oceanivirga, Corynebacterium, Peptoniphilus, Ezakiella, Mobiluncus, Peptostreptococcus, etc. The dominant bacteria were violet Omonas, Filamentous bacillus, Campylobacter and Clostridium. PCA analysis showed that the difference in flora composition on the PC1 axis of HG group, MG group and SG group was 48.54%, and that on the PC2 axis was 15.02%. PCoA method was used to analyze the difference of flora composition on the PC1 axis of HG group, SG group and MG group, and the difference of flora composition on the PC2 axis was 37.76% and 23.44%, respectively. The flora composition of the three groups was similar on the whole, but the community composition of individual samples in MG group differed greatly from that of other samples. The relative abundance of Ciliomycetes in HG group was the highest among the three groups, which was significantly different from the other two groups(P<0.05). The relative abundance of Corynebacterium and Christensenia in SG group was the highest among the three groups, which was significantly different from the other two groups(P<0.05). The relative abundance of Bacteroidetes in MG group was the highest among the three groups, which was significantly different from the other two groups(P<0.05). The results indicated that the diversity of bacterial community in uterus of female donkeys under different health conditions was significantly different.

     

/

返回文章
返回