Genome and metagenome analyses reveal adaptive evolution of the host and interaction with the gut microbiota in the goose

Goose multi-omics database: A comprehensive multi-omics database for goose genomics

Multi-omics has helped elucidate the gene expression patterns and genomic variations closely associated with economically significant traits in geese. Despite the substantial genomic data generated through extensive goose studies, a unified platform for integrating these datasets is lacking. To address this gap, we introduced the Goose Multi-omics Database (GMD), which is accessible at http://goosedb.com/. The GMD is a comprehensive resource enabling streamlined search, analysis, and visualization of genetic information through a unified interface, providing insights into phenotypic traits, gene sequences, structures, expression profiles, genomic variations, gene families, homology, and collinearity. Equipped with robust analytical tools such as GBrowse and BLAST, the GMD facilitates rapid access to target gene information, significantly enhancing the efficiency and productivity of genomic research. By serving as a versatile and intuitive online repository, the GMD offers transformative potential for advancing goose biology, fostering multi-omics investigations, and integrating cutting-edge methodologies such as deep learning to accelerate discoveries in goose genomics.

Effects of rosemary extract and its residue on production, immune performance, and gut microbiota in geese

Introduction

To explore the effects of rosemary extract (RE) and its residue (RR) on the production, immune performance, and gut microbiota of geese.

Methods

We treat 28-day-old Sichuan white geese (n = 180) with three diets: (1) basal diet (control), (2) basal diet supplemented with 0.02% RE, and (3) basal diet supplemented with 15% RR for 42 days.

Results and discussion

On day 70, compared with control treatment, the final body weight, average daily gain and lysozyme levels in the RE treatment increased significantly (p < 0.05). In the RE and RR treatments, there was a significant decrease in alkaline phosphatase, globulin, and high-density lipoprotein levels compared to the control treatment, and there was also a significant increase in aspartate aminotransferase/alanine aminotransferase (p < 0.05). Moreover, for both RE and RR treatments, semi-eviscerated, eviscerated weights, and calcium apparent digestibility increased significantly, along with a decrease in the duodenal index (p < 0.05). Compared with RE treatment, those in the RR treatment had significantly higher duodenal and jejunum relative lengths, aspartate aminotransferase, uric acid, total cholesterol, and low-density lipoprotein levels, and decreased chest depth, chest angle, neck length, semi-eviscerated and eviscerated weights, crude protein digestibility, and levels of globulin, triglyceride, and lysozyme (p < 0.05). There were no differences in gut microbiota α or β diversities among treatments (p > 0.05). Compared to the control treatment, the relative abundance of Turicibacter significantly increased in the RR and RE treatments, and the relative abundance of Sporobacter, Alistipes, and Barnesiella significantly increased in the RR treatment (p < 0.05). Rikenellaceae, Succinivibrionaceae, and Aeromonadales were enriched in the RR treatment, and Lachnospiraceae, Turicibacteraceae, Fusobacteriaceae, and Enterobacteriaceae were enriched in the RE treatment. While we demonstrate the RR diet to be less effective than the RE diet, it did improve production and the gut microbiota of geese to a certain extent.

A novel codominant plumage color pattern of white breast patches in WugangTong geese was controlled by EDNRB2

Two basic plumage color patterns are observed in adult geese: solid grey (G) or colorless white (W). However, a Chinese indigenous breed, the Wugangtong goose (WGT), continues to be subject to selective breeding efforts as it displays segregation of plumage colors, including G, W, and a novel color pattern designated Wb (G with white breast circles). The underlying genetic mechanisms responsible for the Wb phenotype are yet to be determined. The current study employed the population differentiation index (FST) to analyze 90 geese exhibiting diverse plumage colors, identifying the fifth intron of EDNRB2 as a particularly noteworthy region with the highest FST values. Sanger sequencing of the region surrounding the EDNRB2 gene identified a 14-bp insertion within exon 3 as the causal mutation. The heterozygosity of this 14-bp insertion and wild-type alleles was completely associated with the Wb phenotype, thereby substantiating the codominant nature of the G and W phenotypes. An inter-species corroborated this finding cross between the graylag (no 14-bp insertion) and the swan goose (homozygous for the 14-bp insertion) breeds, as hybrids from this cross exhibited the Wb phenotype. Transcriptomes from white breast patches and gray dorsal skins of 4 Wb geese were compared. A significant downregulation of genes involved in melanin synthesis and melanocyte development was observed, including EDRNB2 and MLANA. The downregulation of MLANA indicated that the mutated EDNRB2 resulted in melanocyte loss in specific body regions, as MLANA is a marker gene for melanocytes. The findings were corroborated by melanin staining using the Mansson-Fontana method, which revealed no melanin particles deposited in the white breast patches. In summary, the gray plumage color was codominant to the white color in WGT geese, and plumage color variations were controlled by EDNRB2. The findings of our study offer valuable and practical guidance for the purification of plumage colors among WGT, whether through traditional phenotype selection or molecular breeding methods.

Chromosome-level genome provides insights into evolution and diving adaptability in the vulnerable common pochard (Aythya ferina)

The common pochard (Aythya ferina) is a freshwater diving duck found in the Palearctic region that has been classified as vulnerable by the IUCN due to continuous and rapid population declines across their distribution. To gain a better understanding of its genetic mechanism of adaptive evolution, we successfully sequenced and assembled the first high-quality chromosome-level genome of A. ferina using Illumina, Nanopore and Hi-C sequencing technologies. A total assembly length of 1,130.78 Mbp was obtained, with over 98.81% (1,117.37Mbp) of sequence anchored to 35 pseudo-chromosomes. We predicted 17,232 protein-coding genes, 95.9% of which were functionally annotated. We identified 339 expanded and 937 contracted gene families in the genome of A. ferina, and detected 95 genes that have been positively selected. The significantly enriched Gene Ontology and enriched pathways were related to energy metabolism, immune, nervous, and sensory systems, suggests that these factors likely played an important role in its evolution. Importantly, we recovered signatures of positive selection on genes related to vasoconstriction that may be associated with thermoregulatory adaptations of A. ferina for underwater diving. Overall, the high-quality genome assembly and annotation in this study provides valuable genomic resources for ecological and evolutionary studies, as well as toward the conservation of A. ferina.

Chromosome-level genome sequencing and multi-omics of the Hungarian White Goose (Anser anser domesticus) reveals novel miRNA-mRNA regulation mechanism of waterfowl feather follicle development

The Hungarian White Goose (Anser anser domesticus) is an excellent European goose breed, with high feather and meat production. Despite its importance in the poultry industry, no available genome assembly information has been published. This study aimed to present Chromosome-level and functional genome sequencing of the Hungarian White Goose. The results showed that the genome assembly has a total length of 1115.82 Mb, 39 pairs of chromosomes, 92.98% of the BUSCO index, and contig N50 and scaffold N50 were up to 2.32 Mb and 60.69 Mb, respectively. Annotation of the genome assembly revealed 19550 genes, 286 miRNAs, etc. We identified 235 expanded and 1,167 contracted gene families in this breed compared with the other 16 species. We performed a positive selection analysis between this breed and four species of Anatidae to uncover the genetic information underlying feather follicle development. Further, we detected the function of miR-199-x, miR-143-y, and miR-23-z on goose embryonic skin fibroblast. In summary, we have successfully generated a highly complete genome sequence of the Hungarian white goose, which will provide a great resource to improve our understanding of gene functions and enhance the studies on feather follicle development at the genomic level.

Disposition kinetics and tissue residues of tilmicosin following intravenous, subcutaneous, single and multiple oral dosing in geese (Anser anser domesticus)

Tilmicosin (TMC), a semi-synthetic macrolide antibiotic, is widely used in veterinary medicine due to its broad-spectrum, bacteriostatic properties. Frequently administered in various birds species, it is likely used off-label in geese as well. The study sought to investigate TMC's pharmacokinetics, tissue residues, in geese through in vivo experiments. The study involved longitudinal open studies on 15 healthy adult males, with three phases separated by one-month washout periods. Geese were administered TMC through intravenous (IV, 5 mg/kg), subcutaneous (SC, 10 mg/kg), and oral (PO, 25 mg/kg for five consecutive days) routes, with blood samples drawn at specific intervals. Tissue samples were also collected for subsequent analysis at pre-assigned times. TMC in goose plasma was quantified by a fully validated HPLC method. Plasma concentrations were quantified up to 4 hr for the PO and IV routes, and up to 10 hr in the SC route. Significant variations in bioavailability were observed between SC (87%) and PO (4%) routes. The body extraction ratio was low at 0.03, suggesting minimal ability of the liver and kidneys to eliminate TMC. Multiple oral doses showed no plasma accumulation, but tissue data revealed extensive distribution and prolonged residence, up to 120 h, suggesting a sustained therapeutic effect despite the brief plasma half-life. Regarding the multiple PO doses, provisional withdrawal times of 6, 7.5, and 8 days were suggested for the liver, muscles, and kidneys, respectively, according to the MRL set for these matrices in chickens by EMA. In conclusion, while the practical oral administration is discouraged at the population level, SC administration of TMC may be appropriate for geese, albeit impractical for flock therapy.

Comparative transcriptomic analyses reveal key genes underlying melanin distribution during embryonic development in geese (Anser anser)

1. Melanin distribution typically exhibits a gradient dilution along the dorsal-ventral axis of the body, including in domestic geese. However, the specific genes and molecular mechanisms responsible for this melanin distribution pattern remain incompletely understood.2. The transcriptomic comparisons were conducted at three embryonic stages, specifically on embryonic d 15 (E15), 22 (E22), and 29 (E29), between the pigmented dorsal skin and the depigmented distal foot.3. Differentially expressed genes (DEGs) associated with melanin synthesis were identified, particularly TYR, TYRP1, and EDNRB2, which exhibited significantly higher expression levels in the dorsal skin at E15 and E22. However, expression levels significantly decreased in later stages (E29).4. The ASIP gene showed remarkably high-expression levels in the distal feet compared to the dorsal skin post-E22 stage (log2FC: 5.31/6.88 at E22/E29). Gene Ontology (GO) enrichment analysis detected eight terms associated with melanin synthesis and melanosome formation (p < 0.05), including melanosome membrane (GO: 0033162) and melanin biosynthetic process (GO: 0042438). Additionally, KEGG pathway analysis showed significant enrichment of the melanogenesis pathway (hsa004916) at d 22 (E22).

Identifying candidate genes and biological pathways in muscle development through multi-tissue transcriptome comparisons between male and female geese

Males and females have long shown disparities in body weight and height; yet, the underlying mechanisms influencing growth and development remain unclear. Male and female Zhedong White Geese (ZDW) geese have long been selected for large body size and egg production, respectively. This led to a large difference in body weight between males and females, making them a unique model for studying the effects of sex on growth and development. This study aimed to elucidate these mechanisms by comparing the transcriptomes of muscle and pituitary tissues in male and female ZDW geese to identify the critical genes responsible for the effects of sex on growth performance. Our analysis revealed 1101 differentially expressed genes (DEGs) in leg musculature (507 upregulated, 594 downregulated), 773 DEGs in breast musculature (311 upregulated, 462 downregulated), and 517 DEGs in the pituitary gland (281 upregulated, 236 downregulated) between male and female geese. These DEGs were significantly enriched in gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with endocrine metabolism (e.g., hormonal activities), muscle formation (e.g., sarcomere and myofibril), and bone formation (e.g., bone morphogenesis and cartilage formation). The upregulated genes in males were enriched in KEGG pathways involving nutrient digestion and absorption (vitamin and protein), as well as the secretion of digestive juices (gastric acid and bile). Through protein-protein interaction analyses, we also observed high-density gene networks related to muscle fiber development, calcium ion metabolism, mitochondrial respiratory chain, and bone development. Therefore, our multi-tissue transcriptome analysis provides a deeper understanding of the complex and systematic gender-driven effects on growth and development in geese. IGF1, GHRHR, and NCAPG-LCORL and pathways related to myogenesis might play vital roles in gender differences before hormones exert their effect.

Telomere-to-telomere genome assembly of the goose Anser cygnoides

Our study presents the assembly of a high-quality Taihu goose genome at the Telomere-to-Telomere (T2T) level. By employing advanced sequencing technologies, including Pacific Biosciences HiFi reads, Oxford Nanopore long reads, Illumina short reads, and chromatin conformation capture (Hi-C), we achieved an exceptional assembly. The T2T assembly encompasses a total length of 1,197,991,206 bp, with contigs N50 reaching 33,928,929 bp and scaffold N50 attaining 81,007,908 bp. It consists of 73 scaffolds, including 38 autosomes and one pair of Z/W sex chromosomes. Importantly, 33 autosomes were assembled without any gap, resulting in a contiguous representation. Furthermore, gene annotation efforts identified 34,898 genes, including 436,162 RNA transcripts, encompassing 806,158 exons, 743,910 introns, 651,148 coding sequences (CDS), and 135,622 untranslated regions (UTR). The T2T-level chromosome-scale goose genome assembly provides a vital foundation for future genetic improvement and understanding the genetic mechanisms underlying important traits in geese.

Exploring the dynamic three-dimensional chromatin architecture and transcriptional landscape in goose liver tissues underlying metabolic adaptations induced by a high-fat diet

BACKGROUND

Goose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended.

RESULTS

In this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake.

CONCLUSIONS

We examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.

Comparative analysis of the intestinal microbiota of black-necked cranes (Grus nigricollis) in different wintering areas

Fecal microbiota is essential for host health because it increases digestive effectiveness. The crane species Grus nigricollis (G. nigricollis) is considered to be near threatened. The fecal microbial composition of crane is less understood, particularly in the Tibet, China. This study was performed to investigate the differences in fecal microbial composition and diversity of crane in different wintering areas using third-generation single-molecule real-time sequencing technology in the Tibet, China. According to the findings, 20 samples were used to generate 936 bacterial amplicon sequence variants (ASVs) and 1,800 fungal ASVs, only 4 bacterial ASVs and 20 fungal ASVs were shared in four distinct locations. Firmicutes were the dominant bacterial phylum in all samples, and Ascomycota and Basidiomycota were the dominant fungal phylum. At the genus level, Lactobacillus was the dominant genus in Linzhi City (LZ), Shannan City (SN), and Lasa City (LS), whereas Megamonas was the dominant genus in Rikaze City (RKZ). Naganishia and Mycosphaerella were the dominant fungal genera in SN and RKZ. Mycosphaerella and Tausonia were the dominant fungal genera in LZ. Naganishia and Fusarium were the dominant fungal genera in LS. And the fecal microbial composition varied between the four groups, as shown by the underweighted pair-group method with arithmetic means and principal coordinates analysis. This study offers a theoretical basis for understanding the fecal microbial composition of crane.

Insights into genetic diversity and phenotypic variations in domestic geese through comprehensive population and pan-genome analysis

BACKGROUND

Domestic goose breeds are descended from either the Swan goose (Anser cygnoides) or the Greylag goose (Anser anser), exhibiting variations in body size, reproductive performance, egg production, feather color, and other phenotypic traits. Constructing a pan-genome facilitates a thorough identification of genetic variations, thereby deepening our comprehension of the molecular mechanisms underlying genetic diversity and phenotypic variability.

RESULTS

To comprehensively facilitate population genomic and pan-genomic analyses in geese, we embarked on the task of 659 geese whole genome resequencing data and compiling a database of 155 RNA-seq samples. By constructing the pan-genome for geese, we generated non-reference contigs totaling 612 Mb, unveiling a collection of 2,813 novel genes and pinpointing 15,567 core genes, 1,324 softcore genes, 2,734 shell genes, and 878 cloud genes in goose genomes. Furthermore, we detected an 81.97 Mb genomic region showing signs of genome selection, encompassing the TGFBR2 gene correlated with variations in body weight among geese. Genome-wide association studies utilizing single nucleotide polymorphisms (SNPs) and presence-absence variation revealed significant genomic associations with various goose meat quality, reproductive, and body composition traits. For instance, a gene encoding the SVEP1 protein was linked to carcass oblique length, and a distinct gene-CDS haplotype of the SVEP1 gene exhibited an association with carcass oblique length. Notably, the pan-genome analysis revealed enrichment of variable genes in the "hair follicle maturation" Gene Ontology term, potentially linked to the selection of feather-related traits in geese. A gene presence-absence variation analysis suggested a reduced frequency of genes associated with "regulation of heart contraction" in domesticated geese compared to their wild counterparts. Our study provided novel insights into gene expression features and functions by integrating gene expression patterns across multiple organs and tissues in geese and analyzing population variation.

CONCLUSION

This accomplishment originates from the discernment of a multitude of selection signals and candidate genes associated with a wide array of traits, thereby markedly enhancing our understanding of the processes underlying domestication and breeding in geese. Moreover, assembling the pan-genome for geese has yielded a comprehensive apprehension of the goose genome, establishing it as an indispensable asset poised to offer innovative viewpoints and make substantial contributions to future geese breeding initiatives.

Metagenomic reconstructions of caecal microbiome in Landes, Roman and Zhedong White geese

1. The caecal microbiota in geese play a crucial role in determining the host's health, disease status and behaviour, as evidenced by extensive epidemiological data. The present investigation conducted 10× metagenomic sequencing of caecal content samples obtained from three distinct goose species, namely Landes geese, Roman geese and Zhedong White geese (n = 5), to explore the contribution of the gut microbiome to carbohydrate metabolism.2. In total, 337GB of Illumina data were generated, which identified 1,048,575 complete genes and construction of 331 metagenomic bins, encompassing 78 species from nine phyla. Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria and Bacteria were identified as the dominant phyla while Prevotella, Bacteroides, Streptococcus, and Subdoligranulum were the most abundant genera in the caecum of geese.3. The genes were allocated to 375 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) analysis. The most abundant classes in the caecum of geese were confirmed to be glycoside hydrolases (GHs), glycosyl transferases (GTs), as identified through the carbohydrate-active enzyme (CAZyme) database mapping. Subdoligranulum variabile and Mediterraneibacter glycyrrhizinilyticus were discovered to potentially facilitate carbohydrate digestion in geese.4. Notwithstanding, further investigation and validation are required to establish a connection between these species and CAZymes. Based on binning analysis, Mediterraneibacter glycyrrhizinilyticus and Ruminococcus sp. CAG:177 are potential species in LD geese that contribute to the production of fatty liver.

Chromosomal-level reference genome of a wild North American mallard (Anas platyrhynchos)

The mallard (Anas platyrhynchos) is one of the most common, economically, and socially important birds around the world. Mallards were not only an important food source for early humans but eventually becoming intimately linked with people as they were domesticated over the last 2,000 years. To date, mallard genomes are largely reconstructed from samples of domestic or unknown genetic heritage. Here, we report the first high-quality genome assembly and annotation of a genetically vetted wild mallard from North America (NAwild_v1.0). The genome was assembled using a combination of shotgun libraries, proximity ligation Chicago, and Dovetail Hi-C libraries. The final assembly is ∼1.04 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome-level scaffolds, and with a N50/L50 of 79.1 Mb/4 scaffolds. We used a combination of gene prediction and similarity approaches to annotate a total of 23,584 functional genes, of which 19,242 were associated to GO terms. The genome assembly and the set of annotated genes yielded a 95.4% completeness score when compared with the BUSCO aves_odb10 dataset. Next, we aligned 3 previously published mallard genomes to ours, and demonstrate how runs of homozygosity and nucleotide diversity are substantially higher and lower, respectively, to ours and how these artificially changed genomes resulted in profoundly different and biased demographic histories. Our wild mallard assembly not only provides a valuable resource to shed light onto genome evolution, speciation, and other adaptive processes, but also helping with identifying functional genes that have been significantly altered during the domestication process.

Dietary and Sexual Correlates of Gut Microbiota in the Japanese Gecko, Gekko japonicus (Schlegel, 1836)

Numerous studies have demonstrated that multiple intrinsic and extrinsic factors shape the structure and composition of gut microbiota in a host. The disorder of the gut microbiota may trigger various host diseases. Here, we collected fecal samples from wild-caught Japanese geckos (Gekko japonicus) and captive conspecifics fed with mealworms (mealworm-fed geckos) and fruit flies (fly-fed geckos), aiming to examine the dietary and sexual correlates of the gut microbiota. We used 16S rRNA gene sequencing technology to determine the composition of the gut microbiota. The dominant phyla with a mean relative abundance higher than 10% were Verrucomicrobiota, Bacteroidota, and Firmicutes. Gut microbial community richness and diversity were higher in mealworm-fed geckos than in wild geckos. Neither community evenness nor beta diversity of gut microbiota differed among wild, mealworm-fed, and fly-fed geckos. The beta rather than alpha diversity of gut microbiota was sex dependent. Based on the relative abundance of gut bacteria and their gene functions, we concluded that gut microbiota contributed more significantly to the host's metabolic and immune functions. A higher diversity of gut microbiota in mealworm-fed geckos could result from higher chitin content in insects of the order Coleoptera. This study not only provides basic information about the gut microbiota of G. japonicus but also shows that gut microbiota correlates with dietary habits and sex in the species.

Genome-resolving metagenomics reveals wild western capercaillies (Tetrao urogallus) as avian hosts for antibiotic-resistance bacteria and their interactions with the gut-virome community

The gut microbiome is a critical component of avian health, influencing nutrient uptake and immune functions. While the gut microbiomes of agriculturally important birds have been studied, the microbiomes of wild birds still need to be explored. Filling this knowledge gap could have implications for the microbial rewilding of captive birds and managing avian hosts for antibiotic-resistant bacteria (ARB). Using genome-resolved metagenomics, we recovered 112 metagenome-assembled genomes (MAGs) from the faeces of wild and captive western capercaillies (Tetrao urogallus) (n = 8). Comparisons of bacterial diversity between the wild and captive capercaillies suggest that the reduced diversity in the captive individual could be due to differences in diet. This was further substantiated through the analyses of 517,657 clusters of orthologous groups (COGs), which revealed that gene functions related to amino acids and carbohydrate metabolisms were more abundant in wild capercaillies. Metagenomics mining of resistome identified 751 antibiotic resistance genes (ARGs), of which 40.7 % were specific to wild capercaillies suggesting that capercaillies could be potential reservoirs for hosting ARG-associated bacteria. Additionally, the core resistome shared between wild and captive capercaillies indicates that birds can acquire these ARG-associated bacteria naturally from the environment (43.1 % of ARGs). The association of 26 MAGs with 120 ARGs and 378 virus operational taxonomic units (vOTUs) also suggests a possible interplay between these elements, where putative phages could have roles in modulating the gut microbiota of avian hosts. These findings can have important implications for conservation and human health, such as avian gut microbiota rewilding, identifying the emerging threats or opportunities due to phage-microbe interactions, and monitoring the potential spread of ARG-associated bacteria from wild avian populations.

Bacterial lipopolysaccharide with different administration routes affects intestinal mucosal morphological, immunological, and microbial barrier functions in goslings

The current study was conducted to evaluate the effects of different administration routes of bacterial lipopolysaccharide (LPS) on intestinal mucosal morphological, immunological, and microbial barrier functions in goslings. First, we compared intestinal villi morphology of goslings under intraperitoneal or oral LPS treatment through hematoxylin and eosin staining. Then, we determined the signatures of the microbiome in the ileum mucosa of goslings subjected to oral LPS treatment at 0, 2, 4, and 8 mg/kg BW by 16S sequencing, and analyzed the changes in intestinal barrier functions and permeability, levels of LPS in the ileum mucosa, plasma, and liver tissue, and the induced inflammatory response of Toll-like receptor 4 (TLR4). As a result, intraperitoneal LPS injection resulted in a thicker intestinal wall in the ileum within a short time, whereas villus height was less affected; in contrast, oral LPS treatment exerted a stronger influence on villus height but not on intestinal wall thickness. We also found that oral LPS treatment affected the structure of the intestinal microbiome, reflected by changes in the clustering of intestinal microbiota. The average abundance of Muribaculaceae showed an increasing trend with increasing LPS levels, and that of the genus Bacteroides decreased, compared with the control group. In addition, oral LPS treatment with 8 mg/kg BW affected the intestinal epithelial morphology, damage the mucosal immune barrier, downregulated the expression of tight junction proteins, increased circulating D-lactate levels, and stimulated the secretion of various inflammatory mediators and activation of the TLR4/MyD88/NFκB pathway. This study presented the injuries of intestinal mucosal barrier function induced by LPS challenges in goslings and provided a scientific model for searching the novel strategies to attenuate the immunological stress and gut injury caused by LPS.

The integrated analysis of digestive physiology and gastrointestinal microbiota structure in Changle goose

Changle goose in Fujian, China is a rare genetic resource and in urgent need to be protected. Understanding the characteristics of digestive physiology and spatial variation of gastrointestinal microbiota is crucial for developing nutritional intervention strategies to improve intestinal health and production performance of goose. Hence, histomorphological assay was used for observing development status of proventriculus, jejunum, and cecum in 70-day-old Changle geese, whereas digesta from 6 alimentary canal locations (crop, proventriculus, gizzard, jejunum, cecum, and rectum) were collected for 16S rRNA gene sequencing and short chain fatty acids (SCFAs) quantitative analysis. The histomorphological observation indicated that the jejunum and cecum of Changle goose were well developed. The alpha diversity analysis revealed that, except rectum, microbiota in other noncecum sections were in high diversity as cecum. The Nonmetric MultiDimensional Scaling (NMDS) analysis showed that microbial community of proventriculus, gizzard, and jejunum formed a cluster, which distinctly discrete with the microbiota of the other gastrointestinal locations. Additionally, the proportions of Proteobacteria, Bacteroidota, and Campilobacterota at the phylum level and Lactobacillus, Streptococcus, Helicobacter, and Subdoligranulum at the genus level exhibited tremendous alternations among different gastrointestinal locations. The characteristic bacterial composition in each section was further disclosed by analyzing the core and feature Amplicon Sequence Variants (ASVs) and SCFAs pattern. Importantly, 7 body-weight-associated ASVs and 2 cecum-development-related ASVs were identified via correlation analysis. In a whole, our findings provided the first insights into the specialized digestive physiology of Changle geese and distinctive regional distribution of gastrointestinal microbiota, which laid the important foundation for improving growth performance through microbiota manipulation in geese.

Origins, timing and introgression of domestic geese revealed by whole genome data

BACKGROUND

Geese are among the most important poultry species in the world. The current generally accepted hypothesis is that the European domestic geese originated from greylag geese (Anser anser), and Chinese domestic geese have two origins, most of which originated from swan geese (Anser cygnoides), and the Yili goose originated from greylag geese. To explain the origin and demographic history of geese, we selected 14 goose breeds from Europe and China and wild populations of swan and greylag geese, and whole genome sequencing data were obtained for 74 samples.

RESULTS

Population structure analysis and phylogenetic trees showed that the wild ancestor of Chinese domestic geese, except for Yili, is the swan geese, and the wild ancestor of Chinese Yili and European domestic geese is greylag geese. Analysis of the demographic history suggests that the domestication of Chinese geese occurred ~ 3499 years ago and that of the European geese occurred ~ 7552 years ago. Furthermore, gene flow was observed between domestic geese and their wild ancestors. Analysis of introgression showed that Yili geese had been introgressed by Chinese domestic geese, and the body size of Yili geese may be influenced by introgression events of some growth-related genes, including IGF-1.

CONCLUSIONS

Our study provides evidence for the origin of geese at the genome-wide level and advances the understanding of the history of goose domestication and the traits affected by introgression events.

Highly differentiated loci resolve phylogenetic relationships in the Bean Goose complex

BACKGROUND

Reconstructing phylogenetic relationships with genomic data remains a challenging endeavor. Numerous phylogenomic studies have reported incongruent gene trees when analyzing different genomic regions, complicating the search for a 'true' species tree. Some authors have argued that genomic regions of increased divergence (i.e. differentiation islands) reflect the species tree, although other studies have shown that these regions might produce misleading topologies due to species-specific selective sweeps or ancient introgression events. In this study, we tested the extent to which highly differentiated loci can resolve phylogenetic relationships in the Bean Goose complex, a group of goose taxa that includes the Taiga Bean Goose (Anser fabalis), the Tundra Bean Goose (Anser serrirostris) and the Pink-footed Goose (Anser brachyrhynchus).

RESULTS

First, we show that a random selection of genomic loci-which mainly samples the undifferentiated regions of the genome-results in an unresolved species complex with a monophyletic A. brachyrhynchus embedded within a paraphyletic cluster of A. fabalis and A. serrirostris. Next, phylogenetic analyses of differentiation islands converged upon a topology of three monophyletic clades in which A. brachyrhynchus is sister to A. fabalis, and A. serrirostris is sister to the clade uniting these two species. Close inspection of the locus trees within the differentiated regions revealed that this topology was consistently supported over other phylogenetic arrangements. As it seems unlikely that selection or introgression events have impacted all differentiation islands in the same way, we are convinced that this topology reflects the 'true' species tree. Additional analyses, based on D-statistics, revealed extensive introgression between A. fabalis and A. serrirostris, which partly explains the failure to resolve the species complex with a random selection of genomic loci. Recent introgression between these taxa has probably erased the phylogenetic branching pattern across a large section of the genome, whereas differentiation islands were unaffected by the homogenizing gene flow and maintained the phylogenetic patterns that reflect the species tree.

CONCLUSIONS

The evolution of the Bean Goose complex can be depicted as a simple bifurcating tree, but this would ignore the impact of introgressive hybridization. Hence, we advocate that the evolutionary relationships between these taxa are best represented as a phylogenetic network.

Chromosome-level genome assembly of goose provides insight into the adaptation and growth of local goose breeds

BACKGROUND

Anatidae contains numerous waterfowl species with great economic value, but the genetic diversity basis remains insufficiently investigated. Here, we report a chromosome-level genome assembly of Lion-head goose (Anser cygnoides), a native breed in South China, through the combination of PacBio, Bionano, and Hi-C technologies.

FINDINGS

The assembly had a total genome size of 1.19 Gb, consisting of 1,859 contigs with an N50 length of 20.59 Mb, generating 40 pseudochromosomes, representing 97.27% of the assembled genome, and identifying 21,208 protein-coding genes. Comparative genomic analysis revealed that geese and ducks diverged approximately 28.42 million years ago, and geese have undergone massive gene family expansion and contraction. To identify genetic markers associated with body weight in different geese breeds, including Wuzong goose, Huangzong goose, Magang goose, and Lion-head goose, a genome-wide association study was performed, yielding an average of 1,520.6 Mb of raw data that detected 44,858 single-mucleotide polymorphisms (SNPs). Genome-wide association study showed that 6 SNPs were significantly associated with body weight and 25 were potentially associated. The significantly associated SNPs were annotated as LDLRAD4, GPR180, and OR, enriching in growth factor receptor regulation pathways.

CONCLUSIONS

We present the first chromosome-level assembly of the Lion-head goose genome, which will expand the genomic resources of the Anatidae family, providing a basis for adaptation and evolution. Candidate genes significantly associated with different goose breeds may serve to understand the underlying mechanisms of weight differences.

Using comparative genomics to detect mutations regulating plumage variations in graylag (A. anser) and swan geese (A. cygnoides)

Although graylag geese (A. anser) showed similar plumages of white, grey, and white with grey patches compared to those in swan geese (A. cygnoides), it was believed the substantial molecular mechanism for plumage variations were different. To date, studies on genes responsible for diverse plumages among graylag geese were limited and causal mutations remain unknown. In this study, genomes from 57 individuals belonging to six breeds showing different plumages were sequenced at ∼10X depth. Firstly, the allele frequency differences (AFD) of variants on the scaffold394 (NW_013185915.1) between grey and white goose breeds (A. anser) was calculated and a genomic region between 768,290-779,889 bp was detected to carry candidate variants associated with plumages, including one SNP (g. 775,151G > T, ∼18.6 kb upstream of EDNRB2) found to be fixed in white geese. This region was overlapped with the one detected by the haplotype-based sweep analysis, in which significant signals defined a candidate region of 736,610-820,622 bp on the same scaffold. Results from the transcriptomic data showed that expression levels of EDNRB2 and many other melanogenesis-related genes were significantly decreased among white geese compared to that in grey geese, especially at late embryonic stages (>E15). Modifications at transcriptional levels might result in abnormal melanocyte developments and thus the white plumages when they grow up. In addition, a frameshift mutation (C > -) in exon4 of MLANA gene on scaffold176 (NW_013185876.1) was suggested as the causal mutation for sex-linked dilution phenotype in graylag geese although this requires more demonstration experiments. Together with observed white plumages caused by EDNRB2 mutations in coding regions among swan geese and chicken, our study provided new examples to study the parallel evolution.

Effects of Yeast Culture Supplementation on Growth Performance, Nutrient Digestibility, Blood Metabolites, and Immune Response in Geese

The present study was conducted to investigate the effects of dietary yeast culture (YC) supplementation on growth performance, nutrient digestibility, blood metabolites, and immune functions in geese. One-day-old Sichuan white geese (n = 300) were randomly divided into five groups containing 0 (control), 0.5%, 1.0%, 2.0%, and 4.0% of YC in the diet for 70 days. In general, the dietary supplementation of YC significantly increased the average daily gain and feed conversion ratio (p < 0.05) in which the 1.0% or 2.0% levels were better and significantly reduced the average daily feed intake at the 2.0% level (p < 0.05). YC supplementation increased digestibility of P (quadratic, p = 0.01) and gross energy (quadratic, p = 0.04) from days 23 to 27 and crude protein from days 23 to 27 and days 64 to 68 (quadratic, p ≤ 0.05), with the 2.0% level being the most effective. Serum metabolites were significantly affected by dietary YC (p < 0.05). Supplemental YC increased IL-2 on day 28 (linear, p = 0.01; quadratic, p = 0.04) and lysozyme on day 70 (quadratic, p = 0.04) and decreased complement C4 on day 70 (linear, p = 0.05). Interferon-γ, interleukin-2, and tumor necrosis factor-α genes were mostly up-regulated after YC supplementation, and interferon-γ and interleukin-2 gene expression levels were significantly increased at the 2.0% level (p < 0.05). Taken together, dietary YC supplementation improved growth performance and affected nutrient digestibility, serum metabolites, and immune function in geese, which was optimized at the 2% YC level in the present study.

Environmental fluoride exposure disrupts the intestinal structure and gut microbial composition in ducks

Fluorine (F) and its compounds produced from industrial production and coal combustion can cause air, water and soil contamination, which can accumulate in animals, plants and humans via food chain threatening public health. Fluoride exposure affects liver, kidney, gastrointestinal and reproductive system in humans and animals. Literature regarding fluoride influence on intestinal structure and microbiota composition in ducks is scarce. This study was designed to investigate these effects by using simple and electron microscopy and 16S rRNA sequencing techniques. Results indicated an impaired structure with reduced relative distribution of goblet cells in the fluoride exposed group. Moreover, the gut microbiota showed a significant decrease in alpha diversity. Proteobacteria, Firmicutes and Bacteroidetes were the most abundant phyla in both control and fluoride-exposed groups. Specifically, fluoride exposure resulted in a significant decrease in the relative abundance of 9 bacterial phyla and 15 bacterial genera. Among them, 4 phyla (Latescibacteria, Dependentiae, Zixibacteria and Fibrobacteres) and 4 genera (Thauera, Hydrogenophaga, Reyranella and Arenimonas) weren't even detectable in the gut microbiota of the ducks. In summary, higher fluoride exposure can significantly damage the intestinal structure and gut microbial composition in ducks.

Pooled Sequencing Analysis of Geese (Anser cygnoides) Reveals Genomic Variations Associated With Feather Color

During the domestication of the goose a change in its feather color took place, however, the molecular mechanisms responsible for this change are not completely understood. Here, we performed whole-genome resequencing on three pooled samples of geese (feral and domestic geese), with two distinct feather colors, to identify genes that might regulate feather color. We identified around 8 million SNPs within each of the three pools and validated allele frequencies for a subset of these SNPs using PCR and Sanger sequencing. Several genomic regions with signatures of differential selection were found when we compared the gray and white feather color populations using the F_ST and Hp approaches. When we combined previous functional studies with our genomic analyses we identified 26 genes (KITLG, MITF, TYRO3, KIT, AP3B1, SMARCA2, ROR2, CSNK1G3, CCDC112, VAMP7, SLC16A2, LOC106047519, RLIM, KIAA2022, ST8SIA4, LOC106044163, TRPM6, TICAM2, LOC106038556, LOC106038575, LOC106038574, LOC106038594, LOC106038573, LOC106038604, LOC106047489, and LOC106047492) that potentially regulate feather color in geese. These results substantially expand the catalog of potential feather color regulators in geese and provide a basis for further studies on domestication and avian feather coloration.

Domestication Shapes the Community Structure and Functional Metagenomic Content of the Yak Fecal Microbiota

Domestication is a key factor of genetic variation; however, the mechanism by which domestication alters gut microbiota is poorly understood. Here, to explore the variation in the structure, function, rapidly evolved genes (REGs), and enzyme profiles of cellulase and hemicellulose in fecal microbiota, we studied the fecal microbiota in wild, half-blood, and domestic yaks based on 16S rDNA sequencing, shotgun-metagenomic sequencing, and the measurement of short-chain-fatty-acids (SCFAs) concentration. Results indicated that wild and half-blood yaks harbored an increased abundance of the phylum Firmicutes and reduced abundance of the genus Akkermansia, which are both associated with efficient energy harvesting. The gut microbial diversity decreased in domestic yaks. The results of the shotgun-metagenomic sequencing showed that the wild yak harbored an increased abundance of microbial pathways that play crucial roles in digestion and growth of the host, whereas the domestic yak harbored an increased abundance of methane-metabolism-related pathways. Wild yaks had enriched amounts of REGs in energy and carbohydrate metabolism pathways, and possessed a significantly increased abundance of cellulases and endohemicellulases in the glycoside hydrolase family compared to domestic yaks. The concentrations of acetic, propionic, n-butyric, i-butyric, n-valeric, and i-valeric acid were highest in wild yaks. Our study displayed the domestic effect on the phenotype of composition, function in gut microbiota, and SCFAs associated with gut microbiota, which had a closely association with the growth performance of the livestock. These findings may enlighten the researchers to construct more links between economic characteristics and gut microbiota, and develop new commercial strains in livestock based on the biotechnology of gut microbiota.

Genome-Wide Association Study-Based Identification of SNPs and Haplotypes Associated With Goose Reproductive Performance and Egg Quality

Geese are one of the most economically important waterfowl. However, the low reproductive performance and egg quality of geese hinder the development of the goose industry. The identification and application of genetic markers may improve the accuracy of beneficial trait selection. To identify the genetic markers associated with goose reproductive performance and egg quality traits, we performed a genome-wide association study (GWAS) for body weight at birth (BBW), the number of eggs at 48 weeks of age (EN48), the number of eggs at 60 weeks of age (EN60) and egg yolk color (EYC). The GWAS acquired 2.896 Tb of raw sequencing data with an average depth of 12.44× and identified 9,279,339 SNPs. The results of GWAS showed that 26 SNPs were significantly associated with BBW, EN48, EN60, and EYC. Moreover, five of these SNPs significantly associated with EN48 and EN60 were in a haplotype block on chromosome 35 from 4,512,855 to 4,541,709 bp, oriented to TMEM161A and another five SNPs significantly correlated to EYC were constructed in haplotype block on chromosome 5 from 21,069,009 to 21,363,580, which annotated by TMEM161A, CALCR, TFPI2, and GLP1R. Those genes were enriched in epidermal growth factor-activated receptor activity, regulation of epidermal growth factor receptor signaling pathway. The SNPs, haplotype markers, and candidate genes identified in this study can be used to improve the accuracy of marker-assisted selection for the reproductive performance and egg quality traits of geese. In addition, the candidate genes significantly associated with these traits may provide a foundation for better understanding the mechanisms underlying reproduction and egg quality in geese.

The intestinal microbial composition in Greylag geese differs with steatosis induction mode: spontaneous or induced by overfeeding

Background

Relationships between microbial composition and steatosis are being extensively studied in mammals, and causal relations have been evidenced. In migratory birds the liver can transiently store lipids during pre-migratory and migratory phases, but little is known about the implications of the digestive microbiota in those mechanisms. The Landaise greylag goose (Anser anser) is a good model to study steatosis in migratory birds as it is domesticated, but is still, from a genetic point of view, close to its wild migratory ancestor. It also has a great ingestion capacity and a good predisposition for hepatic steatosis, whether spontaneous or induced by conventional overfeeding. The conventional (overfeeding) and alternative (spontaneous steatosis induction) systems differ considerably in duration and feed intake level and previous studies have shown that aptitudes to spontaneous steatosis are very variable. The present study thus aimed to address two issues: (i) evaluate whether microbial composition differs with steatosis-inducing mode; (ii) elucidate whether a digestive microbial signature could be associated with variable aptitudes to spontaneous liver steatosis.

Results

Performances, biochemical composition of the livers and microbiota differed considerably in response to steatosis stimulation. We namely identified the genus Romboutsia to be overrepresented in birds developing a spontaneous steatosis in comparison to those submitted to conventional overfeeding while the genera Ralstonia, Variovorax and Sphingomonas were underrepresented only in birds that did not develop a spontaneous steatosis compared to conventionally overfed ones, birds developing a spontaneous steatosis having intermediate values. Secondly, no overall differences in microbial composition were evidenced in association with variable aptitudes to spontaneous steatosis, although one OTU, belonging to the Lactobacillus genus, was overrepresented in birds having developed a spontaneous steatosis compared to those that had not.

Conclusions

Our study is the first to evaluate the intestinal microbial composition in association with steatosis, whether spontaneous or induced by overfeeding, in geese. Steatosis induction modes were associated with distinct digestive microbial compositions. However, unlike what can be observed in mammals, no clear microbial signature associated with spontaneous steatosis level was identified.

Pacific Biosciences assembly with Hi-C mapping generates an improved, chromosome-level goose genome

BACKGROUND

The domestic goose is an economically important and scientifically valuable waterfowl; however, a lack of high-quality genomic data has hindered research concerning its genome, genetics, and breeding. As domestic geese breeds derive from both the swan goose (Anser cygnoides) and the graylag goose (Anser anser), we selected a female Tianfu goose for genome sequencing. We generated a chromosome-level goose genome assembly by adopting a hybrid de novo assembly approach that combined Pacific Biosciences single-molecule real-time sequencing, high-throughput chromatin conformation capture mapping, and Illumina short-read sequencing.

FINDINGS

We generated a 1.11-Gb goose genome with contig and scaffold N50 values of 1.85 and 33.12 Mb, respectively. The assembly contains 39 pseudo-chromosomes (2n = 78) accounting for ∼88.36% of the goose genome. Compared with previous goose assemblies, our assembly has more continuity, completeness, and accuracy; the annotation of core eukaryotic genes and universal single-copy orthologs has also been improved. We have identified 17,568 protein-coding genes and a repeat content of 8.67% (96.57 Mb) in this genome assembly. We also explored the spatial organization of chromatin and gene expression in the goose liver tissues, in terms of inter-pseudo-chromosomal interaction patterns, compartments, topologically associating domains, and promoter-enhancer interactions.

CONCLUSIONS

We present the first chromosome-level assembly of the goose genome. This will be a valuable resource for future genetic and genomic studies on geese.

Recent introgression between Taiga Bean Goose and Tundra Bean Goose results in a largely homogeneous landscape of genetic differentiation

Several studies have uncovered a highly heterogeneous landscape of genetic differentiation across the genomes of closely related species. Specifically, genetic differentiation is often concentrated in particular genomic regions (“islands of differentiation”) that might contain barrier loci contributing to reproductive isolation, whereas the rest of the genome is homogenized by introgression. Alternatively, linked selection can produce differentiation islands in allopatry without introgression. We explored the influence of introgression on the landscape of genetic differentiation in two hybridizing goose taxa: the Taiga Bean Goose (Anser fabalis) and the Tundra Bean Goose (A. serrirostris). We re-sequenced the whole genomes of 18 individuals (9 of each taxon) and, using a combination of population genomic summary statistics and demographic modeling, we reconstructed the evolutionary history of these birds. Next, we quantified the impact of introgression on the build-up and maintenance of genetic differentiation. We found evidence for a scenario of allopatric divergence (about 2.5 million years ago) followed by recent secondary contact (about 60,000 years ago). Subsequent introgression events led to high levels of gene flow, mainly from the Tundra Bean Goose into the Taiga Bean Goose. This scenario resulted in a largely undifferentiated genomic landscape (genome-wide F_ST = 0.033) with a few notable differentiation peaks that were scattered across chromosomes. The summary statistics indicated that some peaks might contain barrier loci while others arose in allopatry through linked selection. Finally, based on the low genetic differentiation, considerable morphological variation and incomplete reproductive isolation, we argue that the Taiga and the Tundra Bean Goose should be treated as subspecies.

Molecular cloning of the goose GnRH gene and identification of GnRH polymorphisms associated with laying traits

1. Egg-laying traits are important economic characteristics in goose production (Anser cygnoides). The gene GnRH, which encodes gonadotropin-releasing hormone, is a strong candidate gene for egg-laying traits in avian species. 2. In this study, a 3520 bp genomic sequence and a 279 bp mRNA sequence for GnRH, which encoded 92 amino acids, were determined. The GnRH DNA sequence contains four exons and three introns, and the DNA and deduced amino acid sequences were highly conserved across mammals (human, macaque, cow, and sheep) and avians (chicken, fulmar and quail). 3. Using a direct sequencing method, 46 single nucleotide polymorphisms (SNPs) were identified in the GnRH genomic sequence that were shared between two Sichuan White goose populations (217 and 208 individuals). Furthermore, 44 haplotypes were constructed using a sliding window approach. Association analysis between the SNPs and haplotypes and egg-laying traits showed that 10 SNPs affected the first egg weight, average egg weight, egg number at 48 weeks and egg number at 64 weeks. 4. These results lay the foundation for further studies of the function of GnRH in geese and provide a theoretical basis for marker-assisted selection of egg-laying traits in the Sichuan white goose population.

First de novo whole genome sequencing and assembly of the bar-headed goose

Background

The bar-headed goose (Anser indicus) mainly inhabits the plateau wetlands of Asia. As a specialized high-altitude species, bar-headed geese can migrate between South and Central Asia and annually fly twice over the Himalayan mountains along the central Asian flyway. The physiological, biochemical and behavioral adaptations of bar-headed geese to high-altitude living and flying have raised much interest. However, to date, there is still no genome assembly information publicly available for bar-headed geese.

Methods

In this study, we present the first de novo whole genome sequencing and assembly of the bar-headed goose, along with gene prediction and annotation.

Results

10X Genomics sequencing produced a total of 124 Gb sequencing data, which can cover the estimated genome size of bar-headed goose for 103 times (average coverage). The genome assembly comprised 10,528 scaffolds, with a total length of 1.143 Gb and a scaffold N50 of 10.09 Mb. Annotation of the bar-headed goose genome assembly identified a total of 102 Mb (8.9%) of repetitive sequences, 16,428 protein-coding genes, and 282 tRNAs. In total, we determined that there were 63 expanded and 20 contracted gene families in the bar-headed goose compared with the other 15 vertebrates. We also performed a positive selection analysis between the bar-headed goose and the closely related low-altitude goose, swan goose (Anser cygnoides), to uncover its genetic adaptations to the Qinghai-Tibetan Plateau.

Conclusion

We reported the currently most complete genome sequence of the bar-headed goose. Our assembly will provide a valuable resource to enhance further studies of the gene functions of bar-headed goose. The data will also be valuable for facilitating studies of the evolution, population genetics and high-altitude adaptations of the bar-headed geese at the genomic level.

Gut microbiota correlates with fiber and apparent nutrients digestion in goose

To explore the relationship among the level of fiber, gut microbiota, and nutritional substances, we applied the next generation sequencing technology for the identification of the composition and structure of microbiota in the gastrointestinal tract. In this study, a total of 25 phyla and 298 genera were identified from the gastrointestinal tract; Firmicutes, Bacteroidetes, and Proteobacteria were the predominant phyla. The ability of cecum in carbohydrate metabolism was significantly higher than that of the gizzard and ileum (P < 0.05). The bacterial community structure in various stages of the development of the cecum was different. In the different growth stages of cecum, the increase in the microbiota structure of the fiber level elevates the ability of carbon hydration. Second, the apparent metabolic rates of the other nutrients were affected by the fiber and period except for acid detergent fiber (P < 0.05); the apparent utilization rate of the nutrients increased with time. However, with the increase in the fiber level, the apparent utilization of nutrients was initially increased, followed by a decrease. Therefore, a correlation was established between the fiber level and gastrointestinal microbiota and apparent nutrient utilization rate of the 3 phyla. Our results suggest that the fiber level and growth stages could impact the composition of gut microbiota.

Comparative proteomic analysis of hypothalamus tissue from Huoyan geese between pre-laying period and laying period using an iTRAQ-based approach

The hypothalamus plays a central role in controlling poultry endocrine and reproductive activities. So far there is limited information focused on the proteome profiles of the hypothalamus from geese during different stages of the egg-laying cycle. In order to identify proteins regulating the egg-laying process of Huoyan geese, we investigated the proteome profiles of the hypothalamus from Huoyan geese during the laying period and pre-laying period by applying an isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomic technology. A total number of 3,337 were identified and quantified, of which 18 were significantly up-regulated and 16 were significantly down-regulated. These differentially expressed proteins were subjected to bioinformatics analyses based on the Gene Ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathway. Some of these were revealed to be involved in hormone and neurotransmitter secretion, exocytosis, calcium ion transport and synaptic transmission. Subsequently, excitatory amino acid transporter 2, complexin-1 and inositol 1,4,5-trisphosphate receptor, type 3 were confirmed at the messenger RNA level using quantitative real-time RT-PCR. Then, the abundance change of these proteins was verified further using Western blotting analysis. These data may aid in elucidating the molecular mechanism of higher laying performance in Huoyan geese.

Spatial Organization of the Gastrointestinal Microbiota in Urban Canada Geese

Recent reviews identified the reliance on fecal or cloacal samples as a significant limitation hindering our understanding of the avian gastrointestinal (gut) microbiota and its function. We investigated the microbiota of the esophagus, duodenum, cecum, and colon of a wild urban population of Canada goose (Branta canadensis). From a population sample of 30 individuals, we sequenced the V4 region of the 16S SSU rRNA on an Illumina MiSeq and obtained 8,628,751 sequences with a median of 76,529 per sample. These sequences were assigned to 420 bacterial OTUs and a single archaeon. Firmicutes, Proteobacteria, and Bacteroidetes accounted for 90% of all sequences. Microbiotas from the four gut regions differed significantly in their richness, composition, and variability among individuals. Microbial communities of the esophagus were the most distinctive whereas those of the colon were the least distinctive, reflecting the physical downstream mixing of regional microbiotas. The downstream mixing of regional microbiotas was also responsible for the majority of observed co-occurrence patterns among microbial families. Our results indicate that fecal and cloacal samples inadequately represent the complex patterns of richness, composition, and variability of the gut microbiota and obscure patterns of co-occurrence of microbial lineages.

Differentiated surface fungal communities at point of harvest on apple fruits from rural and peri-urban orchards

The diverse fungal communities that colonize fruit surfaces are closely associated with fruit development, preservation and quality control. However, the overall fungi adhering to the fruit surface and the inference of environmental factors are still unknown. Here, we characterized the fungal signatures on apple surfaces by sequencing internal transcribed spacer 1 (ITS1) region. We collected the surface fungal communities from apple fruits cultivated in rural and peri-urban orchards. A total of 111 fungal genera belonging to 4 phyla were identified, showing remarkable fungal diversity on the apple surface. Comparative analysis of rural samples harboured higher fungal diversity than those from peri-urban orchards. In addition, fungal composition varied significantly across apple samples. At the genus level, the protective genera Coniothyrium, Paraphaeosphaeria and Periconia were enriched in rural samples. The pathogenic genera Acremonium, Aspergillus, Penicillium and Tilletiposis were enriched in peri-urban samples. Our findings indicate that rural samples maintained more diverse fungal communities on apple surfaces, whereas peri-urban-planted apple carried potential pathogenic risks. This study sheds light on ways to improve fruit cultivation and disease prevention practices.

Comparative proteomic analysis of pituitary glands from Huoyan geese between pre-laying and laying periods using an iTRAQ-based approach

In this study, we performed a comprehensive evaluation of the proteomic profile of the pituitary gland of the Huoyan goose during the laying period compared to the pre-laying period using an iTRAQ-based approach. Protein samples were prepared from pituitary gland tissues of nine pre-laying period and nine laying period geese. Then the protein samples from three randomly selected geese within each period were pooled in equal amounts to generate one biological sample pool. We identified 684 differentially expressed proteins, including 418 up-regulated and 266 down-regulated proteins. GO annotation and KEGG pathway analyses of these proteins were conducted. Some of these proteins were found to be associated with hormone and neurotransmitter secretion and transport, neuropeptide signalling and GnRH signalling pathways, among others. Subsequently, the modification of the abundance of three proteins (prolactin, chromogranin-A and ITPR3) was verified using western blotting. Our results will provide a new source for mining genes and gene products related to the egg-laying performance of Huoyan geese, and may provide important information for the conservation and utilization of local goose breeds.

Effects of One-Week Empirical Antibiotic Therapy on the Early Development of Gut Microbiota and Metabolites in Preterm Infants

The early postnatal period is the most dynamic and vulnerable stage in the assembly of intestinal microbiota. Antibiotics are commonly prescribed to newborn preterm babies and are frequently used for a prolonged duration in China. We hypothesized that the prolonged antibiotic therapy would affect the early development of intestinal microbiota and their metabolites. To test this hypothesis, we analyzed the stool microbiota and metabolites in 36 preterm babies with or without antibiotic treatment. These babies were divided into three groups, including two groups treated with the combination of penicillin and moxalactam or piperacillin-tazobactam for 7 days, and the other group was free of antibiotics. Compared to the antibiotic-free group, both antibiotic-treated groups had distinct gut microbial communities and metabolites, including a reduction of bacterial diversity and an enrichment of harmful bacteria such as Streptococcus and Pseudomonas. In addition, there was a significant difference in the composition of gut microbiota and their metabolites between the two antibiotic-treated groups, where the piperacillin-tazobactam treatment group showed an overgrowth of Enterococcus. These findings suggest that prolonged antibiotic therapy affects the early development of gut microbiota in preterm infants, which should be considered when prescribing antibiotics for this population.

Evolutionary Biology Needs Wild Microbiomes

The microbiome is a vital component to the evolution of a host and much of what we know about the microbiome derives from studies on humans and captive animals. But captivity alters the microbiome and mammals have unique biological adaptations that affect their microbiomes (e.g., milk). Birds represent over 30% of known tetrapod diversity and possess their own suite of adaptations relevant to the microbiome. In a previous study, we showed that 59 species of birds displayed immense variation in their microbiomes and host (bird) taxonomy and ecology were most correlated with the gut microbiome. In this Frontiers Focused Review, I put those results in a broader context by discussing how collecting and analyzing wild microbiomes contributes to the main goals of evolutionary biology and the specific ways that birds are unique microbial hosts. Finally, I outline some of the methodological considerations for adding microbiome sampling to the research of wild animals and urge researchers to do so. To truly understand the evolution of a host, we need to understand the millions of microorganisms that inhabit it as well: evolutionary biology needs wild microbiomes.

Metaproteomic strategies and applications for gut microbial research

The human intestine hosts various complex microbial communities that are closely associated with multiple health and disease processes. Determining the composition and function of these microbial communities is critical to unveil disease mechanisms and promote human health. Recently, meta-omic strategies have been developed that use high-throughput techniques to provide a wealth of information, thus accelerating the study of gut microbes. Metaproteomics is a newly emerged analytical approach that aims to identify proteins on a large scale in complex environmental microbial communities (e.g., the gut microbiota). This review introduces the recent analytical strategies and applications of metaproteomics, with a focus on advances in gut microbiota research, including a discussion of the limitations and challenges of these approaches.