Profile of the Spleen Transcriptome in Beef Steers with Variation in Gain and Feed Intake

We have previously identified components of the immune system contributing to feed intake and gain in both the rumen and small intestine of beef steers. In this study, we examined the spleen, a major lymphatic organ near the digestive tract, to determine whether it was also influencing individual feed efficiency status through immune responses. Animals (n=16) that were divergent for gain and intake were selected for tissue sampling. The spleen transcriptomes were evaluated by microarray. A total of 1,216 genes were identified as differentially expressed. Genes were over-represented in Kyoto encyclopedia of genes and genomes (KEGG) pathways including biological regulation, protein folding, cell communication, immune systems process, response to stress and RNA metabolic process. Several stress response or heat shock genes including HSPH1, HSPA1A, HSPA4, DNAJB4, DNAJA4, etc., were identified as a stress response functional gene cluster in the low gain-low intake animals. These genes were up-regulated amongst the low gain-low intake animals compared to all other groups. Canonical pathways associated with the differentially expressed genes included the coagulation system, extrinsic prothrombin activation, protein ubiquitination, unfolded protein response and aldosterone signaling in epithelial cells. An analysis of expressed copy number variable (CNV) genes in the spleen produced some of the same genes and gene families that were differentially expressed. Our data suggests the splenic contribution to some of the underlying variation among gain and intake within this group of animals may be a result of immune function and stress response. In addition, some of the differences in immune response functions may be related to gene copy number

Genes Involved in Feed Efficiency Identified in a Meta-Analysis of Rumen Tissue from Two Populations of Beef Steers

In cattle, the rumen is an important site for the absorption of feed by-products released by bacterial fermentation, and variation in ruminal function plays a role in cattle feed efficiency. Studies evaluating gene expression in the rumen tissue have been performed prior to this. However, validating the expression of genes identified in additional cattle populations has been challenging. The purpose of this study was to perform a meta-analysis of the ruminal transcriptome of two unrelated populations of animals to identify genes that are involved in feed efficiency across populations. RNAseq data from animals with high and low residual feed intake (RFI) from a United States population of cattle (eight high and eight low RFI) and a Canadian population of cattle (nine high and nine low RFI) were analyzed for differences in gene expression. A total of 83 differentially expressed genes were identified. Some of these genes have been previously identified in other feed efficiency studies. These genes included ATP6AP1, BAG6, RHOG, and YPEL3. Differentially expressed genes involved in the Notch signaling pathway and in protein turnover were also identified. This study, combining two unrelated populations of cattle in a meta-analysis, produced several candidate genes for feed efficiency that may be more robust indicators of feed efficiency than those identified from single populations of animals

Evaluation of genotype quality parameters for \u3ci\u3eSowPro90\u3c/i\u3e, a new genotyping array for swine

Understanding early predictors of sow fertility has the potential to improve genomic predictions. A custom SNP array (SowPro90 produced by Affymetrix) was developed to include genetic variants overlapping quantitative trait loci for age at puberty, one of the earliest indicators of sow fertility, as well as variants related to innate and adaptive immunity. The polymorphisms included in the custom genotyping array were identified using multiple genomic approaches including deep genomic and transcriptomic sequencing and genome-wide associations. Animals from research and commercial populations (n = 2,586) were genotyped for 103,476 SNPs included in SowPro90. To assess the quality of data generated, genotype concordance was evaluated between the SowPro90 and Porcine SNP60 BeadArray using a subset of common SNP (n = 44,708) and animals (n = 277). The mean genotype concordance rate per SNP was 98.4%. Differences in distribution of data quality were observed between the platforms indicating the need for platform specific thresholds for quality parameters. The optimal thresholds for SowPro90 (≥97% SNP and ≥93% sample call rate) were obtained by analyzing the data quality distribution and genotype concordance per SNP across platforms. At ≥97% SNP call rate, there were 42,151 SNPs (94.3%) retained with a mean genotype concordance of 98.6% across platforms. Similarly, ≥94% SNPs and ≥85% sample call rates were established as thresholds for Porcine SNP60 BeadArray. At ≥94% SNPs call rate, there were 41,043 SNPs (91.8%) retained with a mean genotype concordance of 98.6% across platforms. Final evaluation of SowPro90 array content (n = 103,476) at ≥97% SNPs and ≥93% sample call rates allowed retention of 89,040 SNPs (86%) for downstream analysis. The findings and strategy for quality control could be helpful in identifying consistent, high-quality genotypes for genomic evaluations, especially when integrating genotype data from different platforms

Neonatal lactocrine deficiency affects the adult porcine endometrial transcriptome at pregnancy day 13

Reproductive performance of female pigs that do not receive sufficient colostrum from birth is permanently impaired. Whether lactocrine deficiency, reflected by low serum immunoglobulin immunocrit (iCrit), affects patterns of endometrial gene expression during the periattachment period of early pregnancy is unknown. Here, objectives were to determine effects of low iCrit at birth on the adult endometrial transcriptome on pregnancy day (PxD) 13. On the first day of postnatal life, gilts were assigned to high or low iCrit groups. Adult high (n = 8) and low (n = 7) iCrit gilts were bred (PxD 0), and humanely slaughtered on PxD 13 when tissues and fluids were collected. The endometrial transcriptome was defined for each group using mRNAseq and microRNAseq. Reads were mapped to the Sus scrofa 11.1 genome build. Mature microRNAs were annotated using miRBase 21. Differential expression was defined based on fold change (≥ ±1.5). Lactocrine deficiency did not affect corpora lutea number, uterine horn length, uterine wet weight, conceptus recovery, or uterine luminal fluid estrogen content on PxD 13. However, mRNAseq revealed 1157 differentially expressed endometrial mRNAs in high versus low iCrit gilts. Differentially expressed genes had functions related to solute transport, endometrial receptivity, and immune response. Six differentially expressed endometrial microRNAs included five predicted to target 62 differentially expressed mRNAs, affecting similar biological processes. Thus, lactocrine deficiency on the first day of postnatal life can alter uterine developmental trajectory with lasting effects on endometrial responses to pregnancy as reflected at the level of the transcriptome on PxD 13

Neonatal lactocrine deficiency affects the adult porcine endometrial transcriptome at pregnancy day 13

Reproductive performance of female pigs that do not receive sufficient colostrum from birth is permanently impaired. Whether lactocrine deficiency, reflected by low serum immunoglobulin immunocrit (iCrit), affects patterns of endometrial gene expression during the periattachment period of early pregnancy is unknown. Here, objectives were to determine effects of low iCrit at birth on the adult endometrial transcriptome on pregnancy day (PxD) 13. On the first day of postnatal life, gilts were assigned to high or low iCrit groups. Adult high (n = 8) and low (n = 7) iCrit gilts were bred (PxD 0), and humanely slaughtered on PxD 13 when tissues and fluids were collected. The endometrial transcriptome was defined for each group using mRNAseq and microRNAseq. Reads were mapped to the Sus scrofa 11.1 genome build. Mature microRNAs were annotated using miRBase 21. Differential expression was defined based on fold change (≥ ±1.5). Lactocrine deficiency did not affect corpora lutea number, uterine horn length, uterine wet weight, conceptus recovery, or uterine luminal fluid estrogen content on PxD 13. However, mRNAseq revealed 1157 differentially expressed endometrial mRNAs in high versus low iCrit gilts. Differentially expressed genes had functions related to solute transport, endometrial receptivity, and immune response. Six differentially expressed endometrial microRNAs included five predicted to target 62 differentially expressed mRNAs, affecting similar biological processes. Thus, lactocrine deficiency on the first day of postnatal life can alter uterine developmental trajectory with lasting effects on endometrial responses to pregnancy as reflected at the level of the transcriptome on PxD 13

Global analysis of differential gene expression within the porcine conceptus transcriptome as it transitions through spherical, ovoid, and tubular morphologies during the initiation of elongation

This study aimed to identify transcriptome differences between distinct or transitional stage spherical, ovoid, and tubular porcine blastocysts throughout the initiation of elongation. We performed a global transcriptome analysis of differential gene expression using RNA‐Seq with high temporal resolution between spherical, ovoid, and tubular stage blastocysts at specific sequential stages of development from litters containing conceptus populations of distinct or transitional blastocysts. After RNA‐Seq analysis, significant differentially expressed genes (DEGs) and pathways were identified between distinct morphologies or sequential development stages. Overall, 1898 significant DEGs were identified between distinct spherical and ovoid morphologies, with 311 total DEGs between developmental stages throughout this first morphological transition, while 15 were identified between distinct ovoid and tubular, with eight total throughout these second morphological transition developmental stages. The high quantity of DEGs and pathways between conceptus stages throughout the spherical to ovoid transition suggests the importance of gene regulation during this first morphological transition for initiating elongation. Further, extensive DEG coverage of known elongation signaling pathways was illustrated from spherical to ovoid, and regulation of lipid signaling and membrane/ECM remodeling across these early conceptus stages were implicated as essential to this process, providing novel insights into potential mechanisms governing this rapid morphological change

GENOMICS SYMPOSIUM: Using genomic approaches to uncover sources of variation in age at puberty and reproductive longevity in sows

Genetic variants associated with traits such as age at puberty and litter size could provide insight into the underlying genetic sources of variation impacting sow reproductive longevity and productivity. Genomewide characterization and gene expression profiling were used using gilts from the University of Nebraska–Lincoln swine resource population (n = 1,644) to identify genetic variants associated with age at puberty and litter size traits. From all reproductive traits studied, the largest fraction of phenotypic variation explained by the Porcine SNP60 BeadArray was for age at puberty (27.3%). In an evaluation data set, the predictive ability of all SNP from highranked 1-Mb windows (1 to 50%), based on genetic variance explained in training, was greater (12.3 to 36.8%) compared with the most informative SNP from these windows (6.5 to 23.7%). In the integrated data set (n = 1,644), the top 1% of the 1-Mb windows explained 6.7% of the genetic variation of age at puberty. One of the high-ranked windows detected (SSC2, 12–12.9 Mb) showed pleiotropic features, affecting both age at puberty and litter size traits. The RNA sequencing of the hypothalami arcuate nucleus uncovered 17 differentially expressed genes (adjusted P \u3c 0.05) between gilts that became pubertal early (180 d of age). Twelve of the differentially expressed genes are upregulated in the late pubertal gilts. One of these genes is involved in energy homeostasis (FFAR2), a function in which the arcuate nucleus plays an important contribution, linking nutrition with reproductive development. Energy restriction during the gilt development period delayed age at puberty by 7 d but increased the probability of a sow to produce up to 3 parities (P \u3c 0.05). Identification of pleotropic functional polymorphisms may improve accuracy of genomic prediction while facilitating a reduction in sow replacement rates and addressing welfare concerns

Transcriptomic Analysis for Pork Color—The Ham Halo Effect in Biceps Femoris

Pork color is a major indicator of product quality that guides consumer purchasing decisions. Recently, industryhas received an increase in consumer complaints about the lightness and nonuniformity of ham color, primarily lighter colorin the periphery termed “ham halo” that is not caused by manufacturing procedures. This effect is seen in fresh and processed hams and the outer lighter muscle is associated with lower myoglobin concentration, pH, and type I fibers. The objective of this study was to identify differences in gene expression profiles between light and normal-colored portions of the biceps femoris muscle from pork hams. RNA-sequencing (RNA-seq) was performed for paired light and normal-colored muscle samples from 10 animals showing the ham halo effect. Over 50 million paired-end reads (2×75 bp) per library were obtained. An average of 99.74% of trimmed high-quality reads was mapped to the Sscrofa 11.1 genome assembly. Differentially expressed genes (DEGs) were identified using both the DESeq2 and GFOLD software packages. A total of 14,049 genes were expressed in the biceps femoris; 13,907 were expressed in both light and normal muscle, while 56 and 86 genes were only expressed in light and normal muscle, respectively. Analysis with DESeq2 identified 392 DEGs with 359 genes being more highly expressed in normal-colored muscle. A total of 61 DEGs were identified in the DESeq2 analysis and identified in at least 7 of the 10 individual animal analyses. All 61 of these DEGs were up-regulated in normal-colored muscle. Gene Ontology enrichment analysis of DEGs identified the transition between fast and slow fibers and skeletal muscle adaptation and contraction as the most significant biological process terms. The evaluation of gene expression by RNA-seq identified DEGs between regions of the biceps femoris with the ham halo effect that are associated with variation in pork color

The roles of age at puberty and energy restriction |in sow reproductive longevity: a genomic perspective

Approximately 50% of sows are culled annually with more than one-third due to poor fertility. Our research demonstrated that age at puberty is an early pre-breeding indicator of reproductive longevity. Age at puberty can be measured early in life, has a moderate heritability, and is negatively correlated with lifetime number of parities. Detection of age at puberty is tedious and time consuming and is therefore not collected by the industry, which limits genetic progress. Genomic prediction is a viable approach to preselect gilts that will express puberty early and have superior reproductive longevity. The hypothesis that genetic variants explaining differences in age at puberty also explain differences in sow reproductive longevity was tested. Phenotypes, genotypes, and tissues from the UNL resource population (n \u3e 1700) were used in genome-wide association analyses, genome, and RNA sequencing to uncover functional polymorphisms that could explain variation in puberty and reproductive longevity. A BeadArray including 56,424 SNP explained 25.2% of the phenotypic variation in age at puberty in a training set (n = 820). Evaluation of major windows and SNPs of subsequent batches of similar genetics (n = 412) showed that if all SNPs located in the major 1-Mb windows were tested, they explained a substantial amount of phenotypic variation (12.3 to 36.8%). Due to differences in linkage disequilibrium status, the most informative SNP from these windows explained a lower proportion of the variation (6.5 to 23.7%). To improve genomic predictive ability, the limited capability of BeadArray was enhanced by potential functional variants uncovered by genome sequencing of selected sires (n = 20; \u3e20X). There were 11.2 mil. SNPs and 2.9 mil. indels discovered across sires and reference genomes. The role of gene expression differences in explaining phenotypic variation in age at puberty was investigated by RNA sequencing of the hypothalamic arcuate nucleus (ARC) in gilts (n = 37) with different pubertal statuses. Seventy genes, including genes involved in reproductive processes, were differentially expressed between gilts with early and late puberty status (Padj \u3c 0.1). Dietary restriction of energy 3 mo before breeding delayed puberty by 7 d but improved the potential of a sow producing up to three parities (P \u3c 0.05). Energy restriction was associated with differential expression in 42 genes in the ARC, including genes involved in energy metabolism. This integrated genomic information will be evaluated in commercial populations to improve the reproductive potential of sows through genomic selection. This project is supported by AFRI Competitive grant no. 2013-68004-20370 from the USDA-NIFA. USDA is an equal opportunity provider and employer

A compendium of genetic regulatory effects across pig tissues

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.</p