Application of proteomic to investigate the different degrees of meat tenderness in Nellore breed.

This study describes the association between meat tenderness and abundance of soluble muscle proteins in Nellore bulls (Bos indicus) using a proteomic approach. We evaluated shear force (SF) of Longissimus thoracis muscle 24 h after slaughter and selected three experimental groups of animals with moderately tender (TE; SF = 3.9 ± 0.7 kg), moderately tough (TO; SF = 5.6 ± 0.7 kg) and very tough meat (TO+; SF = 7.9 ± 1.4 kg). Proteome was investigated by two-dimensional electrophoresis (2D-PAGE) in combination with electrospray ionization-tandem mass spectrometry (ESI-MS/MS). The metabolic proteins triosephosphate isomerase (TPI1) and phosphoglucomutase 1 (PGM1), the structural protein profilin 1 (PFN1), and cytosol aminopeptidase (LAP3) were up-regulated (P < 0.05) in the TE meat group when compared to the TO and TO+ groups. Actin structural proteins (ACTA1, ACTB, and ACTG1), the oxidative stress protein peroxiredoxin (PRDX6, PRDX2, PRDX1, and PARK7), heat shock protein isoforms, and co-chaperones (CDC37 and STIP1) were up-regulated (P < 0.05) in the TO and TO+ meat groups. In addition, we also identified proteins PFN1, LAP3, PRDX1, PRDX2, HSPD1, and ARHGDIA to be associated with beef tenderness. The results reported herein demonstrated that meat tenderness in Nellore cattle depends on the modulation and expression of a set of proteins involved in different biological pathways. SIGNIFICANCE: The manuscript entitled "Application of proteomic to investigate the different degrees of meat tenderness in Nellore breed" describes a classical proteomics work using two-dimensional gel electrophoresis (2D-PAGE), followed by mass spectrometry coupled to electrospray ionization ion trap (ESI-MS/MS) in order to understand the biochemical engineering involved in the process of meat tenderness. We evaluated shear force (SF) of Longissimus thoracis muscle samples of Nellore cattle (n = 90) and select three experimental groups of animals with moderately tender (TE; SF = 3.9 ± 0.7), moderately tough (TO; SF = 5.6 ± 0.7) and very tough meat (TO+; SF = 7.9 ± 1.4). The proteomic approach allowed observing that meat tenderness is influenced by structural proteins (ACTA1, ACTG1, ACTB, MYL1 and PFN1), co-chaperones (CDC37 and STIP1), heat shock proteins (HSP90AA1, HSP90AB1, HSPD1, HSPA1L, HSPA1A and HSPB1), regulatory protein (ARHGDIA), metabolic proteins (TPI1 and PGM1) and oxidative stress proteins (PRDX1, PRDX2, PRDX6, PARK7). Our results suggest that meat tenderness in Nellore depends on the modulation and expression of a set of proteins involved in different biological pathways.

Gene and protein expression of myosin heavy chain in Nellore cattle comparing growth or meat tenderness traits.

The objective was to investigate gene and protein expression of myosin heavy chain (MyHC) in Nellore cattle slaughtered at different weights (BW) or degrees of meat tenderness. Ninety animals with initial BW 370 ± 37 kg, 24 months of age, were slaughtered after 95 days on feed. We evaluated shear force (SF), myofibrillar fragmentation index, ribeye area, backfat thickness, marbling, color, and cooking losses. Subsequently, 24 animals were selected and divided into four contrasting groups, in which light (BW = 504.58 ± 32.36 kg) versus heavy animals (BW = 604.83 ± 42.97 kg) and animals with tender (SF = 3.88 ± 0.57 kg) versus tough meat (SF = 7.95 ± 1.04 kg) were compared. The MYH7, MYH2 and MYH1 genes were analyzed by real-time PCR. The MyHC isoforms (MyHC-I, MyHC-IIa, and MyHC-IIx) were quantified by SDS-PAGE electrophoresis. We found lower expression of MYH2 and MYH1 genes in heavy compared to light animals and a higher amount of MyHC-I isoform in the tough meat group compared to the tender meat group. Protein expression of MyHC-IIa was higher in the tender meat group. A negative correlation was found of this protein and SF (tenderness), suggesting MyHC-IIa as a biomarker of meat quality.

Exome analysis and functional classification of identified variants in racing Quarter Horses.

The main objectives of this study were to identify and functionally classify SNPs and indels by exome sequencing of animals of the racing line of Quarter Horses. Based on the individual genomic estimated breeding values (GEBVs) for maximum speed index (SImax) obtained for 349 animals, two groups of 20 extreme animals were formed. Of these individuals, 20 animals with high GEBVs for SImax and 19 with low GEBVs for SImax had their exons and 5' and 3' UTRs sequenced. Considering SNPs and indels, 105 182 variants were identified in the expressed regions of the Quarter Horse genome. Of these, 72 166 variants were already known and 33 016 are new variants and were deposited in a database. The analysis of the set of gene variants significantly related (Padjusted  < 0.05) to extreme animals in conjunction with the predicted impact of the changes and the physiological role of protein product pointed to two candidate genes potentially related to racing performance: SLC3A1 on ECA15 and CCN6 on ECA10.

Pedigree analysis and inbreeding effects over morphological traits in Campolina horse population.

Genetic improvement, without control of inbreeding, can go to loss of genetic variability, reducing the potential for genetic gains in the domestic populations. The aim of this study was to analyze the population structure and the inbreeding depression in Campolina horses. Phenotype information from 43 465 individuals was analyzed, data provided by the Campolina Breeders Association. A pedigree file containing 107 951 horses was used to connected the phenotyped individuals. The inbreeding coefficient was performed by use of the diagonal of the relationship matrix and the genealogical parameters were computed using proper softwares. The effective population size was estimated based on the rate of inbreeding and census information, and the stratification of the population was verified by the average relationship coefficient between animals born in different regions of Brazil. The effects of inbreeding on morphological traits were made by the use of inbreeding coefficient as a covariate in the model of random regression. The inbreeding coefficient increased from 1990 on, impacting effective population size and, consequently, shrinking genetic variability. The paternal inbreeding was greater than maternal, which may be attributed to the preference for inbred animals in reproduction. The average genetic relationship coefficient of animals born in different states was lower than individuals born within the same state. The increase in the inbreeding coefficient was negatively associated with all studied traits, showing the importance to avoid genetic losses in the long term. Although results do not indicate a severe narrowing of the population until the present date, the average relationship coefficient shows signs of increase, which could cause a drastic reduction in genetic variability if inbred mating is not successfully controlled in the Campolina horse population.

Genetic study of skin thickness and its association with postweaning growth in Nellore cattle: estimation of the genetic parameters.

The objective of the present study was to estimate genetic parameters for skin thickness (ST) and postweaning weight gain (PWG550) in Nellore cattle. Records were obtained from 152,392 Nellore animals born between 2001 and 2011. ST was measured in the posterior region of the animal's scapula with a millimeter caliper. The animals were assigned to different contemporary groups, formed on the basis of farm, year, sex, feeding regimen at weaning, date of weaning, feeding regimen at 450 days of age, and date of weighing at 450 days of age. The genetic parameters were estimated by Bayesian analysis using the GIBBS1F90 program. The mean ST and PWG550 observed were 7.71 ± 2.04 mm and 115.95 ± 36.17 kg, respectively. The posterior mean estimates of heritability (h2) were 0.12 ± 0.02 and 0.29 ± 0.02 for ST and PWG550, respectively. The posterior mean estimates of the phenotypic, genetic, and environmental correlations between the traits were 0.16 ± 0.0, 0.17 ± 0.02, and 0.17 ± 0.09, respectively. The traits ST and PWG550 showed sufficient additive genetic variance to be used as selection criteria in breeding programs. The low genetic correlation obtained indicates that genes favoring the expression of one trait may not influence the other. Consequently, a selection favoring ST would be less efficient in increasing PWG550.

Genetic parameters for earnings in Quarter Horse.

In this study, we estimated the heritability (h(2)) of earnings in the Quarter Horse in order to evaluate the inclusion of this trait in breeding programs. Records from 14,754 races of 2443 horses from 1978-2009 were provided by Sorocaba Hippodrome, São Paulo, Brazil. All ancestors of the registered horses were included in the pedigree file until the 4th generation. Log-transformed performance measures (LPM) were analyzed for animals aged 2, 3, and 4 years and during their entire career. The h(2) estimates were obtained using a multi-trait model and Gibbs sampling that included the effects of sex, year of race, and animal in all analyses. Five analyses were performed: 1 in which LPM was divided by the number of prizes, 1 in which LPM was divided by the number of race starts, and 3 analyses that included the number of prizes, number of race starts, and both (LPM_cNPS) as covariates. Analysis was performed with and without inclusion of the maternal effect. Models were compared based on the deviance information criterion and LPM_cNPS including maternal effects was found to be the best model. The h(2) estimates and standard deviation obtained using model LPM_cNPS were 0.19 ± 0.08, 0.21 ± 0.08, 0.22 ± 0.09, and 0.21 ± 0.07 for earnings at 2, 3, and 4 years of age and total career, respectively. Our analyses indicate that earnings are subject to selection and can be included in breeding programs to improve the racing performance of Quarter Horses.

Prospection of genomic regions divergently selected in racing line of Quarter Horses in relation to cutting line.

Selection of Quarter Horses for different purposes has led to the formation of lines, including racing and cutting horses. The objective of this study was to identify genomic regions divergently selected in racing line of Quarter Horses in relation to cutting line applying relative extended haplotype homozygosity (REHH) analysis, an extension of extended haplotype homozygosity (EHH) analysis, and the fixation index (F ST) statistic. A total of 188 horses of both sexes, born between 1985 and 2009 and registered at the Brazilian Association of Quarter Horse Breeders, including 120 of the racing line and 68 of the cutting line, were genotyped using single nucleotide polymorphism arrays. On the basis of 27 genomic regions identified as selection signatures by REHH and F ST statistics, functional annotations of genes were made in order to identify those that could have been important during formation of the racing line and that could be used subsequently for the development of selection tools. Genes involved in muscle growth (n=8), skeletal growth (n=10), muscle energy metabolism (n=15), cardiovascular system (n=14) and nervous system (n=23) were identified, including the FKTN, INSR, GYS1, CLCN1, MYLK, SYK, ANG, CNTFR and HTR2B.

Characterization of gene polymorphisms related to immune system physiology in Mangalarga horses / Caracterização de polimorfismos gênicos relacionados à fisiologia do sistema imune em equinos da raça Mangalarga

The objectives of this study were to standardize a PCR-RFLP genotyping method for the AY_731081:g.1900T>C SNP of the equine CD14 gene, and to characterize this SNP and two other polymorphisms (AY_005808: c.1530A>G of the TLR4 gene and AX_463789: g.133T>C of the Cε gene) in Mangalarga horses, in order to contribute to future studies investigating the association between DNA markers and traits related to immune system physiology in this breed. A total of 151 Mangalarga horses of both sexes and variable ages, representative of the population of São Paulo State, were used. PCR-RFLP was found to be adequate for genotyping of the AY_731081: g.1900T>C SNP of the equine CD14 gene. However, this polymorphism is probably not present in Mangalarga horses, thus impairing association studies using this marker in the breed. The population genetic parameters obtained for the TLR4 AY_005808:c.1530A>G and Cε AX_463789:g.133T>C polymorphisms suggest the use of these markers in association studies with immune system-related traits in Mangalarga horses.

Os objetivos deste trabalho foram a padronização da metodologia PCR-RFLP para genotipagem do SNP AY_731081:g.1900T>C do gene CD14 equino, bem como a caracterização em equinos da raça brasileira Mangalarga deste e de outros dois polimorfismos, o AY_005808: c.1530A>G do TLR4 e o AX_463789: g.133T>C do Cε, a fim de promover o embasamento necessário para futuras pesquisas visando à associação entre marcadores de DNA e características relacionadas à fisiologia do sistema imune na raça. Para tanto, foram utilizados 151 animais Mangalarga, de ambos os sexos e de idades variadas, representativos da população do estado de São Paulo. O método de PCR-RFLP mostrou-se adequado para a genotipagem do SNP AY_731081: g.1900T>C do gene CD14 equino. Entretanto, tal polimorfismo provavelmente não ocorre em equinos Mangalarga, impossibilitando estudos de associação com o marcador na raça. Os parâmetros genético-populacionais obtidos para os polimorfismos AY_005808:c.1530A>G do gene TLR4 e o AX_463789:g.133T>C do gene Cε demonstraram a possibilidade de realização de pesquisas.

Characterization and transcriptional analysis of the promoter region of the Duffy blood group, chemokine receptor (DARC) gene in cattle.

The Duffy antigen is the only receptor for Plasmodium vivax, a hemoparasite of the phylum Apicomplexa and the cause of vivax malaria in humans. Resistance to this parasite in the majority of black African individuals and their descendents is due to a mutation in the gene promoter region, which blocks its transcription on erythrocytes. Regarding bovine babesiosis, it is known that taurine breeds are more susceptible to parasite infection than zebuine breeds. In order to verify whether the same human resistance occurs in bovine, the 5' flanking region of the DARC gene was isolated and characterized in Bos indicus and Bos taurus. Four single nucleotide polymorphisms were identified and genotyped (SNP1: EF_647729.1:g.91C>T; SNP2: EF_647729.1:g.405C>T; SNP3: EF_647729.1: g.433A>G and SNP4: EF_647729.1:g.588A>G), which showed significant frequency differences among 99 bovines of each species (n=198). Characterization of the isolated region revealed the presence of 6 putative haplotypes, 14 genotypes, which are formed by haplotypes, and numerous putative transcription factor binding sites. Only the thymine presence on SNPs 1 and 2, more common in B. indicus, was observed to alter some of the sites in this region. Despite this fact, analyses through real-time PCR on bovines that present the most common homozygote genotypes of each species, which contrast for all the polymorphism, revealed no difference on the DARC gene transcription. Thus, in principle, it was concluded that the polymorphisms identified would not be useful as molecular markers in an improvement program for resistance to babesiosis.

Effect of single nucleotide polymorphisms of CAPN1 and CAST genes on meat traits in Nellore beef cattle (Bos indicus) and in their crosses with Bos taurus.

The objectives of this work were to study the segregation of single nucleotide polymorphisms of the calpain 1, large subunit (CAPN1) and calpastatin (CAST) genes in Nellore (Bos indicus) and Nellore xBos taurus beef cattle, as well as to evaluate their effects on meat traits. For this, 300 animals, including 114 Nellore, 67 Angus x Nellore, 44 Rubia Gallega x Nellore, 41 Canchim, 19 Brangus three-way crosses and 15 Braunvieh three-way crosses, were genotyped for the CAPN4751 [AF_248054.2:g.6545C>T (GenBank accession AF248054.2)] and CAST/DdeI [AF_159246.1:g.2959A>G (GenBank accession AF159246.1)] polymorphisms and phenotyped for Ribeye Area, Backfat Thickness, Intramuscular Fat, Shear Force (SF) and Myofibrillar Fragmentation Index (MFI). In relation to the CAPN4751 polymorphism, a frequency of 10.5% was observed for the C allele in the Nellore breed. In the total sample of studied animals, a significant association was found between genotypes and meat tenderness, assessed by SF (P = 0.005) and MFI (P = 0.008), with genotype CT being more favourable than TT. For the CAST/DdeI polymorphism, a frequency of 55.7% was found for the A allele in the Nellore breed. In the total sample, a significant association was observed between genotypes and meat tenderness - SF (P = 0.004) and MFI (P = 0.006), with genotype AA being more favourable than AG. The relationship between genotypes and aged meat tenderness in confluence with the distribution of favourable alleles shows great potential for application of the CAPN4751 and CAST/DdeI polymorphisms in the genetic improvement of the Nellore breed, whilst contributing to the validation, in this breed and in its crosses with B. taurus, of the association results previously described in the literature.

Targeted nucleotide exchange in bovine myostatin gene.

The myostatin gene, known as Growth Differentiation Factor 8 (GDF8), located at chromosome 2 (BTA2) in cattle, is specifically expressed during embryo development and in the adult skeletal muscle. Molecular analysis of this gene reveals the presence of three exons and two introns. Several cattle breeds, such as Piedmontese, Belgian Blue, Blond'Aquitaine, among others, show polymorphisms in this gene, which are directly related to double muscling phenotype. Piedmontese cattle shows a nucleotide transition G --> A (G938A) at exon 3, resulting in the substitution of cysteine to tyrosine, leading to a protein structure change, which determines myostatin inactivation and consequent muscular hypertrophy. The objective of this work was to implant the polymorphism G938A, naturally existent in Piedmontese breed, into in vitro propagated foetal myoblasts, from Nellore cattle. Single strand DNA (ssDNA) oligonucleotides were used to direct the same nucleotidic transition (G938A) to exon 3. Two transfection protocols (cationic lipid solution and electroporation) were tested and, 48 hours after transfection, RNA and DNA were extracted from myoblasts. Reverse transcription and polymerase chain reaction (PCR) were performed, using primers flanking the mutation region. The PCR products were cloned and analyzed by DNA sequencing, and it was possible to detect the nucleotidic CT transition at position 938, in the electroporated myoblasts. The existence of a positive signal in the transfection indicates that ssDNA oligonucleotides can be used to introduce this point mutation in specific functional gene sites.

Identification and characterization of polymorphisms within the 5' flanking region, first exon and part of first intron of bovine GH gene.

The aim of the present study was to identify and characterize polymorphisms within the 5' flanking region, first exon and part of first intron of the bovine growth hormone gene among different beef cattle breeds: Nelore (n = 25), Simmental (n = 39), Simbrasil (n = 24), Simmental x Nelore (n = 30), Canchim x Nelore (n = 30) and Angus x Nelore (n = 30). Two DNA fragments (GH1, 464 bp and GH2, 453 bp) were amplified by polymerase chain reaction and then used for polymorphism identification by SSCP. Within the GH1 fragment, five polymorphisms were identified, corresponding to three different alleles: GH1.1, GH1.2 and GH1.3 (GenBank: AY662648, AY662649 and AY662650, respectively). These allele sequences were aligned and compared with bovine GH gene nucleotide sequence (GenBank: M57764 and AF118837), resulting in the identification of five insertion/deletions (INDELs) and five single nucleotide polymorphisms (SNPs). In the GH2 fragment two alleles were identified, GH2.1 and GH2.2 (GenBank: AY662651 and AY662652, respectively). The allele sequences were compared with GenBank sequences (M57764, AF007750 and AH009106) and three INDELs and four SNPs were identified. In conclusion, we were able to identify six new polymorphisms of the bovine GH gene (one INDEL and five SNPs), which can be used as molecular markers in genetic studies.

Effects of polymorphic microsatellites in the regulatory region of IGF1 and GHR on growth and carcass traits in beef cattle.

Growth hormone (GH), insulin-like growth factors 1 and 2 (IGF1 and IGF2) and their associated binding proteins and transmembrane receptors (GHR, IGF1R and IGF2R) play an important role in the physiology of mammalian growth. The objectives of the present study were to estimate the allele and genotype frequencies of microsatellite markers located in the 5'-regulatory region of the IGF1 and GHR genes in beef cattle belonging to different genetic groups and to determine effects of these markers on growth and carcass traits in these animals under an intensive production system. For this purpose, genotyping was performed on 384 bulls including 79 Nellore, 30 Canchim (5/8 Charolais + 3/8 Zebu) and 275 crossbred animals originating from crosses of Simmental (1/2 Simmental, n = 30) and Angus (1/2 Angus, n = 245) sires with Nellore females. The effects of substituting L allele for S allele of GHR microsatellite across Nellore, Canchim and 1/2 Angus were significant for weight gain and body weight (P < 0.05). The IGF1 microsatellite allele substitutions of 229 for 225 within Nellore group and of 225 for 229 within 1/2 Angus were not significant for any of the traits.