An assessment of population structure in eight breeds of cattle using a whole genome SNP panel

Background: Analyses of population structure and breed diversity have provided insight into the origin and evolution of cattle. Previously, these studies have used a low density of microsatellite markers, however, with the large number of single nucleotide polymorphism markers that are now available, it is possible to perform genome wide population genetic analyses in cattle. In this study, we used a high-density panel of SNP markers to examine population structure and diversity among eight cattle breeds sampled from Bos indicus and Bos taurus. Results: Two thousand six hundred and forty one single nucleotide polymorphisms ( SNPs) spanning all of the bovine autosomal genome were genotyped in Angus, Brahman, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black, Limousin and Nelore cattle. Population structure was examined using the linkage model in the program STRUCTURE and Fst estimates were used to construct a neighbor-joining tree to represent the phylogenetic relationship among these breeds. Conclusion: The whole-genome SNP panel identified several levels of population substructure in the set of examined cattle breeds. The greatest level of genetic differentiation was detected between the Bos taurus and Bos indicus breeds. When the Bos indicus breeds were excluded from the analysis, genetic differences among beef versus dairy and European versus Asian breeds were detected among the Bos taurus breeds. Exploration of the number of SNP loci required to differentiate between breeds showed that for 100 SNP loci, individuals could only be correctly clustered into breeds 50% of the time, thus a large number of SNP markers are required to replace the 30 microsatellite markers that are currently commonly used in genetic diversity studies.

European 1: A globally important clonal complex of Mycobacterium bovis

We have identified a globally important clonal complex of Mycobacterium bovis by deletion analysis of over one thousand strains from over 30 countries. We initially show that over 99% of the strains of M. bovis, the cause of bovine tuberculosis, isolated from cattle in the Republic of Ireland and the UK are closely related and are members of a single clonal complex marked by the deletion of chromosomal region RDEu1 and we named this clonal complex European 1 (Eu1). Eu1 strains were present at less than 14% of French, Portuguese and Spanish isolates of M. bovis but are rare in other mainland European countries and Iran. However, strains of the Eu1 clonal complex were found at high frequency in former trading partners of the UK (USA, South Africa, New Zealand, Australia and Canada). The Americas, with the exception of Brazil, are dominated by the Eu1 clonal complex which was at high frequency in Argentina, Chile, Ecuador and Mexico as well as North America. Eu1 was rare or absent in the African countries surveyed except South Africa. A small sample of strains from Taiwan were non-Eu1 but, surprisingly, isolates from Korea and Kazakhstan were members of the Eu1 clonal complex. The simplest explanation for much of the current distribution of the Eu1 clonal complex is that it was spread in infected cattle, such as Herefords, from the UK to former trading partners, although there is evidence of secondary dispersion since. This is the first identification of a globally dispersed clonal complex M. bovis and indicates that much of the current global distribution of this important veterinary pathogen has resulted from relatively recent International trade in cattle. Crown Copyright (C) 2011 Published by Elsevier B.V. All rights reserved.

Antimicrobial susceptibility of coagulase-negative Staphylococci isolated from mastitic cattle in Brazil

Um total de 109 cepas de Staphylococci coagulase-negativa foi isolado de leite de vacas com mastite clínica e subclínica, em 35 fazendas, situadas em nove estados brasileiros, no período de fevereiro a maio de 2005. Os isolados foram investigados em relação a susceptibilidade in vitro a diversos agentes antimicrobianos. A resistência à penicilina foi a observação mais freqüente (93,5%), seguida por sulfonamida (88,9%), novobiocina (88,6%) e ampicilina (85,3%). Todas as cepas examinadas mostraram resistência a pelo menos uma das drogas antimicrobianas testadas. Cepas apresentando resistência múltipla foram extremamente comuns, com 10,0% dos microrganismos isolados apresentando resistência a todas as drogas antimicrobianas. Os resultados obtidos indicaram que as cepas de Staphylococci coagulase-negativas, isoladas no Brasil, apresentaram um alto grau de resistência a antimicrobianos. Estes resultados são, provavelmente, uma conseqüência da pressão devida ao uso intensivo de drogas antimicrobianas.

Effects of dietary energy and nitrogen supplements on rumen fermentation and protozoa population in buffalo and zebu cattle

The effects were assessed of two energy sources in concentrate (ground grain corn vs. citrus pulp) and two nitrogen sources (soybean meal vs. urea) on rumen metabolism in four buffaloes and four zebu cattle (Nellore) with rumen cannula and fed in a 4 × 4 Latin square design with feeds containing 60% sugar cane. Energy supplements had no effect on the rumen ammonia concentration in cattle, but ground grain corn promoted higher ammonia level than citrus pulp in buffalo. Urea produced higher ammonia level than soybean meal in both animal species. On average, the buffaloes maintained a lower rumen ammonia concentration (11.7 mg/dL) than the cattle (14.5 mg/dL). Buffaloes had lower production of acetic acid than cattle (58.7 vs. 61.6 mol/100 mol) and higher of propionic acid (27.4 vs. 23.6 mol/100 mol). There was no difference in the butyric acid production between the buffaloes (13.6 mol/100 mol) and cattle (14.8 mol/100 mol) and neither in the total volatile fatty acids concentration (82.5 vs. 83.6 mM, respectively). The energy or nitrogen sources had no effect on rumen protozoa count in either animal species. The zebu cattle had higher rumen protozoa population (8.8 × 10(5)/mL) than the buffaloes (6.1 × 10(5)/mL). The rumen protozoa population differed between the animal species, except for Dasytricha and Charonina. The buffaloes had a lower Entodinium population than the cattle (61.0 vs 84.9%, respectively) and a greater percentage of species belonging to the Diplodiniinae subfamily than the cattle (28.6 vs. 1.4%, respectively). In cattle, ground corn is a better energy source than citrus pulp for use by Entodinium and Diplodiniinae. In the buffaloes, the Entodinium are favored by urea and Diplodiniinae species by soybean meal.

Association of polymorphisms at the ADIPOR1 regulatory region with type 2 diabetes and body mass index in a Brazilian population with European or African ancestry

Association studies between ADIPOR1 genetic variants and predisposition to type 2 diabetes (DM2) have provided contradictory results. We determined if two single nucleotide polymorphisms (SNP c.-8503G>A and SNP c.10225C>G) in regulatory regions of ADIPOR1 in 567 Brazilian individuals of European (EA; N = 443) or African (AfA; N = 124) ancestry from rural (quilombo remnants; N = 439) and urban (N = 567) areas. We detected a significant effect of ethnicity on the distribution of the allelic frequencies of both SNPs in these populations (EA: -8503A = 0.27; AfA: -8503A = 0.16; P = 0.001 and EA: 10225G = 0.35; AfA: 10225G = 0.51; P < 0.001). Neither of the polymorphisms were associated with DM2 in the case-control study in EA (SNP c.-8503G>A: DM2 group -8503A = 0.26; control group -8503A = 0.30; P = 0.14/SNP 10225C>G: DM2 group 10225G = 0.37; control group 10225G = 0.32; P = 0.40) and AfA populations (SNP c.-8503G>A: DM2 group -8503A = 0.16; control group -8503A = 0.15; P = 0.34/SNP 10225C>G: DM2 group 10225G = 0.51; control group 10225G = 0.52; P = 0.50). Similarly, none of the polymorphisms were associated with metabolic/anthropometric risk factors for DM2 in any of the three populations, except for HDL cholesterol, which was significantly higher in AfA heterozygotes (GC = 53.75 ± 17.26 mg/dL) than in homozygotes. We conclude that ADIPOR1 polymorphisms are unlikely to be major risk factors for DM2 or for metabolic/anthropometric measurements that represent risk factors for DM2 in populations of European and African ancestries.

Genome wide scan for quantitative trait loci affecting tick resistance in cattle (Bos taurus x Bos indicus)

Background : In tropical countries, losses caused by bovine tick Rhipicephalus (Boophilus) microplus infestation have a tremendous economic impact on cattle production systems. Genetic variation between Bos taurus and Bos indicus to tick resistance and molecular biology tools might allow for the identification of molecular markers linked to resistance traits that could be used as an auxiliary tool in selection programs. The objective of this work was to identify QTL associated with tick resistance/susceptibility in a bovine F2 population derived from the Gyr (Bos indicus) x Holstein (Bos taurus) cross. Results: Through a whole genome scan with microsatellite markers, we were able to map six genomic regions associated with bovine tick resistance. For most QTL, we have found that depending on the tick evaluation season (dry and rainy) different sets of genes could be involved in the resistance mechanism. We identified dry season specific QTL on BTA 2 and 10, rainy season specific QTL on BTA 5, 11 and 27. We also found a highly significant genome wide QTL for both dry and rainy seasons in the central region of BTA 23. Conclusions: The experimental F2 population derived from Gyr x Holstein cross successfully allowed the identification of six highly significant QTL associated with tick resistance in cattle. QTL located on BTA 23 might be related with the bovine histocompatibility complex. Further investigation of these QTL will help to isolate candidate genes involved with tick resistance in cattle.

Multi-trait and random regression mature weight heritability and breeding value estimates in Nelore cattle

Mature weight breeding values were estimated using a multi-trait animal model (MM) and a random regression animal model (RRM). Data consisted of 82 064 weight records from 8 145 animals, recorded from birth to eight years of age. Weights at standard ages were considered in the MM. All models included contemporary groups as fixed effects, and age of dam (linear and quadratic effects) and animal age as covariates. In the RRM, mean trends were modelled through a cubic regression on orthogonal polynomials of animal age and genetic maternal and direct and maternal permanent environmental effects were also included as random. Legendre polynomials of orders 4, 3, 6 and 3 were used for animal and maternal genetic and permanent environmental effects, respectively, considering five classes of residual variances. Mature weight (five years) direct heritability estimates were 0.35 (MM) and 0.38 (RRM). Rank correlation between sires' breeding values estimated by MM and RRM was 0.82. However, selecting the top 2% (12) or 10% (62) of the young sires based on the MM predicted breeding values, respectively 71% and 80% of the same sires would be selected if RRM estimates were used instead. The RRM modelled the changes in the (co) variances with age adequately and larger breeding value accuracies can be expected using this model.

Predictive values of aspartate aminotransferase and gamma-glutamyl transferase for the hepatic accumulation of copper in cattle and buffalo

Ten cattle and 10 buffalo were divided into 2 groups (control [n = 8] and experimental [n = 12]) that received daily administration of copper. Three hepatic biopsies and blood samples were performed on days 0, 45, and 105. The concentration of hepatic copper was determined by spectrophotometric atomic absorption, and the activities of aspartate aminotransferase (AST) and gamma-glutamyl transferase (GGT) were analyzed. Regression analyses were done to verify the possible existing relationship between enzymatic activity and concentration of hepatic copper. Sensitivity, specificity, accuracy, and positive and negative predictive values were determined. The serum activities of AST and GGT had coefficients of determination that were excellent predictive indicators of hepatic copper accumulation in cattle, while only GGT serum activity was predictive of hepatic copper accumulation in buffalo. Elevated serum GGT activity may be indicative of increased concentrations of hepatic copper even in cattle and buffalo that appear to be clinically healthy. Thus, prophylactic measures can be implemented to prevent the onset of a hemolytic crisis that is characteristic of copper intoxication.

Allopatric speciation in ticks: genetic and reproductive divergence between geographic strains of Rhipicephalus (Boophilus) microplus

Background: The cattle tick, Rhipicephalus (Boophilus) microplus, economically impact cattle industry in tropical and subtropical regions of the world. The morphological and genetic differences among R. microplus strains have been documented in the literature, suggesting that biogeographical and ecological separation may have resulted in boophilid ticks from America/Africa and those from Australia being different species. To test the hypothesis of the presence of different boophilid species, herein we performed a series of experiments to characterize the reproductive performance of crosses between R. microplus from Australia, Africa and America and the genetic diversity of strains from Australia, Asia, Africa and America. Results: The results showed that the crosses between Australian and Argentinean or Mozambican strains of boophilid ticks are infertile while crosses between Argentinean and Mozambican strains are fertile. These results showed that tick strains from Africa (Mozambique) and America (Argentina) are the same species, while ticks from Australia may actually represent a separate species. The genetic analysis of mitochondrial 12S and 16S rDNA and microsatellite loci were not conclusive when taken separately, but provided evidence that Australian tick strains were genetically different from Asian, African and American strains. Conclusion: The results reported herein support the hypothesis that at least two different species share the name R. microplus. These species could be redefined as R. microplus (Canestrini, 1887) (for American and African strains) and probably the old R. australis Fuller, 1899 (for Australian strains), which needs to be redescribed. However, experiments with a larger number of tick strains from different geographic locations are needed to corroborate these results.

Molecular and geographic analyses of vampire bat-transmitted cattle rabies in central Brazil

Background: Vampire bats are important rabies virus vectors, causing critical problems in both the livestock industry and public health sector in Latin America. In order to assess the epidemiological characteristics of vampire bat-transmitted rabies, the authors conducted phylogenetic and geographical analyses using sequence data of a large number of cattle rabies isolates collected from a wide geographical area in Brazil. Methods: Partial nucleoprotein genes of rabies viruses isolated from 666 cattle and 18 vampire bats between 1987 and 2006 were sequenced and used for phylogenetic analysis. The genetic variants were plotted on topographical maps of Brazil. Results: In this study, 593 samples consisting of 24 genetic variants were analyzed. Regional localization of variants was observed, with the distribution of several variants found to be delimited by mountain ranges which served as geographic boundaries. The geographical distributions of vampire-bat and cattle isolates that were classified as the identical phylogenetic group were found to overlap with high certainty. Most of the samples analyzed in this study were isolated from adjacent areas linked by rivers. Conclusion: This study revealed the existence of several dozen regional variants associated with vampire bats in Brazil, with the distribution patterns of these variants found to be affected by mountain ranges and rivers. These results suggest that epidemiological characteristics of vampire bat-related rabies appear to be associated with the topographical and geographical characteristics of areas where cattle are maintained, and the factors affecting vampire bat ecology.

Characterization of mitochondrial genotypes in the foundation herd of the Canchim beef cattle breed

The Canchim (5/8 Charolais + 3/8 Zebu) beef cattle breed was developed at Southeast-Embrapa Cattle to take advantage of hybrid vigor and to combine the higher growth rate and beef quality of Charolais with tropical adaptations of Zebu. The development of three lineages (old, new, and crossbred) has increased its genetic basis. The genotypic origin (Bos taurus or Bos indicus) of the mitochondrial DNA (mtDNA) of the Canchim breed was unknown. We characterized the mtDNA genotype of this founder herd by allele-specific polymerase chain reaction. The 173 founder Zebu females (62 Indubrasil, 3 Guzerat, and 108 Nellore) and their 6749 offspring were identified. The frequency of B. indicus mtDNA ranged from 1.15 to 2.05% among the descendants (N = 6404) of each maternal line with available DNA, and among animals that were alive (N = 689) in December 2007 among the three lineages. Though mtDNA characterization can be used to direct animal selection, the low frequency of B. indicus mtDNA impairs the evaluation of its effects on production traits in these animals. The high prevalence of B. taurus mtDNA in Canchim proves that the founder Zebu females from the Indubrasil, Guzerat and Nellore breeds were obtained from crosses of Zebu sires with local B. taurus dams.

Genetic parameters for pre-weaning traits in Braunvieh cattle

The objective of this study was to estimate genetic parameters for pre-weaning traits of Braunvieh cattle raised under tropical conditions in Brazil. The weight and weight gain parameters were birth weight (BW, N = 9955), weight at 120 days of age (W120, N = 5901), weaning weight at 205 days (WW, N = 6970), weight gain from birth to 205 days (GAIN205, N = 6013), weight gain from birth to 120 days (GAIN120, N = 5135), and weight gain from 120 to 205 days (GAIN85, N = 4482). Variance components were estimated using the animal model with the MTDFREML software. The relationship matrix included 35,188 animals; phenotypic measures were available for 18,688. Direct and maternal heritability increased from birth to weaning, with estimates of 0.23 +/- 0.037, 0.25 +/- 0.050, 0.41 +/- 0.059 for direct heritability for BW, W120 and WW, respectively, 0.08 +/- 0.012, 0.15 +/- 0.032, 0.22 +/- 0.036 for maternal genetic effects, and 0.18, 0.14 and 0.16 for total heritability estimates. For pre-weaning gains, estimates of heritability were 0.36 +/- 0.059, 0.30 +/- 0.059, 0.12 +/- 0.035 for direct genetic effects of the traits GAIN205, GAIN120 and GAIN85, respectively, 0.23 +/- 0.038, 0.17 +/- 0.037, 0.03 +/- 0.029 for estimates of maternal heritability, and 0.12, 0.13, 0.16 for total heritability, respectively. Genetic correlations between weights were greater between measures taken at shorter intervals. This information can be used to optimize the design of programs for genetic improvement of Braunvieh cattle raised under tropical conditions.

Genetic parameters for productive life traits and reproductive efficiency traits at 6 years in Nellore cattle

The objective of the present study was to estimate (co)variance components for length of productive life (LPL) and some alternative reproductive traits of 6-year-old Nellore cattle. The data set contained 57,410 records for age at first calving from Nellore females and was edited to remove animal records with uncertain paternity and cows with just one piece of calving information. Only animals with age at first calving ranging from 23 to 48 months and calving intervals between 11 and 24 months were kept for analysis. LPL and life production ( LP) were used to describe productive life. LPL was defined as the number of months a cow was kept in the herd until she was 6 years old, given that she was alive at first calving and LP was defined as total number of calves in that time. Four traits were used to describe reproductive traits: two breeding efficiencies on original scale were estimated using Wilcox and Tomar functions (BEW and BET, respectively), and two breeding efficiencies transformed (ASBEW and ASBET, respectively), using the function [arcsine (square root (BEi/100))]. Estimates of heritability for measures of LPL and LP were low and ranged from 0.04 to 0.05. Estimates of heritability for breeding efficiencies on original and transformed scales oscillated from 0.18 to 0.32. Estimates of genetic correlations ranged from -0.57 to 0.79 for LPL and other traits and from 0.28 to 0.63 for LP and other traits.

Association of SNPs on CAPN1 and CAST genes with tenderness in Nellore cattle

We examined whether single-nucleotide polymorphisms (SNPs) in the calpain (CAPN) and calpastatin (CAST) genes, described from Bos primigenius taurus, are polymorphic in Nellore cattle. We also looked for a possible association of linkage disequilibrium of this polymorphism with tenderness of the longissimus dorsi muscle after 7, 14 and 21 days of postmortem aging in 638 purebred Nellore bulls. Meat tenderness was measured as Warner-Bratzler shear force. Additive and dominance effects were tested for SNPs of the three genotypic classes; the substitution effect was tested for SNPs with missing genotypic classes. Genotypic and gene frequencies were also calculated for the different SNPs. An increase in tenderness was observed from 7 to 21 days; the average values for shear force at 7, 14 and 21 days of aging were 5.92 +/- 0.06, 4.92 +/- 0.05, and 4.38 +/- 0.04 kg, respectively. All markers showed polymorphism, but there was no CC genotype for CAPN316, and few animals showed the AA genotype for CAPN530. The alleles CAPN4751, UOGCAST1, and WSUCAST were found to have additive and dominance effects for shear force at 7, 14 and 21 days, while CAPN316 showed a substitution effect for shear force at 7 and 21 days. An additive-by-additive epistatic interaction was observed between CAPN4751 and markers on the CAST gene. In conclusion, these markers should be considered for use in breeding programs.

Effects of polymorphisms of LHR and FSHR genes on sexual precocity in a Bos taurus x Bos indicus beef composite population

The purpose of the present research was to investigate the effects of polymorphisms of luteinizing hormone receptor (LHR) and follicle-stimulating hormone receptor (FSHR) genes, evaluated by polymerase chain reaction-restriction fragment length polymorphism in European-Zebu composite beef heifers from six different breed compositions. The polymorphism site analysis from digestion with HhaI and AluI restriction endonucleases allowed the genotype identification for LHR (TT, CT and CC) and FSHR (GG, CG and CC) genes. A high frequency of heterozygous animals was recorded in all breed compositions for both genes, except in two compositions for LHR. The probability of pregnancy (PP) at first breeding was used to evaluate the polymorphism effect on sexual precocity. The PP was analyzed as a binary trait, with a value of 1 (success) assigned to heifers that were diagnosed pregnant by rectal palpation and a value of 0 (failure) assigned to those that were not pregnant at that time. Heterozygous heifers showed a higher pregnancy rate (67 and 66% for LHR and FSHR genes, respectively), but no significant effects were observed for the genes studied (P=0.9188 and 0.8831 for LHR and FSHR, respectively) on the PP. These results do not justify the inclusion of LHR and FSHR restriction fragment length polymorphism markers in selection programs for sexual precocity in beef heifers. Nevertheless, these markers make possible the genotype characterization and may be used in additional studies to evaluate the genetic structure in other bovine populations.