Interaction between poly(A)-binding protein PABPC4 and nuclear receptor corepressor NCoR1 modulates a metabolic stress response.

Mitochondria are organelles known primarily for generating ATP via the oxidative phosphorylation process. Environmental signals are sensed by whole organisms or cells and markedly affect this process, leading to alterations in gene transcription and, consequently, changes in mitochondrial function and biogenesis. The expression of mitochondrial genes is finely regulated by nuclear transcription factors, including nuclear receptors and their coregulators. Among the best-known coregulators is the nuclear receptor corepressor 1 (NCoR1). Muscle-specific knockout of NCoR1 in mice induces an oxidative phenotype, improving glucose and fatty acid metabolism. However, the mechanism by which NCoR1 is regulated remains elusive. In this work, we identified the poly(A)-binding protein 4 (PABPC4) as a new NCoR1 interactor. Unexpectedly, we found that silencing of PABPC4 induced an oxidative phenotype in both C2C12 and MEF cells, as indicated by increased oxygen consumption, mitochondria content, and reduced lactate production. Mechanistically, we demonstrated that PABPC4 silencing increased the ubiquitination and consequent degradation of NCoR1, leading to the derepression of PPAR-regulated genes. As a consequence, cells with PABPC4 silencing had a greater capacity to metabolize lipids, reduced intracellular lipid droplets, and reduced cell death. Interestingly, in conditions known to induce mitochondrial function and biogenesis, both mRNA expression and PABPC4 protein content were markedly reduced. Our study, therefore, suggests that the lowering of PABPC4 expression may represent an adaptive event required to induce mitochondrial activity in response to metabolic stress in skeletal muscle cells. As such, the NCoR1-PABPC4 interface might be a new road to the treatment of metabolic diseases.

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle.

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (ß = 3.02; P-value = 3.51E-03) and the residual feed intake trait (ß = - 3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore cattle.

Genomewide association study for production and meat quality traits in Canchim beef cattle.

The commercial value of the bovine carcass is determined by a set of traits, such as weight, yield, back fat thickness, and marbling; therefore, the genetic improvement of growth, meat, and carcass quality traits is an important tool to add value to the supply chain. Genomewide association studies (GWAS) enable the identification of loci that control phenotypic expression of quantitative traits (QTL). Therefore, the objective of this work was to perform a GWAS to identify genomic regions and genes associated with growth, carcass traits, and meat quality in Canchim beef cattle. These traits were yearling weight (YW), rib eye area (REA), back fat thickness (BFT), and marbling (MARB). To increase sample size and marker density, genotype imputation was performed, and only markers imputed with greater than 95% accuracy were used. Genomewide association study was performed using a Bayesian approach, by the Bayes B statistical method, incorporating genotypes and phenotypes from 614 animals from both the Canchim breed and the MA genetic group (offspring of Charolais bulls and one-half Canchim + one-half Zebu cows). This investigation identified 1 and 4 genomic regions explaining 0.23 and 7.35% of the genetic variance for REA and YW, respectively. These regions harbor a total of 19 genes, 7 of which were classified for biological functions by functional analysis. Significant associations were not observed for BFT and MARB. The identification of QTL that had been previously described in the literature reinforces associations found in this study.

Genomewide association mapping and pathway analysis of meat tenderness in Polled Nellore cattle.

Brazil is one of the world's largest beef exporters, although the product has a low price due to quality issues. The meat exported by Brazil is considered medium and low quality by international buyers, mainly due to lack of tenderness. The predominant Zebu breeds (80% Nellore) are known for producing tougher beef than taurine breeds. Nonetheless, some studies have shown that there is substantial genetic variability for tenderness within the Nellore breed, although it is a difficult trait to improve by conventional selection methods. Therefore, the aim of this study was to perform a genomewide association study (GWAS) and a gene set enrichment analysis to identify genomic regions and biologically relevant pathways associated with meat tenderness in Polled Nellore cattle. Data consisted of Warner-Bratzler shear force values of LM from 427 Polled Nellore animals divided into 3 experimental slaughters (years 2005, 2008, and 2010). The animals were genotyped with either the Illumina BovineHD BeadChip (777k, on 61 samples) or the GGP Indicus HD chip (77k, on 366 samples). Single nucleotide polymorphisms were excluded when the call rate was <90%, the Hardy-Weinberg proportions -value was <1% (Fisher exact test, Bonferroni adjusted), and the minor allele frequency was <1%. Imputation from the GGP Indicus HD chip to the Illumina BovineHD BeadChip was performed using the FImput program. Genomewide association analysis was performed using the Efficient Mixed Model Association eXpedited (EMMAx) and the population parameters previously determined (P3D) methods. The GWAS was complemented with a gene set enrichment analysis performed using the FatiGO procedure. Significant markers ( < 0.0001) explaining a larger proportion of variation than other significant SNPs were located on chromosomes 3, 13, 17, 20, 21, and 25, indicating QTL associated with meat tenderness throughout the genome. Additionally, gene set analysis identified 22 Gene Ontology functional terms and 2 InterPro entries that showed significant enrichment of genes associated with tenderness. The functional categories included protein tyrosine and serine/threonine kinase activity, calcium ion binding, lipid metabolic process, and growth factors, among others. These results help to elucidate the genetic architecture and metabolic pathways underlying this trait, which is of extreme economic and social importance to Brazil, because Nellore is the dominant beef cattle breed in the country.

Genotypic and allelic frequencies of gene polymorphisms associated with meat tenderness in Nellore beef cattle.

The objectives of this study were to characterize the allelic and genotypic frequencies of polymorphisms in the µ-calpain and calpastatin genes, and to assess their association with meat tenderness and animal growth in Nellore cattle. We evaluated 605 Nellore animals at 24 months of age, on average, at slaughter. The polymorphisms were determined for the molecular markers CAPN316, CAPN530, CAPN4751, CAPN4753, and UOGACAST1. Analyses of meat tenderness at 7, 14, and 21 days of maturation were performed in samples of longissimus thoracis obtained between the 12th and 13th rib and sheared using a Warner Bratzler Shear Force. Significant effects were observed for meat tenderness at days 7, 14, and 21 of maturation for the marker CAPN4751, at day 21 for the marker CAPN4753, and at days 14 and 21 for the marker UOGCAST1. For genotypic combinations of markers, the results were significant for the combination CAPN4751/UOGCAST1 in the three maturation periods and CAPN4753/UOGCAST1 at days 14 and 21 of maturation.

Allele- and parent-of-origin-specific effects on expression of the KCNJ11 gene: A candidate for meat tenderness in cattle.

In contrast to the Mendelian inheritance model, parental alleles can contribute unequally to gene expression, which may result in phenotypic variance among individuals and bias in the predicted additive effect of molecular markers associated with production traits. Given the need to understand the effects of allelic variation and parent-of-origin effects on the expression of genes with a commercial interest in cattle, we analyzed the expression of KCNJ11 (potassium inwardly rectifying channel, subfamily J, member 11), which was previously described as a functional candidate gene for meat tenderness. Allele-specific and parent-of-origin-dependent expression of this gene were assessed in bovine muscle using the rs379610823 single nucleotide polymorphism as a reference. Biallelic expression was observed; however, the T allele was expressed at significantly higher levels than the C allele. Furthermore, increased expression of KCNJ11 was found in animals harboring the maternal T allele. This study is the first to describe the differential allelic expression of bovine KCNJ11. Our findings are important for understanding the mechanisms that underlie the pattern of KCNJ11 expression and its potential impact on the phenotypic variation of meat tenderness in Nelore beef cattle. This reinforces the need for further investigation of allelic- and parent-of-origin expression deviation in genetic markers eligible for the selection of target traits.

Variation in myogenic differentiation 1 mRNA abundance is associated with beef tenderness in Nelore cattle.

The myogenic differentiation 1 gene (MYOD1) has a key role in skeletal muscle differentiation and composition through its regulation of the expression of several muscle-specific genes. We first used a general linear mixed model approach to evaluate the association of MYOD1 expression levels on individual beef tenderness phenotypes. MYOD1mRNA levels measured by quantitative polymerase chain reactions in 136 Nelore steers were significantly associated (P ≤ 0.01) with Warner-Bratzler shear force, measured on the longissimus dorsi muscle after 7 and 14 days of beef aging. Transcript abundance for the muscle regulatory gene MYOD1 was lower in animals with more tender beef. We also performed a co-expression network analysis using whole transcriptome sequence data generated from 30 samples of longissimus muscle tissue to identify genes that are potentially regulated by MYOD1. The effect of MYOD1 gene expression on beef tenderness may emerge from its function as an activator of muscle-specific gene transcription such as for the serum response factor (C-fos serum response element-binding transcription factor) gene (SRF), which determines muscle tissue development, composition, growth and maturation.

Candidate genes for carcass traits in a tropical-adapted Brazilian composite beef breed.

Backfat thickness (BFT) and ribeye area (REA) are important production traits but, because they are measured late in the animal's life, they have not been efficiently included in breeding programs. The aim of this study was to evaluate whether single nucleotide polymorphisms (SNPs) mapped to the leptin, PPARGC1A, PSMC1, CRH, and FABP4 genes, which influence BFT and REA in Canchim cattle, a composite beef breed (5/8 Charolais + 3/8 Zebu). BFT and REA phenotypic records were obtained by ultrasound measurements from 18-month-old animals. All SNP markers were genotyped by restriction fragment length polymorphism-polymerase chain reaction. Restricted maximum likelihood analysis revealed that the non-synonymous SNP located in exon 2 of the FABP4 gene has an additive effect on BFT (P ≤ 0.05). Significant allele substitution effects showed that the substitution of G by A may lead to a decrease of 0.1055 mm in mean BFT. This information can be used for inclusion of this trait-associated marker in commercial SNP panels.

Single nucleotide polymorphisms in candidate genes associated with gastrointestinal nematode infection in goats.

Cytokines are small cell-signaling proteins that play an important role in the immune system, participating in intracellular communication. Four candidate genes of the cytokine family (IL2, IL4, IL13, and IFNG) were selected to identify Single Nucleotide Polymorphisms (SNPs) that might be associated with resistance to gastrointestinal endoparasites in goats. A population of 229 goats, F2 offspring from an F1 intercross was produced by crossing pure Saanen goats, considered as susceptible to gastrointestinal endoparasites, with pure Anglo-Nubian goats, considered resistant. Blood was collected for DNA extraction and fecal samples were also collected for parasite egg count. Polymorphisms were prospected by sequencing animals with extreme phenotype for fecal egg count (FEC) distribution. The association between SNPs and phenotype was determined by using the Fisher exact test with correction for multiple tests. Three of the 10 SNPs were identified as significant (P ≤ 0.03). They were found in intron 1 of IL2 (ENSBTA00000020883), intron 3 of IL13 (ENSBTA00000015953) and exon 3 of IFNG (ENSBTA00000012529), suggesting an association between them and gastrointestinal endoparasite resistance. Further studies will help describe the effects of these markers accurately before implementing them in marker assisted selection. This study is the pioneer in describing such associations in goats.

A genome-wide scan for selection signatures in Nellore cattle.

Brazilian Nellore cattle (Bos indicus) have been selected for growth traits for over more than four decades. In recent years, reproductive and meat quality traits have become more important because of increasing consumption, exports and consumer demand. The identification of genome regions altered by artificial selection can potentially permit a better understanding of the biology of specific phenotypes that are useful for the development of tools designed to increase selection efficiency. Therefore, the aims of this study were to detect evidence of recent selection signatures in Nellore cattle using extended haplotype homozygosity methodology and BovineHD marker genotypes (>777,000 single nucleotide polymorphisms) as well as to identify corresponding genes underlying these signals. Thirty-one significant regions (P < 0.0001) of possible recent selection signatures were detected, and 19 of these overlapped quantitative trait loci related to reproductive traits, growth, feed efficiency, meat quality, fatty acid profiles and immunity. In addition, 545 genes were identified in regions harboring selection signatures. Within this group, 58 genes were associated with growth, muscle and adipose tissue metabolism, reproductive traits or the immune system. Using relative extended haplotype homozygosity to analyze high-density single nucleotide polymorphism marker data allowed for the identification of regions potentially under artificial selection pressure in the Nellore genome, which might be used to better understand autozygosity and the effects of selection on the Nellore genome.

Genome scan for meat quality traits in Nelore beef cattle.

Meat quality traits are economically important because they affect consumers' acceptance, which, in turn, influences the demand for beef. However, selection to improve meat quality is limited by the small numbers of animals on which meat tenderness can be evaluated due to the cost of performing shear force analysis and the resultant damage to the carcass. Genome wide-association studies for Warner-Bratzler shear force measured at different times of meat aging, backfat thickness, ribeye muscle area, scanning parameters [lightness, redness (a*), and yellowness] to ascertain color characteristics of meat and fat, water-holding capacity, cooking loss (CL), and muscle pH were conducted using genotype data from the Illumina BovineHD BeadChip array to identify quantitative trait loci (QTL) in all phenotyped Nelore cattle. Phenotype count for these animals ranged from 430 to 536 across traits. Meat quality traits in Nelore are controlled by numerous QTL of small effect, except for a small number of large-effect QTL identified for a*fat, CL, and pH. Genomic regions harboring these QTL and the pathways in which the genes from these regions act appear to differ from those identified in taurine cattle for meat quality traits. These results will guide future QTL mapping studies and the development of models for the prediction of genetic merit to implement genomic selection for meat quality in Nelore cattle.

Characterization of the exonic regions of the JY-1 gene in zebu cattle and buffaloes.

Protein JY-1 is an oocyte-specific protein that plays an important regulatory role in the granulosa cell layer and during the early embryo development stages. It is the first specific protein of maternal origin discovered in a single-ovulating species. In this study, the exon regions of the JY-1 gene were characterized by sequencing in 20 unrelated cattle (Bos taurus indicus) and 20 unrelated buffaloes (Bubalus bubalis). Eighteen polymorphisms were detected in cattle and 10 polymorphisms in buffaloes. Some of the polymorphisms were identified in codifying regions and caused amino acid changes. The insertion of a thymine was detected in the codifying region of exon 3 of the buffalo sequence when compared to the cattle one. This insertion causes a change in the codons frameshift from this point onwards, modifying the 19 terminal amino acids of the buffalo protein and creating a premature stop codon. This finding may explain reproductive differences between cattle and buffaloes in terms of follicle recruitment, embryo development and incidence of twin pregnancies.

Candidate genes for production traits in Nelore beef cattle.

Meat quality is an important trait for the beef industry. Backfat thickness, ribeye area, and shear force are traits measured late in life, and the investigation of molecular markers associated with these traits can help breeding programs. In cattle, some polymorphisms have been related to production traits. Thus, the purpose of this study was to assess the presence of polymorphisms in the candidate genes insulin-like growth factor 1 (IGF1), fatty acid-binding protein 4 (FABP4), and peroxisome proliferative active receptor gamma coactivator 1 A (PPARGC1A) and associate them with production traits in reference families of Nelore cattle. We used 270 steers descendent from 20 sires that were chosen to represent variability in this breed. The investigation of marker effects on the traits was performed using a mixed model under the restricted maximum likelihood method. A significant allele substitution effect was found for IGF1 and yearling weight (P ≤ 0.017). The mean allele substitution effect was 6.9 kg, with the 229 allele associated with reduced yearling weight in this Nelore population.

A SNP in ASAP1 gene is associated with meat quality and production traits in Nelore breed.

The ASAP1 gene is located in a QTL region for meat production traits and to access the role of the ASAP1 gene, the association between a SNP in this gene and production traits in beef cattle was studied. For this, about 270 steers of reference families of Nelore breed were used. The investigation of marker effects on the traits was performed using a mixed model under the restricted maximum likelihood method. Novel association of a SNP in the ASAP1 gene and shear force measured at 24 h post mortem (P≤0.0083) was described in this population of Nelore cattle. This polymorphism accounted for 1.13% of the total additive variance and 17.51% of total phenotypic variance of the trait, suggesting that this marker could be used in marker assisted selection.

Whole-genome analysis for backfat thickness in a tropically adapted, composite cattle breed from Brazil.

Backfat thickness affects the preservation of the beef carcass after slaughter and confers organoleptic characteristics assessed by the consumer. One of the breeding goals for Canchim, a tropically adapted breed, is to comprehensively increase fat thickness. Our goals were to identify genomic regions associated with backfat in Canchim populations and validate the association of single nucleotide polymorphisms (SNPs) overlapping previously identified QTL regions known to affect fat deposition. Fifteen animals with lower and 15 animals with higher residues for backfat, according to a linear model using the SAS GLM procedure, were selected from a population of 1171 animals and genotyped using the BovineSNP50 BeadChip. Initial analysis revealed more than 100 SNPs that discriminated the tails of phenotypic distribution. One extended region of association included the centromeric region of chromosome (Chr) 14. Because this region overlapped with QTL from previous reports, we developed SNP assays to interrogate two linkage disequilibrium blocks, one in the centromeric region and another in the middle region of Chr 14 to confirm the association. The analysis validated the presence of specific haplotypes affecting fat thickness.

Resistance of cattle of various genetic groups to the tick Rhipicephalus microplus and the relationship with coat traits.

This study evaluated the resistance of cattle of different genetic groups to the tick Rhipicephalus microplus and the relationship with traits of the animals' hair and coat. Cows of the Senepol×Nelore (SN), Angus×Nelore (AN) and Nelore (NX) genetic groups were submitted to four consecutive artificial infestations, at 14-day intervals, each one with approximately 20,000 tick larvae placed on the animals' lumbar region. From the 19th to 23rd day of each infestation five counts of the number of ticks were performed on each animal's left body side. The tick count data (TTC) were transformed into log(10) (n+1), and also into percentage of return (PR), where n is the total number of ticks counted at each infestation. Hair samples were collected 24h after the last infestation with flat-nosed pliers. Measures of the average hair length (HL), coat thickness (CT), number of hairs per cm(2) (NHCM2) and weight of the samples (SW) were obtained. Pearson's correlation coefficients were calculated within genetic group to measure association between PR and the hair and coat data. There was a significant difference among genetic groups for the number of ticks, with the AN group having higher counts than the SN and NX groups. For the hair and coat traits, the NX and SN groups had lower values of HL and SW than did the AN group. The SN genetic group had lower NHCM2 counts than the NX and AN groups. There were positive correlations between TTC and CT (P<0.05) and SW (P<0.05) in the SN group. No significant correlation was found for the AN genetic group (P>0.05).

Candidate gene region for control of rib eye area in Canchim beef cattle.

Investigation of molecular marker effects on production traits is essential to define marker assisted selection strategies in beef cattle. We looked for a possible association of molecular markers and backfat thickness (BFT) and rib eye area (REA) in Canchim (5/8 Charolais + 3/8 Zebu) and MA (offspring of Charolais bulls and 1/2 Canchim + 1/2 Zebu cows) animals raised exclusively on pasture. Traits were measured on 987 individuals from seven herds from two Brazilian States (São Paulo and Goiás), in March and April from 2005 to 2007, when animals were, on average, 19 months of age. Five microsatellite markers lying in QTL regions for BFT and REA (BMS490 and ETH10 on chromosome 5, INRA133 and ILSTS090 on chromosome 6, and BMS2142 on chromosome 19) were genotyped and association analyses were performed under an animal model using the restricted maximum likelihood method. After correction for multiple tests, a significant effect of microsatellite BMS490 on REA was observed, suggesting that at least one QTL affecting carcass traits in this region of the BTA5. No significant effect on BFT was observed for these markers.

mRNA profile of Nellore calves after primary infection with Haemonchus placei.

Haemonchus parasites are responsible for many losses in animal production. However, few studies are available, especially of zebu cattle. In this study, we investigated mRNA differences of immune response genes in naïve Nellore calves infected with Haemonchus placei, relating these differences to patterns of cellular infiltrate. Calves were infected with 15,000 H. placei L3 larvae and after 7 days lymph node and abomasum tissues were collected. IL-2, IL-4, IL-8, IL-12, IL-13, IFN-γ, MCP-1, lysozyme, pepsinogen and TNF-α genes were evaluated by qPCR. Mast cells, eosinophils and globular leukocytes were counted by abomasum histology. In the infected group, IL-4, IL-13 and TNF-α were up-regulated in the abomasal lymph node. In the abomasum, IL-13 increased and TNF-α was down-regulated (p<0.05). No differences were detected for mast cells and eosinophil counts in abomasal tissue (p>0.05). We conclude that for this infection time, there was Th2 polarization but that cellular infiltrate in abomasal tissue takes longer to develop.

Prospecting candidate SNPs for backfat in Canchim beef cattle.

Canchim is a composite cattle breed developed in Brazil for beef production. One of the breeding objectives is to increase fat deposition. QTLs for fat thickness and/or marbling have been reported on BTA4 and BTA14. The IGFBP3 and DDEF1 genes, mapped to BTA4 and BTA14, respectively, affect adipogenesis. We looked for SNPs in the IGFBP3 and DDEF1 genes that could be associated with backfat thickness in Canchim beef cattle. For SNP identification, sires with the highest accuracy were ranked according to expected breeding value for fat thickness; the 12 extremes (six sires with the highest and six with the lowest expected breeding value for the trait) were chosen. Six regions of the IGFBP3 and 14 regions of the DDEF1 were sequenced using the Sanger method. Nine SNPs were identified in IGFBP3 and 76 in the DDEF1. After an initial analysis, two SNPs were selected to be genotyped for the whole population; these were DDEF1g.279401A>G and IGFBP3c.4394T>C(Trp>Arg). We found a significant effect (P ≤ 0.05) of allele substitution on backfat thickness; however, the IGFBP3 SNP did not significantly affect this trait.

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.