Polledness in Argentinean Creole cattle, five centuries surviving.

Polledness has been shown to have autosomal Mendelian inheritance, with the polled locus being dominant to the horned locus. This trait was mapped to the BTA1 centromeric end in several breeds. One of the distinctive attributes of Creole cattle, such as the Argentinean Creole, is the presence of long, lyre-shaped horns. However, polled native animals were reported before the introduction of modern selected European breeds. Here, we studied the origin of the polled mutation, either independent or introgressed, in a Creole line from the Creole cattle founder group at the IIACS-INTA Leales Experimental Station (northwest Argentina). The study sample (65 animals: 26 horned and 39 polled) was genotyped using high-density SNP microarrays and three previously reported genetic markers (P202 ID , P80kb ID and PG ). A genome-wide association study, selection signatures, linkage disequilibrium analysis and copy number variations were used to detect the responsible region and the segregating haplotypes/alleles. The interval mapped in the Leales herd (1.23-2.13 Mb) overlapped with the region previously reported in several European cattle breeds, suggesting that the same locus could be segregating in this population. The previously reported variants PF and PG were not detected, thus dismissing the Holstein-Friesian and Nellore origins of the polled phenotype in this native breed. Conversely, the presence of the Celtic variant PC suggests an almost complete co-segregation. The cluster analysis rejected the hypothesis of recent introgression, which is compatible with the historical record of polled Creole cattle in northwest Argentina.

Beyond genomic selection: The animal model strikes back (one generation)!

Genome inheritance is by segments of DNA rather than by independent loci. We introduce the ancestral regression (AR) as a recursive system of simultaneous equations, with grandparental path coefficients as novel parameters. The information given by the pedigree in the AR is complementary with that provided by a dense set of genomic markers, such that the resulting linear function of grandparental BV is uncorrelated to the average of parental BV in the absence of inbreeding. AR is then connected to segmental inheritance by a causal multivariate Gaussian density for BV. The resulting covariance structure (Σ) is Markovian, meaning that conditional on the BV of parents and grandparents, the BV of the animal is independent of everything else. Thus, an algorithm is presented to invert the resulting covariance structure, with a computing effort that is linear in the number of animals as in the case of the inverse additive relationship matrix.

Characterization of bovine MHC class II DRB3 diversity in South American Holstein cattle populations.

Holstein cattle dominate the global milk production industry because of their outstanding milk production, however, this breed is susceptible to tropical endemic pathogens and suffers from heat stress and thus fewer Holstein populations are raised in tropical areas. The bovine major histocompatibility complex (BoLA)-DRB3 class II gene is used as a marker for disease and immunological traits, and its polymorphism has been studied extensively in Holstein cattle from temperate and cold regions. We studied the genetic diversity of the BoLA-DRB3 gene in South American Holstein populations to determine whether tropical populations have diverged from those bred in temperate and cold regions by selection and/or crossbreeding with local native breeds. We specifically studied Exon 2 of this gene from 855 South American Holstein individuals by a polymerase chain reaction (PCR) sequence-based typing method. We found a high degree of gene diversity at the allelic (Na > 20 and He > 0.87) and molecular (π > 0.080) levels, but a low degree of population structure (FST = 0.009215). A principal components analysis and tree showed that the Bolivian subtropical population had the largest genetic divergence compared with Holsteins bred in temperate or cold regions, and that this population was closely related to Bolivian Creole cattle. Our results suggest that Holstein genetic divergence can be explained by selection and/or gene introgression from local germplasms. This is the first examination of BoLA-DRB3 in Holsteins adapted to tropical environments, and contributes to an ongoing effort to catalog bovine MHC allele frequencies by breed and location.

Assessment of biodiversity in Chilean cattle using the distribution of major histocompatibility complex class II BoLA-DRB3 allele.

Bovine leukocyte antigens (BoLAs) are used extensively as markers for bovine disease and immunological traits. In this study, we estimated BoLA-DRB3 allele frequencies using 888 cattle from 10 groups, including seven cattle breeds and three crossbreeds: 99 Red Angus, 100 Black Angus, 81 Chilean Wagyu, 49 Hereford, 95 Hereford × Angus, 71 Hereford × Jersey, 20 Hereford × Overo Colorado, 113 Holstein, 136 Overo Colorado, and 124 Overo Negro cattle. Forty-six BoLA-DRB3 alleles were identified, and each group had between 12 and 29 different BoLA-DRB3 alleles. Overo Negro had the highest number of alleles (29); this breed is considered in Chile to be an 'Old type' European Holstein Friesian descendant. By contrast, we detected 21 alleles in Holstein cattle, which are considered to be a 'Present type' Holstein Friesian cattle. Chilean cattle groups and four Japanese breeds were compared by neighbor-joining trees and a principal component analysis (PCA). The phylogenetic tree showed that Red Angus and Black Angus cattle were in the same clade, crossbreeds were closely related to their parent breeds, and Holstein cattle from Chile were closely related to Holstein cattle in Japan. Overall, the tree provided a thorough description of breed history. It also showed that the Overo Negro breed was closely related to the Holstein breed, consistent with historical data indicating that Overo Negro is an 'Old type' Holstein Friesian cattle. This allelic information will be important for investigating the relationship between major histocompatibility complex (MHC) and disease.

Development of an ELA-DRA gene typing method based on pyrosequencing technology.

The polymorphism of equine lymphocyte antigen (ELA) class II DRA gene had been detected by polymerase chain reaction-single-strand conformational polymorphism (PCR-SSCP) and reference strand-mediated conformation analysis. These methodologies allowed to identify 11 ELA-DRA exon 2 sequences, three of which are widely distributed among domestic horse breeds. Herein, we describe the development of a pyrosequencing-based method applicable to ELA-DRA typing, by screening samples from eight different horse breeds previously typed by PCR-SSCP. This sequence-based method would be useful in high-throughput genotyping of major histocompatibility complex genes in horses and other animal species, making this system interesting as a rapid screening method for animal genotyping of immune-related genes.

Computational prediction of nsSNPs effects on protein function and structure, a prioritization approach for further in vitro studies applied to bovine GSTP1.

The development of high-throughput technologies in the last decade produced an exponential increase in the amount of biological data available. The case of redox biology and apoptosis is not an exception, and nowadays there is a need to integrate information from multiple "omics" studies. Therefore, validation of proposed discoveries is essential. However, the study in biological systems of the effect of the massive amounts of sequence variation data generated with next-generation sequencing (NGS) technologies can be a very difficult and expensive process. In this context, the present study aimed to demonstrate the advantages of a computational methodology to systematically analyze the structural and functional effects of protein variants, in order to prioritize further studies. This approach stands out for its easy implementation, low costs and low time consumed. First, the possible impact of mutations on protein structure and function was tested by a combination of tools based on evolutionary and structural information. Next, homology modeling was performed to predict and compare the 3D protein structures of unresolved amino acid sequences obtained from genomic resequencing. This analysis applied to the bovine GSTP1 allowed to determine that some of amino acid substitutions may generate important changes in protein structure and function. Moreover, the haplotype analysis highlighted three structure variants worthwhile studying through in vitro or in vivo experiments.

Evidence of positive selection towards Zebuine haplotypes in the BoLA region of Brangus cattle.

The Brangus breed was developed to combine the superior characteristics of both of its founder breeds, Angus and Brahman. It combines the high adaptability to tropical and subtropical environments, disease resistance, and overall hardiness of Zebu cattle with the reproductive potential and carcass quality of Angus. It is known that the major histocompatibility complex (MHC, also known as bovine leucocyte antigen: BoLA), located on chromosome 23, encodes several genes involved in the adaptive immune response and may be responsible for adaptation to harsh environments. The objective of this work was to evaluate whether the local breed ancestry percentages in the BoLA locus of a Brangus population diverged from the estimated genome-wide proportions and to identify signatures of positive selection in this genomic region. For this, 167 animals (100 Brangus, 45 Angus and 22 Brahman) were genotyped using a high-density single nucleotide polymorphism array. The local ancestry analysis showed that more than half of the haplotypes (55.0%) shared a Brahman origin. This value was significantly different from the global genome-wide proportion estimated by cluster analysis (34.7% Brahman), and the proportion expected by pedigree (37.5% Brahman). The analysis of selection signatures by genetic differentiation (F st ) and extended haplotype homozygosity-based methods (iHS and Rsb) revealed 10 and seven candidate regions, respectively. The analysis of the genes located within these candidate regions showed mainly genes involved in immune response-related pathway, while other genes and pathways were also observed (cell surface signalling pathways, membrane proteins and ion-binding proteins). Our results suggest that the BoLA region of Brangus cattle may have been enriched with Brahman haplotypes as a consequence of selection processes to promote adaptation to subtropical environments.

Refining genomewide association for growth and fat deposition traits in an F pig population.

The identification of genomic regions that affect additive genetic variation and contain genes involved in controlling growth and fat deposition has enormous impact in the farm animal industry (e.g., carcass merit and meat quality). Therefore, a genomewide association study was implemented in an F pig population using a 60,000 SNP marker panel for traits related to growth and fat deposition. Estimated genomic EBV were linearly transformed to calculate SNP effects and to identify genomic positions possibly associated with the genetic variability of each trait. Genomic segments were then defined considering the markers included in a region 1 Mb up- and downstream from the SNP with the smallest -value and a false discovery rate < 0.05 for each trait. The significance for each 2-Mb segment was tested using the Bonferroni correction. Significant SNP were detected on SSC2, SSC3, SSC5, and SSC6, but 2-Mb segment significant effects were observed on SSC3 for weight at birth (wt_birth) and on SSC6 for 10th-rib backfat and last-rib backfat measured by ultrasound at different ages. Furthermore, a 6-Mb segment on SSC6 was also considered because the 2-Mb segments for 10 different fat deposition traits were overlapped. Although the segment effects for each trait remain significant, the proportion of additive variance explained by this larger segment was slightly smaller in some traits. In general, the results confirm the presence of genetic variability for wt_birth on SSC3 (18.0-20.2 Mb) and for fat deposition traits on SSC6 (133.8-136.0 Mb). Within these regions, fibrosin () and myosin light chain, phosphorylatable, fast skeletal muscle () genes could be considered as candidates for the wt_birth signal on SSC3, and the SERPINE1 mRNAbinding protein 1 gene () may be a candidate for the fat deposition trait signals on SSC6.

Effectiveness of a 95 SNP panel for the screening of breed label fraud in the Chinese meat market.

Breed assignment has proved to be useful to control meat trade and protect the value of special productions. Meat-related frauds have been detected in China; therefore, 95 SNPs selected from the ISAG core panel were evaluated to develop an automated and technologically updated tool to screen breed label fraud in the Chinese meat market. A total of 271 animals from four Chinese yellow cattle (CYC) populations, six Bos taurus breeds, two Bos indicus and one composite were used. The allocation test distinguished European, Japanese and Zebu breeds, and two Chinese genetic components. It correctly allocated Japanese Black, Zebu and British breeds in 100, 90 and 89% of samples, respectively. CYC evidenced the Zebu, Holstein and Limousin introgression. The test did not detect CYC components in any of the 25 samples from Argentinean butchers. The method could be useful to certify Angus, Hereford and Japanese Black meat, but a modification in the panel would be needed to differentiate other breeds.

Growth, carcass and meat quality traits in beef from Angus, Hereford and cross-breed grazing steers, and their association with SNPs in genes related to fat deposition metabolism.

Grazing steers from Angus and Hereford breeds, their cross-breeds and a three-way cross-breed (Limousin × Angus-Hereford) were measured for growth, carcass and meat quality traits. Breed effects were studied, and the association of SNPs with fat deposition and fatty acid (FA) composition (leptin, melanocortin-4 receptor, stearoyl-CoA desaturase, FA synthase and thyroglobulin) was tested. Limousin cross-breed showed the greatest final body weight, ultrasound rib eye area, dressing percentage, carcass and leg length, and the lowest backfat thickness and intramuscular fat content. Genetic groups had similar pH, shear force, cooking loss, L* and b* and n-6:n-3 ratio. Meat from 1/2-Angus presented greater a* than Limousin cross-breed. Whereas Angus had the highest total SFA content, Hereford had the lowest total SFA and the highest total MUFA. Limousin cross-breed had greater content of several individual PUFAs, total PUFA, n-6 and n-3 FA than Angus and 1/2-Angus. Leptin and FA synthase were associated with some FAs, supporting their influence over fat metabolism for grazing animals.

Study of the influence of genes related to muscle oxidative processes on beef color.

The biochemical bases of meat color are determined by the concentration and redox state of myoglobin, hemoglobin, cytochromes, and other pigments. Post-mortem depletion of cellular oxygen results in oxidative stresses that consume NADH and affects reducing activity, while enzymatic detoxification influences the cellular oxidative processes, both affecting meat color. The aim of this work was to study the influence of several genes related to cellular oxidative processes that could affect CIELAB meat color parameters. The study was performed in steers that received a grass-based diet combined with grain, hays and silages. Results suggest a possible link between colorimetric parameters (a*, b* and chroma) and SNPs in the GSTP1 gene (P<0.05). Although the influence of the enzymes, encoded by GSTP1 gene, on meat color has been proposed previously at biochemical level and protein expression level, further association studies in different populations and functional studies of proteins are needed to confirm the genetic determination of that gene on meat color.

Foreign meat identification by DNA breed assignment for the Chinese market.

Methods for individual identification are usually employed for traceability, whereas breed identification is useful to detect commercial frauds. In this study, Chinese Yellow Cattle (CYC) samples plus data from six Bos taurus breeds, two Bos indicus breeds, and one composite breed were used to develop an allocation test based on 22 microsatellites. The test allowed discriminating all foreign breeds from the CYC, although some CYC individuals were wrongly allocated as Limousin or Holstein, probably due to the recent introduction of these breeds into China. In addition, CYC evidenced a previously reported Zebu cline (south-north) and a possible structure within the B. taurus component that should be confirmed. An independent test performed with meat samples of unknown breed origin from Argentina allocated 92% of them to either Angus, Hereford, or their crossbreed, but none was identified as CYC. We conclude that the test is a suitable tool to certify meat of foreign breed origin and to detect adulterations of CYC beef labeled as imported meat.

Genetic and management factors affecting beef quality in grazing Hereford steers.

Attributes contributing to differences in beef quality of 206 Hereford steers finished on pasture were assessed. Beef quality traits evaluated were: Warner-Bratzler meat tenderness and muscle and fat color at one and seven days after slaughter and trained sensory panel traits (tenderness, juiciness, flavor, and marbling) at seven days. Molecular markers were CAPN1 316 and an SNP in exon 2 on the leptin gene (E2FB). Average daily live weight gain, ultrasound monthly backfat thickness gain and rib-eye area gain were estimated. Molecular markers effects on meat quality traits were analyzed by mixed models. Association of meat quality with post weaning growth traits was analyzed by canonical correlations. Muscle color and marbling were affected by CAPN1 316 and E2FB and Warner-Bratzler meat tenderness by the former. The results confirm that marker assisted selection for tenderness is advisable only when beef aging is a common practice. The most important sources of variation in tenderness and color of meat remained unaccounted for.

Characterization and validation of bovine gonadotripin releasing hormone receptor (GNRHR) polymorphisms.

Gonadotropin releasing hormone and its receptor (GNRHR) play a critical role in sexual differentiation and reproduction. Available evidence shows a strong genetic component in the timing of puberty. In bovines, there are significant differences within and among beef breeds in the time when bulls reach puberty. Despite its economic importance, there are not many SNPs or genetic markers associated with this characteristic. The aims of the study were to identify DNA polymorphism in the bovine GNRHR by re-sequencing analysis, determine haplotype phases, and perform a population study in a selected tag SNP in six breeds. Eight SNPs were detected, including: one in the Upstream Regulatory Region (URR), five in the coding regions, and two in non-coding regions. This polymorphism level corresponds to one variant every 249.4bp and a global nucleotide diversity of 0.385. Two haplogroups comprising nine haplotypes and two linkage blocks were detected. Despite 5 tag SNPs were required to capture all variability, just one SNP allowed to define both haplogroups, and only two SNPs were needed to differentiate the most common haplotypes. An additional taq SNP was necessary to identify both URR variants. Allele-frequency analysis of a selected taq SNP among breeds showed a geographical cline. European Bos taurus breeds had lower frequencies of the C allele than B. indicus type cattle, while Creole cattle and Wagyu breeds had intermediate frequency. There was a significant correlation between frequency profile and timing of puberty among the studied breeds, which seems to suggest that genetic variation within bovine GNRHR gene could explain at least part of the reported variability.

History and selection imprinting on genetic relationships among bovine breeds analyzed trough five genes related with marbling.

Many candidate genes have been suggested as responsible for marbling in beef cattle, for instance diacylglycerol O-acyltransferase 1, thyroglobulin, growth hormone, leptin and stearoyl CoA desaturase. The objective of the present work was to evaluate the polymorphisms of five SNPs of these candidate genes in 389 animals of 18 Bos Taurus and Bos indicus breeds. The obtained results were compared with ones previously obtained with STRs and loci related to milk production in these populations. Moreover we analyzed whether the phylogenies reconstructed using SNPs associated with marbling resulted in the known tree topology. The tree constructed with UPGMA, using genetic distance D(A), exhibit a topology partially consistent with the historical origin of breeds. The result observed in the Correspondence Analysis coincided with the topology of the UPGMA tree. This work allowed us to evaluate the five SNPs genetic diversity and to demonstrate that the grouping of the breeds may be the result of its history, selection process, or both at once.

Effect of consanguinity on Argentinean Angus beef DNA traceability.

Since the 1990s several authors have envisaged the use of DNA to certify meat origin. Two major parameters must be assessed before a DNA based traceability protocol can be implemented in the food chain: (i) the information content of a DNA marker set in a specific livestock breed or group of breeds; (ii) the minimum number of DNA markers needed to obtain a statistically acceptable match probability. The objective of the present work was to establish the effect of different levels of inbreeding in the matching efficiency, and the minimum number of microsatellite markers needed, in a DNA based meat traceability program, starting from an 11-microsatellite marker panel. Samples were obtained from beef production farms in South America, where animals are typically bred under pasture-based extensive conditions. Three groups of animals with different consanguinity rates were sampled. Exclusion power (Q) was higher than 0.999998 and match probability lower than 3.01E-08, for the whole set of markers within each group. Both values were affected by consanguinity. To reach a two mismatch criteria exclusion power (Q(2)) of 99.99, six markers were needed in unrelated animals whereas seven markers were needed in related animals. To reach Q(2)=99.9999, 8 and 10 microsatellite markers, respectively, were needed. In general, one or two more microsatellite markers were needed to identify consanguineous animals. This study proved the DNA marker set used to be suitable for the identification of the meat from all slaughtered animals in Argentina, per week, month, and year.