Inheritance and Establishment of Gut Microbiota in Chickens

In mammals, the microbiota can be transmitted from the placenta, uterus, and vagina of the mother to the infant. Unlike mammals, development of the avian embryo is a process isolated from the mother and thus in the avian embryo the gut microbial developmental process remains elusive. To explore the establishment and inheritance of the gut microbiome in the avian embryo, we used the chicken as the model organism to investigate the gut microbial composition in embryos, chicks, and maternal hens. We observed: (1) 28 phyla and 162 genera of microbes in embryos where the dominated genus was Halomonas (79%). (2) 65 genera were core microbiota in all stages with 42% and 62% gut microbial genera of embryo were found in maternal hen and chick, respectively. There was a moderate correlation (0.40) between the embryo and maternal, and 0.52 between the embryo and chick at the family level. (3) Gut microbes that are involved in substance metabolism, infectious disease, and environmental adaptation are enriched in embryos, chicks, and maternal hens, respectively. (4) 94% genera of gut microbial composition were similar among three different chicken breeds which were maintained under similar conditions. Our findings provide evidence to support the hypothesis that part of the microbial colonizers harbored in early embryos were inherited from maternal hens, and the gut microbial abundance and diversity were influenced by environmental factors and host genetic variation during development.

Introduction

The gastrointestinal tracts of mammals are populated by a dynamic and enormous microbial community that can be viewed as complex, and polygenic traits interacting and coevolving with host genetic and environmental factors (Ley et al., 2006; Sansonetti and Medzhitov, 2009; Yang et al., 2017). The development of metagenomics has allowed for intensive study of the microbial genome and revealed that gut microbiota play important roles in the physiology and immunity of the host (Clemente et al., 2012; Kostic et al., 2015), by enhancing digestive efficiency, promoting immune system development and immune homeostasis, and limiting pathogen colonization (Fraune and Bosch, 2010; Oakley et al., 2014). In mammals, it was previously thought that fetuses lived in a sterile environment and the initial gut microbiota of the infant originated from the mother’s birth canal and body (Tissier, 1900). Recently, more evidence supports the idea that the gut microbiota can be vertically transmitted from the mother to her infant (Blaser, 2006). For example, microorganisms isolated from meconium of healthy newborns by Cesarean section indicated that the microbiota was not derived only postnatally (Jimenez et al., 2008; Moles et al., 2013; Ardissone et al., 2014). Many cultivable microorganisms present in the umbilical cord blood of preterm infants and abundant nonculturable and unclassified microorganisms have been detected in amniotic fluid, suggesting that fetuses are not sterile before delivery (Jimenez et al., 2005; DiGiulio et al., 2008). Morphologic and metagenomic studies have demonstrated that microbes exist in different regions of placenta and microbial DNA can be horizontally transferred from mother to fetus via placenta (Satokari et al., 2009; Stout et al., 2013; Aagaard et al., 2014). Additionally, the construction and succession of gut microbiota were influenced by many factors, such as the mode of delivery, birth environment, and feeding patterns (Dominguez-Bello et al., 2010; Bäckhed et al., 2015).

However, unlike mammals, the avian embryo is an isolated unit whose fertilization occurs in vivo, while development occurs in an extracorporeal egg without an umbilical cord, placenta, and amniotic fluid directly associated with the maternal body. The establishment and inheritance of a microbiome in an avian embryo is unknown. The chicken, domesticized from the red jungle fowl, bridges an evolutionary gap between reptiles and mammals, and is an important economical animal that provides humans with meat, eggs, and feathers. It was the first farm animal to have its genome sequenced (International Chicken Genome Sequencing Consortium, 2004). Chickens also serve as a traditional model for studying embryonic development, since its embryos are packaged in eggshells and can be easily observed and manipulated in vitro (Brown et al., 2003). Here, we used the chicken as a model organism and performed 16S rRNA sequencing to analyze and compare the microbial composition, abundance, and dynamic distribution during different embryonic stages of development, as well as chicks and their maternal hens (Supplementary Figure S1). In this study, we addressed the following: “(1) Can bacteria be detected in the chicken embryo? (2) If so, can we infer if any of these bacteria are maternally inherited? (3) Are there any core bacterial taxa that are shared across host developmental stages? (4) How does the gut microbiome change during development in a variety of host genetic background?”

Results

Profiles of Embryonic Microbiota

Fifty-one samples, including whole embryos that were incubated for 4 days, and intestines from embryos that were incubated for 19 days, were collected in all breeds. The microbial genomic DNA was isolated from ground embryos or intestines. The microbial classifications revealed that 28 phyla, 162 genera, and 76 species were present in the 4-day and the 19-day embryos (Supplementary Table S1). The most abundant phylum was Proteobacteria (86%), followed by Firmicutes (5%), Bacteroidetes (4%), and Actinobacteria (3%) in embryos (Supplementary Figure S2A). The dominant microbial genera were Halomonas (79%) and Ochrobactrum (5%), which belong to the phylum of Proteobacteria (Figures 1A, 2C). The correlation coefficient between the microbiota of the 4-day (E4) and the 19-day (E19) embryos was 0.80 (Table 1). Although there was similarity of microbial species in phylum between 4-day and the 19-day embryos, we observed a decreasing proportion of Proteobacteria, and remarkable increasing tendency of Firmicutes, Actinobacteria, and Bacteroidetes (Figure 1A). Accordingly, Halomonas was significantly different between the two groups, which were 88% and 55%, respectively. Microbial beta diversity of embryos with a Non-metric Multidimensional Scaling (NMDS; unweighted UniFrac distance) plot and heatmap showed the characteristics of the 4-day and 19-day indicating a significant different of microbial communities (Figures 1B,D). Signature genera at 4-day included bacteria found in 19-day, i.e., Halomonas, Bacteroidetes, and Ochrobactrum, in addition, alpha diversity analysis revealed more microbial diversity in the 19-day embryos than the 4-day embryos (Figures 1C, 2A and Supplementary Figure S3).

FIGURE 1

FIGURE 1. Aggregate microbiota characteristics of 4-day (E4) and 19-day (E19) embryos. (A) Dominant taxonomic groups of embryos by phylum. (B) Microbial beta diversity of embryos with a Non-metric Multidimensional Scaling (NMDS) plot showing how distant E4 and E19 communities were. (C) Microbial alpha diversity with a box plot exhibiting the community diversity (The Shannon estimator). (D) Heatmap of hierarchy cluster results for the statistical significant microbial OTUs of two groups at the genus level.

FIGURE 2

FIGURE 2. Phylogenetic tree of the taxa and differences in the microbial communities of different development stages. (A) Phylogenetic tree constructed from the 188 taxa. Colored blocks in the outermost circle indicate phyla, and in the inner circle indicate genera. The heatmap circles show relative abundance of embryos (E) E4 and E19; chicks (L) at 4 (L4), 21 (L21), and 42 days (L42) post hatch; and maternal hens (H). (B) The Venn diagram shows the microbes shared within different host development stages (see Supplementary Table S3 in the Supplementary Material). (C) Microbial composition structure of embryos, chicks, and maternal hens at the genus level.

TABLE 1

TABLE 1. Gut microbial correlation coefficients at different host development stages.

Composition of Gut Microbiome in Different Host Developmental Stages

In order to investigate the establishment and inheritance of gut microbiota during host development, we compared in three breeds the composition and abundance of gut microbiota in maternal hens (H), embryos (E) E4 and E19, and chicks (L) at 4 (L4), 21 (L21), and 42 days (L42) post hatch (Figure 2A and Supplementary Figure S1). One thousand seven hundred and nine microbial operational taxonomic units (OTUs) were detected in at least 75% of the samples, among them, 33 phyla and 196 genera of gut microbes were annotated (Supplementary Table S2). The most abundant phyla in all stages were Proteobacteria, Firmicutes, and Bacteroidetes, followed by Actinobacteria, Cyanobacteria, and Synergistetes (Figure 2A), which were consistent with previous studies (Waite and Taylor, 2014). However, the proportion of these phyla throughout the stages was different. Proteobacteria was mainly present in embryos (86%) and maternal hens (22%). Firmicutes (76% and 44%) and Bacteroidetes (10% and 24%) were the two major phyla in chicks and maternal hens, respectively (Supplementary Figure S2A). At the genus level, Halomonas was the dominant microbe in embryos (79%), and Lactobacillus (67%) and Bacteroides (10%) were the most abundant microbes in chicks (Figure 2C). Fifty-nine genera of microbes with low abundance were detected only in embryos. The declining population of those microbes during embryonic development and an absence in chicks and maternal hens suggested that some microbes were temporarily harbored in embryos and they were influenced by genetic and environmental factors during host development (Figures 2A,B and Supplementary Table S3). Abundance of Enterococcus, Ruminococcus, and Oscillospira was high in chicks and low in embryos and maternal hens. Jonquetella, Sneathia, and Porphyromonas whose function related to physiological and environmental adaptation (Jumas-Bilak et al., 2012; Belibasakis et al., 2013) showed abundant proportions in maternal hens (Figure 2A). Among the 162 genera that were detected in embryos, 68 and 100 of them were detected in maternal hens and chicks, respectively. The correlation coefficients for microbiota in embryos and maternal hens were approximately 0.33 and 0.40 at the genus and family level, respectively, and for embryos and chicks, it decreased initially and subsequently increased during the development of host stages (Table 1) with the correlation coefficient of 0.52 at the family level. Of the132 genera detected in chicks 91 were also observed in the maternal hens (Figure 2B). Correlation coefficients ranging from 0.58 to 0.72 between chicks and maternal hens were found at the family level (Table 1). Principle component analysis (PCA) indicated that the gut microbial composition in embryos was different from that in chicks and maternal hens, while similar between chicks and maternal hens (Supplementary Figure S2B). These results support the hypothesis that early colonizers of embryos were inherited from their maternal hens and a transmission process of gut microbiota from maternal hen to embryo and embryo to chick exists. The higher microbial diversity in embryos than in chicks and maternal hens indicates that the microbial diversity declined during host development, and provides additional evidence for the hypothesis that the harbored microbiome is influenced by the host during development (Figure 2A and Supplementary Figure S2C).

Among 196 annotated genera, 65 were considered as core microbes which existed in embryos, chicks, and maternal hens (Figures 2A,B and Table 2), such as Halomonas, Lactobacillus, Bacteroides, and Enterococcus (Figure 3). Halomonas had a high abundance in embryos and maternal hens, with a decreasing trend in chicks. This suggested that it may be an essential microbe and it may be susceptible to environmental effects. Lactobacillus and Enterococcus exhibited an increasing trend from embryo to chick and stabilized in maternal hens. Bacteroides remained steady in all stages (Figure 3). These results suggest that most of the microbial colonizers harbored in early embryos may originate from the maternal hens, and can be influenced by host genetic and environmental factors during developmental stages.

TABLE 2

TABLE 2. Sixty-five core microbes which were common to embryo, chick, and maternal hen in Venn diagram.

FIGURE 3

FIGURE 3. The dynamic distribution of main core microbes within different host development stages including embryos (E) E4 and E19; chicks (L) at 4 (L4), 21 (L21), and 42 days (L42) post hatch; and maternal hens (H) based on microbial abundance. The boxes have lines at the lower quartile, median, and upper quartile values.

Functional Performance of Microbiota in Different Host Developmental Stages

Comparing the functional capacity of the gut microbiota, in different host developmental stages including all breeds, could help us understand the relationships between the gut microbiota and its host. Twenty-nine abundant secondary metabolism pathways were identified while 25 pathways showed significant differences among embryos, chicks, and maternal hens in all breeds (Supplementary Table S4). As the gut microbiota of embryos evolve, they may develop more complex functional pathways such as substance metabolism, cell motility, signal transduction, cellular processes and signaling, cofactors, and vitamins metabolism. Most of these pathways are related to growth and development. Pathways allocated for translation, replication and repair, nucleotide metabolism, genetic information processing, carbohydrate metabolism, and transcription were more enriched in the microbiome harbored in chicks and maternal hens (Figure 4). The microbiota in chicks have a high abundance of functional capacities involved in metabolic pathways such as infection diseases, energy metabolism, and the digestive system (Figure 4C). Fewer differences in metabolism pathways between chicks and maternal hens indicated that, although they have a different gut microbial composition, their gut microbiome performs similar functions. Photosynthesis-antenna proteins usually exist in phycobilisomes of cyanobacteria and act as peripheral antenna systems enabling more efficient absorption of light energy. In our study, photosynthesis-antenna proteins, which fall into the category of energy metabolism, were enriched in embryos but not in chicks nor maternal hens (Supplementary Figure S4). Calcium signaling pathways are known to participate in bone formation and brain development in fetuses and in increasing bone mineral content of elders (Maas et al., 2005; Chen et al., 2010). We believe that the high abundance of microorganisms with calcium signaling pathways in the microbiota from embryos may be associated with growth and development. Microbiota whose function related to the replication and repair systems, such as nucleotide excision repair and non-homologous end-joining, which are involved in DNA repair to resist infectivity damage, were significantly enriched in chicks and maternal hens. The phosphotransferase system (PTS) is a major mechanism used by bacteria for uptake of carbohydrates, particularly hexoses, hexitols, and disaccharides, to enable the use of phosphoenolpyruvate as a source of energy. Microbiota related to this system were also enriched in chicks and maternal hens (Supplementary Figure S4). This phenomenon may be caused by the food intake of chicks and maternal hens and suggest that there is an inevitable relationship between growth and development of the host and microbial functional capacity.

FIGURE 4

FIGURE 4. Significant differences in microbial metabolic pathways at different stages. Extended error bar plot indicates the difference in mean proportion of microbial metabolic pathways between the two groups along with the associated confidence interval of this effect size and the p-value of Welch’s t-test (p < 0.05). (A) Embryo vs. chick. (B) Embryo vs. maternal hen. (C) Chick vs. maternal hen.

Host Genetic Background Influences the Gut Microbiota

In the above results, microbial community and functional capacity were analyzed in different host developmental stages including all breeds. To investigate the extent of the influence of host genetic factors on gut microbiota, we also studied microbial composition and abundance in three native Chinese breeds: Beijing Fatty (B), Shiqiza (C), and Xianju (X) chickens. 16S rRNA annotation results revealed that, under similar feeding conditions, the three breeds had similar dominant microorganism communities at the genus level: Lactobacillus, Halomonas, Bacteroides, and Enterococcus (Supplementary Figures S5A,B). Alpha diversity results indicated that there was no difference in the richness and diversity of the gut microbiota among the Beijing Fatty, Shiqiza, and Xianju chickens (Supplementary Figure S5D). Ninety-four percent genera (185/196) from the gut microbiota were common in Beijing Fatty, Shiqiza, and Xianju chickens (Supplementary Figure S5C). Although the dominant microorganisms were uniform among the breeds, 109 genera were significantly different (p < 0.05) among Beijing Fatty, Shiqiza, and Xianju chickens (Supplementary Table S5). Principal component analysis showed that the samples clustered according to breeds within different host developmental stages, and demonstrated that the Beijing Fatty, Shiqiza, and Xianju separated clear in E4 and E19 embryos (Figure 5). However, these three breeds were more dispersed and similar in chicks and maternal hens compared to embryos indicating the similar feeding conditions influence the construction of host microbiota (Figure 5). The microbiota were still influenced by breed, albeit with more similarities in chicks and maternal hens.

FIGURE 5

FIGURE 5. Microbiota compositions for Beijing Fatty (B), Shiqiza (C), and Xianju (X) of chickens within different host development stages including embryos (E) E4 and E19; chicks (L) at 4 (L4), 21 (L21), and 42 days (L42) post hatch; and maternal hens (H).

Discussion

Studies of gut microbial establishment have been performed most in adult animals or newborn offspring, such as human neonates (Mackie et al., 1999), insects (Tang et al., 2012), zebra fish (Davis et al., 2016), chickens (Ding et al., 2016), snakes (Costello et al., 2010; McLaughlin et al., 2015), alligators (Keenan et al., 2013; Keenan and Elsey, 2015), and lizards (Ren et al., 2016; Kohl et al., 2017), while few studies involved the fetus or embryo. In our study, microbiota were detected qualitatively in the embryo. The formation of the egg is an open process where large amounts of lipids and proteins are transported to the oocyte in the maternal hen ovary. When the oocyte is ovulated, fertilization occurs in the funnel of the infundibulum, and then the oosperm is packaged with secretory proteins in the oviduct to form a fertilized egg. Although there is no umbilical cord, placenta, or amniotic fluid directly associated with the maternal hen, avian embryos obtain abundant nutrients from maternal hen. Mammalian studies have suggested that the upper female reproductive tract is an essential place for the beginning of life. In a recent study, researchers found different types of bacteria in the ovaries and fallopian tubes which suggested that the upper female reproductive tract was not free from microorganisms (Keku, 2016). Additional evidence indicated that vaginal tract contains a remarkably complex microbial community which can be transmitted from mother to fetus (Mackie et al., 1999; Schultz et al., 2004; White et al., 2011; Martin, 2012; Prince et al., 2014). The oviduct of chickens joins the urinary and digestive tracts in the cloaca which may result in the transmission of microbiota from the cloacae through the oviduct and even the body cavity (Gross and Siegel, 1959; Stanley et al., 2015; Zhao et al., 2015; Pedroso et al., 2016). So, it is a normal phenomenon to detect microbes in the embryo. Also, there is possible that the microbes were from the environment, chemical reagent, or the process of sample collection (Salter et al., 2014). However, in order to verify whether the microbiota were contaminated by the environment or experimenter, we performed the preliminary experiments. Seven embryos were selected as the experimental group (B) and five embryo-free samples as the control group (C) for 16S rRNA sequencing. Comparing the microbial OTUs between the two groups, we found that there are much more abundant microbial taxa (Supplementary Figure S6A) in the experimental than the control group. The result of rarefaction curve also indicated the deficient sequencing depth in control group (Supplementary Figure S6B). Further analysis revealed significant differences in microbial diversity and abundance between the experimental and control groups (Supplementary Figure S6C). As suggested above, we speculated that the establishment and inheritance of the embryo gut microbiota were derived from the maternal hen by the process of fertilization and egg formation in the oviduct. Further studies are still needed to determine the microbial composition of avian reproductive system.

In order to adapt to their rapidly changing developmental environments, patterns of gut microbial diversity and abundance varied across different host development stages. However, some of the gut microbial membership in embryos, chicks, and maternal hens was strikingly similar, with a shared core microbiota. The definition of core microbiota has been addressed in fish and mammals (Turnbaugh et al., 2009; Roeselers et al., 2011; Tang et al., 2012), although the producing mechanisms of core microbiota remain unclear. Our study suggests that core gut microbiota also exists in chickens where 65 genera were detected and considered as core microbes in embryos, chicks, and maternal hens. Interestingly, Halomonas, which belongs to Halomonadaceae family and Proteobacteria phylum, was the dominant genus in embryos and maternal hens but not chicks (Figure 2A). Halomonadaceae are Gram-negative, chemoorganotrophic, aerobic, or facultative anaerobic, and moderately halophilic, haloalkaliphilic, halotolerant, or non-halophilic. Certain Halomonas species are also known to have a temperature-dependent salt requirement or tolerance and are able to utilize cellulose as carbon sources, although Borsodi et al. (2005) did not observe direct cellulolytic activity. Halomonas venusta was identified as the abundant species in the embryos and maternal hens. Berlanga et al. (2014) reported that H. venusta retained the capacity of sessile cells to adhere to polystyrene and form a biofilm. Ectoine, which is synthesized and released by H. venusta, can balance the cell’s osmotic pressure and protect the structures of enzymes, DNA, and the cytomembrane, implying that H. venusta may be involved in cell motility, amino acid, carbohydrate metabolism, and environmental adaptation in embryos and maternal hens (Lippert and Galinski, 1992; Zheng et al., 2010; Zhou et al., 2014). A strong correlation between bacterial growth and microbial utilization of phenol was discovered, and H. venusta was able to consume phenol as the sole carbon source indicating it can degrade toxic and recalcitrant compounds (Munoz et al., 2001). Compared with chicks and maternal hens, xenobiotics biodegradation and metabolism were significantly enriched in embryos of our study (Figures 4A,B and Supplementary Figure S7). As embryos produce poisonous and harmful substances during development, we speculate that the existence of H. venusta can assist the host by decomposing these toxic compounds. Therefore, the abundance of specific genera within the chicken core microbiota may be due, in part, to the host distinct growth and development, as well as diet in each developmental stage, i.e., the transition from embryos to chicks which utilizing the external feeds that are rich in carbohydrates (Noy et al., 2001; Zhao et al., 2013). Our results provide evidence for the hypothesis that the composition of embryo microbiome was inherited part from maternal hens and adjusted by host genetic and environmental factors during different developmental stages.

Materials and Methods

Study Design and Animal Sampling

The intent of this study was to characterize the gut microbial establishment and inheritance from maternal hen to embryo and to chick (Supplementary Figure S1). We selected 12 individual maternal hens at 35 weeks of age from three breeds (Beijing Fatty, Xianju, and Shiqiza) from the Animal Husbandry and Veterinary Research Institute in the Shanghai Academy of Agricultural Science for reproduction. Fertilized eggs from the selected maternal hens were used to hatch embryos and chicks which also include three different breeds (Supplementary Figure S1). The incubation was performed in Ova-Easy Advance Series II Digital Cabinet Egg Incubator (Brinsea, United Kingdom) within a sterilized room at 37.8°C and 55–65% humidity. To investigate embryo gut microbiota, we aseptically collected 51 samples including whole embryos that had incubated to day 4, and intestines from embryos that had incubated to day 19. The collected embryo samples were immediately frozen at -80°C. The procedures of embryonic sample collection and subsequent operation were carried out on a clean bench under aseptic conditions. The newly hatched chicks were maintained in same husbandry conditions with their maternal hen. Fecal samples from the maternal hens (n = 12) at reproduction and chicks (n = 113) aged 4, 21, and 42 days were collected (Supplementary Table S6). All of the fecal specimens were kept on ice when collecting, transported to the laboratory and frozen at -80°C immediately until further analysis. The procedures of embryonic sample collection, microbial genomic DNA extraction, and PCR amplification were performed with sterile procedures. Protocols used for this experiment were approved by ethics committee for the Care and Use of Laboratory Animals in Shanghai Jiao Tong University, China.

DNA Extraction and 16S rRNA Sequencing

Genomic DNA was isolated from intestinal and fecal samples using TIANGEN DNA Stool Mini Kit (TIANGEN, cat#DP328) following the manufacturer’s instructions. Harvested DNA was quantified on a Nanodrop spectrophotometer (Thermo scientific, United States) to assess the quantity and quality of the DNA. The V4 hypervariable region of the 16S rRNA gene was amplified by PCR using sample-specific sequence barcoded fusion primers: forward 5′-AYTGGGYDTAAAGNG-3′, reverse 5′-TACNVGGGTATCTAATCC-3′. The PCR reaction conditions were: 94°C for 5 min; 94°C for 30 s, 50°C for 30 s, and 72°C for extension; repeated for 27 cycles; with a final 72°C for 7 min (Zhao et al., 2013). The PCR products were excised from a 1.5% agarose gel and purified using a QIAGEN quick Gel Extraction Kit (QIAGEN, cat#28706). Purified PCR products from 176 samples were used to construct a sequencing library using Illumina TruSeq (Illumina, United States) following the manufacturer’s suggested protocols. 16S rRNA sequencing was performed at the Shanghai Personal Biotechnology Limited Company, Shanghai, China, using the Illumina MiSeq (Illumina, United States) sequencing platform. Barcoded V4 amplicons were sequenced using the pair-end method by Illumina Miseq. Original pair-end sequences with mean quality lower than 30, containing ambiguous bases, sequence length shorter than 150 bp, chimera, adaptor contamination, or host contaminating reads were removed. The original pair-end sequence reads that passed our quality control criteria and contained a sequence overlap of at least 10 bp without any mismatch were assembled according to their overlap sequence. A total of 9,887,520 sequences from the V4 region of the 16S rRNA sequence from 176 samples that passed our quality filters were used for this study (Supplementary Table S7), with an average length of 225 bp for each sequence (Supplementary Figure S8). Trimmed sequences were uploaded to QIIME for further analysis. The metagenome sequences used in this paper are publicly available from Metagenome Rapid Annotation using Subsystem Technology¹.

Taxonomic Assignment and Statistical Analysis

Microbial OTUs were derived from the trimmed sequences of the PCR amplicon for the V4 hypervariable region of the 16S rRNA gene which were compared to the Greengene (DeSantis et al., 2006) databases using the uclust and blast functions in QIIME (Caporaso et al., 2010). 76,840 OTUs were annotated from the 9,887,520 amplicons and classified at the domain, phylum, class, order, family, genus, and species levels (Supplementary Table S8). Then, OTUs that were present in at least 44 samples were used for statistical analysis. The OTU abundance count was log2 transformed and then normalized as follows: from each log-transformed measure, the arithmetic mean of all transformed values was subtracted; the difference was divided by the standard deviation of all log-transformed values for the given sample. After this procedure, the abundance profiles for all samples exhibited a mean of 0 and a standard deviation of 1 (Zhao et al., 2013). Normalized abundance was used to generate a heatmap by Cluster3.0 and Java Treeview (Eisen, 2002) and to draw a 3dPCA figure by R (Dessau and Pipper, 2008). Alpha-diversity analysis was performed in maternal hens (Schloss et al., 2009) with the alpha-diversity.py script to calculate the chao1 (Pitta et al., 2014) and Shannon (Mahaffee and Kloepper, 1997) metrics. The Venn diagrams were generated using Venny (Oliveros, 2007). Box plots and bar charts were performed by SigmaPlot (Kornbrot, 2000) and STAMP (Parks and Beiko, 2010), and all p-values were adjusted by the Benjamini–Hochberg FDR procedure (Benjamini and Hochberg, 1995). The phylogenetic tree was generated by GraPhlAn, which is a software tool for producing high-quality circular representations of taxonomic and phylogenetic trees (Asnicar et al., 2015). Correlation analysis was conducted to identify the association of microbiome among different developmental stages using Microsoft Excel (Pearson, 1901).

Functional Profile Analysis

Microbial functional profile was predicted using PICRUSt (Langille et al., 2013). The OTUs were mapped to gg13.5 database at 97% similarity by QIIME’s command “pick_closed_otus.” The OTU abundance was normalized automatically using 16S rRNA gene copy numbers from known bacterial genomes in Integrated Microbial Genomes (IMG). The predicted genes and their functions were aligned to Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2014) database and the differences among groups were compared using the software STAMP. Two-sided Welch’s t-test (Welch, 1947) and Benjamini–Hochberg FDR correction were used in group analysis.

Ethics Statement

Protocols used for this experiment were approved by ethics committee for the Care and Use of Laboratory Animals in Shanghai Jiao Tong University, China.

Author Contributions

JD, YZ, and HM wrote the manuscript. JD and HM conceived and designed the experimental procedure, and supervised the study. JD performed statistical analysis, sample collection, DNA extraction, and data processing. RD, LY, and CH participated in sample collection, statistical analysis, and figure design. KX, LX, and LL coordinated sample collection. All authors read and approved the final manuscript.

Funding

This study was supported by the National key Research and Development Program of China (grant no. 2017YFD0500506) and the National Science Foundation of China (grant no. 31572384).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.01967/full#supplementary-material

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Footnotes

1. ^ http://metagenomics.anl.gov/mgmain.html?mgpage=project&project=20686

References

Aagaard, K., Ma, J., Antony, K. M., Ganu, R., Petrosino, J., and Versalovic, J. (2014). The placenta harbors a unique microbiome. Sci. Transl. Med. 6, 237ra65. doi: 10.1126/scitranslmed.3008599

PubMed Abstract | CrossRef Full Text | Google Scholar

Ardissone, A. N., De La Cruz, D. M., Davis-Richardson, A. G., Rechcigl, K. T., Li, N., Drew, J. C., et al. (2014). Meconium microbiome analysis identifies bacteria correlated with premature birth. PLOS ONE 9:e90784. doi: 10.1371/journal.pone.0090784

PubMed Abstract | CrossRef Full Text | Google Scholar

Asnicar, F., Weingart, G., Tickle, T. L., Huttenhower, C., and Segata, N. (2015). Compact graphical representation of phylogenetic data and metadata with GraPhlAn. PeerJ 3:e1029. doi: 10.7717/peerj.1029

PubMed Abstract | CrossRef Full Text | Google Scholar

Bäckhed, F., Roswall, J., Peng, Y., Feng, Q., Jia, H., Kovatcheva-Datchary, P., et al. (2015). Dynamics and stabilization of the human gut microbiome during the first year of Life. Cell Host Microbe 17, 690–703. doi: 10.1016/j.chom.2015.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Belibasakis, G. N., Ozturk, V. O., Emingil, G., and Bostanci, N. (2013). Synergistetes cluster A in saliva is associated with periodontitis. J. Periodontal Res. 48, 727–732. doi: 10.1111/jre.12061

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57, 289–300.

Google Scholar

Berlanga, M., Domenech, O., and Guerrero, R. (2014). Biofilm formation on polystyrene in detached vs. planktonic cells of polyhydroxyalkanoate-accumulating Halomonas venusta. Int. Microbiol. 17, 205–212. doi: 10.2436/20.1501.01.223

PubMed Abstract | CrossRef Full Text | Google Scholar

Blaser, M. J. (2006). Who are we? Indigenous microbes and the ecology of human diseases. EMBO Rep. 7, 956–960. doi: 10.1038/sj.embor.7400812

PubMed Abstract | CrossRef Full Text | Google Scholar

Borsodi, A. K., Micsinai, A., Rusznyak, A., Vladar, P., Kovacs, G., Toth, E. M., et al. (2005). Diversity of alkaliphilic and alkalitolerant bacteria cultivated from decomposing reed rhizomes in a Hungarian soda lake. Microb. Ecol. 50, 9–18. doi: 10.1007/s00248-004-0063-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, W. R., Hubbard, S. J., Tickle, C., and Wilson, S. A. (2003). The chicken as a model for large-scale analysis of vertebrate gene function. Nat. Rev. Genet. 4, 87–98. doi: 10.1038/nrg998

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). Qiime allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. H., Liu, C., You, L., and Simmons, C. A. (2010). Boning up on Wolff’s Law: mechanical regulation of the cells that make and maintain bone. J. Biomech. 43, 108–118. doi: 10.1016/j.jbiomech.2009.09.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Clemente, J. C., Ursell, L. K., Parfrey, L. W., and Knight, R. (2012). The impact of the gut microbiota on human health: an integrative view. Cell 148, 1258–1270. doi: 10.1016/j.cell.2012.01.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Costello, E. K., Gordon, J. I., Secor, S. M., and Knight, R. (2010). Postprandial remodeling of the gut microbiota in Burmese pythons. ISME J. 4, 1375–1385. doi: 10.1038/ismej.2010.71

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, D. J., Bryda, E. C., Gillespie, C. H., and Ericsson, A. C. (2016). 16S rRNA amplicon sequencing dataset for conventionalized and conventionally raised zebra fish larvae. Data Brief. 5, 938–943. doi: 10.1016/j.dib.2016.06.057

PubMed Abstract | CrossRef Full Text | Google Scholar

DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05

PubMed Abstract | CrossRef Full Text | Google Scholar

Dessau, R. B., and Pipper, C. B. (2008). ”R”–project for statistical computing. Ugeskr. Laeger. 170, 328–330.

Google Scholar

DiGiulio, D. B., Romero, R., Amogan, H. P., Kusanovic, J. P., Bik, E. M., Gotsch, F., et al. (2008). Microbial prevalence, diversity and abundance in amniotic fluid during preterm labor: a molecular and culture-based investigation. PLOS ONE 3:e3056. doi: 10.1371/journal.pone.0003056

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, J., Zhao, L., Wang, L., Zhao, W., Zhai, Z., Leng, L., et al. (2016). Divergent selection-induced obesity alters the composition and functional pathways of chicken gut microbiota. Genet. Sel. Evol. 48, 93. doi: 10.1186/s12711-016-0270-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Dominguez-Bello, M. G., Costello, E. K., Contreras, M., Magris, M., Hidalgo, G., Fierer, N., et al. (2010). Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc. Natl. Acad. Sci. U.S.A. 107, 11971–11975. doi: 10.1073/pnas.1002601107

PubMed Abstract | CrossRef Full Text | Google Scholar

Eisen, M. (2002). Gene Cluster Treeview. Stanford, CA: Stanford University.

Google Scholar

Fraune, S., and Bosch, T. C. (2010). Why bacteria matter in animal development and evolution. Bioessays 32, 571–580. doi: 10.1002/bies.200900192

PubMed Abstract | CrossRef Full Text | Google Scholar

Gross, W. B., and Siegel, P. B. (1959). Coliform peritonitis of chickens. Avian Dis. 3, 370–373. doi: 10.2307/1587575

CrossRef Full Text | Google Scholar

International Chicken Genome Sequencing Consortium (2004). Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature 432, 695–716.

Google Scholar

Jimenez, E., Fernandez, L., Marin, M. L., Martin, R., Odriozola, J. M., Nueno-Palop, C., et al. (2005). Isolation of commensal bacteria from umbilical cord blood of healthy neonates. Curr. Microbiol. 51, 270–274. doi: 10.1007/s00284-005-0020-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Jimenez, E., Marin, M. L., Martin, R., Odriozola, J. M., Olivares, M., Xaus, J., et al. (2008). Is meconium from healthy newborns actually sterile? Res. Microbiol. 159, 187–193. doi: 10.1016/j.resmic.2007.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Jumas-Bilak, E., Carlier, J. P., Jean-Pierre, H., Citron, D., Bernard, K., Damay, A., et al. (2012). Jonquetella anthropi gen. nov., sp. nov., the first member of the candidate phylum ’Synergistetes’ isolated from man. Int. J. Syst. Evol. Microbiol. 57, 2743–2748. doi: 10.1099/ijs.0.65213-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Goto, S., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2014). Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 42, D199–D205. doi: 10.1093/nar/gkt1076

PubMed Abstract | CrossRef Full Text | Google Scholar

Keenan, S. W., and Elsey, R. M. (2015). The good, the bad, and the unknown: microbial symbioses of the American alligator. Integr. Comp. Biol. 55, 972–985. doi: 10.1093/icb/icv006

PubMed Abstract | CrossRef Full Text | Google Scholar

Keenan, S. W., Engel, A. S., and Elsey, R. M. (2013). The alligator gut microbiome and implications for archosaur symbioses. Sci. Rep. 3:2877. doi: 10.1038/srep02877

PubMed Abstract | CrossRef Full Text | Google Scholar

Keku, T. O. (2016). “Bacteria found in female upper reproductive tract, once thought sterile,” in American Society of Clinical Oncology Annual Meeting, (Chicago, IL: University of North Carolina).

PubMed Abstract

Kohl, K. D., Brun, A., Magallanes, M., Brinkerhoff, J., Laspiur, A., Acosta, J. C., et al. (2017). Gut microbial ecology of lizards: insights into diversity in the wild, effects of captivity, variation across gut regions, and transmission. Mol. Ecol. 26, 1175–1189. doi: 10.1111/mec.13921

PubMed Abstract | CrossRef Full Text | Google Scholar

Kornbrot, D. (2000). Statistical software for microcomputers: SigmaPlot 2000 and SigmaStat2. Br. J. Math. Stat. Psychol. 53, 335–337. doi: 10.1348/000711000159268

PubMed Abstract | CrossRef Full Text | Google Scholar

Kostic, A. D., Gevers, D., Siljander, H., Vatanen, T., Hyotylainen, T., Hamalainen, A. M., et al. (2015). The dynamics of the human infant gut microbiome in development and in progression toward type 1 diabetes. Cell Host Microbe 17, 260–273. doi: 10.1016/j.chom.2015.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Langille, M. G. I., Zaneveld, J., Caporaso, J. G., Mcdonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676

PubMed Abstract | CrossRef Full Text | Google Scholar

Ley, R. E., Peterson, D. A., and Gordon, J. I. (2006). Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell 124, 837–848. doi: 10.1016/j.cell.2006.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Lippert, K., and Galinski, E. A. (1992). Enzyme stabilization be ectoine-type compatible solutes: protection against heating, freezing and drying. Appl. Microbiol. Biotechnol. 37, 61–65. doi: 10.1007/BF00174204

CrossRef Full Text | Google Scholar

Maas, J. W. J., Indacochea, R. A., Muglia, L. M., Tran, T. T., Vogt, S. K., West, T., et al. (2005). Calcium-stimulated adenylyl cyclases modulate ethanol-induced neurodegeneration in the neonatal brain. J. Neurosci. 25, 2376–2385. doi: 10.1523/JNEUROSCI.4940-04.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mackie, R. I., Sghir, A., and Gaskins, H. R. (1999). Developmental microbial ecology of the neonatal gastrointestinal tract. Am. J. Clin. Nutr. 69, 1035S–1045S.

Google Scholar

Mahaffee, W. F., and Kloepper, J. W. (1997). Temporal changes in the bacterial communities of soil, rhizosphere, and endorhiza associated with field-grown cucumber (Cucumis sativus L.). Microb Ecol. 34, 210–223. doi: 10.1007/s002489900050

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, D. H. (2012). The microbiota of the vagina and its influence on women’s health and disease. Am. J. Med. Sci. 343, 2–9. doi: 10.1097/MAJ.0b013e31823ea228

PubMed Abstract | CrossRef Full Text | Google Scholar

McLaughlin, R. W., Cochran, P. A., and Dowd, S. E. (2015). Metagenomic analysis of the gut microbiota of the Timber Rattlesnake, Crotalus horridus. Mol. Biol. Rep. 42, 1187–1195. doi: 10.1007/s11033-015-3854-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Moles, L., Gómez, M., Heilig, H., Bustos, G., Fuentes, S., De Vos, W., et al. (2013). Bacterial diversity in meconium of preterm neonates and evolution of their fecal microbiota during the first month of life. PLOS ONE 8:e66986. doi: 10.1371/journal.pone.0066986

PubMed Abstract | CrossRef Full Text | Google Scholar

Munoz, J. A., Perez-Esteban, B., Esteban, M., De La Escalera, S., Gomez, M. A., Martinez-Toledo, M. V., et al. (2001). Growth of moderately halophilic bacteria isolated from sea water using phenol as the sole carbon source. Folia Microbiol. 46, 297–302. doi: 10.1007/BF02815617

PubMed Abstract | CrossRef Full Text | Google Scholar

Noy, Y., Geyra, A., and Sklan, D. (2001). The effect of early feeding on growth and small intestinal development in the posthatch poult. Poult. Sci. 80, 912–919. doi: 10.1093/ps/80.7.912

CrossRef Full Text | Google Scholar

Oakley, B. B., Lillehoj, H. S., Kogut, M. H., Kim, W. K., Maurer, J. J., Pedroso, A., et al. (2014). The chicken gastrointestinal microbiome. FEMS Microbiol. Lett. 360, 100–112. doi: 10.1111/1574-6968.12608

PubMed Abstract | CrossRef Full Text | Google Scholar

Oliveros, J. C. (2007). An Interactive Tool for Comparing Lists with Venn’s Diagrams. Available at: http://bioinfogp.cnb.csic.es/tools/venny/index.html

Google Scholar

Parks, D. H., and Beiko, R. G. (2010). Identifying biologically relevant differences between metagenomic communities. Bioinformatics 26, 715–721. doi: 10.1093/bioinformatics/btq041

PubMed Abstract | CrossRef Full Text | Google Scholar

Pearson, K. (1901). Mathematical contributions to the theory of evolution. Parts 13 and 14. 13. On the theory of contingency and its relation to association and normal correlation. 14. On the general theory of skew correlation and non linear regression. Drapers Comp. Res. Mem. Biom. 191, 229–311.

Pedroso, A. A., Batal, A. B., and Lee, M. D. (2016). Effect of in ovo administration of an adult-derived microbiota on establishment of the intestinal microbiome in chickens. Am. J. Vet. Res. 77, 514–526. doi: 10.2460/ajvr.77.5.514

PubMed Abstract | CrossRef Full Text | Google Scholar

Pitta, D. W., Parmar, N., Patel, A. K., Indugu, N., Kumar, S., Prajapathi, K. B., et al. (2014). Bacterial diversity dynamics associated with different diets and different primer. PLOS ONE 9:e111710. doi: 10.1371/journal.pone.0111710

PubMed Abstract | CrossRef Full Text | Google Scholar

Prince, A. L., Antony, K. M., Ma, J., and Aagaard, K. M. (2014). The microbiome and development: a mother’s perspective. Semin. Reprod. Med. 32, 14–22. doi: 10.1055/s-0033-1361818

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, T., Kahrl, A. F., Wu, M., and Cox, R. M. (2016). Does adaptive radiation of a host lineage promote ecological diversity of its bacterial communities? A test using gut microbiota of Anolis lizards. Mol. Ecol. 25, 4793–4804. doi: 10.1111/mec.13796

PubMed Abstract | CrossRef Full Text | Google Scholar

Roeselers, G., Mittge, E. K., Stephens, W. Z., Parichy, D. M., Cavanaugh, C. M., Guillemin, K., et al. (2011). Evidence for a core gut microbiota in the zebrafish. ISME J. 5, 1595–1608. doi: 10.1038/ismej.2011.38

PubMed Abstract | CrossRef Full Text | Google Scholar

Salter, S. J., Cox, M. J., Turek, E. M., Calus, S. T., Cookson, W. O., Moffatt, M. F., et al. (2014). Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 12:87. doi: 10.1186/s12915-014-0087-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Sansonetti, P. J., and Medzhitov, R. (2009). Learning tolerance while fighting ignorance. Cell 138, 416–420. doi: 10.1016/j.cell.2009.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Satokari, R. J., Gronroos, T., Laitinen, K., Salminen, S., and Isolauri, E. (2009). Bifidobacterium and Lactobacillus DNA in the human placenta. Lett. Appl. Microbiol. 48, 8–12. doi: 10.1111/j.1472-765X.2008.02475.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Schultz, M., Gottl, C., Young, R. J., Iwen, P., and Vanderhoof, J. A. (2004). Administration of oral probiotic bacteria to pregnant women causes temporary infantile colonization. J. Pediatr. Gastroenterol. Nutr. 38, 293–297. doi: 10.1097/00005176-200403000-00012

PubMed Abstract | CrossRef Full Text | Google Scholar

Stanley, D., Geier, M. S., Chen, H., Hughes, R. J., and Moore, R. J. (2015). Comparison of fecal and cecal microbiotas reveals qualitative similarities but quantitative differences. BMC Microbiol. 15:51. doi: 10.1186/s12866-015-0388-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Stout, M. J., Conlon, B., Landeau, M., Lee, I., Bower, C., Zhao, Q., et al. (2013). Identification of intracellular bacteria in the basal plate of the human placenta in term and preterm gestations. Am. J. Obstet. Gynecol. 208, 226.e1–226.e7. doi: 10.1016/j.ajog.2013.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, X., Freitak, D., Vogel, H., Ping, L., Shao, Y., Cordero, E. A., et al. (2012). Complexity and variability of gut commensal microbiota in polyphagous lepidopteran larvae. PLOS ONE 7:e36978. doi: 10.1371/journal.pone.0036978

PubMed Abstract | CrossRef Full Text | Google Scholar

Tissier, H. (1900). Recherches sur la Flore Intestinale des Nourrissons (E’tatnormal et Pathologique). Paris: G. Carre and C. Naud.

Turnbaugh, P. J., Hamady, M., Yatsunenko, T., Cantarel, B. L., Duncan, A., Ley, R. E., et al. (2009). A core gut microbiome in obese and lean twins. Nature 457, 480–484. doi: 10.1038/nature07540

PubMed Abstract | CrossRef Full Text | Google Scholar

Waite, D. W., and Taylor, M. W. (2014). Characterizing the avian gut microbiota: membership, driving influences, and potential function. Front. Microbiol. 5:223. doi: 10.3389/fmicb.2014.00223

PubMed Abstract | CrossRef Full Text | Google Scholar

Welch, B. L. (1947). The generalisation of student’s problems when several different population variances are involved. Biometrika 34, 28–35.

Google Scholar

White, B. A., Creedon, D. J., Nelson, K. E., and Wilson, B. A. (2011). The vaginal microbiome in health and disease. Trends Endocrinol. Metab. 22, 389–393. doi: 10.1016/j.tem.2011.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Liu, S., Ding, J., Dai, R., He, C., Xu, K., et al. (2017). Gut microbiota co-microevolution with selection for host humoral immunity. Front. Microbiol. 8:1243. doi: 10.3389/fmicb.2017.01243

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, L. L., Wang, G., Siegel, P., He, C., Wang, H., Zhao, W. J., et al. (2013). Quantitative genetic background of the host influences gut microbiomes in chickens. Sci. Rep. 3:1163. doi: 10.1038/srep01163

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, W., Wang, Y., Liu, S., Huang, J., Zhai, Z., He, C., et al. (2015). The dynamic distribution of porcine microbiota across different ages and gastrointestinal tract segments. PLOS ONE 10:e0117441. doi: 10.1371/journal.pone.0117441

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, X., Ma, H., Yan, X. W., and Zhang, L. H. (2010). Ectoine synthesis and release under osmotic shock in halomonas venusta DSM4743. Microbiol. China 37, 1090–1096.

Google Scholar

Zhou, Y., Li, R., Gao, X. Y., Lapidus, A., Han, J., Haynes, M., et al. (2014). High quality draft genome sequence of the slightly halophilic bacterium Halomonas zhanjiangensis type strain JSM 078169(T) (DSM 21076(T)) from a sea urchin in southern China. Stand. Genomic Sci. 9, 1020–1030. doi: 10.4056/sigs.5449586

PubMed Abstract | CrossRef Full Text | Google Scholar