Genomic Comparison of Indigenous African and Northern European Chickens Reveals Putative Mechanisms of Stress Tolerance Related to Environmental Selection Pressure

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.

Climate-driven mitochondrial selection in lacertid lizards

The mitochondrion, which is an intracellular organelle responsible for most of the energy-producing pathways, can have its genome targeted for climate-driven selection. However, climate-driven mitochondrial selection remains a sparsely studied area in reptiles. Here, we reported the complete mitochondrial genome sequence of a lacertid lizard (Takydromus intermedius) and used mitogenomes from 54 species of lacertid lizards to study their phylogenetic relationships and to identify the mitochondrial genes under positive selection by climate. The length of the complete mitochondrial genome sequence of T. intermedius was 17,713 bp, which was within the range of lengths (17,224-18,943) ever reported for Takydromus species. The arrangement of mitochondrial genes in T. intermedius was the same as in other congeneric species. The 54 lacertid species could be divided into three geographically and climatically different clades. We identified three mitochondrial genes (ATP6, ATP8, and ND3) under positive selection by climate, and found that isothermality, temperature seasonality, precipitation of wettest month, and precipitation seasonality were the most important climatic variables contributing to the gene selection.

Genome-wide scan for runs of homozygosity in South American Camelids

BACKGROUND

Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection.

RESULTS

A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively.

CONCLUSIONS

The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources

Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.

Gene expression networks and functionally enriched pathways involved in the response of domestic chicken to acute heat stress

Heat stress in poultry houses, especially in warm areas, is one of the main environmental factors that restrict the growth of broilers or laying performance of layers, suppresses the immune system, and deteriorates egg quality and feed conversion ratio. The molecular mechanisms underlying the response of chicken to acute heat stress (AHS) have not been comprehensively elucidated. Therefore, the main object of the current work was to investigate the liver gene expression profile of chickens under AHS in comparison with their corresponding control groups, using four RNA-seq datasets. The meta-analysis, GO and KEGG pathway enrichment, WGCNA, machine-learning, and eGWAS analyses were performed. The results revealed 77 meta-genes that were mainly related to protein biosynthesis, protein folding, and protein transport between cellular organelles. In other words, under AHS, the expression of genes involving in the structure of rough reticulum membrane and in the process of protein folding was adversely influenced. In addition, genes related to biological processes such as "response to unfolded proteins," "response to reticulum stress" and "ERAD pathway" were differentially regulated. We introduce here a couple of genes such as HSPA5, SSR1, SDF2L1, and SEC23B, as the most significantly differentiated under AHS, which could be used as bio-signatures of AHS. Besides the mentioned genes, the main findings of the current work may shed light to the identification of the effects of AHS on gene expression profiling of domestic chicken as well as the adaptive response of chicken to environmental stresses.

Three-Dimensional Organization of Chicken Genome Provides Insights into Genetic Adaptation to Extreme Environments

The high-throughput chromosome conformation capture (Hi-C) technique is widely used to study the functional roles of the three-dimensional (3D) architecture of genomes. However, the knowledge of the 3D genome structure and its dynamics during extreme environmental adaptations remains poor. Here, we characterized 3D genome architectures using the Hi-C technique for chicken liver cells. Upon comparing Lindian chicken (LDC) liver cells with Wenchang chicken (WCC) liver cells, we discovered that environmental adaptation contributed to the switching of A/B compartments, the reorganization of topologically associated domains (TADs), and TAD boundaries in both liver cells. In addition, the analysis of the switching of A/B compartments revealed that the switched compartmental genes (SCGs) were strongly associated with extreme environment adaption-related pathways, including tight junction, notch signaling pathway, vascular smooth muscle contraction, and the RIG-I-like receptor signaling pathway. The findings of this study advanced our understanding of the evolutionary role of chicken 3D genome architecture and its significance in genome activity and transcriptional regulation.

Whole Genome Analyses Reveal Novel Genes Associated with Chicken Adaptation to Tropical and Frigid Environments

INTRODUCTION

Investigating the genetic footprints of historical temperature selection can get insights to the local adaptation and feasible influences of climate change on long-term population dynamics.

OBJECT

Chicken is a significative species to study genetic adaptation on account of its similar domestication track related to human activity with the most diversified varieties. Yet, few studies have demonstrated the genetic signatures of its adaptation to naturally tropical and frigid environments.

METHOD

Here, we generated whole genome resequencing of 119 domesticated chickens in China including the following breeds which are in order of breeding environmental temperature from more tropical to more frigid: Wenchang chicken (WCC), green-shell chicken (GSC), Tibetan chicken (TBC), and Lindian chicken (LDC).

RESULTS

Our results showed WCC branched off earlier than LDC with an evident genetic admixture between WCC and LDC, suggesting their closer genetic relationship. Further comparative genomic analyses solute carrier family 33 member 1 (SLC33A1) and thyroid stimulating hormone receptor (TSHR) genes exhibited stronger signatures for positive selection in the genome of the more tropical WCC. Furthermore, genotype data from about 3,000 African local ecotypes confirmed that allele frequencies of single nucleotide polymorphisms (SNPs) in these 2 genes appeared strongly associated with tropical environment adaptation. In addition, the NADH:ubiquinone oxidoreductase subunit S4 (NDUFS4) gene exhibited a strong signature for positive selection in the LDC genome, and SNPs with marked allele frequency differences indicated a significant relationship with frigid environment adaptation.

CONCLUSION

Our findings partially clarify how selection footprints from environmental temperature stress can lead to advantageous genomic adaptions to tropical and frigid environments in poultry and provide a valuable resource for selective breeding of chickens.

Genetic diversity and signatures of selection for heat tolerance and immune response in Iranian native chickens

Background

Understanding how evolutionary forces relating to climate have shaped the patterns of genetic variation within and between species is a fundamental pursuit in biology. Iranian indigenous chickens have evolved genetic adaptations to their local environmental conditions, such as hot and arid regions. In the present study, we provide a population genome landscape of genetic variations in 72 chickens representing nine Iranian indigenous ecotypes (Creeper, Isfahan, Lari, Marand, Mashhad, Naked neck, Sari, Shiraz and Yazd) and two commercial lines (White Leghorn and Arian). We further performed comparative population genomics to evaluate the genetic basis underlying variation in the adaptation to hot climate and immune response in indigenous chicken ecotypes. To detect genomic signatures of adaptation, we applied nucleotide diversity (θπ) and F_ST statistical measurements, and further analyzed the results to find genomic regions under selection for hot adaptation and immune response-related traits.

Results

By generating whole-genome data, we assessed the relationship between the genetic diversity of indigenous chicken ecotypes and their genetic distances to two different commercial lines. The results of genetic structure analysis revealed clustering of indigenous chickens in agreement with their geographic origin. Among all studied chicken groups, the highest level of linkage disequilibrium (LD) (~ 0.70) was observed in White Leghorn group at marker pairs distance of 1 Kb. The results from admixture analysis demonstrated evidence of shared ancestry between Arian individuals and indigenous chickens, especially those from the north of the country. Our search for potential genomic regions under selection in indigenous chicken ecotypes revealed several immune response and heat shock protein-related genes, such as HSP70, HSPA9, HSPH1, HSP90AB1 and PLCB4 that have been previously unknown to be involved in environmental-adaptive traits. In addition, we found some other candidate loci on different chromosomes probably related with hot adaptation and immune response-related traits.

Conclusions

The work provides crucial insights into the structural variation in the genome of Iranian indigenous chicken ecotypes, which up to now has not been genetically investigated. Several genes were identified as candidates for drought, heat tolerance, immune response and other phenotypic traits. These candidate genes may be helpful targets for understanding of the molecular basis of adaptation to hot environmental climate and as such they should be used in chicken breeding programs to select more efficient breeds for desert climate.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08434-7.

Selection signatures for heat tolerance in Brazilian horse breeds

Since domestication, horse breeds have adapted to their environments and differentiated from one another. This paper uses two methods to detect selection signatures in 23 horse breeds, eight of which are Brazilian (610 animals), both cold-blooded and warm-blooded, from temperate and tropical regions. These animals were genotyped using the GGP Equine BeadChip and we analysed the data by Principal Component Analysis (PCA). The samples were separated into groups based on their geographical area of origin and PCA results studied. The genomic regions under selection were detected by hapFLK and PCAdapt methodologies, identifying six regions under selection with at least one Brazilian horse breed. These regions contain genes associated with heat tolerance, skin colour, body size, energy production/metabolism, genes involved in protein degradation/turnover/DNA repair, genes reducing the impact of oxidative stress/cellular repair, and transcriptional regulation. This work confirmed LCORL and NCAPG gene regions in previous studies associated with body size on Equine Chromosome Autosome 3 (ECA3). On the same ECA3, a region implicating genes linked to coat colour was identified, also previously related to heat stress. Regions with genes coding heat shock proteins were found on ECA1 and 2, and many candidate genes for oxidation-reduction which are a natural response to heat stress. However, a larger sample size and whole-genome SNPs are needed to understand better and identify new candidate regions as well as their functional relation with heat tolerance.

Runs of homozygosity analysis reveals population characteristics of yellow-feathered chickens using re-sequencing data

1. To evaluate the inbreeding of yellow-feathered chickens (YFCs) and identify genes related to their unique characteristics, whole-genome re-sequencing data were applied to detect runs of homozygosity (ROH) in the genomes of ten YFC breeds. The number, length, distribution of ROH, and inbreeding coefficient in different YFC populations were calculated. Genomic regions with high frequency in ROH were annotated.2. In total, 25,547 ROH with an average length of 335 kb were detected, with most being <1 Mb. The domination of short ROH reflected the long breeding history of this chicken. The number, length, frequency, and distribution of ROH varied among chicken populations, and high genetic diversity was maintained.3. Numerous genes related to YFC characteristics were identified in the high-frequency ROH regions. Among these, IFNA, IFNB, IL11RA, IL22RA1, IFNLR1, and TRIF genes were involved in disease resistance. The AMY, G6PC, SDHB, GCNT4, and ACO genes were associated with energy material metabolism; and FABPL, AQP7, ACAA2, and RYR2 were related to meat quality and flavour. The KITLG, CREB3, RYR2, and LGR4 genes, related to pigmentation, were detected.4. This ROH-based inbreeding evaluation laid the foundation for breeding and conservation of YFC populations, and the candidate genes identified can be used for marker-assisted selection.

A genome-wide scan to identify signatures of selection in two Iranian indigenous chicken ecotypes

Background

Various regions of the chicken genome have been under natural and artificial selection for thousands of years. The substantial diversity that exits among chickens from different geographic regions provides an excellent opportunity to investigate the genomic regions under selection which, in turn, will increase our knowledge about the mechanisms that underlie chicken diversity and adaptation. Several statistics have been developed to detect genomic regions that are under selection. In this study, we applied approaches based on differences in allele or haplotype frequencies (F_ST and hapFLK, respectively) between populations, differences in long stretches of consecutive homozygous sequences (ROH), and differences in allele frequencies within populations (composite likelihood ratio (CLR)) to identify inter- and intra-populations traces of selection in two Iranian indigenous chicken ecotypes, the Lari fighting chicken and the Khazak or creeper (short-leg) chicken.

Results

Using whole-genome resequencing data of 32 individuals from the two chicken ecotypes, approximately 11.9 million single nucleotide polymorphisms (SNPs) were detected and used in genomic analyses after quality processing. Examination of the distribution of ROH in the two populations indicated short to long ROH, ranging from 0.3 to 5.4 Mb. We found 90 genes that were detected by at least two of the four applied methods. Gene annotation of the detected putative regions under selection revealed candidate genes associated with growth (DCN, MEOX2 and CACNB1), reproduction (ESR1 and CALCR), disease resistance (S1PR1, ALPK1 and MHC-B), behavior pattern (AGMO, GNAO1 and PSEN1), and morphological traits (IHH and NHEJ1).

Conclusions

Our findings show that these two phenotypically different indigenous chicken populations have been under selection for reproduction, immune, behavioral, and morphology traits. The results illustrate that selection can play an important role in shaping signatures of differentiation across the genomic landscape of two chicken populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00664-9.

Species and Phenotypic Distribution Models Reveal Population Differentiation in Ethiopian Indigenous Chickens

Smallholder poultry production dominated by indigenous chickens is an important source of livelihoods for most rural households in Ethiopia. The long history of domestication and the presence of diverse agroecologies in Ethiopia create unique opportunities to study the effect of environmental selective pressures. Species distribution models (SDMs) and Phenotypic distribution models (PDMs) can be applied to investigate the relationship between environmental variation and phenotypic differentiation in wild animals and domestic populations. In the present study we used SDMs and PDMs to detect environmental variables related with habitat suitability and phenotypic differentiation among nondescript Ethiopian indigenous chicken populations. 34 environmental variables (climatic, soil, and vegetation) and 19 quantitative traits were analyzed for 513 adult chickens from 26 populations. To have high variation in the dataset for phenotypic and ecological parameters, animals were sampled from four spatial gradients (each represented by six to seven populations), located in different climatic zones and geographies. Three different ecotypes are proposed based on correlation test between habitat suitability maps and phenotypic clustering of sample populations. These specific ecotypes show phenotypic differentiation, likely in response to environmental selective pressures. Nine environmental variables with the highest contribution to habitat suitability are identified. The relationship between quantitative traits and a few of the environmental variables associated with habitat suitability is non-linear. Our results highlight the benefits of integrating species and phenotypic distribution modeling approaches in characterization of livestock populations, delineation of suitable habitats for specific breeds, and understanding of the relationship between ecological variables and quantitative traits, and underlying evolutionary processes.

Genome-Wide Analysis Reveals Selection Signatures Involved in Meat Traits and Local Adaptation in Semi-Feral Maremmana Cattle

The Maremmana cattle is an ancient Podolian-derived Italian breed raised in semi-wild conditions with distinctive morphological and adaptive traits. The aim of this study was to detect potential selection signatures in Maremmana using medium-density single nucleotide polymorphism array. Putative selection signatures were investigated combining three statistical approaches designed to quantify the excess of haplotype homozygosity either within (integrated haplotype score, iHS) or among pairs of populations (Rsb and XP-EHH), and contrasting the Maremmana with a single reference population composed of a pool of seven Podolian-derived Italian breeds. Overall, the three haplotype-based analyses revealed selection signatures distributed over 19 genomic regions. Of these, six relevant candidate regions were identified by at least two approaches. We found genomic signatures of selective sweeps spanning genes related to mitochondrial function, muscle development, growth, and meat traits (SCIN, THSD7A, ETV1, UCHL1, and MYOD1), which reflects the different breeding schemes between Maremmana (semi-wild conditions) and the other Podolian-derived Italian breeds (semi-extensive). We also identified several genes linked to Maremmana adaptation to the environment of the western-central part of Italy, known to be hyperendemic for malaria and other tick-borne diseases. These include several chemokine (C-C motif) ligand genes crucially involved in both innate and adaptive immune responses to intracellular parasite infections and other genes playing key roles in pulmonary disease (HEATR9, MMP28, and ASIC2) or strongly associated with malaria resistance/susceptibility (AP2B1). Our results provide a glimpse into diverse selection signatures in Maremmana cattle and can be used to enhance our understanding of the genomic basis of environmental adaptation in cattle.

Emerging Genetic Tools to Investigate Molecular Pathways Related to Heat Stress in Chickens: A Review

Simple Summary

New genomic tools have been used as an instrument in order to assess the molecular pathway involved in heat stress resistance. Local chicken breeds have a better attitude to face heat stress. This review aims to summarize studies linked to chickens, heat stress, and heat shock protein.

Abstract

Chicken products are the most consumed animal-sourced foods at a global level across greatly diverse cultures, traditions, and religions. The consumption of chicken meat has increased rapidly in the past few decades and chicken meat is the main animal protein source in developing countries. Heat stress is one of the environmental factors which decreases the productive performance of poultry and meat quality. Heat stress produces the over-expression of heat shock factors and heat shock proteins in chicken tissues. Heat shock proteins regulate several molecular pathways in cells in response to stress conditions, changing the homeostasis of cells and tissues. These changes can affect the physiology of the tissue and hence the production ability of chickens. Indeed, commercial chicken strains can reach a high production level, but their body metabolism, being comparatively accelerated, has poor thermoregulation. In contrast, native backyard chickens are more adapted to the environments in which they live, with a robustness that allows them to survive and reproduce constantly. In the past few years, new molecular tools have been developed, such as RNA-Seq, Single Nucleotide Polymorphisms (SNPs), and bioinformatics approaches such as Genome-Wide Association Study (GWAS). Based on these genetic tools, many studies have detected the main pathways involved in cellular response mechanisms. In this context, it is necessary to clarify all the genetic and molecular mechanisms involved in heat stress response. Hence, this paper aims to review the ability of the new generation of genetic tools to clarify the molecular pathways associated with heat stress in chickens, offering new perspectives for the use of these findings in the animal breeding field.

Epigenetics of Animal Personality: DNA Methylation Cannot Explain the Heritability of Exploratory Behavior in a Songbird

The search for the hereditary mechanisms underlying quantitative traits traditionally focused on the identification of underlying genomic polymorphisms such as single-nucleotide polymorphisms. It has now become clear that epigenetic mechanisms, such as DNA methylation, can consistently alter gene expression over multiple generations. It is unclear, however, if and how DNA methylation can stably be transferred from one generation to the next and can thereby be a component of the heritable variation of a trait. In this study, we explore whether DNA methylation responds to phenotypic selection using whole-genome and genome-wide bisulfite approaches. We assessed differential erythrocyte DNA methylation patterns between extreme personality types in the Great Tit (Parus major). For this, we used individuals from a four-generation artificial bi-directional selection experiment and siblings from eight F2 inter-cross families. We find no differentially methylated sites when comparing the selected personality lines, providing no evidence for the so-called epialleles associated with exploratory behavior. Using a pair-wise sibling design in the F2 intercrosses, we show that the genome-wide DNA methylation profiles of individuals are mainly explained by family structure, indicating that the majority of variation in DNA methylation in CpG sites between individuals can be explained by genetic differences. Although we found some candidates explaining behavioral differences between F2 siblings, we could not confirm this with a whole-genome approach, thereby confirming the absence of epialleles in these F2 intercrosses. We conclude that while epigenetic variation may underlie phenotypic variation in behavioral traits, we were not able to find evidence that DNA methylation can explain heritable variation in personality traits in Great Tits.

Analyses of mitogenomic markers shed light on the divergence, population dynamics, and demographic history of Pakistani chickens

Pakistan is one of a few sites, associated with the earliest known independent domestication event in the evolutionary history of chicken, which is socio-economically and historically the most important poultry bird in the country. However, the divergence, past population dynamics, and demographic history of Pakistani chickens have not been addressed so far. Therefore, we herein investigated the indigenous Pakistani chickens using mitogenomic markers. We first prepared individual DNA samples from the chicken feathers, and generated nucleotide sequence data, which was then subjected to various population genetics analyses. In molecular phylogenetic analysis, the Pakistani chickens were clustered under nine different clades. Among the wild fowls, the Indian red jungle fowl (IRJF) shared very close affinities to Pakistani chickens. The Bayesian skyline plot showed an increase in the effective population size of Pakistani chickens during the last 50 years. Finally, a time-calibrated phylogeny inferred molecular divergence of the Pakistani chickens. A molecular rate of 3.6 × 10^-6 mutations/site/year (95% HPD interval: 2.28 × 10^-8 to 9.32 × 10^-6) was estimated for the data set. In a rooted tree with root-age of 12058 years (95% HPD interval: 1161-38411), the Pakistani chicken haplotypes showed divergence from IRJF haplotypes around 6987 years (95% HPD interval: 1132-20746) ago, and they shared their most recent common ancestor with Gallus gallus spadiceus, and G. g. jabouillei at the root of the tree. Overall, these results suggest that Pakistani chicken haplotypes share their ancestral gene pool with the IRJF as compared to other red jungle fowl subspecies.

Avian Physiology: Are Birds Simply Feathered Mammals?

There are marked differences between the physiology of birds and mammals. These reflect the evolutionary distance between the two classes with the last common ancestor estimated as existing 318 million years ago. There are analogous organ systems in birds and mammals. However, marked differences exist. For instance, in the avian gastro-intestinal tract, there is a crop at the lower end of the esophagus. This functions both to store feed and for microbial action. The avian immune system lacks lymph nodes and has a distinct organ producing B-lymphocytes, namely the bursa Fabricius. The important of spleen has been largely dismissed until recently. However, its importance in both innate and specific immunity is increasingly recognized. There is a major difference between birds and mammals is the female reproductive system as birds produce large yolk filled eggs. The precursors of the yolk are synthesized by the liver. Another difference is that there is a single ovary and oviduct in birds.

Phenotypic variability and population structure analysis of Tanzanian free-range local chickens

BACKGROUND

Free-range local chickens (FRLC) farming is an important activity in Tanzania, however, they have not been well-characterized. This study aimed to phenotypically characterize three Tanzanian FRLCs and to determine their population structure. A total of 389 mature breeder chickens (324 females and 65 males) from three popular Tanzanian FRLC ecotypes (Kuchi, Morogoro-medium and Ching'wekwe) were used for the phenotypic characterization. Progenies of these chickens were utilized to assess population structure. The ecotypes were collected from four geographical zones across Tanzania: Lake, Central, Northern and Coastal zones. Body weights and linear measurements were obtained from the mature breeders, including body, neck, shanks, wingspan, chest girth, and shank girth. Descriptive statistics were utilized to characterize the chickens. Correlations between the linear measurements and differences among the means of measured linear traits between ecotypes and between sexes were assessed. A total of 1399 progeny chicks were genotyped using a chicken 600 K high density single nucleotide polymorphism (SNP) panel for determination of population structure.

RESULTS

The means for most traits were significantly higher in Kuchi relative to Ching'wekwe and Morogoro-medium. However, shank length and shank girth were similar between Kuchi and Morogoro-medium females. All traits were correlated with the exception of shank girth in Morogoro-medium. Admixture analyses revealed that Morogoro-medium and Ching'wekwe clustered together as one population, separate from Kuchi.

CONCLUSIONS

Phenotypic traits could be used to characterize FRLCs, however, there were variations in traits among individuals within ecotypes; therefore, complementary genomic methods should be considered to improve the characterization for selective breeding.

Genetic Diversity and Selection Signatures Within Diannan Small-Ear Pigs Revealed by Next-Generation Sequencing

Genetic characterization of Chinese indigenous pig breeds is essential to promote scientific conservation and sustainable development of pigs. Here, we systematically surveyed the genomes of 75 unrelated Diannan small-ear (DSE) pigs from three diverse regions (Yingjiang County, Jinping County, and Sipsongpanna in Yunnan Province) to describe their population structures, genetic diversity, inbreeding coefficients, and selection signatures. First, these individuals were sequenced and genotyped using the genome reducing and sequencing (GGRS) protocol. A total of 438,038 autosomal single-nucleotide polymorphisms (SNPs) were obtained and used for subsequent statistical analysis. The results showed that these DSE pigs were clearly differentiated into three separate clades revealed by the population structure and principal component analysis, which is consistent with their geographical origins. Diannan small-ear pigs owned lower genetic diversity when compared with some other pig breeds, which demonstrated the need to strengthen the conservation strategies for DSE pigs. In addition, the inbreeding coefficients based on runs of homozygosity (ROH) length (F _ROH) were calculated in each ROH length categories, respectively. And the results indicated that the ancient (up to 50 generations ago) inbreeding had greater impacts than recent (within the last five generations) inbreeding within DSE pigs. Some candidate selection signatures within the DSE pig population were detected through the ROH islands and integrated haplotype homozygosity score (iHS) methods. And genes associated with meat quality (COL15A1, RPL3L, and SLC9A3R2), body size (PALM2-AKAP2, NANS, TRAF7, and PACSIN1), adaptability (CLDN9 and E4F1), and appetite (GRM4) were identified. These findings can help to understand the genetic characteristics and provide insights into the molecular background of special phenotypes of DSE pigs to promote conservation and sustainability of the breed.

Genetic Basis of Response of Ghanaian Local Chickens to Infection With a Lentogenic Newcastle Disease Virus

Newcastle disease (ND) is a global threat to domestic poultry, especially in rural areas of Africa and Asia, where the loss of entire backyard local chicken flocks often threatens household food security and income. To investigate the genetics of Ghanaian local chicken ecotypes to Newcastle disease virus (NDV), in this study, three popular Ghanaian chicken ecotypes (regional populations) were challenged with a lentogenic NDV strain at 28 days of age. This study was conducted in parallel with a similar study that used three popular Tanzanian local chicken ecotypes and after two companion studies in the United States, using Hy-line Brown commercial laying birds. In addition to growth rate, NDV response traits were measured following infection, including anti-NDV antibody levels [pre-infection and 10 days post-infection (dpi)], and viral load (2 and 6 dpi). Genetic parameters were estimated, and two genome-wide association study analysis methods were used on data from 1,440 Ghanaian chickens that were genotyped on a chicken 600K Single Nucleotide Polymorphism (SNP) chip. Both Ghana and Tanzania NDV challenge studies revealed moderate to high (0.18 – 0.55) estimates of heritability for all traits, except viral clearance where the heritability estimate was not different from zero for the Tanzanian ecotypes. For the Ghana study, 12 quantitative trait loci (QTL) for growth and/or response to NDV from single-SNP analyses and 20 genomic regions that explained more than 1% of genetic variance using the Bayes B method were identified. Seven of these windows were also identified as having at least one significant SNP in the single SNP analyses for growth rate, anti-NDV antibody levels, and viral load at 2 and 6 dpi. An important gene for growth during stress, CHORDC1 associated with post-infection growth rate was identified as a positional candidate gene, as well as other immune related genes, including VAV2, IL12B, DUSP1, and IL17B. The QTL identified in the Ghana study did not overlap with those identified in the Tanzania study. However, both studies revealed QTL with genes vital for growth and immune response during NDV challenge. The Tanzania parallel study revealed an overlapping QTL on chromosome 24 for viral load at 6 dpi with the US NDV study in which birds were challenged with NDV under heat stress. This QTL region includes genes related to immune response, including TIRAP, ETS1, and KIRREL3. The moderate to high estimates of heritability and the identified QTL suggest that host response to NDV of local African chicken ecotypes can be improved through selective breeding to enhance increased NDV resistance and vaccine efficacy.

Genomic regions and pathways associated with thermotolerance in layer-type strain Taiwan indigenous chickens

This study aimed to investigate genetic markers and candidate genes associated with thermotolerance in a layer-type strain Taiwan indigenous chickens exposed to acute heat stress. One hundred and ninety-two 30-week-old roosters were subjected to acute heat stress. Changes in body temperature (BT, ΔT) were calculated by measuring the difference between the initial BT and the highest BT during heat stress and the results were categorized into dead, susceptible, tolerant, and intermediate groups depending on their survival and ΔT values at the end of the experiment. A genome-wide association study on survival and ΔT values was conducted using the Cochran-Armitage trend test and Fisher's exact test. Association analyses identified 80 significant SNPs being annotated to 23 candidate genes, 440 SNPs to 71 candidate genes, 64 SNPs to 25 candidate genes, and 378 SNPs to 78 candidate genes in the dead versus survivor, tolerant versus susceptible, intermediate versus tolerant, and intermediate versus susceptible groups, respectively. The annotated genes were associated with apoptosis, cellular stress responses, DNA repair, and metabolic oxidative stress. In conclusion, the identified SNPs of candidate genes provide insights into the potential mechanisms underlying physiological responses to acute heat stress in chickens.

Cardiac and Skeletal Muscle Transcriptome Response to Heat Stress in Kenyan Chicken Ecotypes Adapted to Low and High Altitudes Reveal Differences in Thermal Tolerance and Stress Response

Heat stress (HS) negatively affects chicken performance. Agricultural expansion will happen in regions that experience high ambient temperatures, where fast-growing commercial chickens are vulnerable. Indigenous chickens of such regions, due to generations of exposure to environmental challenges, might have higher thermal tolerance. In this study, two indigenous chicken ecotypes, from the hot and humid Mombasa (lowland) and the colder Naivasha (highland) regions, were used to investigate the effects of acute (5 h, 35°C) and chronic (3 days of 35°C for 8 h/day) HS on the cardiac and skeletal muscle, through RNA sequencing. The rectal temperature gain and the number of differentially expressed genes (DEGs) [False Discovery Rate (FDR) < 0.05] were two times higher in the acute stage than in the chronic stage in both ecotypes, suggesting that cyclic exposure to HS can lead to adaptation. A tissue- and stage-specific difference in response to HS was observed, with peroxisome proliferator-activated-receptor (PPAR) signaling and mitogen-activate protein kinase (MAPK) signaling pathways, enriched in heart and skeletal muscle, respectively, and the p53 pathway enriched only in the acute stage in both tissues. The acute and chronic stage DEGs were integrated by a region-specific gene coexpression network (GCN), and genes with the highest number of connections (hub genes) were identified. The hub genes in the lowland network were CCNB2, Crb2, CHST9, SESN1, and NR4A3, while COMMD4, TTC32, H1F0, ACYP1, and RPS28 were the hub genes in the highland network. Pathway analysis of genes in the GCN showed that p53 and PPAR signaling pathways were enriched in both low and highland networks, while MAPK signaling and protein processing in endoplasmic reticulum were enriched only in the gene network of highland chickens. This shows that to dissipate the accumulated heat, to reduce heat induced apoptosis, and to promote DNA damage repair, the ecotypes activated or suppressed different genes, indicating the differences in thermal tolerance and HS response mechanisms between the ecotypes. This study provides information on the HS response of chickens, adapted to two different agro climatic environments, extending our understanding of the mechanisms of HS response and the effect of adaptation in counteracting HS.

Natural Selection Footprints Among African Chicken Breeds and Village Ecotypes

Natural selection is likely a major factor in shaping genomic variation of the African indigenous rural chicken, driving the development of genetic footprints. Selection footprints are expected to be associated with adaptation to locally prevailing environmental stressors, which may include diverse factors as high altitude, disease resistance, poor nutrition, oxidative and heat stresses. To determine the existence of a selection footprint, 268 birds were randomly sampled from three indigenous ecotypes from East Africa (Rwanda and Uganda) and North Africa (Baladi), and two registered Egyptian breeds (Dandarawi and Fayoumi). Samples were genotyped using the chicken Affymetrix 600K Axiom^® Array. A total of 494,332 SNPs were utilized in the downstream analysis after implementing quality control measures. The intra-population runs of homozygosity (ROH) that occurred in >50% of individuals of an ecotype or in >75% of a breed were studied. To identify inter-population differentiation due to genetic structure, F_ST was calculated for North- vs. East-African populations and Baladi and Fayoumi vs. Dandarawi for overlapping windows (500 kb with a step-size of 250 kb). The ROH and F_ST mapping detected several selective sweeps on different autosomes. Results reflected selection footprints of the environmental stresses, breed behavior, and management. Intra-population ROH of the Egyptian chickens showed selection footprints bearing genes for adaptation to heat, solar radiation, ion transport and immunity. The high-altitude-adapted East-African populations' ROH showed a selection signature with genes for angiogenesis, oxygen-heme binding and transport. The neuroglobin gene (GO:0019825 and GO:0015671) was detected on a Chromosome 5 ROH of Rwanda-Uganda ecotypes. The sodium-dependent noradrenaline transporter, SLC6A2 on a Chromosome 11 ROH in Fayoumi breed may reflect its active behavior. Inter-population F_ST among Egyptian populations reflected genetic mechanisms for the Fayoumi resistance to Newcastle Disease Virus (NDV), while F_ST between Egyptian and Rwanda-Uganda populations indicated the Secreted frizzled related protein 2, SFRP2, (GO:0009314) on Chromosome 4, that contributes to melanogenic activity and most likely enhances the Dandarawi chicken adaptation to high-intensity of solar radiation in Southern Egypt. These results enhance our understanding of the natural selection forces role in shaping genomic structure for adaptation to the stressful African conditions.

Leveraging Available Resources and Stakeholder Involvement for Improved Productivity of African Livestock in the Era of Genomic Breeding

The African continent is home to diverse populations of livestock breeds adapted to harsh environmental conditions with more than 70% under traditional systems of management. Animal productivity is less than optimal in most cases and is faced with numerous challenges including limited access to adequate nutrition and disease management, poor institutional capacities and lack of adequate government policies and funding to develop the livestock sector. Africa is home to about 1.3 billion people and with increasing demand for animal proteins by an ever growing human population, the current state of livestock productivity creates a significant yield gap for animal products. Although a greater section of the population, especially those living in rural areas depend largely on livestock for their livelihoods; the potential of the sector remains underutilized and therefore unable to contribute significantly to economic development and social wellbeing of the people. With current advances in livestock management practices, breeding technologies and health management, and with inclusion of all stakeholders, African livestock populations can be sustainably developed to close the animal protein gap that exists in the continent. In particular, advances in gene technologies, and application of genomic breeding in many Western countries has resulted in tremendous gains in traits like milk production with the potential that, implementation of genomic selection and other improved practices (nutrition, healthcare, etc.) can lead to rapid improvement in traits of economic importance in African livestock populations. The African livestock populations in the context of this review are limited to cattle, goat, pig, poultry, and sheep, which are mainly exploited for meat, milk, and eggs. This review examines the current state of livestock productivity in Africa, the main challenges faced by the sector, the role of various stakeholders and discusses in-depth strategies that can enable the application of genomic technologies for rapid improvement of livestock traits of economic importance.

Detection of Selection Signatures Among Brazilian, Sri Lankan, and Egyptian Chicken Populations Under Different Environmental Conditions

Extreme environmental conditions are a major challenge in livestock production. Changes in climate, particularly those that contribute to weather extremes like drought or excessive humidity, may result in reduced performance and reproduction and could compromise the animal's immune function. Animal survival within extreme environmental conditions could be in response to natural selection and to artificial selection for production traits that over time together may leave selection signatures in the genome. The aim of this study was to identify selection signatures that may be involved in the adaptation of indigenous chickens from two different climatic regions (Sri Lanka = Tropical; Egypt = Arid) and in non-indigenous chickens that derived from human migration events to the generally tropical State of São Paulo, Brazil. To do so, analyses were conducted using fixation index (Fst) and hapFLK analyses. Chickens from Brazil (n = 156), Sri Lanka (n = 92), and Egypt (n = 96) were genotyped using the Affymetrix Axiom^®600k Chicken Genotyping Array. Pairwise Fst analyses among countries did not detect major regions of divergence between chickens from Sri Lanka and Brazil, with ecotypes/breeds from Brazil appearing to be genetically related to Asian-Indian (Sri Lanka) ecotypes. However, several differences were detected in comparisons of Egyptian with either Sri Lankan or Brazilian populations, and common regions of difference on chromosomes 2, 3 and 8 were detected. The hapFLK analyses for the three separate countries suggested unique regions that are potentially under selection on chromosome 1 for all three countries, on chromosome 4 for Sri Lankan, and on chromosomes 3, 5, and 11 for the Egyptian populations. Some of identified regions under selection with hapFLK analyses contained genes such as TLR3, SOCS2, EOMES, and NFAT5 whose biological functions could provide insights in understanding adaptation mechanisms in response to arid and tropical environments.

Looking at genetic structure and selection signatures of the Mexican chicken population using single nucleotide polymorphism markers

Genetic variation enables both adaptive evolutionary changes and artificial selection. Genetic makeup of populations is the result of a long-term process of selection and adaptation to specific environments and ecosystems. The aim of this study was to characterize the genetic variability of México's chicken population to reveal any underlying population structure. A total of 213 chickens were sampled in different rural production units located in 25 states of México. Genotypes were obtained using the Affymetrix Axiom® 600 K Chicken Genotyping Array. The Identity by Descent (IBD) and the principal components analysis (PCA) were performed by SVS software on pruned single nucleotide polymorphisms (SNPs).ADMIXTURE analyses identified 3 ancestors and the proportion of the genetic contribution of each of them has been determined in each individual. The results of the Neighbor-Joining (NJ) analysis resulted consistent with those obtained by the PCA. All methods utilized in this study did not allow a classification of Mexican chicken in distinct clusters or groups. A total of 3,059 run of homozygosity (ROH) were identified and, being mainly short in length (<4 Mb), these regions are indicative of a low inbreeding level in the population. Finally, findings from the ROH analysis indicated the presence of natural selective pressure in the population of Mexican chicken.The study indicates that the Mexican chicken clearly appear to be a unique creole chicken population that was not subjected to a specific artificial selection. Results provide a genetic knowledge that can be used as a basis for the genetic management of a unique and very large creole population, especially in the view of using it in production of hybrids to increase the productivity and economic revenue of family farming agriculture, which is widely present in México.