In vitro- and in vivo-produced male dairy calves show molecular differences in the hepatic and muscular energy regulation†.

In cattle, the in vitro production (IVP) of embryos is becoming more relevant than embryos produced in vivo, i.e. after multiple ovulation and embryo transfer (MOET). However, the effects of IVP on the developmental programming of specific organs in the postnatal calves are yet unknown. Previously, we reported an epigenomic and transcriptomic profile of the hypothalamus-pituitary-testicular axis compatible with its earlier activation in IVP calves compared to MOET animals. Here, we studied the hepatic and muscular epigenome and transcriptome of those same male dairy calves (n = 4 per group). Tissue samples from liver and semitendinosus muscle were obtained at 3 months of age, and the extracted gDNA and RNA were sequenced through whole-genome bisulfite sequencing and RNA-sequencing, respectively. Next, bioinformatic analyses determined differentially methylated cytosines or differentially expressed genes [false discovery rate (FDR) < 0.05] for each Omic dataset; and nonparametrically combined genes (NPCG) for both integrated omics (P < 0.05). KEGG pathways enrichment analysis showed that NPCG upregulated in the liver and the muscle of the IVP calves were involved in oxidative phosphorylation and the tricarboxylic acid cycle. In contrast, ribosome and translation were upregulated in the liver but downregulated in the muscle of the IVP calves compared to the MOET calves (FDR < 0.05). A model considering the effect of the methylation levels and the group on the expression of all the genes involved in these pathways confirmed these findings. In conclusion, the multiomics data integration approach indicated an altered hepatic and muscular energy regulation in phenotypically normal IVP calves compared to MOET calves.

Epigenomic and transcriptomic analyses reveal early activation of the HPG axis in in vitro-produced male dairy calves.

In cattle, several calves born after IVP ("in vitro" embryo production) present similar birthweight to those generated after MOET (multiple ovulation and embryo transfer). However, the underlying molecular patterns in organs involved in the developmental process are unknown and could indicate physiological programming. The objectives of this study were: (1) to compare epigenomic and transcriptomic modifications in the hypothalamus, pituitary, gonadal and adrenal organs between 3 months old ovum pick-up-IVP and MOET male calves (n = 4 per group) and (2) to use blood epigenomic data to proxy methylation of the inner organs. Extracted gDNA and RNA were sequenced through whole-genome bisulfite sequencing and RNA sequencing, respectively. Next, bioinformatic analyses determined differentially methylated cytosines (DMC) and differentially expressed genes (DEG) (FDR < 0.05) in IVP versus MOET samples and the KEGG pathways that were overrepresented by genes associated with DMC or DEG (FDR < 0.1). Pathways related to hypothalamus, pituitary, gonadal (HPG) axis activation (GnRH secretion in the hypothalamus, GnRH signaling in the pituitary, and steroidogenesis in the testicle) were enriched in IVP calves. Modeling the effect of the methylation levels and the group on the expression of all the genes involved in these pathways confirmed their upregulation in HPG organs in IVP calves. The application of the DIABLO method allowed the identification of 15 epigenetic and five transcriptomic biomarkers, which were able to predict the embryo origin using the epigenomic data from the blood. In conclusion, the use of an integrated epigenomic-transcriptomic approach suggested an early activation of the HPG axis in male IVP calves compared to MOET counterparts, and the identification of potential biomarkers allowed the use of blood samples to proxy methylation levels of the relevant internal organs.

Impact of IL-17F 7488T/C Functional Polymorphism on Progressive Rheumatoid Arthritis: Novel Insight from the Molecular Dynamic Simulations.

Resorption of bones and cartilage coupled with structural changes in the inflamed joints are the major hallmark of rheumatoid arthritis (RA). Genetic polymorphisms in pro-inflammatory interleukins (ILs) appear to play an important role in the susceptibility towards progressive RA. We therefore aimed to investigate the association of single nucleotide polymorphisms (SNP), present in the hotspot coding/promoter regions of IL-6, -17 and -18, with RA susceptibility or severity in a larger study cohort from Pakistan together with finding clues as to how a functional SNP impacts the predisposition towards RA. TaqMan SNP genotyping approach was first used to assess IL-6 (rs1800795), IL-17 F (rs763780), IL-17A (rs2275913), and IL-18 (rs1946518) polymorphisms in 310 subjects (150 RA and 160 control). Molecular dynamic simulations (MDS) of wild- and mutant-type IL-17A with corresponding receptor were thereafter performed using AMBER-16; Chimera 1.13 was used for analyses. Our results showed the association of two SNPs, namely IL-6 - 174 G/C [allelic (OR = 0.960, 95% CI = 0.929-0.992, p = .009)] and IL-17 F 7488 T/C [allelic (OR = 0.907, 95%CI = 0.861-0.954, p = .000)] with increased RA risk in Pakistani subjects. When mapped, IL-17 F 7488 T/C was found involved in His161→Arg161 change near the C-terminus of IL-17 F. Comparative MDS revealed enhanced stability of the mutant hence advocating a potential role of IL-17F functional SNP in RA susceptibility and/or severity. This study provides a novel structural insight for SNP-derived functional mutation and its overall impact on binding with heterotrimeric receptor complex of IL-17 receptor thereby opening new avenues for understanding the biochemical basis of the disease.

Characterization of the endometrial transcriptome in early diestrus influencing pregnancy status in dairy cattle after transfer of in vitro-produced embryos.

Modifications of the endometrial transcriptome at day 7 of the estrus cycle are crucial to maintain gestation after transfer of in vitro-produced (IVP) embryos, although these changes are still largely unknown. The aim of this study was to identify genes, and their related biological mechanisms, important for pregnancy establishment based on the endometrial transcriptome of recipient lactating dairy cows that become pregnant in the subsequent estrus cycle, upon transfer of IVP embryos. Endometrial biopsies were taken from Holstein Friesian cows on day 6-8 of the estrus cycle followed by embryo transfer in the following cycle. Animals were classified retrospectively as pregnant (PR, n = 8) or nonpregnant (non-PR, n = 11) cows, according to pregnancy status at 26-47 days. Extracted mRNAs from endometrial samples were sequenced with an Illumina platform to determine differentially expressed genes (DEG) between the endometrial transcriptome from PR and non-PR cows. There were 111 DEG (false discovery rate < 0.05), which were mainly related to extracellular matrix interaction, histotroph metabolic composition, prostaglandin synthesis, transforming growth factor-ß signaling as well as inflammation and leukocyte activation. Comparison of these DEG with DEG identified in two public external data sets confirmed the more fertile endometrial molecular profile of PR cows. In conclusion, this study provides insights into the key early endometrial mechanisms for pregnancy establishment, after IVP embryo transfer in dairy cows.

Interplay among miR-29 family, mineral metabolism, and gene regulation in Bos indicus muscle.

An interplay between gene expression, mineral concentration, and beef quality traits in Bos indicus muscle has been reported previously under a network approach. However, growing evidence suggested that miRNAs not only modulate gene expression but are also involved with mineral homeostasis. To our knowledge, understanding of the miRNA-gene expression-mineral concentration relationship in mammals is still minimal. Therefore, we carried out a miRNA co-expression and multi-level miRNA-mRNA integration analyses to predict the putative drivers (miRNAs and genes) associated with muscle mineral concentration in Nelore steers. In this study, we identified calcium and iron to be the pivotal minerals associated with miRNAs and gene targets. Furthermore, we identified the miR-29 family (miR-29a, -29b, -29c, -29d-3p, and -29e) as the putative key regulators modulating mineral homeostasis. The miR-29 family targets genes involved with AMPK, insulin, mTOR, and thyroid hormone signaling pathways. Finally, we reported an interplay between miRNAs and minerals acting cooperatively to modulate co-expressed genes and signaling pathways both involved with mineral and energy homeostasis in Nelore muscle. Although we provided some evidence to understand this complex relationship, future work should determine the functional implications of minerals for miRNA levels and their feedback regulation system.\\An interplay between gene expression, mineral concentration, and beef quality traits in Bos indicus muscle has been reported previously under a network approach. However, growing evidence suggested that miRNAs not only modulate gene expression but are also involved with mineral homeostasis. To our knowledge, understanding of the miRNA-gene expression-mineral concentration relationship in mammals is still minimal. Therefore, we carried out a miRNA co-expression and multi-level miRNA-mRNA integration analyses to predict the putative drivers (miRNAs and genes) associated with muscle mineral concentration in Nelore steers. In this study, we identified calcium and iron to be the pivotal minerals associated with miRNAs and gene targets. Furthermore, we identified the miR-29 family (miR-29a, -29b, -29c, -29d-3p, and -29e) as the putative key regulators modulating mineral homeostasis. The miR-29 family targets genes involved with AMPK, insulin, mTOR, and thyroid hormone signaling pathways. Finally, we reported an interplay between miRNAs and minerals acting cooperatively to modulate co-expressed genes and signaling pathways both involved with mineral and energy homeostasis in Nelore muscle. Although we provided some evidence to understand this complex relationship, future work should determine the functional implications of minerals for miRNA levels and their feedback regulation system.

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.

Simulated physiological oocyte maturation has side effects on bovine oocytes and embryos.

PURPOSE: Oocyte maturation is a complex process involving nuclear and cytoplasmic modulations, during which oocytes acquire their ability to become fertilized and support embryonic development. The oocyte is apparently "primed" for maturation during its development in the dominant follicle. As bovine oocytes immediately resume meiosis when cultured, it was hypothesized that delaying resumption of meiosis with cyclic nucleotide modulators before in vitro maturation (IVM) would allow the oocytes to acquire improved developmental competence. METHODS: We tested the Simulated Physiological Oocyte Maturation (SPOM) system that uses forskolin and 3-isobutyl-1-methylxanthine for 2 h prior to IVM against two different systems of conventional IVM (Con-IVM). We evaluated the ultrastructure of matured oocytes and blastocysts and also assessed the expression of 96 genes related to embryo quality in the blastocysts. RESULTS: In summary, the SPOM system resulted in lower blastocyst rates than both Con-IVM systems (30 ± 9.1 vs. 35 ± 8.7; 29 ± 2.6 vs. 38 ± 2.8). Mature SPOM oocytes had significantly increased volume and number of vesicles, reduced volume and surface density of large smooth endoplasmic reticulum clusters, and lower number of mitochondria than Con-IVM oocytes. SPOM blastocysts showed only subtle differences with parallel undulations of adjacent trophectoderm plasma membranes and peripherally localized ribosomes in cells of the inner cell mass compared with Con-IVM blastocysts. SPOM blastocysts, however, displayed significant downregulation of genes related to embryonic developmental potential when compared to Con-IVM blastocysts. CONCLUSIONS: Our results show that the use of the current version of the SPOM system may have adverse effects on oocytes and blastocysts calling for optimized protocols for improving oocyte competence.

WISH-R- a fast and efficient tool for construction of epistatic networks for complex traits and diseases.

BACKGROUND: Genetic epistasis is an often-overlooked area in the study of the genomics of complex traits. Genome-wide association studies are a useful tool for revealing potential causal genetic variants, but in this context, epistasis is generally ignored. Data complexity and interpretation issues make it difficult to process and interpret epistasis. As the number of interaction grows exponentially with the number of variants, computational limitation is a bottleneck. Gene Network based strategies have been successful in integrating biological data and identifying relevant hub genes and pathways related to complex traits. In this study, epistatic interactions and network-based analysis are combined in the Weighted Interaction SNP hub (WISH) method and implemented in an efficient and easy to use R package. RESULTS: The WISH R package (WISH-R) was developed to calculate epistatic interactions on a genome-wide level based on genomic data. It is easy to use and install, and works on regular genomic data. The package filters data based on linkage disequilibrium and calculates epistatic interaction coefficients between SNP pairs based on a parallelized efficient linear model and generalized linear model implementations. Normalized epistatic coefficients are analyzed in a network framework, alleviating multiple testing issues and integrating biological signal to identify modules and pathways related to complex traits. Functions for visualizing results and testing runtimes are also provided. CONCLUSION: The WISH-R package is an efficient implementation for analyzing genome-wide epistasis for complex diseases and traits. It includes methods and strategies for analyzing epistasis from initial data filtering until final data interpretation. WISH offers a new way to analyze genomic data by combining epistasis and network based analysis in one method and provides options for visualizations. This alleviates many of the existing hurdles in the analysis of genomic interactions.

An adult-based insulin resistance genetic risk score associates with insulin resistance, metabolic traits and altered fat distribution in Danish children and adolescents who are overweight or obese.

AIMS/HYPOTHESIS: A genetic risk score (GRS) consisting of 53 insulin resistance variants (GRS53) was recently demonstrated to associate with insulin resistance in adults. We speculated that the GRS53 might already associate with insulin resistance during childhood, and we therefore aimed to investigate this in populations of Danish children and adolescents. Furthermore, we aimed to address whether the GRS associates with components of the metabolic syndrome and altered body composition in children and adolescents. METHODS: We examined a total of 689 children and adolescents who were overweight or obese and 675 children and adolescents from a population-based study. Anthropometric data, dual-energy x-ray absorptiometry scans, BP, fasting plasma glucose, fasting serum insulin and fasting plasma lipid measurements were obtained, and HOMA-IR was calculated. The GRS53 was examined for association with metabolic traits in children by linear regressions using an additive genetic model. RESULTS: In overweight/obese children and adolescents, the GRS53 associated with higher HOMA-IR (ß = 0.109 ± 0.050 (SE); p = 2.73 × 10-2), fasting plasma glucose (ß = 0.010 ± 0.005 mmol/l; p = 2.51 × 10-2) and systolic BP SD score (ß = 0.026 ± 0.012; p = 3.32 × 10-2) as well as lower HDL-cholesterol (ß = -0.008 ± 0.003 mmol/l; p = 1.23 × 10-3), total fat-mass percentage (ß = -0.143 ± 0.054%; p = 9.15 × 10-3) and fat-mass percentage in the legs (ß = -0.197 ± 0.055%; p = 4.09 × 10-4). In the population-based sample of children, the GRS53 only associated with lower HDL-cholesterol concentrations (ß = -0.007 ± 0.003 mmol/l; p = 1.79 × 10-2). CONCLUSIONS/INTERPRETATION: An adult-based GRS comprising 53 insulin resistance susceptibility SNPs associates with insulin resistance, markers of the metabolic syndrome and altered fat distribution in a sample of Danish children and adolescents who were overweight or obese.

Gene co-expression networks in liver and muscle transcriptome reveal sex-specific gene expression in lambs fed with a mix of essential oils.

BACKGROUND: Essential oil (EO) dietary supplementation is a new strategy to improve animal health. EO compounds have antiparasitic, antimicrobial, antiviral, antimycotic, antioxidant and anti-inflammatory proprieties. Nutrigenomics investigations represent innovative approaches in understanding the relation between diet effect and gene expression related to the animal performance. Few nutrigenomics studies have used a high-throughput RNA-Sequencing (RNA-Seq) approach, despite great potential of RNA-Seq data in gene expression quantification and in co-expression network analyses. Our aim is to use the potential of RNA-Sequencing data in order to evaluate the effect of an EO supplementary diet on gene expression in both lamb liver and muscle. RESULTS: Using a treatment and sex interaction model, 13 and 4 differentially expressed genes were identified in liver and muscle respectively. Sex-specific differentially expressed (DE) genes were identified in both sexes. Using network based analysis, different clusters of co-expressed genes that were highly correlated to the diet were detected in males vs. females, in agreement with DE analysis. A total of five regulatory genes in liver tissue associated to EO diet were identified: DNAJB9, MANF, UFM1, CTNNLA1 and NFX1. Our study reveals a sex-dependent effect of EO diet in both tissues, and an influence on the expression of genes mainly involved in immune, inflammatory and stress pathway. CONCLUSION: Our analysis suggests a sex-dependent effect of the EO dietary supplementation on the expression profile of both liver and muscle tissues. We hypothesize that the presence of EOs could have beneficial effects on wellness of male lamb and further analyses are needed to understand the biological mechanisms behind the different effect of EO metabolites based on sex. Using lamb as a model for nutrigenomics studies, it could be interesting to investigate the effects of EO diets in other species and in humans.

Genome-wide association analysis identifies three new susceptibility loci for childhood body mass index.

A large number of genetic loci are associated with adult body mass index. However, the genetics of childhood body mass index are largely unknown. We performed a meta-analysis of genome-wide association studies of childhood body mass index, using sex- and age-adjusted standard deviation scores. We included 35 668 children from 20 studies in the discovery phase and 11 873 children from 13 studies in the replication phase. In total, 15 loci reached genome-wide significance (P-value < 5 × 10(-8)) in the joint discovery and replication analysis, of which 12 are previously identified loci in or close to ADCY3, GNPDA2, TMEM18, SEC16B, FAIM2, FTO, TFAP2B, TNNI3K, MC4R, GPR61, LMX1B and OLFM4 associated with adult body mass index or childhood obesity. We identified three novel loci: rs13253111 near ELP3, rs8092503 near RAB27B and rs13387838 near ADAM23. Per additional risk allele, body mass index increased 0.04 Standard Deviation Score (SDS) [Standard Error (SE) 0.007], 0.05 SDS (SE 0.008) and 0.14 SDS (SE 0.025), for rs13253111, rs8092503 and rs13387838, respectively. A genetic risk score combining all 15 SNPs showed that each additional average risk allele was associated with a 0.073 SDS (SE 0.011, P-value = 3.12 × 10(-10)) increase in childhood body mass index in a population of 1955 children. This risk score explained 2% of the variance in childhood body mass index. This study highlights the shared genetic background between childhood and adult body mass index and adds three novel loci. These loci likely represent age-related differences in strength of the associations with body mass index.

Transcription Factor Co-expression Networks of Adipose RNA-Seq Data Reveal Regulatory Mechanisms of Obesity.

BACKGROUND: Transcription Factors (TFs) control actuation of genes in the genome and are key mediators of complex processes such as obesity. Master Regulators (MRs) are the genes at the top of a regulation hierarchy which regulate other genes. OBJECTIVE: To elucidate clusters of highly co-expressed TFs (modules), involved pathways, highly inter-connected TFs (hub-TFs) and MRs leading to obesity and leanness, using porcine model for human obesity. METHODS: We identified 817 expressed TFs in RNA-Sequencing dataset representing extreme degrees of obesity (DO; lean, obese). We built a single Weighted Transcription Factor Co-expression Network (WTFCN) and TF sub-networks (based on the DO). Hub-TFs and MRs (using iRegulon) were identi-fied in biologically relevant WTFCNs modules. RESULTS: Single WTFCN detected the Red module significantly associated with DO (P < 0.03). This module was enriched for regulation processes in the immune system, e.g.: Immune system process (Padj = 2.50E-06) and metabolic lifestyle disorders, e.g. Circadian rhythm - mammal pathway (Padj = 2.33E-11). Detected MR, hub-TF SPI1 was involved in obesity, immunity and osteoporosis. Within the obese sub-network, the Red module suggested possible associations with immunity, e.g. TGF-beta signaling pathway (Padj = 1.73E-02) and osteoporosis, e.g. Osteoclast differentiation (Padj = 1.94E-02). Within the lean sub-network, the Magenta module displayed associations with type 2 diabetes, obesity and os-teoporosis e.g. Notch signaling pathway (Padj = 2.40E-03), osteoporosis e.g. hub-TF VDR (a prime candidate gene for osteoporosis). CONCLUSION: Our results provide insights into the regulatory network of TFs and biologically relevant hub TFs in obesity.

Systematic in vitro and in vivo characterization of Leukemia-inhibiting factor- and Fibroblast growth factor-derived porcine induced pluripotent stem cells.

Derivation and stable maintenance of porcine induced pluripotent stem cells (piPSCs) is challenging. We herein systematically analyzed two piPSC lines, derived by lentiviral transduction and cultured under either leukemia inhibitory factor (LIF) or fibroblast growth factor (FGF) conditions, to shed more light on the underlying biological mechanisms of porcine pluripotency. LIF-derived piPSCs were more successful than their FGF-derived counterparts in the generation of in vitro chimeras and in teratoma formation. When LIF piPSCs chimeras were transferred into surrogate sows and allowed to develop, only their prescence within the embryonic membranes could be detected. Whole-transcriptome analysis of the piPSCs and porcine neonatal fibroblasts showed that they clustered together, but apart from the two pluripotent cell populations of early porcine embryos, indicating incomplete reprogramming. Indeed, bioinformatic analysis of the pluripotency-related gene network of the LIF- versus FGF-derived piPSCs revealed that ZFP42 (REX1) expression was absent in both piPSC-like cells, whereas it was expressed in the porcine inner cell mass at Day 7/8. A second striking difference was the expression of ATOH1 in piPSC-like cells, which was absent in the inner cell mass. Moreover, our gene expression analyses plus correlation analyses of known pluripotency genes identified unique relationships between pluripotency genes in the inner cell mass, which are to some extent, in the piPSC-like cells. This deficiency in downstream gene activation and divergent gene expression may be underlie the inability to derive germ line-transmitting piPSCs, and provides unique insight into which genes are necessary to achieve fully reprogrammed piPSCs. 84: 229-245, 2017. © 2016 Wiley Periodicals, Inc.

Multi-omic data integration and analysis using systems genomics approaches: methods and applications in animal production, health and welfare.

In the past years, there has been a remarkable development of high-throughput omics (HTO) technologies such as genomics, epigenomics, transcriptomics, proteomics and metabolomics across all facets of biology. This has spearheaded the progress of the systems biology era, including applications on animal production and health traits. However, notwithstanding these new HTO technologies, there remains an emerging challenge in data analysis. On the one hand, different HTO technologies judged on their own merit are appropriate for the identification of disease-causing genes, biomarkers for prevention and drug targets for the treatment of diseases and for individualized genomic predictions of performance or disease risks. On the other hand, integration of multi-omic data and joint modelling and analyses are very powerful and accurate to understand the systems biology of healthy and sustainable production of animals. We present an overview of current and emerging HTO technologies each with a focus on their applications in animal and veterinary sciences before introducing an integrative systems genomics framework for analysing and integrating multi-omic data towards improved animal production, health and welfare. We conclude that there are big challenges in multi-omic data integration, modelling and systems-level analyses, particularly with the fast emerging HTO technologies. We highlight existing and emerging systems genomics approaches and discuss how they contribute to our understanding of the biology of complex traits or diseases and holistic improvement of production performance, disease resistance and welfare.

Liver transcriptomic networks reveal main biological processes associated with feed efficiency in beef cattle.

BACKGROUND: The selection of beef cattle for feed efficiency (FE) traits is very important not only for productive and economic efficiency but also for reduced environmental impact of livestock. Considering that FE is multifactorial and expensive to measure, the aim of this study was to identify biological functions and regulatory genes associated with this phenotype. RESULTS: Eight genes were differentially expressed between high and low feed efficient animals (HFE and LFE, respectively). Co-expression analyses identified 34 gene modules of which 4 were strongly associated with FE traits. They were mainly enriched for inflammatory response or inflammation-related terms. We also identified 463 differentially co-expressed genes which were functionally enriched for immune response and lipid metabolism. A total of 8 key regulators of gene expression profiles affecting FE were found. The LFE animals had higher feed intake and increased subcutaneous and visceral fat deposition. In addition, LFE animals showed higher levels of serum cholesterol and liver injury biomarker GGT. Histopathology of the liver showed higher percentage of periportal inflammation with mononuclear infiltrate. CONCLUSION: Liver transcriptomic network analysis coupled with other results demonstrated that LFE animals present altered lipid metabolism and increased hepatic periportal lesions associated with an inflammatory response composed mainly by mononuclear cells. We are now focusing to identify the causes of increased liver lesions in LFE animals.

Genome-wide association and systems genetic analyses of residual feed intake, daily feed consumption, backfat and weight gain in pigs.

BACKGROUND: Feed efficiency is one of the major components determining costs of animal production. Residual feed intake (RFI) is defined as the difference between the observed and the expected feed intake given a certain production. Residual feed intake 1 (RFI1) was calculated based on regression of individual daily feed intake (DFI) on initial test weight and average daily gain. Residual feed intake 2 (RFI2) was as RFI1 except it was also regressed with respect to backfat (BF). It has been shown to be a sensitive and accurate measure for feed efficiency in livestock but knowledge of the genomic regions and mechanisms affecting RFI in pigs is lacking. The study aimed to identify genetic markers and candidate genes for RFI and its component traits as well as pathways associated with RFI in Danish Duroc boars by genome-wide associations and systems genetic analyses. RESULTS: Phenotypic and genotypic records (using the Illumina Porcine SNP60 BeadChip) were available on 1,272 boars. Fifteen and 12 loci were significantly associated (p < 1.52 × 10-6) with RFI1 and RFI2, respectively. Among them, 10 SNPs were significantly associated with both RFI1 and RFI2 implying the existence of common mechanisms controlling the two RFI measures. Significant QTL regions for component traits of RFI (DFI and BF) were detected on pig chromosome (SSC) 1 (for DFI) and 2 for (BF). The SNPs within MAP3K5 and PEX7 on SSC 1, ENSSSCG00000022338 on SSC 9, and DSCAM on SSC 13 might be interesting markers for both RFI measures. Functional annotation of genes in 0.5 Mb size flanking significant SNPs indicated regulation of protein and lipid metabolic process, gap junction, inositol phosphate metabolism and insulin signaling pathway are significant biological processes and pathways for RFI, respectively. CONCLUSIONS: The study detected novel genetic variants and QTLs on SSC 1, 8, 9, 13 and 18 for RFI and indicated significant biological processes and metabolic pathways involved in RFI. The study also detected novel QTLs for component traits of RFI. These results improve our knowledge of the genetic architecture and potential biological pathways underlying RFI; which would be useful for further investigations of key candidate genes for RFI and for development of biomarkers.

GEDI: An R Package for Integration of Transcriptomic Data from Multiple Platforms for Bioinformatics Applications.

Transcriptomic data is often expensive and difficult to generate in large cohorts relative to genomic data; therefore, it is often important to integrate multiple transcriptomic datasets from both microarray- and next generation sequencing (NGS)-based transcriptomic data across similar experiments or clinical trials to improve analytical power and discovery of novel transcripts and genes. However, transcriptomic data integration presents a few challenges including reannotation and batch effect removal. We developed the Gene Expression Data Integration (GEDI) R package to enable transcriptomic data integration by combining existing R packages. With just four functions, the GEDI R package makes constructing a transcriptomic data integration pipeline straightforward. Together, the functions overcome the complications in transcriptomic data integration by automatically reannotating the data and removing the batch effect. The removal of the batch effect is verified with principal component analysis and the data integration is verified using a logistic regression model with forward stepwise feature selection. To demonstrate the functionalities of the GEDI package, we integrated five bovine endometrial transcriptomic datasets from the NCBI Gene Expression Omnibus. These transcriptomic datasets were from multiple high-throughput platforms, namely, array-based Affymetrix and Agilent platforms, and NGS-based Illumina paired-end RNA-seq platform. Furthermore, we compared the GEDI package to existing tools and found that GEDI is the only tool that provides a full transcriptomic data integration pipeline including verification of both batch effect removal and data integration for downstream genomic and bioinformatics applications. © 2024 The Author(s). Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: ReadGE, a function to import gene expression datasets Basic Protocol 2: GEDI, a function to reannotate and merge gene expression datasets Basic Protocol 3: BatchCorrection, a function to remove batch effects from gene expression data Basic Protocol 4: VerifyGEDI, a function to confirm successful integration of gene expression data.

Genetic parameters of infectious bovine keratoconjunctivitis and its relationship with weight and parasite infestations in Australian tropical Bos taurus cattle.

BACKGROUND: Infectious bovine keratoconjunctivitis (IBK) or 'pinkeye' is an economically important ocular disease that significantly impacts animal performance. Genetic parameters for IBK infection and its genetic and phenotypic correlations with cattle tick counts, number of helminth (unspecified species) eggs per gram of faeces and growth traits in Australian tropically adapted Bos taurus cattle were estimated. METHODS: Animals were clinically examined for the presence of IBK infection before and after weaning when the calves were 3 to 6 months and 15 to 18 months old, respectively and were also recorded for tick counts, helminth eggs counts as an indicator of intestinal parasites and live weights at several ages including 18 months. RESULTS: Negative genetic correlations were estimated between IBK incidence and weight traits for animals in pre-weaning and post-weaning datasets. Genetic correlations among weight measurements were positive, with moderate to high values. Genetic correlations of IBK incidence with tick counts were positive for the pre-weaning and negative for the post-weaning datasets but negative with helminth eggs counts for the pre-weaning dataset and slightly positive for the post-weaning dataset. Genetic correlations between tick and helminth eggs counts were moderate and positive for both datasets. Phenotypic correlations of IBK incidence with helminth eggs per gram of faeces were moderate and positive for both datasets, but were close to zero for both datasets with tick counts. CONCLUSIONS: Our results suggest that genetic selection against IBK incidence in tropical cattle is feasible and that calves genetically prone to acquire IBK infection could also be genetically prone to have a slower growth. The positive genetic correlations among weight traits and between tick and helminth eggs counts suggest that they are controlled by common genes (with pleiotropic effects). Genetic correlations between IBK incidence and tick and helminth egg counts were moderate and opposite between pre-weaning and post-weaning datasets, suggesting that the environmental and (or) maternal effects differ between these two growth phases. This preliminary study provides estimated genetic parameters for IBK incidence, which could be used to design selection and breeding programs for tropical adaptation in beef cattle.

Host Genome-Metagenome Analyses Using Combinatorial Network Methods Reveal Key Metagenomic and Host Genetic Features for Methane Emission and Feed Efficiency in Cattle.

Cattle production is one of the key contributors to global warming due to methane emission, which is a by-product of converting feed stuff into milk and meat for human consumption. Rumen hosts numerous microbial communities that are involved in the digestive process, leading to notable amounts of methane emission. The key factors underlying differences in methane emission between individual animals are due to, among other factors, both specific enrichments of certain microbial communities and host genetic factors that influence the microbial abundances. The detection of such factors involves various biostatistical and bioinformatics methods. In this study, our main objective was to reanalyze a publicly available data set using our proprietary Synomics Insights platform that is based on novel combinatorial network and machine learning methods to detect key metagenomic and host genetic features for methane emission and residual feed intake (RFI) in dairy cattle. The other objective was to compare the results with publicly available standard tools, such as those found in the microbiome bioinformatics platform QIIME2 and classic GWAS analysis. The data set used was publicly available and comprised 1,016 dairy cows with 16S short read sequencing data from two dairy cow breeds: Holstein and Nordic Reds. Host genomic data consisted of both 50 k and 150 k SNP arrays. Although several traits were analyzed by the original authors, here, we considered only methane emission as key phenotype for associating microbial communities and host genetic factors. The Synomics Insights platform is based on combinatorial methods that can identify taxa that are differentially abundant between animals showing high or low methane emission or RFI. Focusing exclusively on enriched taxa, for methane emission, the study identified 26 order-level taxa that combinatorial networks reported as significantly enriched either in high or low emitters. Additionally, a Z-test on proportions found 21/26 (81%) of these taxa were differentially enriched between high and low emitters (p value <.05). In particular, the phylum of Proteobacteria and the order Desulfovibrionales were found enriched in high emitters while the order Veillonellales was found to be more abundant in low emitters as previously reported for cattle (Wallace et al., 2015). In comparison, using the publicly available tool ANCOM only the order Methanosarcinales could be identified as differentially abundant between the two groups. We also investigated a link between host genome and rumen microbiome by applying our Synomics Insights platform and comparing it with an industry standard GWAS method. This resulted in the identification of genetic determinants in cows that are associated with changes in heritable components of the rumen microbiome. Only four key SNPs were found by both our platform and GWAS, whereas the Synomics Insights platform identified 1,290 significant SNPs that were not found by GWAS. Gene Ontology (GO) analysis found transcription factor as the dominant biological function. We estimated heritability of a core 73 taxa from the original set of 150 core order-level taxonomies and showed that some species are medium to highly heritable (0.25-0.62), paving the way for selective breeding of animals with desirable core microbiome characteristics. We identified a set of 113 key SNPs associated with >90% of these core heritable taxonomies. Finally, we have characterized a small set (<10) of SNPs strongly associated with key heritable bacterial orders with known role in methanogenesis, such as Desulfobacterales and Methanobacteriales.

Genetic architecture of gene expression in ovine skeletal muscle.

BACKGROUND: In livestock populations the genetic contribution to muscling is intensively monitored in the progeny of industry sires and used as a tool in selective breeding programs. The genes and pathways conferring this genetic merit are largely undefined. Genetic variation within a population has potential, amongst other mechanisms, to alter gene expression via cis- or trans-acting mechanisms in a manner that impacts the functional activities of specific pathways that contribute to muscling traits. By integrating sire-based genetic merit information for a muscling trait with progeny-based gene expression data we directly tested the hypothesis that there is genetic structure in the gene expression program in ovine skeletal muscle. RESULTS: The genetic performance of six sires for a well defined muscling trait, longissimus lumborum muscle depth, was measured using extensive progeny testing and expressed as an Estimated Breeding Value by comparison with contemporary sires. Microarray gene expression data were obtained for longissimus lumborum samples taken from forty progeny of the six sires (4-8 progeny/sire). Initial unsupervised hierarchical clustering analysis revealed strong genetic architecture to the gene expression data, which also discriminated the sire-based Estimated Breeding Value for the trait. An integrated systems biology approach was then used to identify the major functional pathways contributing to the genetics of enhanced muscling by using both Estimated Breeding Value weighted gene co-expression network analysis and a differential gene co-expression network analysis. The modules of genes revealed by these analyses were enriched for a number of functional terms summarised as muscle sarcomere organisation and development, protein catabolism (proteosome), RNA processing, mitochondrial function and transcriptional regulation. CONCLUSIONS: This study has revealed strong genetic structure in the gene expression program within ovine longissimus lumborum muscle. The balance between muscle protein synthesis, at the levels of both transcription and translation control, and protein catabolism mediated by regulated proteolysis is likely to be the primary determinant of the genetic merit for the muscling trait in this sheep population. There is also evidence that high genetic merit for muscling is associated with a fibre type shift toward fast glycolytic fibres. This study provides insight into mechanisms, presumably subject to strong artificial selection, that underpin enhanced muscling in sheep populations.