Search

Abstract Background The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses. Results Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants. Conclusions Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

Abstract 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.

Abstract 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.

Abstract 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.

Abstract Background The aim of this study was to analyse the mechanisms that underlie phenotypic quantitative trait loci (QTL) in overfed mule ducks by identifying co-localized proteomic QTL (pQTL). The QTL design consisted of three families of common ducks that were progeny-tested by using 294 male mule ducks. This population of common ducks was genotyped using a genetic map that included 334 genetic markers located across 28 APL chromosomes (APL for Anas platyrhynchos). Mule ducks were phenotyped for 49 traits related to growth, metabolism, overfeeding ability and meat and fatty liver quality, and 326 soluble fatty liver proteins were quantified. Results One hundred and seventy-six pQTL and 80 phenotypic QTL were detected at the 5% chromosome-wide significance threshold. The great majority of the identified pQTL were trans-acting and localized on a chromosome other than that carrying the coding gene. The most significant pQTL (1% genome-wide significance) were found for alpha-enolase on APL18 and fatty acid synthase on APL24. Some proteins were associated with numerous pQTL (for example, 17 and 14 pQTL were detected for alpha-enolase and apolipoprotein A1, respectively) and pQTL hotspots were observed on some chromosomes (APL18, 24, 25 and 29). We detected 66 co-localized phenotypic QTL and pQTL for which the significance of the two-trait QTL (2t-QTL) analysis was higher than that of the strongest QTL using a single-trait approach. Among these, 16 2t-QTL were pleiotropic. For example, on APL15, melting rate and abundance of two alpha-enolase spots appeared to be impacted by a single locus that is involved in the glycolytic process. On APLZ, we identified a pleiotropic QTL that modified both the blood level of glucose at the beginning of the force-feeding period and the concentration of glutamate dehydrogenase, which, in humans, is involved in increased glucose absorption by the liver when the glutamate dehydrogenase 1 gene is mutated. Conclusions We identified pleiotropic loci that affect metabolic pathways linked to glycolysis or lipogenesis, and in the end to fatty liver quality. Further investigation, via transcriptomics and metabolomics approaches, is required to confirm the biomarkers that were found to impact the genetic variability of these phenotypic traits.

Abstract 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.

Abstract 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.

Abstract Selection for increased resistance to Salmonella colonisation and excretion could reduce the risk of foodborne Salmonella infection. In order to identify potential loci affecting resistance, differences in resistance were identified between the N and 61 inbred lines and two QTL research performed. In an F2 cross, the animals were inoculated at one week of age with Salmonella enteritidis and cloacal swabs were carried out 4 and 5 wk post inoculation (thereafter called CSW4F2 and CSW4F2) and caecal contamination (CAECF2) was assessed 1 week later. The animals from the (N × 61) × N backcross were inoculated at six weeks of age with Salmonella typhimurium and cloacal swabs were studied from wk 1 to 4 (thereafter called CSW1BC to CSW4BC). A total of 33 F2 and 46 backcross progeny were selectively genotyped for 103 and 135 microsatellite markers respectively. The analysis used least-squares-based and non-parametric interval mapping. Two genome-wise significant QTL were observed on Chromosome 1 for CSW2BC and on Chromosome 2 for CSW4F2, and four suggestive QTL for CSW5F2 on Chromosome 2, for CSW5F2 and CSW2BC on chromosome 5 and for CAECF2 on chromosome 16. These results suggest new regions of interest and the putative role of SAL1.

Abstract 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 cR6000 on the RH map.

Abstract An attempt was made to assign five genes, CPS1, OTC, ASS, CRYD2, and ARG2, to chicken chromosomes (GGA) by radiation-hybrid mapping. OTC was assigned to GGA1; ARG2 to GGA5; CPS1 to GGA7; and CRYD2 to GGA19. ASS was not, however, assigned to a specific chromosomal position.

Abstract The quail is a valuable farm and laboratory animal. Yet molecular information about this species remains scarce. We present here the first genetic linkage map of the Japanese quail. This comprehensive map is based solely on amplified fragment length polymorphism (AFLP) markers. These markers were developed and genotyped in an F2 progeny from a cross between two lines of quail differing in stress reactivity. A total of 432 polymorphic AFLP markers were detected with 24 TaqI/EcoRI primer combinations. On average, 18 markers were produced per primer combination. Two hundred and fifty eight of the polymorphic markers were assigned to 39 autosomal linkage groups plus the ZW sex chromosome linkage groups. The linkage groups range from 2 to 28 markers and from 0.0 to 195.5 cM. The AFLP map covers a total length of 1516 cM, with an average genetic distance between two consecutive markers of 7.6 cM. This AFLP map can be enriched with other marker types, especially mapped chicken genes that will enable to link the maps of both species and make use of the powerful comparative mapping approach. This AFLP map of the Japanese quail already provides an efficient tool for quantitative trait loci (QTL) mapping.

Abstract In a project on the biodiversity of chickens funded by the European Commission (EC), eight laboratories collaborated to assess the genetic variation within and between 52 populations from a wide range of chicken types. Twenty-two di-nucleotide microsatellite markers were used to genotype DNA pools of 50 birds from each population. The polymorphism measures for the average, the least polymorphic population (inbred C line) and the most polymorphic population (Gallus gallus spadiceus) were, respectively, as follows: number of alleles per locus, per population: 3.5, 1.3 and 5.2; average gene diversity across markers: 0.47, 0.05 and 0.64; and proportion of polymorphic markers: 0.91, 0.25 and 1.0. These were in good agreement with the breeding history of the populations. For instance, unselected populations were found to be more polymorphic than selected breeds such as layers. Thus DNA pools are effective in the preliminary assessment of genetic variation of populations and markers. Mean genetic distance indicates the extent to which a given population shares its genetic diversity with that of the whole tested gene pool and is a useful criterion for conservation of diversity. The distribution of population-specific (private) alleles and the amount of genetic variation shared among populations supports the hypothesis that the red jungle fowl is the main progenitor of the domesticated chicken.

Abstract As a first step towards the development of radiation hybrid maps, we have produced a radiation hybrid panel in the chicken by fusing female embryonic diploid fibroblasts irradiated at 6 000 rads with HPRT-deficient hamster Wg3hCl2 cells. Due to the low retention frequency of the chicken fragments, a high number of clones was produced from which the best ones were selected. Thus, 452 fusion clones were tested for retention frequencies with a panel of 46 markers. Based on these results, 103 clones with a mean marker retention of 23.8% were selected for large scale culture to produce DNA in sufficient quantities for the genotyping of numerous markers. Retention frequency was tested again with the same 46 markers and the 90 best clones, with a final mean retention frequency of 21.9%, were selected for the final panel. This panel will be a valuable resource for fine mapping of markers and genes in the chicken, and will also help in building BAC contigs.

Abstract During the last ten years, the use of molecular markers, revealing polymorphism at the DNA level, has been playing an increasing part in animal genetics studies. Amongst others, the microsatellite DNA marker has been the most widely used, due to its easy use by simple PCR, followed by a denaturing gel electrophoresis for allele size determination, and to the high degree of information provided by its large number of alleles per locus. Despite this, a new marker type, named SNP, for Single Nucleotide Polymorphism, is now on the scene and has gained high popularity, even though it is only a bi-allelic type of marker. In this review, we will discuss the reasons for this apparent step backwards, and the pertinence of the use of SNPs in animal genetics, in comparison with other marker types.

Abstract The bulked segregant analysis methodology has been used to map, with microsatellite markers, two morphological mutations in the chicken: polydactyly (PO) and naked neck (NA). These autosomal mutations show partial dominance for NA, and dominance with incomplete penetrance for PO. They were mapped previously to different linkage groups of the classical map, PO to the linkage group IV and NA being linked to the erythrocyte antigen CPPP. An informative family of 70 offspring was produced by mating a sire, heterozygous for each of the mutations, to 7 dams homozygous recessive for each locus. Three DNA pools were prepared, pool PO included 20 chicks exhibiting at least one extra-toe, pool NA included 20 non-polydactyly chicks showing the typical phenotype associated with heterozygosity for the naked neck mutation, and pool NP included 20 chicks exhibiting neither of the mutant phenotypes. Typings were done on an ABI-373 automatic sequencer with 147 microsatellite markers covering most of the genome. An unbalanced distribution of sire marker alleles were detected between pool PO, and pools NA and NP, for two markers of chromosome 2p, MCW0082 and MCW0247. A linkage analysis taking into account the incomplete penetrance of polydactyly (80%) was performed with additional markers of this region and showed that the closest marker to the PO locus was MCW0071 (5 cM, lod score = 9). MCW0071 lies within the engrailed gene EN2 in the chicken. In the mouse, the homologous gene maps on chromosome 5, close to the hemimelic extra-toes mutation Hx. In the case of the NA locus, markers of chromosome 3 were selected because CPPP was mapped on this chromosome. Analysis of individual typings showed a linkage of 5.7 cM (lod score = 13) between the NA locus and ADL0237 in the distal region of chromosome 3q. These results contribute to connecting the former classical map to the molecular genetic map of the chicken, and open the way to the identification of the molecular nature of two developmental mutations of the chicken that are known to occur in many breeds of chickens.

The transformation of Egypt as a destination for mass tourism has profoundly reshaped the national territory. This article examines how the revolution of Egyptian tourism has led to the development of previously desert areas – the shores of the Red Sea and of the Mediterranean Sea – in an inversion of the previously Nile Valley-centred tourism maps. It shows that those new tourism territories are produced by a combination of processes acting on different geographical scales, from the local to the global. The article examines in particular the combined role of four players: the tourist, the state, the local entrepreneur, and the global tourism industry.

Abstract 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.