Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions.

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep.

BACKGROUND: The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control. METHODS: Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content. RESULTS: The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep. CONCLUSIONS: The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Brain fatty acid and transcriptome profiles of pig fed diets with different levels of soybean oil.

BACKGROUND: The high similarity in anatomical and neurophysiological processes between pigs and humans make pigs an excellent model for metabolic diseases and neurological disorders. Lipids are essential for brain structure and function, and the polyunsaturated fatty acids (PUFA) have anti-inflammatory and positive effects against cognitive dysfunction in neurodegenerative diseases. Nutrigenomics studies involving pigs and fatty acids (FA) may help us in better understanding important biological processes. In this study, the main goal was to evaluate the effect of different levels of dietary soybean oil on the lipid profile and transcriptome in pigs' brain tissue. RESULTS: Thirty-six male Large White pigs were used in a 98-day study using two experimental diets corn-soybean meal diet containing 1.5% soybean oil (SOY1.5) and corn-soybean meal diet containing 3.0% soybean oil (SOY3.0). No differences were found for the brain total lipid content and FA profile between the different levels of soybean oil. For differential expression analysis, using the DESeq2 statistical package, a total of 34 differentially expressed genes (DEG, FDR-corrected p-value < 0.05) were identified. Of these 34 DEG, 25 are known-genes, of which 11 were up-regulated (log2 fold change ranging from + 0.25 to + 2.93) and 14 were down-regulated (log2 fold change ranging from - 3.43 to -0.36) for the SOY1.5 group compared to SOY3.0. For the functional enrichment analysis performed using MetaCore with the 34 DEG, four pathway maps were identified (p-value < 0.05), related to the ALOX15B (log2 fold change - 1.489), CALB1 (log2 fold change - 3.431) and CAST (log2 fold change + 0.421) genes. A "calcium transport" network (p-value = 2.303e-2), related to the CAST and CALB1 genes, was also identified. CONCLUSION: The results found in this study contribute to understanding the pathways and networks associated with processes involved in intracellular calcium, lipid metabolism, and oxidative processes in the brain tissue. Moreover, these results may help a better comprehension of the modulating effects of soybean oil and its FA composition on processes and diseases affecting the brain tissue.

Rumen and fecal microbiomes are related to diet and production traits in Bos indicus beef cattle.

Background: Ruminants harbor a complex microbial community within their gastrointestinal tract, which plays major roles in their health and physiology. Brazil is one of the largest producers of beef in the world and more than 90% of the beef cattle herds are composed of pure and crossbred Nelore (Bos indicus). Despite its importance to the Brazilian economy and human feeding, few studies have characterized the Nelore microbiome. Therefore, using shotgun metagenomics, we investigated the impact of diet on the composition and functionality of the Nelore microbiome, and explored the associations between specific microbial taxa and their functionality with feed efficiency and methane emission. Results: The ruminal microbiome exhibited significantly higher microbial diversity, distinctive taxonomic profile and variations in microbial functionality compared to the fecal microbiome, highlighting the distinct contributions of the microbiomes of these environments. Animals subjected to different dietary treatments exhibited significant differences in their microbiomes' archaeal diversity and in the abundance of 89 genera, as well as in the functions associated with the metabolism of components of each diet. Moreover, depending on the diet, feed-efficient animals and low methane emitters displayed higher microbial diversity in their fecal microbiome. Multiple genera were associated with an increase or decrease of the phenotypes. Upon analyzing the functions attributed to these taxa, we observed significant differences on the ruminal taxa associated with feed efficient and inefficient cattle. The ruminal taxa that characterized feed efficient cattle stood out for having significantly more functions related to carbohydrate metabolism, such as monosaccharides, di-/oligosaccharides and amino acids. The taxa associated with methane emission had functions associated with methanogenesis and the production of substrates that may influence methane production, such as hydrogen and formate. Conclusion: Our findings highlight the significant role of diet in shaping Nelore microbiomes and how its composition and functionality may affect production traits such as feed efficiency and methane emission. These insights provide valuable support for the implementation of novel feeding and biotechnological strategies.

Multi-omic data integration for the study of production, carcass, and meat quality traits in Nellore cattle.

Data integration using hierarchical analysis based on the central dogma or common pathway enrichment analysis may not reveal non-obvious relationships among omic data. Here, we applied factor analysis (FA) and Bayesian network (BN) modeling to integrate different omic data and complex traits by latent variables (production, carcass, and meat quality traits). A total of 14 latent variables were identified: five for phenotype, three for miRNA, four for protein, and two for mRNA data. Pearson correlation coefficients showed negative correlations between latent variables miRNA 1 (mirna1) and miRNA 2 (mirna2) (-0.47), ribeye area (REA) and protein 4 (prot4) (-0.33), REA and protein 2 (prot2) (-0.3), carcass and prot4 (-0.31), carcass and prot2 (-0.28), and backfat thickness (BFT) and miRNA 3 (mirna3) (-0.25). Positive correlations were observed among the four protein factors (0.45-0.83): between meat quality and fat content (0.71), fat content and carcass (0.74), fat content and REA (0.76), and REA and carcass (0.99). BN presented arcs from the carcass, meat quality, prot2, and prot4 latent variables to REA; from meat quality, REA, mirna2, and gene expression mRNA1 to fat content; from protein 1 (prot1) and mirna2 to protein 5 (prot5); and from prot5 and carcass to prot2. The relations of protein latent variables suggest new hypotheses about the impact of these proteins on REA. The network also showed relationships among miRNAs and nebulin proteins. REA seems to be the central node in the network, influencing carcass, prot2, prot4, mRNA1, and meat quality, suggesting that REA is a good indicator of meat quality. The connection among miRNA latent variables, BFT, and fat content relates to the influence of miRNAs on lipid metabolism. The relationship between mirna1 and prot5 composed of isoforms of nebulin needs further investigation. The FA identified latent variables, decreasing the dimensionality and complexity of the data. The BN was capable of generating interrelationships among latent variables from different types of data, allowing the integration of omics and complex traits and identifying conditional independencies. Our framework based on FA and BN is capable of generating new hypotheses for molecular research, by integrating different types of data and exploring non-obvious relationships.

Integrative Analysis Between Genome-Wide Association Study and Expression Quantitative Trait Loci Reveals Bovine Muscle Gene Expression Regulatory Polymorphisms Associated With Intramuscular Fat and Backfat Thickness.

Understanding the architecture of gene expression is fundamental to unravel the molecular mechanisms regulating complex traits in bovine, such as intramuscular fat content (IMF) and backfat thickness (BFT). These traits are economically important for the beef industry since they affect carcass and meat quality. Our main goal was to identify gene expression regulatory polymorphisms within genomic regions (QTL) associated with IMF and BFT in Nellore cattle. For that, we used RNA-Seq data from 193 Nellore steers to perform SNP calling analysis. Then, we combined the RNA-Seq SNP and a high-density SNP panel to obtain a new dataset for further genome-wide association analysis (GWAS), totaling 534,928 SNPs. GWAS was performed using the Bayes B model. Twenty-one relevant QTL were associated with our target traits. The expression quantitative trait loci (eQTL) analysis was performed using Matrix eQTL with the complete SNP dataset and 12,991 genes, revealing a total of 71,033 cis and 36,497 trans-eQTL (FDR < 0.05). Intersecting with QTL for IMF, we found 231 eQTL regulating the expression levels of 117 genes. Within those eQTL, three predicted deleterious SNPs were identified. We also identified 109 eQTL associated with BFT and affecting the expression of 54 genes. This study revealed genomic regions and regulatory SNPs associated with fat deposition in Nellore cattle. We highlight the transcription factors FOXP4, FOXO3, ZSCAN2, and EBF4, involved in lipid metabolism-related pathways. These results helped us to improve our knowledge about the genetic architecture behind important traits in cattle.

Cross-population selection signatures in Canchim composite beef cattle.

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes.

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Genetic regulators of mineral amount in Nelore cattle muscle predicted by a new co-expression and regulatory impact factor approach.

Mineral contents in bovine muscle can affect meat quality, growth, health, and reproductive traits. To better understand the genetic basis of this phenotype in Nelore (Bos indicus) cattle, we analysed genome-wide mRNA and miRNA expression data from 114 muscle samples. The analysis implemented a new application for two complementary algorithms: the partial correlation and information theory (PCIT) and the regulatory impact factor (RIF), in which we included the estimated genomic breeding values (GEBVs) for the phenotypes additionally to the expression levels, originally proposed for these methods. We used PCIT to determine putative regulatory relationships based on significant associations between gene expression and GEBVs for each mineral amount. Then, RIF was adopted to determine the regulatory impact of genes and miRNAs expression over the GEBVs for the mineral amounts. We also investigated over-represented pathways, as well as pieces of evidences from previous studies carried in the same population and in the literature, to determine regulatory genes for the mineral amounts. For example, NOX1 expression level was positively correlated to Zinc and has been described as Zinc-regulated in humans. Based on our approach, we were able to identify genes, miRNAs and pathways not yet described as underlying mineral amount. The results support the hypothesis that extracellular matrix interactions are the core regulator of mineral amount in muscle cells. Putative regulators described here add information to this hypothesis, expanding the knowledge on molecular relationships between gene expression and minerals.

Correction: Tournaments between markers as a strategy to enhance genomic predictions.

[This corrects the article DOI: 10.1371/journal.pone.0217283.].

Tournaments between markers as a strategy to enhance genomic predictions.

Analysis of a large number of markers is crucial in both genome-wide association studies (GWAS) and genome-wide selection (GWS). However there are two methodological issues that restrict statistical analysis: high dimensionality (p≫n) and multicollinearity. Although there are methodologies that can be used to fit models for data with high dimensionality (eg, the Bayesian Lasso), a big problem that can occurs in this cases is that the predictive ability of the model should perform well for the individuals used to fit the model, but should not perform well for other individuals, restricting the applicability of the model. This problem can be circumvent by applying some selection methodology to reduce the number of markers (but keeping the markers associated with the phenotypic trait) before adjusting a model to predict GBVs. We revisit a tournament-based strategy between marker samples, where each sample has good statistical properties for estimation: n>p and low collinearity. Such tournaments are elaborated using multiple linear regression to eliminate markers. This method is adapted from previous works found in the literature. We used simulated data as well as real data derived from a study with SNPs in beef cattle. Tournament strategies not only circumvent the p≫n issue, but also minimize spurious associations. For real data, when we selected a few more than 20 markers, we obtained correlations greater than 0.70 between predicted Genomic Breeding Values (GBVs) and phenotypes in validation groups of a cross-validation scheme; and when we selected a larger number of markers (more than 100), the correlations exceeded 0.90, showing the efficiency in identifying relevant SNPs (or segregations) for both GWAS and GWS. In the simulation study, we obtained similar results.

Widespread modulation of gene expression by copy number variation in skeletal muscle.

Copy number variation (CNV) is a frequently observed deviation from the diploid state due to duplication or deletion of genomic regions. Although intensively analyzed for association with diseases and production traits, the specific mechanisms and extent by which such variations affect the phenotype are incompletely understood. We present an integrative study on CNV and genome-wide gene expression in Brazilian Bos indicus cattle. We analyzed CNVs inferred from SNP-chip data for effects on gene expression measured with RNA-seq in skeletal muscle samples of 183 steers. Local effects, where expression changes coincided with CNVs in the respective genes, were restricted to immune genes. Distal effects were attributable to several high-impact CNVs that modulated remote expression in an orchestrated and intertwined fashion. These CNVs were located in the vicinity of major skeletal muscle pathway regulators and associated genes were enriched for proteolysis, autophagy, and muscle structure development. From association analysis between CNVs and several meat quality and production traits, we found CNV-associated expression effects to also manifest at the phenotype level. Based on genome sequences of the population founders, we further demonstrate that CNVs with impact on expression and phenotype are passed on from one generation to another.

Candidate genes for male and female reproductive traits in Canchim beef cattle.

BACKGROUND: Beef cattle breeding programs in Brazil have placed greater emphasis on the genomic study of reproductive traits of males and females due to their economic importance. In this study, genome-wide associations were assessed for scrotal circumference at 210 d of age, scrotal circumference at 420 d of age, age at first calving, and age at second calving, in Canchim beef cattle. Data quality control was conducted resulting in 672,778 SNPs and 392 animals. RESULTS: Associated SNPs were observed for scrotal circumference at 420 d of age (435 SNPs), followed by scrotal circumference at 210 d of age (12 SNPs), age at first calving (six SNPs), and age at second calving (four SNPs). We investigated whether significant SNPs were within genic or surrounding regions. Biological processes of genes were associated with immune system, multicellular organismal process, response to stimulus, apoptotic process, cellular component organization or biogenesis, biological adhesion, and reproduction. CONCLUSIONS: Few associations were observed for scrotal circumference at 210 d of age, age at first calving, and age at second calving, reinforcing their polygenic inheritance and the complexity of understanding the genetic architecture of reproductive traits. Finding many associations for scrotal circumference at 420 d of age in various regions of the Canchim genome also reveals the difficulty of targeting specific candidate genes that could act on fertility; nonetheless, the high linkage disequilibrium between loci herein estimated could aid to overcome this issue. Therefore, all relevant information about genomic regions influencing reproductive traits may contribute to target candidate genes for further investigation of causal mutations and aid in future genomic studies in Canchim cattle to improve the breeding program.

Study on the introgression of beef breeds in Canchim cattle using single nucleotide polymorphism markers.

The aim of this study was to evaluate the level of introgression of breeds in the Canchim (CA: 62.5% Charolais-37.5% Zebu) and MA genetic group (MA: 65.6% Charolais-34.4% Zebu) cattle using genomic information on Charolais (CH), Nelore (NE), and Indubrasil (IB) breeds. The number of animals used was 395 (CA and MA), 763 (NE), 338 (CH), and 37 (IB). The Bovine50SNP BeadChip from Illumina panel was used to estimate the levels of introgression of breeds considering the Maximum likelihood, Bayesian, and Single Regression method. After genotype quality control, 32,308 SNPs were considered in the analysis. Furthermore, three thresholds to prune out SNPs in linkage disequilibrium higher than 0.10, 0.05, and 0.01 were considered, resulting in 15,286, 7,652, and 1,582 SNPs, respectively. For k = 2, the proportion of taurine and indicine varied from the expected proportion based on pedigree for all methods studied. For k = 3, the Regression method was able to differentiate the animals in three main clusters assigned to each purebred breed, showing more reasonable according to its biological viewpoint. Analyzing the data considering k = 2 seems to be more appropriate for Canchim-MA animals due to its biological interpretation. The usage of 32,308 SNPs in the analyses resulted in similar findings between the estimated and expected breed proportions. Using the Regression approach, a contribution of Indubrasil was observed in Canchim-MA when k = 3 was considered. Genetic parameter estimation could account for this breed composition information as a source of variation in order to improve the accuracy of genetic models. Our findings may help assemble appropriate reference populations for genomic prediction for Canchim-MA in order to improve prediction accuracy. Using the information on the level of introgression in each individual could also be useful in breeding or crossing design to improve individual heterosis in crossbred cattle.

Iron Content Affects Lipogenic Gene Expression in the Muscle of Nelore Beef Cattle.

Iron (Fe) is an essential mineral for metabolism and plays a central role in a range of biochemical processes. Therefore, this study aimed to identify differentially expressed (DE) genes and metabolic pathways in Longissimus dorsi (LD) muscle from cattle with divergent iron content, as well as to investigate the likely role of these DE genes in biological processes underlying beef quality parameters. Samples for RNA extraction for sequencing and iron, copper, manganese, and zinc determination were collected from LD muscles at slaughter. Eight Nelore steers, with extreme genomic estimated breeding values for iron content (Fe-GEBV), were selected from a reference population of 373 animals. From the 49 annotated DE genes (FDR<0.05) found between the two groups, 18 were up-regulated and 31 down-regulated for the animals in the low Fe-GEBV group. The functional enrichment analyses identified several biological processes, such as lipid transport and metabolism, and cell growth. Lipid metabolism was the main pathway observed in the analysis of metabolic and canonical signaling pathways for the genes identified as DE, including the genes FASN, FABP4, and THRSP, which are functional candidates for beef quality, suggesting reduced lipogenic activities with lower iron content. Our results indicate metabolic pathways that are partially influenced by iron, contributing to a better understanding of its participation in skeletal muscle physiology.

Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle.

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

Selection signatures in Canchim beef cattle.

BACKGROUND: Recent technological advances in genomics have allowed the genotyping of cattle through single nucleotide polymorphism (SNP) panels. High-density SNP panels possess greater genome coverage and are useful for the identification of conserved regions of the genome due to selection, known as selection signatures (SS). The SS are detectable by different methods, such as the extended haplotype homozygosity (EHH); and the integrated haplotype score (iHS), which is derived from the EHH. The aim of this study was to identify SS regions in Canchim cattle (composite breed), genotyped with high-density SNP panel. RESULTS: A total of 687,655 SNP markers and 396 samples remained for SS analysis after the genotype quality control. The iHS statistic for each marker was transformed into piHS for better interpretation of the results. Chromosomes BTA5 and BTA14 showed piHS > 5, with 39 and nine statistically significant SNPs (P < 0.00001), respectively. For the candidate selection regions, iHS values were computed across the genome and averaged within non-overlapping windows of 500 Kb. We have identified genes that play an important role in metabolism, melanin biosynthesis (pigmentation), and embryonic and bone development. CONCLUSIONS: The observation of SS indicates that the selection processes performed in Canchim, as well as in the founder breeds (i.e. Charolais), are maintaining specific genomic regions, particularly on BTA5 and BTA14. These selection signatures regions could be associated with Canchim characterization.

Effect of IGF1, GH, and PIT1 markers on the genetic parameters of growth and reproduction traits in Canchim cattle.

The availability of dense genomic information has increased genome-wide association studies for the bovine species; however research to assess the effect of single genes on production traits is still important to elucidate the genes functions. On this study the association of IGF1, GH, and PIT1 markers with growth and reproductive traits (birth weight, weaning weight, weight at 12 and 18 months of age, preweaning average daily weight gain, age and weight at first calving, and scrotal circumference at 12 and 18 months of age) were assessed by means of the variance component approach. The phenotypes were adjusted and then analyzed under two animal models, one which considered the polygenic and genotype (IGF1, GH or PIT1 markers) effects (Model 1), and the other which considers only the polygenic effect (Model 2). When the likelihood ratio test and the Bonferroni correction was applied at 5 % significance level, the genetic markers for the IGF1, GH, and PIT1 genes did not influence significantly the traits (p > 0.002). However, evidence of association of IGF1 with birth weight (p = 0.06) and GH with weight at first calving (p = 0.03) and with weight at 12 months of age (p = 0.08) was observed. In conclusion we could not confirm the associations between IGF1, GH, and PIT1 and growth traits that were previously reported in Canchim cattle, and no association was observed between these genes and reproductive traits. Future studies involving functional markers of IGF1, GH and PIT1 genes may help to clarify the role of these genes in growth and reproductive processes.

Heat shock and structural proteins associated with meat tenderness in Nellore beef cattle, a Bos indicus breed.

Nellore beef cattle, a Bos indicus (Zebu) breed, is well adapted to tropical conditions and has allowed Brazil to become one of the largest producers of red meat. Nevertheless, B. indicus breeds are reported to have less tender meat than Bos taurus. This study was designed to identify genes associated with meat tenderness and thus provides important information for breeding programs. A group of 138 animals was evaluated for longissimus thoracis muscle shear force (SF). Animals with the highest and lowest SF values (six animals each) were then selected for protein abundance studies. Samples were subjected to two-dimensional gel electrophoresis (2-DE) followed by peptide sequencing through mass spectrometry (MS) to identify differentially expressed proteins associated with SF values. Seventeen differentially expressed spots were observed (p<0.05) between the two groups. The 13 proteins identified included structural proteins (alpha actin-1, MLC1, MLC3, MLC2F and tropomyosin), related to cell organization (HSPB1 and HSP70), metabolism (beta-LG, ACBD6 and Complex III subunit I) and some uncharacterized proteins. Results confirm the existence of differentially expressed proteins associated with SF, which can lead to a better understanding of mechanisms involved in meat tenderness.

Association between JY-1 gene polymorphisms and reproductive traits in beef cattle.

Reproductive traits have a high economic value and it is interesting to include them in the selection objectives of an animal breeding program. These traits generally show low heritability and molecular markers may therefore be used in genetic evaluations to improve the accuracy of predictions. The JY-1 gene is expressed in the oocyte and it is associated with folliculogenesis and early embryo development. It has been suggested to affect reproductive traits. In this study, exons 1 and 2 of the JY-1 gene were studied in 385 Nellore females by PCR-sequencing. Seventeen polymorphisms were identified. After analysis of linkage disequilibrium, association tests were performed between eight SNPs and the occurrence of early pregnancy, age at first calving, days to calving, and reconception of primiparous heifers. Seven SNPs were significant for three traits. The most significant was chr29:12,999 T/A (p=0.003) which was associated with the occurrence of early pregnancy. This SNP might be involved in protein translation inhibition since it affects the initial methionine codon. The JY-1, an oocyte specific gene, influences reproductive traits; further studies investigating other regions of the gene or other genes expressed in tissues of the female reproductive system would be interesting to be performed.