Metagenomic Analysis Identifies Sex-Related Cecal Microbial Gene Functions and Bacterial Taxa in the Quail.

Background: Japanese quail (Coturnix japonica) are important and widely distributed poultry in China. Researchers continue to pursue genetic selection for heavier quail. The intestinal microbiota plays a substantial role in growth promotion; however, the mechanisms involved in growth promotion remain unclear. Results: We generated 107.3 Gb of cecal microbiome data from ten Japanese quail, providing a series of quail gut microbial gene catalogs (1.25 million genes). We identified a total of 606 main microbial species from 1,033,311 annotated genes distributed among the ten quail. Seventeen microbial species from the genera Anaerobiospirillum, Alistipes, Barnesiella, and Butyricimonas differed significantly in their abundances between the female and male gut microbiotas. Most of the functional gut microbial genes were involved in metabolism, primarily in carbohydrate transport and metabolism, as well as some active carbohydrate-degrading enzymes. We also identified 308 antibiotic-resistance genes (ARGs) from the phyla Bacteroidetes, Firmicutes and Euryarchaeota. Studies of the differential gene functions between sexes indicated that abundances of the gut microbes that produce carbohydrate-active enzymes varied between female and male quail. Bacteroidetes was the predominant ARG-containing phylum in female quail; Euryarchaeota was the predominant ARG-containing phylum in male quail. Conclusion: This article provides the first description of the gene catalog of the cecal bacteria in Japanese quail as well as insights into the bacterial taxa and predictive metagenomic functions between male and female quail to provide a better understanding of the microbial genes in the quail ceca.

Impact of GC content on gene expression pattern in chicken.

BACKGROUND: GC content varies greatly between different genomic regions in many eukaryotes. In order to determine whether this organization named isochore organization influences gene expression patterns, the relationship between GC content and gene expression has been investigated in man and mouse. However, to date, this question is still a matter for debate. Among the avian species, chicken (Gallus gallus) is the best studied representative with a complete genome sequence. The distinctive features and organization of its sequence make it a good model to explore important issues in genome structure and evolution. METHODS: Only nuclear genes with complete information on protein-coding sequence with no evidence of multiple-splicing forms were included in this study. Chicken protein coding sequences, complete mRNA sequences (or full length cDNA sequences), and 5' untranslated region sequences (5' UTR) were downloaded from Ensembl and chicken expression data originated from a previous work. Three indices i.e. expression level, expression breadth and maximum expression level were used to measure the expression pattern of a given gene. CpG islands were identified using hgTables of the UCSC Genome Browser. Correlation analysis between variables was performed by SAS Proprietary Software Release 8.1. RESULTS: In chicken, the GC content of 5' UTR is significantly and positively correlated with expression level, expression breadth, and maximum expression level, whereas that of coding sequences and introns and at the third coding position are negatively correlated with expression level and expression breadth, and not correlated with maximum expression level. These significant trends are independent of recombination rate, chromosome size and gene density. Furthermore, multiple linear regression analysis indicated that GC content in genes could explain approximately 10% of the variation in gene expression. CONCLUSIONS: GC content is significantly associated with gene expression pattern and could be one of the important regulation factors in the chicken genome.

Indel segregating within introns in the chicken genome are positively correlated with the recombination rates.

Insertions and deletions (Indel) are important sources of genetic diversity and phenotypic divergence. Many factors such as mutation, recombination, selection and genetic drift can jointly affect the indel distribution across the genome. Studies of the relationship between recombination and indel density can, to a certain extent, reflect the selective constrain on indel. Based on the improved genetic map, genome sequence assembly and the partial (0.25X) shotgun sequencing of three breeds of domestic chicken, we calculated the recombination rates and the indel density segregating within introns and intergenic for 4 Mb windows (n = 210). Regression analyses demonstrated that recombination rates are significantly correlated with intron indel density, but not with the intergenic indel density. After adjusted regional effect, the significant trend was remained. This implies that selection is an important factor to influence the indel distribution within introns in chicken genome. By contrast, the intergenic indel seem to be neutral. Since the intron indel density on Z chromosome is less than half of that on autosomes, we preliminarily deduced that genetic hitchhiking might be more important than background selection in producing the observed correlation. As these two processes are not mutually exclusive, it is most likely that both contribute somewhat to the observed pattern. In result similar to previous study, we also found SNP density is highly correlated with indel density. Based on this characteristic, a hypothesis suggested that there are common effects of mutation and/or selection on the occurrence of indel and point mutations. This hypothesis can not explain our observations.

Extent of linkage disequilibrium in wild and domestic chicken populations.

Linkage disequilibrium (LD) analyses play a fundamental role in gene mapping, both as a tool for fine mapping of complex trait gene and in genome-wide association studies. The use of LD analyses in practice depends crucially on the understanding of the patterns of LD in the genome. In the present study, a total of 36 SNP were selected initially in a region (200 kb) of Contig.060226.1 on GGA1 based on the average physical distance. After verifying their level of polymorphism, 21 SNP were selected finally to genotype one wild and two domestic chicken populations, Red Jungle Fowl (RJF), Taihe Silkie chicken (TS) and White Recessive Rock chicken (WRR). Two distinct measures of linkage disequilibrium, D' and r(2), between marker pairs were used. The D'map of RJF is spurious. Many marker pairs showed complete LD. TS and WRR showed distinct characteristic of decreasing D' value over increasing physical distance. The r(2) showed much less inflation than did D' in RJF population, and also showed a characteristic decreasing value over increasing physical distance. In TS and WRR populations, although the two measures differed in scale, their decay profiles were similar. The data in the present study suggested that the extent of LD in this region is about 150 kb, corresponding to 0.4 cM. Our results imply that a very dense map of SNP markers will be required for LD mapping methods. Thus, association studies based on polymorphisms within all known genes, and candidate QTL mapping, may ultimately prove to be a more effective strategy.

[Comparative study of SNP diversity and calculation of the effective size of population in chicken].

A region (200 kb) of Contig. 060226.1 on GGA1 was selected to study the average diversity of Red Jungle Fowl (RJF), Taihe Silk chicken (TS), and White Recessive Rock (WRR), and estimate the effective population size (Ne) of chicken. The mean heterozygosities of RJF, TS and WRR were 0.28533+/-0.034747, 0.32926+/-0.039191 and 0.30168+/-0.040382, respectively. Significant test indicted that there is not significant difference between them (P=0.2368>0.05). The initial chicken effective population size was estimated to be about 20 000-150 000. Chicken has undergone serious population bottleneck effect during the first stage of domestication. However, this bottleneck effect did not result in a substantial loss of diversity as revealed by SNP. The possible explanations for similar SNP diversity between domesticated chicken and Red Jungle Fowl might due to many factors, such as a big Ne in earlier phase of domestication, population expending in breed differentiation, abroad crossing between breeds (especially crossing with RJF), together with high recombination rate in chicken genome.

[Single nucleotide polymorphisms and fine mapping QTL in chickens].

As an important economically animal and a model animal, a great of progress of SNP diversity and QTL mapping in chickens have been made in recent years. The present paper not only summarized these achievements but also commanded an extensive view about the fine mapping of QTL by the use of SNP marker.