Mitochondrial DNA polymorphism in Moroccan goats

Insights into the genomic homogeneity of Moroccan indigenous sheep breeds though the lens of runs of homozygosity

Numerous studies have indicated that Morocco's indigenous sheep breeds are genetically homogenous, posing a risk to their survival in the challenging harsh climate conditions where they predominantly inhabit. To understand the genetic behind genetic homogeneity through the lens of runs of homozygosity (ROH), we analyzed the whole genome sequences of five indigenous sheep breeds (Beni Guil, Ouled Djellal, D'man, Sardi, Timahdite and Admixed).The results from principal component, admixture, Fst, and neighbour joining tree analyses consistently showed a homogenous genetic structure. This structure was characterized by an average length of 1.83 Mb for runs of homozygosity (ROH) segments, with a limited number of long ROH segments (24-48 Mb and > 48 Mb). The most common ROH segments were those ranging from 1-6 Mb. The most significant regions of homozygosity (ROH Islands) were mostly observed in two chromosomes, namely Chr1 and Chr5. Specifically, ROH Islands were exclusively discovered in the Ouled Djellal breed on Chr1, whereas Chr5 exhibited ROH Islands in all breeds. The analysis of ROH Island and iHS technique was employed to detect signatures of selection on Chr1 and Chr5. The results indicate that Chr5 had a high level of homogeneity, with the same genes being discovered across all breeds. In contrast, Chr1 displays some genetic variances between breeds. Genes identified on Chr5 included SLC39A1, IL23A, CAST, IL5, IL13, and IL4 which are responsible for immune response while genes identified on Chr1 include SOD1, SLAMF9, RTP4, CLDN1, and PRKAA2. ROH segment profile and effective population sizes patterns suggests that the genetic uniformity of studied breeds is the outcome of events that transpired between 250 and 300 generations ago. This research not only contributes to the understanding of ROH distribution across breeds but helps design and implement native sheep breeding and conservation strategies in Morocco. Future research, incorporating a broader sample size and utilizing the pangenome for reference, is recommended to further elucidate these breeds' genomic landscapes and adaptive mechanisms.

Mitochondrial DNA D-loop sequence analysis reveals high variation and multiple maternal origins of indigenous Tanzanian goat populations

The Small East African (SEA) goat are widely distributed in different agro-ecological zones of Tanzania. We report the genetic diversity, maternal origin, and phylogenetic relationship among the 12 Tanzanian indigenous goat populations, namely Fipa, Songwe, Tanga, Pwani, Iringa, Newala, Lindi, Gogo, Pare, Maasai, Sukuma, and Ujiji, based on the mitochondrial DNA (mtDNA) D-loop. High haplotype (H d = 0.9619-0.9945) and nucleotide (π = 0.0120-0.0162) diversities were observed from a total of 389 haplotypes. The majority of the haplotypes (n = 334) belonged to Haplogroup A which was consistent with the global scenario on the genetic pattern of maternal origin of all goat breeds in the world. Haplogroup G comprised of 45 haplotypes drawn from all populations except the Ujiji goat population while Haplogroup B with 10 haplotypes was dominated by Ujiji goats (41%). Tanzanian goats shared four haplotypes with the Kenyan goats and two with goats from South Africa, Namibia, and Mozambique. There was no sharing of haplotypes observed between individuals from Tanzanian goat populations with individuals from North or West Africa. The indigenous goats in Tanzania have high genetic diversity defined by 389 haplotypes and multiple maternal origins of haplogroup A, B, and G. There is a lot of intermixing and high genetic variation within populations which represent an abundant resource for selective breeding in the different agro-ecological regions of the country.

Characterizing the Mitochondrial Diversity of Arbi Goats from Tunisia

Arbi is one of the main local goat breeds in Tunisia, representing an important economic resource in arid and hot areas where cattle and sheep cannot thrive successfully. In the current work, we have characterized the mitochondrial diversity of 26 Arbi goats by partially sequencing the mitochondrial D-loop region. These sequences plus 10 retrieved from GenBank were analyzed with the DnaSP v.5.10.1, evidencing the existence of 12 different haplotypes. Nucleotide and haplotype diversities were 0.02 and 0.96. Moreover, median-joining network analysis showed that all D-loop sequences from Arbi goats correspond to haplogroup A and that in general they do not cluster with sequences from other goat breeds. The high diversity that has been observed in North African goats is compatible with the maritime diffusion of the Neolithic package 10,000-7000 YBP. Moreover, there are evidences that local Tunisian breeds have been extensively crossed with highly productive transboundary breeds in order to improve meat and milk yields. These uncontrolled crossing practices may lead to the loss of alleles that play key roles in the adaptation of Tunisian local breeds to a harsh environment.

Genome-wide diversity and demographic dynamics of Cameroon goats and their divergence from east African, north African, and Asian conspecifics

Indigenous goats make significant contributions to Cameroon’s national and local economy, but little effort has been devoted to identifying the populations. Here, we assessed the genetic diversity and demographic dynamics of Cameroon goat populations using mitochondrial DNA (two populations) and autosomal markers (four populations) generated with the Caprine 50K SNP chip. To infer genetic relationships at continental and global level, genotype data on six goat populations from Ethiopia and one population each from Egypt, Morocco, Iran, and China were included in the analysis. The mtDNA analysis revealed 83 haplotypes, all belonging to haplogroup A, in Cameroon goats. Four haplotypes were shared between goats found in Cameroon, Mozambique, Namibia, Zimbabwe, Kenya, and Ethiopia. Analysis of autosomal SNPs in Cameroon goats revealed the lowest H_O (0.335±0.13) and H_E (0.352±0.15) in the North-west Highland and Central Highland populations, respectively. Overall, the highest H_O (0.401±0.12) and H_E (0.422±0.12) were found for Barki and Iranian goats, respectively. Barki goats had the highest average MAF, while Central Highland Cameroon goats had the lowest. Overall, Cameroon goats demonstrated high F_IS. AMOVA revealed that 13.29% of the variation was explained by genetic differences between the six population groups. Low average F_ST (0.01) suggests intermixing among Cameroon goats. All measures indicated that Cameroon goats are closer to Moroccan goats than to other goat populations. PCA and STRUCTURE analyses poorly differentiated the Cameroon goats, as did genetic distance, Neighbor-Net network, and neighbor-joining tree analyses. The haplotype analysis of mtDNA showed the initial dispersion of goats to Cameroon and central Africa from north-east Africa following the Nile Delta. Whereas, the approximate Bayesian computation indicated Cameroon goats were separated from Moroccan goats after 506 generations in later times (~1518 YA), as supported by the phylogenetic net-work and admixture outputs. Overall, indigenous goats in Cameroon show weak phylogenetic structure, suggesting either extensive intermixing.

Genetic homogeneity of North-African goats

North Africa represents a rich and early reservoir of goat genetic diversity, from which the main African breeds have been derived. In this study, the genetic diversity of four indigenous Algerian goat breeds (i.e., Arabia, Makatia, M’Zabite and Kabyle, with n = 12 for each breed) has been investigated for the first time by genome-wide SNP genotyping; moreover in a broader context, genetic structuration of Algerian and Moroccan goats was explored (via F_ST, MDS, STRUCTURE, FineSTRUCTURE, BAPS, sPCA and DAPC analyses). At national level, the study revealed high level of genetic diversity and a significant phenomenon of admixture affecting all the Algerian breeds. At broader scale, clear global genetic homogeneity appeared considering both Algerian and Moroccan stocks. Indeed, genetic structuration was almost nonexistent among Arabia (from Algeria), Draa, Black and Nord (from Morocco), while the ancestral Kabyle and M’Zabite breeds, reared by Berber peoples, showed genetic distinctness. The study highlighted the threat to the Maghrebin stock, probably induced by unsupervised cross-breeding practices which have intensified in recent centuries. Moreover, it underlined the necessity to deepen our understanding of the genetic resources represented by the resilient North-African goat stock.

Characterizing neutral genomic diversity and selection signatures in indigenous populations of Moroccan goats (Capra hircus) using WGS data

Since the time of their domestication, goats (Capra hircus) have evolved in a large variety of locally adapted populations in response to different human and environmental pressures. In the present era, many indigenous populations are threatened with extinction due to their substitution by cosmopolitan breeds, while they might represent highly valuable genomic resources. It is thus crucial to characterize the neutral and adaptive genetic diversity of indigenous populations. A fine characterization of whole genome variation in farm animals is now possible by using new sequencing technologies. We sequenced the complete genome at 12× coverage of 44 goats geographically representative of the three phenotypically distinct indigenous populations in Morocco. The study of mitochondrial genomes showed a high diversity exclusively restricted to the haplogroup A. The 44 nuclear genomes showed a very high diversity (24 million variants) associated with low linkage disequilibrium. The overall genetic diversity was weakly structured according to geography and phenotypes. When looking for signals of positive selection in each population we identified many candidate genes, several of which gave insights into the metabolic pathways or biological processes involved in the adaptation to local conditions (e.g., panting in warm/desert conditions). This study highlights the interest of WGS data to characterize livestock genomic diversity. It illustrates the valuable genetic richness present in indigenous populations that have to be sustainably managed and may represent valuable genetic resources for the long-term preservation of the species.