The Anti-Müllerian Hormone as Endocrine and Molecular Marker Associated with Reproductive Performance in Holstein Dairy Cows Exposed to Heat Stress.

Anti-Müllerian hormone (AMH) is proposed as a biomarker for fertility in cattle, yet this associative relationship appears to be influenced by heat stress (HS). The objective was to test serum AMH and AMH-related single nucleotide polymorphisms (SNPs) as markers potentially predictive of reproductive traits in dairy cows experiencing HS. The study included 300 Holstein cows that were genotyped using BovineSNP50 (54,000 SNP). A genome-wide association study was then executed. Nine intragenic SNPs within the pathways that influence the AMH gene were found important with multiple comparisons adjustment tests (p < 1.09 × 10-6). A further validation study was performed in an independent Holstein cattle population, which was divided into moderate (MH; n = 152) and severe heat-stressed (SH; n = 128) groups and then subjected to a summer reproductive management program. Serum AMH was confirmed as a predictor of fertility measures (p < 0.05) in MH but not in the SH group. Cows were genotyped, which revealed four SNPs as predictive markers for serum AMH (p < 0.01), reproductive traits (p < 0.01), and additional physiological variables (p < 0.05). These SNPs were in the genes AMH, IGFBP1, LGR5, and TLR4. In conclusion, serum AMH concentrations and AMH polymorphisms are proposed as predictive markers that can be used in conjunction with genomic breeding value approaches to improve reproductive performance in Holstein cows exposed to summer HS conditions.

Novel lncRNA regulatory elements in milk somatic cells of Holstein dairy cows associated with mastitis.

Despite regulatory elements such as long non - coding RNAs representing most of the transcriptome, the functional understanding of long non - coding RNAs in relation to major health conditions including bovine mastitis is limited. This study examined the milk somatic cell transcriptome from udder quarters of 6 Holstein dairy cows to identify differentially expressed long non - coding RNAs using RNA - Sequencing. Ninety - four differentially expressed long non - coding RNAs are identified, 5 of which are previously annotated for gene name and length, 11 are annotated for gene name and 78 are novel, having no gene name or length previously annotated. Significant inflammatory response and regulation of immune response pathways (false discovery rate < 0.05) are associated with the differentially expressed long non - coding RNAs. QTL annotation analysis revealed 31 QTL previously annotated in the genomic regions of the 94 differentially expressed long non - coding RNAs, and the majority are associated with milk traits. This research provides a better understanding of long non - coding RNAs regulatory elements in milk somatic cells, which may enhance current breeding strategies for more adaptable or high mastitis resistant cattle.

Genetic Markers Associated with Milk Production and Thermotolerance in Holstein Dairy Cows Managed in a Heat-Stressed Environment.

Dairy production in Holstein cows in a semiarid environment is challenging due to heat stress. Under such conditions, genetic selection for heat tolerance appears to be a useful strategy. The objective was to validate molecular markers associated with milk production and thermotolerance traits in Holstein cows managed in a hot and humid environment. Lactating cows (n = 300) exposed to a heat stress environment were genotyped using a medium-density array including 53,218 SNPs. A genome-wide association study (GWAS) detected six SNPs associated with total milk yield (MY305) that surpassed multiple testing (p < 1.14 × 10-6). These SNPs were further validated in 216 Holstein cows from two independent populations that were genotyped using the TaqMan bi-allelic discrimination method and qPCR. In these cows, only the SNPs rs8193046, rs43410971, and rs382039214, within the genes TLR4, GRM8, and SMAD3, respectively, were associated (p < 0.05) with MY305, rectal temperature (RT), and respiratory rate. Interestingly, these variables improved as the number of favorable genotypes of the SNPs increased from 0 to 3. In addition, a regression analysis detected RT as a significant predictor (R2 = 0.362) for MY305 in cows with >1 favorable genotype, suggesting this close relationship was influenced by genetic markers. In conclusion, SNPs in the genes TLR4, GRM8, and SMAD3 appear to be involved in the molecular mechanism that regulates milk production in cows under heat-stressed conditions. These SNPs are proposed as thermotolerance genetic markers for a selection program to improve the milk performance of lactating Holstein cows managed in a semiarid environment.

Single-Step Genome Editing of Small Ruminant Embryos by Electroporation.

We investigated the possibility of single-step genome editing in small ruminants by CRISPR-Cas9 zygote electroporation. We targeted SOCS2 and PDX1 in sheep embryos and OTX2 in goat embryos, utilizing a dual sgRNA approach. Gene editing efficiency was compared between microinjection and three different electroporation settings performed at four different times of embryo development. Electroporation of sheep zygotes 6 h after fertilization with settings that included short high-voltage (poring) and long low-voltage (transfer) pulses was efficient at producing SOCS2 knock-out blastocysts. The mutation rate after CRISPR/Cas9 electroporation was 95.6% ± 8%, including 95.4% ± 9% biallelic mutations; which compared favorably to 82.3% ± 8% and 25% ± 10%, respectively, when using microinjection. We also successfully disrupted the PDX1 gene in sheep and the OTX2 gene in goat embryos. The biallelic mutation rate was 81 ± 5% for PDX1 and 85% ± 6% for OTX2. In conclusion, using single-step CRISPR-Cas9 zygote electroporation, we successfully introduced biallelic deletions in the genome of small ruminant embryos.

Effect of calving season on the parameters and components of the lactation curve in Holstein dairy cows managed in a semi-desert climate.

The lactation curve in dairy cows is influenced by the calving season, which is highly dependent on the warm climate in semi-arid regions. Objective herein was to evaluate effects of calving season on the parameters and components of the lactation curve in Holstein cows. The study included 278,317 milk records collected from 1086 cows from the 1st to 5th calving and good body condition score. The cows were grouped according to the season in which they calved: winter-calved (CS1), summer-calved (CS2), and autumn-calved cows (CS3). Ambient temperature and humidity data were used to calculate the temperature-humidity index (THI). The NLIN procedure was used to estimate the parameters of the lactation curve that served to calculate the components. The mixed procedure was executed to analyze the fixed effect of calving season. Associations between lactation curve traits were tested using correlation and regression analyses. A univariate model was utilized to calculate heritability. Average THI values during the lactation period were 73.5, 68.5, and 69.5 units for CS1, CS2, and CS3 groups, respectively. Initial milk production and increasing rate to the maximum milk yield in CS1 and CS3 groups were higher (P < 0.05) than CS2 cows. However, persistency and total milk yield during the entire lactation period were superior (P < 0.05) for CS2 and CS3 cows compared to CS1 cows, probably due to the moderate heat stress during the lactation period in the CS1 group. In cows from CS2 and CS3 groups, total milk production at 305 days was moderately correlated with initial milk production (r = 0.47; P < 0.05), and highly correlated with milk yield at peak day (r = 0.91; P < 0.05) which resulted as reliable predictor for total milk yield during the entire lactation (R2 = 0.83). In conclusion, the THI prevailing during the different calving seasons appeared to be an important factor influencing the performance of the lactation curve.

Large-Scale Multiplexing Permits Full-Length Transcriptome Annotation of 32 Bovine Tissues From a Single Nanopore Flow Cell.

A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome - or any mammalian transcriptome - appears significantly more feasible.

A SNP within the PMCH gene as a molecular marker associated with fertility traits in Angus and Brangus beef heifers raised under a desert environment.

Components of the GH/IGF1 endocrine axis regulate growth and reproductive traits in cattle. The pro-melanin-concentrating hormone (PMCH) gene located within chromosome 5 belongs to this axis. Objective herein was to evaluate PMCH single-nucleotide polymorphisms (SNPs) as molecular markers associated with age at first calving, calving interval, and age at second calving in Angus and Brangus beef heifers raised in desert conditions. Five SNPs within the PMCH gene were included in the study. Three of these SNPs had minor allele frequency > 10% and only one SNP did not deviate from Hardy-Weinberg equilibrium. A genotype to phenotype association analyses was performed using a mixed-effects model which included phenotype as the response variable, SNP genotype, breed, year of birth and age of dam as fixed terms, and sire as a random effect. Genotypes from the SNP rs135033882 were found to be associated (P < 0.05) with all evaluated fertility traits, and the term breed resulted as a significant source of variation only for age at second calving. The allele A was the favorable allele because it decreased the age at first calving 98.6 days, the calving interval 85.3 days, and the age at second calving 183.1 days, in Angus and Brangus heifers. In conclusion, we proposed a SNP within the PMCH gene as a potential candidate marker associated with reproductive performance in Angus and Brangus beef heifers raised in a desert climate.

Functional annotations of three domestic animal genomes provide vital resources for comparative and agricultural research.

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources.

A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues.

BACKGROUND: Although considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues. RESULTS: Overall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions. CONCLUSIONS: The lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

De novo assembly of the cattle reference genome with single-molecule sequencing.

BACKGROUND: Major advances in selection progress for cattle have been made following the introduction of genomic tools over the past 10-12 years. These tools depend upon the Bos taurus reference genome (UMD3.1.1), which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. RESULTS: We present the new reference genome for cattle, ARS-UCD1.2, based on the same animal as the original to facilitate transfer and interpretation of results obtained from the earlier version, but applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly includes 2.7 Gb and is >250× more continuous than the original assembly, with contig N50 >25 Mb and L50 of 32. We also greatly expanded supporting RNA-based data for annotation that identifies 30,396 total genes (21,039 protein coding). The new reference assembly is accessible in annotated form for public use. CONCLUSIONS: We demonstrate that improved continuity of assembled sequence warrants the adoption of ARS-UCD1.2 as the new cattle reference genome and that increased assembly accuracy will benefit future research on this species.

Deducing signaling pathways from parallel actions of arsenite and antimonite in human epidermal keratinocytes.

Inorganic arsenic oxides have been identified as carcinogens in several human tissues, including epidermis. Due to the chemical similarity between trivalent inorganic arsenic (arsenite) and antimony (antimonite), we hypothesized that common intracellular targets lead to similarities in cellular responses. Indeed, transcriptional and proteomic profiling revealed remarkable similarities in differentially expressed genes and proteins resulting from exposure of cultured human epidermal keratinocytes to arsenite and antimonite in contrast to comparisons of arsenite with other metal compounds. These data were analyzed to predict upstream regulators and affected signaling pathways following arsenite and antimonite treatments. A majority of the top findings in each category were identical after treatment with either compound. Inspection of the predicted upstream regulators led to previously unsuspected roles for oncostatin M, corticosteroids and ephrins in mediating cellular response. The influence of these predicted mediators was then experimentally verified. Together with predictions of transcription factor effects more generally, the analysis has led to model signaling networks largely accounting for arsenite and antimonite action. The striking parallels between responses to arsenite and antimonite indicate the skin carcinogenic risk of exposure to antimonite merits close scrutiny.

Transcriptomic Profiles of Monocyte-Derived Macrophages in Response to Escherichia coli is Associated with the Host Genetics.

Reactive Nitrogen Species (RNS) are a group of bactericidal molecules produced by macrophages in response to pathogens in a process called oxidative burst. Nitric oxide (NO-) is a member of RNS produced from arginine by inducible Nitric Oxide Synthase (iNOS) enzyme. The activity of iNOS and production of NO- by macrophages following stimulation is one of the indicators of macrophage polarization towards M1/proinflammatory. Production of NO- by bovine monocyte-derived macrophage (MDM) and mouse peritoneal macrophages has been shown to be strongly associated with host genetic with the heritability of 0.776 in bovine MDM and 0.8 in mouse peritoneal macrophages. However, the mechanism of genetic regulation of macrophage response has remained less explored. In the current study, the transcriptome of bovine MDMs was compared between two extreme phenotypes that had been classified as high and low responder based on NO- production. The results showed that 179 and 392 genes were differentially expressed (DE) between high and low responder groups at 3 and 18 hours after exposure to Escherichia coli, respectively. A set of 11 Transcription Factors (TFs) (STAT1, IRF7, SPI1, STAT4, IRF1, HIF1A, FOXO3, REL, NFAT5, HIC1, and IRF4) at 3 hours and a set of 13 TFs (STAT1, IRF1, HIF1A, STAT4, ATF4, TP63, EGR1, CDKN2A, RBL1, E2F1, PRDM1, GATA3, and IRF4) at 18 hours after exposure to E. coli were identified to be differentially regulated between the high and low responder phenotypes. These TFs were found to be divided into two clusters of inflammatory- and hypoxia-related TFs. Functional analysis revealed that some key canonical pathways such as phagocytosis, chemotaxis, antigen presentation, and cell-to-cell signalling are enriched among the over-expressed genes by high responder phenotype. Based on the results of this study, it was inferred that the functional characteristics of bovine MDMs are associated with NO-based classification. Since NO- production is strongly associated with host genetics, this study for the first time shows the distinct proinflammatory profiles of macrophages are controlled by the natural genetic polymorphism in an outbred population. In addition, the results suggest that genetics can be considered as a new dimension in the current model of macrophage polarization which is currently described by the combination of stimulants, only.

Non-Coding RNA Sequencing of Equine Endometrium During Maternal Recognition of Pregnancy.

Maternal recognition of pregnancy (MRP) in the mare is not well defined. In a non-pregnant mare, prostaglandin F2α (PGF) is released on day 14 post-ovulation (PO) to cause luteal regression, resulting in loss of progesterone production. Equine MRP occurs prior to day 14 to halt PGF production. Studies have failed to identify a gene candidate for MRP, so attention has turned to small, non-coding RNAs. The objective of this study was to evaluate small RNA (<200 nucleotides) content in endometrium during MRP. Mares were used in a cross-over design with each having a pregnant and non-mated cycle. Each mare was randomly assigned to collection day 11 or 13 PO (n = 3/day) and endometrial biopsies were obtained. Total RNA was isolated and sequencing libraries were prepared using a small RNA library preparation kit and sequenced on a HiSeq 2000. EquCab3 was used as the reference genome and DESeq2 was used for statistical analysis. On day 11, 419 ncRNAs, representing miRNA, snRNA, snoRNA, scaRNA, and vaultRNA, were different between pregnancy statuses, but none on day 13. Equine endometrial ncRNAs with unknown structure and function were also identified. This study is the first to describe ncRNA transcriptome in equine endometrium. Identifying targets of these ncRNAs could lead to determining MRP.

Coding RNA Sequencing of Equine Endometrium during Maternal Recognition of Pregnancy.

Equine maternal recognition of pregnancy (MRP) is a process whose signal remains unknown. During MRP the conceptus and endometrium communicate to attenuate prostaglandin F2α (PGF) secretion, sparing the corpus luteum and maintaining progesterone production. Recognition of a mobile conceptus by the endometrium is critical by days 14-16 post-ovulation (PO), when endometrium produces PGF, initiating luteolysis. The objective of this study was to evaluate endometrial gene expression changes based upon pregnancy status via RNA sequencing. This experiment utilized a cross-over design with each mare serving as both a pregnant and non-mated control on days nine, 11, and 13 PO (n = 3/status/day). Mares were randomly assigned to collection day and pregnancy confirmed by terminal uterine lavage at the time of endometrial biopsy. Total RNA was isolated and libraries prepared using Illumina TruSeq RNA sample preparation kit. Reads were mapped and annotated using HISAT2 and Stringtie. Expression values were evaluated with DESEQ2 (P ≤ 0.05 indicated significance). On day nine, 11, and 13 there were 1435, 1435 and 916 significant transcripts, respectively. Multiple genes with splice variants had different expression patterns within the same day. These are the first data to evaluate the endometrial transcriptome during MRP on days nine, 11, and 13.

Transcriptome profiles of isolated murine Achilles tendon proper- and peritenon-derived progenitor cells.

Progenitor cells of the tendon proper and peritenon have unique properties that could impact their utilization in tendon repair strategies. While a few markers have been found to aid in distinguishing progenitors cells from each region, there is great value in identifying more markers. In this study, we hypothesized that RNAseq could be used to improve our understanding of those markers that define these cell types. Transcriptome profiles were generated for pools of mouse Achilles tendon progenitor cells from both regions and catalogues of potential markers were generated. Moreover, common (e.g., glycoprotein, signaling, and proteinaceous extracellular matrix) and unique (e.g., cartilage development versus angiogenesis and muscle contraction) biological processes and molecular functions were described for progenitors from each region. Real-time quantitative PCR of a subset of genes was used to gain insight into the heterogeneity amongst individual progenitor colonies from each region. Markers like Scx, Mkx, Thbs4, and Wnt10a were consistently able to distinguish tendon proper progenitors from peritenon progenitors; expression variability for other genes suggested greater cell type complexity for potential peritenon progenitor markers. This is the first effort to define Achilles tendon progenitor markers by region. Further efforts to investigate the value of these cataloged markers are required by screening more individual colonies of progenitors for more markers. Clinical Significance: Findings from this study advance efforts in the discernment of cell type specific markers for tendon proper and peritenon progenitor cells; insight into marker sets could improve tracking and sorting strategies for these cells for future therapeutic strategies. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:1409-1418, 2019.

In vitro breeding: application of embryonic stem cells to animal production†.

Embryonic stem cells (ESCs) are derived from the inner cell mass of preimplantation blastocysts. For decades, attempts to efficiently derive ESCs in animal livestock species have been unsuccessful, but this goal has recently been achieved in cattle. Together with the recent reconstitution of the germ cell differentiation processes from ESCs in mice, these achievements open new avenues for the development of promising technologies oriented toward improving health, animal production, and the environment. In this article, we present a strategy that will notably accelerate genetic improvement in livestock populations by reducing the generational interval, namely in vitro breeding (IVB). IVB combines genomic selection, a widely used strategy for genetically improving livestock, with ESC derivation and in vitro differentiation of germ cells from pluripotent stem cells. We also review the most recent findings in the fields on which IVB is based. Evidence suggests this strategy will be soon within reach.

Combining multi-OMICs information to identify key-regulator genes for pleiotropic effect on fertility and production traits in beef cattle.

The identification of biological processes related to the regulation of complex traits is a difficult task. Commonly, complex traits are regulated through a multitude of genes contributing each to a small part of the total genetic variance. Additionally, some loci can simultaneously regulate several complex traits, a phenomenon defined as pleiotropy. The lack of understanding on the biological processes responsible for the regulation of these traits results in the decrease of selection efficiency and the selection of undesirable hitchhiking effects. The identification of pleiotropic key-regulator genes can assist in developing important tools for investigating biological processes underlying complex traits. A multi-breed and multi-OMICs approach was applied to study the pleiotropic effects of key-regulator genes using three independent beef cattle populations evaluated for fertility traits. A pleiotropic map for 32 traits related to growth, feed efficiency, carcass and meat quality, and reproduction was used to identify genes shared among the different populations and breeds in pleiotropic regions. Furthermore, data-mining analyses were performed using the Cattle QTL database (CattleQTLdb) to identify the QTL category annotated in the regions around the genes shared among breeds. This approach allowed the identification of a main gene network (composed of 38 genes) shared among breeds. This gene network was significantly associated with thyroid activity, among other biological processes, and displayed a high regulatory potential. In addition, it was possible to identify genes with pleiotropic effects related to crucial biological processes that regulate economically relevant traits associated with fertility, production and health, such as MYC, PPARG, GSK3B, TG and IYD genes. These genes will be further investigated to better understand the biological processes involved in the expression of complex traits and assist in the identification of functional variants associated with undesirable phenotypes, such as decreased fertility, poor feed efficiency and negative energetic balance.

Genome-wide identification of tissue-specific long non-coding RNA in three farm animal species.

BACKGROUND: Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs. RESULTS: Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen. CONCLUSIONS: LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Polymorphisms within the prolactin and growth hormone/insulin-like growth factor-1 functional pathways associated with fertility traits in Holstein cows raised in a hot-humid climate.

Prolactin (PRL), growth hormone (GH), and insulin-like growth factor-1 (IGF-1) are in hormone-response pathways involved in energy metabolism during thermoregulation processes in cattle. Objective herein was to study the association between single nucleotide polymorphisms (SNP) within genes of the PRL and GH/IGF-1 pathways with fertility traits such as services per conception (SPC) and days open (DO) in Holstein cattle lactating under a hot-humid climate. Ambient temperature and relative humidity were used to calculate the temperature-humidity index (THI) which revealed that the cows were exposed to heat stress conditions from June to November of 2012 in southern Sonora, Mexico. Individual blood samples from all cows were collected, spotted on FTA cards, and used to genotype a 179 tag SNP panel within 44 genes from the PRL and GH/IGF-1 pathways. The associative analyses among SNP genotypes and fertility traits were performed using mixed-effect models. Allele substitution effects were calculated using a regression model that included the genotype term as covariate. Single-SNP association analyses indicated that eight SNP within the genes IGF-1, IGF-1R, IGFBP5, PAPPA1, PMCH, PRLR, SOCS5, and SSTR2 were associated with SPC (P < 0.05), whereas four SNP in the genes GHR, PAPPA2, PRLR, and SOCS4 were associated with DO (P < 0.05). In conclusion, SNP within genes of the PRL and GH/IGF-1 pathways resulted as predictors of reproductive phenotypes in heat-stressed Holstein cows, and these SNP are proposed as candidates for a marker-assisted selection program intended to improve fertility of dairy cattle raised in warm climates.

STAT6, PBX2, and PBRM1 Emerge as Predicted Regulators of 452 Differentially Expressed Genes Associated With Puberty in Brahman Heifers.

The liver plays a central role in metabolism and produces important hormones. Hepatic estrogen receptors and the release of insulin-like growth factor 1 (IGF1) are critical links between liver function and the reproductive system. However, the role of liver in pubertal development is not fully understood. To explore this question, we applied transcriptomic analyses to liver samples of pre- and post-pubertal Brahman heifers and identified differentially expressed (DE) genes and genes encoding transcription factors (TFs). Differential expression of genes suggests potential biological mechanisms and pathways linking liver function to puberty. The analyses identified 452 DE genes and 82 TF with significant contribution to differential gene expression by using a regulatory impact factor metric. Brain-derived neurotrophic factor was observed as the most down-regulated gene (P = 0.003) in post-pubertal heifers and we propose this gene influences pubertal development in Brahman heifers. Additionally, co-expression network analysis provided evidence for three TF as key regulators of liver function during pubertal development: the signal transducer and activator of transcription 6, PBX homeobox 2, and polybromo 1. Pathway enrichment analysis identified transforming growth factor-beta and Wnt signaling pathways as significant annotation terms for the list of DE genes and TF in the co-expression network. Molecular information regarding genes and pathways described in this work are important to further our understanding of puberty onset in Brahman heifers.