Host genetics drives differences in cecal microbiota composition and immune traits of laying hens raised in the same environment.

Vaccination is one of the most effective strategies for preventing infectious diseases but individual vaccine responses are highly heterogeneous. Host genetics and gut microbiota composition are 2 likely drivers of this heterogeneity. We studied 94 animals belonging to 4 lines of laying hens: a White Leghorn experimental line genetically selected for a high antibody response against the Newcastle Disease Virus (NDV) vaccine (ND3) and its unselected control line (CTR), and 2 commercial lines (White Leghorn [LEG] and Rhode Island Red [RIR]). Animals were reared in the same conditions from hatching to 42 d of age, and animals from different genetic lines were mixed. Animals were vaccinated at 22 d of age and their humoral vaccine response against NDV was assessed by hemagglutination inhibition assay and ELISA from blood samples collected at 15, 19, and 21 d after vaccination. The immune parameters studied were the 3 immunoglobulins subtypes A, M, and Y and the blood cell composition was assessed by flow cytometry. The composition of the cecal microbiota was assessed at the end of the experiment by analyzing amplified 16S rRNA gene sequences to obtain amplicon sequence variants (ASV). The 4 lines showed significantly different levels of NDV vaccine response at the 3 measured points, with, logically, a higher response of the genetically selected ND3 line, and intermediate and low responses for the unselected CTR control line and for the 2 commercial lines, respectively. The ND3 line displayed also a higher proportion of immunoglobulins (IgA, IgM, and IgY). The RIR line showed the most different blood cell composition. The 4 lines showed significantly different microbiota characteristics: composition, abundances at all taxonomic levels, and correlations between genera and vaccine response. The tested genetic lines differ for immune parameters and gut microbiota composition and functions. These phenotypic differences can be attributed to genetic differences between lines. Causal relationships between both types of parameters are discussed and will be investigated in further studies.

A frame-shift mutation in COMTD1 is associated with impaired pheomelanin pigmentation in chicken.

The biochemical pathway regulating the synthesis of yellow/red pheomelanin is less well characterized than the synthesis of black/brown eumelanin. Inhibitor of gold (IG phenotype) is a plumage colour variant in chicken that provides an opportunity to further explore this pathway since the recessive allele (IG) at this locus is associated with a defect in the production of pheomelanin. IG/IG homozygotes display a marked dilution of red pheomelanin pigmentation, whilst black pigmentation (eumelanin) is only slightly affected. Here we show that a 2-base pair insertion (frame-shift mutation) in the 5th exon of the Catechol-O-methyltransferase containing domain 1 gene (COMTD1), expected to cause a complete or partial loss-of-function of the COMTD1 enzyme, shows complete concordance with the IG phenotype within and across breeds. We show that the COMTD1 protein is localized to mitochondria in pigment cells. Knockout of Comtd1 in a mouse melanocytic cell line results in a reduction in pheomelanin metabolites and significant alterations in metabolites of glutamate/glutathione, riboflavin, and the tricarboxylic acid cycle. Furthermore, COMTD1 overexpression enhanced cellular proliferation following chemical-induced transfection, a potential inducer of oxidative stress. These observations suggest that COMTD1 plays a protective role for melanocytes against oxidative stress and that this supports their ability to produce pheomelanin.

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.

The feather pattern autosomal barring in chicken is strongly associated with segregation at the MC1R locus.

Color patterns within individual feathers are common in birds but little is known about the genetic mechanisms causing such patterns. Here, we investigate the genetic basis for autosomal barring in chicken, a horizontal striping pattern on individual feathers. Using an informative backcross, we demonstrate that the MC1R locus is strongly associated with this phenotype. A deletion at SOX10, underlying the dark brown phenotype on its own, affects the manifestation of the barring pattern. The coding variant L133Q in MC1R is the most likely causal mutation for autosomal barring in this pedigree. Furthermore, a genetic screen across six different breeds showing different patterning phenotypes revealed that the most striking shared characteristics among these breeds were that they all carried the MC1R alleles Birchen or brown. Our data suggest that the presence of activating MC1R mutations enhancing pigment synthesis is an important mechanism underlying pigmentation patterns on individual feathers in chicken. We propose that MC1R and its antagonist ASIP play a critical role for determining within-feather pigmentation patterns in birds by acting as activator and inhibitor possibly in a Turing reaction-diffusion model.

Chickens can durably clear herpesvirus vaccine infection in feathers while still carrying vaccine-induced antibodies.

Marek's disease (MD) is a major disease of chickens induced by Marek's disease virus (MDV) associated to lethal lymphomas. Current MD vaccines protect against lymphomas, but fail to prevent infection and shedding. The control of MDV shedding is crucial in order to eradicate this highly contagious virus. Like pathogenic MDV, MD vaccines infect the feather follicles of the skin before being shed into the environment. MD vaccines constitute excellent models to study virus interaction with feathers, the unique excretion source of these viruses. Herein we studied the viral persistence in feathers of a MD vaccine, the recombinant turkey herpesvirus (rHVT-ND). We report that most of the birds showed a persistent HVT infection of feathers over 41 weeks with moderate viral loads. Interestingly, 20% of the birds were identified as low HVT producers, among which six birds cleared the infection. Indeed, after week 14-26, these birds named controllers had undetectable HVT DNA in their feathers through week 41. All vaccinated birds developed antibodies to NDV, which lasted until week 41 in 95% of the birds, including the controllers. No correlation was found between HVT loads in feathers and NDV antibody titers over time. Interestingly, no HVT DNA was detected in the spleens of four controllers. This is the first description of chickens that durably cleared MD vaccine infection of feathers suggesting that control of Mardivirus shedding is achievable by the host.

The quail genome: insights into social behaviour, seasonal biology and infectious disease response.

BACKGROUND: The Japanese quail (Coturnix japonica) is a popular domestic poultry species and an increasingly significant model species in avian developmental, behavioural and disease research. RESULTS: We have produced a high-quality quail genome sequence, spanning 0.93 Gb assigned to 33 chromosomes. In terms of contiguity, assembly statistics, gene content and chromosomal organisation, the quail genome shows high similarity to the chicken genome. We demonstrate the utility of this genome through three diverse applications. First, we identify selection signatures and candidate genes associated with social behaviour in the quail genome, an important agricultural and domestication trait. Second, we investigate the effects and interaction of photoperiod and temperature on the transcriptome of the quail medial basal hypothalamus, revealing key mechanisms of photoperiodism. Finally, we investigate the response of quail to H5N1 influenza infection. In quail lung, many critical immune genes and pathways were downregulated after H5N1 infection, and this may be key to the susceptibility of quail to H5N1. CONCLUSIONS: We have produced a high-quality genome of the quail which will facilitate further studies into diverse research questions using the quail as a model avian species.

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.

A missense mutation in TYRP1 causes the chocolate plumage color in chicken and alters melanosome structure.

The chocolate plumage color in chickens is due to a sex-linked recessive mutation, choc, which dilutes eumelanin pigmentation. Because TYRP1 is sex-linked in chickens, and TYRP1 mutations determine brown coat color in mammals, TYRP1 appeared as the obvious candidate gene for the choc mutation. By combining gene mapping with gene capture, a complete association was identified between the chocolate phenotype and a missense mutation leading to a His214Asn change in the ZnA zinc-binding domain of the protein. A diagnostic test confirmed complete association by screening 428 non-chocolate chickens of various origins. This is the first TYRP1 mutation described in the chicken. Electron microscopy analysis showed that melanosomes were more numerous in feather follicles of chocolate chickens but exhibited an abnormal structure characterized by a granular content and an irregular shape. A similar altered morphology was observed on melanosomes of another TYRP1 mutant in birds, the roux mutation of the quail.

Accounting for linkage disequilibrium in genome scans for selection without individual genotypes: The local score approach.

Detecting genomic footprints of selection is an important step in the understanding of evolution. Accounting for linkage disequilibrium in genome scans increases detection power, but haplotype-based methods require individual genotypes and are not applicable on pool-sequenced samples. We propose to take advantage of the local score approach to account for linkage disequilibrium in genome scans for selection, cumulating (possibly small) signals from single markers over a genomic segment, to clearly pinpoint a selection signal. Using computer simulations, we demonstrate that this approach detects selection with higher power than several state-of-the-art single-marker, windowing or haplotype-based approaches. We illustrate this on two benchmark data sets including individual genotypes, for which we obtain similar results with the local score and one haplotype-based approach. Finally, we apply the local score approach to Pool-Seq data obtained from a divergent selection experiment on behaviour in quail and obtain precise and biologically coherent selection signals: while competing methods fail to highlight any clear selection signature, our method detects several regions involving genes known to act on social responsiveness or autistic traits. Although we focus here on the detection of positive selection from multiple population data, the local score approach is general and can be applied to other genome scans for selection or other genomewide analyses such as GWAS.

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.

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.

Chicken skin virome analyzed by high-throughput sequencing shows a composition highly different from human skin.

Recent studies show that human skin at homeostasis is a complex ecosystem whose virome include circular DNA viruses, especially papillomaviruses and polyomaviruses. To determine the chicken skin virome in comparison with human skin virome, a chicken swabs pool sample from fifteen indoor healthy chickens of five genetic backgrounds was examined for the presence of DNA viruses by high-throughput sequencing (HTS). The results indicate a predominance of herpesviruses from the Mardivirus genus, coming from either vaccinal origin or presumably asymptomatic infection. Despite the high sensitivity of the HTS method used herein to detect small circular DNA viruses, we did not detect any papillomaviruses, polyomaviruses, or circoviruses, indicating that these viruses may not be resident of the chicken skin. The results suggest that the turkey herpesvirus is a resident of chicken skin in vaccinated chickens. This study indicates major differences between the skin viromes of chickens and humans. The origin of this difference remains to be further studied in relation with skin physiology, environment, or virus population dynamics.

Genome-Wide Characterization of RNA Editing in Chicken Embryos Reveals Common Features among Vertebrates.

RNA editing results in a post-transcriptional nucleotide change in the RNA sequence that creates an alternative nucleotide not present in the DNA sequence. This leads to a diversification of transcription products with potential functional consequences. Two nucleotide substitutions are mainly described in animals, from adenosine to inosine (A-to-I) and from cytidine to uridine (C-to-U). This phenomenon is described in more details in mammals, notably since the availability of next generation sequencing technologies allowing whole genome screening of RNA-DNA differences. The number of studies recording RNA editing in other vertebrates like chicken is still limited. We chose to use high throughput sequencing technologies to search for RNA editing in chicken, and to extend the knowledge of its conservation among vertebrates. We performed sequencing of RNA and DNA from 8 embryos. Being aware of common pitfalls inherent to sequence analyses that lead to false positive discovery, we stringently filtered our datasets and found fewer than 40 reliable candidates. Conservation of particular sites of RNA editing was attested by the presence of 3 edited sites previously detected in mammals. We then characterized editing levels for selected candidates in several tissues and at different time points, from 4.5 days of embryonic development to adults, and observed a clear tissue-specificity and a gradual increase of editing level with time. By characterizing the RNA editing landscape in chicken, our results highlight the extent of evolutionary conservation of this phenomenon within vertebrates, attest to its tissue and stage specificity and provide support of the absence of non A-to-I events from the chicken transcriptome.

Quantitative effect of a CNV on a morphological trait in chickens.

Copy Number Variation has been associated with morphological traits, developmental defects or disease susceptibility. The autosomal dominant Pea-comb mutation in chickens is due to the massive amplification of a CNV in intron 1 of SOX5 and provides a unique opportunity to assess the effect of variation in the number of repeats on quantitative traits such as comb size and comb mass in Pea-comb chickens. The quantitative variation of comb size was estimated by 2D morphometry and the number of repeats (RQ) was estimated by qPCR, in a total of 178 chickens from 3 experimental lines, two of them showing segregation for the Pea-comb mutation. This study included only Pea-comb chickens. Analysis of variance showed highly significant effects of line and sex on comb measurements. Adult body weight (BW) and RQ were handled as covariates. BW significantly influenced comb mass but not comb size. RQ values significantly influenced comb size, and the linear regression coefficient was highest for heterozygous carriers: the higher the number of repeats, the smaller the comb size. A similar trend was observed for comb mass. The CNV contributed to 3.4% of the phenotypic variance of comb size in heterozygous carriers of the CNV, an order of magnitude frequently encountered for QTLs. Surprisingly, there was no such relationship between RQ values and comb size in the homozygous line. It may be concluded that heterozygosity for a CNV in a non-coding region may contribute to phenotypic plasticity.

A genomic duplication is associated with ectopic eomesodermin expression in the embryonic chicken comb and two duplex-comb phenotypes.

Duplex-comb (D) is one of three major loci affecting comb morphology in the domestic chicken. Here we show that the two Duplex-comb alleles, V-shaped (D*V) and Buttercup (D*C), are both associated with a 20 Kb tandem duplication containing several conserved putative regulatory elements located 200 Kb upstream of the eomesodermin gene (EOMES). EOMES is a T-box transcription factor that is involved in mesoderm specification during gastrulation. In D*V and D*C chicken embryos we find that EOMES is ectopically expressed in the ectoderm of the comb-developing region as compared to wild-type embryos. The confinement of the ectopic expression of EOMES to the ectoderm is in stark contrast to the causal mechanisms underlying the two other major comb loci in the chicken (Rose-comb and Pea-comb) in which the transcription factors MNR2 and SOX5 are ectopically expressed strictly in the mesenchyme. Interestingly, the causal mutations of all three major comb loci in the chicken are now known to be composed of large-scale structural genomic variants that each result in ectopic expression of transcription factors. The Duplex-comb locus also illustrates the evolution of alleles in domestic animals, which means that alleles evolve by the accumulation of two or more consecutive mutations affecting the phenotype. We do not yet know whether the V-shaped or Buttercup allele correspond to the second mutation that occurred on the haplotype of the original duplication event.

A medium density genetic map and QTL for behavioral and production traits in Japanese quail.

BACKGROUND: Behavioral traits such as sociability, emotional reactivity and aggressiveness are major factors in animal adaptation to breeding conditions. In order to investigate the genetic control of these traits as well as their relationships with production traits, a study was undertaken on a large second generation cross (F2) between two lines of Japanese Quail divergently selected on their social reinstatement behavior. All the birds were measured for several social behaviors (social reinstatement, response to social isolation, sexual motivation, aggression), behaviors measuring the emotional reactivity of the birds (reaction to an unknown object, tonic immobility reaction), and production traits (body weight and egg production). RESULTS: We report the results of the first genome-wide QTL detection based on a medium density SNP panel obtained from whole genome sequencing of a pool of individuals from each divergent line. A genetic map was constructed using 2145 markers among which 1479 could be positioned on 28 different linkage groups. The sex-averaged linkage map spanned a total of 3057 cM with an average marker spacing of 2.1 cM. With the exception of a few regions, the marker order was the same in Japanese Quail and the chicken, which confirmed a well conserved synteny between the two species. The linkage analyses performed using QTLMAP software revealed a total of 45 QTLs related either to behavioral (23) or production (22) traits. The most numerous QTLs (15) concerned social motivation traits. Interestingly, our results pinpointed putative pleiotropic regions which controlled emotional reactivity and body-weight of birds (on CJA5 and CJA8) or their social motivation and the onset of egg laying (on CJA19). CONCLUSION: This study identified several QTL regions for social and emotional behaviors in the Quail. Further research will be needed to refine the QTL and confirm or refute the role of candidate genes, which were suggested by bioinformatics analysis. It can be hoped that the identification of genes and polymorphisms related to behavioral traits in the quail will have further applications for other poultry species (especially the chicken) and will contribute to solving animal welfare issues in poultry production.

QTL detection for coccidiosis (Eimeria tenella) resistance in a Fayoumi × Leghorn F2 cross, using a medium-density SNP panel.

BACKGROUND: Coccidiosis is a major parasitic disease that causes huge economic losses to the poultry industry. Its pathogenicity leads to depression of body weight gain, lesions and, in the most serious cases, death in affected animals. Genetic variability for resistance to coccidiosis in the chicken has been demonstrated and if this natural resistance could be exploited, it would reduce the costs of the disease. Previously, a design to characterize the genetic regulation of Eimeria tenella resistance was set up in a Fayoumi × Leghorn F2 cross. The 860 F2 animals of this design were phenotyped for weight gain, plasma coloration, hematocrit level, intestinal lesion score and body temperature. In the work reported here, the 860 animals were genotyped for a panel of 1393 (157 microsatellites and 1236 single nucleotide polymorphism (SNP) markers that cover the sequenced genome (i.e. the 28 first autosomes and the Z chromosome). In addition, with the aim of finding an index capable of explaining a large amount of the variance associated with resistance to coccidiosis, a composite factor was derived by combining the variables of all these traits in a single variable. QTL detection was performed by linkage analysis using GridQTL and QTLMap. Single and multi-QTL models were applied. RESULTS: Thirty-one QTL were identified i.e. 27 with the single-QTL model and four with the multi-QTL model and the average confidence interval was 5.9 cM. Only a few QTL were common with the previous study that used the same design but focused on the 260 more extreme animals that were genotyped with the 157 microsatellites only. Major differences were also found between results obtained with QTLMap and GridQTL. CONCLUSIONS: The medium-density SNP panel made it possible to genotype new regions of the chicken genome (including micro-chromosomes) that were involved in the genetic control of the traits investigated. This study also highlights the strong variations in QTL detection between different models and marker densities.

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.

Evidence of phenotypic and genetic relationships between sociality, emotional reactivity and production traits in Japanese quail.

The social behavior of animals, which is partially controlled by genetics, is one of the factors involved in their adaptation to large breeding groups. To understand better the relationships between different social behaviors, fear behaviors and production traits, we analyzed the phenotypic and genetic correlations of these traits in Japanese quail by a second generation crossing of two lines divergently selected for their social reinstatement behavior. Analyses of results for 900 individuals showed that the phenotypic correlations between behavioral traits were low with the exception of significant correlations between sexual behavior and aggressive pecks both at phenotypic (0.51) and genetic (0.90) levels. Significant positive genetic correlations were observed between emotional reactivity toward a novel object and sexual (0.89) or aggressive (0.63) behaviors. The other genetic correlations were observed mainly between behavioral and production traits. Thus, the level of emotional reactivity, estimated by the duration of tonic immobility, was positively correlated with weight at 17 and 65 days of age (0.76 and 0.79, respectively) and with delayed egg laying onset (0.74). In contrast, a higher level of social reinstatement behavior was associated with an earlier egg laying onset (-0.71). In addition, a strong sexual motivation was correlated with an earlier laying onset (-0.68) and a higher number of eggs laid (0.82). A low level of emotional reactivity toward a novel object and also a higher aggressive behavior were genetically correlated with a higher number of eggs laid (0.61 and 0.58, respectively). These results bring new insights into the complex determinism of social and emotional reactivity behaviors in birds and their relationships with production traits. Furthermore, they highlight the need to combine animal welfare and production traits in selection programs by taking into account traits of sociability and emotional reactivity.