Genetic diversity of 10 indigenous chicken ecotypes from Southern Highlands of Tanzania based on Major Histocompatibility Complex-linked microsatellite LEI0258 marker typing

Effectiveness of DArTseq markers application in genetic diversity and population structure of indigenous chickens in Eastern Province of Rwanda

BACKGROUND

The application of biotechnologies which make use of genetic markers in chicken breeding is developing rapidly. Diversity Array Technology (DArT) is one of the current Genotyping-By-Sequencing techniques allowing the discovery of whole genome sequencing. In livestock, DArT has been applied in cattle, sheep, and horses. Currently, there is no study on the application of DArT markers in chickens. The aim was to study the effectiveness of DArTSeq markers in the genetic diversity and population structure of indigenous chickens (IC) and SASSO in the Eastern Province of Rwanda.

METHODS

In total 87 blood samples were randomly collected from 37 males and 40 females of indigenous chickens and 10 females of SASSO chickens purposively selected from 5 sites located in two districts of the Eastern Province of Rwanda. Genotyping by Sequencing (GBS) using DArTseq technology was employed. This involved the complexity reduction method through digestion of genomic DNA and ligation of barcoded adapters followed by PCR amplification of adapter-ligated fragments.

RESULTS

From 45,677 DArTseq SNPs and 25,444 SilicoDArTs generated, only 8,715 and 6,817 respectively remained for further analysis after quality control. The average call rates observed, 0.99 and 0.98 for DArTseq SNPs and SilicoDArTs respectively were quite similar. The polymorphic information content (PIC) from SilicoDArTs (0.33) was higher than that from DArTseq SNPs (0.22). DArTseq SNPs and SilicoDArTs had 34.4% and 34% of the loci respectively mapped on chromosome 1. DArTseq SNPs revealed distance averages of 0.17 and 0.15 within IC and SASSO chickens respectively while the respective averages observed with SilicoDArTs were 0.42 and 0.36. The average genetic distance between IC and SASSO chickens was moderate for SilicoDArTs (0.120) compared to that of DArTseq SNPs (0.048). The PCoA and population structure clustered the chicken samples into two subpopulations (1 and 2); 1 is composed of IC and 2 by SASSO chickens. An admixture was observed in subpopulation 2 with 12 chickens from subpopulation 1.

CONCLUSIONS

The application of DArTseq markers have been proven to be effective and efficient for genetic relationship between IC and separated IC from exotic breed used which indicate their suitability in genomic studies. However, further studies using all chicken genetic resources available and large big sample sizes are required.

Application of next-generation sequencing for the high-resolution typing of MHC-B in Korean native chicken

The major histocompatibility complex-B (MHC-B) region of chicken is crucially important in their immunogenesis and highly diverse among different breeds, lines, and even populations. Because it determines the resistance/susceptibility to numerous infectious diseases, it is important to analyze this genomic region, particularly classical class I and II genes, to determine the variation and diversity that ultimately affect antigen presentation. This study investigated five lines of indigenous Korean native chicken (KNC) and the Ogye breed using next-generation sequencing (NGS) data with Geneious Prime-based assembly and variant calling with the Genome Analysis Toolkit (GATK) best practices pipeline. The consensus sequences of MHC-B (BG1-BF2) were obtained for each chicken line/breed and their variants were analyzed. All of the Korean native chicken lines possessed an excessive number of variants, including an ample amount of high-impact variants that provided useful information regarding modified major histocompatibility complex molecules. The study confirmed that next-generation sequencing techniques can effectively be used to detect MHC variabilities and the KNC lines are highly diverse for the MHC-B region, suggesting a substantial divergence from red junglefowl.

The diversity of smallholder chicken farming in the Southern Highlands of Tanzania reveals a range of underlying production constraints

The poultry industry in Tanzania has grown steadily over the past decade. We surveyed 121 chicken farming households along an intensification gradient from backyard to semi-intensive and intensive production systems based on rearing system and assumed purpose and poultry breed in the Iringa region. About 30% of households had more than one breed and/or rearing system combination. The subdivision of poultry systems was refined by adding the size of the flocks to highlight variation in scale of operations. On this basis we distinguished 3 main types: 1) subsistence small-scale free-range chicken production; 2) market-oriented small to medium scale semi-intensive and 3) small to medium-large scale intensive systems. ‘Intensification’ involves the transition from keeping indigenous chickens to improved dual-purpose and exotic breeds driven by greater productivity and potential for income generation. The more intensive the production system, the more the intensity and diversity of diseases identified by farmers as their main problem, which was partly attributed to the greater sensitivity of the improved breeds, poor veterinary measures, and the high chicken density facilitating disease spread. Based on the survey we constructed a problem tree to classify the underlying constraints and their interrelations, and to identify common root causes, based on which we propose practical solutions to improve chicken production. Development of medium-large scale systems is particularly constrained by a limited supply of 1-day-old chicks and theft. By contrast, intensification of small-scale systems is constrained by limited access to quality feed, vaccines and medicines, capital, and lack of a reliable market, partly due to the absence of farmer organization. These constraints can be addressed through formation of producer groups and promotion of outgrower and enterprise development models. Enterprise development appears to be the most promising business model for smallholder chicken farmers given that it allows farmers more freedom in decision-making and management while strengthening linkages with input suppliers and output markets to ensure a viable and profitable business.

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.

Comparison of major histocompatibility complex-B variability in Sri Lankan indigenous chickens with five global chicken populations using MHC-B SNP panel

In the present study, we investigated the major histocompatibility complex (MHC)-B haplotypes diversity of Sri Lankan indigenous chickens from three different geographical sites consisting of highly mixed populations using 90 SNPs in the MHC-B region. A total of 48 haplotypes were identified. Those included 37 novel haplotypes and 11 previously identified 'standard' haplotypes. The MHC-linked marker, LEI0258, had 23 alleles showing less diversity than defined by MHC-B SNP haplotypes. Among those identified haplotypes, five standard haplotypes-BSNP-O02, BSNP-M01, BSNP-A04, BSNP-K03, BSNP-T04-were most commonly observed, suggesting past introgression of imported breeds. Comparison of the MHC-B haplotypes of Sri Lankan and four other global populations with previously defined haplotypes indicated the sharing of 23 standard haplotypes with common origins. Novel haplotypes are population-specific and not shared among the geographical boundaries. Backyard indigenous chickens are unselected, highly crossbred, and generally thrive under dynamic environmental conditions. Hence free-range production systems may be responsible for maintaining high diversity in the MHC-B region with novel haplotypes.

Comparative Analysis of SLA-1, SLA-2, and DQB1 Genetic Diversity in Locally-Adapted Kenyan Pigs and Their Wild Relatives, Warthogs

Swine leukocyte antigen (SLA) plays a central role in controlling the immune response by discriminating self and foreign antigens and initiating an immune response. Studies on SLA polymorphism have demonstrated associations between SLA allelic variants, immune response, and disease resistance. The SLA polymorphism is due to host-pathogen co-evolution resulting in improved adaptation to diverse environments making SLA a crucial genomic region for comparative diversity studies. Although locally-adapted African pigs have small body sizes, they possess increased resilience under harsh environmental conditions and robust immune systems with reported tolerance to some diseases, including African swine fever. However, data on the SLA diversity in these pigs are not available. We characterized the SLA of unrelated locally-adapted domestic pigs from Homa Bay, Kenya, alongside exotic pigs and warthogs. We undertook SLA comparative diversity of the functionally expressed SLA class I (SLA-1, SLA-2) and II (DQB1) repertoires in these three suids using the reverse transcription polymerase chain reaction (RT-PCR) sequence-based typing (SBT) method. Our data revealed higher genetic diversity in the locally-adapted pigs and warthogs compared to the exotic pigs. The nucleotide substitution rates were higher in the peptide-binding regions of the SLA-1, SLA-2, and DQB1 loci, indicative of adaptive evolution. We obtained high allele frequencies in the three SLA loci, including some breed-specific private alleles, which could guide breeders to increase their frequency through selection if confirmed to be associated with enhanced resilience. Our study contributes to the growing body of knowledge on genetic diversity in free-ranging animal populations in their natural environment, availing the first DQB1 gene data from locally-adapted Kenyan pigs.

MHC-B variation in maternal and paternal synthetic lines of the Argentinian Campero INTA chicken

The Campero-INTA chicken of Argentina was developed to provide a robust bird that can survive under Argentinian pasture conditions with no significant additional nutrition, producing a source of animal protein for small producers or low-income families. In previous work, we described the AH paternal line of Campero and its Major Histocompatibility Complex B region (MHC-B) variation. In this work we analyzed the three remaining synthetic lines used to produce the Campero-INTA production bird: lines AS, A, and E. Because of the association between variation within the MHC of chickens and disease resistance, MHC variation within this breed is of particular interest. MHC variability within the lines used to produce the Campero-INTA chicken was examined using a 90 SNP panel encompassing the chicken MHC-B region plus the VNTR, LEI0258, located within the chicken MHC. Across all 4lines 12 haplotypes were found, with 7 of these being previously reported in North America/European breeds, reflecting the original breed sources for these birds. Three Campero unique haplotypes were found, 2 of which likely originated from MHC recombination events. MHC-B variation for all lines involved with production of the final Campero-INTA bird has now been determined.

High Levels of Genetic Variation in MHC-Linked Microsatellite Markers from Native Chicken Breeds

The major histocompatibility complex (MHC) is a highly polymorphic gene region that regulates cellular communication in all specific immune responses. In this study, we investigated 11 microsatellite (MS) markers in the MHC-B region of chicken populations from four countries: Sri Lanka, Bangladesh, South Korea, and Nigeria. The MS markers were divided into two sets. Set 1 included five novel MS markers, which we assessed using 192 samples from 21 populations. Set 2 included six previously reported markers, which we assessed using 881 samples from 29 populations. The Set 1 MS markers had lower polymorphism (polymorphic information content (PIC) < 0.5) than the Set 2 markers (PIC = 0.4-0.9). In all populations, the LEI0258 marker was the most polymorphic, with a total of 38 alleles (PIC = 0.912, expected heterozygosity (He) = 0.918). Local populations from Sri Lanka, Bangladesh, and Nigeria had higher allele diversity and more haplotypes for Set 2 MS markers than Korean and commercial populations. The Sri Lankan Karuwalagaswewa village population had the highest MHC diversity (mean allele number = 8.17, He = 0.657), whereas the white leghorn population had the lowest (mean allele number = 2.33, He = 0.342). A total of 409 haplotypes (89 shared and 320 unique), with a range of 4 (Rhode Island red) to 46 (Karuwalagaswewa village (TA)), were identified. Among the shared haplotypes, the B21-like haplotype was identified in 15 populations. The genetic relationship observed in a neighbour-joining tree based on the DA distance agreed with the breeding histories and geographic separations. The results indicated high MHC diversity in the local chicken populations. The difference in the allelic pattern among populations presumably reflects the effects of different genotypes, environments, geographic variation, and breeding policies in each country. The selection of MHC allele in domestic poultry can vary due to intensification of poultry production. Preserved MHC diversity in local chicken provides a great opportunity for future studies that address the relationships between MHC polymorphisms and differential immune responses.

Major histocompatibility complex B variability in Korean native chicken breeds

Adaptive genetic variations have direct influences on the fitness traits of the animal. The major histocompatibility complex B (MHC-B) region is responsible for adaptive and innate immune responses in chickens. In native Korean chicken breeds, no information on serologically defined B haplotypes is available. We investigated the MHC-B diversity in 5 restored lines of Korean native chicken and Ogye chicken breeds using a recently described MHC-B single-nucleotide polymorphism (SNP) panel and the MHC-linked LEI0258 variable number of tandem repeat marker. High SNP haplotype diversity was observed in Korean native chicken breeds with an average of 9.7 MHC-B SNP haplotypes per line. The total number of haplotypes ranged from 6 to 12 per line, and population-specific haplotypes ranged from 3 to 4. A total of 41 BSNP haplotypes, including 26 novel population-specific haplotypes and 15 common haplotypes, were reported over all populations. The 15 common haplotypes included 7 novel and 8 previously reported standard haplotypes. Selection and breeding evidence supports the observation of common haplotypes between the Korean native chicken and exotic breeds. Similarly, the LEI0258 marker showed allele variation, between 193 bp and 474 bp having 5 to 8 alleles per population. Some of these alleles (193, 249, 309, and 443 bp) were shared and more frequently observed. Comparison between SNP haplotypes and LEI0258 allele sizes for the same samples showed that some LEI0258 allele sizes correspond to more than one BSNP haplotype. The use of the MHC-B SNP panel greatly enhances the identification of MHC diversity compared with the sole use of the LEI0258 marker in native chicken populations.

Genetic diversity at major histocompatibility complex and its effect on production and immune traits in indigenous chicken breeds of India

The genetic diversity at major histocompatibility complex (MHC) in indigenous chicken breeds of India (Ghagus and Nicobari) in comparison with the White Leghorn (WLH) breed was investigated by genotyping the MHC-linked LEI0258 marker. Altogether 38 alleles and 96 genotypes were observed among three breeds. The observed and effective alleles were highest in Ghagus (23, 8.3) followed by Nicobari (14, 3.2) and WLH (10 and 2.2) breeds. The size of alleles ranged from 193 to 489 bp in Ghagus, 193 to 552 bp in Nicobari and 241 to 565 bp in the WLH breed. The number of private alleles was also highest in Ghagus (18) followed by Nicobari (8) and WLH (5) breeds. The most frequent allele was 261 bp in WLH (66 %), 343 bp in Nicobari (50.4 %) and 309 bp in the Ghagus (28.15 %) breed. Observed and expected heterozygosities were highest in Ghagus (0.83, 0.88) followed by Nicobari (0.58, 0.68) and WLH (0.53, 0.54). The genetic distance (Nei) between Ghagus and Nicobari breeds (2.24) was higher as compared to that of Ghagus and WLH (1.23) and that between Nicobari and WLH breeds (0.89). Association analysis revealed significant influence of MHC alleles on body weight, egg production in Ghagus and WLH breeds and antibody titres to Newcastle disease vaccine in the Nicobari breed.

Association of LEI0258 Marker Alleles and Susceptibility to Virulent Newcastle Disease Virus Infection in Kuroiler, Sasso, and Local Tanzanian Chicken Embryos

Newcastle disease (ND) control by vaccination and an institution of biosecurity measures is less feasible in backyard chicken in developing countries. Therefore, an alternative disease control strategy like the genetic selection of less susceptible chicken genotypes is a promising option. In the present study, genetic polymorphism of LEIO258 marker and association with susceptibility to virulent Newcastle disease virus (NDV) infection in Kuroilers, Sasso, and local Tanzanian chicken embryos were investigated. Samples from high (15%) and less (15%) susceptible cohorts were genotyped by sequencing of LEI0258 marker. A total of 75 DNA sequences comprised of 29 Kuroiler, 29 local Tanzanian chickens, and 17 Sasso were analyzed. Neighbor-joining phylogenetic trees were constructed to depict the clustering of LEI0258 marker alleles and relationship with susceptibility. Alleles with frequency ≥3 were considered for association with susceptibility by the use of the inference technique. The present findings suggest that some LEI0258 marker genetic polymorphisms apart from LEI0258 marker allelic based on sizes may be linked with chicken MHC-B haplotypes that confer chickens variability in resistance or susceptibility to infections. Furthermore, these results demonstrate the presence of relationship between LEI0258 marker polymorphisms and variations in chicken susceptibility to NDV infection, which could be utilized in breeding programs designed to improve chicken disease resistance.

Genetic diversity and population structure of indigenous chicken in Rwanda using microsatellite markers

Rwanda has about 4.5 million of indigenous chicken (IC) that are very low in productivity. To initiate any genetic improvement programme, IC needs to be accurately characterized. The key purpose of this study was to ascertain the genetic diversity of IC in Rwanda using microsatellite markers. Blood samples of IC sampled from 5 agro-ecological zones were collected from which DNA was extracted, amplified by PCR and genotyped using 28 microsatellite markers. A total of 325 (313 indigenous and 12 exotic) chickens were genotyped and revealed a total number of 305 alleles varying between 2 and 22 with a mean of 10.89 per locus. One hundred eighty-six (186) distinct alleles and 60 private alleles were also observed. The frequency of private alleles was highest in samples from the Eastern region, whereas those from the North West had the lowest. The influx of genes was lower in the Eastern agro-ecological zone than the North West. The mean observed heterozygosity was 0.6155, whereas the average expected heterozygosity was 0.688. The overall inbreeding coefficient among the population was 0.040. Divergence from the Hardy-Weinberg equilibrium was significant (p<0.05) in 90% of loci in all the populations. The analysis of molecular variance revealed that about 92% of the total variation originated from variation within populations. Additionally, the study demonstrated that IC in Rwanda could be clustered into four gene groups. In conclusion, there was considerable genetic diversity in IC in Rwanda, which represents a crucial genetic resource that can be conserved or optimized through genetic improvement.

Advances in methodologies for detecting MHC-B variability in chickens

The chicken major histocompatibility B complex (MHC-B) region is of great interest owing to its very strong association with resistance to many diseases. Variation in the MHC-B was initially identified by hemagglutination of red blood cells with specific alloantisera. New technologies, developed to identify variation in biological materials, have been applied to the chicken MHC. Protein variation encoded by the MHC genes was examined by immunoprecipitation and 2-dimensional gel electrophoresis. Increased availability of DNA probes, PCR, and sequencing resulted in the application of DNA-based methods for MHC detection. The chicken reference genome, completed in 2004, allowed further refinements in DNA methods that enabled more rapid examination of MHC variation and extended such analyses to include very diverse chicken populations. This review progresses from the inception of MHC-B identification to the present, describing multiple methods, plus their advantages and disadvantages.

Innate Immune Genes Associated With Newcastle Disease Virus Load in Chick Embryos From Inbred and Outbred Lines

Newcastle disease virus (NDV) causes substantial economic losses to smallholder farmers in low- and middle-income countries with high levels of morbidity and mortality in poultry flocks. Previous investigations have suggested differing levels of susceptibility to NDV between specific inbred lines and amongst breeds of chickens, however, the mechanisms contributing to this remain poorly understood. Studies have shown that some of these differences in levels of susceptibility to NDV infection may be accounted for by variability in the innate immune response amongst various breeds of poultry to NDV infection. Recent studies, in inbred Fayoumi and Leghorn lines, uncovered conserved, breed-dependent, and subline-dependent responses. To better understand the role of innate immune genes in engendering a protective immune response, we assessed the transcriptional responses to NDV of three highly outbred Tanzanian local chicken ecotypes, the Kuchi, the Morogoro Medium, and the Ching’wekwe. Hierarchical clustering and principal coordinate analysis of the gene expression profiles of 21-day old chick embryos infected with NDV clustered in an ecotype-dependent manner and was consistent with the relative viral loads for each of the three ecotypes. The Kuchi and Morogoro Medium exhibit significantly higher viral loads than the Ching’wekwe. The results show that the outbred ecotypes with increased levels of expression of CCL4, NOS2, and SOCS1 also had higher viral loads. The higher expression of SOCS1 is inconsistent with the expression in inbred lines. These differences may uncover new mechanisms or pathways in these populations that may have otherwise been overlooked when examining the response in highly inbred lines. Taken together, our findings provide insights on the specific conserved and differentially expressed innate immune-related genes involved the response of highly outbred chicken lines to NDV. This also suggests that several of the specific innate immunity related genes identified in the current investigation may serve as markers for the selection of chickens with reduced susceptibility to NDV.