Long-term reproductive effects in male rats prenatally exposed to sodium arsenite.

Arsenic is an environmental toxicant known to be a carcinogen and endocrine disruptor. Maternal exposure to arsenic has been associated with fetus malformation and reproductive disorders in male offspring. However, it is unclear the extent to which those effects remain during postnatal development and adulthood. Therefore, this study aimed to investigate the long-term effects of prenatal arsenic exposure on reproductive parameters of male offspring at peripubertal and adult periods. Pregnant female Wistar rats were exposed to 0 or 10 mg/L sodium arsenite in drinking water from gestational day 1 (GD 1) until GD 21 and male pups were analyzed at postnatal day 44 (PND 44) and PND 70. We observed that some reproductive parameters were affected differently by arsenic exposure at each age evaluated. The body and reproductive organs weights, as well as testicular and epididymal morphology were strongly affected in peripubertal animals and recovered at adult period. On the other hand, the antioxidant genes expression (SOD1, SOD2, CAT and GSTK1) and the endogenous antioxidant system were affected in the testes and epididymides from both peripubertal and adult rats. Finally, an impairment in daily sperm production and in sperm parameters was observed in adult animals. Taken together, our findings show that prenatal arsenic exposure affected reproductive parameters of peripubertal and adult male rats mainly due to oxidative stress. Collectively, those alterations may be affecting fertility potential of adult animals.

Fine Mapping of a Major Backfat QTL Reveals a Causal Regulatory Variant Affecting the CCND2 Gene.

Backfat is an important trait in pork production, and it has been included in the breeding objectives of genetic companies for decades. Although adipose tissue is a good energy storage, excessive fat results in reduced efficiency and economical losses. A large QTL for backfat thickness on chromosome 5 is still segregating in different commercial pig breeds. We fine mapped this QTL region using a genome-wide association analysis (GWAS) with 133,358 genotyped animals from five commercial populations (Landrace, Pietrain, Large White, Synthetic, and Duroc) imputed to the porcine 660K SNP chip. The lead SNP was located at 5:66103958 (G/A) within the third intron of the CCND2 gene, with the G allele associated with more backfat, while the A allele is associated with less backfat. We further phased the QTL region to discover a core haplotype of five SNPs associated with low backfat across three breeds. Linkage disequilibrium analysis using whole-genome sequence data revealed three candidate causal variants within intronic regions and downstream of the CCND2 gene, including the lead SNP. We evaluated the association of the lead SNP with the expression of the genes in the QTL region (including CCND2) in a large cohort of 100 crossbred samples, sequenced in four different tissues (lung, spleen, liver, muscle). Results show that the A allele increases the expression of CCND2 in an additive way in three out of four tissues. Our findings indicate that the causal variant for this QTL region is a regulatory variant within the third intron of the CCND2 gene affecting the expression of CCND2.

Transcription Landscape of the Early Developmental Biology in Pigs.

Since pre- and postnatal development are programmed during early prenatal life, studies addressing the complete transcriptional landscape during organogenesis are needed. Therefore, we aimed to disentangle differentially expressed (DE) genes between fetuses (at 35 days old) and embryos (at 25 days old) through RNA-sequencing analysis using the pig as model. In total, 1705 genes were DE, including the top DE IBSP, COL6A6, HBE1, HBZ, HBB, and NEUROD6 genes, which are associated with developmental transition from embryos to fetuses, such as ossification, skeletal muscle development, extracellular matrix organization, cardiovascular system, erythrocyte differentiation, and neuronal system. In pathway analysis, embryonic development highlighted those mainly related to morphogenic signaling and cell interactions, which are crucial for transcriptional control during the establishment of the main organs in early prenatal development, while pathways related to myogenesis, neuronal development, and cardiac and striated muscle contraction were enriched for fetal development, according to the greater complexity of organs and body structures at this developmental stage. Our findings provide an exploratory and informative transcriptional landscape of pig organogenesis, which might contribute to further studies addressing specific developmental events in pigs and in other mammals.

Applying an association weight matrix in weighted genomic prediction of boar taint compounds.

Biological information regarding markers and gene association may be used to attribute different weights for single nucleotide polymorphism (SNP) in genome-wide selection. Therefore, we aimed to evaluate the predictive ability and the bias of genomic prediction using models that allow SNP weighting in the genomic relationship matrix (G) building, with and without incorporating biological information to obtain the weights. Firstly, we performed a genome-wide association studies (GWAS) in data set containing single- (SL) or a multi-line (ML) pig population for androstenone, skatole and indole levels. Secondly, 1%, 2%, 5%, 10%, 30% and 50% of the markers explaining the highest proportions of the genetic variance for each trait were selected to build gene networks through the association weight matrix (AWM) approach. The number of edges in the network was computed and used to derive weights for G (AWM-WssGBLUP). The single-step GBLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used as standard scenarios. All scenarios presented predictive abilities different from zero; however, the great overlap in their confidences interval suggests no differences among scenarios. Most of scenarios of based on AWM provide overestimations for skatole in both SL and ML populations. On the other hand, the skatole and indole prediction were no biased in the ssGBLUP (S1) in both SL and ML populations. Most of scenarios based on AWM provide no biased predictions for indole in both SL and ML populations. In summary, using biological information through AWM matrix and gene networks to derive weights for genomic prediction resulted in no increase in predictive ability for boar taint compounds. In addition, this approach increased the number of analyses steps. Thus, we can conclude that ssGBLUP is most appropriate for the analysis of boar taint compounds in comparison with the weighted strategies used in the present work.

Single-step genome-wide association studies (GWAS) and post-GWAS analyses to identify genomic regions and candidate genes for milk yield in Brazilian Girolando cattle.

Milk production is economically important to the Brazilian agribusiness, and the majority of the country's milk production derives from Girolando (Gir × Holstein) cows. This study aimed to identify quantitative trait loci (QTL) and candidate genes associated with 305-d milk yield (305MY) in Girolando cattle. In addition, we investigated the SNP-specific variances for Holstein and Gir breeds of origin within the sequence of candidate genes. A single-step genomic BLUP procedure was used to identify QTL associated with 305MY, and the most likely candidate genes were identified through follow-up analyses. Genomic breeding values specific for Holstein and Gir were estimated in the Girolando animals using a model that uses breed-specific partial relationship matrices, which were converted to breed of origin SNP effects. Differences between breed of origin were evaluated by comparing estimated SNP variances between breeds. From 10 genome regions explaining most additive genetic variance for 305MY in Girolando cattle, 7 candidate genes were identified on chromosomes 1, 4, 6, and 26. Within the sequence of these 7 candidate genes, Gir breed of origin SNP alleles showed the highest genetic variance. These results indicated QTL regions that could be further explored in genomic selection panels and which may also help in understanding the gene mechanisms involved in milk production in the Girolando breed.

l-Arginine supplementation of gilts during early gestation modulates energy sensitive pathways in pig conceptuses.

Dietary l-arginine (ARG) supplementation has been studied as a nutritional strategy to improve reproductive performance of pregnant sows, since arginine is a conditionally essential amino acid. However, reports addressing the molecular mechanisms that mediate supplementation effects on embryos and fetuses development are still scarce. Therefore, we aimed to evaluate the effects of 1.0% ARG supplementation of commercial pregnant gilts on genes and proteins from energy metabolism and antioxidant defense pathways in embryos and fetuses. We also analyzed the global transcriptome profile of 25- and 35-day-old conceptuses. At Day 25, we observed a lower abundance of phospho-AMP-activated protein kinase (phospho-AMPK) protein and downregulation of oxidative phosphorylation system genes in ARG embryos. On the other hand, ARG fetuses showed greater expression of MLST8 and lower expression of MTOR genes, in addition to lower abundance of phospho-AMPK and phospho-mammalian target of rapamycin (phospho-mTOR) proteins. Transcriptome analysis at Day 35 did not present differentially expressed genes. For the antioxidant defense pathway, no differences were found between CON and ARG conceptuses, only trends. In general, supplementation of gilts with 1.0% ARG during early gestation affects energy sensitive pathways in 25- and 35-day conceptuses; however, no effects of supplementation were found on the antioxidative defense pathway in 25-day embryos.

Proteomic analysis reveals changes in energy metabolism of skeletal muscle in beef cattle supplemented with vitamin A.

BACKGROUND: Vitamin A has been reported as a factor influencing marbling deposition in meat from animals. Although the mechanisms by which vitamin A regulates lipid metabolism in mature adipocytes are already well-established, information regarding molecular mechanisms underlying the effects of vitamin A on the regulation of intramuscular fat deposition in beef cattle still remains limited. The present study aimed to assess the molecular mechanisms involved in the intramuscular fat deposition in beef cattle supplemented with vitamin A during the fattening phase using a proteomic approach. RESULTS: Vitamin A supplementation during the fattening phase decreased intramuscular fat deposition in beef cattle. Proteome and phospho-proteome analysis together with biological and networking analysis of the protein differentially abundant between treatments indicated that Vitamin A supplementation affects the overall energy metabolism of skeletal muscle, impairing lipid biosynthesis in skeletal muscle. CONCLUSION: Vitamin A supplementation at fattening phase impairs intramuscular fat deposition in beef cattle likely by changing the energy metabolism of skeletal muscle. The interaction of retinoic acid and heat shock 70-kDa protein may play a pivotal role in intramuscular fat deposition as a consequence of vitamin A supplementation by impairing de novo fatty acid synthesis as a result of a possible decrease in insulin sensitivity in the skeletal muscle. © 2020 Society of Chemical Industry.

Sex Determination Using RNA-Sequencing Analyses in Early Prenatal Pig Development.

Sexual dimorphism is a relevant factor in animal science, since it can affect the gene expression of economically important traits. Eventually, the interest in the prenatal phase in a transcriptome study may not comprise the period of development in which male and female conceptuses are phenotypically divergent. Therefore, it would be interesting if sex differentiation could be performed using transcriptome data, with no need for extra techniques. In this study, the sex of pig conceptuses (embryos at 25 days-old and fetuses at 35 days-old) was determined by reads counts per million (CPM) of Y chromosome-linked genes that were discrepant among samples. Thus, ten genes were used: DDX3Y, KDM5D, ZFY, EIF2S3Y, EIF1AY, LOC110255320, LOC110257894, LOC396706, LOC100625207, and LOC110255257. Conceptuses that presented reads CPM sum for these genes (ΣCPMchrY) greater than 400 were classified as males and those with ΣCPMchrY below 2 were classified as females. It was demonstrated that the sex identification can be performed at early stages of pig development from RNA-sequencing analysis of genes mapped on Y chromosome. Additionally, these results reinforce that sex determination is a mechanism conserved across mammals, highlighting the importance of using pigs as an animal model to study sex determination during human prenatal development.

Genome-wide association studies for heat stress response in Bos taurus × Bos indicus crossbred cattle.

Heat stress is an important issue in the global dairy industry. In tropical areas, an alternative to overcome heat stress is the use of crossbred animals or synthetic breeds, such as the Girolando. In this study, we performed a genome-wide association study (GWAS) and post-GWAS analyses for heat stress in an experimental Gir × Holstein F2 population. Rectal temperature (RT) was measured in heat-stressed F2 animals, and the variation between 2 consecutive RT measurements (ΔRT) was used as the dependent variable. Illumina BovineSNP50v1 BeadChip (Illumina Inc., San Diego, CA) and single-SNP approach were used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene-transcription factor (TF) networks, generated from enriched TF. The breed origin of marker alleles in the F2 population was assigned using the breed of origin of alleles (BOA) approach. Heritability and repeatability estimates (± standard error) for ΔRT were 0.13 ± 0.08 and 0.29 ± 0.06, respectively. Association analysis revealed 6 SNP significantly associated with ΔRT. Genes involved with biological processes in response to heat stress effects (LIF, OSM, TXNRD2, and DGCR8) were identified as putative candidate genes. After performing the BOA approach, the 10% of F2 animals with the lowest breeding values for ΔRT were classified as low-ΔRT, and the 10% with the highest breeding values for ΔRT were classified as high-ΔRT. On average, 49.4% of low-ΔRT animals had 2 alleles from the Holstein breed (HH), and 39% had both alleles from the Gir breed (GG). In high-ΔRT animals, the average proportion of animals for HH and GG were 1.4 and 50.2%, respectively. This study allowed the identification of candidate genes for ΔRT in Gir × Holstein crossbred animals. According to the BOA approach, Holstein breed alleles could be associated with better response to heat stress effects, which could be explained by the fact that Holstein animals are more affected by heat stress than Gir animals and thus require a genetic architecture to defend the body from the deleterious effects of heat stress. Future studies can provide further knowledge to uncover the genetic architecture underlying heat stress in crossbred cattle.

Differentially expressed mRNAs, proteins and miRNAs associated to energy metabolism in skeletal muscle of beef cattle identified for low and high residual feed intake.

BACKGROUND: Feed efficiency is one of the most important parameters that affect beef production costs. The energy metabolism of skeletal muscle greatly contributes to variations in feed efficiency. However, information regarding differences in proteins involved in the energy metabolism of the skeletal muscle in beef cattle divergently identified for feed efficiency is scarce. In this study, we aimed to investigate energy metabolism of skeletal muscle of Nellore beef cattle, identified for low and high residual feed intake using a proteomics approach. We further assessed the expression of candidate microRNAs as a one of the possible mechanisms controlling the biosynthesis of the proteins involved in energy metabolism that were differentially abundant between high and low residual feed intake animals. RESULTS: A greater abundance of 14-3-3 protein epsilon (P = 0.01) was observed in skeletal muscle of residual feed intake (RFI) high animals (RFI-High). Conversely, a greater abundance of Heat Shock Protein Beta 1 (P < 0.01) was observed in the skeletal muscle of RFI-Low cattle. A greater mRNA expression of YWHAE, which encodes the 14-3-3 protein epsilon, was also observed in the skeletal muscle of RFI-High animals (P = 0.01). A lower mRNA expression of HSPB1, which encodes the Heat Shock Protein Beta 1, was observed in the skeletal muscle of RFI-High animals (P = 0.01). The miR-665 was identified as a potential regulator of the 14-3-3 protein epsilon, and its expression was greater in RFI-Low animals (P < .001). A greater expression of miR-34a (P = 0.01) and miR-2899 (P < .001) was observed in the skeletal muscle of RFI-High animals, as both miRNAs were identified as potential regulators of HSPB1 expression. CONCLUSION: Our results show that Nellore cattle divergently identified for feed efficiency by RFI present changes in the abundance of proteins involved in energy expenditure in skeletal muscle. Moreover, our data point towards that miR-665, miR34a and miR-2899 are likely involved in controlling both 14-3-3 epsilon and HSPB1 proteins identified as differentially abundant in the skeletal muscle of RFI-High and RFI-Low Nellore cattle.

In vitro evaluation of cryopreserved bovine sperm and its relation to field fertility in fixed-time artificial insemination.

This study aimed to assess characteristics of bovine cryopreserved sperm and evaluate its relation to field fertility in fixed-time artificial insemination (FTAI). Semen samples of 16 bulls were used to inseminate 811 Nellore cows, and four of these bulls were also used to inseminate 101 Nellore heifers. Samples of the same ejaculate used for FTAI from each bull were analysed in the laboratory after thawing. Sperm motility and vigour were subjectively assessed by light microscope, and integrity of the plasma and acrosome membranes, and H2 O2 production were evaluated by flow cytometer. Relation among sperm characteristics and pregnancy rate of cows and heifers were evaluated by univariate and multivariate logistic regression. Subjective sperm motility and vigour did not affect the probability of pregnancy in cows or heifers. In univariate analysis for pregnancy in cows, sperm traits related to acrosome injury positively affected probability of pregnancy mainly when associated with plasma membrane integrity; H2 O2 production seems to be less important than plasma membrane integrity in affecting probability of pregnancy. In multivariate analysis, sperm traits related to injured acrosome positively affected probability of cow and heifer pregnancies while intact acrosome was negatively related to cow pregnancy. Intact plasma membrane and high H2 O2 production were positively related to cow pregnancy but negatively related to heifer pregnancy. Results suggest that a capacitation-like status of the acrosome may benefit probability of pregnancy in cows.

Effects of alleles in crossbred pigs estimated for genomic prediction depend on their breed-of-origin.

BACKGROUND: This study investigated if the allele effect of a given single nucleotide polymorphism (SNP) for crossbred performance in pigs estimated in a genomic prediction model differs depending on its breed-of-origin, and how these are related to estimated effects for purebred performance. RESULTS: SNP-allele substitution effects were estimated for a commonly used SNP panel using a genomic best linear unbiased prediction model with breed-specific partial relationship matrices. Estimated breeding values for purebred and crossbred performance were converted to SNP-allele effects by breed-of-origin. Differences between purebred and crossbred, and between breeds-of-origin were evaluated by comparing percentage of variance explained by genomic regions for back fat thickness (BF), average daily gain (ADG), and residual feed intake (RFI). From ten regions explaining most additive genetic variance for crossbred performance, 1 to 5 regions also appeared in the top ten for purebred performance. The proportion of genetic variance explained by a genomic region and the estimated effect of a haplotype in such a region were different depending upon the breed-of-origin. To illustrate underlying mechanisms, we evaluated the estimated effects across breeds-of-origin for haplotypes associated to the melanocortin 4 receptor (MC4R) gene, and for the MC4Rsnp itself which is a missense mutation with a known effect on BF and ADG. Although estimated allele substitution effects of the MC4Rsnp mutation were very similar across breeds, explained genetic variance of haplotypes associated to the MC4R gene using a SNP panel that does not include the mutation, was considerably lower in one of the breeds where the allele frequency of the mutation was the lowest. CONCLUSIONS: Similar regions explaining similar additive genetic variance were observed across purebred and crossbred performance. Moreover, there was some overlap across breeds-of-origin between regions that explained relatively large proportions of genetic variance for crossbred performance; albeit that the actual proportion of variance deviated across breeds-of-origin. Results based on a missense mutation in MC4R confirmed that even if a causal locus has similar effects across breeds-of-origin, estimated effects and explained variance in its region using a commonly used SNP panel can strongly depend on the allele frequency of the underlying causal mutation.

Genome-wide association studies for tick resistance in Bos taurus × Bos indicus crossbred cattle: A deeper look into this intricate mechanism.

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F2 population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F2 population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.

Weighted single-step GWAS and gene network analysis reveal new candidate genes for semen traits in pigs.

BACKGROUND: In recent years, there has been increased interest in the study of the molecular processes that affect semen traits. In this study, our aim was to identify quantitative trait loci (QTL) regions associated with four semen traits (motility, progressive motility, number of sperm cells per ejaculate and total morphological defects) in two commercial pig lines (L1: Large White type and L2: Landrace type). Since the number of animals with both phenotypes and genotypes was relatively small in our dataset, we conducted a weighted single-step genome-wide association study, which also allows unequal variances for single nucleotide polymorphisms. In addition, our aim was also to identify candidate genes within QTL regions that explained the highest proportions of genetic variance. Subsequently, we performed gene network analyses to investigate the biological processes shared by genes that were identified for the same semen traits across lines. RESULTS: We identified QTL regions that explained up to 10.8% of the genetic variance of the semen traits on 12 chromosomes in L1 and 11 chromosomes in L2. Sixteen QTL regions in L1 and six QTL regions in L2 were associated with two or more traits within the population. Candidate genes SCN8A, PTGS2, PLA2G4A, DNAI2, IQCG and LOC102167830 were identified in L1 and NME5, AZIN2, SPATA7, METTL3 and HPGDS in L2. No regions overlapped between these two lines. However, the gene network analysis for progressive motility revealed two genes in L1 (PLA2G4A and PTGS2) and one gene in L2 (HPGDS) that were involved in two biological processes i.e. eicosanoid biosynthesis and arachidonic acid metabolism. PTGS2 and HPGDS were also involved in the cyclooxygenase pathway. CONCLUSIONS: We identified several QTL regions associated with semen traits in two pig lines, which confirms the assumption of a complex genetic determinism for these traits. A large part of the genetic variance of the semen traits under study was explained by different genes in the two evaluated lines. Nevertheless, the gene network analysis revealed candidate genes that are involved in shared biological pathways that occur in mammalian testes, in both lines.

Morphological and molecular differences in corpus luteum of pregnant sows from divergent genetic groups.

Comprehending mechanisms controlling corpus luteum (CL) angiogenesis and apoptosis in pregnant sows is essential to understand the physiological role of these processes in CL function, progesterone production and consequently in conceptus development and prenatal mortality. CL from 54 sows from two genetic groups, a commercial line (COM) and the local Piau breed (LPB), were obtained for gene expression (n = 3 COM; n = 6 LPB), histological and protein analysis (n = 3 COM; n = 3 LPB), divided in six gestational ages (seven, 15, 30, 45, 60 and 90 days). We observed differences between gestational ages in CL morphology, in which the average number of blood vessels/capillaries at 90-days was greater than at the seventh day by Tukey test. RT-qPCR analysis revealed that apoptotic genes (BAX, BCL2 and CASP3) were differentially expressed between genetic groups and gestational ages in each group. Angiogenesis genes also presented differences between genetic groups (ANGPT1) and gestational ages (MMP9, VEGFA and ANGPT1). No differences in protein abundance of steroidogenic enzymes (CYP11A1 and HSD3B1) were observed. Our findings indicate that despite the differences in gene expression, differences in corpus luteum vascularization were observed only across gestational ages, with no dissimilarities between genetic groups.

Gene networks for total number born in pigs across divergent environments.

For reproductive traits such as total number born (TNB), variance due to different environments is highly relevant in animal breeding. In this study, we aimed to perform a gene-network analysis for TNB in pigs across different environments using genomic reaction norm models. Thus, based on relevant single-nucleotide polymorphisms and linkage disequilibrium blocks across environments obtained from GWAS, different sets of candidate genes having biological roles linked to TNB were identified. Network analysis across environment levels resulted in gene interactions consistent with known mammal's fertility biology, captured relevant transcription factors for TNB biology and pointing out different sets of candidate genes for TNB in different environments. These findings may have important implication for animal production, as optimal breeding may vary depending on later environments. Based on these results, genomic diversity was identified and inferred across environments highlighting differential genetic control in each scenario.

Revealing new candidate genes for reproductive traits in pigs: combining Bayesian GWAS and functional pathways.

BACKGROUND: Reproductive traits such as number of stillborn piglets (SB) and number of teats (NT) have been evaluated in many genome-wide association studies (GWAS). Most of these GWAS were performed under the assumption that these traits were normally distributed. However, both SB and NT are discrete (e.g. count) variables. Therefore, it is necessary to test for better fit of other appropriate statistical models based on discrete distributions. In addition, although many GWAS have been performed, the biological meaning of the identified candidate genes, as well as their functional relationships still need to be better understood. Here, we performed and tested a Bayesian treatment of a GWAS model assuming a Poisson distribution for SB and NT in a commercial pig line. To explore the biological role of the genes that underlie SB and NT and identify the most likely candidate genes, we used the most significant single nucleotide polymorphisms (SNPs), to collect related genes and generated gene-transcription factor (TF) networks. RESULTS: Comparisons of the Poisson and Gaussian distributions showed that the Poisson model was appropriate for SB, while the Gaussian was appropriate for NT. The fitted GWAS models indicated 18 and 65 significant SNPs with one and nine quantitative trait locus (QTL) regions within which 18 and 57 related genes were identified for SB and NT, respectively. Based on the related TF, we selected the most representative TF for each trait and constructed a gene-TF network of gene-gene interactions and identified new candidate genes. CONCLUSIONS: Our comparative analyses showed that the Poisson model presented the best fit for SB. Thus, to increase the accuracy of GWAS, counting models should be considered for this kind of trait. We identified multiple candidate genes (e.g. PTP4A2, NPHP1, and CYP24A1 for SB and YLPM1, SYNDIG1L, TGFB3, and VRTN for NT) and TF (e.g. NF-κB and KLF4 for SB and SOX9 and ELF5 for NT), which were consistent with known newborn survival traits (e.g. congenital heart disease in fetuses and kidney diseases and diabetes in the mother) and mammary gland biology (e.g. mammary gland development and body length).

Molecular Factors Underlying the Deposition of Intramuscular Fat and Collagen in Skeletal Muscle of Nellore and Angus Cattle.

Studies have shown that intramuscular adipogenesis and fibrogenesis may concomitantly occur in skeletal muscle of beef cattle. Thus, we hypothesized that the discrepancy of intramuscular fat content in beef from Nellore and Angus was associated with differences in intramuscular adipogenesis and fibrogenesis during the finishing phase. To test our hypothesis, longissimus muscle samples of Nellore (n = 6; BW = 372.5 ± 37.3 kg) and Angus (n = 6; BW = 382.8 ± 23.9 kg) cattle were collected for analysis of gene and protein expression, and quantification of intramuscular fat and collagen. Least-squares means were estimated for the effect of Breed and differences were considered at P ≤ 0.05. A greater intramuscular fat content was observed in skeletal muscle of Angus compared to Nellore cattle (P≤0.05). No differences were observed for mRNA expression of lipogenic and lipolytic markers ACC, FAS, FABP4, SERBP-1, CPT-2, LPL, and ACOX (P > 0.05) in skeletal muscle of Nellore and Angus cattle. Similarly, no differences were observed in mRNA expression of adipogenic markers Zfp423, PPARγ, and C/EBPα (P>0.05) However, a greater PPARγ protein content was observed in skeletal muscle of Angus compared to Nellore cattle (P≤0.05). A greater abundance of adipo/fibrogenic cells, evaluated by the PDGFRα content, was observed in skeletal muscle of Angus than Nellore cattle (P≤0.05). No differences in fibrogenesis were observed in skeletal muscle of Angus and Nellore cattle, which is in accordance with the lack of differences in intramuscular collagen content in beef from both breeds (P>0.05). These findings demonstrate that difference in intramuscular fat content is associated with a slightly enhanced adipogenesis in skeletal muscle of Angus compared to Nellore cattle, while no difference in fibrogenesis.

Identification of Suitable Reference Genes for Real Time Quantitative Polymerase Chain Reaction Assays on Pectoralis major Muscle in Chicken (Gallus gallus ).

Thirteen reference genes were investigated to determine their stability to be used as a housekeeping in gene expression studies in skeletal muscle of chickens. Five different algorithms were used for ranking of reference genes and results suggested that individual rankings of the genes differed among them. The stability of the expression of reference genes were validated using samples obtained from the Pectoralis major muscle in chicken. Samples were obtained from chickens in different development periods post hatch and under different nutritional diets. For gene expression calculation the ΔΔCt approach was applied to compare relative expression of pairs of genes within each of 52 samples when normalized to mitochondrially encoded cytochrome c oxidase II (MT-CO2) target gene. Our findings showed that hydroxymethylbilane synthase (HMBS) and hypoxanthine phosphoribosyl transferase 1 (HPRT1) are the most stable reference genes while transferrin receptor (TFRC) and beta-2-microglobulin (B2M) ranked as the least stable genes in the Pectoralis major muscle of chickens. Moreover, our results revealed that HMBS and HPRT1 gene expression did not change due to dietary variations and thus it is recommended for accurate normalization of RT-qPCR data in chicken Pectoralis major muscle.

Linkage disequilibrium patterns and persistence of phase in purebred and crossbred pig (Sus scrofa) populations.

BACKGROUND: Genomic selection and genomic wide association studies are widely used methods that aim to exploit the linkage disequilibrium (LD) between markers and quantitative trait loci (QTL). Securing a sufficiently large set of genotypes and phenotypes can be a limiting factor that may be overcome by combining data from multiple breeds or using crossbred information. However, the estimated effect of a marker in one breed or a crossbred can only be useful for the selection of animals in another breed if there is a correspondence of the phase between the marker and the QTL across breeds. Using data of five pure pig (Sus scrofa) lines (SL1, SL2, SL3, DL1, DL2), one F1 cross (DLF1) and two commercial finishing crosses (TER1 and TER2), the objectives of this study were: (i) to compare the equality of LD decay curves of different pig populations; and (ii) to evaluate the persistence of the LD phase across lines or final crosses. RESULTS: Almost all of the lines presented different extents of LD, except for the SL2 and DL3, both of which exhibited the same extent of LD. Similar levels of LD over large distances were found in crossbred and pure lines. The crossbred animals (DLF1, TER1 and TER2) presented a high persistence of phase with their parental lines, suggesting that the available porcine single nucleotide polymorphism (SNP) chip should be dense enough to include markers that have the same LD phase with QTL across crossbred and parental pure lines. The persistence of phase across pure lines varied considerably between the different line comparisons; however, correlations were above 0.8 for all line comparisons when marker distances were smaller than 50 kb. CONCLUSIONS: This study showed that crossbred populations could be very useful as a reference for the selection of pure lines by means of the available SNP chip panel. Here, we also pinpoint pure lines that could be combined in a multiline training population. However, if multiline reference populations are used for genomic selection, the required density of SNP panels should be higher compared with a single breed reference population.