Analysis of genetic diversity and conservation priorities for six north Ethiopian cattle breeds

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01 Mbp to 163.88 Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.

Genetic variability and biodiversity of Ukrainian Gray cattle by the BoLA-DRB3 gene

At the current stage of genetic studies of cattle, more and more attention is being drawn to autochthonous breeds. Native cattle have a number of prominent phenotypic traits and have preserved unique genes and their combinations lost by modern commercial breeds, which would be valuable to use in selective programs. We surveyed polymorphism of the Ukrainian autochthonous Gray breed according to alleles of exon 2 of the BoLA-DRB3 gene. The uniqueness of the gene lies in the broad variability of its allele variants. Significant informativeness at DNA level is quite important for genetic studies. We surveyed allele polymorphism using the PCR-RLFP method on DNA isolated from 88 samples of blood of cows and 5 samples of sperm. We identified 28 alleles, of which 23 variants were nomenclature ones and 5 (jba, *jab, *jbb, *nad and *nda) were “without established nomenclature”, their share accounting for 8.9%. Four alleles *06, *12, *16 and *jba had a frequency above 5% and occupied 69.9% of the breed’s allele fund overall. The commonest allele was BoLA-DRB3.2*16 (44.1%). In total, we found 40 genotypes. Considering the significant dominance of variant *16, as expected, 5 genotypes with its inclusion occurred: *16/*16, *12/*16, *06/*16, *16/*24 and *jba/*16. It was present in the genotype of two out three studied animals. Parameters of heterozygosity, effective number of alleles, Shannon and Pielou indices indicate that Ukrainian Gray cattle are characterized by lowest level of genetic variability and biodiversity according to the BoLA-DRB3 gene compared with other breeds. Due to significant dominance of allele *16, the breed has no inbred motifs. We noted deviation toward increase in homozygosity without deviations from the norm of the distribution according to Hardy-Weinberg equilibrium. The obtained results will be used for genetic-populational programs with the purpose of improving the genetic potential of cattle breeds in terms of economically beneficial traits and diseases of cattle.

Marginal diversity analysis of conservation of Chinese domestic duck breeds

The present study aimed to systematically evaluate the genetic diversity of Chinese domestic duck breeds and ensure the most effective allocation and usage of conservation funds. We first performed an analysis of DNA genetic distance in 21 duck breeds by measuring short tandem repeats. Then, we calculated the extinction probability, contribution rate, and marginal diversity for each breed. The results showed that the extinction rate of the Zhongshan duck, Guangxi duck, and Ji'an duck were the highest at 0.67, 0.59, and 0.59, respectively, and that of the Linwu duck, Jinding duck, and Gaoyou duck were the lowest at 0.15, 0.18, and 0.19, respectively. The current diversity of populations was 7.72 and the expected diversity in five hundred years is 5.14 ± 1.15. The marginal diversity of the Chinese Muscovy duck was the largest (-2.20), accounting for 42.61% of the expected diversity, followed by the Guangxi duck (-0.49, 9.44%), whereas the Jinding duck was the smallest (-0.12; 2.32%). The protection potency of the Chinese Muscovy duck was the largest (0.61), followed by Guangxi duck (0.29), whereas the Jinding duck was the smallest (0.02). This study provides a reference for determining the conservation priority of Chinese domestic duck breeds or genetic resources.

Prevalence and risk factors analysis of bovine tuberculosis in cattle raised in mixed crop-livestock farming system in Tigray region, Ethiopia

Bovine tuberculosis (BTB) is a disease of animal and public health importance in developing countries. In rural Ethiopia, there is potential for a shift in the epidemiologic of this disease driven by transformation of dairy industry. This includes gradual change from the traditional mixed crop-livestock husbandry practice to a semi-intensification system. It is therefore, essential to document the prevalence and risk factors of BTB to continuously update the designing and implementation of control and prevention strategies. Here, we present findings of a cross-sectional study on the prevalence and associated risk factors of BTB among cattle reared under mixed crop-livestock farming system in Tigray region, Ethiopia. A multistage purposive sampling approach was used to select districts, villages, herds and individual cattle. A total of 1,357 cattle from 310 herds were examined for BTB infection using a comparative intradermal tuberculin skin test (CIDT). Questionnaires were used to gather data on herd structure and herd management practices. A multilevel logistic mixed effect model was used to determine risk factors after accounting for clustering effect at three levels (village, herd and individual animal). Overall prevalence of BTB was 4.3% (95% CI = 3.4-5.6), with the highest prevalence recorded in Alamata district (5.6%) and lowest in Korem (1.6%). Multilevel logistic mixed effect model analysis identified exotic breed (OR = 3, p = 0.014), closed barn (OR = 2.6, p = 0.018), large herd size (OR = 2.6, p = 0.05) and purchase of cattle (OR = 2.1, p = 0.027) as important risk factors for BTB. Taken together, these findings suggest that the current dairy development program centred on the introduction of exotic and or crossed animals could have contributed to changing epidemiological situations of BTB in the study area.

Genetic diversity and admixture of indigenous cattle from North Ethiopia: implications of historical introgressions in the gateway region to Africa

Microsatellite variation was surveyed to determine the genetic diversity, population structure and admixture of seven North Ethiopian cattle breeds by combining multiple microsatellite data sets of Indian and West African zebu, and European, African and Near-Eastern taurine in genetic analyses. Based on allelic distribution, we identified four diagnostic alleles (HEL1-123 bp, CSSM66-201 bp, BM2113-150 bp and ILSTS6-285 bp) specific to the Near-Eastern taurine. Results of genetic relationship and population structure analyses confirmed the previously established marked genetic distinction between taurine and zebu, and indicated further divergence among the bio-geographical groupings of breeds such as North Ethiopian, Indian and West African zebu, and African, European and Near-Eastern taurine. Using the diagnostic alleles for bio-geographical groupings and a Bayesian method for population structure inference, we estimated the genetic influences of major historical introgressions in North Ethiopian cattle. The breeds have been heavily (>90%) influenced by zebu, followed by African, European and the Near-Eastern taurine. Overall, North Ethiopian cattle show a high level of within-population genetic variation (e.g. observed heterozygosity = 0.659-0.687), which is in the upper range of that reported for domestic cattle and indicates their potential for future breeding applications, even in a global context. Rather low but significant population differentiation (F(ST) = 1.1%, P < 0.05) was recorded as a result of multiple introgression events and strong genetic exchanges among the North Ethiopian breeds.

Y chromosome genetic diversity in the Lidia bovine breed: a highly fragmented population

To assess the paternal gene pool in the Lidia bovine breed (or fighting bull), a total of 603 animals belonging to 81 herds classified in 33 lineages were genotyped for six Y chromosome microsatellites, one single nucleotide polymorphism and one indel. A total of 10 haplotypes were determined with a high level of frequency variation between them, ranging from 0.2 to 74%. All the haplotypes identified belong to two previously defined major haplogroups (Y1 and Y2). Two major paternal influences were identified, corresponding to the two most common haplotypes (H1Y1 and H3Y2) with frequencies of 74 and 18%, respectively. The detection of the INRA189-104 allele evidenced an African influence in the Lidia bovine breed. Low levels of haplotype diversity have been achieved and only eight lineages showed more than one haplotype. Analysis of molecular variance showed a high level of interlineage variance (F(ST) = 86%). Network results evidenced two main clusters made for those haplotypes belonging to Y1 and Y2 haplogroups, respectively. The findings support a high level of genetic structure together with a low level of genetic diversity in the Lidia bovine breed.

Objectives, criteria and methods for using molecular genetic data in priority setting for conservation of animal genetic resources

The genetic diversity of the world's livestock populations is decreasing, both within and across breeds. A wide variety of factors has contributed to the loss, replacement or genetic dilution of many local breeds. Genetic variability within the more common commercial breeds has been greatly decreased by selectively intense breeding programmes. Conservation of livestock genetic variability is thus important, especially when considering possible future changes in production environments. The world has more than 7500 livestock breeds and conservation of all of them is not feasible. Therefore, prioritization is needed. The objective of this article is to review the state of the art in approaches for prioritization of breeds for conservation, particularly those approaches that consider molecular genetic information, and to identify any shortcomings that may restrict their application. The Weitzman method was among the first and most well-known approaches for utilization of molecular genetic information in conservation prioritization. This approach balances diversity and extinction probability to yield an objective measure of conservation potential. However, this approach was designed for decision making across species and measures diversity as distinctiveness. For livestock, prioritization will most commonly be performed among breeds within species, so alternatives that measure diversity as co-ancestry (i.e. also within-breed variability) have been proposed. Although these methods are technically sound, their application has generally been limited to research studies; most existing conservation programmes have effectively primarily based decisions on extinction risk. The development of user-friendly software incorporating these approaches may increase their rate of utilization.