Evaluation of breeding strategies to reduce the inbreeding rate in the Friesian horse population: Looking back and moving forward.

In the past, small population sizes and unequal ancestor contributions have resulted in high inbreeding rates (ΔF) in the Friesian horse. Two decades ago, the studbook implemented a mating quota and started publishing individual kinships and reduced ΔF below 1% per generation. However, since then, the breeding population size has decreased and this raises the question whether current breeding strategies are sufficient to keep ΔF below desired rates. The aim of this study was to (1) reflect on past inbreeding trends and their main determinants, using pedigree analysis and (2) evaluate the effectiveness of the current and additional breeding strategies using stochastic simulations. We estimated the current ΔF (2013-2022) at 0.72% per generation. While the total contribution of the top 10 sires to the number of offspring per year has decreased from 75% in 1980 to 35% in 2022, this was mainly due to an increased number of approved studbook sires, and not due to more equalized contributions among sires. Of the simulated breeding strategies, selecting only breeding stallions with a below average mean kinship (i.e., "mean kinship selection") was most effective to decrease ΔF (from 0.66% to 0.33%). Increasing the number of breeding sires only had an effect when also a mating quota was applied. However, its effect remained limited. For example, a ~1.5 fold increase, combined with a mating quota of 80 offspring per sire per year, reduced ΔF from 0.55% to 0.51%. When increasing the number of breeding mares, a practically unfeasible large increase was needed for a meaningful reduction in ΔF (e.g. twice as many mares were needed to reduce ΔF from 0.66% to 0.56%). Stratified mating quotas, a novel approach in which we assigned each sire a mating quota (of 60, 80, 100 or 120 offspring per year) based on its mean kinship to recently born foals, resulted in a lower ΔF (0.43%) than a general mating quota of 90 offspring per sire per year (0.55%). Overall, while the current ΔF is below 1%, we recommend to implement additional strategies to further reduce ΔF below 0.5% in the Friesian horse population. For this breed and similar populations, we recommend to focus on breeding strategies based on kinship levels to effectively reduce ΔF.

Landscape genomics reveals regions associated with adaptive phenotypic and genetic variation in Ethiopian indigenous chickens.

Climate change is a threat to sustainable livestock production and livelihoods in the tropics. It has adverse impacts on feed and water availability, disease prevalence, production, environmental temperature, and biodiversity. Unravelling the drivers of local adaptation and understanding the underlying genetic variation in random mating indigenous livestock populations informs the design of genetic improvement programmes that aim to increase productivity and resilience. In the present study, we combined environmental, genomic, and phenotypic information of Ethiopian indigenous chickens to investigate their environmental adaptability. Through a hybrid sampling strategy, we captured wide biological and ecological variabilities across the country. Our environmental dataset comprised mean values of 34 climatic, vegetation and soil variables collected over a thirty-year period for 260 geolocations. Our biological dataset included whole genome sequences and quantitative measurements (on eight traits) from 513 individuals, representing 26 chicken populations spread along 4 elevational gradients (6-7 populations per gradient). We performed signatures of selection analyses ([Formula: see text] and XP-EHH) to detect footprints of natural selection, and redundancy analyses (RDA) to determine genotype-environment and genotype-phenotype-associations. RDA identified 1909 outlier SNPs linked with six environmental predictors, which have the highest contributions as ecological drivers of adaptive phenotypic variation. The same method detected 2430 outlier SNPs that are associated with five traits. A large overlap has been observed between signatures of selection identified by[Formula: see text]and XP-EHH showing that both methods target similar selective sweep regions. Average genetic differences measured by [Formula: see text] are low between gradients, but XP-EHH signals are the strongest between agroecologies. Genes in the calcium signalling pathway, those associated with the hypoxia-inducible factor (HIF) transcription factors, and sports performance (GALNTL6) are under selection in high-altitude populations. Our study underscores the relevance of landscape genomics as a powerful interdisciplinary approach to dissect adaptive phenotypic and genetic variation in random mating indigenous livestock populations.

Large scale sequence-based screen for recessive variants allows for identification and monitoring of rare deleterious variants in pigs.

Most deleterious variants are recessive and segregate at relatively low frequency. Therefore, high sample sizes are required to identify these variants. In this study we report a large-scale sequence based genome-wide association study (GWAS) in pigs, with a total of 120,000 Large White and 80,000 Synthetic breed animals imputed to sequence using a reference population of approximately 1,100 whole genome sequenced pigs. We imputed over 20 million variants with high accuracies (R2>0.9) even for low frequency variants (1-5% minor allele frequency). This sequence-based analysis revealed a total of 14 additive and 9 non-additive significant quantitative trait loci (QTLs) for growth rate and backfat thickness. With the non-additive (recessive) model, we identified a deleterious missense SNP in the CDHR2 gene reducing growth rate and backfat in homozygous Large White animals. For the Synthetic breed, we revealed a QTL on chromosome 15 with a frameshift variant in the OBSL1 gene. This QTL has a major impact on both growth rate and backfat, resembling human 3M-syndrome 2 which is related to the same gene. With the additive model, we confirmed known QTLs on chromosomes 1 and 5 for both breeds, including variants in the MC4R and CCND2 genes. On chromosome 1, we disentangled a complex QTL region with multiple variants affecting both traits, harboring 4 independent QTLs in the span of 5 Mb. Together we present a large scale sequence-based association study that provides a key resource to scan for novel variants at high resolution for breeding and to further reduce the frequency of deleterious alleles at an early stage in the breeding program.

Identification of the genetic basis of sow pelvic organ prolapse.

Introduction: Pelvic organ prolapse (POP) is one contributor to recent increases in sow mortality that have been observed in some populations and environments, leading to financial losses and welfare concerns. Methods: With inconsistent previous reports, the objective here was to investigate the role of genetics on susceptibility to POP, using data on 30,429 purebred sows, of which 14,186 were genotyped (25K), collected from 2012 to 2022 in two US multiplier farms with a high POP incidence of 7.1% among culled and dead sows and ranging from 2% to 4% of all sows present by parity. Given the low incidence of POP for parities 1 and >6, only data from parities 2 to 6 were retained for analyses. Genetic analyses were conducted both across parities, using cull data (culled for POP versus another reason), and by parity, using farrowing data. (culled for POP versus culled for another reason or not culled). Results and Discussion: Estimates of heritability from univariate logit models on the underlying scale were 0.35 ± 0.02 for the across-parity analysis and ranged from 0.41 ± 0.03 in parity 2 to 0.15 ± 0.07 in parity 6 for the by-parity analyses. Estimates of genetic correlations of POP between parities based on bivariate linear models indicated a similar genetic basis of POP across parities but less similar with increasing distance between parities. Genome wide association analyses revealed six 1 Mb windows that explained more than 1% of the genetic variance in the across-parity data. Most regions were confirmed in several by-parity analyses. Functional analyses of the identified genomic regions showed a potential role of several genes on chromosomes 1, 3, 7, 10, 12, and 14 in susceptibility to POP, including the Estrogen Receptor gene. Gene set enrichment analyses showed that genomic regions that explained more variation for POP were enriched for several terms from custom transcriptome and gene ontology libraries. Conclusion: The influence of genetics on susceptibility to POP in this population and environment was confirmed and several candidate genes and biological processes were identified that can be targeted to better understand and mitigate the incidence of POP.

Detailed molecular and epigenetic characterization of the pig IPEC-J2 and chicken SL-29 cell lines.

The pig IPEC-J2 and chicken SL-29 cell lines are of interest because of their untransformed nature and wide use in functional studies. Molecular characterization of these cell lines is important to gain insight into possible molecular aberrations. The aim of this paper is to provide a molecular and epigenetic characterization of the IPEC-J2 and SL-29 cell lines, a cell-line reference for the FAANG community, and future biomedical research. Whole genome sequencing, gene expression, DNA methylation, chromatin accessibility, and ChIP-seq of four histone marks (H3K4me1, H3K4me3, H3K27ac, H3K27me3) and an insulator (CTCF) are used to achieve these aims. Heteroploidy (aneuploidy) of various chromosomes was observed from whole genome sequencing analysis in both cell lines. Furthermore, higher gene expression for genes located on chromosomes with aneuploidy in comparison to diploid chromosomes was observed. Regulatory complexity of gene expression, DNA methylation, and chromatin accessibility was investigated through an integrative approach.

Preimplantation Genetic Testing for Aneuploidy (PGT-A) Reveals High Levels of Chromosomal Errors in In Vivo-Derived Pig Embryos, with an Increased Incidence When Produced In Vitro.

Preimplantation genetic testing for aneuploidy (PGT-A) is widespread, but controversial, in humans and improves pregnancy and live birth rates in cattle. In pigs, it presents a possible solution to improve in vitro embryo production (IVP), however, the incidence and origin of chromosomal errors remains under-explored. To address this, we used single nucleotide polymorphism (SNP)-based PGT-A algorithms in 101 in vivo-derived (IVD) and 64 IVP porcine embryos. More errors were observed in IVP vs. IVD blastocysts (79.7% vs. 13.6% p < 0.001). In IVD embryos, fewer errors were found at blastocyst stage compared to cleavage (4-cell) stage (13.6% vs. 40%, p = 0.056). One androgenetic and two parthenogenetic embryos were also identified. Triploidy was the most common error in IVD embryos (15.8%), but only observed at cleavage, not blastocyst stage, followed by whole chromosome aneuploidy (9.9%). In IVP blastocysts, 32.8% were parthenogenetic, 25.0% (hypo-)triploid, 12.5% aneuploid, and 9.4% haploid. Parthenogenetic blastocysts arose from just three out of ten sows, suggesting a possible donor effect. The high incidence of chromosomal abnormalities in general, but in IVP embryos in particular, suggests an explanation for the low success of porcine IVP. The approaches described provide a means of monitoring technical improvements and suggest future application of PGT-A might improve embryo transfer success.

Assessing the relationship between the in silico predicted consequences of 97 missense mutations mapping to 68 genes related to lipid metabolism and their association with porcine fatness traits.

In general, the relationship between the predicted functional consequences of missense mutations mapping to genes known to be involved in human diseases and the severity of disease manifestations is weak. In this study, we tested in pigs whether missense single nucleotide polymorphisms (SNPs), predicted to have consequences on the function of genes related to lipid metabolism are associated with lipid phenotypes. Association analysis demonstrated that nine out of 72 nominally associated SNPs were classified as "highly" or "very highly consistent" in silico-predicted functional mutations and did not show association with lipid traits expected to be affected by inactivation of the corresponding gene. Although the lack of endophenotypes and the limited sample size of certain genotypic classes might have limited to some extent the reach of the current study, our data indicate that present-day bioinformatic tools have a modest ability to predict the impact of missense mutations on complex phenotypes.

Corrigendum: A novel loss-of-function variant in transmembrane protein 263 (TMEM263) of autosomal dwarfism in chicken.

[This corrects the article DOI: 10.3389/fgene.2018.00193.].

Meta-Analysis of SNPs Determining Litter Traits in Pigs.

Nearly 2000 SNPs associated with pig litter size traits have been reported based on genome-wide association studies (GWASs). The aims of this study were to gather and integrate previously reported associations between SNPs and five litter traits: total number born (TNB), number born alive (NBA), number of stillborn (SB), litter birth weight (LWT), and corpus luteum number (CLN), in order to evaluate their common genetic background and to perform a meta-analysis (MA) of GWASs for total number born (TNB) recorded for animals from five pig populations. In this study, the genes with the largest number of associations with evaluated litter traits were GABRG3, RBP7, PRKD1, and STXBP6. Only 21 genes out of 233 associated with the evaluated litter traits were reported in more than one population or for more than one trait. Based on this evaluation, the most interesting candidate gene is PRKD1, which has an association with SB and TNB traits. Based on GO term analysis, PRKD1 was shown to be involved in angiogenesis as well. As a result of the MA, two new genomic regions, which have not been previously reported, were found to be associated with the TNB trait. One SNP was located on Sus scrofa chromosome (SSC) 14 in the intron of the FAM13C gene. The second SNP was located on SSC9 within the intron of the AGMO gene. Functional analysis revealed a strong candidate causal gene underlying the QTL on SSC9. The third best hit and the most promising candidate gene for litter size was found within the SOSTDC1 gene, associated with lower male fertility in rats. We showed that litter traits studied across pig populations have only a few genomic regions in common based on candidate gene comparison. PRKD1 could be an interesting candidate gene with a wider association with fertility. The MA identified new genomic regions on SSC9 and SSC14 associated with TNB. Further functional analysis indicated the most promising gene was SOSTDC1, which was confirmed to affect male fertility in other mammals. This is an important finding, as litter traits are by default linked with females rather than males.

Fine Mapping of a Major Backfat QTL Reveals a Causal Regulatory Variant Affecting the CCND2 Gene.

Backfat is an important trait in pork production, and it has been included in the breeding objectives of genetic companies for decades. Although adipose tissue is a good energy storage, excessive fat results in reduced efficiency and economical losses. A large QTL for backfat thickness on chromosome 5 is still segregating in different commercial pig breeds. We fine mapped this QTL region using a genome-wide association analysis (GWAS) with 133,358 genotyped animals from five commercial populations (Landrace, Pietrain, Large White, Synthetic, and Duroc) imputed to the porcine 660K SNP chip. The lead SNP was located at 5:66103958 (G/A) within the third intron of the CCND2 gene, with the G allele associated with more backfat, while the A allele is associated with less backfat. We further phased the QTL region to discover a core haplotype of five SNPs associated with low backfat across three breeds. Linkage disequilibrium analysis using whole-genome sequence data revealed three candidate causal variants within intronic regions and downstream of the CCND2 gene, including the lead SNP. We evaluated the association of the lead SNP with the expression of the genes in the QTL region (including CCND2) in a large cohort of 100 crossbred samples, sequenced in four different tissues (lung, spleen, liver, muscle). Results show that the A allele increases the expression of CCND2 in an additive way in three out of four tissues. Our findings indicate that the causal variant for this QTL region is a regulatory variant within the third intron of the CCND2 gene affecting the expression of CCND2.

A new haplotype-resolved turkey genome to enable turkey genetics and genomics research.

BACKGROUND: The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry. RESULTS: By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding. CONCLUSIONS: We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.

Review: Balancing Selection for Deleterious Alleles in Livestock.

Harmful alleles can be under balancing selection due to an interplay of artificial selection for the variant in heterozygotes and purifying selection against the variant in homozygotes. These pleiotropic variants can remain at moderate to high frequency expressing an advantage for favorable traits in heterozygotes, while harmful in homozygotes. The impact on the population and selection strength depends on the consequence of the variant both in heterozygotes and homozygotes. The deleterious phenotype expressed in homozygotes can range from early lethality to a slightly lower fitness in the population. In this review, we explore a range of causative variants under balancing selection including loss-of-function variation (i.e., frameshift, stop-gained variants) and regulatory variation (affecting gene expression). We report that harmful alleles often affect orthologous genes in different species, often influencing analogous traits. The recent discoveries are mainly driven by the increasing genomic and phenotypic resources in livestock populations. However, the low frequency and sometimes subtle effects in homozygotes prevent accurate mapping of such pleiotropic variants, which requires novel strategies to discover. After discovery, the selection strategy for deleterious variants under balancing selection is under debate, as variants can contribute to the heterosis effect in crossbred animals in various livestock species, compensating for the loss in purebred animals. Nevertheless, gene-assisted selection is a useful tool to decrease the frequency of the harmful allele in the population, if desired. Together, this review marks various deleterious variants under balancing selection and describing the functional consequences at the molecular, phenotypic, and population level, providing a resource for further study.

Investigating inbreeding in the turkey (Meleagris gallopavo) genome.

The detrimental effects of increased homozygosity due to inbreeding have prompted the development of methods to reduce inbreeding. The detection of runs of homozygosity (ROH), or contiguous stretches of homozygous marker genotypes, can be used to describe and quantify the level of inbreeding in an individual. The estimation of inbreeding coefficients can be calculated based on pedigree information, ROH, or the genomic relationship matrix. The aim of this study was to detect and describe ROH in the turkey genome and compare estimates of pedigree-based inbreeding coefficients (FPED) with genomic-based inbreeding coefficients estimated from ROH (FROH) and the genomic relationship matrix (FGRM). A total of 2,616,890 pedigree records were available. Of these records, 6,371 genotyped animals from three purebred turkey (Meleagris gallopavo) lines between 2013 and 2019 were available, and these were obtained using a dense single nucleotide polymorphism array (56,452 SNPs). The overall mean length of detected ROH was 2.87 ± 0.29 Mb with a mean number of 84.87 ± 8.79 ROH per animal. Short ROH with lengths of 1 to 2 Mb long were the most abundant throughout the genome. Mean ROH coverage differed greatly between chromosomes and lines. Considering inbreeding coefficient means across all lines, genomic derived inbreeding coefficients (FROH = 0.27; FGRM = 0.32) were higher than coefficients estimated from pedigree records (FPED = 0.14). Correlations between FROH and FPED, FROH and FGRM, and FPED and FGRM ranged between 0.19 to 0.31, 0.68 to 0.73, and 0.17 to 0.30, respectively. Additionally, correlations between FROH from different lengths and FPED substantially increased with ROH length from -0.06 to 0.33. Results of the current research, including the distribution of ROH throughout the genome and ROH-derived inbreeding estimates, can provide a more comprehensive description of inbreeding in the turkey genome. This knowledge can be used to evaluate genetic diversity, a requirement for genetic improvement, and develop methods to minimize inbreeding in turkey breeding programs.

Introgression contributes to distribution of structural variations in cattle.

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.

Accelerated discovery of functional genomic variation in pigs.

The genotype-phenotype link is a major research topic in the life sciences but remains highly complex to disentangle. Part of the complexity arises from the number of genes contributing to the observed phenotype. Despite the vast increase of molecular data, pinpointing the causal variant underlying a phenotype of interest is still challenging. In this study, we present an approach to map causal variation and molecular pathways underlying important phenotypes in pigs. We prioritize variation by utilizing and integrating predicted variant impact scores (pCADD), functional genomic information, and associated phenotypes in other mammalian species. We demonstrate the efficacy of our approach by reporting known and novel causal variants, of which many affect non-coding sequences. Our approach allows the disentangling of the biology behind important phenotypes by accelerating the discovery of novel causal variants and molecular mechanisms affecting important phenotypes in pigs. This information on molecular mechanisms could be applicable in other mammalian species, including humans.

Heterogeneity of a dwarf phenotype in Dutch traditional chicken breeds revealed by genomic analyses.

The growth of animals is a complex trait, in chicken resulting in a diverse variety of forms, caused by a heterogeneous genetic basis. Bantam chicken, known as an exquisite form of dwarfism, has been used for crossbreeding to create corresponding dwarf counterparts for native fowls in the Dutch populations. Here, we demonstrate the heterogeneity of the bantam trait in Dutch chickens and reveal the underlying genetic causes, using whole-genome sequence data from matching pairs of bantam and normal-sized breeds. During the bantam-oriented crossbreeding, various bantam origins were used to introduce the bantam phenotype, and three major bantam sources were identified and clustered. The genome-wide association studies revealed multiple genetic variants and genes associated with bantam phenotype, including HMGA2 and PRDM16, genes involved in body growth and stature. The comparison of associated variants among studies illustrated differences related to divergent bantam origins, suggesting a clear heterogeneity among bantam breeds. We show that in neo-bantam breeds, the bantam-related regions underwent a strong haplotype introgression from the bantam source, outcompeting haplotypes from the normal-sized counterpart. The bantam heterogeneity is further confirmed by the presence of multiple haplotypes comprising associated alleles, which suggests the selection of the bantam phenotype is likely subject to a convergent direction across populations. Our study demonstrates that the diverse history of human-mediated crossbreeding has contributed to the complexity and heterogeneity of the bantam phenotype.

Phenotypic and Genomic Analysis of Cystic Hygroma in Pigs.

Cystic hygroma is a malformation of the lymphatic and vascular system and is recognized as a benign congenital tumor that affects humans and animals in the perinatal period. This congeni-tal disorder is rarely described in animals, and until today, cystic hygroma in pigs has not been described in the literature. In a purebred Piètrain litter with twelve live-born piglets, cystic hy-groma was noticed on the rump of two male pigs within the first week of life. In addition, a third case of a crossbred weaner (Large White × Landrace) was detected during a herd examina-tion. To rule out common differential diagnoses, e.g., abscess or hematoma, further clinical and pathological investigations were conducted. During clinical examination, a painless and soft mass, which was compressible, was detected on the rump of all affected animals. The ultra-sonographic examination revealed a fluid-filled and cavernous subcutaneous structure. In addi-tion, a puncture of the cyst was conducted, revealing a serosanguinous fluid with negative bacte-riological culture. In all cases, a necropsy was performed, showing that the animals had fluid-filled cysts lined by well-differentiated lymphatic endothelium. Based on the clinicopathological examination, cystic hygroma was diagnosed. Furthermore, SNP array genotyping and whole-genome sequencing was performed and provided no evidence for a chromosomal disorder. In the Piètrain family, several genome regions were homozygous in both affected piglets. None-theless, a dominant acting de novo germline variant could not be ruled out, and therefore differ-ent filtering strategies were used to find pathogenic variants. The herein presented lists of pri-vate variants after filtering against hundreds of control genomes provide no plausible candidate and no shared variants among the two sequenced cases. Therefore, further studies are needed to evaluate possible genetic etiology. In general, systematic surveillance is needed to identify ge-netic defects as early as possible and to avoid the occurrence of losses in the pig population.

Using short read sequencing to characterise balanced reciprocal translocations in pigs.

BACKGROUND: A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions. RESULTS: Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes. CONCLUSIONS: We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.

Parallel Genetic Origin of Foot Feathering in Birds.

Understanding the genetic basis of similar phenotypes shared between lineages is a long-lasting research interest. Even though animal evolution offers many examples of parallelism, for many phenotypes little is known about the underlying genes and mutations. We here use a combination of whole-genome sequencing, expression analyses, and comparative genomics to study the parallel genetic origin of ptilopody (Pti) in chicken. Ptilopody (or foot feathering) is a polygenic trait that can be observed in domesticated and wild avian species and is characterized by the partial or complete development of feathers on the ankle and feet. In domesticated birds, ptilopody is easily selected to fixation, though extensive variation in the type and level of feather development is often observed. By means of a genome-wide association analysis, we identified two genomic regions associated with ptilopody. At one of the loci, we identified a 17-kb deletion affecting PITX1 expression, a gene known to encode a transcription regulator of hindlimb identity and development. Similarly to pigeon, at the second loci, we observed ectopic expression of TBX5, a gene involved in forelimb identity and a key determinant of foot feather development. We also observed that the trait evolved only once as foot-feathered birds share the same haplotype upstream TBX5. Our findings indicate that in chicken and pigeon ptilopody is determined by the same set of genes that affect similar molecular pathways. Our study confirms that ptilopody has evolved through parallel evolution in chicken and pigeon.