Differential expression and co-expression gene networks reveal candidate biomarkers of boar taint in non-castrated pigs

Differential Regulation of Male-Hormones-Related Enhancers Revealed by Chromatin Accessibility and Transcriptional Profiles in Pig Liver

Surgical castration can effectively avoid boar taint and improve pork quality by removing the synthesis of androstenone in the testis, thereby reducing its deposition in adipose tissue. The expression of genes involved in testis-derived hormone metabolism was altered following surgical castration, but the upstream regulatory factors and underlying mechanism remain unclear. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics in liver tissue of castrated and intact full-sibling Yorkshire pigs. First, we identified 897 differentially expressed genes and 6864 differential accessible regions (DARs) using RNA- and ATAC-seq. By integrating the RNA- and ATAC-seq results, 227 genes were identified, and a significant positive correlation was revealed between differential gene expression and the ATAC-seq signal. We constructed a transcription factor regulatory network after motif analysis of DARs and identified a candidate transcription factor (TF) SP1 that targeted the HSD3B1 gene, which was responsible for the metabolism of androstenone. Subsequently, we annotated DARs by incorporating H3K27ac ChIP-seq data, marking 2234 typical enhancers and 245 super enhancers involved in the regulation of all testis-derived hormones. Among these, four typical enhancers associated with HSD3B1 were identified. Furthermore, an in-depth investigation was conducted on the androstenone-related enhancers, and an androstenone-related mutation was identified in a newfound candidatetypical enhancer (andEN) with dual-luciferase assays. These findings provide further insights into how enhancers function as links between phenotypic and non-coding area variations. The discovery of upstream TF and enhancers of HSD3B1 contributes to understanding the regulatory networks of androstenone metabolism and provides an important foundation for improving pork quality.

GEO dataset mining analysis reveals novel Staphylococcus aureus virulence gene regulatory networks and diagnostic targets in mice

Staphylococcus (S.) aureus infection is a serious, worldwide health concern, particularly in many communities and hospitals. Understanding the S. aureus pathogenetic regulatory network will provide significant insights into diagnostic target screening to improve clinical treatment of diseases caused by S. aureus. We screened differentially expressed genes between normal mice and S. aureus-infected mice. We used the Gene Expression Omnibus (GEO) DataSets database for functional analysis (GO-analysis) and the DAVID and KEGG databases for signaling pathway analyses. We next integrated the gene and pathway analyses with Transcriptional Regulatory Element Database (TRED) to build an antimicrobial resistance gene regulatory network of S. aureus. We performed association analysis of network genes and diseases using DAVID online annotation tools. We identified a total of 437 virulence genes and 15 transcription factors (TFs), as well as 444 corresponding target genes, in the S. aureus TF regulatory network. We screened seven key network nodes (Met, Mmp13, Il12b, Il4, Tnf, Ptgs2, and Ctsl), four key transcription factors (Jun, C3, Spil, and Il6) and an important signaling pathway (TNF). We hypothesized that the cytokine activity and growth factor activity of S. aureus are combinatorically cross-regulated by Met, Mmp13, Il12b, Il4, Tnf, Ptgs2, and Ctsl genes, the TFs Jun, C3, Spi1, and Il6, as well as the immune response, cellular response to lipopolysaccharide, and inflammatory response. Our study provides information and reference values for the molecular understanding of the S. aureus pathogenetic gene regulatory network.

Nuclear Receptor Pathways Mediating the Development of Boar Taint

The nuclear receptors PXR, CAR, and FXR are activated by various ligands and function as transcription factors to control the expression of genes that regulate the synthesis and metabolism of androstenone and skatole. These compounds are produced in entire male pigs and accumulate in the fat to cause the development of a meat quality issue known as boar taint. The extent of this accumulation is influenced by the synthesis and hepatic clearance of androstenone and skatole. For this reason, PXR, CAR, and FXR-mediated signaling pathways have garnered interest as potential targets for specialized treatments designed to reduce the development of boar taint. Recent research has also identified several metabolites produced by gut microbes that act as ligands for these nuclear receptors (e.g., tryptophan metabolites, short-chain fatty acids, bile acids); however, the connection between the gut microbiome and boar taint development is not clear. In this review, we describe the nuclear receptor signaling pathways that regulate the synthesis and metabolism of boar taint compounds and outline the genes involved. We also discuss several microbial-derived metabolites and dietary additives that are known or suspected nuclear receptor ligands and suggest how these compounds could be used to develop novel treatments for boar taint.

Novel role of COX6c in the regulation of oxidative phosphorylation and diseases

Cytochrome c oxidase subunit VIc (COX6c) is one of the most important subunits of the terminal enzyme of the respiratory chain in mitochondria. Numerous studies have demonstrated that COX6c plays a critical role in the regulation of oxidative phosphorylation (OXPHOS) and energy production. The release of COX6c from the mitochondria may be a hallmark of the intrinsic apoptosis pathway. Moreover, The changes in COX6c expression are widespread in a variety of diseases and can be chosen as a potential biomarker for diagnosis and treatment. In light of its exclusive effects, we present the elaborate roles that COX6c plays in various diseases. In this review, we first introduced basic knowledge regarding COX6c and its functions in the OXPHOS and apoptosis pathways. Subsequently, we described the regulation of COX6c expression and activity in both positive and negative ways. Furthermore, we summarized the elaborate roles that COX6c plays in various diseases, including cardiovascular disease, kidney disease, brain injury, skeletal muscle injury, and tumors. This review highlights recent advances and provides a comprehensive summary of COX6c in the regulation of OXPHOS in multiple diseases and may be helpful for drug design and the prediction, diagnosis, treatment, and prognosis of diseases.

Mechanism of Phellodendron and Anemarrhena Drug Pair on the Treatment of Liver Cancer Based on Network Pharmacology and Bioinformatics

Background

This study aims to explore the key targets and signaling pathways of the traditional Chinese medicine Phellodendron and Anemarrhena drug pair (PADP) for the treatment of liver cancer.

Methods

Firstly, bioinformatics technology was used to analyze GSE62232 gene chip to obtain the differential genes of liver cancer. A network pharmacology technology was used to find the active components of PADP and their targets. Secondly, the differential genes were imported into STRING database to draw a PPI network, and network topology structure map combined with Cytoscape software. And the R language was used to identify differential gene targets and pathways through GO and KEGG pathway enrichment analysis. In addition, AutoDock Vina was used for molecular docking of core targets and core compounds. Moreover, GEPIA online analysis tool was used to perform survival analysis of the core target genes. Finally, RT-PCR was used to verify the changes of key target genes. CCK−8 assay was performed to detect cell proliferation. Flow cytometry was performed to detect the cell cycle and apoptotic. Transwell invasion assay was performed to detect cell invasion.

Results

Firstly, a total of 21,654 genes were obtained. After screening, 1019 differential genes were obtained, including 614 down-regulated genes and 405 up-regulated genes. Furthermore, after screening by ADME standards, 52 active ingredients were obtained, of which 37 were Phellodendron and 15 were Anemarrhena. And a total of 36 differential genes have been identified, including 13 up-regulated genes and 23 down-regulated genes. Moreover, through enrichment analysis, we found that PADP may treat liver cancer through multiple channels and multiple pathways including the p53 signaling pathway, IL-17 signaling pathway, TNF signaling pathway, Toll-like receptor signaling pathway and so on. Secondly, the molecular docking results showed that there was certain affinity between the core compounds and core target genes. In addition, GEPIA online analysis showed that ESR1, AR, CCNB1, CDK1, AKR1C3 and CCNA2 might become potential target genes for the survival and prognosis of PADP for the treatment of liver cancer. Finally, it was found that PADP could up regulate genes ESR1 and AR, down regulate genes CCNB1, CDK1, AKR1C3, and CCNA2. PADP could promote the apoptosis of liver cancer cells, shorten the cell cycle, and inhibit the proliferation and invasion of liver cancer cells.

Conclusion

PADP may treat liver cancer through multiple targets, multiple channels, and multiple pathways, thereby suppressing cancer cells and improving the living quality of patients.

Gene Co-Expression Network Analysis Identifies Vitamin D-Associated Gene Modules in Adult Normal Rectal Epithelium Following Supplementation

Colorectal cancer (CRC) is a common, multifactorial disease. While observational studies have identified an association between lower vitamin D and higher CRC risk, supplementation trials have been inconclusive and the mechanisms by which vitamin D may modulate CRC risk are not well understood. We sought to perform a weighted gene co-expression network analysis (WGCNA) to identify modules present after vitamin D supplementation (when plasma vitamin D level was sufficient) which were absent before supplementation, and then to identify influential genes in those modules. The transcriptome from normal rectal mucosa biopsies of 49 individuals free from CRC were assessed before and after 12 weeks of 3200IU/day vitamin D (Fultium-D3) supplementation using paired-end total RNAseq. While the effects on expression patterns following vitamin D supplementation were subtle, WGCNA identified highly correlated genes forming gene modules. Four of the 17 modules identified in the post-vitamin D network were not preserved in the pre-vitamin D network, shedding new light on the biochemical impact of supplementation. These modules were enriched for GO terms related to the immune system, hormone metabolism, cell growth and RNA metabolism. Across the four treatment-associated modules, 51 hub genes were identified, with enrichment of 40 different transcription factor motifs in promoter regions of those genes, including VDR:RXR. Six of the hub genes were nominally differentially expressed in studies of vitamin D effects on adult normal mucosa organoids: LCN2, HLA-C, AIF1L, PTPRU, PDE4B and IFI6. By taking a gene-correlation network approach, we have described vitamin D induced changes to gene modules in normal human rectal epithelium in vivo, the target tissue from which CRC develops.

SEaCorAl: Identifying and contrasting the regulation-correlation bias in RNA-Seq paired expression data of patient groups

The Cancer Genome Atlas database offers the possibility of analyzing genome-wide expression RNA-Seq cancer data using paired counts, that is, studies where expression data are collected in pairs of normal and cancer cells, by taking samples from the same individual. Correlation of gene expression profiles is the most common analysis to study co-expression groups, which is used to find biological interpretation of -omics big data. The aim of the paper is threefold: firstly we show for the first time, the presence of a "regulation-correlation bias" in RNA-Seq paired expression data, that is an artifactual link between the expression status (up- or down-regulation) of a gene pair and the sign of the corresponding correlation coefficient. Secondly, we provide a statistical model able to theoretically explain the reasons for the presence of such a bias. Thirdly, we present a bias-removal algorithm, called SEaCorAl, able to effectively reduce bias effects and improve the biological significance of correlation analysis. Validation of the SEaCorAl algorithm is performed by showing a significant increase in the ability to detect biologically meaningful associations of positive correlations and a significant increase of the modularity of the resulting unbiased correlation network.

Differential expression and clinical significance of COX6C in human diseases

Mitochondria, independent double-membrane organelles, are intracellular power plants that feed most eukaryotic cells with the ATP produced via the oxidative phosphorylation (OXPHOS). Consistently, cytochrome c oxidase (COX) catalyzes the electron transfer chain's final step. Electrons are transferred from reduced cytochrome c to molecular oxygen and play an indispensable role in oxidative phosphorylation of cells. Cytochrome c oxidase subunit 6c (COX6C) is encoded by the nuclear genome in the ribosome after translation and is transported to mitochondria via different pathways, and eventually forms the COX complex. In recent years, many studies have shown the abnormal level of COX6C in familial hypercholesterolemia, chronic kidney disease, diabetes, breast cancer, prostate cancer, uterine leiomyoma, follicular thyroid cancer, melanoma tissues, and other conditions. Its underlying mechanism may be related to the cellular oxidative phosphorylation pathway in tissue injury disease. Here reviews the varied function of COX6C in non-tumor and tumor diseases.

Recent genetic advances on boar taint reduction as an alternative to castration: a review

Boar taint is an unpleasant odor in male pig meat, mainly caused by androstenone, skatole, and indole, which are deposited in the fat tissue. Piglet castration is the most common practice to prevent boar taint. However, castration is likely to be banished in a few years due to animal welfare concerns. Alternatives to castration, such as genetic selection, have been assessed. Androstenone and skatole have moderate to high heritability, which makes it feasible to select against these compounds. This review presents the latest results obtained on genetic selection against boar taint, on correlation with other traits, on differences in breeds, and on candidate genes related to boar taint. QTLs for androstenone and skatole have been reported mainly on chromosomes 6, 7, and 14. These chromosomes were reported to contain genes responsible for synthesis and degradation of androstenone and skatole. A myriad of work has been done to find markers or genes that can be used to select animals with lower boar taint. The selection against boar taint could decrease performance of some reproduction traits. However, a favorable response on production traits has been observed by selecting against boar taint. Selection results have shown that it is possible to reduce boar taint in few generations. In addition, modifications in diet and environment conditions could be associated with genetic selection to reduce boar taint. Nevertheless, costs to measure and select against boar taint should be rewarded with incentives from the market; otherwise, it would be difficult to implement genetic selection.

(Differential) Co-Expression Analysis of Gene Expression: A Survey of Best Practices

Analysis of gene expression data is widely used in transcriptomic studies to understand functions of molecules inside a cell and interactions among molecules. Differential co-expression analysis studies diseases and phenotypic variations by finding modules of genes whose co-expression patterns vary across conditions. We review the best practices in gene expression data analysis in terms of analysis of (differential) co-expression, co-expression network, differential networking, and differential connectivity considering both microarray and RNA-seq data along with comparisons. We highlight hurdles in RNA-seq data analysis using methods developed for microarrays. We include discussion of necessary tools for gene expression analysis throughout the paper. In addition, we shed light on scRNA-seq data analysis by including preprocessing and scRNA-seq in co-expression analysis along with useful tools specific to scRNA-seq. To get insights, biological interpretation and functional profiling is included. Finally, we provide guidelines for the analyst, along with research issues and challenges which should be addressed.

Characterization of Global DNA Methylation in Different Gene Regions Reveals Candidate Biomarkers in Pigs with High and Low Levels of Boar Taint

DNA methylation of different gene components, including different exons and introns, or different lengths of exons and introns is associated with differences in gene expression. To investigate the methylation of porcine gene components associated with the boar taint (BT) trait, this study used reduced representation bisulfite sequencing (RRBS) data from nine porcine testis samples in three BT groups (low, medium and high BT). The results showed that the methylation levels of the first exons and first introns were lower than those of the other exons and introns. The first exons/introns of CpG island regions had even lower levels of methylation. A total of 123 differentially methylated promoters (DMPs), 194 differentially methylated exons (DMEs) and 402 differentially methylated introns (DMIs) were identified, of which 80 DMPs (DMP-CpGis), 112 DMEs (DME-CpGis) and 166 DMIs (DMI-CpGis) were discovered in CpG islands. Importantly, GPX1 contained one each of DMP, DME, DMI, DMP-CpGi, DME-CpGi and DMI-CpGi. Gene-GO term relationships and pathways analysis showed DMP-CpGi-related genes are mainly involved in methylation-related biological functions. In addition, gene-gene interaction networks consisted of nodes that were hypo-methylated GPX1, hypo-methylated APP, hypo-methylated ATOX1, hyper-methylated ADRB2, hyper-methylated RPS6KA1 and hyper-methylated PNMT. They could be used as candidate biomarkers for reducing boar taint in pigs, after further validation in large cohorts.

An ectomycorrhizal fungus alters sensitivity to jasmonate, salicylate, gibberellin, and ethylene in host roots

The phytohormones jasmonate, gibberellin, salicylate, and ethylene regulate an interconnected reprogramming network integrating root development with plant responses against microbes. The establishment of mutualistic ectomycorrhizal symbiosis requires the suppression of plant defense responses against fungi as well as the modification of root architecture and cortical cell wall properties. Here, we investigated the contribution of phytohormones and their crosstalk to the ontogenesis of ectomycorrhizae (ECM) between grey poplar (Populus tremula x alba) roots and the fungus Laccaria bicolor. To obtain the hormonal blueprint of developing ECM, we quantified the concentrations of jasmonates, gibberellins, and salicylate via liquid chromatography-tandem mass spectrometry. Subsequently, we assessed root architecture, mycorrhizal morphology, and gene expression levels (RNA sequencing) in phytohormone-treated poplar lateral roots in the presence or absence of L. bicolor. Salicylic acid accumulated in mid-stage ECM. Exogenous phytohormone treatment affected the fungal colonization rate and/or frequency of Hartig net formation. Colonized lateral roots displayed diminished responsiveness to jasmonate but regulated some genes, implicated in defense and cell wall remodelling, that were specifically differentially expressed after jasmonate treatment. Responses to salicylate, gibberellin, and ethylene were enhanced in ECM. The dynamics of phytohormone accumulation and response suggest that jasmonate, gibberellin, salicylate, and ethylene signalling play multifaceted roles in poplar L. bicolor ectomycorrhizal development.

Functional Annotation of the Transcriptome of the Pig, Sus scrofa, Based Upon Network Analysis of an RNAseq Transcriptional Atlas

The domestic pig (Sus scrofa) is both an economically important livestock species and a model for biomedical research. Two highly contiguous pig reference genomes have recently been released. To support functional annotation of the pig genomes and comparative analysis with large human transcriptomic data sets, we aimed to create a pig gene expression atlas. To achieve this objective, we extended a previous approach developed for the chicken. We downloaded RNAseq data sets from public repositories, down-sampled to a common depth, and quantified expression against a reference transcriptome using the mRNA quantitation tool, Kallisto. We then used the network analysis tool Graphia to identify clusters of transcripts that were coexpressed across the merged data set. Consistent with the principle of guilt-by-association, we identified coexpression clusters that were highly tissue or cell-type restricted and contained transcription factors that have previously been implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes, such as the cell cycle and oxidative phosphorylation. The same approach was used to identify coexpression clusters within RNAseq data from multiple individual liver and brain samples, highlighting cell type, process, and region-specific gene expression. Evidence of conserved expression can add confidence to assignment of orthology between pig and human genes. Many transcripts currently identified as novel genes with ENSSSCG or LOC IDs were found to be coexpressed with annotated neighbouring transcripts in the same orientation, indicating they may be products of the same transcriptional unit. The meta-analytic approach to utilising public RNAseq data is extendable to include new data sets and new species and provides a framework to support the Functional Annotation of Animals Genomes (FAANG) initiative.

Metabolomic networks and pathways associated with feed efficiency and related-traits in Duroc and Landrace pigs

Improving feed efficiency (FE) is a major goal of pig breeding, reducing production costs and providing sustainability to the pig industry. Reliable predictors for FE could assist pig producers. We carried out untargeted blood metabolite profiling in uncastrated males from Danbred Duroc (n = 59) and Danbred Landrace (n = 50) pigs at the beginning and end of a FE testing phase to identify biomarkers and biological processes underlying FE and related traits. By applying linear modeling and clustering analyses coupled with WGCNA framework, we identified 102 and 73 relevant metabolites in Duroc and Landrace based on two sampling time points. Among them, choline and pyridoxamine were hub metabolites in Duroc in early testing phase, while, acetoacetate, cholesterol sulfate, xanthine, and deoxyuridine were identified in the end of testing. In Landrace, cholesterol sulfate, thiamine, L-methionine, chenodeoxycholate were identified at early testing phase, while, D-glutamate, pyridoxamine, deoxycytidine, and L-2-aminoadipate were found at the end of testing. Validation of these results in larger populations could establish FE prediction using metabolomics biomarkers. We conclude that it is possible to identify a link between blood metabolite profiles and FE. These results could lead to improved nutrient utilization, reduced production costs, and increased FE.

Whole blood transcriptome analysis in ewes fed with hemp seed supplemented diet

The hemp plant (Cannabis sativa L.) has a long tradition of being used for many different purposes such as industry, medicine and nutrition. In particular, because hemp seed (HS) is rich in oil protein and considerable amounts of dietary fiber, vitamins and minerals that are particularly suitable also for animal nutrition. Different studies have evaluated HS on qualitative and quantitative properties of livestock products but as of today, nobody has investigated the molecular pathway behind HS supplementation in farm animals. Thus, in this study, we will report the first RNA sequencing of the whole-blood transcriptome of ewes fed either with a controlled diet (CTR, n = 5) or with a diet supplemented with 5% of hemp seed (HSG, n = 5). Applying a false discovery rate (FDR) <0.05 and a log₂FC either higher than 0.5 or lower than −0.5, we identified 314 differentially regulated genes in the HS-supplemented group compared to the CTR group. Several genes encoding for different subunits belonging to the complex I, II, III, IV and ATP-synthase were up-regulated making oxidative phosphorylation (FDR: 3.05e-19) and thermogenesis (FDR: 2.17e-16) the highest up-regulated pathways in our study. Moreover, we found up-regulation in different genes involved in lactose biosyntheses such as GALK1 and PGM1 and, as a result, we observed a statistically higher lactose percentage in the HSG group (p < 0.05). These results indicate that HS supplementation positively affects the energy production pathway in lactating ewes conferring them also more resistance to adverse climatic conditions such as low temperature. Finally, the higher milk lactose content makes the derived dairy products more profitable.

Whole blood transcriptome analysis in ewes fed with hemp seed supplemented diet

The hemp plant (Cannabis sativa L.) has a long tradition of being used for many different purposes such as industry, medicine and nutrition. In particular, because hemp seed (HS) is rich in oil protein and considerable amounts of dietary fiber, vitamins and minerals that are particularly suitable also for animal nutrition. Different studies have evaluated HS on qualitative and quantitative properties of livestock products but as of today, nobody has investigated the molecular pathway behind HS supplementation in farm animals. Thus, in this study, we will report the first RNA sequencing of the whole-blood transcriptome of ewes fed either with a controlled diet (CTR, n = 5) or with a diet supplemented with 5% of hemp seed (HSG, n = 5). Applying a false discovery rate (FDR) <0.05 and a log₂FC either higher than 0.5 or lower than -0.5, we identified 314 differentially regulated genes in the HS-supplemented group compared to the CTR group. Several genes encoding for different subunits belonging to the complex I, II, III, IV and ATP-synthase were up-regulated making oxidative phosphorylation (FDR: 3.05e-19) and thermogenesis (FDR: 2.17e-16) the highest up-regulated pathways in our study. Moreover, we found up-regulation in different genes involved in lactose biosyntheses such as GALK1 and PGM1 and, as a result, we observed a statistically higher lactose percentage in the HSG group (p < 0.05). These results indicate that HS supplementation positively affects the energy production pathway in lactating ewes conferring them also more resistance to adverse climatic conditions such as low temperature. Finally, the higher milk lactose content makes the derived dairy products more profitable.

Characterization of eQTLs associated with androstenone by RNA sequencing in porcine testis

Characterization of genetic variants affecting genome-wide gene expression levels (expression quantitative trait loci or eQTLs) in pig testes may improve our understanding of genetic architecture of boar taint (an animal welfare trait) and helps in genome-assisted or genomic selection programs. The aims of this study were to identify eQTLs associated with androstenone, to find candidate eQTLs for low androstenone, and to validate the top eQTL by reverse transcriptase quantitative PCR (RT-qPCR). Gene expression profiles were obtained by RNA sequencing in testis from Danish cross-bred pigs and genotype data by 80K single nucleotide polymorphism panel. A total of 262 eQTLs [false discovery rate (FDR) < 0.05] were identified by using two software packages: Matrix eQTL and Krux eQTL. Of these, 149 cis-acting eQTLs were significantly associated with androstenone concentrations and gene expression (FDR < 0.05). The eQTLs were associated with several genes of boar taint relevance including CYP1A2, CYB5D1, and SPHK2. One eQTL gene, AMPH, was differentially expressed (FDR < 0.05) and affected by chicory. Five candidate eQTLs associated with low androstenone concentrations were discovered, including the top eQTL associated with CYP1A2. RT-qPCR confirmed target gene expression to be significantly (P < 0.05) different based on eQTL genotypes. Furthermore, eQTLs were enriched as QTLs for 15 boar taint related traits from the PigQTLdb. This is the first study to report eQTLs in testes of commercial crossbred pigs used in pork production and to reveal genetic architecture of boar taint. Potential applications include development of a DNA test and in advanced genomic selection models for boar taint.

Genome-wide DNA methylation analysis using next-generation sequencing to reveal candidate genes responsible for boar taint in pigs

Boar taint (BT) is an offensive flavor observed in non-castrated male pigs that reduces the carcass price. Surgical castration effectively avoids the taint but is associated with animal welfare concerns. The functional annotation of farm animal genomes for understanding the biology of complex traits can be used in the selection of breeding animals to achieve favorable phenotypic outcomes. The characterization of pig epigenomes/methylation changes between animals with high and low BT and genome-wide epigenetic markers that can predict BT are lacking. Reduced representation bisulfite sequencing of DNA methylation patterns based on next-generation sequencing is an efficient technology to identify candidate epigenetic biomarkers associated with BT. Three different BT levels were analyzed using reduced representation bisulfite sequencing data to calculate the methylation levels of cytosine and guanine dinucleotide (CpG) sites. The co-analysis of differentially methylated CpG sites identified by this study and differentially expressed genes identified by a previous study found 32 significant co-located genes. The joint analysis of GO terms and pathways revealed that methylation and gene expression of seven candidate genes were associated with BT; in particular, FASN plays a key role in fatty acid biosynthesis, and PEMT might be involved in estrogen regulation and the development of BT. This study is the first to report the genome-wide DNA methylation profiles of BT in pigs using next-generation sequencing and summarize candidate genes associated with epigenetic markers of BT, which could contribute to the understanding of the functional biology of BT traits and selective breeding of pigs against BT based on epigenetic biomarkers.

Transcription Profiling of Monocyte-Derived Macrophages Infected In Vitro With Two Strains of Streptococcus agalactiae Reveals Candidate Pathways Affecting Subclinical Mastitis in Cattle

Macrophages are key cells of innate immune response and serve as the first line of defense against bacteria. Transcription profiling of bacteria-infected macrophages could provide important insights on the pathogenicity and host defense mechanisms during infection. We have examined transcription profiles of bovine monocyte-derived macrophages (bMDMs) isolated from the blood of 12 animals and infected in vitro with two strains of Streptococcus agalactiae. Illumina sequencing of RNA from 36 bMDMs cultures exposed in vitro to either one of two sequence types of S. agalactiae (ST103 or ST12) for 6 h and unchallenged controls was performed. Analyses of over 1,656 million high-quality paired-end sequence reads revealed 5,936 and 6,443 differentially expressed genes (p < 0.05) in bMDMs infected with ST103 and ST12, respectively, versus unchallenged controls. Moreover, 588 genes differentially expressed between bMDMs infected with ST103 versus ST12 were identified. Ingenuity pathway analysis of the differentially up-regulated genes in the bMDMs infected with ST103 revealed significant enrichment for granulocyte adhesion and diapedesis, while significant enrichment for the phagosome formation pathway was found among down-regulated genes. Moreover, Ingenuity pathway analysis of the differentially up-regulated genes in the bMDMs infected with ST12 showed significant enrichment for type 1/type 2 T helper cell activation, while the complement activation pathway was overrepresented in the down-regulated genes. Our study identified pathogen-induced regulation of key genes and pathways involved in the immune response of macrophages against infection but also likely involved in bacterial evasion of the host immune system. These results may contribute to better understanding of the mechanisms underlying subclinical infection such as bovine streptococcal mastitis.

An Epigenome-Wide DNA Methylation Map of Testis in Pigs for Study of Complex Traits

Epigenetic changes are important for understanding complex trait variation and inheritance in pigs that are also a valuable biomedical model for human health research. Testis is the main organ for reproduction and boar taint in pigs; however, there have been no studies to-date on adult pig testis epigenome. The main objective of this study was to establish a genome-wide DNA methylation map of pig testis that would help identify candidate epigenetic biomarkers and methylated genes for complex traits such as male reproduction, fertility or boar taint. Reduced Representation Bisulfite Sequencing (RRBS) was used to study methylation levels of cytosine in nine pig testis samples. The results showed that genome-wide methylation status of nine samples overlapped greatly and their variation among pigs were low. The methylation levels of promoter, exon, intron, cytosine and guanine dinucleotide (CpG) islands and CpG island shores regions were 0.15, 0.47, 0.55, 0.39, and 0.53, respectively. Cytosines binding to CpG islands showed different methylation levels between exon and intron regions. All methylation levels of CpG islands were lower than CpG island shores in different genic features. The distribution of 12,738 differentially methylated cytosines (DMCs) within CpG islands, CpG island shores and other regions was 36.86, 21.65, and 41.49%, respectively, and was 0.33, 1.71, 5.95, and 92.01% in promoter, exon, intron and intergenic regions, respectively. Methylation levels of DMCs in promoter, exon and intron regions were significantly different between CpG islands and CpG island shores (P < 0.05). A total of 898 genes with 2089 DMCs were enriched in 112 Gene Ontology (GO) terms. Fifteen methylated genes from our study were associated with fertility or boar taint traits. Our analysis revealed the methylation patterns in different genic features and CpG island regions of testis in pigs, and summarized several candidate genes associated with DMCs and the involved GO terms. These findings are helpful to understand the relationship between DNA methylation and genic CpG islands, to provide candidate epigenetic regions or biomarkers for pig production and welfare and for translational epigenomic studies that use pigs as an animal model for human research.

Gene Co-expression Analysis Indicates Potential Pathways and Regulators of Beef Tenderness in Nellore Cattle

Beef tenderness, a complex trait affected by many factors, is economically important to beef quality, industry, and consumer’s palatability. In this study, RNA-Seq was used in network analysis to better understand the biological processes that lead to differences in beef tenderness. Skeletal muscle transcriptional profiles from 24 Nellore steers, selected by extreme estimated breeding values (EBVs) for shear force after 14 days of aging, were analyzed and 22 differentially expressed transcripts were identified. Among these were genes encoding ribosomal proteins, glutathione transporter ATP-binding cassette, sub-family C (CFTR/MRP), member 4 (ABCC4), and synaptotagmin IV (SYT4). Complementary co-expression analyses using Partial Correlation with Information Theory (PCIT), Phenotypic Impact Factor (PIF) and the Regulatory Impact Factor (RIF) methods identified candidate regulators and related pathways. The PCIT analysis identified ubiquitin specific peptidase 2 (USP2), growth factor receptor-bound protein 10 (GBR10), anoctamin 1 (ANO1), and transmembrane BAX inhibitor motif containing 4 (TMBIM4) as the most differentially hubbed (DH) transcripts. The transcripts that had a significant correlation with USP2, GBR10, ANO1, and TMBIM4 enriched for proteasome KEGG pathway. RIF analysis identified microRNAs as candidate regulators of variation in tenderness, including bta-mir-133a-2 and bta-mir-22. Both microRNAs have target genes present in the calcium signaling pathway and apoptosis. PIF analysis identified myoglobin (MB), enolase 3 (ENO3), and carbonic anhydrase 3 (CA3) as potentially having fundamental roles in tenderness. Pathways identified in our study impacted in beef tenderness included: calcium signaling, apoptosis, and proteolysis. These findings underscore some of the complex molecular mechanisms that control beef tenderness in Nellore cattle.

Systems genomics study reveals expression quantitative trait loci, regulator genes and pathways associated with boar taint in pigs

Boar taint is an offensive odour and/or taste from a proportion of non-castrated male pigs caused by skatole and androstenone accumulation during sexual maturity. Castration is widely used to avoid boar taint but is currently under debate because of animal welfare concerns. This study aimed to identify expression quantitative trait loci (eQTLs) with potential effects on boar taint compounds to improve breeding possibilities for reduced boar taint. Danish Landrace male boars with low, medium and high genetic merit for skatole and human nose score (HNS) were slaughtered at ~100 kg. Gene expression profiles were obtained by RNA-Seq, and genotype data were obtained by an Illumina 60K Porcine SNP chip. Following quality control and filtering, 10,545 and 12,731 genes from liver and testis were included in the eQTL analysis, together with 20,827 SNP variants. A total of 205 and 109 single-tissue eQTLs associated with 102 and 58 unique genes were identified in liver and testis, respectively. By employing a multivariate Bayesian hierarchical model, 26 eQTLs were identified as significant multi-tissue eQTLs. The highest densities of eQTLs were found on pig chromosomes SSC12, SSC1, SSC13, SSC9 and SSC14. Functional characterisation of eQTLs revealed functions within regulation of androgen and the intracellular steroid hormone receptor signalling pathway and of xenobiotic metabolism by cytochrome P450 system and cellular response to oestradiol. A QTL enrichment test revealed 89 QTL traits curated by the Animal Genome PigQTL database to be significantly overlapped by the genomic coordinates of cis-acting eQTLs. Finally, a subset of 35 cis-acting eQTLs overlapped with known boar taint QTL traits. These eQTLs could be useful in the development of a DNA test for boar taint but careful monitoring of other overlapping QTL traits should be performed to avoid any negative consequences of selection.