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

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

Reintroducing genetic diversity in populations from cryopreserved material: the case of Abondance, a French local dairy cattle breed.

BACKGROUND: Genetic diversity is a necessary condition for populations to evolve under natural adaptation, artificial selection, or both. However, genetic diversity is often threatened, in particular in domestic animal populations where artificial selection, genetic drift and inbreeding are strong. In this context, cryopreserved genetic resources are a promising option to reintroduce lost variants and to limit inbreeding. However, while the use of ancient genetic resources is more common in plant breeding, it is less documented in animals due to a longer generation interval, making it difficult to fill the gap in performance due to continuous selection. This study investigates one of the only concrete cases available in animals, for which cryopreserved semen from a bull born in 1977 in a lost lineage was introduced into the breeding scheme of a French local dairy cattle breed, the Abondance breed, more than 20 years later. RESULTS: We found that this re-introduced bull was genetically distinct with respect to the current population and thus allowed part of the genetic diversity lost over time to be restored. The expected negative gap in milk production due to continuous selection was absorbed in a few years by preferential mating with elite cows. Moreover, the re-use of this bull more than two decades later did not increase the level of inbreeding, and even tended to reduce it by avoiding mating with relatives. Finally, the reintroduction of a bull from a lost lineage in the breeding scheme allowed for improved performance for reproductive abilities, a trait that was less subject to selection in the past. CONCLUSIONS: The use of cryopreserved material is an efficient way to manage the genetic diversity of an animal population, by mitigating the effects of both inbreeding and strong selection. However, attention should be paid to mating of animals to limit the disadvantages associated with incorporating original genetic material, notably a discrepancy in the breeding values for selected traits or an increase in inbreeding. Therefore, careful characterization of the genetic resources available in cryobanks could help to ensure the sustainable management of populations, in particular local or small populations. These results could also be transferred to the conservation of wild threatened populations.

Managing genetic diversity in breeding programs of small populations: the case of French local chicken breeds.

BACKGROUND: On-going climate change will drastically modify agriculture in the future, with a need for more sustainable systems, in particular regarding animal production. In this context, genetic diversity is a key factor for adaptation to new conditions: local breeds likely harbor unique adaptive features and represent a key component of diversity to reach resilience. However, local breeds often suffer from small population sizes, which puts these valuable resources at risk of extinction. In chickens, population management programs were initiated a few decades ago in France, relying on a particular niche market that aims at promoting and protecting local breeds. We conducted a unique comprehensive study of 22 French local breeds, along with four commercial lines, to evaluate their genetic conservation status and the efficiency of the population management programs. RESULTS: Using a 57K single nucleotide polymorphism (SNP) chip, we demonstrated that both the between- and within-breed genetic diversity levels are high in the French local chicken populations. Diversity is mainly structured according to the breeds' selection and history. Nevertheless, we observed a prominent sub-structuring of breeds according to farmers' practices in terms of exchange, leading to more or less isolated flocks. By analysing demographic parameters and molecular information, we showed that consistent management programs are efficient in conserving genetic diversity, since breeds that integrated such programs earlier had older inbreeding. CONCLUSIONS: Management programs of French local chicken breeds have maintained their genetic diversity at a good level. We recommend that future programs sample as many individuals as possible, with emphasis on both males and females from the start, and focus on a quick and strong increase of population size while conserving as many families as possible. We also stress the usefulness of molecular tools to monitor small populations for which pedigrees are not always available. Finally, the breed appears to be an appropriate operational unit for the conservation of genetic diversity, even for local breeds, for which varieties, if present, could also be taken into account.

A research agenda for scaling up agroecology in European countries.

A profound transformation of agricultural production methods has become unavoidable due to the increase in the world's population, and environmental and climatic challenges. Agroecology is now recognized as a challenging model for agricultural systems, promoting their diversification and adaptation to environmental and socio-economic contexts, with consequences for the entire agri-food system and the development of rural and urban areas. Through a prospective exercise performed at a large interdisciplinary institute, INRAE, a research agenda for agroecology was built that filled a gap through its ambition and interdisciplinarity. It concerned six topics. For genetics, there is a need to study genetic aspects of complex systems (e.g., mixtures of genotypes) and to develop breeding methods for them. For landscapes, challenges lie in effects of heterogeneity at multiple scales, in multifunctionality and in the design of agroecological landscapes. Agricultural equipment and digital technologies show high potential for monitoring dynamics of agroecosystems. For modeling, challenges include approaches to complexity, consideration of spatial and temporal dimensions and representation of the cascade from cropping practices to ecosystem services. The agroecological transition of farms calls for modeling and observational approaches as well as for creating new design methods. Integration of agroecology into food systems raises the issues of product specificity, consumer behavior and organization of markets, standards and public policies. In addition, transversal priorities were identified: (i) generating sets of biological data, through research and participatory mechanisms, that are appropriate for designing agroecological systems and (ii) collecting and using coherent sets of data to enable assessment of vulnerability, resilience and risk in order to evaluate the performance of agroecological systems and to contribute to scaling up. The main lessons learned from this collective exercise can be useful for the entire scientific community engaged in research into agroecology.

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.

The evolution of Sex-linked barring alleles in chickens involves both regulatory and coding changes in CDKN2A.

Sex-linked barring is a fascinating plumage pattern in chickens recently shown to be associated with two non-coding and two missense mutations affecting the ARF transcript at the CDKN2A tumor suppressor locus. It however remained a mystery whether all four mutations are indeed causative and how they contribute to the barring phenotype. Here, we show that Sex-linked barring is genetically heterogeneous, and that the mutations form three functionally different variant alleles. The B0 allele carries only the two non-coding changes and is associated with the most dilute barring pattern, whereas the B1 and B2 alleles carry both the two non-coding changes and one each of the two missense mutations causing the Sex-linked barring and Sex-linked dilution phenotypes, respectively. The data are consistent with evolution of alleles where the non-coding changes occurred first followed by the two missense mutations that resulted in a phenotype more appealing to humans. We show that one or both of the non-coding changes are cis-regulatory mutations causing a higher CDKN2A expression, whereas the missense mutations reduce the ability of ARF to interact with MDM2. Caspase assays for all genotypes revealed no apoptotic events and our results are consistent with a recent study indicating that the loss of melanocyte progenitors in Sex-linked barring in chicken is caused by premature differentiation and not apoptosis. Our results show that CDKN2A is a major locus driving the differentiation of avian melanocytes in a temporal and spatial manner.

Multi-species annotation of transcriptome and chromatin structure in domesticated animals.

BACKGROUND: Comparative genomics studies are central in identifying the coding and non-coding elements associated with complex traits, and the functional annotation of genomes is a critical step to decipher the genotype-to-phenotype relationships in livestock animals. As part of the Functional Annotation of Animal Genomes (FAANG) action, the FR-AgENCODE project aimed to create reference functional maps of domesticated animals by profiling the landscape of transcription (RNA-seq), chromatin accessibility (ATAC-seq) and conformation (Hi-C) in species representing ruminants (cattle, goat), monogastrics (pig) and birds (chicken), using three target samples related to metabolism (liver) and immunity (CD4+ and CD8+ T cells). RESULTS: RNA-seq assays considerably extended the available catalog of annotated transcripts and identified differentially expressed genes with unknown function, including new syntenic lncRNAs. ATAC-seq highlighted an enrichment for transcription factor binding sites in differentially accessible regions of the chromatin. Comparative analyses revealed a core set of conserved regulatory regions across species. Topologically associating domains (TADs) and epigenetic A/B compartments annotated from Hi-C data were consistent with RNA-seq and ATAC-seq data. Multi-species comparisons showed that conserved TAD boundaries had stronger insulation properties than species-specific ones and that the genomic distribution of orthologous genes in A/B compartments was significantly conserved across species. CONCLUSIONS: We report the first multi-species and multi-assay genome annotation results obtained by a FAANG project. Beyond the generation of reference annotations and the confirmation of previous findings on model animals, the integrative analysis of data from multiple assays and species sheds a new light on the multi-scale selective pressure shaping genome organization from birds to mammals. Overall, these results emphasize the value of FAANG for research on domesticated animals and reinforces the importance of future meta-analyses of the reference datasets being generated by this community on different species.

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

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

Detection of QTL for traits related to adaptation to sub-optimal climatic conditions in chickens.

BACKGROUND: Growth traits can be used as indicators of adaptation to sub-optimal conditions. The current study aimed at identifying quantitative trait loci (QTL) that control performance under variable temperature conditions in chickens. METHODS: An F2 population was produced by crossing the Taiwan Country chicken L2 line (selected for body weight, comb area, and egg production) with an experimental line of Rhode Island Red layer R- (selected for low residual feed consumption). A total of 844 animals were genotyped with the 60 K Illumina single nucleotide polymorphism (SNP) chip. Whole-genome interval linkage mapping and a genome-wide association study (GWAS) were performed for body weight at 0, 4, 8, 12, and 16 weeks of age, shank length at 8 weeks of age, size of comb area at 16 weeks of age, and antibody response to sheep red blood cells at 11 weeks of age (7 and 14 days after primary immunization). Relevant genes were identified based on functional annotation of candidate genes and potentially relevant SNPs were detected by comparing whole-genome sequences of several birds between the parental lines. RESULTS: Whole-genome QTL analysis revealed 47 QTL and 714 effects associated with 178 SNPs were identified by GWAS with 5% Bonferroni genome-wide significance. Little overlap was observed between the QTL and GWAS results, with only two chromosomal regions detected by both approaches, i.e. one on GGA24 (GGA for Gallus gallus chromosome) for BW04 and one on GGAZ for six growth-related traits. Based on whole-genome sequence, differences between the parental lines based on several birds were screened in the genome-wide QTL regions and in a region detected by both methods, resulting in the identification of 106 putative candidate genes with a total of 15,443 SNPs, of which 41 were missense and 1698 were not described in the dbSNP archive. CONCLUSIONS: The QTL detected in this study for growth and morphological traits likely influence adaptation of chickens to sub-tropical climate. Using whole-genome sequence data, we identified candidate SNPs for further confirmation of QTL in the F2 design. A strong QTL effect found on GGAZ underlines the importance of sex-linked inheritance for growth traits in chickens.

Whole-genome resequencing reveals loci under selection during chicken domestication.

Domestic animals are excellent models for genetic studies of phenotypic evolution. They have evolved genetic adaptations to a new environment, the farm, and have been subjected to strong human-driven selection leading to remarkable phenotypic changes in morphology, physiology and behaviour. Identifying the genetic changes underlying these developments provides new insight into general mechanisms by which genetic variation shapes phenotypic diversity. Here we describe the use of massively parallel sequencing to identify selective sweeps of favourable alleles and candidate mutations that have had a prominent role in the domestication of chickens (Gallus gallus domesticus) and their subsequent specialization into broiler (meat-producing) and layer (egg-producing) chickens. We have generated 44.5-fold coverage of the chicken genome using pools of genomic DNA representing eight different populations of domestic chickens as well as red jungle fowl (Gallus gallus), the major wild ancestor. We report more than 7,000,000 single nucleotide polymorphisms, almost 1,300 deletions and a number of putative selective sweeps. One of the most striking selective sweeps found in all domestic chickens occurred at the locus for thyroid stimulating hormone receptor (TSHR), which has a pivotal role in metabolic regulation and photoperiod control of reproduction in vertebrates. Several of the selective sweeps detected in broilers overlapped genes associated with growth, appetite and metabolic regulation. We found little evidence that selection for loss-of-function mutations had a prominent role in chicken domestication, but we detected two deletions in coding sequences that we suggest are functionally important. This study has direct application to animal breeding and enhances the importance of the domestic chicken as a model organism for biomedical research.

The Rose-comb mutation in chickens constitutes a structural rearrangement causing both altered comb morphology and defective sperm motility.

Rose-comb, a classical monogenic trait of chickens, is characterized by a drastically altered comb morphology compared to the single-combed wild-type. Here we show that Rose-comb is caused by a 7.4 Mb inversion on chromosome 7 and that a second Rose-comb allele arose by unequal crossing over between a Rose-comb and wild-type chromosome. The comb phenotype is caused by the relocalization of the MNR2 homeodomain protein gene leading to transient ectopic expression of MNR2 during comb development. We also provide a molecular explanation for the first example of epistatic interaction reported by Bateson and Punnett 104 years ago, namely that walnut-comb is caused by the combined effects of the Rose-comb and Pea-comb alleles. Transient ectopic expression of MNR2 and SOX5 (causing the Pea-comb phenotype) occurs in the same population of mesenchymal cells and with at least partially overlapping expression in individual cells in the comb primordium. Rose-comb has pleiotropic effects, as homozygosity in males has been associated with poor sperm motility. We postulate that this is caused by the disruption of the CCDC108 gene located at one of the inversion breakpoints. CCDC108 is a poorly characterized protein, but it contains a MSP (major sperm protein) domain and is expressed in testis. The study illustrates several characteristic features of the genetic diversity present in domestic animals, including the evolution of alleles by two or more consecutive mutations and the fact that structural changes have contributed to fast phenotypic evolution.

Cryptic patterning of avian skin confers a developmental facility for loss of neck feathering.

Vertebrate skin is characterized by its patterned array of appendages, whether feathers, hairs, or scales. In avian skin the distribution of feathers occurs on two distinct spatial levels. Grouping of feathers within discrete tracts, with bare skin lying between the tracts, is termed the macropattern, while the smaller scale periodic spacing between individual feathers is referred to as the micropattern. The degree of integration between the patterning mechanisms that operate on these two scales during development and the mechanisms underlying the remarkable evolvability of skin macropatterns are unknown. A striking example of macropattern variation is the convergent loss of neck feathering in multiple species, a trait associated with heat tolerance in both wild and domestic birds. In chicken, a mutation called Naked neck is characterized by a reduction of body feathering and completely bare neck. Here we perform genetic fine mapping of the causative region and identify a large insertion associated with the Naked neck trait. A strong candidate gene in the critical interval, BMP12/GDF7, displays markedly elevated expression in Naked neck embryonic skin due to a cis-regulatory effect of the causative mutation. BMP family members inhibit embryonic feather formation by acting in a reaction-diffusion mechanism, and we find that selective production of retinoic acid by neck skin potentiates BMP signaling, making neck skin more sensitive than body skin to suppression of feather development. This selective production of retinoic acid by neck skin constitutes a cryptic pattern as its effects on feathering are not revealed until gross BMP levels are altered. This developmental modularity of neck and body skin allows simple quantitative changes in BMP levels to produce a sparsely feathered or bare neck while maintaining robust feather patterning on the body.

Diversity and evolution of the highly polymorphic tandem repeat LEI0258 in the chicken MHC-B region.

The chicken major histocompatibility complex (MHC) is located on the microchromosome 16 and is described as the most variable region in the genome. The genes of the MHC play a central role in the immune system. Particularly, genes encoding proteins involved in the antigen presentation to T cells. Therefore, describing the genetic polymorphism of this region is crucial in understanding host-pathogen interactions. The tandem repeat LEI0258 is located within the core area of the B region of the chicken MHC (MHC-B region) and its genotypes correlate with serology. This marker was used to provide a picture of the worldwide diversity of the chicken MHC-B region and to categorize chicken MHC haplotypes. More than 1,600 animals from 80 different populations or lines of chickens from Africa, Asia, and Europe, including wild fowl species, were genotyped at the LEI0258 locus. Fifty novel alleles were described after sequencing. The resulting 79 alleles were classified into 12 clusters, based on the SNPs and indels found within the sequences flanking the repeats. Furthermore, hypotheses were formulated on the evolutionary dynamics of the region. This study constitutes the largest variability report for the chicken MHC and establishes a framework for future diversity or association studies.

Gene diversity, agroecological structure and introgression patterns among village chicken populations across North, West and Central Africa.

BACKGROUND: Chickens represent an important animal genetic resource for improving farmers' income in Africa. The present study provides a comparative analysis of the genetic diversity of village chickens across a subset of African countries. Four hundred seventy-two chickens were sampled in 23 administrative provinces across Cameroon, Benin, Ghana, Côte d'Ivoire, and Morocco. Geographical coordinates were recorded to analyze the relationships between geographic distribution and genetic diversity. Molecular characterization was performed with a set of 22 microsatellite markers. Five commercial lines, broilers and layers, were also genotyped to investigate potential gene flow. A genetic diversity analysis was conducted both within and between populations. RESULTS: High heterozygosity levels, ranging from 0.51 to 0.67, were reported for all local populations, corresponding to the values usually found in scavenging populations worldwide. Allelic richness varied from 2.04 for a commercial line to 4.84 for one population from Côte d'Ivoire. Evidence of gene flow between commercial and local populations was observed in Morocco and in Cameroon, which could be related to long-term improvement programs with the distribution of crossbred chicks. The impact of such introgressions seemed rather limited, probably because of poor adaptation of exotic birds to village conditions, and because of the consumers' preference for local chickens. No such gene flow was observed in Benin, Ghana, and Côte d'Ivoire, where improvement programs are also less developed. The clustering approach revealed an interesting similarity between local populations found in regions sharing high levels of precipitation, from Cameroon to Côte d'Ivoire. Restricting the study to Benin, Ghana, and Côte d'Ivoire, did not result in a typical breed structure but a south-west to north-east gradient was observed. Three genetically differentiated areas (P<0.01) were identified, matching with Major Farming Systems (namely Tree Crop, Cereal-Root Crop, and Root Crop) described by the FAO. CONCLUSIONS: Local chickens form a highly variable gene pool constituting a valuable resource for human populations. Climatic conditions, farming systems, and cultural practices may influence the genetic diversity of village chickens in Africa. A higher density of markers would be needed to identify more precisely the relative importance of these factors.

Copy number variation in intron 1 of SOX5 causes the Pea-comb phenotype in chickens.

Pea-comb is a dominant mutation in chickens that drastically reduces the size of the comb and wattles. It is an adaptive trait in cold climates as it reduces heat loss and makes the chicken less susceptible to frost lesions. Here we report that Pea-comb is caused by a massive amplification of a duplicated sequence located near evolutionary conserved non-coding sequences in intron 1 of the gene encoding the SOX5 transcription factor. This must be the causative mutation since all other polymorphisms associated with the Pea-comb allele were excluded by genetic analysis. SOX5 controls cell fate and differentiation and is essential for skeletal development, chondrocyte differentiation, and extracellular matrix production. Immunostaining in early embryos demonstrated that Pea-comb is associated with ectopic expression of SOX5 in mesenchymal cells located just beneath the surface ectoderm where the comb and wattles will subsequently develop. The results imply that the duplication expansion interferes with the regulation of SOX5 expression during the differentiation of cells crucial for the development of comb and wattles. The study provides novel insight into the nature of mutations that contribute to phenotypic evolution and is the first description of a spontaneous and fully viable mutation in this developmentally important gene.

Identification of the yellow skin gene reveals a hybrid origin of the domestic chicken.

Yellow skin is an abundant phenotype among domestic chickens and is caused by a recessive allele (W*Y) that allows deposition of yellow carotenoids in the skin. Here we show that yellow skin is caused by one or more cis-acting and tissue-specific regulatory mutation(s) that inhibit expression of BCDO2 (beta-carotene dioxygenase 2) in skin. Our data imply that carotenoids are taken up from the circulation in both genotypes but are degraded by BCDO2 in skin from animals carrying the white skin allele (W*W). Surprisingly, our results demonstrate that yellow skin does not originate from the red junglefowl (Gallus gallus), the presumed sole wild ancestor of the domestic chicken, but most likely from the closely related grey junglefowl (Gallus sonneratii). This is the first conclusive evidence for a hybrid origin of the domestic chicken, and it has important implications for our views of the domestication process.

Ex situ and in situ data for endangered livestock breeds in Spain.

Improvements in ex situ storage of genetic and reproductive materials offer an alternative for endangered livestock breed conservation. This paper presents a dataset for current ex situ collections and in situ population for 179 Spanish livestock breeds of seven species, cattle, sheep, pig, chicken, goat, horse and donkey. Ex situ data was obtained via survey administered to 18 functioning gene banks in Spain and relates to the reproductive genetic materials (semen doses) of 210 livestock breeds distributed across the gene banks. In situ data combines CENSUS information with linear regression techniques and relates to the geographic distribution of 179 Spanish autochthonous livestock breeds (2009-2018), and in situ population projections and extinction probabilities (2019-2060). We use a decision variable defining an "acceptable level of risk" that allows decision makers to specify tolerable levels of in situ breed endangerment when taking ex situ collection and storage decisions.

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

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

Tissue Resources for the Functional Annotation of Animal Genomes.

In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

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

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