Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Anim. Sci., 21 January 2026

Sec. Animal Breeding and Genetics

Volume 6 - 2025 | https://doi.org/10.3389/fanim.2025.1732717

Maternal lineage diversity and genetic relationships of Sukuma chicken ecotype with other Tanzanian indigenous chickens based on mitochondrial DNA D-loop region

Zidihery Maquis Mhando,*&#x;Zidihery Maquis Mhando1,2*†Said Hemed Mbaga&#x;Said Hemed Mbaga1†Athumani Shabani Nguluma&#x;Athumani Shabani Nguluma1†Elisa Daniel Mwega&#x;Elisa Daniel Mwega3†Charles Moses Lyimo&#x;Charles Moses Lyimo1†
  • 1Department of Animal, Aquaculture and Range Sciences, College of Agriculture, Sokoine University of Agriculture (SUA), Morogoro, Tanzania
  • 2Department of Veterinary Microbiology, Parasitology and Biotechnology, Ministry of Livestock and Fisheries, Dodoma, Tanzania
  • 3Department of Veterinary Microbiology, Parasitology, and Biotechnology, College of Veterinary Medicine and Biomedical Sciences (CVMBS), Sokoine University of Agriculture, Morogoro, Tanzania

Understanding the genetic diversity, phylogenetic relationships, and demographic history of chickens is crucial for elucidating domestication processes, conserving indigenous genetic resources, and improving productivity and resilience in sustainable poultry production systems. This study investigates the genetic diversity, phylogenetic relationships, and demographic history of six Tanzanian indigenous chicken ecotypes, namely, Sukuma, Kuchi, Unguja, Pemba, Morogoro Medium, and Ching’wekwe, using mitochondrial DNA (mtDNA) D-loop sequences. A total of 100 mitochondrial DNA (mtDNA) D-loop sequences of 310 base pairs, representing six Tanzanian chicken ecotypes (Sukuma, Kuchi, Unguja, Pemba, Morogoro Medium, and Ching’wekwe), were analyzed. Genetic variation among the six Tanzanian chicken ecotypes was evaluated using several statistical approaches. Molecular diversity indices revealed substantial heterogeneity, with Sukuma chickens exhibiting the highest haplotype diversity (Hd = 0.95789) and nucleotide diversity (π = 0.30299), reflecting multiple maternal origins and deep evolutionary history. Conversely, Kuchi chickens showed the lowest diversity (Hd = 0.36316; π = 0.00158), consistent with founder effects. Analysis of molecular variance (AMOVA) indicated significant genetic differentiation (FST = 0.35452, p < 0.001), with 35.45% variation among populations and 64.55% within populations. Pairwise FST values and phylogenetic analyses revealed clear clustering, with Sukuma, Kuchi, and Morogoro Medium forming distinct clades, while coastal and Island ecotypes showed closer relationships. Neutrality tests suggested demographic stability in most ecotypes, except Kuchi, which showed signals of past expansion. These findings highlight Sukuma chickens as a critical genetic resource for breeding programs, while the low diversity of Kuchi chickens underscores the need for genetic management. The results provide a foundation for conservation strategies aimed at preserving the adaptive potential and genetic heritage of Tanzania’s indigenous chickens.

1 Introduction

Indigenous chickens are a critical genetic resource in sub-Saharan Africa, providing meat, eggs, income, and cultural value to rural households (Mtileni et al., 2011; FAO, 2015; Lyimo, 2025). In Tanzania, indigenous chicken ecotypes such as Sukuma, Kuchi, Pemba, Unguja, Morogoro Medium, and Ching’wekwe are reared under low-input, free-range systems and have adapted to diverse agro-ecological conditions, endemic diseases, and feed scarcity (Msoffe et al., 2004; Lyimo et al., 2013). These populations exhibit substantial phenotypic and genetic variation, making them vital for sustainable poultry improvement (Mtileni et al., 2011). Despite these potential, indigenous chicken populations are increasingly threatened by genetic erosion driven by uncontrolled crossbreeding, habitat changes, and the introduction of exotic breeds, which can undermine their unique adaptive traits. This necessitates a study on maternal lineage diversity and genetic relationships of the Sukuma chicken ecotype with other Tanzanian chickens. Critical insights based on the genetic diversity, phylogenetic relationships, and demographic history of chickens are vital for informative domestication processes, conserving indigenous genetic resources and improving productivity and resilience in sustainable poultry production systems. Historical and phylogeographic evidence indicates that domestic chickens reached Africa through at least two major introduction waves: one via Egypt during the Ptolemaic period (ca. 300 BC) and another through the Indian Ocean trade between the early and mid–1st millennium AD (MacDonald and Edward, 1993; Mwacharo et al., 2013; Lyimo, 2025). Archaeological remains from Egypt, Zanzibar, and the Horn of Africa, alongside linguistic and ethnographic records, point to complex dispersal patterns shaped by trade, migration, and cultural exchange. Genetic analyses support this by revealing a mixture of haplogroups from Southeast Asia, the Indian subcontinent, and possibly the Middle East (Williamson, 2000; Lyimo and Al-Qamashoui, 2022). Such multi-origin gene flow has contributed to the rich genetic landscape observed in Tanzanian ecotypes, including Sukuma chickens, which may have benefited from historical trade connections between inland Tanzania and coastal markets.

Sukuma chickens are also known as “Horashi” or “Horasi,” predominantly found in the Western and Lake Zone regions of Tanzania, which are characterized by a semiarid agro-ecological zone. They are highly valued for their distinctive features, including a larger body size, tall stature, long legs, elongated neck, abundant plumage, and docile temperament (Esatu et al., 2022). Oral histories from Sukuma livestock-keeping communities affirm that these chickens have long been integral to traditional farming systems. The name “Horashi,” derived from the Sukuma Bantu language, reflects this chicken’s physical attributes of larger body size, heavier weight, taller stature, and extended neck. Genetically, Sukuma chickens predominantly belong to haplogroup E, as reported by Liu et al. (2006), and they show close phenotypic resemblance to Asian Aseel chickens from the Indian subcontinent, which also cluster primarily within the same haplogroup (Kanakachari et al., 2023).

Mitochondrial DNA (mtDNA), particularly the hypervariable D-loop region, has proven to be a powerful marker for tracing maternal genetic diversity, population structure, and domestication history because of its rapid mutation rate and strict maternal inheritance (Saccone et al., 1987; Liu et al., 2006). Previous studies on African village chickens revealed multiple maternal origins, highlighting significant contributions from Asian lineages that reached the continent through Indian Ocean maritime exchange, Mediterranean trade, and inland dispersal routes (Mwacharo et al., 2013; Osman et al., 2016; Lawal and Hanotte, 2021; Lyimo, 2025). These diverse introductions produced distinct haplogroups that not only reflect the complex history of chicken domestication and dispersal but also strengthen the adaptive potential, resilience, and genetic richness of Tanzanian indigenous ecotypes.

Previous molecular studies on Tanzanian indigenous chickens using microsatellites and mtDNA D-loop sequences demonstrated clear genetic structuring among ecotypes, often clustering populations by geographical and cultural origin (Lyimo et al., 2013). For example, island populations (Unguja and Pemba) grouped together, mainland Bankiva-type ecotypes (Morogoro Medium and Ching’wekwe) formed another cluster, and game-type Kuchi remained genetically distinct patterns reflecting limited interbreeding and unique evolutionary histories. In this study, the inclusion of Sukuma chickens in such comparative frameworks allows for a deeper understanding of their maternal lineage diversity, potential adaptive traits, and their role in the broader genetic landscape of Tanzanian poultry. This knowledge is essential for designing conservation and breeding programs that will maintain productivity while safeguarding genetic heritage.

2 Materials and methods

2.1 Study area

The study involved investigating the maternal lineage diversity and genetic relationships of indigenous chicken ecotypes in Tanzania. Six indigenous chicken ecotypes, namely, Mwanza, Geita, Shinyanga, Tabora, Singida, Dodoma, and Morogoro, from five regions of mainland Tanzania were used in this study. Likewise, two indigenous chicken ecotypes from Unguja and Pemba Island were used (Figure 1).

Figure 1
Map of Tanzania showing regional boundaries with red dots marking study areas. Regions highlighted include Mwanza, Shinyanga, Singida, Morogoro, Pemba North, Pemba South, Unguja North and South, Mjini Magharibi, and Dar es Salaam. A legend explains symbols. North is indicated with an arrow.

Figure 1. Map showing the study areas.

2.2 Materials

2.2.1 mitochondrial DNA sequences

A total of 100 mtDNA sequences were used. Eighty sequences from Sukuma chickens were generated after the extraction of genomic mitochondrial DNA, and 20 sequences from five Tanzanian chicken ecotypes (Kuchi, Pemba, Unguja, Morogoro Medium, and Ching’wekwe), each being represented by four sequences, were retrieved from the NCBI GenBank database. Consequently, a total of 100 mtDNA D-loop sequences of 310 base pairs, representing six Tanzanian chicken ecotypes, were analyzed. Two of these, Pemba and Unguja, are named after their primary islands in the Zanzibar archipelago (Msoffe et al., 2004; Lyimo et al., 2013). The remaining four ecotypes originate from mainland Tanzania and include Sukuma from Mwanza; Shinyanga Kuchi from Mwanza, Shinyanga, Tabora, and Geita; Morogoro Medium from Morogoro; and Ching’wekwe chickens from Morogoro and Tanga. These ecotypes represent diverse geographic and cultural backgrounds, providing valuable insights into the genetic diversity, evolutionary history, and potential adaptive traits of indigenous chickens in the region.

2.2.2 Data collection

Fresh blood samples were collected from 100 Sukuma chicken ecotypes. Among these, 50 blood samples were collected from the Mwanza region, while the remaining 50 samples were collected from the Shinyanga region. Whole blood was obtained by bleeding the wing vein using a small sterile needle and then transferred immediately into vacutainer tubes containing ethylenediaminetetraacetic acid (EDTA), which was used as an anticoagulant. All samples were subsequently stored at 8°C to preserve DNA integrity for downstream molecular analysis.

2.2.3 DNA extraction

Genomic DNA was extracted from the blood using the Quick-DNA MiniPrep plus Kit (ZYMO Research, Co., USA, Cat. No. D4068) following the manufacturer’s protocol for nucleated blood samples. The extracted DNA was stored at −20°C until further analysis.

2.2.4 Amplification (polymerase chain reaction)

The primers L16750 (5′AGGACTACGGCTTGAAAAGC-3′) as forward primer and H547 (5′ATGTGCCTGACCGAGGAACCAG-3′) as reverse primer were used to amplify a fragment of 550 bp from the HV1 region of the Sukuma chicken mitochondrial genome. All polymerase chain reactions (PCRs) were carried out in a ProFlex™ 3 × 32-well PCR System (Thermo Fisher Scientific, Carlsbad, CA, USA) thermocycler. Each amplification reaction mix constituted a final volume of 25 μL, which included 12.5 μL of Quick-Load®, oneTaq2X Master Mix with standard buffer (New England Biolabs, UK), 2 μL of BSA (New England Biolabs, UK), 0.5 μL of forward and reverse primers, 5.5 μL of nuclease-free water, and 4 μL of DNA template. The reaction profile was as follows: initial denaturation at 94°C for 2 min, followed by 35 cycles of 94°C for 30 s, 58°C for 30 s, and 68°C for 1 min, and a final extension step at 68°C for 10 min.

2.2.5 Gel electrophoresis and visualization

Agarose gel electrophoresis was performed using the HU10 Mini-Plus (10 cm) horizontal gel electrophoresis unit (Scie-Plas, Camlab, UK) together with the nanoPAC-300P Electrophoresis Power Supply (SLS Flowgen, UK). The amplicons/PCR products, along with Quick-Load 100 bp DNA ladder (New England Biolabs, UK), were electrophoresed on a 1.5% agarose gel pre-stained with ethidium bromide DNA staining dye in 1× TAE buffer. The electrophoresis setup was set at 80 V for 45 min. Gel visualization was done using the T 2201 Ultraviolet (UV) transilluminator (Sigma Chemical Co., USA).

2.2.6 Polymerase chain reaction results

Eighty samples showed a successful amplification of the 550-bp fragment of the HV1 region of the chicken mitochondrial genome. DNA bands were clearly visualized under UV light, confirming successful amplification of target fragments (Figure 2).

Figure 2
Gel electrophoresis image showing thirteen lanes labeled M and 1 to 13. Lane M contains a DNA ladder with visible bands at 100 base pairs, 500 base pairs, 1000 base pairs, and 1500 base pairs. Lanes 1 to 13 show bands around 500 base pairs, indicating DNA fragments of similar size.

Figure 2. Gel image of representative amplicons produced by the primers used. M is the maker used (100 bp ladder); lanes 1 to 13 are positive samples.

2.2.7 mtDNA sequencing

The Sanger sequencing technique (the chain termination method) was employed for the determination of nucleotide sequences. A volume of 20 μL of each reaction was sent to the Macrogen Europe, Netherlands laboratory for Sanger dideoxynucleotide (ddNTP) sequencing, where each sample was sequenced in both directions, preceded by purification of the PCR products. PCR products from 80 samples, which are representatives from five districts of the study area, were sent for sequencing. Primers CR-forward 5′TCTATATTCCACATTTCTC3′ and CR-reverse 5′GCGAGCATAACCAAATGG-3′ were used for sequencing. Forward and reverse reads were manually edited and aligned to obtain consensus sequences using BioEdit v7.0.9.0 (Hall and Carlsbad, 2011).

2.2.8 Data analysis

Genetic variation among the six Tanzanian chicken ecotypes was evaluated using several statistical approaches. AMOVA and pairwise genetic differentiation (FST) were performed in Arlequin v3.5 to assess population structure and gene flow (Excoffier and Lischer, 2010). Pairwise FST genetic distance matrices were used to infer phylogenetic relationships, visualized in SplitsTree4 Version 4.14.2 (Huson and Bryant, 2006). A neighbor-joining phylogenetic tree was constructed in R using pairwise FST distances to visualize ecotype clustering patterns (Paradis and Schliep, 2019). Molecular diversity indices, including haplotype diversity (Hd), nucleotide diversity (Pi), number of segregating sites (S), and the average pairwise differences (k), were calculated using Arlequin and DNASP v6 (Rozas et al., 2017), with significance tested through 1,000 permutations. Demographic history and neutrality tests, including Tajima’s D, Fu’s Fs, Harpending’s raggedness index (r), and sum of squared deviation (SSD), were also computed to infer historical population dynamics and test for deviations from neutrality, providing insight into demographic stability or expansion trends.

3 Results

3.1 Population structure of Tanzanian indigenous chickens

AMOVA revealed significant genetic differentiation among the six ecotypes (FST = 0.35452, p < 0.001), partitioning 35.45% of variation among populations and 64.55% within populations (Table 1). The high and significant FST indicates strong population structure, while the larger within-population variance suggests considerable genetic diversity within each ecotype, reflecting rich genetic resources despite pronounced differentiation across groups.

Table 1
www.frontiersin.org

Table 1. AMOVA results showing genetic variation partitioning among and within populations.

3.2 Genetic distance metrics and gene flow estimates in Tanzania chicken populations

Pairwise FST values reveal varying levels of genetic differentiation among the six Tanzanian chicken ecotypes (Table 2), reflecting both genetic distance and historical gene flow. The greatest genetic distances occur between Morogoro Medium and Kuchi (0.42325), Kuchi and Sukuma (0.38158), and Morogoro Medium and Sukuma (0.37996), indicating limited gene flow. In contrast, minimal distances between Morogoro Medium and Unguja (0.01601) and between Pemba and Ching’wekwe (0.01101) suggest close genetic relationships and high recent gene flow. Overall, the results highlight pronounced genetic structuring among most ecotypes, with a few pairs maintaining substantial connectivity or sharing recent common ancestry.

Table 2
www.frontiersin.org

Table 2. Pairwise FST values indicating genetic differentiation among the six Tanzanian indigenous chicken ecotypes.

The phylogenetic tree constructed from pairwise FST distances (Figure 3) reveals distinct clustering patterns among the six Tanzanian chicken ecotypes, highlighting their genetic differentiation. The Sukuma, Kuchi, and Morogoro Medium ecotypes form clearly separated clades, indicating limited gene flow, suggesting long-term maternal isolation likely shaped by geographic barriers, cultural preferences, and localized breeding practices. In contrast, Pemba, Unguja, and Ching’wekwe cluster more closely, suggesting recent genetic exchange or shared ancestry. These relationships align with pairwise FST values, where high differentiation was observed among inland ecotypes, while lower values between island and coastal populations indicate recent gene flow. These patterns align with observed genetic distance metrics, highlighting both historical isolation and connectivity among Tanzania’s indigenous chicken populations. The distinct positioning of Sukuma chickens highlights their deep maternal lineage structure. The distinct positioning of Sukuma chickens highlights their deep maternal lineage structure, consistent with their high haplotype and nucleotide diversity.

Figure 3
Phylogenetic tree diagram showing relationships between six entities: Sukuma, Unguja, Morogoro, Pemba, Chingwekwe, and Kuchi. Sukuma is distinct, while Unguja and Morogoro are closely related. Pemba and Chingwekwe are also closely related, with both groups related to Kuchi. Scale bar at the top indicates distance.

Figure 3. Phylogenetic tree generated from FST pairwise difference, showing the clustering relationships among the six Tanzanian indigenous chicken populations.

3.3 Genetic diversity of Tanzanian indigenous chickens

The molecular diversity indices reveal marked genetic variation among the six Tanzanian chicken ecotypes (Table 3). Sukuma chickens exhibited the highest haplotype diversity (Hd = 0.95789) and nucleotide diversity (π = 0.30299), with 204 segregating sites (S) and an average pairwise difference (k) of 93.926, indicating rich maternal lineages and deep genetic structure. In contrast, Kuchi chickens showed the lowest diversity (Hd = 0.36316; π = 0.00158), with only 4 segregating sites and k = 0.489, suggesting past bottlenecks or founder effects. Unguja, Morogoro, Pemba, and Ching’wekwe displayed moderate variability (S = 10–11; k ≈ 3.9–4.6), likely reflecting admixture or localized selection. The pooled population diversity (Hd = 0.86765; π = 0.08756; S = 310; k = 27.144) highlights significant overall variation, essential for conservation and breeding programs aimed at maintaining adaptability and resilience.

Table 3
www.frontiersin.org

Table 3. Molecular diversity indices of the six Tanzanian chicken ecotypes.

3.4 Demographic history and neutrality analysis

The demographic history and neutrality test results in Table 4 suggest varied evolutionary dynamics among the six Tanzanian chicken ecotypes. Sukuma chickens show high Tajima’s D (2.63095, p = 0.993) and Fu’s Fs (5.20711, p = 0.981) values, both non-significant, indicating no strong deviation from neutrality and suggesting demographic stability or balancing selection. Unguja, Morogoro, and Pemba exhibit similarly high, non-significant positive values, reflecting potential population stability. Kuchi has negative Tajima’s D (−2.09760) and Fu’s Fs (−0.03325), also non-significant, possibly hinting at past expansion or purifying selection. Harpending’s raggedness and SSD p-values across ecotypes indicate generally good model fit for population stability rather than recent bottlenecks or expansions.

Table 4
www.frontiersin.org

Table 4. Demographic history and neutrality test parameters of six Tanzanian chicken ecotypes.

4 Discussion

The phylogenetic and molecular diversity analyses conducted in this study provide compelling evidence of substantial genetic structuring among six Tanzanian indigenous chicken ecotypes, namely, Sukuma, Kuchi, Unguja, Pemba, Morogoro Medium, and Ching’wekwe. This structuring reflects their diverse origins, evolutionary trajectories, and historical patterns of gene flow shaped by geographic isolation, human-mediated dispersal, and localized selection pressures (Lyimo et al., 2013; Lyimo, 2025). The observed FST value of 0.35452 is highly significant (p < 0.001), indicating strong genetic differentiation among the ecotypes. Such differentiation reinforces the role of maternal lineage diversity in preserving unique genetic resources essential for conservation, breeding programs, and sustainable improvement of Tanzania’s poultry sector (Mwacharo et al., 2013; Osman et al., 2016).

4.1 Genetic differentiation and population structure

The genetic differentiation observed among Tanzanian indigenous chicken ecotypes highlights both localized isolation and historical connections. AMOVA revealed that 35.45% of variation occurs among populations, with a high and significant fixation index (FST = 0.35452; p < 0.001), confirming strong structuring. Nevertheless, the substantial within-population variation (64.55%) highlights the richness of maternal lineages preserved within each ecotype, a pattern consistent with previous findings in African village chickens (MuChadeyi et al., 2008; Mtileni et al., 2011; Lyimo et al., 2013).

Pairwise FST comparisons indicate the greatest genetic differentiation between Morogoro Medium and Kuchi (0.42325), as well as between Kuchi and Sukuma (0.38158), suggesting long-term maternal isolation possibly reinforced by cultural preferences and geographic separation. Conversely, the lowest values, observed between Morogoro Medium and Unguja (0.01601) and between Pemba and Ching’wekwe (0.01101), reflect recent gene flow or shared ancestry likely mediated by interisland poultry exchange and coastal trade routes (Msoffe et al., 2004; Osman et al., 2016). These results align with the broader understanding that African village chickens exhibit a mosaic of differentiation shaped by both ancient introductions and contemporary management practices (Groeneveld et al., 2010). Such population structure highlights the need for conservation and breeding programs that safeguard both inter-ecotype distinctiveness and intrapopulation diversity critical for adaptation and resilience.

4.2 Phylogenetic relationships and maternal origins

The phylogenetic reconstruction based on FST values reveals distinct clades corresponding to Sukuma, Kuchi, and Morogoro Medium chickens, supporting their maternal divergence and relative isolation. These ecotypes likely retained unique lineages shaped by localized breeding systems and adaptation to inland environments. In contrast, the coastal and island ecotypes, Unguja and Pemba, cluster more closely, reflecting shared haplotypes and higher levels of maternal connectivity. Such relationships suggest stronger historical and contemporary exchanges facilitated by maritime and interisland trade networks (Mwacharo et al., 2013; Osman et al., 2016). Such clustering patterns align with the observed FST values, which indicate strong differentiation among most ecotypes yet closer relationships within the coastal and island groups. The Sukuma ecotype is particularly notable for its deep maternal structure and exceptionally high haplotype diversity, consistent with multiple maternal origins introduced through inland trade routes connecting western regions of Tanzania to coastal markets (Lyimo, 2025).

The observed phylogenetic topology reflects not only ecological and cultural selection but also the broader history of domestic chicken introduction and dispersal across Africa. Archaeological and historical evidence suggests that chickens entered the continent through multiple waves, including introductions via the Mediterranean and the Indian Ocean, with significant contributions from Southeast Asia and the Indian subcontinent (MacDonald and Edward, 1993; Mwacharo et al., 2013; Lyimo et al., 2013). These repeated introductions created a mosaic of maternal lineages, explaining the admixture found among coastal populations and the persistence of unique haplotypes within inland ecotypes. In particular, island and coastal groups show strong genetic connectivity, shaped by maritime trade and interisland poultry exchange, whereas inland populations such as the Sukuma ecotype retained distinct maternal signatures due to limited exposure to such external gene flow. These findings emphasize the importance of integrating phylogenetic evidence with archaeological and ethnographic records to trace the maternal history of African poultry and inform conservation strategies.

4.3 Molecular diversity patterns

Molecular diversity indices demonstrate substantial heterogeneity among the six Tanzanian indigenous chicken ecotypes. Sukuma chickens display the highest haplotype diversity (Hd = 0.95789), nucleotide diversity (π = 0.30299), segregating sites (S = 249), and average pairwise differences (k = 114.329), indicating multiple maternal lineages, long evolutionary history, and minimal bottleneck effects. Such high variability likely reflects historical trade connections and livestock exchanges, particularly in the western part of Tanzania and the Lake Zone, that introduced diverse maternal haplotypes (Lyimo et al., 2013).

The high haplotype diversity observed in Sukuma chickens may partly reflect breeding autonomy at the household level, supported by strong informal networks that enable the exchange of breeding stock. Such exchanges, though beneficial for genetic variation, can also risk genetic dilution if not complemented by formal conservation programs. In Tanzania, conservation and sustainable use of indigenous chicken ecotypes are influenced not only by genetic and ecological factors but also by the interplay between formal and informal institutions (Kapella et al., 2020). Formal institutions provide technical guidance, breeding policies, and research support, whereas informal structures such as traditional knowledge systems, local norms, and community networks play a critical role in maintaining breeding practices and genetic integrity at the grassroots level.

Kuchi chickens exhibited the lowest genetic diversity (Hd = 0.36316; π = 0.00158; S = 12; k = 0.873), reflecting founder effects, genetic drift, or strong directional selection. Their ancestry links to South Asian game fowl, traditionally bred for combat under rigorous selection (Lyimo et al., 2013). Mitochondrial analyses further showed that most Kuchi individuals clustered within a single haplogroup previously identified in Shamo game birds, underlining their narrow maternal base. This restricted variability parallels patterns observed in recently introduced or intensively selected populations, highlighting potential vulnerability to inbreeding and reduced adaptability (Liu et al., 2006).

Unguja, Pemba, Morogoro Medium, and Ching’wekwe populations exhibit intermediate diversity (Hd = 0.76842–0.90526; π = 0.00378–0.16249), reflecting a balance between localized adaptation and sporadic introgression. Such patterns indicate moderate gene flow among coastal and island ecotypes, likely mediated by interisland poultry trade and cultural exchange (Osman et al., 2016). This exchange has promoted genetic connectivity while retaining ecotype-specific variation. For these populations, sustaining current levels of gene flow while protecting distinct haplotypes is crucial for maintaining adaptive potential and ensuring long-term resilience under changing production environments (Grant and Bowen, 1998; Tajima, 1989; Fu, 1997).

4.4 Demographic history and neutrality tests

Neutrality test results provide valuable insights into the demographic history of Tanzanian indigenous chicken ecotypes. Most ecotypes, including Sukuma, Unguja, Morogoro Medium, and Pemba, exhibited high, non-significant positive Tajima’s D and Fu’s Fs values, suggesting demographic stability or balancing selection rather than recent population expansion. In particular, Sukuma’s high Tajima’s D (2.63095) and Fu’s Fs (5.20711) reflect the long-term maintenance of diverse haplotypes, possibly reinforced by cultural preferences for phenotypic variation and the retention of multiple maternal lineages (Mwacharo et al., 2013; Mtileni et al., 2011; Lyimo et al., 2013). Such stability may also be linked to their geographic isolation and adaptive significance in local environments.

In contrast, Kuchi chickens showed a negative Tajima’s D (−2.09760) and Fu’s Fs (−0.03325), though non-significant, which may suggest a historical expansion following a founder event, a scenario consistent with their low diversity and known derivation from Southeast Asian game fowl lineages (Lyimo et al., 2013; Liu et al., 2006). The narrow genetic base of the Kuchi ecotype signals vulnerability to disease and environmental change, necessitating strategies to broaden their genetic pool. Complementing these demographic signals, neutrality and diversity analyses offer deeper insights into the evolutionary forces shaping Tanzanian chicken ecotypes. Harpending’s raggedness index and SSD p-values across ecotypes further support a stable demographic model, with no strong evidence of recent bottlenecks. These patterns mirror observations in other African village chickens, where localized selection, limited long-distance dispersal, and occasional admixture shape long-term population structure (MuChadeyi et al., 2008).

4.5 Implications for conservation and breeding

The high genetic diversity in Sukuma (Hd = 0.95789; π = 0.30299) and Ching’wekwe chickens highlights their significance as reservoirs of adaptive traits valuable for breeding programs. This diversity enhances resilience to emerging diseases, climatic variability, and production challenges, aligning with the FAO (2015) guidelines on safeguarding livestock genetic resources. In Sukuma chickens, the deep maternal lineage structure suggests long-term adaptation to the Lake Zone’s agro-ecological conditions, making them ideal candidates for conservation-oriented breeding to maintain traits for resistance and environmental adaptability (Lyimo et al., 2013).

Conversely, the Kuchi ecotype possessing low diversity (Hd = 0.36316; π = 0.00158) reflects a narrow genetic base, increasing vulnerability to inbreeding depression and reduced adaptability. Targeted management, including genetic monitoring and controlled introgression from genetically compatible indigenous populations, could help mitigate these risks (Liu et al., 2006; Mtileni et al., 2011). Conservation strategies should combine in situ approaches by safeguarding traditional free-range management systems, with ex situ measures such as cryopreservation of germplasm. Breeding programs must integrate productivity goals with preservation of unique genetic signatures to avoid erosion of adaptive variation, ensuring both food security and cultural heritage preservation (FAO, 2015; Mwacharo et al., 2013). Selective breeding programs should aim to balance productivity traits with the conservation of unique genetic lineages to avoid erosion of adaptive diversity.

5 Conclusion

This study demonstrates marked genetic structuring among Tanzanian chicken ecotypes, shaped by historical trade routes, cultural selection, and geographical isolation. While certain populations exhibit signs of shared ancestry and recent gene flow, others, such as the Sukuma and Kuchi ecotypes, retain distinct genetic signatures that reflect their unique evolutionary paths. The exceptional diversity of Sukuma chickens positions them as a crucial genetic resource for breeding programs that may focus on resilience and adaptability. Contrarywise, Kuchi’s narrow genetic base emphasizes the need for strategies to broaden genetic variation and mitigate inbreeding risks. For ecotypes with intermediate diversity, maintaining existing gene flow while conserving unique haplotypes is essential to safeguarding the adaptive potential and long-term productivity of Tanzania’s indigenous chicken genetic heritage.

Data availability statement

The data presented in the study are deposited in the NCBI GenBank repository under accession numbers KP067445 – KP067545.

Ethics statement

The animal study was approved by Sokoine University of Agriculture. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

ZM: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. SM: Conceptualization, Data curation, Project administration, Supervision, Validation, Visualization, Writing – review & editing. AN: Data curation, Methodology, Software, Supervision, Validation, Visualization, Writing – review & editing. EM: Methodology, Project administration, Supervision, Visualization, Writing – review & editing. CL: Conceptualization, Data curation, Methodology, Project administration, Software, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. The authors sincerely thank the Ministry of Livestock and Fisheries for funding this research.

Acknowledgments

ZM thanks the Ministry of Livestock and Fisheries for permitting him to pursue MSc studies at SUA. The authors are also grateful to the Genome Science Center Laboratory at the College of Veterinary Medicine and Biomedical Sciences from SUA for providing facilities, guidance, and support during molecular analyses, including DNA extraction and sequencing of the chicken samples.

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

References

Esatu W., Goromela E. H., Kafuku S. H., Mpemba C. H., Chando M., Andrew O., et al. (2022). A guide to setting up a selective indigenous chicken improvement program: The Horasi breed in Tanzania. ILRI Manual No. 58. Nairobi, Kenya: International Livestock Research Institute (ILRI).

Google Scholar

Excoffier L. and Lischer H. E. L. (2010). Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x

PubMed Abstract | Crossref Full Text | Google Scholar

FAO (2015). The Second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture. Eds. Scherf B. D. and Pilling D. (Rome: Food and Agriculture Organization of the United Nations).

Google Scholar

Fu Y. X. (1997). Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147, 915–925. doi: 10.1093/genetics/147.2.915

PubMed Abstract | Crossref Full Text | Google Scholar

Grant W. S. and Bowen B. W. (1998). Shallow population histories in deep evolutionary lineages of marine fishes: Insights from sardines and anchovies and lessons for conservation. J. Heredity 89, 415–426. doi: 10.1093/jhered/89.5.415

Crossref Full Text | Google Scholar

Groeneveld L. F., Lenstra J. A., Eding H., Toro M. A., Scherf B., Pilling D., et al. (2010). Genetic diversity in farm animals – A review. Anim. Genet. 41, 6–31. doi: 10.1111/j.1365-2052.2010.02038.x

PubMed Abstract | Crossref Full Text | Google Scholar

Hall T. I. and Carlsbad C. (2011). BioEdit: an important software for molecular biology. GERF Bull. Biosciences. 2, 60–61.

Google Scholar

Huson D. H. and Bryant D. (2006). Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23, 254–267. doi: 10.1093/molbev/msj030

PubMed Abstract | Crossref Full Text | Google Scholar

Kanakachari M., Chatterjee R. N., Reddy M. R., Dange M., and Bhattacharya T. K. (2023). Indian Red Jungle fowl reveals a genetic relationship with South East Asian Red Jungle fowl and Indian native chicken breeds as evidenced through whole mitochondrial genome sequences. Front. Genet. 14, 1083976. doi: 10.3389/fgene.2023.1083976

PubMed Abstract | Crossref Full Text | Google Scholar

Kapella L. E., Nyanda S. S., and Mahonge C. P. (2020). The relevance of formal and informal institutions in local chicken genetic resource conservation: A case of Igunga District, Tanzania. Tanzania J. Agric. Sci. 21, 134–149. https://www.ajol.info/index.php/tjags/article/view/234424 (Accessed August 06, 2025).

Google Scholar

Lawal R. A. and Hanotte O. (2021). Domestic chicken diversity: Origin, distribution, and adaptation. Anim. Genet. 52, 385–394. doi: 10.1111/age.13091

PubMed Abstract | Crossref Full Text | Google Scholar

Liu Y. P., Wu G. S., Yao Y. G., Miao Y. W., Luikart G., Baig M., et al. (2006). Multiple maternal origins of chickens: out of the Asian jungles. Mol. Phylogenet. Evol. 38, 12–19. doi: 10.1016/j.ympev.2005.09.014

PubMed Abstract | Crossref Full Text | Google Scholar

Lyimo C. M. (2025). The origins and spread of domestic chickens in Africa: A synthesis of archaeological, ethnographic, and genetic perspectives. Int. J. Anim. Sci. Technol. 9, 123–139. doi: 10.11648/j.ijast.20250903.11

Crossref Full Text | Google Scholar

Lyimo C. M. and Al-Qamashoui B. (2022). Chicken maternal lineage retained a long historical relationship between Zanzibar and Oman. Tanzania J. Agric. Sci. 21, 277–287. Available online at: https://www.ajol.info/index.php/tjags/article/view/234448 (Accessed March 15, 2025).

Google Scholar

Lyimo C. M., Weigend A., Janßen-Tapken U., Msoffe P. L., Simianer H., and Weigend S. (2013). Assessing the genetic diversity of five Tanzanian chicken ecotypes using molecular tools. South Afr. J. Anim. Sci. 43, 499–510. doi: 10.4314/sajas.v43i4.7

Crossref Full Text | Google Scholar

MacDonald K. C. and Edward D. N. (1993). Chicken in Africa: the importance of Qasr Ibrim. Antiquity 67, 584–590. doi: 10.1017/S0003598X00045786

Crossref Full Text | Google Scholar

Msoffe P. L. M., Mtambo M. M. A., Minga U. M., Olsen J. E., Juul-Madsen H. R., Gwakisa P. S., et al. (2004). Productivity and reproductive performance of the free-range local domestic fowl ecotypes in Tanzania. Livest. Res. Rural Dev. 16. Available online at: http://www.lrrd.org/lrrd16/9/msof16067.htm (Accessed April 11, 2025).

Google Scholar

Mtileni B. J., MuChadeyi F. C., Maiwashe A., Chimonyo M., Groeneveld E., Weigend S., et al. (2011). Diversity and origins of South African chickens. Poult. Sci. 90, 2189–2194. doi: 10.3382/ps.2011-01505

PubMed Abstract | Crossref Full Text | Google Scholar

MuChadeyi F. C., Eding. H., Simianer H., Wolliny C. B. A., Groeneveld E., and Weigend S. (2008). Mitochondrial DNA D-loop sequences suggest a Southeast Asian and Indian origin of Zimbabwean village chickens. Anim. Genet. 39, 615–622. doi: 10.1111/j.1365-2052.2008.01785.x

PubMed Abstract | Crossref Full Text | Google Scholar

Mwacharo J. M., Bjørnstad G., Han J. L., and Hanotte O. (2013). The history of African village chickens: an archeological and molecular perspective. Afr. Archaeological Rev. 30, 97–114. doi: 10.1007/s10437-013-9128-1

PubMed Abstract | Crossref Full Text | Google Scholar

NCBI. (2025). National Center for Biotechnology Information: Retrieval and BLAST analysis of chicken mitochondrial DNA nucleotide sequences. GenBank. Available online at: https://www.ncbi.nlm.nih.gov/genbank/ (Accessed July 7–August 22, 2025).

Google Scholar

Osman S. A. M., Yonezawa T., and Nishibori M. (2016). Origin and genetic diversity of Egyptian native chickens based on complete sequence of mitochondrial DNA D-loop region. Poultry Sci. 95, 1248–1256. doi: 10.3382/ps/pew029

PubMed Abstract | Crossref Full Text | Google Scholar

Paradis E. and Schliep K. (2019). ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528. doi: 10.1093/bioinformatics/bty633

PubMed Abstract | Crossref Full Text | Google Scholar

Rozas J., Ferrer-Mata A., Sánchez-DelBarrio J. C., Guirao-Rico S., Librado P., Ramos-Onsins S. E., et al. (2017). DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299–3302. doi: 10.1093/molbev/msx248

PubMed Abstract | Crossref Full Text | Google Scholar

Saccone C., Attimonelli M., and Sbisa E. (1987). Structural elements highly preserved during the evolution of the D-loop-containing region in vertebrate mitochondrial DNA. J. Mol. Evol. 26, 205–211. doi: 10.1007/BF02099853

PubMed Abstract | Crossref Full Text | Google Scholar

Tajima F. (1989). Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123, 585–595. doi: 10.1093/genetics/123.3.585

PubMed Abstract | Crossref Full Text | Google Scholar

Williamson K. (2000). “Did chickens go west?,” in The origins and development of African livestock. Eds. Blench R. M. and MacDonald K. C. (UCL Press, London), 368–448. Chapter 23.

Google Scholar

Keywords: genetic diversity, mtDNA D-loop, phylogenetics, population dynamics, Tanzanian indigenous chickens

Citation: Mhando ZM, Mbaga SH, Nguluma AS, Mwega ED and Lyimo CM (2026) Maternal lineage diversity and genetic relationships of Sukuma chicken ecotype with other Tanzanian indigenous chickens based on mitochondrial DNA D-loop region. Front. Anim. Sci. 6:1732717. doi: 10.3389/fanim.2025.1732717

Received: 26 October 2025; Accepted: 29 December 2025; Revised: 29 December 2025;
Published: 21 January 2026.

Edited by:

Francisco Javier Navas González, University of Cordoba, Spain

Reviewed by:

Khanyisile Hadebe, Agricultural Research Council of South Africa (ARC-SA), South Africa
Boko Michel Orounladji, University of Abomey-Calavi, Benin

Copyright © 2026 Mhando, Mbaga, Nguluma, Mwega and Lyimo. 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: Zidihery Maquis Mhando, emlkaWhlcnlAZ21haWwuY29t

ORCID: Zidihery Maquis Mhando, orcid.org/0009-0008-0930-869X
Said Hemed Mbaga, orcid.org/0000-0003-4043-6629
Athumani Shabani Nguluma, orcid.org/0000-0002-1453-982X
Elisa Daniel Mwega, orcid.org/0000-0002-6192-4974
Charles Moses Lyimo, orcid.org/0000-0002-8739-0253

Disclaimer: 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.