Natural clines and human management impact the genetic structure of Algerian honey bee populations

BACKGROUND

The Algerian honey bee population is composed of two described subspecies A. m. intermissa and A. m. sahariensis, of which little is known regarding population genomics, both in terms of genetic differentiation and of possible contamination by exogenous stock. Moreover, the phenotypic differences between the two subspecies are expected to translate into genetic differences and possible adaptation to heat and drought in A. m. sahariensis. To shed light on the structure of this population and to integrate these two subspecies in the growing dataset of available haploid drone sequences, we performed whole-genome sequencing of 151 haploid drones.

RESULTS

Integrated analysis of our drone sequences with a similar dataset of European reference populations did not detect any significant admixture in the Algerian honey bees. Interestingly, most of the genetic variation was not found between the A. m. intermissa and A. m. sahariensis subspecies; instead, two main genetic clusters were found along an East-West axis. We found that the correlation between genetic and geographic distances was higher in the Western cluster and that close-family relationships were mostly detected in the Eastern cluster, sometimes at long distances. In addition, we selected a panel of 96 ancestry-informative markers to decide whether a sampled bee is Algerian or not, and tested this panel in simulated cases of admixture.

CONCLUSIONS

The differences between the two main genetic clusters suggest differential breeding management between eastern and western Algeria, with greater exchange of genetic material over long distances in the east. The lack of detected admixture events suggests that, unlike what is seen in many places worldwide, imports of queens from foreign countries do not seem to have occurred on a large scale in Algeria, a finding that is relevant for conservation purposes. In addition, the proposed panel of 96 markers was found effective to distinguish Algerian from European honey bees. Therefore, we conclude that applying this approach to other taxa is promising, in particular when genetic differentiation is difficult to capture.

AmelHap: Leveraging drone whole-genome sequence data to create a honey bee HapMap

Honey bee, Apis mellifera, drones are typically haploid, developing from an unfertilized egg, inheriting only their queen's alleles and none from the many drones she mated with. Thus the ordered combination or 'phase' of alleles is known, making drones a valuable haplotype resource. We collated whole-genome sequence data for 1,407 drones, including 45 newly sequenced Scottish drones, collectively representing 19 countries, 8 subspecies and various hybrids. Following alignment to Amel_HAv3.1, variant calling and quality filtering, we retained 17.4 M high quality variants across 1,328 samples with a genotyping rate of 98.7%. We demonstrate the utility of this haplotype resource, AmelHap, for genotype imputation, returning >95% concordance when up to 61% of data is missing in haploids and up to 12% of data is missing in diploids. AmelHap will serve as a useful resource for the community for imputation from low-depth sequencing or SNP chip data, accurate phasing of diploids for association studies, and as a comprehensive reference panel for population genetic and evolutionary analyses.

Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole-genome sequencing data

Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole-genome queen genotypes, from a pool-seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well-defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F ROH) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F ROH (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F ROH (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F ROH of 2.33%. Within the breed-specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh-BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.

Reconstructing queen genotypes by pool sequencing colonies in eusocial insects: Statistical Methods and their application to honeybee

Eusocial insects are crucial to many ecosystems, and particularly the honeybee (Apis mellifera). One approach to facilitate their study in molecular genetics, is to consider whole-colony genotyping by combining DNA of multiple individuals in a single pool sequencing experiment. Cheap and fast, this technique comes with the drawback of producing data requiring dedicated methods to be fully exploited. Despite this limitation, pool sequencing data have been shown to be informative and cost-effective when working on random mating populations. Here, we present new statistical methods for exploiting pool sequencing of eusocial colonies in order to reconstruct the genotypes of the queen of such colony. This leverages the possibility to monitor genetic diversity, perform genomic-based studies or implement selective breeding. Using simulations and honeybee real data, we show that the new methods allow for a fast and accurate estimation of the queen's genetic ancestry, with correlations of about 0.9 to that obtained from individual genotyping. Also, it allows for an accurate reconstruction of the queen genotypes, with about 2% genotyping error. We further validate these inferences using experimental data on colonies with both pool sequencing and individual genotyping of drones. In brief, in this study we present statistical models to accurately estimate the genetic ancestry and reconstruct the genotypes of the queen from pool sequencing data from workers of an eusocial colony. Such information allows to exploit pool sequencing for traditional population genetics analyses, association studies and for selective breeding. While validated in Apis mellifera, these methods are applicable to other eusocial hymenopterans.

Complex population structure and haplotype patterns in the Western European honey bee from sequencing a large panel of haploid drones

Honey bee subspecies originate from specific geographical areas in Africa, Europe and the Middle East, and beekeepers interested in specific phenotypes have imported genetic material to regions outside of the bees' original range for use either in pure lines or controlled crosses. Moreover, imported drones are present in the environment and mate naturally with queens from the local subspecies. The resulting admixture complicates population genetics analyses, and population stratification can be a major problem for association studies. To better understand Western European honey bee populations, we produced a whole genome sequence and single nucleotide polymorphism (SNP) genotype data set from 870 haploid drones and demonstrate its utility for the identification of nine genetic backgrounds and various degrees of admixture in a subset of 629 samples. Five backgrounds identified correspond to subspecies, two to isolated populations on islands and two to managed populations. We also highlight several large haplotype blocks, some of which coincide with the position of centromeres. The largest is 3.6 Mb long and represents 21% of chromosome 11, with two major haplotypes corresponding to the two dominant genetic backgrounds identified. This large naturally phased data set is available as a single vcf file that can now serve as a reference for subsequent populations genomics studies in the honey bee, such as (i) selecting individuals of verified homogeneous genetic backgrounds as references, (ii) imputing genotypes from a lower-density data set generated by an SNP-chip or by low-pass sequencing, or (iii) selecting SNPs compatible with the requirements of genotyping chips.

Two quantitative trait loci are associated with recapping of Varroa destructor-infested brood cells in Apis mellifera mellifera

Recapping of Varroa destructor-infested brood cells is a trait that has recently attracted interest in honey bee breeding to select mite-resistant Apis mellifera colonies. To investigate the genetic architecture of this trait, we evaluated a sample of A. mellifera mellifera colonies (N = 155) from Switzerland and France and performed a genome-wide association study, using a pool of 500 workers per colony for next-generation sequencing. The results revealed that two QTL were significantly (P < 0.05) associated with recapping of V. destructor-infested brood cells. The best-associated QTL is located on chromosome 5 in a region previously found to be associated with grooming behaviour, a resistance trait against V. destructor, in A. mellifera and Apis cerana. The second best-associated QTL is located on chromosome 4 in an intron of the Dscam gene, which is involved in neuronal wiring. Previous research demonstrated that genes involved in neuronal wiring are associated with recapping and varroa sensitive hygiene. Therefore, our study confirms the role of a gene region on chromosome 5 in social immunity and simultaneously provides novel insights into genetic interactions between common mite resistance traits in honey bees.

Identification of quantitative trait loci associated with calmness and gentleness in honey bees using whole-genome sequences

The identification of quantitative trait loci (QTL) through genome-wide association studies (GWAS) is a powerful method for unravelling the genetic background of selected traits and improving early-stage predictions. In honey bees (Apis mellifera), past genetic analyses have particularly focused on individual queens and workers. In this study, we used pooled whole-genome sequences to ascertain the genetic variation of the entire colony. In total, we sampled 216 Apis mellifera mellifera and 28 Apis mellifera carnica colonies. Different experts subjectively assessed the gentleness and calmness of the colonies using a standardised protocol. Conducting a GWAS for calmness on 211 purebred A. m. mellifera colonies, we identified three QTL, on chromosomes 8, 6, and 12. The two first QTL correspond to LOC409692 gene, coding for a disintegrin and metalloproteinase domain-containing protein 10, and to Abscam gene, coding for a Dscam family member Abscam protein, respectively. The last gene has been reported to be involved in the domestication of A. mellifera. The third QTL is located 13 kb upstream of LOC102655631, coding for a trehalose transporter. For gentleness, two QTL were identified on chromosomes 4 and 3. They are located within gene LOC413669, coding for a lap4 protein, and gene LOC413416, coding for a bicaudal C homolog 1-B protein, respectively. The identified positional candidate genes of both traits mainly affect the olfaction and nervous system of honey bees. Further research is needed to confirm the results and to better understand the genetic and phenotypic basis of calmness and gentleness.

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.

Descriptive Analysis of the Varroa Non-Reproduction Trait in Honey Bee Colonies and Association with Other Traits Related to Varroa Resistance

In the current context of worldwide honey bee colony losses, among which the varroa mite plays a major role, the hope to improve honey bee health lies in part in the breeding of varroa resistant colonies. To do so, methods used to evaluate varroa resistance need better understanding. Repeatability and correlations between traits such as mite non-reproduction (MNR), varroa sensitive hygiene (VSH), and hygienic behavior are poorly known, due to practical limitations and to their underlying complexity. We investigate (i) the variability, (ii) the repeatability of the MNR score, and (iii) its correlation with other resistance traits. To reduce the inherent variability of MNR scores, we propose to apply an empirical Bayes correction. In the short-term (ten days), MNR had a modest repeatability of 0.4, whereas in the long-term (a month), it had a low repeatability of 0.2, similar to other resistance traits. Within our dataset, there was no correlation between MNR and VSH. Although MNR is amongst the most popular varroa resistance estimates in field studies, its underlying complex mechanism is not fully understood. Its lack of correlation with better described resistance traits and low repeatability suggest that MNR needs to be interpreted cautiously, especially when used for selection.

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.

Structure of the intergenic spacers in chicken ribosomal DNA

BACKGROUND

Ribosomal DNA (rDNA) repeats are situated in the nucleolus organizer regions (NOR) of chromosomes and transcribed into rRNA for ribosome biogenesis. Thus, they are an essential component of eukaryotic genomes. rDNA repeat units consist of rRNA gene clusters that are transcribed into single pre-rRNA molecules, each separated by intergenic spacers (IGS) that contain regulatory elements for rRNA gene cluster transcription. Because of their high repeat content, rDNA sequences are usually absent from genome assemblies. In this work, we used the long-read sequencing technology to describe the chicken IGS and fill the knowledge gap on rDNA sequences of one of the key domesticated animals.

METHODS

We used the long-read PacBio RSII technique to sequence the BAC clone WAG137G04 (Wageningen BAC library) known to contain chicken NOR elements and the HGAP workflow software suit to assemble the PacBio RSII reads. Whole-genome sequence contigs homologous to the chicken rDNA repetitive unit were identified based on the Gallus_gallus-5.0 assembly with BLAST. We used the Geneious 9.0.5 and Mega software, maximum likelihood method and Chickspress project for sequence evolution analysis, phylogenetic tree construction and analysis of the raw transcriptome data.

RESULTS

Three complete IGS sequences in the White Leghorn chicken genome and one IGS sequence in the red junglefowl contig AADN04001305.1 (Gallus_gallus-5.0) were detected. They had various lengths and contained three groups of tandem repeats (some of them being very GC rich) that form highly organized arrays. Initiation and termination sites of rDNA transcription were located within small and large unique regions (SUR and LUR), respectively. No functionally significant sites were detected within the tandem repeat sequences.

CONCLUSIONS

Due to the highly organized GC-rich repeats, the structure of the chicken IGS differs from that of IGS in human, apes, Xenopus or fish rDNA. However, the chicken IGS shares some molecular organization features with that of the turtles, which are other representatives of the Sauropsida clade that includes birds and reptiles. Our current results on the structure of chicken IGS together with the previously reported ribosomal gene cluster sequence provide sufficient data to consider that the complete chicken rDNA sequence is assembled with confidence in terms of molecular DNA organization.

A guinea fowl genome assembly provides new evidence on evolution following domestication and selection in galliformes

The helmeted guinea fowl Numida meleagris belongs to the order Galliformes. Its natural range includes a large part of sub‐Saharan Africa, from Senegal to Eritrea and from Chad to South Africa. Archaeozoological and artistic evidence suggest domestication of this species may have occurred about 2,000 years BP in Mali and Sudan primarily as a food resource, although villagers also benefit from its capacity to give loud alarm calls in case of danger, of its ability to consume parasites such as ticks and to hunt snakes, thus suggesting its domestication may have resulted from a commensal association process. Today, it is still farmed in Africa, mainly as a traditional village poultry, and is also bred more intensively in other countries, mainly France and Italy. The lack of available molecular genetic markers has limited the genetic studies conducted to date on guinea fowl. We present here a first‐generation whole‐genome sequence draft assembly used as a reference for a study by a Pool‐seq approach of wild and domestic populations from Europe and Africa. We show that the domestic populations share a higher genetic similarity between each other than they do to wild populations living in the same geographical area. Several genomic regions showing selection signatures putatively related to domestication or importation to Europe were detected, containing candidate genes, most notably EDNRB2, possibly explaining losses in plumage coloration phenotypes in domesticated populations.

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.

Autosomal and Mitochondrial Adaptation Following Admixture: A Case Study on the Honeybees of Reunion Island

The honeybee population of the tropical Reunion Island is a genetic admixture of the Apis mellifera unicolor subspecies, originally described in Madagascar, and of European subspecies, mainly A. m. carnica and A. m. ligustica, regularly imported to the island since the late 19th century. We took advantage of this population to study genetic admixing of the tropical-adapted indigenous and temperate-adapted European genetic backgrounds. Whole genome sequencing of 30 workers and 6 males from Reunion, compared with samples from Europe, Madagascar, Mauritius, Rodrigues, and the Seychelles, revealed the Reunion honeybee population to be composed on an average of 53.2 ± 5.9% A. m. unicolor nuclear genomic background, the rest being mainly composed of A. m. carnica and to a lesser extent A. m. ligustica. In striking contrast to this, only 1 out of the 36 honeybees from Reunion had a mitochondrial genome of European origin, suggesting selection has favored the A. m. unicolor mitotype, which is possibly better adapted to the island’s bioclimate. Local ancestry was determined along the chromosomes for all Reunion samples, and a test for preferential selection for the A. m. unicolor or European background revealed 15 regions significantly associated with the A. m. unicolor lineage and 9 regions with the European lineage. Our results provide insights into the long-term consequences of introducing exotic specimen on the nuclear and mitochondrial genomes of locally adapted populations.

Developing reduced SNP assays from whole-genome sequence data to estimate introgression in an organism with complex genetic patterns, the Iberian honeybee (Apis mellifera iberiensis)

The most important managed pollinator, the honeybee (Apis mellifera L.), has been subject to a growing number of threats. In western Europe, one such threat is large-scale introductions of commercial strains (C-lineage ancestry), which is leading to introgressive hybridization and even the local extinction of native honeybee populations (M-lineage ancestry). Here, we developed reduced assays of highly informative SNPs from 176 whole genomes to estimate C-lineage introgression in the most diverse and evolutionarily complex subspecies in Europe, the Iberian honeybee (Apis mellifera iberiensis). We started by evaluating the effects of sample size and sampling a geographically restricted area on the number of highly informative SNPs. We demonstrated that a bias in the number of fixed SNPs (FST = 1) is introduced when the sample size is small (N ≤ 10) and when sampling only captures a small fraction of a population's genetic diversity. These results underscore the importance of having a representative sample when developing reliable reduced SNP assays for organisms with complex genetic patterns. We used a training data set to design four independent SNP assays selected from pairwise FST between the Iberian and C-lineage honeybees. The designed assays, which were validated in holdout and simulated hybrid data sets, proved to be highly accurate and can be readily used for monitoring populations not only in the native range of A. m. iberiensis in Iberia but also in the introduced range in the Balearic islands, Macaronesia and South America, in a time- and cost-effective manner. While our approach used the Iberian honeybee as model system, it has a high value in a wide range of scenarios for the monitoring and conservation of potentially hybridized domestic and wildlife populations.

Genome-wide scans between two honeybee populations reveal putative signatures of human-mediated selection

Human‐mediated selection has left signatures in the genomes of many domesticated animals, including the European dark honeybee, Apis mellifera mellifera, which has been selected by apiculturists for centuries. Using whole‐genome sequence information, we investigated selection signatures in spatially separated honeybee subpopulations (Switzerland, n = 39 and France, n = 17). Three different test statistics were calculated in windows of 2 kb (fixation index, cross‐population extended haplotype homozygosity and cross‐population composite likelihood ratio) and combined into a recently developed composite selection score. Applying a stringent false discovery rate of 0.01, we identified six significant selective sweeps distributed across five chromosomes covering eight genes. These genes are associated with multiple molecular and biological functions, including regulation of transcription, receptor binding and signal transduction. Of particular interest is a selection signature on chromosome 1, which corresponds to the WNT4 gene, the family of which is conserved across the animal kingdom with a variety of functions. In Drosophila melanogaster, WNT4 alleles have been associated with differential wing, cross vein and abdominal phenotypes. Defining phenotypic characteristics of different Apis mellifera ssp., which are typically used as selection criteria, include colour and wing venation pattern. This signal is therefore likely to be a good candidate for human mediated‐selection arising from different applied breeding practices in the two managed populations.

GWAS analyses reveal QTL in egg layers that differ in response to diet differences

BACKGROUND

The genetic architecture of egg production and egg quality traits, i.e. the quantitative trait loci (QTL) that influence these traits, is still poorly known. To date, 33 studies have focused on the detection of QTL for laying traits in chickens, but less than 10 genes have been identified. The availability of a high-density SNP (single nucleotide polymorphism) chicken array developed by Affymetrix, i.e. the 600K Affymetrix(®) Axiom(®) HD genotyping array offers the possibility to narrow down the localization of previously detected QTL and to detect new QTL. This high-density array is also anticipated to take research beyond the classical hypothesis of additivity of QTL effects or of QTL and environmental effects. The aim of our study was to search for QTL that influence laying traits using the 600K SNP chip and to investigate whether the effects of these QTL differed between diets and age at egg collection.

RESULTS

One hundred and thirty-one QTL were detected for 16 laying traits and were spread across all marked chromosomes, except chromosomes 16 and 25. The percentage of variance explained by a QTL varied from 2 to 10 % for the various traits, depending on diet and age at egg collection. Chromosomes 3, 9, 10 and Z were overrepresented, with more than eight QTL on each one. Among the 131 QTL, 60 had a significantly different effect, depending on diet or age at egg collection. For egg production traits, when the QTL × environment interaction was significant, numerous inversions of sign of the SNP effects were observed, whereas for egg quality traits, the QTL × environment interaction was mostly due to a difference of magnitude of the SNP effects.

CONCLUSIONS

Our results show that numerous QTL influence egg production and egg quality traits and that the genomic regions, which are involved in shaping the ability of layer chickens to adapt to their environment for egg production, vary depending on the environmental conditions. The next question will be to address what the impact of these genotype × environment interactions is on selection.

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.

Mule duck "foie gras" shows different metabolic states according to its quality phenotype by using a proteomic approach

This study aimed at identifying the mechanisms implicated in "foie gras" quality variability through the study of the relationships between liver protein compositions and four liver quality phenotypes: liver weight, melting rate, and protein contents on crude or dry matter. Spots of soluble proteins were separated by bidimensional electrophoresis, and the relative abundance of proteins according to quality traits values was investigated. Twenty-three protein spots (19 unique identified proteins) showed different levels of abundance according to one or more of the traits' values. These abundance differences highlighted two groups of livers with opposite trends of abundance levels. Proteins of the first group, associated with low liver weight and melting rate, are involved in synthesis and anabolism processes, whereas proteins of the second group, associated with high liver weight and melting rate, are proteins involved in stress response. Altogether, these results highlight the variations in metabolic states underlying foie gras quality traits.

Mapping and genotypic analysis of the NK-lysin gene in chicken

Background

Antimicrobial peptides (AMP) are important elements of the first line of defence against pathogens in animals. NK-lysin is a cationic AMP that plays a critical role in innate immunity. The chicken NK-lysin gene has been cloned and its antimicrobial and anticancer activity has been described but its location in the chicken genome remains unknown. Here, we mapped the NK-lysin gene and examined the distribution of a functionally significant single nucleotide polymorphism (SNP) among different chicken inbred lines and heritage breeds.

Results

A 6000 rad radiation hybrid panel (ChickRH6) was used to map the NK-lysin gene to the distal end of chromosome 22. Two additional genes, the adipocyte enhancer-binding protein 1-like gene (AEBP1) and the DNA polymerase delta subunit 2-like (POLD2) gene, are located in the same NW_003779909 contig as NK-lysin, and were thus indirectly mapped to chromosome 22 as well. Previously, we reported a functionally significant SNP at position 271 of the NK-lysin coding sequence in two different chicken breeds. Here, we examined this SNP and found that the A allele appears to be more common than the G allele in these heritage breeds and inbred lines.

Conclusions

The chicken NK-lysin gene mapped to the distal end of chromosome 22. Two additional genes, AEBP1 and POLD2, were indirectly mapped to chromosome 22 also. SNP analyses revealed that the A allele, which encodes a peptide with a higher antimicrobial activity, is more common than the G allele in our tested inbred lines and heritage breeds.

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.

The duck genome and transcriptome provide insight into an avian influenza virus reservoir species

The duck (Anas platyrhynchos) is one of the principal natural hosts of influenza A viruses. We present the duck genome sequence and perform deep transcriptome analyses to investigate immune-related genes. Our data indicate that the duck possesses a contractive immune gene repertoire, as in chicken and zebra finch, and this repertoire has been shaped through lineage-specific duplications. We identify genes that are responsive to influenza A viruses using the lung transcriptomes of control ducks and ones that were infected with either a highly pathogenic (A/duck/Hubei/49/05) or a weakly pathogenic (A/goose/Hubei/65/05) H5N1 virus. Further, we show how the duck's defense mechanisms against influenza infection have been optimized through the diversification of its β-defensin and butyrophilin-like repertoires. These analyses, in combination with the genomic and transcriptomic data, provide a resource for characterizing the interaction between host and influenza viruses.

Detection of QTL controlling metabolism, meat quality, and liver quality traits of the overfed interspecific hybrid mule duck

The mule duck, an interspecific hybrid obtained by crossing common duck (Anas platyrhynchos) females with Muscovy (Cairina moschata) drakes, is widely used for fatty liver production. The purpose of the present study was to detect and map single and pleiotropic QTL that segregate in the common duck species, and influence the expression of traits in their overfed mule duck offspring. To this end, we generated a common duck backcross (BC) population by crossing Kaiya and heavy Pekin experimental lines, which differ notably in regard to the BW and overfeeding ability of their mule progeny. The BC females were mated to Muscovy drakes and, on average, 4 male mule ducks hatched per BC female (1600 in total) and were measured for growth, metabolism during growth and the overfeeding period, overfeeding ability, and the quality of their breast meat and fatty liver. The phenotypic value of BC females was estimated for each trait by assigning to each female the mean value of the phenotypes of her offspring. Estimations allowed for variance, which depended on the number of male offspring per BC and the heritability of the trait considered. The genetic map used for QTL detection consisted of 91 microsatellite markers aggregated into 16 linkage groups (LG) covering a total of 778 cM. Twenty-two QTL were found to be significant at the 1% chromosome-wide threshold level using the single-trait detection option of the QTLMap software. Most of the QTL detected were related to the quality of breast meat and fatty liver: QTL for meat pH 20 min post mortem were mapped to LG4 (at the 1% genome-wide significance level), and QTL for meat lipid content and cooking losses were mapped to LG2a. The QTL related to fatty liver weight and liver protein and lipid content were for the most part detected on LG2c and LG9. Multitrait analysis highlighted the pleiotropic effects of QTL in these chromosome regions. Apart from the strong QTL for plasma triglyceride content at the end of the overfeeding period mapped to chromosome Z using single-trait analysis, all metabolic trait QTL were detected with the multitrait approach: the QTL mapped to LG14 and LG21 affected the plasma cholesterol and triglyceride contents, whereas the QTL mapped to LG2a seemed to impact glycemia and the basal plasma corticosterone content. A greater density genetic map will be needed to further fine map the QTL.

Fine mapping of complex traits in non-model species: using next generation sequencing and advanced intercross lines in Japanese quail

Background

As for other non-model species, genetic analyses in quail will benefit greatly from a higher marker density, now attainable thanks to the evolution of sequencing and genotyping technologies. Our objective was to obtain the first genome wide panel of Japanese quail SNP (Single Nucleotide Polymorphism) and to use it for the fine mapping of a QTL for a fear-related behaviour, namely tonic immobility, previously localized on Coturnix japonica chromosome 1. To this aim, two reduced representations of the genome were analysed through high-throughput 454 sequencing: AFLP (Amplified Fragment Length Polymorphism) fragments as representatives of genomic DNA, and EST (Expressed Sequence Tag) as representatives of the transcriptome.

Results

The sequencing runs produced 399,189 and 1,106,762 sequence reads from cDNA and genomic fragments, respectively. They covered over 434 Mb of sequence in total and allowed us to detect 17,433 putative SNP. Among them, 384 were used to genotype two Advanced Intercross Lines (AIL) obtained from three quail lines differing for duration of tonic immobility. Despite the absence of genotyping for founder individuals in the analysis, the previously identified candidate region on chromosome 1 was refined and led to the identification of a candidate gene.

Conclusions

These data confirm the efficiency of transcript and AFLP-sequencing for SNP discovery in a non-model species, and its application to the fine mapping of a complex trait. Our results reveal a significant association of duration of tonic immobility with a genomic region comprising the DMD (dystrophin) gene. Further characterization of this candidate gene is needed to decipher its putative role in tonic immobility in Coturnix.

A duck RH panel and its potential for assisting NGS genome assembly

Background

Owing to the low cost of the high throughput Next Generation Sequencing (NGS) technology, more and more species have been and will be sequenced. However, de novo assemblies of large eukaryotic genomes thus produced are composed of a large number of contigs and scaffolds of medium to small size, having no chromosomal assignment. Radiation hybrid (RH) mapping is a powerful tool for building whole genome maps and has been used for several animal species, to help assign sequence scaffolds to chromosomes and determining their order.

Results

We report here a duck whole genome RH panel obtained by fusing female duck embryonic fibroblasts irradiated at a dose of 6,000 rads, with HPRT-deficient Wg3hCl₂ hamster cells. The ninety best hybrids, having an average retention of 23.6% of the duck genome, were selected for the final panel. To allow the genotyping of large numbers of markers, as required for whole genome mapping, without having to cultivate the hybrid clones on a large scale, three different methods involving Whole Genome Amplification (WGA) and/or scaling down PCR volumes by using the Fluidigm BioMark^TM Integrated Fluidic Circuits (IFC) Dynamic Array^TM for genotyping were tested. RH maps of APL12 and APL22 were built, allowing the detection of intrachromosomal rearrangements when compared to chicken. Finally, the panel proved useful for checking the assembly of sequence scaffolds and for mapping EST located on one of the smallest microchromosomes.

Conclusion

The Fluidigm BioMark^TM Integrated Fluidic Circuits (IFC) Dynamic Array^TM genotyping by quantitative PCR provides a rapid and cost-effective method for building RH linkage groups. Although the vast majority of genotyped markers exhibited a picture coherent with their associated scaffolds, a few of them were discordant, pinpointing potential assembly errors. Comparative mapping with chicken chromosomes GGA21 and GGA11 allowed the detection of the first chromosome rearrangements on microchromosomes between duck and chicken. As in chicken, the smallest duck microchromosomes appear missing in the assembly and more EST data will be needed for mapping them. Altogether, this underlines the added value of RH mapping to improve genome assemblies.

Epilepsy caused by an abnormal alternative splicing with dosage effect of the SV2A gene in a chicken model

Photosensitive reflex epilepsy is caused by the combination of an individual's enhanced sensitivity with relevant light stimuli, such as stroboscopic lights or video games. This is the most common reflex epilepsy in humans; it is characterized by the photoparoxysmal response, which is an abnormal electroencephalographic reaction, and seizures triggered by intermittent light stimulation. Here, by using genetic mapping, sequencing and functional analyses, we report that a mutation in the acceptor site of the second intron of SV2A (the gene encoding synaptic vesicle glycoprotein 2A) is causing photosensitive reflex epilepsy in a unique vertebrate model, the Fepi chicken strain, a spontaneous model where the neurological disorder is inherited as an autosomal recessive mutation. This mutation causes an aberrant splicing event and significantly reduces the level of SV2A mRNA in homozygous carriers. Levetiracetam, a second generation antiepileptic drug, is known to bind SV2A, and SV2A knock-out mice develop seizures soon after birth and usually die within three weeks. The Fepi chicken survives to adulthood and responds to levetiracetam, suggesting that the low-level expression of SV2A in these animals is sufficient to allow survival, but does not protect against seizures. Thus, the Fepi chicken model shows that the role of the SV2A pathway in the brain is conserved between birds and mammals, in spite of a large phylogenetic distance. The Fepi model appears particularly useful for further studies of physiopathology of reflex epilepsy, in comparison with induced models of epilepsy in rodents. Consequently, SV2A is a very attractive candidate gene for analysis in the context of both mono- and polygenic generalized epilepsies in humans.

Integrative mapping analysis of chicken microchromosome 16 organization

BACKGROUND

The chicken karyotype is composed of 39 chromosome pairs, of which 9 still remain totally absent from the current genome sequence assembly, despite international efforts towards complete coverage. Some others are only very partially sequenced, amongst which microchromosome 16 (GGA16), particularly under-represented, with only 433 kb assembled for a full estimated size of 9 to 11 Mb. Besides the obvious need of full genome coverage with genetic markers for QTL (Quantitative Trait Loci) mapping and major genes identification studies, there is a major interest in the detailed study of this chromosome because it carries the two genetically independent MHC complexes B and Y. In addition, GGA16 carries the ribosomal RNA (rRNA) genes cluster, also known as the NOR (nucleolus organizer region). The purpose of the present study is to construct and present high resolution integrated maps of GGA16 to refine its organization and improve its coverage with genetic markers.

RESULTS

We developed 79 STS (Sequence Tagged Site) markers to build a physical RH (radiation hybrid) map and 34 genetic markers to extend the genetic map of GGA16. We screened a BAC (Bacterial Artificial Chromosome) library with markers for the MHC-B, MHC-Y and rRNA complexes. Selected clones were used to perform high resolution FISH (Fluorescent In Situ Hybridization) mapping on giant meiotic lampbrush chromosomes, allowing meiotic mapping in addition to the confirmation of the order of the three clusters along the chromosome. A region with high recombination rates and containing PO41 repeated elements separates the two MHC complexes.

CONCLUSIONS

The three complementary mapping strategies used refine greatly our knowledge of chicken microchromosome 16 organisation. The characterisation of the recombination hotspots separating the two MHC complexes demonstrates the presence of PO41 repetitive sequences both in tandem and inverted orientation. However, this region still needs to be studied in more detail.

Non PCR-amplified Transcripts and AFLP fragments as reduced representations of the quail genome for 454 Titanium sequencing

BACKGROUND

SNP (Single Nucleotide Polymorphism) discovery is now routinely performed using high-throughput sequencing of reduced representation libraries. Our objective was to adapt 454 GS FLX based sequencing methodologies in order to obtain the largest possible dataset from two reduced representations libraries, produced by AFLP (Amplified Fragment Length Polymorphism) for genomic DNA, and EST (Expressed Sequence Tag) for the transcribed fraction of the genome.

FINDINGS

The expressed fraction was obtained by preparing cDNA libraries without PCR amplification from quail embryo and brain. To optimize the information content for SNP analyses, libraries were prepared from individuals selected in three quail lines and each individual in the AFLP library was tagged. Sequencing runs produced 399,189 sequence reads from cDNA and 373,484 from genomic fragments, covering close to 250 Mb of sequence in total.

CONCLUSIONS

Both methods used to obtain reduced representations for high-throughput sequencing were successful after several improvements.The protocols may be used for several sequencing applications, such as de novo sequencing, tagged PCR fragments or long fragment sequencing of cDNA.

Genetic control of resistance to salmonellosis and to Salmonella carrier-state in fowl: a review

Salmonellosis is a frequent disease in poultry stocks, caused by several serotypes of the bacterial species Salmonella enterica and sometimes transmitted to humans through the consumption of contaminated meat or eggs. Symptom-free carriers of the bacteria contribute greatly to the propagation of the disease in poultry stocks. So far, several candidate genes and quantitative trait loci (QTL) for resistance to carrier state or to acute disease have been identified using artificial infection of S. enterica serovar Enteritidis or S. enterica serovar Typhimurium strains in diverse genetic backgrounds, with several different infection procedures and phenotypic assessment protocols. This diversity in experimental conditions has led to a complex sum of results, but allows a more complete description of the disease. Comparisons among studies show that genes controlling resistance to Salmonella differ according to the chicken line studied, the trait assessed and the chicken's age. The loci identified are located on 25 of the 38 chicken autosomal chromosomes. Some of these loci are clustered in several genomic regions, indicating the possibility of a common genetic control for different models. In particular, the genomic regions carrying the candidate genes TLR4 and SLC11A1, the Major Histocompatibility Complex (MHC) and the QTL SAL1 are interesting for more in-depth studies. This article reviews the main Salmonella infection models and chicken lines studied under a historical perspective and then the candidate genes and QTL identified so far.

Mapping main, epistatic and sex-specific QTL for body composition in a chicken population divergently selected for low or high growth rate

Background

Delineating the genetic basis of body composition is important to agriculture and medicine. In addition, the incorporation of gene-gene interactions in the statistical model provides further insight into the genetic factors that underlie body composition traits. We used Bayesian model selection to comprehensively map main, epistatic and sex-specific QTL in an F₂ reciprocal intercross between two chicken lines divergently selected for high or low growth rate.

Results

We identified 17 QTL with main effects across 13 chromosomes and several sex-specific and sex-antagonistic QTL for breast meat yield, thigh + drumstick yield and abdominal fatness. Different sets of QTL were found for both breast muscles [Pectoralis (P) major and P. minor], which suggests that they could be controlled by different regulatory mechanisms. Significant interactions of QTL by sex allowed detection of sex-specific and sex-antagonistic QTL for body composition and abdominal fat. We found several female-specific P. major QTL and sex-antagonistic P. minor and abdominal fatness QTL. Also, several QTL on different chromosomes interact with each other to affect body composition and abdominal fatness.

Conclusions

The detection of main effects, epistasis and sex-dimorphic QTL suggest complex genetic regulation of somatic growth. An understanding of such regulatory mechanisms is key to mapping specific genes that underlie QTL controlling somatic growth in an avian model.

QTL for several metabolic traits map to loci controlling growth and body composition in an F2 intercross between high- and low-growth chicken lines

Quantitative trait loci (QTL) for metabolic and body composition traits were mapped at 7 and 9 wk, respectively, in an F2 intercross between high-growth and low-growth chicken lines. These lines also diverged for abdominal fat percentage (AFP) and plasma insulin-like growth factor-I (IGF-I), insulin, and glucose levels. Genotypings were performed with 129 microsatellite markers covering 21 chromosomes. A total of 21 QTL with genomewide level of significance were detected by single-trait analyses for body weight (BW), breast muscle weight (BMW) and percentage (BMP), AF weight (AFW) and percentage (AFP), shank length (ShL) and diameter (ShD), fasting plasma glucose level (Gluc), and body temperature (Tb). Other suggestive QTL were identified for these parameters and for plasma IGF-I and nonesterified fatty acid levels. QTL controlling adiposity and Gluc were colocalized on GGA3 and GGA5 and QTL for BW, ShL and ShD, adiposity, and Tb on GGA4. Multitrait analyses revealed two QTL controlling Gluc and AFP on GGA5 and Gluc and Tb on GGA26. Significant effects of the reciprocal cross were observed on BW, ShD, BMW, and Gluc, which may result from mtDNA and/or maternal effects. Most QTL regions for Gluc and adiposity harbor genes for which alleles have been associated with increased susceptibility to diabetes and/or obesity in humans. Identification of genes responsible for these metabolic QTL will increase our understanding of the constitutive “hyperglycemia” found in chickens. Furthermore, a comparative approach could provide new information on the genetic causes of diabetes and obesity in humans.

QTL for resistance to Salmonella carrier state confirmed in both experimental and commercial chicken lines

The ability of chickens to carry Salmonella without displaying disease symptoms is responsible for Salmonella propagation in poultry stocks and for subsequent human contamination through the consumption of contaminated eggs or meat. The selection of animals more resistant to carrier state might be a way to decrease the propagation of Salmonella in poultry stocks and its transmission to humans. Five QTL controlling variation for resistance to carrier state in a chicken F(2) progeny derived from the White Leghorn inbred lines N and 6(1) had been previously identified using a selective genotyping approach. Here, a second analysis on the whole progeny was performed, which led to the confirmation of two QTL on chromosomes 2 and 16. To assess the utility of these genomic regions for selection in commercial lines, we tested them together with other QTL identified in an [Nx6(1)] x N backcross progeny and with the candidate genes SLC11A1 and TLR4. We used a commercial line divergently selected for either low or high carrier-state resistance both in young chicks and in adult hens. In divergent chick lines, one QTL on chromosome 1 and one in the SLC11A1 region were significantly associated with carrier-state resistance variations; in divergent adult lines, one QTL located in the major histocompatibility complex on chromosome 16 and one in the SLC11A1 region were involved in these variations. Genetic studies conducted on experimental lines can therefore be of potential interest for marker-assisted selection in commercial lines.

Whole genome comparative studies between chicken and turkey and their implications for avian genome evolution

Background

Comparative genomics is a powerful means of establishing inter-specific relationships between gene function/location and allows insight into genomic rearrangements, conservation and evolutionary phylogeny. The availability of the complete sequence of the chicken genome has initiated the development of detailed genomic information in other birds including turkey, an agriculturally important species where mapping has hitherto focused on linkage with limited physical information. No molecular study has yet examined conservation of avian microchromosomes, nor differences in copy number variants (CNVs) between birds.

Results

We present a detailed comparative cytogenetic map between chicken and turkey based on reciprocal chromosome painting and mapping of 338 chicken BACs to turkey metaphases. Two inter-chromosomal changes (both involving centromeres) and three pericentric inversions have been identified between chicken and turkey; and array CGH identified 16 inter-specific CNVs.

Conclusion

This is the first study to combine the modalities of zoo-FISH and array CGH between different avian species. The first insight into the conservation of microchromosomes, the first comparative cytogenetic map of any bird and the first appraisal of CNVs between birds is provided. Results suggest that avian genomes have remained relatively stable during evolution compared to mammalian equivalents.

FISH mapping of 57 BAC clones reveals strong conservation of synteny between Galliformes and Anseriformes

Karyotypes of chicken (Gallus gallus domesticus; 2n = 78) and mallard duck (Anas platyrhynchos; 2n = 80) share the typical organization of avian karyotypes including a few macrochromosome pairs, numerous indistinguishable microchromosomes, and Z and W sex chromosomes. Previous banding studies revealed great similarities between chickens and ducks, but it was not possible to use comparative banding for the microchromosomes. In order to establish precise chromosome correspondences between these two species, particularly for microchromosomes, we hybridized 57 BAC clones previously assigned to the chicken genome to duck metaphase spreads. Although most of the clones showed similar localizations, we found a few intrachromosomal rearrangements of the macrochromosomes and an additional microchromosome pair in ducks. BAC clones specific for chicken microchromosomes were localized to separate duck microchromosomes and clones mapping to the same chicken microchromosome hybridized to the same duck microchromosome, demonstrating a high conservation of synteny. These results demonstrate that the evolution of karyotypes in avian species is the result of fusion and/or fission processes and not translocations.

An integrated approach of genetic resistance to Salmonella carrier state in fowls: from genetics to genomics and modelling

Increasing resistance to acute Salmonellosis (that is, contamination level shortly after infection) is not sufficient to reduce the risk for consumers to be contaminated by Salmonella. Indeed, animals may remain contaminated at a low level for weeks or months. Increased resistance to the Salmonella carrier state, i.e., animals' ability to clear bacteria, is needed; it involves measuring bacterial contamination several weeks after inoculation with a low dose. To study such resistance traits, three convergent approaches were used. A quantitative trait loci (QTL) study was performed, taking advantage of inbred lines differing in resistance. Several QTLs controlling resistance at a younger age were identified and are currently being confirmed in a new cross before finer mapping, using advanced intercross lines. These inbred lines are also presently being compared using functional genomics. In parallel, a selection experiment for increased or decreased resistance at a younger and a later age was undertaken. Besides providing genetic models differing in their levels of resistance, it underlined the importance of the choice of selection criterion, whether marker assisted or not. Indeed, genes controlling resistance are strongly dependant on age; selecting for resistance at a younger age might result in increased susceptibility at an older age. Finally, the results of this experiment were used in a model of the intra-flock propagation of Salmonella. It showed that introducing a proportion of resistant animals within a flock of susceptible hens could dramatically change the evolution of contamination. Moreover, it demonstrated the magnitude of synergy between selection and vaccination, which should enhance the interest of increased resistance. The results show that selection for increased resistance to the Salmonella carrier state may be efficient, providing that the appropriate criteria of selection are used.

The chicken RH map: current state of progress and microchromosome mapping

The ChickRH6 radiation hybrid panel has been used to construct consensus chromosome radiation hybrid (RH) maps of the chicken genome. Markers genotyped were either from throughout the genome or targeted to specific chromosomes and a large proportion (one third) of data was the result of collaborative efforts. Altogether, 2,531 markers were genotyped, allowing the construction of RH reference maps for 20 chromosomes and linkage groups for four other chromosomes. Amongst the markers, 581 belong to the framework maps, while 1,721 are on the comprehensive maps. Around 800 markers still have to be assigned to linkage groups. Our attempt to assign the supercontigs from the chrun (virtual chromosome containing all the genome sequence that could not be attributed to a chromosome) as well as EST (Expressed Sequence Tag) contigs that do not have a BLAST hit in the genome assembly led to the construction of new maps for microchromosomes either absent or for which very little data is present in the genome assembly. RH data is presented through our ChickRH webserver (http://chickrh.toulouse.inra.fr/), which is a mapping tool as well as the official repository RH database for genotypes. It also displays the RH reference maps and comparison charts with the sequence thus highlighting the possible discrepancies. Future improvements of the RH maps include complete coverage of the sequence assigned to chromosomes, further mapping of the chrun and mapping of EST contigs absent from the assembly. This will help finish the mapping of the smallest gene-rich microchromosomes.

Female-specific DNA sequences in the chicken genome

Eight in silico W-specific sequences from the WASHUC1 chicken genome assembly gave female-specific PCR products using chicken DNA. Some of these fragments gave female-specific products with turkey and peacock DNA. Sequence analysis of these 8 fragments (3077 bp total) failed to detect any polymorphisms among 10 divergent chickens. In contrast, comparison of the DNA sequences of chicken with those of turkey and peacock revealed a nucleotide difference every 25 and 28 bp, respectively. Radiation hybrid mapping verified that these amplicons exist only on chromosome W. The homology of 6 W-specific fragments with chromo-helicase-DNA-binding gene and expressed sequenced tags from chicken and other species indicate that these fragments may have or have had a biological function. These fragments may be used for early sexing in commercial chicken and turkey flocks.

Cloning of Ovocalyxin-36, a Novel Chicken Eggshell Protein Related to Lipopolysaccharide-binding Proteins, Bactericidal Permeability-increasing Proteins, and Plunc Family Proteins

The avian eggshell is a composite biomaterial composed of noncalcifying eggshell membranes and the overlying calcified shell matrix. The shell is deposited in a uterine fluid where the concentration of different protein species varies at different stages of its formation. The role of avian eggshell proteins during shell formation remains poorly understood, and we have sought to identify and characterize the individual components in order to gain insight into their function during elaboration of the eggshell. In this study, we have used direct sequencing, immunochemistry, expression screening, and EST data base mining to clone and characterize a 1995-bp full-length cDNA sequence corresponding to a novel chicken eggshell protein that we have named Ovocalyxin-36 (OCX-36). Ovocalyxin-36 protein was only detected in the regions of the oviduct where egg-shell formation takes place; uterine OCX-36 message was strongly up-regulated during eggshell calcification. OCX-36 localized to the calcified eggshell predominantly in the inner part of the shell, and to the shell membranes. BlastN data base searching indicates that there is no mammalian version of OCX-36; however, the protein sequence is 20-25% homologous to proteins associated with the innate immune response as follows: lipopolysaccharide-binding proteins, bactericidal permeability-increasing proteins, and Plunc family proteins. Moreover, the genomic organization of these proteins and OCX-36 appears to be highly conserved. These observations suggest that OCX-36 is a novel and specific chicken eggshell protein related to the superfamily of lipopolysaccharide-binding proteins/bactericidal permeability-increasing proteins and Plunc proteins. OCX-36 may therefore participate in natural defense mechanisms that keep the egg free of pathogens.

FISH on avian lampbrush chromosomes produces higher resolution gene mapping

Giant lampbrush chromosomes, which are characteristic of the diplotene stage of prophase I during avian oogenesis, represent a very promising system for precise physical gene mapping. We applied 35 chicken BAC and 4 PAC clones to both mitotic metaphase chromosomes and meiotic lampbrush chromosomes of chicken (Gallus gallus domesticus) and Japanese quail (Coturnix coturnix japonica). Fluorescence in situ hybridization (FISH) mapping on lampbrush chromosomes allowed us to distinguish closely located probes and revealed gene order more precisely. Our data extended the data earlier obtained using FISH to chicken and quail metaphase chromosomes 1-6 and Z. Extremely low levels of inter- and intra-chromosomal rearrangements in the chicken and Japanese quail were demonstrated again. Moreover, we did not confirm the presence of a pericentric inversion in Japanese quail chromosome 4 as compared to chicken chromosome 4. Twelve BAC clones specific for chicken chromosome 4p and 4q showed the same order in quail as in chicken when FISH was performed on lampbrush chromosomes. The centromeres of chicken and quail chromosomes 4 seem to have formed independently after centric fusion of ancestral chromosome 4 and a microchromosome.

Search for QTL affecting the shape of the egg laying curve of the Japanese quail

Background

Egg production is of critical importance in birds not only for their reproduction but also for human consumption as the egg is a highly nutritive and balanced food. Consequently, laying in poultry has been improved through selection to increase the total number of eggs laid per hen. This number is the cumulative result of the oviposition, a cyclic and repeated process which leads to a pattern over time (the egg laying curve) which can be modelled and described individually. Unlike the total egg number which compounds all variations, the shape of the curve gives information on the different phases of egg laying, and its genetic analysis using molecular markers might contribute to understand better the underlying mechanisms. The purpose of this study was to perform the first QTL search for traits involved in shaping the egg laying curve, in an F₂ experiment with 359 female Japanese quail.

Results

Eight QTL were found on five autosomes, and six of them could be directly associated with egg production traits, although none was significant at the genome-wide level. One of them (on CJA13) had an effect on the first part of the laying curve, before the production peak. Another one (on CJA06) was related to the central part of the curve when laying is maintained at a high level, and the four others (on CJA05, CJA10 and CJA14) acted on the last part of the curve where persistency is determinant. The QTL for the central part of the curve was mapped at the same position on CJA06 than a genome-wide significant QTL for total egg number detected previously in the same F₂.

Conclusion

Despite its limited scope (number of microsatellites, size of the phenotypic data set), this work has shown that it was possible to use the individual egg laying data collected daily to find new QTL which affect the shape of the egg laying curve. Beyond the present results, this new approach could also be applied to longitudinal traits in other species, like growth and lactation in ruminants, for which good marker coverage of the genome and theoretical models with a biological significance are available.

Integrated maps in quail (Coturnix japonica) confirm the high degree of synteny conservation with chicken (Gallus gallus) despite 35 million years of divergence

BACKGROUND

By comparing the quail genome with that of chicken, chromosome rearrangements that have occurred in these two galliform species over 35 million years of evolution can be detected. From a more practical point of view, the definition of conserved syntenies helps to predict the position of genes in quail, based on information taken from the chicken sequence, thus enhancing the utility of this species in biological studies through a better knowledge of its genome structure. A microsatellite and an Amplified Fragment Length Polymorphism (AFLP) genetic map were previously published for quail, as well as comparative cytogenetic data with chicken for macrochromosomes. Quail genomics will benefit from the extension and the integration of these maps.

RESULTS

The integrated linkage map presented here is based on segregation analysis of both anonymous markers and functional gene loci in 1,050 quail from three independent F2 populations. Ninety-two loci are resolved into 14 autosomal linkage groups and a Z chromosome-specific linkage group, aligned with the quail AFLP map. The size of linkage groups ranges from 7.8 cM to 274.8 cM. The total map distance covers 904.3 cM with an average spacing of 9.7 cM between loci. The coverage is not complete, as macrochromosome CJA08, the gonosome CJAW and 23 microchromosomes have no marker assigned yet. Significant sequence identities of quail markers with chicken enabled the alignment of the quail linkage groups on the chicken genome sequence assembly. This, together with interspecific Fluorescence In Situ Hybridization (FISH), revealed very high similarities in marker order between the two species for the eight macrochromosomes and the 14 microchromosomes studied.

CONCLUSION

Integrating the two microsatellite and the AFLP quail genetic maps greatly enhances the quality of the resulting information and will thus facilitate the identification of Quantitative Trait Loci (QTL). The alignment with the chicken chromosomes confirms the high conservation of gene order that was expected between the two species for macrochromosomes. By extending the comparative study to the microchromosomes, we suggest that a wealth of information can be mined in chicken, to be used for genome analyses in quail.

Fatness QTL on chicken chromosome 5 and interaction with sex

Quantitative trait loci (QTL) affecting fatness in male chickens were previously identified on chromosome 5 (GGA5) in a three-generation design derived from two experimental chicken lines divergently selected for abdominal fat weight. A new design, established from the same pure lines, produced 407 F2 progenies (males and females) from 4 F1-sire families. Body weight and abdominal fat were measured on the F2 at 9 wk of age. In each sire family, selective genotyping was carried out for 48 extreme individuals for abdominal fat using seven microsatellite markers from GGA5. QTL analyses confirmed the presence of QTL for fatness on GGA5 and identified a QTL by sex interaction. By crossing one F1 sire heterozygous at the QTL with lean line dams, three recombinant backcross 1 (BC1) males were produced and their QTL genotypes were assessed in backcross 2 (BC2) progenies. These results confirmed the QTL by sex interaction identified in the F2 generation and they allow mapping of the female QTL to less than 8 Mb at the distal part of the GGA5. They also indicate that fat QTL alleles were segregating in both fat and lean lines.

Mapping quantitative trait loci affecting fatness and breast muscle weight in meat-type chicken lines divergently selected on abdominal fatness

Quantitative trait loci (QTL) for abdominal fatness and breast muscle weight were investigated in a three-generation design performed by inter-crossing two experimental meat-type chicken lines that were divergently selected on abdominal fatness. A total of 585 F2 male offspring from 5 F1 sires and 38 F1 dams were recorded at 8 weeks of age for live body, abdominal fat and breast muscle weights. One hundred-twenty nine microsatellite markers, evenly located throughout the genome and heterozygous for most of the F1 sires, were used for genotyping the F2 birds. In each sire family, those offspring exhibiting the most extreme values for each trait were genotyped. Multipoint QTL analyses using maximum likelihood methods were performed for abdominal fat and breast muscle weights, which were corrected for the effects of 8-week body weight, dam and hatching group. Isolated markers were assessed by analyses of variance. Two significant QTL were identified on chromosomes 1 and 5 with effects of about one within-family residual standard deviation. One breast muscle QTL was identified on GGA1 with an effect of 2.0 within-family residual standard deviation.

A genome scan with AFLP markers to detect fearfulness-related QTLs in Japanese quail

A quantitative trait loci (QTL) study was undertaken to identify genome regions involved in the control of fearfulness in Japanese quail (Coturnix japonica). An F2 cross was made between two quail lines divergently selected over 29 generations on duration of tonic immobility (DTI), a catatonic-like state of reduced responsiveness to a stressful stimulation. A total of 1065 animals were measured for the logarithm of DTI (LOGTI), the number of inductions (NI) necessary to induce the immobility reaction, open-field behaviour including locomotor activity (MOVE), latency before first movement (LAT), number of jumps (JUMP), dejections (DEJ) and shouts (SHOUT), corticosterone level after a contention stress (LOGCORT) and body weight at 2 weeks of age (BW2). A total of 310 animals were included in a genome scan using selective genotyping with 248 AFLP markers. A total of 21 suggestive or genome-wide significant QTL were observed. Two highly significant QTL were identified on linkage group 1 (GL1), one for LOGTI and one for NI. In the vicinity of the QTL for LOGTI, a nearly significant QTL for SHOUT and a suggestive QTL for LAT were also identified. On GL3, genome-wide significant QTL were observed for JUMP and DEJ as well as suggestive QTL for LOGTI, MOVE, SHOUT and LAT. A significant QTL for BW2 was observed on GL2 and a nearly significant one on GL1. These results may be useful in the understanding of fearfulness in quail and related species provided that fearfulness has the same genetic basis.

A gene-based radiation hybrid map of chicken microchromosome 14: comparison to human and alignment to the assembled chicken sequence

We present a gene-based RH map of the chicken microchromosome GGA14, known to have synteny conservations with human chromosomal regions HSA16p13.3 and HSA17p11.2. Microsatellite markers from the genetic map were used to check the validity of the RH map and additional markers were developed from chicken EST data to yield comparative mapping data. A high rate of intra-chromosomal rearrangements was detected by comparison to the assembled human sequence. Finally, the alignment of the RH map to the assembled chicken sequence showed a small number of discordances, most of which involved the same region of the chromosome spanning between 40.5 and 75.9 cR(6000) on the RH map.