- 1Department of Animal Science, University of Ghana, Accra, Ghana
- 2United States Agency for International Development Feed the Future Innovation Lab for Genomics to Improve Poultry, University of California, Davis, CA, United States
- 3School of Veterinary Medicine, University of Ghana, Accra, Ghana
- 4One Health Institute, University of California, Davis, CA, United States
- 5School of Veterinary Medicine, University of California, Davis, CA, United States
- 6Department of Animal Science, Iowa State University, Ames, IA, United States
- 7Department of Animal Science, University of California, Davis, CA, United States
Information on the genetic architecture of the production traits of indigenous African chicken is limited. We performed a genome-wide association study using imputed Affymetrix Axiom® 600K SNP-chip genotypes on 1,113 chickens from three agroecological zones of Ghana. After quality control, a total of 382,240 SNPs remained. Variance components and heritabilities for some growth, carcass and internal organ traits were estimated. The genetic and phenotypic correlations among these traits were also estimated. The estimated heritabilities of body weight at week 22 (BW22), average daily gain (ADG), dressed weight, breast weight, thigh weight, wing weight, drumstick weight, and neck weight were high and ranged from 0.50 to 0.69. Estimates of heritabilities for head weight, shank weight, and gizzard weight were moderate (0.31–0.35) while those of liver weight, back weight, dressing percentage, and heart weight were low (0.13–0.21). The estimated heritabilities of dressed weight, breast weight, wing weight, drumstick weight, neck weight, shank weight, and gizzard weight, corrected for BW22, were moderate (0.29–0.38), while the remaining traits had low heritability estimates (0.13–0.21). A total of 58 1-Mb SNP windows on chromosomes 1, 2, 4, 5, 6, 7, 8, 9, 13, 18, and 33 each explained more than 1% of the genetic variance for at least one of these traits. These genomic regions contained many genes previously reported to have effects on growth, carcass, and internal organ traits of chickens, including EMX2, CALCUL1, ACVR1B, CACNB1, RB1, MLNR, FOXO1, NCARPG, LCORL, LAP3, LDB2, KPNA3, and CAB39L. The moderate to high heritability estimates and high positive genetic correlations suggest that BW22, ADG, dressed weight, breast weight, thigh weight, wing weight, drumstick weight, and neck weight could be improved through selective breeding.
1 Introduction
Several indigenous African chicken ecotypes, including the Forest (FO), Interior savanna (IS) and Costal savanna (CS) ecotypes of Ghana have been reported (Walugembe et al., 2020). These chickens are hardy and thrive quite well in severe climates and environments (Pius et al., 2021). They are a major source of protein and play very important roles in sustaining the livelihoods of many households in Africa. Furthermore, there is also a perception that the meat of indigenous chicken ecotypes is very tasty, thus contributing in a large part, to a high demand for the meat of indigenous chickens in Ghana and many other parts of Africa (Asante-Addo and Weible, 2020; Ragasa et al., 2020).
The ability of indigenous chickens to thrive in different agroecological zones of Africa can in part be attributed to the variety of adaptive traits they possess, including thermotolerance, ability of escape predation, resistance to several endemic diseases (Mpenda et al., 2019), and a capacity to thrive under conditions of feed and water scarcity. Notwithstanding these important adaptive traits, indigenous African chicken ecotypes tend to have comparatively lower growth rates and body sizes (Munisi et al., 2015; Birteeb et al., 2016). As a result, many subsistence farmers tend to breed them with other breeds of chicken with the objective of increasing their body weights, a situation that can occasion the loss of their adaptive traits.
Some studies on the production traits of indigenous African chicken are available (Osei-Amponsah et al., 2013; Dekkers et al., 2018) but very few comprehensive genome-wide association studies (GWAS) on their production traits have been carried out. It is therefore imperative to unravel the genetic architecture of the production traits of indigenous chicken populations of Ghana to provide better insights for the genetic improvement of these traits in future. This GWAS therefore sought to examine the genetic architecture of the growth, carcass, and internal organ traits of the Forest, Interior and Coastal Savanna chicken ecotypes of Ghana.
2 Materials and methods
2.1 Experimental design
A total of 1,113 chickens, made up of the CS, the IS, and FO chicken ecotypes were used in this study. These are chickens whose parents have been described in Walugembe et al., 2020. Each ecotype was housed separately in deep litter pens. The dimensions of each pen were 2.54 m × 2.2 m × 2.2 m and housed a maximum of 40 birds. From day 1 to week 8, all birds were fed a standard chick starter mash, while from week 9 to week 22 they were fed a standard chick grower mash. Water was available on an ad libitum basis. Vaccination, feeding, and all other management practices were the same for all the chickens in the study.
At hatch, the body weight of every bird was measured and thereafter measured fortnightly until 22 weeks of age. From this data, average daily gain (ADG) was calculated as the linear regression of body weight on days of age. At week 23, the birds were euthanized and several carcass and internal organ traits including breast, thigh, wing, drumstick, neck, back, shank, head, gizzard, heart, liver and dressed weights were measured. Except for the gizzard, heart, and liver, the rest of the parts contained some skin.
2.2 Genotyping
Blood samples were collected from the wing veins of the chicks at 5 weeks of age using Whatman FTA cards (Sigma-Aldrich, St. Louis, MO, United States). Genomic DNA was isolated from the FTA cards for genotyping by sequencing (GBS) using a 5K GBS panel which was developed specifically for local Ghanaian and Tanzanian chicken ecotypes. A total of 5,238 SNPs were included in the SNP panel. Details on the development of the GBS panel are given in Walugembe et al. (2022). The genome sequences obtained were subjected to a customized SNP-pipeline that resulted in 5K SNP genotypes of each bird. These genotypes were then imputed to 382,240 SNPs that remained after quality control of high-density genotype data of relatives using Affymetrix Axiom® 600K SNP chip [the high-density genotype data are described in Walugembe et al. (2020)]. Imputation was performed using Fimpute (Sargolzaei et al., 2014).
2.3 Population structure
The FO, CS, and IS chicken ecotypes of Ghana that were used in this study are reported to originate from three ancestral populations (Walugembe et al., 2020). To deduce the proportion of ancestral subpopulations in each chicken, we carried out admixture analyses on the imputed genotypes using the Admixture software (Alexander et al., 2009), with the number of sub-populations set to three. These ancestral subpopulation proportions were used as covariates in the downstream genetic analyses.
2.4 Genetic parameters
Variance components and heritabilities were estimated using the following univariate linear model:
Body weight at week 22 (BW22) was also fixed as a covariate (Model 2) for some of the traits, i.e., breast, thigh, wing, drumstick, neck, back, shank, head, gizzard, heart, liver and dressed weight. The covariate explains out some of the variation in these traits due to body weight.
The genetic and phenotypic correlations between traits were estimated by fitting pairwise bivariate models with the same effects as in the univariate linear models. All models were implemented in ASReml 4 (Gilmour et al., 2015).
2.5 Genome-wide association and bioinformatics analyses
Genome-wide association analysis was performed using Bayes B (Meuwissen et al., 2001; Cheng et al., 2015), as implemented in the JWAS package (Cheng et al., 2015), to estimate the genetic variance accounted for by each 1-megabase (Mb) SNP window across Gallus gallus 6 genome build. Both Models 1 and 2 were used. 1-Mb SNP regions that explained more than 1% of the genetic variance in a trait were considered significant. To identify genes within significant 1-Mb SNP windows, we resorted to the Genome Data Viewer in NCBI—https://www.ncbi.nlm.nih.gov/genome/gdv/browser/genome/?id=GCF_000002315.6.
3 Results
3.1 Population structure
The admixture analysis based on identity by state as shown in Figure 1, indicates that notwithstanding the evidence of admixture, all the three ecotypes appeared to have come from three distinct ancestral populations. The IS ecotype had a high proportion of subpopulation 1 (0.83) but with lower proportions of subpopulations 2 (0.11) and 3 (0.05). The CS ecotype had a higher proportion of subpopulation 2 (0.74) and lower proportions of subpopulations 1 (0.09) and 3 (0.16), while the FO ecotype had a higher proportion of subpopulation 3 (0.63) and lower proportions of subpopulation 1 (0.22) and subpopulation 2 (0.15).
FIGURE 1. Admixture plot showing mixed ancestry among birds from the three Ghanaian chicken ecotypes.
3.2 Genetic parameters
The heritabilities and correlations of the growth, carcass, and internal organ traits of the FO, CS, and IS chicken ecotypes of Ghana, i.e., body weight at 22 weeks of age (BW22), average daily gain (ADG), breast weight (BrW), drumstick weight (DW), thigh weight (TW), wing weight (WW), dressed weight (DrW), and dressing percentage (DP) were estimated. These are presented in Tables 1, 2, while the estimated heritabilities and correlations of other body parts and internal organs, i.e., head weight (HW), neck weight (NeW), shank weight (ShW), back weight (BaW), gizzard weight (GzW), liver weight (LiW), and heart weight (HeW), are presented in Tables 3, 4.
TABLE 1. Estimates of heritabilities (along diagonal) and of genetic (above diagonal) and phenotypic (below diagonal) correlations for growth and carcass traits from Model 1 (without BW22 as a covariate).
TABLE 2. Estimates of heritabilities (along diagonal) and of genetic (above diagonal) and phenotypic (below diagonal) correlations for growth and carcass traits from Model 2 (with BW22 covariate).
TABLE 3. Estimates of heritabilities (along diagonal) and of genetic (above diagonal) and phenotypic (below diagonal) correlations for carcass traits and internal organs from Model 1 (without BW22 as a covariate).
TABLE 4. Estimates of heritabilities (along diagonal) and of genetic (above diagonal) and phenotypic (below diagonal) correlations for carcass traits and internal organs from Model 2 (with BW22 as a covariate).
For Model 1 (without BW22 as a covariate), as shown in Tables 1, 3, the estimated heritabilities of DP, LiW, HeW and BaW were low (0.13–0.21), while those for the other traits ranged from medium to high (0.31–0.69). Estimates of the genetic correlation between DP and BW22 and between DP and ADG were low. DP also had low phenotypic correlations with BW22, ADG, and TW. HW had low phenotypic correlations with GzW and LiW, while GzW also had a weak correlation with HeW. The rest of the traits had positive medium to high genetic and phenotypic correlations with each other.
For Model 2 (with BW22 as a covariate), as shown in Tables 2, 4, the estimated heritabilities for TW, LiW, and BaW were low (0.13–0.21), while the estimated heritabilities for some other traits were moderate and ranged from 0.29 for BrW to 0.38 for WW. GzW had negative estimates of genetic correlations with HW and NeW. LiW also had negative genetic correlation estimates with HW, NeW and ShW. The phenotypic correlation between BaW and ShW, and between GzW and LiW was negative. In addition, HW also had negative phenotypic correlations with GzW and LiW. Among the traits, ShW and NeW, BaW and NeW, GzW and LiW, and HeW, and NeW had high genetic correlation estimates, while HW and NeW, and HW, and BaW had high phenotypic correlations. The rest of the traits had low estimates of genetic and phenotypic correlations.
The effects of the ancestral subpopulation proportions as covariates in the downstream genetic analysis were statistically not significant.
3.3 Genome-wide association study
After quality control, a total of 1,113 birds and 382,240 SNPs were used for the GWAS. The percentage of genetic variance explained by 1-Mb genomic regions that are associated with the growth, carcass, and internal organ traits, with and without BW22 as covariate, are shown in Table 5 and the genes that are within significant 1-Mb windows are shown in Tables 6, 7.
TABLE 5. Percentage of genetic variance explained by 1-Mb genomic regions that are associated with growth, carcass, and internal organ traits (≥1.0% of genetic variance) based on the Bayes-B method, using Model 1 (without BW22 as covariate) and Model 2 (with BW22 as covariate).
TABLE 6. Positions and genes located in 1-Mb windows that explained ≥ 1% of genetic variance for growth, carcass, and internal organ traits (Model 1: without BW22 as a covariate).
TABLE 7. Positions and genes located in 1-Mb windows which explain ≥1% of genetic variance for growth, carcass, and internal organ traits (Model 2, with BW22 as a covariate).
Eight 1-Mb windows of SNPs that explained more than 1% of the genetic variance of body weight at week 22 were found. Two of these windows on chromosomes 1 and 4 explained more than 10% of the genetic variance of this trait (see Supplementary Figure S1A). These regions contained several annotated genes including SLIT2, LCORL, NCARPG, LAP3, MED28, INTS6, DLEU7, CKAP2, KPNA3, and CAB39L (Table 6).
Five 1-Mb windows of SNPs explained more than 1% of the genetic variance of ADG (Supplementary Figure S1B). One of these windows located on chromosome 4 explained more than 10% of the genetic variance of ADG. This window contained several genes including SLIT2, LCORL, FAM184B, NCARPG, QDPR, LAP3, MED28, and MIR218-1 (Table 6).
For breast weight, seven 1-Mb SNP windows on chromosomes 1, 4, and 33 reached the level of significance (see Supplementary Figure S2C). These windows together explained about 20.5% of the genetic variance of this trait. Many genes including RNASEH2B, INTS6, SERPINE3, DLEU7, WDFY2, NEK3, NEK5, CKAP2, TPTE 2, SLC25A15, FOXO1, SLIT2, LCORL, NCARPG, LAP3, SCN8A, FIGNL2, ANKRD33, ACVRL1, and ACVR1B were found in these regions (Table 6). With BW22 included as a covariate (Model 2), the number of significant SNP windows decreased from seven to three and explained only about 6.63% of the genetic variance of this trait. Within these SNP windows, several genes, including SCN8A, FIGNL2, ANKRD33, ACVRL1, ACVR1B, MYH10, ELAC2, DNAH9, TRNAM-CAU, PIRT, SHISA A6, RNF222, NDEL1, CCDC42, PIK3R5, PIK3R6 were co-located on chromosome 1, chromosome 18 and chromosome 33 (Table 7).
Two SNP windows on chromosomes 1 and 4 explained 32.5% of the genetic variance of drumstick weight (Supplementary Figure S2D). These windows contained 34 annotated genes including SLIT2, LCORL, NCARPG, LAP3, MED28, CLRN2, MIR218-1, KPNA3, CAB39L, CDADC1, RCBTB1, ARL11, SPRYD7, TRIM13, KCNRG, MIR15-A, SETDB2, LPAR6, and MED4 (Table 6). However, with BW22 as a covariate, the number of significant SNP windows increased from two to three, but they explained only 4.43% of the genetic variance of this trait. Seventy-two genes on chromosomes 1, 7, and 15 were observed. These include MBD5, ACVR2A, ORC4, EPC2, KIF5C, LYPD6, MMADHC, MIR1C, DERL3, SLC2A11, MYO7A and CAPN5 (Table 7).
One SNP window on chromosome 4 and three SNP windows on chromosome 1 explained about 40.9% of the genetic variance of WW (Supplementary Figure S3E). These regions contained many genes some of which are UVRAG, LRRC32, GUCY2F, ENSY, THAP12, TRNAP-AGG, TRNAP-UGG, WNT11, ART1, ART7B, ART7C, MADPRT1, IL18BP, RNF121, RNF169, TRPC2L, NUMA1, LAMTOR1, LRTOMT, ANAPC15, WDR73, ADAM15, SLCO2B1, TPBGL, PGM2L1, KCNE3, LIPT2, POLD3, CHRDL2, and XRRA1 (Table 6). With BW22 as a covariate, four SNP windows on chromosomes 1, 2, 4 and Z were significant, and together they explained 10.21% of the genetic variance of WW (Table 7).
Four SNP windows on chromosomes 1 and 4 explained 31.8% of the genetic variance of thigh weight (Supplementary Figure S3F). Many annotated genes, including LSAMP, EPHAS, CADM2, ROBO1, DSCAM, DMD, GPC5, PCDH9, NBEA, CNTNS, FAT3, DLG2, TENM4, IL1RAPL1, FAM155A, RNASEH2B, INTS6, FAM124A, SERPINE3, DLEU7, WDFY2, DHRS12, TMEM272, ATP7B, ALG11, NEK3, NEK5, CKAP2, VPS36, THSD1, FGL1L, TPTE2, SLC25A15, MRPS31, and FOXO1 were observed within these regions (Table 6). With the inclusion of BW22 as a covariate, only 1 SNP window on chromosome 1 was significant and explained about 1.01% of the genetic variances for thigh weight. This genomic region contained OLFM4 (Table 7).
One SNP window on chromosome 4 explained about 14% of the genetic variance of dressed weight (Supplementary Figure S4G). This region contained SLIT2, LCORL, FAM184B, NCARPG, QDPR, LAP3, MED28, MIR218 (Table 6).
Three SNP windows, two on chromosome 4 and one on chromosome 1, explained about 11% of the genetic variance of back weight (Supplementary Figure S4H). These genomic regions contained many genes as shown in Table 6. With BW22 as a covariate, only 2% of the genetic variance was explained by one SNP window on chromosome 4.
A total of three SNP windows on chromosomes 1, 4, and 8 explained about 5.5% of the genetic variance of head weight (Supplementary Figure S5I). The genes in these genomic regions can be seen in Table 6. The inclusion of BW22 as a covariate had two SNP windows on chromosomes 3 and 10 explaining 4.45% of the genetic variance of this trait and the genes in these genomic regions are shown in Table 7.
About 14.6% of the genetic variance of neck weight was explained by four SNP windows. Two of these SNPs were on chromosome 1 while the rest were on chromosomes 4 and 9 (Supplementary Figure S5J). The inclusion of BW22 as a covariate increased the amount of genetic variance explained by the same SNP windows to 19.45%. The annotated genes located in these genomic regions are shown in Tables 6, 7.
Two SNP windows on chromosome 1 and one SNP window on chromosome 4 explained about 44.5% of the genetic variance of shank weight (Supplementary Figure S6K). The inclusion of BW22 as a covariate increased the amount of genetic variance explained by the same SNP windows to 49.1%. Tables 6, 7 show the annotated genes that are found in these genomic regions.
Two SNP windows on chromosome 4 and one SNP window on chromosome 1 explained about 23.7% of the genetic variance of gizzard weight (Supplementary Figure S6L). The inclusion of BW22 as a covariate decreased the amount of genetic variance explained by the same SNP windows to 20.95%. Tables 6, 7 show the annotated genes that are found in these genomic regions.
One SNP window on chromosome 4 explained 20.5% of the genetic variance of liver weight (Supplementary Figure S7M). The inclusion of BW22 as a covariate increased the amount of genetic variance explained by this SNP windows to 24.43%. This window contained several genes (see Tables 6, 7).
About 12% of the genetic variance of heart weight was explained by seven SNP windows. Two of these were on chromosome 1 while the rest were on chromosomes 4, 6, 7, 13, and 33 (Supplementary Figure S7N). With the inclusion of BW22 as a covariate, the above-mentioned SNP windows, except for the SNP window on chromosome 13, explained 9.44% of the genetic variance of this trait. The annotated genes located in these genomic regions can also be seen in Tables 6, 7.
4 Discussion
4.1 Population structure
The admixture analysis based on identity by state (Figure 1), indicates that notwithstanding the evidence of admixture, all the three chicken ecotypes of Ghana appeared to have come from three distinct ancestral populations. Similar observations have been reported by Osei-Amponsah et al. (2010b). In another study involving the same chicken ecotypes, some of which were related to those used in this study, Walugembe et al. (2020) also arrived at the same conclusion. The IS ecotype appeared distinct from the CS and FO ecotypes, which, on the other hand, appeared to have a somewhat similar ancestry. The admixture of the CS and FO ecotypes could be a result of significant gene flows between the forest and coastal agroecological zones due to their proximity to each other.
4.2 Genetic parameters
Estimates of heritabilities for the growth and carcass traits of the three chicken ecotypes of Ghana, without BW22 as a covariate, as shown in Tables 1, 3 ranged from a low of 0.17 for BaW to a high of 0.69 for ADG. These estimates generally agree with the findings of several authors, including Rance et al. (2002), Venturini et al. (2014), El-Attroun et al. (2017), and El-Attroun et al. (2021). On the other hand, estimates of heritability for the internal organs ranged from a low of 0.13 for LiW to a medium of 0.35 for GzW. Heritability estimates for HeW and GzW are similar to those reported by Gaya et al. (2006), Venturini et al. (2014), and Dou et al. (2019) while the estimate of heritability for LiW was also similar to the findings of Moriera et al. (2019) but different from those of Venturini et al. (2014) and Dou et al. (2019). Growth traits of unselected chicken populations tend to have relatively high heritability. For example, Walugembe et al. (2020) also found a heritability for growth rate even after a challenge with La Sota Newcastle Disease Virus strain of above 0.4, and a pre-challenge heritability of 0.55.
Inclusion of BW22 as a covariate did not significantly affect heritability estimates of any of the internal organ traits (Table 4). However, with BW22 as a covariate, estimates of heritabilities of the carcass traits (Table 2) ranged from 0.21 for TW to 0.38 for WW. Direct selection could therefore be effective in improving some of these traits in Ghanaian local chicken populations.
High positive genetic and phenotypic correlations were found between most of the traits, except between DP and BW22 and between DP and ADG, which had low genetic and phenotypic correlation estimates. HW also had low phenotypic correlations with GzW and LiW. Estimates of the genetic correlation of BW22 with ADG, DrW, BrW, TW, WW, and DW were positive and high, suggesting that these traits could make indirect genetic gains when selection is directed at increasing BW22. Furthermore, the estimate of the genetic correlation between LiW and HeW was also high and similar to the findings of Rance et al. (2002) and Gaya et al. (2006). Inclusion of BW22 as covariate did not change the magnitude and direction of estimates of the phenotypic and genetic correlations for most carcass traits but led to a reduction in correlation estimates between the internal organ traits. Moderate to strong positive genetic and phenotypic correlations of body weight with carcass traits of chicken have also been reported by Venturini et al. (2014), Bungsrisawat et al. (2018), and El-Attrouny et al. (2021) but these conclusions are at variance with those of Osei-Amponsah (2010a), who reported weak to moderate negative phenotypic correlations between live weight and most carcass traits of Forest and Savannah chicken populations of Ghana. The small differences between estimates of the phenotypic versus the genetic correlations suggest that the environmental correlation was of similar magnitude as the genetic correlation. Furthermore, the medium to strong positive genetic correlation estimates between some of the growth, carcass, and internal organ traits of local chicken ecotype populations in Ghana suggests that selection based on body weight could enhance some of the carcass traits.
4.3 Positional candidate genes for growth traits of local chicken
Body weight is a polygenic trait, and chromosomes 1 and 4 of the chicken genome have been widely reported to harbour QTL for growth (Podisi, et al., 2013; Mebratie et al., 2019; Wang et al., 2022). In this study, a 1-Mb SNP window on chromosome 4 explained 19.7% and 11.6% of the genetic variance for BW22 and ADG, respectively, while another SNP window on chromosome 1 also explained 10.4% and 9.6% of the genetic variance of BW22 and ADG. These two chromosomal regions contain many genes, some of which have previously been reported to be associated with growth and carcass traits in chicken and other farm animals (Yang et al., 2021; Wang et al., 2022). The genes in these genomic regions include ligand dependent nuclear receptor corepressor like (LCORL) and non-SMC condensin I complex subunit G (NCAPG), which play an important role in arginine metabolism and are linked with growth in animals (Wu et al., 2009; Tetens et al., 2013; Tiensuu et al., 2019); leucine aminopeptidase 3 (LAP3) and LIM domain binding 2 (LDB2) genes which have an influence on growth traits of chicken (Gu et al., 2011). SNPs in karyopherin subunit alpha 3 (KPNA3) and RCBTB1 genes are also associated with growth in chicken (Wang et al., 2022; Zhu et al., 2023). Calcium binding protein 39 like (CAB39L) which is on chromosome 1 plays an important role in the regulation of food intake by activating AMP-activated protein kinase through the process of phosphorylation (Proszkowiec et al., 2006) and regulates body weight in chicken (Li et al., 2021; Zhang et al., 2021; Zhu et al., 2023). Some SNPs in the deleted lymphocytic leukemia 7 (DLEU7) gene have also been reported by Abdalhag et al. (2015) to be associated with growth traits in Jinghai yellow chickens. Forkhead box O1 (FOXO1) is another gene that has also been widely reported to influence average daily intake and the formation of adipose tissue and skeletal muscle of chickens (Xie et al., 2012). Xie et al. (2012) also observed that some SNPs in INTS6 are significantly associated with body weight of chicken at 90 days of age.
4.4 Positional candidate genes for carcass traits of local chicken
Carcass traits in chickens are also influenced by many genes with small individual effects. Genes associated with breast muscle weight, drumstick weight, thigh weight, wing weight, dressed weight, head weight, back weight, neck weight, and shank weight were mainly located on chromosomes 1 and 4. Some of these genes have been reported in the literature, including FOXO1, which plays an important role in muscle development by mediating PI3K-AKT-MAPK and PI3K-AKT-mTOR pathways (Xie et al., 2012; Jia et al., 2017). LCORL, a gene on chromosome 4, is reported to be expressed at higher levels in the breast muscle of high-muscle-weight chickens than in low-muscle-weight chickens (Liu et al., 2015). SLIT2 plays a regulatory role in the differentiation of osteoblast (Sun et al., 2009) and the inhibition of bone resorption (Park et al., 2019) and KPNA3 influences growth and muscle quality in chicken (Pértille et al., 2015; Li et al., 2022). LAP3 and FAM184B have been associated with organ weight in cattle and sheep (An et al., 2018; La et al., 2019). SERPINE3, one of the serine proteinase inhibitor (serpin) gene family members, and INTS6 are associated with bone quality (Guo et al., 2017) while Mediator Complex Subunit 4 (MED4) regulating vitamin D metabolism also affects development and maintenance of mineral ion homeostasis and skeletal integrity (Sutton and MacDonald, 2003). MLNR gene encodes a motilin receptor that promotes the release of growth hormone. In chicken, motilin receptor is largely involved in gastrointestinal functions including increments in Ca+2 levels and is associated with bone traits (Takahashi et al., 2014; Li et al., 2021). FNDC3A is also associated with bone traits (Li et al., 2021).
4.5 Positional candidate genes for internal organ traits of local chicken
The internal organs of animals are highly nutritious and contain high levels of bioavailable protein, amino acids, vitamins, and micronutrients (Fayemi et al., 2018). They are relatively cheaper than other meats and are easily available. As a result, their consumption among low-income and food insecure households in developing countries is quite high. Selective breeding aimed at improving the internal organ traits of chicken could therefore play a significant role in improving nutritional outcomes amongst children and low income-households in developing countries.
Chromosomes 1 and 4 contained several genes that exhibited a pleiotropic effect in gizzard, liver, and heart weights. These include SLIT2, LCORL, NCARPG, QDPR, LAP3, MED28, KPNA3, CAB39L, SPRYD7, TRIM13, KCNRG, SETDB2, MLNR, CRSLTR2, LPAR6, RB1, ITM2B, FNDC3A, and MED4. Other significant genes that are associated with internal organ traits of chicken include Empty spiracles homeobox 2 (EMX2) which is on chromosome 6 and is associated with heart weight. This gene plays a major role in transcriptional regulation of the slow myosin heavy chain 2 (MyHC2) gene in fast/slow embryonic muscle fibers (Weimer et al., 2013). Other genes that were found to be associated with heart weight include: CDK2 Associated Cullin Domain (CALCUL1) which is on chromosome 6 and is implicated in positive regulation of cell population proliferation and protein kinase activity (Kong et al., 2009; Zhang et al., 2021); Activin A receptor type 1B (ACVR1B) gene which encodes an activin A type IB receptor. Activins are dimeric growth and differentiation factors which belong to the transforming growth factor-beta (TGF-beta). SNPs in (TGF)-β2, 3, and 4 have been reported by Hosnedlova et al. (2020) to be associated with growth, skeletal and body composition traits of chicken; SLAIN motif family member 2 (SLAIN2) is involved in cytoplasmic microtubule organization (van der Vaart et al., 2011); Retinoblastoma 1 (RB1) is associated with body weight and bone traits in chicken (Zhang et al., 2011); Motilin receptor (MLNR) gene which encodes a motilin receptor and is also associated with growth and bone traits in chicken (Takahashi et al., 2014). Some significant positional genes for gizzard weight include the Follistatin like 5 (FSTL5) which is predicted to facilitate calcium ion binding activity and cell differentiation (Zhang et al., 2017); Calcium voltage-gated channel auxiliary subunit beta 1 (CACNB1) which affects skeletal muscle development in mice (Chen et al., 2011); SLAIN2; RB1; MLNR; SAMD9 and FNDC3A. A SNP window on chromosome 4 which explained about 21% of the genetic variance of liver weight contained SLIT2, LCORL, FAM184B, NCARPG, QDPR, LAP3, MED28, MIR218-1 genes.
5 Conclusion
We estimated genetic parameters and performed GWAS for several growth, carcass, and internal organ traits in local Ghanaian chicken ecotypes. The results show that heritabilities for growth and carcass traits were moderate to high, while the genetic correlations between these traits were generally positively high. The moderate to high heritabilities of BW22, ADG, dressed weight, drumstick weight, thigh weight, breast weight, wing weight, head weight, neck weight, shank weight, and gizzard weight indicates that these traits could be improved in these populations through selective breeding.
A total of 58 1-Mb SNP windows each of which explained more than 1% of the genetic variance of the growth, carcass, and internal organ traits studied contained many genes including EMX2, CALCUL1, ACVR1B, CACNB1, RB1, MLNR, FOXO1, NCARPG, LCORL, LAP3, LDB2, KPNA3, DLEU7 and CAB39L. These genes, which are reported to be associated with growth, carcass, and internal organ traits of chickens, could play important roles in future genetic improvement efforts targeted at the chicken ecotypes of Ghana.
Data availability statement
The data presented in the study are deposited in the European Variation Archive (EVA) at EMBL-EBI under accession number PRJEB69292 (https://www.ebi.ac.uk/eva/?eva-study=PRJEB69292).
Ethics statement
The animal study was approved by the Council for Scientific and Industrial Research Institutional Animal Care and Use Committee No. RPN 001/CSIR-IACUC/2014. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
RK: Formal Analysis, Investigation, Writing–original draft. EA-A: Data Curation, Formal Analysis, Software, Writing–review and editing. AN: Investigation, Writing–review and editing. HO: Investigation, Writing–review and editing. TK: Investigation, Validation, Writing–review and editing. RG: Investigation, Validation, Writing–review and editing. SL: Conceptualization, Methodology, Validation, Writing–review and editing. HZ: Conceptualization, Funding, Methodology, Resources, Supervision, Project Administration, Writing–review and editing. JD: Conceptualization, Methodology, Formal Analysis, Supervision, Writing–review and editing. BK: Conceptualization, Investigation, Methodology, Resources, Supervision, Project Administration, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was made possible by the generous support of the American people through the United States Agency for International Development (USAID) Feed the Future Innovation Lab for Genomics to Improve Poultry (cooperative agreement number AID-OAA-A-13-00080). The contents are the responsibility of the Feed the Future Innovation Lab for Genomics to Improve Poultry and do not necessarily reflect the views of USAID or the United States Government.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1297034/full#supplementary-material
References
Abdalhag, M. A., Zhang, T., Fan, Q. C., Zhang, X. Q., Zhang, G. X., Wang, J. Y., et al. (2015). Single nucleotide polymorphisms associated with growth traits in Jinghai yellow chickens. Genet. Mol. Res. 14 (4), 16169–16177. doi:10.4238/2015.December.8.6
Alexander, D. H., Novembre, J., and Lange, K. (2009). Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19 (9), 1655–1664. doi:10.1101/gr.094052.109
An, B., Xia, J., Chang, T., Wang, X., Miao, J., Xu, L., et al. (2018). Genome-wide association study identifies loci and candidate genes for internal organ weights in Simmental beef cattle. Physiol. genomics 50 (7), 523–531. doi:10.1152/physiolgenomics.00022.2018
Asante-Addo, C., and Weible, D. (2020). Imported versus domestic chicken consumption in Ghana: do attitudes and perceptions matter? J. Int. Food Agribus. Mark. 32 (5), 503–526. doi:10.1080/08974438.2020.1751767
Birteeb, P. T., Essuman, A. K., and Adzitey, F. (2016). Variations in morphometric traits of local chicken in Gomoa West district, southern Ghana. J. World’s Poult. Res. 6, 153–160.
Bungsrisawat, P., Tumwasorn, S., Loongyai, W., Nakthong, S., and Sopannarath, P. (2018). Genetic parameters of some carcass and meat quality traits in Betong chicken (KU line). Agric. Nat. Resour. 52 (3), 274–279. doi:10.1016/j.anres.2018.09.010
Chen, F., Liu, Y., Sugiura, Y., Allen, P. D., Gregg, R. G., and Lin, W. (2011). Neuromuscular synaptic patterning requires the function of skeletal muscle dihydropyridine receptors. Nat. Neurosci. 14 (5), 570–577. doi:10.1038/nn.2792
Cheng, H., Garrick, D., and Fernando, R. (2015). A fast and efficient Gibbs sampler for BayesB in whole-genome analyses. Genet. Sel. Evol. 47, 80. doi:10.1186/s12711-015-0157-x
Dekkers, J., Botchway, P. K., Amuzu-Aweh, E. N., Naazie, A., Aning, K. G., Zhou, H., et al. (2018). Genotypic and phenotypic characterisation of three local chicken ecotypes of Ghana based on principal component analysis and body measurements. Proc. World Congr. Genet. Appl. Livest. Prod. 11, 812.
Dou, T., Shen, M., Ma, M., Qu, L., Li, Y., Hu, Y., et al. (2019). Genetic architecture and candidate genes detected for chicken internal organ weight with a 600 K single nucleotide polymorphism array. Asian-Australasian J. Animal Sci. 32 (3), 341–349. doi:10.5713/ajas.18.0274
El-Attrouny, M. M., Iraqi, M. M., Khalil, M. H., and El-Moghazy, G. M. (2017). Genetic and phenotypic evaluation of growth traits in selection experiment performed in synthesized Benha chickens. Ann. Agric. Sci., Moshtohor 51 (1), 33–42.
El-Attrouny, M. M., Iraqi, M. M., and Mohamed, S. A. H. (2021). The estimation of genetic parameters for body weight, body dimension, and carcass traits in four Egyptian chickens strains. J. World’s Poult. Res. 11 (2), 230–240. doi:10.36380/jwpr.2021.28
Fayemi, P. O., Muchenje, V., Yetim, H., and Ahhmed, A. (2018). Targeting the pains of food insecurity and malnutrition among internally displaced persons with nutrient synergy and analgesics in organ meat. Food Res. Int. 104, 48–58. doi:10.1016/j.foodres.2016.11.038
Gaya, L. D. G., Ferraz, J. B. S., Rezende, F. M., Mourão, G. B., Mattos, E. C., Eler, J. P., et al. (2006). Heritability and genetic correlation estimates for performance and carcass and body composition traits in a male broiler line. Poult. Sci. 85 (5), 837–843. doi:10.1093/ps/85.5.837
Gilmour, A. R., Gogel, B. J., Cullis, B. R., Welham, S., and Thompson, R. (2015). ASReml user guide release 4.1 structural specification. Hemel Hempstead, United Kingdom: Hemel hempstead: VSN international ltd.
Gu, X., Feng, C., Ma, L., Song, C., Wang, Y., Da, Y., et al. (2011). Genome-wide association study of body weight in chicken F2 resource population. PloS one 6 (7), e21872. doi:10.1371/journal.pone.0021872
Guo, J., Sun, C., Qu, L., Shen, M., Dou, T., Ma, M., et al. (2017). Genetic architecture of bone quality variation in layer chickens revealed by a genome-wide association study. Sci. Rep. 7, 45317. doi:10.1038/srep45317
Hosnedlova, B., Vernerova, K., Kizek, R., Bozzi, R., Kadlec, J., Curn, V., et al. (2020). Associations between IGF1, IGFBP2 and TGFß3 genes polymorphisms and growth performance of broiler chicken lines. Animals 10 (5), 800. doi:10.3390/ani10050800
Jia, X., Ouyang, H., Abdalla, B. A., Xu, H., Nie, Q., and Zhang, X. (2017). miR-16 controls myoblast proliferation and apoptosis through directly suppressing Bcl2 and FOXO1 activities. Biochimica Biophysica Acta (BBA)-Gene Regul. Mech. 1860 (6), 674–684. doi:10.1016/j.bbagrm.2017.02.010
Kong, Y., Kejun, N., and Yin, Y. (2009). Identification and characterization of CAC1 as a novel CDK2-associated cullin. Cell Cycle 8 (21), 3552–3561. doi:10.4161/cc.8.21.9955
La, Y., Zhang, X., Li, F., Zhang, D., Li, C., Mo, F., et al. (2019). Molecular characterization and expression of SPP1, LAP3 and LCORL and their association with growth traits in sheep. Genes 10 (8), 616. doi:10.3390/genes10080616
Li, Y. D., Liu, X., Li, Z. W., Wang, W. J., Li, Y. M., Cao, Z. P., et al. (2021). A combination of genome-wide association study and selection signature analysis dissects the genetic architecture underlying bone traits in chickens. Animal 15 (8), 100322. doi:10.1016/j.animal.2021.100322
Li, Y. D., Xue, B. A. I., Xin, L. I. U., Wang, W. J., Li, Z. W., Ning, W. A. N. G., et al. (2022). Integration of genome-wide association study and selection signatures reveals genetic determinants for skeletal muscle production traits in an F2 chicken population. J. Integr. Agric. 21 (7), 2065–2075. doi:10.1016/s2095-3119(21)63805-4
Liu, R., Sun, Y., Zhao, G., Wang, F., Wu, D., Zheng, M., et al. (2013). Genome-wide association study identifies loci and candidate genes for body composition and meat quality traits in Beijing-You chickens. PLoS One 8 (4), e61172. doi:10.1371/journal.pone.0061172
Liu, R., Sun, Y., Zhao, G., Wang, H., Zheng, M., Li, P., et al. (2015). Identification of loci and genes for growth related traits from a genome-wide association study in a slow-× fast-growing broiler chicken cross. Genes and Genomics 37, 829–836. doi:10.1007/s13258-015-0314-1
Mebratie, W., Reyer, H., Wimmers, K., Bovenhuis, H., and Jensen, J. (2019). Genome wide association study of body weight and feed efficiency traits in a commercial broiler chicken population, a re-visitation. Sci. Rep. 9 (1), 922–1010. doi:10.1038/s41598-018-37216-z
Meuwissen, T. H., Hayes, B. J., and Goddard, M. (2001). Prediction of total genetic value using genome-wide dense marker maps. genetics 157 (4), 1819–1829. doi:10.1093/genetics/157.4.1819
Moreira, G. C. M., Salvian, M., Boschiero, C., Cesar, A. S. M., Reecy, J. M., Godoy, T. F., et al. (2019). Genome-wide association scan for QTL and their positional candidate genes associated with internal organ traits in chickens. BMC genomics 20 (1), 669–715. doi:10.1186/s12864-019-6040-3
Mpenda, F. N., Schilling, M. A., Campbell, Z., Mngumi, E. B., and Buza, J. (2019). The genetic diversity of local african chickens: a potential for selection of chickens resistant to viral infections. J. Appl. Poult. Res. 28 (1), 1–12. doi:10.3382/japr/pfy063
Munisi, W. G., Katule, A. M., and Mbaga, S. H. (2015). Comparative growth and livability performance of exotic, indigenous chickens and their crosses in Tanzania. Livest. Res. rural Dev. 27 (4).
Osei-Amponsah, R. (2010a). Phenotypic and genetic characterization of local chicken (Gallus gallus) ecotypes of Ghana. PhD Thesis. Legon: Department of Animal Science. University of Ghana.
Osei-Amponsah, R., Kayang, B. B., and Naazie, A. (2013). Phenotypic and genetic parameters for production traits of local chickens in Ghana. Anim. Genet. Resources/Resources Génétiques animales/Recursos Genéticos Anim. 53, 45–50. doi:10.1017/s2078633613000271
Osei-Amponsah, R., Kayang, B. B., Naazie, A., Osei, Y. D., Youssao, I. A., Yapi-Gnaorey, V. C., et al. (2010b). Genetic diversity of Forest and Savannah chicken populations of Ghana as estimated by microsatellite markers. Animal Sci. J. 81 (3), 297–303. doi:10.1111/j.1740-0929.2010.00749.x
Park, S. J., Lee, J. Y., Lee, S. H., Koh, J. M., and Kim, B. J. (2019). SLIT2 inhibits osteoclastogenesis and bone resorption by suppression of Cdc42 activity. Biochem. biophysical Res. Commun. 514 (3), 868–874. doi:10.1016/j.bbrc.2019.05.046
Pértille, F., Zanella, R., Felício, A. M., Ledur, M. C., Peixoto, J. O., and Coutinho, L. L. (2015). Identification of polymorphisms associated with production traits on chicken (Gallus gallus) chromosome 4. Genet. Mol. Res. 14 (3), 10717–10728. doi:10.4238/2015.September.9.11
Pius, L. O., Strausz, P., and Kusza, S. (2021). Overview of poultry management as a key factor for solving food and nutritional security with a special focus on chicken breeding in East African countries. Biology 10 (8), 810. doi:10.3390/biology10080810
Podisi, B. K., Knott, S. A., Burt, D. W., and Hocking, P. M. (2013). Comparative analysis of quantitative trait loci for body weight, growth rate and growth curve parameters from 3 to 72 weeks of age in female chickens of a broiler–layer cross. Bmc Genet. 14, 22–11. doi:10.1186/1471-2156-14-22
Proszkowiec-Weglarz, M., Richards, M. P., Ramachandran, R., and McMurtry, J. P. (2006). Characterization of the AMP-activated protein kinase pathway in chickens. Comp. Biochem. Physiology Part B Biochem. Mol. Biol. 143 (1), 92–106. doi:10.1016/j.cbpb.2005.10.009
Ragasa, C., Andam, K. S., Asante, S. B., and Amewu, S. (2020). Can local products compete against imports in West Africa? Supply-and demand-side perspectives on chicken, rice, and tilapia in Ghana. Glob. Food Secur. 26, 100448. doi:10.1016/j.gfs.2020.100448
Rance, K. A., McEntee, G. M., and McDevitt, R. M. (2002). Genetic and phenotypic relationships between and within support and demand tissues in a single line of broiler chicken. Br. Poult. Sci. 43 (4), 518–527. doi:10.1080/0007166022000004426
Sargolzaei, M., Chesnais, J. P., and Schenkel, F. S. (2014). A new approach for efficient genotype imputation using information from relatives. BMC genomics 15 (1), 478–512. doi:10.1186/1471-2164-15-478
Sun, H., Dai, K., Tang, T., and Zhang, X. (2009). Regulation of osteoblast differentiation by slit2 in osteoblastic cells. Cells Tissues Organs 190 (2), 69–80. doi:10.1159/000178020
Sutton, A. L., and MacDonald, P. N. (2003). Vitamin D: more than a “bone-a-fide” hormone. Mol. Endocrinol. 17 (5), 777–791. doi:10.1210/me.2002-0363
Takahashi, H., Rikimaru, K., Komatsu, M., Uemoto, Y., and Suzuki, K. (2014). Association between motilin receptor gene haplotypes and growth traits in Japanese Hinai-dori crossbred chickens. Asian-Australasian J. Animal Sci. 27 (3), 316–323. doi:10.5713/ajas.2013.13500
Tarsani, E., Kranis, A., Maniatis, G., Avendano, S., Hager-Theodorides, A. L., and Kominakis, A. (2019). Discovery and characterization of functional modules associated with body weight in broilers. Sci. Rep. 9 (1), 9125. doi:10.1038/s41598-019-45520-5
Tetens, J., Widmann, P., Kühn, C., and Thaller, G. (2013). A genome-wide association study indicates LCORL/NCAPG as a candidate locus for withers height in German Warmblood horses. Anim. Genet. 44 (4), 467–471. doi:10.1111/age.12031
Tiensuu, H., Haapalainen, A. M., Karjalainen, M. K., Pasanen, A., Huusko, J. M., Marttila, R., et al. (2019). Risk of spontaneous preterm birth and fetal growth associates with fetal SLIT2. PLoS Genet. 15 (6), e1008107. doi:10.1371/journal.pgen.1008107
van der Vaart, B., Manatschal, C., Grigoriev, I., Olieric, V., Gouveia, S. M., Bjelić, S., et al. (2011). SLAIN2 links microtubule plus end–tracking proteins and controls microtubule growth in interphase. J. Cell Biol. 193 (6), 1083–1099. doi:10.1083/jcb.201012179
Venturini, G. C., da Cruz, V. A. R., Rosa, J. O., Baldi, F., El Faro, L., Munari, D. P., et al. (2014). Genetic and phenotypic parameters of carcass and organ traits of broiler chickens.
Walugembe, M., Amuzu-Aweh, E. N., Botchway, P. K., Naazie, A., Aning, G., Wang, Y., et al. (2020). Genetic basis of response of Ghanaian local chickens to infection with a lentogenic Newcastle disease virus. Front. Genet. 11, 739. doi:10.3389/fgene.2020.00739
Walugembe, M., Amuzu-Aweh, E. N., Lim, K. S., Wang, Y., Chouicha, N., Kelly, T., et al. (2022). “Development of a low-density SNP panel for local Ghanaian and Tanzanian chicken ecotypes,” in Proceedings of 12th World Congress on Genetics Applied to Livestock Production (WCGALP) Technical and species orientated innovations in animal breeding, and contribution of genetics to solving societal challenges (Wageningen Academic Publishers), 1917–1920.
Wang, S., Wang, Y., Li, Y., Xiao, F., Guo, H., Gao, H., et al. (2022). Genome-wide association study and selective sweep analysis reveal the genetic architecture of body weights in a chicken F2 resource population. Front. Veterinary Sci. 9. doi:10.3389/fvets.2022.875454
Weimer, K., Theobald, J., Campbell, K. S., Esser, K. A., and DiMario, J. X. (2013). Genome-wide expression analysis and EMX2 gene expression in embryonic myoblasts committed to diverse skeletal muscle fiber type fates. Dev. Dyn. 242 (8), 1001–1020. doi:10.1002/dvdy.23988
Wu, G., Bazer, F. W., Davis, T. A., Kim, S. W., Li, P., Marc Rhoads, J., et al. (2009). Arginine metabolism and nutrition in growth, health, and disease. Amino acids 37, 153–168. doi:10.1007/s00726-008-0210-y
Xie, L., Luo, C., Zhang, C., Zhang, R., Tang, J., Nie, Q., et al. (2012). Genome-wide association study identified a narrow chromosome 1 region associated with chicken growth traits. PloS one 7 (2), e30910. doi:10.1371/journal.pone.0030910
Yang, R., Xu, Z., Wang, Q., Zhu, D., Bian, C., Ren, J., et al. (2021). Genome-wide association study and genomic prediction for growth traits in yellow-plumage chicken using genotyping-by-sequencing. Genet. Sel. Evol. 53, 82–15. doi:10.1186/s12711-021-00672-9
Zhang, D. Y., Sun, W. L., Ma, X., Zhang, P., Wu, W., Wu, H., et al. (2017). Up-regulated FSTL5 inhibits invasion of hepatocellular carcinoma through the Wnt/β-catenin/YAP pathway. Int. J. Clin. Exp. Pathology 10 (10), 10325–10333.
Zhang, H., Liu, S. H., Zhang, Q., Zhang, Y. D., Wang, S. Z., Wang, Q. G., et al. (2011). Fine-mapping of quantitative trait loci for body weight and bone traits and positional cloning of the RB1 gene in chicken. J. Animal Breed. Genet. 128 (5), 366–375. doi:10.1111/j.1439-0388.2011.00927.x
Zhang, Y., Wang, Y., Li, Y., Wu, J., Wang, X., Bian, C., et al. (2021). Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population. Heredity 126 (2), 293–307. doi:10.1038/s41437-020-00365-x
Keywords: GWAS, Ghanaian chicken ecotypes, carcass traits, internal organ traits, growth traits
Citation: Kanlisi RA, Amuzu-Aweh EN, Naazie A, Otsyina HR, Kelly TR, Gallardo RA, Lamont SJ, Zhou H, Dekkers J and Kayang BB (2024) Genetic architecture of body weight, carcass, and internal organs traits of Ghanaian local chickens. Front. Genet. 15:1297034. doi: 10.3389/fgene.2024.1297034
Received: 19 September 2023; Accepted: 26 February 2024;
Published: 13 March 2024.
Edited by:
Ntanganedzeni Mapholi, University of South Africa, South AfricaReviewed by:
Paolo Zambonelli, University of Bologna, ItalyIsidore Houaga, University of Edinburgh, United Kingdom
Copyright © 2024 Kanlisi, Amuzu-Aweh, Naazie, Otsyina, Kelly, Gallardo, Lamont, Zhou, Dekkers and Kayang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Boniface B. Kayang, bbkayang@ug.edu.gh