Brittney N. Keel

Genome-wide copy number variation in the bovine genome detected using low coverage sequence of popular beef breeds.

Copy number variations (CNVs) are large insertions, deletions or duplications in the genome that vary between members of a species and are known to affect a wide variety of phenotypic traits. In this study, we identified CNVs in a population of bulls using low coverage next-generation sequence data. First, in order to determine a suitable strategy for CNV detection in our data, we compared the performance of three distinct CNV detection algorithms on benchmark CNV datasets and concluded that using the multiple sample read depth approach was the best method for identifying CNVs in our sequences. Using this technique, we identified a total of 1341 copy number variable regions (CNVRs) from genome sequences of 154 purebred sires used in Cycle VII of the USMARC Germplasm Evaluation Project. These bulls represented the seven most popular beef breeds in the United States: Hereford, Charolais, Angus, Red Angus, Simmental, Gelbvieh and Limousin. The CNVRs covered 6.7% of the bovine genome and spanned 2465 protein-coding genes and many known quantitative trait loci (QTL). Genes harbored in the CNVRs were further analyzed to determine their function as well as to find any breed-specific differences that may shed light on breed differences in adaptation, health and production.

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 1216 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.

Evolutionary and Functional Features of Copy Number Variation in the Cattle Genome

Genomic structural variations are an important source of genetic diversity. Copy number variations (CNVs), gains and losses of large regions of genomic sequence between individuals of a species, have been associated with a wide variety of phenotypic traits. However, in cattle, as well as many other species, relatively little is understood about CNV, including frequency of CNVs in the genome, sizes, and locations, chromosomal properties, and evolutionary processes acting to shape CNV. In this work, we focused on copy number variation in the bovine genome, with the aim to detect CNVs in Bos taurus coding sequence and explore potential evolutionary mechanisms shaping these CNV. We identified and characterized CNV regions by utilizing exome sequence from 175 influential sires used in the Germplasm Evaluation project, representing 10 breeds. We examined various evolutionary and functional aspects of these CNVs, including selective constraint on CNV-overlapped genes, centrality of CNV genes in protein-protein interaction networks, and tissue-specific expression of CNV genes. Patterns of CNV in the Bos taurus genome reveal that reduced functional constraint and mutational bias may play a prominent role in shaping this type of structural variation.

Genome-Wide Association of Myoglobin Concentrations in Pork Loins

Lean color is a major focus for identifying pork loins for export markets, and myoglobin is the primary pigment driving pork color. Thus, increasing myoglobin concentration should increase redness of pork products and the number of loins acceptable for exportation. Therefore, understanding genetic variation and parameters affecting myoglobin concentration is critical for improving pork color. The objective of this study was to identify genetic markers associated with myoglobin concentration in pork loin muscle. Ultimate pH and myoglobin concentrations were measured in longissimus thoracis et lumborum samples of pigs (n = 599) from two different commercial finishing swine facilities. A Bayes-C model implemented in GenSel identified regions within 7 chromosomes that explained greater than 63% of the genetic variance in myoglobin concentration. Chromosome 7 had 1 significant region which accounted for 37% of the genetic variance, while chromosome 14 had 4 significant regions accounting for 9.8% of the genetic variance. Candidate genes in the region on chromosome 7 were involved in iron homeostasis, and genes in the significant regions on chromosome 14 were involved in calcium regulation. Genes identified in this study represent potential biomarkers that could be used to select for higher myoglobin concentrations in pork, which may improve lean meat color.

The Wnt signaling pathway is differentially expressed during the bovine herpesvirus 1 latency-reactivation cycle: evidence that two protein kinases associated with neuronal survival (Akt3 and bone morphogenetic protein receptor 2) are expressed at higher levels during latency.

ABSTRACT Sensory neurons in trigeminal ganglia (TG) of calves latently infected with bovine herpesvirus 1 (BoHV-1) abundantly express latency-related (LR) gene products, including a protein (ORF2) and two micro-RNAs. Recent studies in mouse neuroblastoma cells (Neuro-2A) demonstrated ORF2 interacts with β-catenin and a β-catenin coactivator, high-mobility group AT-hook 1 (HMGA1) protein, which correlates with increased β-catenin-dependent transcription and cell survival. β-Catenin and HMGA1 are readily detected in a subset of latently infected TG neurons but not TG neurons from uninfected calves or reactivation from latency. Consequently, we hypothesized that the Wnt/β-catenin signaling pathway is differentially expressed during the latency and reactivation cycle and an active Wnt pathway promotes latency. RNA-sequencing studies revealed that 102 genes associated with the Wnt/β-catenin signaling pathway were differentially expressed in TG during the latency-reactivation cycle in calves. Wnt agonists were generally expressed at higher levels during latency, but these levels decreased during dexamethasone-induced reactivation. The Wnt agonist bone morphogenetic protein receptor 2 (BMPR2) was intriguing because it encodes a serine/threonine receptor kinase that promotes neuronal differentiation and inhibits cell death. Another differentially expressed gene encodes a protein kinase (Akt3), which is significant because Akt activity enhances cell survival and is linked to herpes simplex virus 1 latency and neuronal survival. Additional studies demonstrated ORF2 increased Akt3 steady-state protein levels and interacted with Akt3 in transfected Neuro-2A cells, which correlated with Akt3 activation. Conversely, expression of Wnt antagonists increased during reactivation from latency. Collectively, these studies suggest Wnt signaling cooperates with LR gene products, in particular ORF2, to promote latency. IMPORTANCE Lifelong BoHV-1 latency primarily occurs in sensory neurons. The synthetic corticosteroid dexamethasone consistently induces reactivation from latency in calves. RNA sequencing studies revealed 102 genes associated with the Wnt/β-catenin signaling pathway are differentially regulated during the latency-reactivation cycle. Two protein kinases associated with the Wnt pathway, Akt3 and BMPR2, were expressed at higher levels during latency but were repressed during reactivation. Furthermore, five genes encoding soluble Wnt antagonists and β-catenin-dependent transcription inhibitors were induced during reactivation from latency. These findings are important because Wnt, BMPR2, and Akt3 promote neurogenesis and cell survival, processes crucial for lifelong viral latency. In transfected neuroblastoma cells, a viral protein expressed during latency (ORF2) interacts with and enhances Akt3 protein kinase activity. These findings provide insight into how cellular factors associated with the Wnt signaling pathway cooperate with LR gene products to regulate the BoHV-1 latency-reactivation cycle.

Genomewide association study of liver abscess in beef cattle.

Fourteen percent of U.S. cattle slaughtered in 2011 had liver abscesses, resulting in reduced carcass weight, quality, and value. Liver abscesses can result from a common bacterial cause, , which inhabits rumen lesions caused by acidosis and subsequently escapes into the blood stream, is filtered by the liver, and causes abscesses in the liver. Our aim was to identify SNP associated with liver abscesses in beef cattle. We used lung samples as a DNA source because they have low economic value, they have abundant DNA, and we had unrestricted access to sample them. We collected 2,304 lung samples from a beef processing plant: 1,152 from animals with liver abscess and 1,152 from animals without liver abscess. Lung tissue from pairs of animals, 1 with abscesses and another without, were collected from near one another on the viscera table to ensure that pairs of phenotypically extreme animals came from the same lot. Within each phenotype (abscess or no abscess), cattle were pooled by slaughter sequence into 12 pools of 96 cattle for each phenotype for a total of 24 pools. The pools were constructed by equal volume of frozen lung tissue from each animal. The DNA needed to allelotype each pool was then extracted from pooled lung tissue and the BovineHD Bead Array (777,962 SNP) was run on all 24 pools. Total intensity (TI), an indicator of copy number variants, was the sum of intensities from red and green dyes. Pooling allele frequency (PAF) was red dye intensity divided TI. Total intensity and PAF were weighted by the inverse of their respective genomic covariance matrices computed over all SNP across the genome. A false discovery rate ≤ 5% was achieved for 15 SNP for PAF and 20 SNP for TI. Genes within 50 kbp from significant SNP were in diverse pathways including maintenance of pH homeostasis in the gastrointestinal tract, maintain immune defenses in the liver, migration of leukocytes from the blood into infected tissues, transport of glutamine into the kidney in response to acidosis to facilitate production of bicarbonate to increase pH, aggregate platelets to liver injury to facilitate liver repair, and facilitate axon guidance. Evidence from the 35 detected SNP associations combined with evidence of polygenic variation indicate that there is adequate genetic variation in incidence rate of liver abscesses, which could be exploited to select sires for reduced susceptibility to subacute acidosis and associated liver abscess.

The effects of Capn1 gene inactivation on the differential expression of genes in skeletal muscle

Protein turnover is required for muscle growth and regeneration and several proteolytic enzymes, including the calpains, degrade myofibrillar proteins during this process. In a previous experiment, phenotypic differences were observed between μ-calpain knockout (KO) and wild type (WT) mice, including nutrient accretion and fiber type differences. These changes were particularly evident as the animals aged. Thus, we utilized 18 mice (9 KO and 9 WT) to compare transcript abundance to identify differentially expressed genes (DEGs) at 52 wk of age. A total of 55 genes were differentially expressed, including adiponectin, phosphoenolpyruvate carboxykinase 1, uncoupling protein 1, and lysine deficient protein kinase 2. These genes were analyzed for over- and underrepresented gene ontology (GO) terms. Several GO terms, including response to cytokine, response to interferon-beta, regulation of protein phosphorylation, and hydrolase activity, were identified as overrepresented. Pathways related to taurine biosynthesis, nitric oxide synthase signaling, amyloid processing, and L-cysteine degradation were also identified. Our results are consistent with previous experiments, in that identified DEGs may explain, at least in part, some of the phenotypic differences between μ-calpain KO and WT mice. Clearly muscle growth and maintenance are complex, multifaceted processes. Genes affected by the silencing of the μ-calpain gene have been identified, but the relationship between μ-calpain and these pathways requires further investigation.

Comparison of Burrows-Wheeler Transform-Based Mapping Algorithms Used in High-Throughput Whole-Genome Sequencing: Application to Illumina Data for Livestock Genomes1

Ongoing developments and cost decreases in next-generation sequencing (NGS) technologies have led to an increase in their application, which has greatly enhanced the fields of genetics and genomics. Mapping sequence reads onto a reference genome is a fundamental step in the analysis of NGS data. Efficient alignment of the reads onto the reference genome with high accuracy is very important because it determines the global quality of downstream analyses. In this study, we evaluate the performance of three Burrows-Wheeler transform-based mappers, BWA, Bowtie2, and HISAT2, in the context of paired-end Illumina whole-genome sequencing of livestock, using simulated sequence data sets with varying sequence read lengths, insert sizes, and levels of genomic coverage, as well as five real data sets. The mappers were evaluated based on two criteria, computational resource/time requirements and robustness of mapping. Our results show that BWA and Bowtie2 tend to be more robust than HISAT2, while HISAT2 was significantly faster and used less memory than both BWA and Bowtie2. We conclude that there is not a single mapper that is ideal in all scenarios but rather the choice of alignment tool should be driven by the application and sequencing technology.

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

Abstract 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. Summary Sentence Lactocrine deficiency on the first day of postnatal life alters the uterine developmental program with long-term effects on patterns of porcine endometrial gene expression during the periattachment period of early pregnancy.

RNA-Seq Meta-analysis identifies genes in skeletal muscle associated with gain and intake across a multi-season study of crossbred beef steers

BackgroundFeed intake and body weight gain are economically important inputs and outputs of beef production systems. The purpose of this study was to discover differentially expressed genes that will be robust for feed intake and gain across a large segment of the cattle industry. Transcriptomic studies often suffer from issues with reproducibility and cross-validation. One way to improve reproducibility is by integrating multiple datasets via meta-analysis. RNA sequencing (RNA-Seq) was performed on longissimus dorsi muscle from 80 steers (5 cohorts, each with 16 animals) selected from the outside fringe of a bivariate gain and feed intake distribution to understand the genes and pathways involved in feed efficiency. In each cohort, 16 steers were selected from one of four gain and feed intake phenotypes (n = 4 per phenotype) in a 2 × 2 factorial arrangement with gain and feed intake as main effect variables. Each cohort was analyzed as a single experiment using a generalized linear model and results from the 5 cohort analyses were combined in a meta-analysis to identify differentially expressed genes (DEG) across the cohorts.ResultsA total of 51 genes were differentially expressed for the main effect of gain, 109 genes for the intake main effect, and 11 genes for the gain x intake interaction (Pcorrected < 0.05). A jackknife sensitivity analysis showed that, in general, the meta-analysis produced robust DEGs for the two main effects and their interaction. Pathways identified from over-represented genes included mitochondrial energy production and oxidative stress pathways for the main effect of gain due to DEG including GPD1, NDUFA6, UQCRQ, ACTC1, and MGST3. For intake, metabolic pathways including amino acid biosynthesis and degradation were identified, and for the interaction analysis the pathways identified included GADD45, pyridoxal 5’phosphate salvage, and caveolar mediated endocytosis signaling.ConclusionsVariation among DEG identified by cohort suggests that environment and breed may play large roles in the expression of genes associated with feed efficiency in the muscle of beef cattle. Meta-analyses of transcriptome data from groups of animals over multiple cohorts may be necessary to elucidate the genetics contributing these types of biological phenotypes.