Microphthalmia-associated transcription factor mutations are associated with white-spotted coat color in swamp buffalo

A candidate gene analysis of the microphthalmia-associated transcription factor (MITF) gene was used in an attempt to identify the genetic basis for a white-spotted coat color phenotype in the Asian swamp buffalo (Bubalus bubalis carabanensis). Ninety-three buffaloes-32 solid, 38 spotted and 23 white individuals-were Sanger-sequenced for all MITF exons as well as highly conserved intronic and flanking regions. MITF cDNA representing skin and iris tissue from six spotted, nine solid and one white buffaloes was also Sanger-sequenced to confirm detected mutations. Two independent loss-of-function mutations, a premature stop codon (c.328C>T, p.Arg110*) and a donor splice-site mutation (c.840+2T>A, p.Glu281_Leu282Ins8), both of which cause white-spotted coat color in swamp buffaloes, were identified. The nonsense mutation leads to a premature stop codon in exon 3, and likely removal of the resulting mRNA via nonsense-mediated decay pathway, whereas the donor splice-site mutation leads to aberrant splicing of exon 8 that encodes part of a highly conserved region of MITF. The resulting insertion of eight amino acid residues is expected to perturb the leucine zipper part in the basic helix-loop-helix leucine zipper (bHLH-Zip) domain and will most likely influence dimerization and DNA binding capacity. Electrophoretic mobility shift assay was performed using mutant and wild-type MITF proteins and showed that the mutant MITF protein resulting from the splice-site mutation decreased in vitro DNA binding capacity compared to wild-type MITF. White-spotted buffalo bulls are sacrificed in funeral ceremonies in Tana Toraja, Indonesia, because they are considered holy, and our results show that genetic variation causes a tie to the cultural use of these buffaloes.

MARKER ASSISTED EVALUATION OF MORPHOLOGICAL AND GENETIC ATTRIBUTES OF SUB-POPULATIONS OF NILI-RAVI BUFFALO: A VULNERABLE DAIRY TYPE RIVERINE BREED OF INDIA

In the present study, we report the distribution of true to type and atypical Nili-Ravi buffalo, a vulnerable dairy type riverine breed of North India and its underlying genetic structure. Out of total investigated buffaloes 73.5% had bilateral wall eyes while 5.4% had unilateral wall eyes and 21.1% had no wall eyes. 41.15% of Nili-Ravi buffaloes maintained in the breeding farm were having typical true to the type characteristics (both eyes walled, white markings in forehead, muzzle/chin, all the four legs and tail) while only 28.5% of Nili-Ravi buffaloes were true to the type under field conditions. Genotypic data were generated in four groups of Nili-Ravi buffalo (FMTNR--Typical Nili-Ravi from farm; FMANR--Atypical Nili-Ravi from farm; FDTNR--Typical Nili-Ravi from field; FDANR--Atypical Nili-Ravi from field) at 16 microsatellite loci. Comparative genetic analysis of various groups of Nili-Ravi buffaloes with Murrah revealed significant between group differences with an estimated global F(ST) of 0.063. Pair-wise F(ST) values ranged from 0.003 (between FDTNR and FDANR) to 0.112 (between FMTNR and FDTNR). Phylogenetic analysis and multi-dimensional scaling revealed clustering of FDTNR and FDANR together while FMTNR and FMANR clustered separately with Murrah in between farm and field Nili-Ravi buffaloes. Based on the results, the paper also proposes three pronged strategy for conservation and sustainable genetic improvement of Nili-Ravi buffalo in India.

Genetic Evaluation of Dual-Purpose Buffaloes (Bubalus bubalis) in Colombia Using Principal Component Analysis

Genealogy and productive information of 48621 dual-purpose buffaloes born in Colombia between years 1996 and 2014 was used. The following traits were assessed using one-trait models: milk yield at 270 days (MY270), age at first calving (AFC), weaning weight (WW), and weights at the following ages: first year (W12), 18 months (W18), and 2 years (W24). Direct additive genetic and residual random effects were included in all the traits. Maternal permanent environmental and maternal additive genetic effects were included for WW and W12. The fixed effects were: contemporary group (for all traits), sex (for WW, W12, W18, and W24), parity (for WW, W12, and MY270). Age was included as covariate for WW, W12, W18 and W24. Principal component analysis (PCA) was conducted using the genetic values of 133 breeding males whose breeding-value reliability was higher than 50% for all the traits in order to define the number of principal components (PC) which would explain most of the variation. The highest heritabilities were for W18 and MY270, and the lowest for AFC; with 0.53, 0.23, and 0.17, respectively. The first three PCs represented 66% of the total variance. Correlation of the first PC with meat production traits was higher than 0.73, and it was -0.38 with AFC. Correlations of the second PC with maternal genetic component traits for WW and W12 were above 0.75. The third PC had 0.84 correlation with MY270. PCA is an alternative approach for analyzing traits in dual-purpose buffaloes and reduces the dimension of the traits.

The Differential Expression of Immune Genes between Water Buffalo and Yellow Cattle Determines Species-Specific Susceptibility to Schistosoma japonicum Infection

Water buffalo are less susceptible to Schistosoma japonicum infection than yellow cattle. The factors that affect such differences in susceptibility remain unknown. A Bos taurus genome-wide gene chip was used to analyze gene expression profiles in the peripheral blood of water buffalo and yellow cattle pre- and post-infection with S. japonicum. This study showed that most of the identified differentially expressed genes (DEGs) between water buffalo and yellow cattle pre- and post-infection were involved in immune-related processes, and the expression level of immune genes was lower in water buffalo. The unique DEGs (390) in yellow cattle were mainly associated with inflammation pathways, while the unique DEGs (2,114) in water buffalo were mainly associated with immune-related factors. The 83 common DEGs may be the essential response genes during S. japonicum infection, the highest two gene ontology (GO) functions were associated with the regulation of fibrinolysis. The pathway enrichment analysis showed that the DEGs constituted similar immune-related pathways pre- and post-infection between the two hosts. This first analysis of the transcriptional profiles of natural hosts has enabled us to gain new insights into the mechanisms that govern their susceptibility or resistance to S. japonicum infections.

Molecular characterization and expression of buffalo (Bubalus bubalis) DEAD-box family VASA gene and mRNA transcript variants isolated from testis tissue

VASA is a member of the DEAD-box protein family that plays an indispensable role in mammalian spermatogenesis, particularly during meiosis. In the present study, we isolated, sequenced, and characterized VASA gene in buffalo testis. Here, we demonstrated that VASA mRNA is expressed as multiple isoforms and uses four alternative transcriptional start sites (TSSs) and four different polyadenylation sites. The TSSs identified by 5'-RNA ligase-mediated rapid amplification of cDNA ends (RLM-5'-RACE) were positioned at 48, 53, 85, and 88 nucleotides upstream relative to the translation initiation codon. 3'-RACE experiment revealed the presence of tandem polyadenylation signals, which lead to the expression of at least four different 3'-untranslated regions (209, 233, 239 and 605 nucleotides). The full-length coding region of VASA was 2190 bp, which encodes a 729 amino acid (aa) protein containing nine consensus regions of the DEAD box protein family. VASA variants are highly expressed in testis of adult buffalo. We found five variants, one full length VASA (729 aa) and four splice variants VASA 2, 4, 5, 6 (683, 685, 679, 703 aa). The expression level of VASA 1 was significantly higher than rest of all (P < 0.05) except VASA 6. The relative ratio for VASA 1:2:4:5:6 was 100:1.0:1.6:0.9:48.

Establishment of Trophectoderm Cell Lines from Buffalo (Bubalus bubalis) Embryos of Different Sources and Examination of In Vitro Developmental Competence, Quality, Epigenetic Status and Gene Expression in Cloned Embryos Derived from Them

Despite being successfully used to produce live offspring in many species, somatic cell nuclear transfer (NT) has had a limited applicability due to very low (>1%) live birth rate because of a high incidence of pregnancy failure, which is mainly due to placental dysfunction. Since this may be due to abnormalities in the trophectoderm (TE) cell lineage, TE cells can be a model to understand the placental growth disorders seen after NT. We isolated and characterized buffalo TE cells from blastocysts produced by in vitro fertilization (TE-IVF) and Hand-made cloning (TE-HMC), and compared their growth characteristics and gene expression, and developed a feeder-free culture system for their long-term culture. The TE-IVF cells were then used as donor cells to produce HMC embryos following which their developmental competence, quality, epigenetic status and gene expression were compared with those of HMC embryos produced using fetal or adult fibroblasts as donor cells. We found that although TE-HMC and TE-IVF cells have a similar capability to grow in culture, significant differences exist in gene expression levels between them and between IVF and HMC embryos from which they are derived, which may have a role in the placental abnormalities associated with NT pregnancies. Although TE cells can be used as donor cells for producing HMC blastocysts, their developmental competence and quality is lower than that of blastocysts produced from fetal or adult fibroblasts. The epigenetic status and expression level of many important genes is different in HMC blastocysts produced using TE cells or fetal or adult fibroblasts or those produced by IVF.

Splice variants and seasonal expression of buffalo HSF genes

In eukaryotes, the heat shock factors (HSFs) are recognized as the master regulator of the heat shock response. In this respect, the genes encoding the heat shock factors seem to be important for adaptation to thermal stress in organisms. Despite this, only few mammalian HSFs has been characterized. In this study, four major heat shock factor genes viz. HSF-1, 2, 4, and 5 were studied. The main objective of the present study was to characterize the cDNA encoding using conserved gene specific primers and to investigate the expression status of these buffalo HSF genes. Our RT-PCR analysis uncovered two distinct variants of buffalo HSF-1 and HSF-2 gene transcripts. In addition, we identified a variant of the HSF5 transcript in buffalo lacking a DNA-binding domain. In silico analysis of deduced amino acid sequences for buffalo HSF genes showed domain architecture similar to other mammalian species. Changes in the gene expression profile were noted by quantitative real-time PCR (qRT-PCR) analysis. We detected the transcript of buffalo HSF genes in different tissues. We also evaluated the seasonal changes in the expression of HSF genes. Interestingly, the transcript level of HSF-1 gene was found upregulated in months of high and low ambient temperatures. In contrast, the expression of the HSF-4 and 5 genes was found to be downregulated in months of high ambient temperature. This suggests that the intricate balance of different HSFs is adjusted to minimize the effect of seasonal changes in environmental conditions. These findings advance our understanding of the complex, context-dependent regulation of HSF gene expression under normal and stressful conditions.

Expression and localization of ghrelin and its receptor in ovarian follicles during different stages of development and the modulatory effect of ghrelin on granulosa cells function in buffalo

Ghrelin, a hormone predominantly found in the stomach, was recently described as a factor that controls female reproductive function. The aim of our study was to investigate the expression and localization of ghrelin and its active receptor, growth hormone secretagogue receptor type 1a (GHS-R1a) in buffalo ovarian follicles of different follicular size and to investigate role of ghrelin on estradiol (E2) secretion, aromatase (CYP19A1), proliferating cell nuclear antigen (PCNA) and apoptosis regulator Bax gene expression on granulosa cell culture. Using real time PCR and western blot, we measured gene and protein expression of examined factors. Localization was done with immunofluorescence method. Expression of ghrelin increased with follicle size with significantly highest in dominant or pre-ovulatory follicle (P<0.05). Expression of GHS-R1a was comparable in medium and large follicle but was higher than small follicles (P<0.05). Both the factors were localized in granulosa and theca cells. Pattern of intensity of immunofluorescence was similar with mRNA and protein expression. In the in vitro study granulosa cells (GCs) were cultured and treated with ghrelin each at 1, 10 and 100ng/ml concentrations for two days after obtaining 75-80 per cent confluence. Ghrelin treatment significantly (P<0.05) inhibited E2 secretion, CYP19A1 expression, apoptosis and promoted cell proliferation. In conclusion, this study provides novel evidence for the presence of ghrelin and receptor GHS-R1a in ovarian follilcles and modulatory role of ghrelin on granulosa cell function in buffalo.

Positive selection of deleterious alleles through interaction with a sex-ratio suppressor gene in African Buffalo: a plausible new mechanism for a high frequency anomaly

Although generally rare, deleterious alleles can become common through genetic drift, hitchhiking or reductions in selective constraints. Here we present a possible new mechanism that explains the attainment of high frequencies of deleterious alleles in the African buffalo (Syncerus caffer) population of Kruger National Park, through positive selection of these alleles that is ultimately driven by a sex-ratio suppressor. We have previously shown that one in four Kruger buffalo has a Y-chromosome profile that, despite being associated with low body condition, appears to impart a relative reproductive advantage, and which is stably maintained through a sex-ratio suppressor. Apparently, this sex-ratio suppressor prevents fertility reduction that generally accompanies sex-ratio distortion. We hypothesize that this body-condition-associated reproductive advantage increases the fitness of alleles that negatively affect male body condition, causing genome-wide positive selection of these alleles. To investigate this we genotyped 459 buffalo using 17 autosomal microsatellites. By correlating heterozygosity with body condition (heterozygosity-fitness correlations), we found that most microsatellites were associated with one of two gene types: one with elevated frequencies of deleterious alleles that have a negative effect on body condition, irrespective of sex; the other with elevated frequencies of sexually antagonistic alleles that are negative for male body condition but positive for female body condition. Positive selection and a direct association with a Y-chromosomal sex-ratio suppressor are indicated, respectively, by allele clines and by relatively high numbers of homozygous deleterious alleles among sex-ratio suppressor carriers. This study, which employs novel statistical techniques to analyse heterozygosity-fitness correlations, is the first to demonstrate the abundance of sexually-antagonistic genes in a natural mammal population. It also has important implications for our understanding not only of the evolutionary and ecological dynamics of sex-ratio distorters and suppressors, but also of the functioning of deleterious and sexually-antagonistic alleles, and their impact on population viability.

Genetic structure of fragmented southern populations of African Cape buffalo (Syncerus caffer caffer)

BACKGROUND

African wildlife experienced a reduction in population size and geographical distribution over the last millennium, particularly since the 19th century as a result of human demographic expansion, wildlife overexploitation, habitat degradation and cattle-borne diseases. In many areas, ungulate populations are now largely confined within a network of loosely connected protected areas. These metapopulations face gene flow restriction and run the risk of genetic diversity erosion. In this context, we assessed the "genetic health" of free ranging southern African Cape buffalo populations (S.c. caffer) and investigated the origins of their current genetic structure. The analyses were based on 264 samples from 6 southern African countries that were genotyped for 14 autosomal and 3 Y-chromosomal microsatellites.

RESULTS

The analyses differentiated three significant genetic clusters, hereafter referred to as Northern (N), Central (C) and Southern (S) clusters. The results suggest that splitting of the N and C clusters occurred around 6000 to 8400 years ago. Both N and C clusters displayed high genetic diversity (mean allelic richness (A r ) of 7.217, average genetic diversity over loci of 0.594, mean private alleles (P a ) of 11), low differentiation, and an absence of an inbreeding depression signal (mean F IS = 0.037). The third (S) cluster, a tiny population enclosed within a small isolated protected area, likely originated from a more recent isolation and experienced genetic drift (F IS = 0.062, mean A r = 6.160, P a = 2). This study also highlighted the impact of translocations between clusters on the genetic structure of several African buffalo populations. Lower differentiation estimates were observed between C and N sampling localities that experienced translocation over the last century.

CONCLUSIONS

We showed that the current genetic structure of southern African Cape buffalo populations results from both ancient and recent processes. The splitting time of N and C clusters suggests that the current pattern results from human-induced factors and/or from the aridification process that occurred during the Holocene period. The more recent S cluster genetic drift probably results of processes that occurred over the last centuries (habitat fragmentation, diseases). Management practices of African buffalo populations should consider the micro-evolutionary changes highlighted in the present study.

Developmental competence and expression pattern of bubaline (Bubalus bubalis) oocytes subjected to elevated temperatures during meiotic maturation in vitro

OBJECTIVE

To determine the direct effect of physiologically relevant high temperatures (40.5 and 41.5 °C) for two time periods (12 and 24 h) on bubaline oocytes during in vitro maturation.

METHOD

The control group oocytes were cultured at 38.5 °C for 24 h. The treatment 1 (T1) and 3 (T3) group oocytes were cultured at 40.5 and 41.5 °C respectively, for the first 12 h and at 38.5 °C for rest of the 12 h. However, treatment 2 (T2) and 4 (T4) group oocytes were cultured at 40.5 and 41.5 °C for complete 24 h.

RESULTS

Development of oocytes to blastocyst was severely compromised (p < 0.001) when matured at 40.5 and 41.5 °C for both exposure periods (12 h and 24 h). It was found that the cleavage rates, blastocyst yield and mean cell number decreased remarkably (p < 0.001) in the treatment groups compared to control. The relative mRNA expression of heat shock protein (Hsp 70.1, 70.2, 70.8, 60, 10 and HSF1), pro-apoptotic (caspases-3, -7, -8, Bid and Bax) and oxidative stress (iNOS) related genes was significantly higher (p < 0.05) in all the treatment groups compared to control. However, mRNA abundance of anti-apoptotic (Bcl-2, Mcl-1, Bcl-xl), glucose transport (Glut1, Glut3 and IGF1R), developmental competence (ZAR1 and BMP15) and oxidative stress (MnSOD) related genes was significantly decreased (p < 0.05) in the treatment groups compared to control.

CONCLUSION

The present study clearly establishes that physiologically relevant elevated temperatures during in vitro meiotic maturation reduce developmental competence of bubaline oocytes.

Global transcriptome profiles of Italian Mediterranean buffalo embryos with normal and retarded growth

The transcriptome profiles were compared for buffalo embryos with normal growth and embryos with retarded growth on Day 25 after mating. Embryos with retarded growth on Day 25 after mating have a reduced likelihood of undergoing attachment to the uterine endometrium and establishing a pregnancy. Italian Mediterranean buffaloes were mated by AI and on Day 25 underwent trans-rectal ultrasonography to ascertain embryo development. Embryos with an embryonic width (EW)>2.7 mm were classed as normal embryos and embryos with an EW<2.7 mm were classed as retarded embryos. Three buffaloes with embryos of the largest EW (3.7, 3.7 and 3.9 mm) and three buffaloes with embryos of the smallest EW (1.5, 1.6 and 1.9 mm) were slaughtered on Day 27 to recover embryos for transcriptome analysis using a bovine custom designed oligo array. A total of 1,047 transcripts were differentially expressed between embryos with normal growth and embryos with retarded growth. Retarded embryos showed 773/1,047 (74%) transcripts that were down-regulated and 274/1,047 (26%) transcripts that were up-regulated relative to normal embryos; in silico analyses focused on 680/1,047 (65%) of the differentially expressed transcripts. The most altered transcripts observed in retarded embryos were associated with membrane structure and function and with metabolic and homeostasis maintenance functions. Other notable functions altered in retarded embryos were developmental processes and in particular nervous system differentiation and function. Specific biochemical pathways such as the complement cascade and coagulation were also altered in retarded embryos. It was concluded from the findings that buffalo embryos with retarded growth on Day 25 after mating show altered gene expression compared with normal embryos, and some de-regulated functions are associated with attachment to the uterine endometrium.

Transcriptional profile of endometrial TLR4 and 5 genes during the estrous cycle and uterine infection in the buffalo (Bubalus bubalis)

Endometritis is one of the leading causes of infertility in the cattle and buffalo and innate immune mechanism plays an important role in clearing the infection. In this regard, endometrial expression and function of Toll Like Receptors (TLR) are focus of investigation in the recent years. In this study, we report the transcriptional profiles of TLR4 and 5 in the buffalo endometrium during the follicular, early, mid and late luteal phases of estrous cycle and 'subclinical and clinical endometritis' and also at true anestrus (n = 10 for each stage) using RT-PCR and qRT-PCR as they are the ligands for the lipopolysaccharide and flagellin components of E.coli, the most common cause of postpartum endometritis. We found a significant positive correlation between TLR4 and 5 in all the groups (r = 0.696-0.803; P < 0.05) except late luteal phase (r = 0.522; P > 0.05). Chi-square analysis showed that the qualitative expression of endometrial TLR4 and 5 transcripts was significantly associated with the phase of estrous cycle and also with uterine infection (P < 0.05). Further, using true anestrus category as a calibrator group, relative quantitation of TLR4 and 5 revealed that the transcriptional expression of TLR4 and 5 genes were highly upregulated (24.6-83.3 folds) during endometritis conditions and moderately upregulated during mid-luteal phase (6.8-16.2) of the estrous cycle (P < 0.05). The results suggested a role of progesterone in the expression of TLR4 and 5.

Molecular characterization and expression profile of ghrelin gene during different reproductive phases in buffalo (Bubalus bubalis)

Ghrelin, a novel motilin-related endogenous ligand for growth hormone secretagouge receptor, is implicated in various biological functions, including regulation of female reproduction. But the presence of ghrelin and its role in reproductive functions in buffalo, a species with poor reproductive efficiency, is not known. In the present study full-length ghrelin cDNA was isolated from bubaline abomasum, which encodes the entire prepropeptide of 116 amino acids. The deduced amino acid sequence of ghrelin of buffalo showed >95% and 31% identity with that of ruminants (cattle, sheep, and goat) and humans, respectively. Analysis of synonymous and nonsynonymous nucleotide substitutions in the coding region of ghrelin indicated that these sequences of different species have been under purifying selection. The 3995-bp amplicon of ghrelin gene consisting of 4 exons and 3 introns was cloned with genomic DNA from buffalo. Further, ghrelin expression was determined by quantitative real-time PCR, in situ hybridization, and immunohistochemistry in bubaline endometrial tissues at different stages of the estrous cycle and early pregnancy. Our results indicated the persistent expression of ghrelin mRNA and protein in the endometrium during stage I (day 3-5), stage II (day 6-15), and stage III (day 16-21) of the estrous cycle and also during early (~day 30-40) pregnancy. Immunohistochemistry and quantitative real-time PCR experiments indicated the relatively higher expression of ghrelin in the endometrium during stage II (day 6-15) of the estrous cycle and early pregnancy than during stage I (day 3-5) and stage III (day 16-21) of the estrous cycle, but no statistically significant difference in ghrelin expression was observed among stages. To conclude, the results of the present study indicate the persistent expression of ghrelin in the uterine endometrium throughout the estrous cycle and in early pregnancy which might be helpful in determining its role in buffalo reproduction.

Novel polymorphisms of the IGF1R gene and their association with average daily gain in Egyptian buffalo (Bubalus bubalis)

The objective of this study was to detect insulin-like growth factor 1 receptor (IGF1R) polymorphisms, their allele, and genotype frequencies and to determine associations between these polymorphisms and growth traits in Egyptian water buffalo. Three loci of the IGF1R coding region were amplified by RT-PCR and, subsequently, subjected to sequence analysis, followed by single-strand conformation polymorphism to identify different allelic patterns. A total of 11 novel polymorphisms were detected; 6 SNPs among Egyptian water buffaloes and 5 polymorphisms compared with Indian buffalo (Y12700). Three of those polymorphisms; GAG Indel polymorphism, C261G, and G263C SNPs, were nonsynonymous mutations. The GAG Indel polymorphism led to deletion of E (glutamic) amino acid (aa) in the IGF1R of Egyptian water buffaloes compared with Indian buffalo. However, C261G SNP, which replaced A (alanine) by G (glycine) aa, and G263C SNP, which changed A (alanine) to P (proline) aa, were detected among Egyptian water buffaloes. Three different single-strand conformation polymorphism patterns were observed in exon 21: CC/CC, GG/GG, and CG/GC with frequencies of 0.291, 0.253, and 0.556, respectively. The heterozygous animals (CG/GC) had a higher ADG than homozygous animals (CC/CC and GG/GG) from birth to 6 mo of age. We conclude that the heterozygous haplotype, C261G/G263C, in exon 21 of the IGF1R gene is associated with the ADG during the early stages of life (from birth to 6 mo of age) and could be used as a genetic marker for selection of growth traits in Egyptian buffalo.

Identification of three novel myeloid cathelicidin cDNAs and their predicted peptides in buffalo (Bubalus bubalis)

Antimicrobial peptides (AMPs) are broad spectrum antibiotics, which mostly act without specific receptors. Identification of AMPs is important in the current scenario of emerging multi-drug resistant bacteria. In the present study, in an attempt to identify new AMPs, myeloid cathelicidin cDNAs were synthesized from buffalo (Bubalus bubalis) bone marrow and were amplified using specific primers. Sequence analysis of cloned cDNAs revealed three novel myeloid cathelicidins. They were named based on the number of active amino acids in the C-terminal region of their predicted peptide sequences as BuMAP-28 (having an additional Gly at position 22nd), BuMAP-29 (having an additional IIe at position 27) and BuMAP-34, compared to BMAP-27, BMAP-28 and BMAP-34 of cattle. The BuMAPs showed 93%, 95% and 87% homology respectively with that of its cattle counterpart. Predicted number of amino acids of the cDNAs was 159, 155 and 157 residues, with cationic C-terminal sequences of 28, 29 and 34, respectively, which correspond to putative antimicrobial domains. Several amino acid substitutions were observed in all the three cathelicidins. The conformation of the peptides was predicted to be alpha helical, having total net positive charge and hydrophobicity, similar to that of BMAPs in cattle. Comparative analysis of the predicted peptides suggested potential antimicrobial activity and the sequence variations detected might enable the peptides to act as effective broad spectrum antimicrobial agents.

Novel polymorphisms in UTR and coding region of inducible heat shock protein 70.1 gene in tropically adapted Indian zebu cattle (Bos indicus) and riverine buffalo (Bubalus bubalis)

Due to evolutionary divergence, cattle (taurine, and indicine) and buffalo are speculated to have different responses to heat stress condition. Variation in candidate genes associated with a heat-shock response may provide an insight into the dissimilarity and suggest targets for intervention. The present work was undertaken to characterize one of the inducible heat shock protein genes promoter and coding regions in diverse breeds of Indian zebu cattle and buffaloes. The genomic DNA from a panel of 117 unrelated animals representing 14 diversified native cattle breeds and 6 buffalo breeds were utilized to determine the complete sequence and gene diversity of HSP70.1 gene. The coding region of HSP70.1 gene in Indian zebu cattle, Bos taurus and buffalo was similar in length (1,926 bp) encoding a HSP70 protein of 641 amino acids with a calculated molecular weight (Mw) of 70.26 kDa. However buffalo had a longer 5' and 3' untranslated region (UTR) of 204 and 293 nucleotides respectively, in comparison to Indian zebu cattle and Bos taurus wherein length of 5' and 3'-UTR was 172 and 286 nucleotides, respectively. The increased length of buffalo HSP70.1 gene compared to indicine and taurine gene was due to two insertions each in 5' and 3'-UTR. Comparative sequence analysis of cattle (taurine and indicine) and buffalo HSP70.1 gene revealed a total of 54 gene variations (50 SNPs and 4 INDELs) among the three species in the HSP70.1 gene. The minor allele frequencies of these nucleotide variations varied from 0.03 to 0.5 with an average of 0.26. Among the 14 B. indicus cattle breeds studied, a total of 19 polymorphic sites were identified: 4 in the 5'-UTR and 15 in the coding region (of these 2 were non-synonymous). Analysis among buffalo breeds revealed 15 SNPs throughout the gene: 6 at the 5' flanking region and 9 in the coding region. In bubaline 5'-UTR, 2 additional putative transcription factor binding sites (Elk-1 and C-Re1) were identified, other than three common sites (CP2, HSE and Pax-4) observed across all the analyzed animals. No polymorphism was found within the 3'-UTR of Indian cattle or buffalo as it was found to be monomorphic. The promoter sequences generated in 117 individuals showed a rich array of sequence elements known to be involved in transcription regulation. A total of 11 nucleotide changes were observed in the promoter sequence across the analyzed species, 3 of these changes were located within the potential transcription factor binding domains. We also identified 4 microsatellite markers within the buffalo HSP70.1 gene and 3 microsatellites within bovine HSP70.1. The present study identified several distinct changes across indicine, taurine and bubaline HSP70.1 genes that could further be evaluated as molecular markers for thermotolerance.

Gene polymorphisms in African buffalo associated with susceptibility to bovine tuberculosis infection

Bovine tuberculosis (BTB) is a chronic, highly infectious disease that affects humans, cattle and numerous species of wildlife. In developing countries such as South Africa, the existence of extensive wildlife-human-livestock interfaces poses a significant risk of Mycobacterium bovis transmission between these groups, and has far-reaching ecological, economic and public health impacts. The African buffalo (Syncerus caffer), acts as a maintenance host for Mycobacterium bovis, and maintains and transmits the disease within the buffalo and to other species. In this study we aimed to investigate genetic susceptibility of buffalo for Mycobacterium bovis infection. Samples from 868 African buffalo of the Cape buffalo subspecies were used in this study. SNPs (n = 69), with predicted functional consequences in genes related to the immune system, were genotyped in this buffalo population by competitive allele-specific SNP genotyping. Case-control association testing and statistical analyses identified three SNPs associated with BTB status in buffalo. These SNPs, SNP41, SNP137 and SNP144, are located in the SLC7A13, DMBT1 and IL1α genes, respectively. SNP137 remained significantly associated after permutation testing. The three genetic polymorphisms identified are located in promising candidate genes for further exploration into genetic susceptibility to BTB in buffalo and other bovids, such as the domestic cow. These polymorphisms/genes may also hold potential for marker-assisted breeding programmes, with the aim of breeding more BTB-resistant animals and herds within both the national parks and the private sector.

Effects of a single nucleotide polymorphism in the leptin gene on the productive traits of dairy buffaloes (Bubalus bubalis)

The gene responsible for coding the leptin hormone has been associated with productive and reproductive traits in cattle. In dairy cattle, different polymorphisms found in the leptin gene have been associated with several traits of economic interest, such as energy balance, milk yield and composition, live weight, fertility and dry matter consumption. The aim of this study was to detect genetic variability in the leptin gene of buffaloes and to test possible associations with milk yield, fat and protein percentages, age at first calving and first calving interval. Three genotypes (AA, AG and GG) were identified by polymerase chain reaction-restriction fragment length polymorphism, which presented genotypic frequencies of 0.30, 0.54 and 0.16, respectively. The allele frequencies were 0.57 for the A allele and 0.43 for the G allele. No significant effects were found in the present study, but there is an indicative that leptin gene affects lipid metabolism.

The confounding effect of population structure on Bayesian skyline plot inferences of demographic history

Many coalescent-based methods aiming to infer the demographic history of populations assume a single, isolated and panmictic population (i.e. a Wright-Fisher model). While this assumption may be reasonable under many conditions, several recent studies have shown that the results can be misleading when it is violated. Among the most widely applied demographic inference methods are Bayesian skyline plots (BSPs), which are used across a range of biological fields. Violations of the panmixia assumption are to be expected in many biological systems, but the consequences for skyline plot inferences have so far not been addressed and quantified. We simulated DNA sequence data under a variety of scenarios involving structured populations with variable levels of gene flow and analysed them using BSPs as implemented in the software package BEAST. Results revealed that BSPs can show false signals of population decline under biologically plausible combinations of population structure and sampling strategy, suggesting that the interpretation of several previous studies may need to be re-evaluated. We found that a balanced sampling strategy whereby samples are distributed on several populations provides the best scheme for inferring demographic change over a typical time scale. Analyses of data from a structured African buffalo population demonstrate how BSP results can be strengthened by simulations. We recommend that sample selection should be carefully considered in relation to population structure previous to BSP analyses, and that alternative scenarios should be evaluated when interpreting signals of population size change.

Phylogeography and domestication of Chinese swamp buffalo

To further probe into whether swamp buffaloes were domesticated once or multiple times in China, this survey examined the mitochondrial DNA (mtDNA) Control Region (D-loop) diversity of 471 individuals representing 22 populations of 455 Chinese swamp buffaloes and 16 river buffaloes. Phylogenetic analysis revealed that Chinese swamp buffaloes could be divided into two distinct lineages, A and B, which were defined previously. Of the two lineages, lineage A was predominant across all populations. For predominant lineage A, Southwestern buffalo populations possess the highest genetic diversity among the three hypothesized domestication centers (Southeastern, Central, and Southwestern China), suggesting Southwestern China as the most likely location for the domestication of lineage A. However, a complex pattern of diversity is detected for the lineage B, preventing the unambiguous pinpointing of the exact place of domestication center and suggesting the presence of a long-term, strong gene flow among swamp buffalo populations caused by extensive migrations of buffaloes and frequent human movements along the Yangtze River throughout history. Our current study suggests that Southwestern China is the most likely domestication center for lineage A, and may have been a primary center of swamp buffalo domestication. More archaeological and genetic evidence is needed to show the process of domestication.

Novel SNP Discovery in African Buffalo, Syncerus caffer, using high-throughput Sequencing

The African buffalo, Syncerus caffer, is one of the most abundant and ecologically important species of megafauna in the savannah ecosystem. It is an important prey species, as well as a host for a vast array of nematodes, pathogens and infectious diseases, such as bovine tuberculosis and corridor disease. Large-scale SNP discovery in this species would greatly facilitate further research into the area of host genetics and disease susceptibility, as well as provide a wealth of sequence information for other conservation and genomics studies. We sequenced pools of Cape buffalo DNA from a total of 9 animals, on an ABI SOLiD4 sequencer. The resulting short reads were mapped to the UMD3.1 Bos taurus genome assembly using both BWA and Bowtie software packages. A mean depth of 2.7× coverage over the mapped regions was obtained. Btau4 gene annotation was added to all SNPs identified within gene regions. Bowtie and BWA identified a maximum of 2,222,665 and 276,847 SNPs within the buffalo respectively, depending on analysis method. A panel of 173 SNPs was validated by fluorescent genotyping in 87 individuals. 27 SNPs failed to amplify, and of the remaining 146 SNPs, 43-54% of the Bowtie SNPs and 57-58% of the BWA SNPs were confirmed as polymorphic. dN/dS ratios found no evidence of positive selection, and although there were genes that appeared to be under negative selection, these were more likely to be slowly evolving house-keeping genes.

Comparative analysis of the Shadoo gene between cattle and buffalo reveals significant differences

Background

While prions play a central role in the pathogenesis of transmissible spongiform encephalopathies, the biology of these proteins and the pathophysiology of these diseases remain largely unknown. Since no case of bovine spongiform encephalopathy (BSE) has ever been reported in buffalo despite their phylogenetic proximity to cattle, genetic differences may be driving the different susceptibilities of these two species to BSE. We thus hypothesized that differences in expression of the most recently identified member of the prion family or Shadoo (SPRN) gene may relate to these species-specific differences.

Principal Findings

We first analyzed and compared the polymorphisms of the SPRN gene (∼4.4 kb), including the putative promoter, coding and 3′ regions, and further verified the entire ORF and putative promoter. This yielded a total of 117 fixed differences, remarkably: 1) a 12-bp insertion/deletion polymorphism in the hydrophobic domain of the cattle but not buffalo gene, introducing a four amino acid expansion/contraction in a series of 5 tandem Ala/Gly-containing repeats; 2) two fixed missense mutations (102Ser→Gly and 119Thr→Ala), and three missense mutations (92Pro>Thr/Met, 122Thr>Ile and 139Arg>Trp) in the coding region presenting different (P<0.05) genotypic and allelic frequency distributions between cattle and buffalo; and, 3) functional luciferase-reporter experiments for the predicted promoter region, consistent with a significantly higher activity in buffalo than cattle. Supporting these findings, immunoblotting revealed higher relative expression levels of Sho protein in cerebrum from buffalo than from cattle. In addition, for cattle, highest Sho expression was detected in obex, as compared to cerebrum or cerebellum.

Significance

Our findings support Sho as a non-PrP specific marker for prion infections, with obex as the best tissue source for the detection of Sho in TSE rapid tests. Moreover, these discoveries may prove advantageous for further understanding the biology of prion diseases.

Analysis of rumen methanogen diversity in water buffaloes (Bubalus bubalis) under three different diets

The water buffalo (Bubalus bubalis) is a prominent livestock species for the production of milk and meat in many countries. We investigated the diversity of rumen methanogens in Mediterranean water buffaloes maintained in Brazil under different diets: corn silage, grazing pasture, or sugar cane. A total of 467 clones were isolated from three methanogen 16S rRNA gene clone libraries that each represented a distinct feed type. The 467 clones were assigned to 19 species-level operational taxonomic units (OTUs). Four OTUs were represented in all three libraries, eight OTUs were library-specific, six OTUs were found in only the corn silage and pasture grazing libraries, and one OTU was shared only between pasture grazing and sugar cane libraries. We found that Methanobrevibacter-related sequences were the most abundant in the water buffaloes sampled for our analysis, in contrast to previously reported studies showing that Methanomicrobium mobile-like methanogens were the most abundant methanogens in water buffaloes of Murrah and Surti breeds sampled in India. Considering the worldwide distribution of water buffaloes and the likely wide variety of diets provided, our results combined with studies from other groups support that larger scope analyses of microbiomes for this livestock species would provide great insight into the contribution of geographical location, breed, and diet in determining the population structure of rumen microorganisms.