Relationship between feed efficiency and gut microbiota in laying chickens under contrasting feeding conditions.

The gut microbiota is known to play an important role in energy harvest and is likely to affect feed efficiency. In this study, we used 16S metabarcoding sequencing to analyse the caecal microbiota of laying hens from feed-efficient and non-efficient lines obtained by divergent selection for residual feed intake. The two lines were fed either a commercial wheat-soybean based diet (CTR) or a low-energy, high-fibre corn-sunflower diet (LE). The analysis revealed a significant line x diet interaction, highlighting distinct differences in microbial community composition between the two lines when hens were fed the CTR diet, and more muted differences when hens were fed the LE diet. Our results are consistent with the hypothesis that a richer and more diverse microbiota may play a role in enhancing feed efficiency, albeit in a diet-dependent manner. The taxonomic differences observed in the microbial composition seem to correlate with alterations in starch and fibre digestion as well as in the production of short-chain fatty acids. As a result, we hypothesise that efficient hens are able to optimise nutrient absorption through the activity of fibrolytic bacteria such as Alistipes or Anaerosporobacter, which, via their production of propionate, influence various aspects of host metabolism.

RNA-Seq Data for Reliable SNP Detection and Genotype Calling: Interest for Coding Variant Characterization and Cis-Regulation Analysis by Allele-Specific Expression in Livestock Species.

In addition to their common usages to study gene expression, RNA-seq data accumulated over the last 10 years are a yet-unexploited resource of SNPs in numerous individuals from different populations. SNP detection by RNA-seq is particularly interesting for livestock species since whole genome sequencing is expensive and exome sequencing tools are unavailable. These SNPs detected in expressed regions can be used to characterize variants affecting protein functions, and to study cis-regulated genes by analyzing allele-specific expression (ASE) in the tissue of interest. However, gene expression can be highly variable, and filters for SNP detection using the popular GATK toolkit are not yet standardized, making SNP detection and genotype calling by RNA-seq a challenging endeavor. We compared SNP calling results using GATK suggested filters, on two chicken populations for which both RNA-seq and DNA-seq data were available for the same samples of the same tissue. We showed, in expressed regions, a RNA-seq precision of 91% (SNPs detected by RNA-seq and shared by DNA-seq) and we characterized the remaining 9% of SNPs. We then studied the genotype (GT) obtained by RNA-seq and the impact of two factors (GT call-rate and read number per GT) on the concordance of GT with DNA-seq; we proposed thresholds for them leading to a 95% concordance. Applying these thresholds to 767 multi-tissue RNA-seq of 382 birds of 11 chicken populations, we found 9.5 M SNPs in total, of which ∼550,000 SNPs per tissue and population with a reliable GT (call rate ≥ 50%) and among them, ∼340,000 with a MAF ≥ 10%. We showed that such RNA-seq data from one tissue can be used to (i) detect SNPs with a strong predicted impact on proteins, despite their scarcity in each population (16,307 SIFT deleterious missenses and 590 stop-gained), (ii) study, on a large scale, cis-regulations of gene expression, with ∼81% of protein-coding and 68% of long non-coding genes (TPM ≥ 1) that can be analyzed for ASE, and with ∼29% of them that were cis-regulated, and (iii) analyze population genetic using such SNPs located in expressed regions. This work shows that RNA-seq data can be used with good confidence to detect SNPs and associated GT within various populations and used them for different analyses as GTEx studies.

Assessment of trade-offs between feed efficiency, growth-related traits, and immune activity in experimental lines of layer chickens.

BACKGROUND: In all organisms, life-history traits are constrained by trade-offs, which may represent physiological limitations or be related to energy resource management. To detect trade-offs within a population, one promising approach is the use of artificial selection, because intensive selection on one trait can induce unplanned changes in others. In chickens, the breeding industry has achieved remarkable genetic progress in production and feed efficiency over the last 60 years. However, this may have been accomplished at the expense of other important biological functions, such as immunity. In the present study, we used three experimental lines of layer chicken-two that have been divergently selected for feed efficiency and one that has been selected for increased antibody response to inactivated Newcastle disease virus (ND3)-to explore the impact of improved feed efficiency on animals' immunocompetence and, vice versa, the impact of improved antibody response on animals' growth and feed efficiency. RESULTS: There were detectable differences between the low (R+) and high (R-) feed-efficiency lines with respect to vaccine-specific antibody responses and counts of monocytes, heterophils, and/or T cell population. The ND3 line presented reduced body weight and feed intake compared to the control line. ND3 chickens also demonstrated an improved antibody response against a set of commercial viral vaccines, but lower blood leucocyte counts. CONCLUSIONS: This study demonstrates the value of using experimental chicken lines that are divergently selected for RFI or for a high antibody production, to investigate the modulation of immune parameters in relation to growth and feed efficiency. Our results provide further evidence that long-term selection for the improvement of one trait may have consequences on other important biological functions. Hence, strategies to ensure optimal trade-offs among competing functions will ultimately be required in multi-trait selection programs in livestock.

Unraveling the history of the genus Gallus through whole genome sequencing.

The genus Gallus is distributed across a large part of Southeast Asia and has received special interest because the domestic chicken, Gallus gallus domesticus, has spread all over the world and is a major protein source for humans. There are four species: the red junglefowl (G. gallus), the green junglefowl (G. varius), the Lafayette's junglefowl (G. lafayettii) and the grey junglefowl (G. sonneratii). The aim of this study is to reconstruct the history of these species by a whole genome sequencing approach and resolve inconsistencies between well supported topologies inferred using different data and methods. Using deep sequencing, we identified over 35 million SNPs and reconstructed the phylogeny of the Gallus genus using both distance (BioNJ) and maximum likelihood (ML) methods. We observed discrepancies according to reconstruction methods and genomic components. The two most supported topologies were previously reported and were discriminated by using phylogenetic and gene flow analyses, based on ABBA statistics. Terminology fix requested by the deputy editor led to support a scenario with G. gallus as the earliest branching lineage of the Gallus genus, instead of G. varius. We discuss the probable causes for the discrepancy. A likely one is that G. sonneratii samples from parks or private collections are all recent hybrids, with roughly 10% of their autosomal genome originating from G. gallus. The removal of those regions is needed to provide reliable data, which was not done in previous studies. We took care of this and additionally included two wild G. sonneratii samples from India, showing no trace of introgression. This reinforces the importance of carefully selecting and validating samples and genomic components in phylogenomics.

An integrative atlas of chicken long non-coding genes and their annotations across 25 tissues.

Long non-coding RNAs (LNC) regulate numerous biological processes. In contrast to human, the identification of LNC in farm species, like chicken, is still lacunar. We propose a catalogue of 52,075 chicken genes enriched in LNC ( http://www.fragencode.org/ ), built from the Ensembl reference extended using novel LNC modelled here from 364 RNA-seq and LNC from four public databases. The Ensembl reference grew from 4,643 to 30,084 LNC, of which 59% and 41% with expression ≥ 0.5 and ≥ 1 TPM respectively. Characterization of these LNC relatively to the closest protein coding genes (PCG) revealed that 79% of LNC are in intergenic regions, as in other species. Expression analysis across 25 tissues revealed an enrichment of co-expressed LNC:PCG pairs, suggesting co-regulation and/or co-function. As expected LNC were more tissue-specific than PCG (25% vs. 10%). Similarly to human, 16% of chicken LNC hosted one or more miRNA. We highlighted a new chicken LNC, hosting miR155, conserved in human, highly expressed in immune tissues like miR155, and correlated with immunity-related PCG in both species. Among LNC:PCG pairs tissue-specific in the same tissue, we revealed an enrichment of divergent pairs with the PCG coding transcription factors, as for example LHX5, HXD3 and TBX4, in both human and chicken.

Embryonic thermal manipulation has short and long-term effects on the development and the physiology of the Japanese quail.

In vertebrates, the embryonic environment is known to affect the development and the health of individuals. In broiler chickens, the thermal-manipulation (TM) of eggs during the incubation period was shown to improve heat tolerance at slaughter age (35 days of age) in association with several modifications at the molecular, metabolic and physiological levels. However, little is known about the Japanese quail (Coturnix japonica), a closely related avian species widely used as a laboratory animal model and farmed for its meat and eggs. Here we developed and characterized a TM procedure (39.5°C and 65% relative humidity, 12 h/d, from days 0 to 13 of incubation) in quails by analyzing its short and long-term effects on zootechnical, physiological and metabolic parameters. Heat-tolerance was tested by a heat challenge (36°C for 7h) at 35 days of age. TM significantly reduced the hatching rate of the animals and increased mortality during the first four weeks of life. At hatching, TM animals were heavier than controls, but lighter at 25 days of age for both sexes. Thirty-five days after hatching, TM decreased the surface temperature of the shank in females, suggesting a modulation of the blood flow to maintain the internal temperature. TM also increased blood partial pressure and oxygen saturation percentage at 35 days of age in females, suggesting a long-term modulation of the respiration physiology. Quails physiologically responded to the heat challenge, with a modification of several hematologic and metabolic parameters, including an increase in plasma corticosterone concentration. Several physiological parameters such as beak surface temperature and blood sodium concentration revealed that TM birds responded differently to the heat challenge compared to controls. Altogether, this first comprehensive characterization of TM in Japanese quail showed durable effects that may affect the response of TM quails to heat.

Long noncoding RNAs in lipid metabolism: literature review and conservation analysis across species.

BACKGROUND: Lipids are important for the cell and organism life since they are major components of membranes, energy reserves and are also signal molecules. The main organs for the energy synthesis and storage are the liver and adipose tissue, both in humans and in more distant species such as chicken. Long noncoding RNAs (lncRNAs) are known to be involved in many biological processes including lipid metabolism. RESULTS: In this context, this paper provides the most exhaustive list of lncRNAs involved in lipid metabolism with 60 genes identified after an in-depth analysis of the bibliography, while all "review" type articles list a total of 27 genes. These 60 lncRNAs are mainly described in human or mice and only a few of them have a precise described mode-of-action. Because these genes are still named in a non-standard way making such a study tedious, we propose a standard name for this list according to the rules dictated by the HUGO consortium. Moreover, we identified about 10% of lncRNAs which are conserved between mammals and chicken and 2% between mammals and fishes. Finally, we demonstrated that two lncRNA were wrongly considered as lncRNAs in the literature since they are 3' extensions of the closest coding gene. CONCLUSIONS: Such a lncRNAs catalogue can participate to the understanding of the lipid metabolism regulators; it can be useful to better understand the genetic regulation of some human diseases (obesity, hepatic steatosis) or traits of economic interest in livestock species (meat quality, carcass composition). We have no doubt that this first set will be rapidly enriched in coming years.

Two new structural mutations in the 5' region of the ASIP gene cause diluted feather color phenotypes in Japanese quail.

BACKGROUND: In quail, two feather colour phenotypes i.e. fawn-2/beige and yellow are associated with the ASIP locus. The aim of our study was to characterize the structural modifications within this locus that explain the yellow mutation (large deletion) and the fawn-2/beige mutation (assumed to be caused by a different structural modification). RESULTS: For the yellow phenotype, we identified a complex mutation that involves a 141,162-bp long deletion. For the fawn-2/beige phenotype, we identified a 71-kb tandem duplication that comprises one unchanged copy of ASIP and one copy present in the ITCH-ASIP fusion gene, which leads to a transcript coding for a normal ASIP protein. Although this agrees with previous reports that reported an increased level of ASIP transcripts in the skin of mutant animals, we show that in the skin from fawn-2/beige embryos, this level is higher than expected with a simple duplication of the ASIP gene. Thus, we hypothesize that the 5' region of the ITCH-ASIP fusion gene leads to a higher transcription level than the 5' region of the ASIP gene. CONCLUSIONS: We were able to conclude that the fawn-2 and beige phenotypes are caused by the same allele at the ASIP locus. Both of the associated mutations fawn-2/beige and yellow lead to the formation of a fusion gene, which encodes a transcript for the ASIP protein. In both cases, transcription of ASIP depends on the promoter of a different gene, which includes alternative up-regulating sequences. However, we cannot exclude the possibility that the loss of the 5' region of the ASIP gene itself has additional impacts, especially for the fawn-2/beige mutation. In addition, in several other species including mammals, the existence of other dominant gain-of-function structural modifications that are localized upstream of the ASIP coding sequences has been reported, which supports our hypothesis that repressors in the 5' region of ASIP are absent in the fawn-2/beige mutant.

Chicken semen cryopreservation and use for the restoration of rare genetic resources.

For the past 50 yr, practices for ex situ preservation of endangered breeds have been extended. Semen and primordial germ cells, gonadic tissues have been frozen to create genetic stocks of chicken genetic diversity in cryobanks. Semen cryopreservation stays the preferred method since it is not invasive. Many protocols have been developed to cryopreserve chicken semen, but they give highly variable success rate. The aim of the present study was to standardize and prove the effectiveness of semen long-term storage for the restitution of lost families. We showed that semen straws stored for 18 yr in liquid nitrogen did not lose their fertilizing ability. We demonstrated the usefulness of cryopreserved semen stored in the French National Cryobank for the recovery of families of a subfertile experimental chicken line. In order to highlight the standardization of the cryopreserved method, different cryoprotectant protocols were also tested on a rare breed, freezing/thawing and insemination conditions were controlled. The best results were obtained using glycerol protocol, a sperm dilution of 1:4 (semen:extender). The insemination dose of 200 million sperm/female was as efficient as 400 million of sperm. Altogether, these results demonstrated the effectiveness of chicken semen long-term storage for the restoration of lost genetic resources and highlighted the importance of standardized chicken semen cryopreservation using procedures combining biophysical (cryoprotectants, freezing/thawing conditions) and zootechnical (artificial insemination) features.

Cancer- and behavior-related genes are targeted by selection in the Tasmanian devil (Sarcophilus harrisii).

Devil Facial Tumor Disease (DFTD) is an aggressive cancer notorious for its rare etiology and its impact on Tasmanian devil populations. Two regions underlying an evolutionary response to this cancer were recently identified using genomic time-series pre- and post-DTFD arrival. Here, we support that DFTD shaped the genome of the Tasmanian devil in an even more extensive way than previously reported. We detected 97 signatures of selection, including 148 protein coding genes having a human orthologue, linked to DFTD. Most candidate genes are associated with cancer progression, and an important subset of candidate genes has additional influence on social behavior. This confirms the influence of cancer on the ecology and evolution of the Tasmanian devil. Our work also demonstrates the possibility to detect highly polygenic footprints of short-term selection in very small populations.

Genomics for Ruminants in Developing Countries: From Principles to Practice.

Using genomic information, local ruminant populations can be better characterized and compared to selected ones. Genetic relationships between animals can be established even without systematic pedigree recording, provided a budget is available for genotyping. Genomic selection (GS) can rely on a subset of the total population and does not require a costly national infrastructure, e.g., based on progeny testing. Yet, the use of genomic tools for animal breeding in developing countries is still limited. We identify three main reasons for this: (i) the instruments for cheap recording of phenotypes and data management are still limiting. (ii) many developing countries are recurrently exposed to unfavorable conditions (heat, diseases, poor nutrition) requiring special attention to fitness traits, (iii) a high level of expertise in quantitative genetics, modeling, and data manipulation is needed to perform genomic analyses. Yet, the potential outcomes go much beyond genetic improvements and can improve the resilience of the whole farming system. They include a better management of genetic diversity of local populations, a more balanced genetic progress and the possibility to unravel the genetic basis of adaptation of local breeds through whole genome approaches. A GS program being developed by BAIF, a large Indian NGO, is analyzed as a pilot case. It relies on the creation of a female reference population of Bos indicus and crossbreds, recorded with modern technology (e.g., smartphones) to collect performances at low cost in tiny herds on production and fertility. Finally, recommendations for the implementation of GS in developing countries are proposed.

Identification of marginal causal relationships in gene networks from observational and interventional expression data.

Causal network inference is an important methodological challenge in biology as well as other areas of application. Although several causal network inference methods have been proposed in recent years, they are typically applicable for only a small number of genes, due to the large number of parameters to be estimated and the limited number of biological replicates available. In this work, we consider the specific case of transcriptomic studies made up of both observational and interventional data in which a single gene of biological interest is knocked out. We focus on a marginal causal estimation approach, based on the framework of Gaussian directed acyclic graphs, to infer causal relationships between the knocked-out gene and a large set of other genes. In a simulation study, we found that our proposed method accurately differentiates between downstream causal relationships and those that are upstream or simply associative. It also enables an estimation of the total causal effects between the gene of interest and the remaining genes. Our method performed very similarly to a classical differential analysis for experiments with a relatively large number of biological replicates, but has the advantage of providing a formal causal interpretation. Our proposed marginal causal approach is computationally efficient and may be applied to several thousands of genes simultaneously. In addition, it may help highlight subsets of genes of interest for a more thorough subsequent causal network inference. The method is implemented in an R package called MarginalCausality (available on GitHub).

Embryonic environment and transgenerational effects in quail.

BACKGROUND: Environmental exposures, for instance to chemicals, are known to impact plant and animal phenotypes on the long term, sometimes across several generations. Such transgenerational phenotypes were shown to be promoted by epigenetic alterations such as DNA methylation, an epigenetic mark involved in the regulation of gene expression. However, it is yet unknown whether transgenerational epigenetic inheritance of altered phenotypes exists in birds. The purpose of this study was to develop an avian model to investigate whether changes to the embryonic environment had a transgenerational effect that could alter the phenotypes of third-generation offspring. Given its impact on the mammalian epigenome and the reproductive system in birds, genistein was used as an environment stressor. RESULTS: We compared several third-generation phenotypes of two quail "epilines", which were obtained from genistein-injected eggs (Epi+) or from untreated eggs (Epi-) from the same founders. A "mirrored" crossing strategy was used to minimize between-line genetic variability by maintaining similar ancestor contributions across generations in each line. Three generations after genistein treatment, a significant difference in the sexual maturity of the females, which, after three generations, could not be attributed to direct maternal effects, was observed between the lines, with Epi+ females starting to lay eggs later. Adult body weight was significantly affected by genistein treatment applied in a previous generation, and a significant interaction between line and sex was observed for body weight at 3 weeks. Behavioral traits, such as evaluating the birds' reaction to social isolation, were also significantly affected by genistein treatment. Yet, global methylation analyses revealed no significant difference between the epilines. CONCLUSIONS: These findings demonstrate that embryonic environment affects the phenotype of offspring three generations later in quail. While one cannot rule out the existence of some initial genetic variability between the lines, the mirrored animal design should have minimized its effects, and thus, the observed differences in animals of the third generation may be attributed, at least partly, to transgenerational epigenetic phenomena.

Transcriptomes of whole blood and PBMC in chickens.

Global transcriptome analysis of chicken whole blood to discover biomarkers of different phenotypes or physiological disorders has never been investigated so far. Whole blood provides significant advantages, allowing large scale and non-invasive sampling. However, generation of gene expression data from the blood of non-mammalian species remains a challenge, notably due to the nucleated red blood cells, hindering the use of well-established protocols. The aim of this study was to analyze the relevance of using whole blood cells (WB) to find biomarkers, instead of Peripheral Blood Mononuclear Cells (PBMC), usually chosen for immune challenges. RNA sources from WB and PBMC was characterized by microarray analysis. Our results show that the quality and quantity of RNA obtained from WB was suitable for further analyses, although the quality was lower than that from PBMC. The transcriptome profiling comparison revealed that the majority of genes were expressed in both WB and PBMC. Hemoglobin subunits were the major transcripts in WB, whereas the most enriched biological process was related to protein catabolic process. Most of the over-represented transcripts in PBMC were implicated in functions specific to thrombocytes, like coagulation and platelet activation, probably due to the large proportion of this nucleated cell type in chicken PBMC. Functions related to B and T cells and to other immune functions were also enriched in the PBMC subset. We conclude that WB is more suitable for large scale immunity oriented studies and other biological processes that have been poorly investigated so far.

The Extent of mRNA Editing Is Limited in Chicken Liver and Adipose, but Impacted by Tissular Context, Genotype, Age, and Feeding as Exemplified with a Conserved Edited Site in COG3.

RNA editing is a posttranscriptional process leading to differences between genomic DNA and transcript sequences, potentially enhancing transcriptome diversity. With recent advances in high-throughput sequencing, many efforts have been made to describe mRNA editing at the transcriptome scale, especially in mammals, yielding contradictory conclusions regarding the extent of this phenomenon. We show, by detailed description of the 25 studies focusing so far on mRNA editing at the whole-transcriptome scale, that systematic sequencing artifacts are considered in most studies whereas biological replication is often neglected and multi-alignment not properly evaluated, which ultimately impairs the legitimacy of results. We recently developed a rigorous strategy to identify mRNA editing using mRNA and genomic DNA sequencing, taking into account sequencing and mapping artifacts, and biological replicates. We applied this method to screen for mRNA editing in liver and white adipose tissue from eight chickens and confirm the small extent of mRNA recoding in this species. Among the 25 unique edited sites identified, three events were previously described in mammals, attesting that this phenomenon is conserved throughout evolution. Deeper investigations on five sites revealed the impact of tissular context, genotype, age, feeding conditions, and sex on mRNA editing levels. More specifically, this analysis highlighted that the editing level at the site located on COG3 was strongly regulated by four of these factors. By comprehensively characterizing the mRNA editing landscape in chickens, our results highlight how this phenomenon is limited and suggest regulation of editing levels by various genetic and environmental factors.

Impact of the genetic background on the composition of the chicken plasma MiRNome in response to a stress.

Circulating extra-cellular microRNAs (miRNAs) have emerged as promising minimally invasive markers in human medicine. We evaluated miRNAs isolated from total plasma as biomarker candidates of a response to an abiotic stress (feed deprivation) in a livestock species. Two chicken lines selected for high (R+) and low (R-) residual feed intake were chosen as an experimental model because of their extreme divergence in feed intake and energy metabolism. Adult R+ and R- cocks were sampled after 16 hours of feed deprivation and again four hours after re-feeding. More than 292 million sequence reads were generated by small RNA-seq of total plasma RNA. A total of 649 mature miRNAs were identified; after quality filtering, 148 miRNAs were retained for further analyses. We identified 23 and 19 differentially abundant miRNAs between feeding conditions and between lines respectively, with only two miRNAs identified in both comparisons. We validated a panel of six differentially abundant miRNAs by RT-qPCR on a larger number of plasma samples and checked their response to feed deprivation in liver. Finally, we evaluated the conservation and tissue distribution of differentially abundant miRNAs in plasma across a variety of red jungle fowl tissues. We show that the chicken plasma miRNome reacts promptly to the alteration of the animal physiological condition driven by a feed deprivation stress. The plasma content of stress-responsive miRNAs is strongly influenced by the genetic background, with differences reflecting the phenotypic divergence acquired through long-term selection, as evidenced by the profiles of conserved miRNAs with a regulatory role in energy metabolism (gga-miR-204, gga-miR-let-7f-5p and gga-miR-122-5p). These results reinforce the emerging view in human medicine that even small genetic differences can have a considerable impact on the resolution of biomarker studies, and provide support for the emerging interest in miRNAs as potential novel and minimally invasive biomarkers for livestock species.

Transcriptome-wide investigation of genomic imprinting in chicken.

Genomic imprinting is an epigenetic mechanism by which alleles of some specific genes are expressed in a parent-of-origin manner. It has been observed in mammals and marsupials, but not in birds. Until now, only a few genes orthologous to mammalian imprinted ones have been analyzed in chicken and did not demonstrate any evidence of imprinting in this species. However, several published observations such as imprinted-like QTL in poultry or reciprocal effects keep the question open. Our main objective was thus to screen the entire chicken genome for parental-allele-specific differential expression on whole embryonic transcriptomes, using high-throughput sequencing. To identify the parental origin of each observed haplotype, two chicken experimental populations were used, as inbred and as genetically distant as possible. Two families were produced from two reciprocal crosses. Transcripts from 20 embryos were sequenced using NGS technology, producing ∼200 Gb of sequences. This allowed the detection of 79 potentially imprinted SNPs, through an analysis method that we validated by detecting imprinting from mouse data already published. However, out of 23 candidates tested by pyrosequencing, none could be confirmed. These results come together, without a priori, with previous statements and phylogenetic considerations assessing the absence of genomic imprinting in chicken.