Comparison of SNP variation and distribution in indigenous ethiopian and korean cattle (hanwoo) populations

Genetic diversity of DGAT1 gene linked to milk production in cattle populations of Ethiopia

BACKGROUND

Diacylglycerol acyl-CoA acyltransferase 1 (DGAT1) has become a promising candidate gene for milk production traits because of its important role as a key enzyme in catalyzing the final step of triglyceride synthesis. Thus use of bovine DGAT1 gene as milk production markers in cattle is well established. However, there is no report on polymorphism of the DGAT1 gene in Ethiopian cattle breeds. The present study is the first comprehensive report on diversity, evolution, neutrality evaluation and genetic differentiation of DGAT1 gene in Ethiopian cattle population. The aim of this study was to characterize the genetic variability of exon 8 region of DGAT1 gene in Ethiopian cattle breeds.

RESULTS

Analysis of the level of genetic variability at the population and sequence levels with genetic distance in the breeds considered revealed that studied breeds had 11, 0.615 and 0.010 haplotypes, haplotype diversity and nucleotide diversity respectively. Boran-Holstein showed low minor allele frequency and heterozygosity, while Horro showed low nucleotide and haplotype diversities. The studied cattle DGAT1 genes were under purifying selection. The neutrality test statistics in most populations were negative and statistically non-significant (p > 0.10) and consistent with a populations in genetic equilibrium or in expansion. Analysis for heterozygosity, polymorphic information content and inbreeding coefficient revealed sufficient genetic variation in DGAT1 gene. The pairwise FST values indicated significant differentiation among all the breeds (FST = 0.13; p ≤ 0.05), besides the rooting from the evolutionary or domestication history of the cattle inferred from the phylogenetic tree based on the neighbourhood joining method. There was four separated cluster among the studied cattle breeds, and they shared a common node from the constructed tree.

CONCLUSION

The cattle populations studied were polymorphic for DGAT1 locus. The DGAT1 gene locus is extremely crucial and may provide baseline information for in-depth understanding, exploitation of milk gene variation and could be used as a marker in selection programmes to enhance the production potential and to accelerate the rate of genetic gain in Ethiopian cattle populations exposed to different agro ecology condition.

Genome analyses revealed genetic admixture and selection signatures in Bos indicus

The genomic diversity and relationship among seven diverse cattle breeds viz. Sahiwal, Tharparkar, Gir, Vechur, Ongole, Kangayam and Hariana were investigated in 132 random samples based on high density SNP array comprising > 777 K SNPs. A total of 1993 SNPs (0.25% of the total) having greater power (F_ST ≥ 0.20) to differentiate these cattle populations were identified, and utilized to partition genome of each animal into a predefined number of clusters. The structure of these cattle indicated shared ancestry of dairy breeds viz. Gir, Tharparkar and Sahiwal. Most of the animals (> 76%) of different populations under study except Vechur clustered into their own group of animals called breed. Vechur population retained highest rate of admixture, consistent with its crossing with other breeds. Ongole, Kangayam and Hariana shared comparatively less of their genome (≤ 15%) with other breeds. The study indicated that all seven breeds evolved from their independent ancestry but there was intermixing of these breeds in the recent past. The selection signatures identified between draft (Kangayam) and dairy breeds included several genes like FAM19A2, RAB31P, BEST3, DGKA, AHCY, PIGU and PFKP which are involved in immune response, metabolic pathway, transportation of glucose and sugars, signaling pathways, cellular processes, cell division and glycolysis regulation, respectively. Moreover, these genomic regions also harbour QTLs affecting milk performance traits. The signatures were also identified even between the dairy breeds. In comparison to large-sized cattle, there were significant differences in the number of QTLs affecting production (body weight, growth rate etc.) and morphological traits (height) in short-statured Vechur breed. The presence of HMGA2 gene in the selection signature on chromosome 5 may explain the variations in stature between these cattle.

Single nucleotide polymorphisms of leptin gene in five Ethiopian indigenous cattle breeds and the Korean Hanwoo breed

Considering the escalating number of scientific reports on the association between the leptin gene and diverse physiological traits and performance of cattle populations, this study was directed towards identifying SNPs in the leptin gene among five indigenous cattle breeds of Ethiopia. DNA samples were extracted from the nasal swabs of the Ethiopian indigenous cattle breeds: Arsi (n = 18), Horro (n = 20), Begait (n = 21), Boran (n = 19), and Fogera (n = 17) and the Korean Hanwoo (a representative taurine breed) (n = 20), followed by PCR amplification of exon 2 and exon 3 regions of the leptin gene and sequence analysis of the PCR products. Five SNPs, two (generating missense mutations) on exon 2 and three (generating silent mutations) on exon 3 regions, were explicated in this study. Allele frequency and genotype frequency distribution pertaining to the SNPs were recorded for the studied cattle breeds besides the minor allele frequency and deviation from the Hardy-Weinberg equilibrium. Positive F_IS index values were recorded for all the markers except SNP2, illustrative of heterozygote deficiency. MEGA X software-based evolutionary divergence analysis of the phylogenetic tree based on the SNP data revealed that the large-sized breeds, Hanwoo, Begait, Boran, and Fogera, were more closely clustered compared to the small-sized Arsi breed. Among the seven haplotypes documented from the various breeds, sequence analysis was suggestive of haplotypes 1 and 2 to be ancestral haplotypes for the leptin gene. This study is envisaged to accelerate molecular breeding programs for the genetic improvement of the Ethiopian cattle breeds.

Genetic diversity and population structure of three native cattle populations in Mozambique

In the present work, the population diversity and structure of three populations of native Mozambican cattle were studied, to develop knowledge that is required for sound conservation and genetic improvement programs of these genetic resources. A total of 228 animals (Landim, Angone, and Tete) were genotyped using the International Dairy and Beef version three (IDBV3) SNP BeadChip array. Population parameters varied within a limited scope, with the average minor allele frequency (MAF) ranging from 0.228 ± 0.154 in the Angone to 0.245 ± 0.145 in the Tete population, while estimates of expected heterozygosities varied from 0.304 ± 0.166 in the Angone to 0.329 ± 0.148 in the Tete population. Low positive (0.065 ± 0.109) inbreeding rates were detected in the three cattle groups. Population structure and admixture analyses indicated low genetic differentiation and various degrees of admixture among the populations. The effective population size has decreased over time and at 12 generations ago ranged between 349 (Tete) and 929 (Landim). The average linkage disequilibrium (LD) of the studied populations ranged from 0.400 ± 0.213 (Tete) to 0.434 ± 0.232 (Landim). The findings of this study will be valuable for formulating management and conservation strategies for indigenous Mozambican cattle populations.

Single nucleotide polymorphisms at heat shock protein 90 gene and their association with thermo-tolerance potential in selected indigenous Nigerian cattle

Heat shock protein (HSP) 90 gene provides protection and adaptation to thermal assault and certain polymorphisms have been associated to heat tolerance in humans and animals. Single nucleotide polymorphisms (SNPs) of HSP 90 gene were used to evaluate the scientific basis of heat tolerance in four zebu breeds of Nigeria. The DNA was extracted from skin tissue of 90 adult bulls representing White Fulani (WF), Sokoto Gudali (SG), Red Bororo (RB), and Ambala (AM). The SNPs were determined in DNAs using PCR, sequencing, and visualization and bio-editing by chromatogram in SeqMan Ngen tool. Subsequently, respective genotypes were constructed and genotypic and allelic frequencies were computed. Also, body parameters related to heat stress (HS) including body temperature (BT), rectal temperature (RT), and respiratory rates (RR) were taken for each animal before biological sampling and heat tolerance coefficient (HTC) was calculated. We detected four SNPs distinct/specific for each breed as follows: change from thymine (T) to guanine (G) at position 116 (T116G) in RB, G to cytosine (C) at 220 (G220C) in SG, G to adenine (A) at two positions, 346 (G346A) and 390 (G390A) in AM and WF, respectively. Heterozygous SNPs showed significantly lower values (P < 0.0001) for BT, RT, RR, and HTC than homozygous genotypes at all positions. We hypothesize that animals with heterozygous SNPs in exon 3 of HSP 90 may be tolerant to HS. These SNPs can be used as bio-markers for screening large populations of cattle for tolerance to hot tropical conditions in Nigeria and other sub-humid places.

Evaluation of two bovine SNP genotyping arrays for breed clustering and stratification analysis in well-known taurine and indicine breeds

The present study aimed to evaluate the efficiency of two Bovine SNP genotyping arrays (i.e., 50 K and HD) for breed clustering and stratification related studies in taurine and indicine breeds. The whole-genome SNP data at two densities were assembled into three datasets (A, B and C). Dataset A (N = 213) included 50 K genotypic data for five taurine (Holstein-Friesian, Guernsey, Brown Swiss, Angus and Jersey) and two indicine (Gir and Nellore) breeds. Dataset B (N = 241) included the same breeds with HD density data. Dataset C (N = 299) included 50 K SNP genotypic data for six taurine (Holstein-Friesian, Jersey, Guernsey, Brown Swiss, Angus and Hereford) and six indicine (Hariana, Kankrej, Brahman, Nellore, Sahiwal and Gir) breeds. The analysis was done using ADMIXTURE program (bioinformatics-based) and cross-validation errors and Principal Component Analysis (statistical analysis). The proportion of polymorphic markers and minor allele frequencies were assessed for each breed. The proportion of markers polymorphic was consistently higher in taurine breeds when compared with breeds from indicine group. Minor allele frequency estimates and ADMIXTURE results showed differential patterns for both the lineages. However, no significant increase in the accuracy of genomic clustering was found on moving from 50 K to HD density data.

Genetic diversities and population structures of four popular Myanmar local cattle breeds

Local cattle breeds are popular in Myanmar because they make excellent draught animals, so it is of fundamental importance that their genetic diversities and population structures are elucidated for breeding and conservation purposes. In this study, we characterized the genetic diversities and population structures of four popular Myanmar local cattle breeds and five native cattle populations from neighboring countries (Bangladesh, Laos, Cambodia, and Vietnam) using the GGP Bovine 50K array. The mean observed and expected heterozygosity estimates using pruned datasets ranged from 0.317 and 0.322, respectively, in Cambodia to 0.448 and 0.421, respectively, in Vietnam South. The four Myanmar local breeds exhibited similar levels of genetic diversity. However, AMOVA revealed significant genetic differentiation among the nine populations (p < 0.00001), and Neighbor-Net analysis showed that Shan Ngwar Pu was distinct from the other Myanmar local breeds. Furthermore, PCA and population structure analyses revealed that the native cattle from neighboring countries genetically influenced the Myanmar local breeds to some extent and that the genetic origins could also be observed in the local breeds. These findings highlight the importance of Myanmar local breeds as genetic resources and provide useful information for the future development of breeding strategies and conservation management plans.

Evaluation of Bovine High-Density SNP Genotyping Array in Indigenous Dairy Cattle Breeds

In total 52 samples of Sahiwal ( 19 ), Tharparkar ( 17 ), and Gir ( 16 ) were genotyped by using BovineHD SNP chip to analyze minor allele frequency (MAF), genetic diversity, and linkage disequilibrium among these cattle. The common SNPs of BovineHD and 54K SNP Chips were also extracted and evaluated for their performance. Only 40%-50% SNPs of these arrays was found informative for genetic analysis in these cattle breeds. The overall mean of MAF for SNPs of BovineHD SNPChip was 0.248 ± 0.006, 0.241 ± 0.007, and 0.242 ± 0.009 in Sahiwal, Tharparkar and Gir, respectively, while that for 54K SNPs was on lower side. The average Reynold's genetic distance between breeds ranged from 0.042 to 0.055 based on BovineHD Beadchip, and from 0.052 to 0.084 based on 54K SNP Chip. The estimates of genetic diversity based on HD and 54K chips were almost same and, hence, low density chip seems to be good enough to decipher genetic diversity of these cattle breeds. The linkage disequilibrium started decaying (r² < 0.2) at 140 kb inter-marker distance and, hence, a 20K low density customized SNP array from HD chip could be designed for genomic selection in these cattle else the 54K Bead Chip as such will be useful.

Determination of Genetic Structure and Signatures of Selection in Three Strains of Tanzania Shorthorn Zebu, Boran and Friesian Cattle by Genome-Wide SNP Analyses

Background

More than 90 percent of cattle in Tanzania belong to the indigenous Tanzania Short Horn Zebu (TSZ) population which has been classified into 12 strains based on historical evidence, morphological characteristics, and geographic distribution. However, specific genetic information of each TSZ population has been lacking and has caused difficulties in designing programs such as selection, crossbreeding, breed improvement or conservation. This study was designed to evaluate the genetic structure, assess genetic relationships, and to identify signatures of selection among cattle of Tanzania with the main goal of understanding genetic relationship, variation and uniqueness among them.

Methodology/Principal findings

The Illumina Bos indicus SNP 80K BeadChip was used to genotype genome wide SNPs in 168 DNA samples obtained from three strains of TSZ cattle namely Maasai, Tarime and Sukuma as well as two comparative breeds; Boran and Friesian. Population structure and signatures of selection were examined using principal component analysis (PCA), admixture analysis, pairwise distances (F_ST), integrated haplotype score (iHS), identical by state (IBS) and runs of homozygosity (ROH). There was a low level of inbreeding (F~0.01) in the TSZ population compared to the Boran and Friesian breeds. The analyses of F_ST, IBS and admixture identified no considerable differentiation between TSZ trains. Importantly, common ancestry in Boran and TSZ were revealed based on admixture and IBD, implying gene flow between two populations. In addition, Friesian ancestry was found in Boran. A few common significant iHS were detected, which may reflect influence of recent selection in each breed or strain.

Conclusions

Population admixture and selection signatures could be applied to develop conservation plan of TSZ cattle as well as future breeding programs in East African cattle.

Genome-wide Single Nucleotide Polymorphism Analyses Reveal Genetic Diversity and Structure of Wild and Domestic Cattle in Bangladesh

In spite of variation in coat color, size, and production traits among indigenous Bangladeshi cattle populations, genetic differences among most of the populations have not been investigated or exploited. In this study, we used a high-density bovine single nucleotide polymorphism (SNP) 80K Bead Chip derived from Bos indicus breeds to assess genetic diversity and population structure of 2 Bangladeshi zebu cattle populations (red Chittagong, n = 28 and non-descript deshi, n = 28) and a semi-domesticated population (gayal, n = 17). Overall, 95% and 58% of the total SNPs (69,804) showed polymorphisms in the zebu and gayal populations, respectively. Similarly, the average minor allele frequency value was as high 0.29 in zebu and as low as 0.09 in gayal. The mean expected heterozygosity varied from 0.42±0.14 in zebu to 0.148±0.14 in gayal with significant heterozygosity deficiency of 0.06 (F_IS) in the latter. Coancestry estimations revealed that the two zebu populations are weakly differentiated, with over 99% of the total genetic variation retained within populations and less than 1% accounted for between populations. Conversely, strong genetic differentiation (F_ST = 0.33) was observed between zebu and gayal populations. Results of population structure and principal component analyses suggest that gayal is distinct from Bos indicus and that the two zebu populations were weakly structured. This study provides basic information about the genetic diversity and structure of Bangladeshi cattle and the semi-domesticated gayal population that can be used for future appraisal of breed utilization and management strategies.

Genetic diversity, population structure and relationships in indigenous cattle populations of Ethiopia and Korean Hanwoo breeds using SNP markers

In total, 166 individuals from five indigenous Ethiopian cattle populations – Ambo (n = 27), Borana (n = 35), Arsi (n = 30), Horro (n = 36), and Danakil (n = 38) – were genotyped for 8773 single nucleotide polymorphism (SNP) markers to assess genetic diversity, population structure, and relationships. As a representative of taurine breeds, Hanwoo cattle (n = 40) were also included in the study for reference. Among Ethiopian cattle populations, the proportion of SNPs with minor allele frequencies (MAFs) ≥0.05 ranged from 81.63% in Borana to 85.30% in Ambo, with a mean of 83.96% across all populations. The Hanwoo breed showed the highest proportion of polymorphism, with MAFs ≥0.05, accounting for 95.21% of total SNPs. The mean expected heterozygosity varied from 0.370 in Danakil to 0.410 in Hanwoo. The mean genetic differentiation (F_ST; 1%) in Ethiopian cattle revealed that within individual variation accounted for approximately 99% of the total genetic variation. As expected, F_ST and Reynold genetic distance were greatest between Hanwoo and Ethiopian cattle populations, with average values of 17.62 and 18.50, respectively. The first and second principal components explained approximately 78.33% of the total variation and supported the clustering of the populations according to their historical origins. At K = 2 and 3, a considerable source of variation among cattle is the clustering of the populations into Hanwoo (taurine) and Ethiopian cattle populations. The low estimate of genetic differentiation (F_ST) among Ethiopian cattle populations indicated that differentiation among these populations is low, possibly owing to a common historical origin and high gene flow. Genetic distance, phylogenic tree, principal component analysis, and population structure analyses clearly differentiated the cattle population according to their historical origins, and confirmed that Ethiopian cattle populations are genetically distinct from the Hanwoo breed.