The present study is the first comprehensive report on diversity, population structure, genetic admixture and mitochondrial DNA variation in South Indian draught type zebu cattle. The diversity of ...South Indian cattle was moderately high. A significantly strong negative correlation coefficient of -0.674 (P6.25%) was observed in Punganur, Vechur, Umblachery and Pulikulam cattle breeds. Two major maternal haplogroups, I1 and I2, typical of zebu cattle were observed, with the former being predominant than the later. The pairwise differences among the I2 haplotypes of South Indian cattle were relatively higher than West Indian (Indus valley site) zebu cattle. The results indicated the need for additional sampling and comprehensive analysis of mtDNA control region variations to unravel the probable location of origin and domestication of I2 zebu lineage. The present study also revealed major concerns on South Indian zebu cattle (i) risk of endangerment due to small effective population size and high rate of inbreeding (ii) lack of sufficient purebred zebu bulls for breeding and (iii) increasing level of taurine admixture in zebu cattle. Availability of purebred semen for artificial insemination, incorporation of genomic/molecular information to identify purebred animals and increased awareness among farmers will help to maintain breed purity, conserve and improve these important draught cattle germplasms of South India.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on ...the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.
The cow mammary transcriptome was explored at -30, -15, 1, 15, 30, 60, 120, 240, and 300 d relative to parturition. A total of 6,382 differentially expressed genes (DEG) at a false discovery rate ≤ ...0.001 were found throughout lactation. The greatest number of DEG (>3,500 DEG) was observed at 60 and 120 d vs. -30 d with the largest change between consecutive time points observed at -15 vs. 1 d and 120 vs. 240 d. Functional analysis of microarray data was performed using the Dynamic Impact Approach (DIA). The DIA analysis of KEGG pathways uncovered as the most impacted and induced 'Galactose metabolism', 'Glycosylphosphatidylinositol (GPI)-anchor biosynthesis', and 'PPAR signaling'; whereas, 'Antigen processing and presentation' was among the most inhibited. The integrated interpretation of the results suggested an overall increase in metabolism during lactation, particularly synthesis of carbohydrates and lipid. A marked degree of utilization of amino acids as energy source, an increase of protein export, and a decrease of the protein synthesis machinery as well cell cycle also were suggested by the DIA analysis. The DIA analysis of Gene Ontology and other databases uncovered an induction of Golgi apparatus and angiogenesis, and the inhibition of both immune cell activity/migration and chromosome modifications during lactation. All of the highly-impacted and activated functions during lactation were evidently activated at the onset of lactation and inhibited when milk production declined. The overall analysis indicated that the bovine mammary gland relies heavily on a coordinated transcriptional regulation to begin and end lactation. The functional analysis using DIA underscored the importance of genes associated with lactose synthesis, lipid metabolism, protein synthesis, Golgi, transport, cell cycle/death, epigenetic regulation, angiogenesis, and immune function during lactation.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The present study reports the population structure, genetic admixture and phylogeography of cattle breeds of Sri Lanka viz. Batu Harak, Thawalam and White cattle. Moderately high level of genetic ...diversity was observed in all the three Sri Lankan zebu cattle breeds. Estimates of inbreeding for Thawalam and White cattle breeds were relatively high with 6.1% and 7.2% respectively. Genetic differentiation of Sri Lankan Zebu (Batu Harak and White cattle) was lowest with Red Sindhi among Indus Valley Zebu while it was lowest with Hallikar among the South Indian cattle. Global F statistics showed 6.5% differences among all the investigated Zebu cattle breeds and 1.9% differences among Sri Lankan Zebu breeds. The Sri Lankan Zebu cattle breeds showed strong genetic relationships with Hallikar cattle, an ancient breed considered to be ancestor for most of the Mysore type draught cattle breeds of South India. Genetic admixture analysis revealed high levels of breed purity in Lanka White cattle with >97% Zebu ancestry. However, significant taurine admixture was observed in Batu Harak and Thawalam cattle. Two major Zebu haplogroups, I1 and I2 were observed in Sri Lankan Zebu with the former predominating the later in all the three breeds. A total of 112 haplotypes were observed in the studied breeds, of which 50 haplotypes were found in Sri Lankan Zebu cattle. Mismatch analysis revealed unimodal distribution in all the three breeds indicating population expansion. The sum of squared deviations (SSD) and raggedness index were non-significant in both the lineages of all the three breeds except for I1 lineage of Thawalam cattle (P<0.01) and I2 lineage of Batu Harak cattle (P<0.05). The results of neutrality tests revealed negative Tajima’s D values for both the lineages of Batu Harak (P>0.05) and White cattle (P>0.05) indicating an excess of low frequency polymorphisms and demographic expansion. Genetic dilution of native Zebu cattle germplasm observed in the study is a cause for concern. Hence, it is imperative that national breeding organizations consider establishing conservation units for the three native cattle breeds to maintain breed purity and initiate genetic improvement programs.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The overrepresented approach (ORA) is the most widely-accepted method for functional analysis of microarray datasets. The ORA is computationally-efficient and robust; however, it suffers from the ...inability of comparing results from multiple gene lists particularly with time-course experiments or those involving multiple treatments. To overcome such limitation a novel method termed Dynamic Impact Approach (DIA) is proposed. The DIA provides an estimate of the biological impact of the experimental conditions and the direction of the impact. The impact is obtained by combining the proportion of differentially expressed genes (DEG) with the log2 mean fold change and mean -log P-value of genes associated with the biological term. The direction of the impact is calculated as the difference of the impact of up-regulated DEG and down-regulated DEG associated with the biological term. The DIA was validated using microarray data from a time-course experiment of bovine mammary gland across the lactation cycle. Several annotation databases were analyzed with DIA and compared to the same analysis performed by the ORA. The DIA highlighted that during lactation both BTA6 and BTA14 were the most impacted chromosomes; among Uniprot tissues those related with lactating mammary gland were the most positively-impacted; within KEGG pathways 'Galactose metabolism' and several metabolism categories related to lipid synthesis were among the most impacted and induced; within Gene Ontology "lactose biosynthesis" among Biological processes and "Lactose synthase activity" and "Stearoyl-CoA 9-desaturase activity" among Molecular processes were the most impacted and induced. With the exception of the terms 'Milk', 'Milk protein' and 'Mammary gland' among Uniprot tissues and SP_PIR_Keyword, the use of ORA failed to capture as significantly-enriched (i.e., biologically relevant) any term known to be associated with lactating mammary gland. Results indicate the DIA is a biologically-sound approach for analysis of time-course experiments. This tool represents an alternative to ORA for functional analysis.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
•12 Genotypes of three Haemonchus species circulate among ruminants in Pakistan.•High rate of gene flow among Haemonchus parasites infecting different hosts.•H. placei isolates from Pakistan were ...found to be genetically distinct from others.•No significant deviation of Pak Haemonchus populations from selective neutrality.•7.86% frequency of F200Y β-tubulin allele associated with benzimidazole resistance.
Haemonchus species are major gastro-intestinal parasites affecting ruminants across the world. The present study aimed to assess the sympatric species distribution, genetic diversity, population structure and frequency of β-tubulin isotype 1 alleles associated with benzimidazole resistance. Internal transcribed spacer 2 (ITS2) sequences revealed three sympatric species of Haemonchus, H. contortus, H. placei and H. longistipes with 12 distinct genotypes circulating among ruminant hosts in Pakistan. High genetic variability was observed in Pakistani Haemonchus isolates at nicotine amide dehydrogenase subunit 4 (ND4) and cytochrome oxidase subunit 1 (COI) gene loci. Intra-population diversity parameters were higher in H. contortus isolates than H. placei. Phylogenetic analysis of ND4 and COI sequences did not reveal clustering of haplotypes originating from a particular host indicating high rate of gene flow among Haemonchus parasites infecting sheep, goat and cattle in Pakistan. ND4 and COI haplotypes from Pakistan were compared to sequences of Haemonchus isolates from 11 countries to elucidate the population structure. Multidimensional scaling (MDS) plot of pairwise FST derived from 531 ND4 haplotypes revealed clustering together of H. contortus from Pakistan, China, Malaysia and Italy while the isolates from Yemen and United States were found to be genetically distinct. With respect to H. placei, isolates from Pakistan were found to be genetically differentiated from isolates of other countries. The tests for selective neutrality revealed negative D statistics and did not reveal significant deviations in Pakistani Haemonchus populations while significant deviation (P<0.05) was observed in Brazilian and Chinese H. contortus populations. Median Joining (MJ) network of ND4 haplotypes revealed Yemenese H. contortus being closer to H. placei cluster. β-tubulin isotype 1 genotyping revealed 7.86% frequency of Y allele associated with benzimidazole resistance at F200Y locus in Pakistani Haemonchus isolates.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
Sheep chromosome 3 (Oar3) has the largest number of QTLs reported to be significantly associated with resistance to gastro-intestinal nematodes. This study aimed to identify single nucleotide ...polymorphisms (SNPs) within candidate genes located in sheep chromosome 3 as well as genes involved in major immune pathways. A total of 41 SNPs were identified across 38 candidate genes in a panel of unrelated sheep and genotyped in 713 animals belonging to 22 breeds across Asia, Europe and South America. The variations and evolution of immune pathway genes were assessed in sheep populations across these macro-environmental regions that significantly differ in the diversity and load of pathogens. The mean minor allele frequency (MAF) did not vary between Asian and European sheep reflecting the absence of ascertainment bias. Phylogenetic analysis revealed two major clusters with most of South Asian, South East Asian and South West Asian breeds clustering together while European and South American sheep breeds clustered together distinctly. Analysis of molecular variance revealed strong phylogeographic structure at loci located in immune pathway genes, unlike microsatellite and genome wide SNP markers. To understand the influence of natural selection processes, SNP loci located in chromosome 3 were utilized to reconstruct haplotypes, the diversity of which showed significant deviations from selective neutrality. Reduced Median network of reconstructed haplotypes showed balancing selection in force at these loci. Preliminary association of SNP genotypes with phenotypes recorded 42 days post challenge revealed significant differences (P<0.05) in fecal egg count, body weight change and packed cell volume at two, four and six SNP loci respectively. In conclusion, the present study reports strong phylogeographic structure and balancing selection operating at SNP loci located within immune pathway genes. Further, SNP loci identified in the study were found to have potential for future large scale association studies in naturally exposed sheep populations.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The diversity of goats in Myanmar is represented by three indigenous breeds, Jade Ni, Nyaung Oo and Waithar Li. This study aimed at characterizing the genetic diversity and relationship of Myanmar ...goat breeds using microsatellite and mitochondrial DNA variations. A total of 147 goats from all three indigenous breeds were genotyped at 27 microsatellite loci. Genetic diversity in terms of allelic polymorphisms, observed and expected heterozygosities were moderately high. The mean observed heterozygosity within breeds varied between 0.566±0.183 (Nyaung Oo) and 0.595±0.182 (Waithar Li) while the expected heterozygosity varied from 0.605±0.181 (Jade Ni) to 0.647±0.176 (Waithar Li). Considerable heterozygosity deficit ranging from 5.5% to 8.2% was observed in Myanmar goat breeds. Wright’s F statistics revealed most of the variations within breeds and only 1.9% of the total observed variation was explained by between breed differences. Principal components and Bayesian clustering analyses showed complete admixture of Nyaung Oo and Waithar Li goats indicating high rate of gene flow among these populations. Population stratification was observed in Jade Ni with a subset of individuals clustering distinctly. Variations in mitochondrial DNA control region revealed 22 distinct haplotypes belonging to two major haplogroups A and B. Haplogroup A was found to predominate Myanmar goats similar to other goat populations in Asia. Comparative analysis of mtDNA variations indicated possible Chinese origin of the maternal haplotypic lineages of Myanmar goats.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP
Haemonchosis is a major economic problem in goat production in humid, tropical and subtropical regions. The disease is caused by an abomasal nematode,
, which is highly pathogenic in small ruminants
...The aim of this study was to identifying single-nucleotide polymorphisms (SNP) that were associated with fecal egg counts (FEC) and could be used as markers to identify resistance to
in goats.
Ten novel variants in the
, and
genes were associated with FEC in goats with a nominal significance level of
< 0.05. Two missense mutation in the exon region of the caprine
gene resulted in replacement of arginine with cysteine at position 9473550 (R9473550C) and aspartic acid with glutamic acid at position 9473870 (D9473870E). Chinese goat breeds had significantly higher FEC than Bangladeshi goat breeds within their respective genotypes. Polymorphism information content (
), effective allele number (
), and heterozygosity (
) were greatest for the
_197_A > G SNP locus in all goat breeds. Pairwise coefficients of linkage disequilibrium (D´,
) revealed complete LD (
= 1) between significant SNP polymorphisms in
and
and strong LD (
= 0.93 and 0.98) between polymorphisms in
and
, respectively. Correlation coefficient (
) between FEC and body weight (BW) was significantly positive (
= 0.56***,
< 0.001) but that between FEC and packed cell volume (PCV) was negatively significant (
= - 0.47**,
< 0.01) in the total population of goats. On the other hand, correlation coefficient (
) between BW and PCV was not significant in total population of goats. Association analysis revealed that haplotypes within
, and
were significantly associated with FEC. Quantitative real-time PCR revealed that the relative expression of mRNA was higher (
< 0.001) for resistant, compared to susceptible, groups of goats for all candidate genes except
.
This study identified SNP markers that can potentially be used in marker-assisted selection programs to develop goat breeds that are resistant to
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
During the last few decades, the effective population size of indigenous zebu cattle breeds has declined drastically, resulting in the classification of some of them into the vulnerable, endangered, ...or critically endangered category. Drastic reductions in the effective size of a population may result in genetic bottlenecks and can affect within-breed genetic variability and its viability. The present study was undertaken with the objective of evaluating South Indian zebu cattle populations for mutation drift equilibrium and to detect the occurrence of recent genetic bottleneck events. A total of 293 cattle from eight indigenous breeds were genotyped at 27 FAO/ISAG-recommended microsatellite marker loci. Three different statistical tests, viz., the sign test, standardized differences test, and Wilcoxon sign rank test were performed using allele frequency data to detect loci with heterozygosity excess under the infinite alleles, stepwise, and two-phase mutation models. Under the infinite alleles model, the observed number of loci with heterozygosity excess (He > Heq) ranged between 10 and 19 among the investigated cattle breeds. However, the observed heterozygosity excess was not statistically significant (p > 0.05) in any of the studied breeds. Similarly, the standardized differences test and Wilcoxon sign rank test revealed no concrete evidence for the occurrence of a recent genetic bottleneck in South Indian zebu cattle breeds. The qualitative test for mode-shift distortion revealed a normal L-shaped distribution of allele frequencies, suggesting a lack of evidence for the loss of low-frequency alleles in all the investigated South Indian zebu cattle breeds.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK