Search

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.

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.

Summary: Progesterone signaling and uterine function are crucial in terms of pregnancy establishment. To investigate how the uterine tissue and its secretion changes in relation to puberty, we sampled tissue and uterine fluid from six pre‐ and six post‐pubertal Brahman heifers. Post‐pubertal heifers were sampled in the luteal phase. Gene expression of the uterine tissue was investigated with RNA‐sequencing, whereas the uterine fluid was used for protein profiling with mass spectrometry. A total of 4034 genes were differentially expressed (DE) at a nominal P‐value of 0.05, and 26 genes were significantly DE after Bonferroni correction (P < 3.1 × 10−6). We also identified 79 proteins (out of 230 proteins) that were DE (P < 1 × 10−5) in the uterine fluid. When we compared proteomics and transcriptome results, four DE proteins were identified as being encoded by DE genes: OVGP1, GRP, CAP1 and HBA. Except for CAP1, the other three had lower expression post‐puberty. The function of these four genes hypothetically related to preparation of the uterus for a potential pregnancy is discussed in the context of puberty. All DE genes and proteins were also used in pathway and ontology enrichment analyses to investigate overall function. The DE genes were enriched for terms related to ribosomal activity. Transcription factors that were deemed key regulators of DE genes are also reported. Transcription factors ZNF567,ZNF775,RELA,PIAS2,LHX4,SOX2,MEF2C,ZNF354C,HMG20A,TCF7L2,ZNF420,HIC1,GTF3A and two novel genes had the highest regulatory impact factor scores. These data can help to understand how puberty influences uterine function. [ABSTRACT FROM AUTHOR]

To understand genes, pathways, and networks related to puberty, we characterized the transcriptome of two tissues: the pituitary gland and ovaries. Samples were harvested from pre- and postpubertal Brahman heifers (same age group). Brahman heifers (Bos indicus) are older at puberty compared with Bos taurus, a productivity issue. With RNA sequencing, we identified differentially expressed (DEx) genes and important transcription factors (TF) and predicted coexpression networks. The number of DEx genes detected in the pituitary gland was 284 (P < 0.05), and VWC2L was the most DEx gene (fold change = 4.12, P = 0.01). The gene VWC2L promotes bone mineralization through transforming growth factor-β (TGFβ) signaling. Further studies of the link between bone mineralization and puberty could target VWC2L. In ovaries, 3,871 genes were DEx (P < 0.05). Four highly DEx genes were noteworthy for their function: SLC6A13 (a γ-aminobutyric acid [GABA] transporter), OXT (oxytocin), and NPY (neuropeptide Y) and its receptor NPY2R. These genes had higher ovarian expression in postpubertal heifers. The GABA and its receptors and transporters were expressed in the ovaries of many mammals, suggesting a role for this pathway beyond the brain. The OXT pathway has been known to influence the timing of puberty in rats, via modulation of GnRH. The effects of NPY at the hypothalamus, pituitary gland, and ovaries have been documented. Neuropeptide Y and its receptors are known factors in the release of GnRH, similar to OXT and GABA, although their roles in ovarian tissue are less clear. Pathways previously related to puberty such as TGFβ signaling (P = 6.71 × 10-5), Wnt signaling (P = 4.1 × 10-2), and peroxisome proliferator-activated receptor (PPAR) signaling (P = 4.84 × 10-2) were enriched in our data set. Seven genes were identified as key TF in both tissues: HIC2, ZIC4, ZNF219, ZSCAN26, LHX1, OLIG1, and a novel gene. An ovarian subnetwork created with TF and significant ovarian DEx genes revealed five zinc fingers as regulators: ZNF507, ZNF12, ZNF512, ZNF184, and ZNF432. Recent work of hypothalamic gene expression also pointed to zinc fingers as TF for bovine puberty. Although some zinc fingers may be ubiquitously expressed, the identification of DEx genes in common across tissues points to key regulators of puberty. The hypothalamus and pituitary gland had eight DEx genes in common. The hypothalamus and ovaries had 89 DEx genes in common. The pituitary gland and ovaries had 48 DEx genes in common. Our study confirmed the complexity of puberty and suggested further investigation on genes that code zinc fingers. [ABSTRACT FROM AUTHOR]

Puberty onset is a developmental process influenced by genetic determinants, environment, and nutrition. Mutations and regulatory gene networks constitute the molecular basis for the genetic determinants of puberty onset. The emerging knowledge of these genetic determinants presents opportunities for innovation in the breeding of early pubertal cattle. This paper presents new data on hypothalamic gene expression related to puberty in Bos indicus (Brahman) in age- and weight-matched heifers. Six postpubertal heifers were compared with 6 prepubertal heifers using whole-genome RNA sequencing methodology for quantification of global gene expression in the hypothalamus. Five transcription factors (TF) with potential regulatory roles in the hypothalamus were identified in this experiment: E2F8, NFAT5, SIX5, ZBTB38, and ZNF605. These TF genes were significantly differentially expressed in the hypothalamus of postpubertal versus prepubertal heifers and were also identified as significant according to the applied regulatory impact factor metric (P < 0.05). Two of these 5 TF, ZBTB38 and ZNF605, were zinc fingers, belonging to a gene family previously reported to have a central regulatory role in mammalian puberty. The SIX5 gene belongs to the family of homologues of Drosophila sine oculis (SIX) genes implicated in transcriptional regulation of gonadotrope gene expression. Tumor-related genes such as E2F8 and NFAT5 are known to affect basic cellular processes that are relevant in both cancer and developmental processes. Mutations in NFAT5 were associated with puberty in humans. Mutations in these TF, together with other genetic determinants previously discovered, could be used in genomic selection to predict the genetic merit of cattle (i.e., the likelihood of the offspring presenting earlier than average puberty for Brahman). Knowledge of key mutations involved in genetic traits is an advantage for genomic prediction because it can increase its accuracy. [ABSTRACT FROM AUTHOR]

The objective of this study was to examine five different sources of RNA, namely mammary gland tissue (MGT), milk somatic cells (SC), laser microdissected mammary epithelial cells (LCMEC), milk fat globules (MFG) and antibody-captured milk mammary epithelial cells (mMEC) to analyze the bovine mammary gland transcriptome using RNA-Sequencing. Our results provide a comparison between different sampling methods (invasive and non-invasive) to define the transcriptome of mammary gland tissue and milk cells. This information will be of value to investigators in choosing the most appropriate sampling method for different research applications to study specific physiological states during lactation. One of the simplest procedures to study the transcriptome associated with milk appears to be the isolation of total RNA directly from SC or MFG released into milk during lactation. Our results indicate that the SC and MFG transcriptome are representative of MGT and LCMEC and can be used as effective and alternative samples to study mammary gland expression without the need to perform a tissue biopsy. [ABSTRACT FROM AUTHOR]

The technological properties of milk have significant importance for the dairy industry. Citrate, a normal constituent of milk, forms one of the main buffer systems that regulate the equilibrium between Ca2+ and H+ ions. Higher-than-normal citrate content is associated with poor coagulation properties of milk. To identify the genes responsible for the variation of citrate content in milk in dairy cattle, the metabolic steps involved in citrate and fatty acid synthesis pathways in ruminant mammary tissue using RNA sequencing were studied. Genetic markers that could influence milk citrate content in Holstein cows were used in a marker-trait association study to establish the relationship between 74 single nucleotide polymorphisms (SNP) in 20 candidate genes and citrate content in 250 Holstein cows. This analysis revealed 6 SNP in key metabolic pathway genes [isocitrate dehydrogenase 1 (NADP+), soluble (IDH1); pyruvate dehydrogenase (lipoamide) β (PDHB); pyruvate kinase (PKM2); and solute carrier family 25 (mitochondrial carrier; citrate transporter), member 1 (SLC25A1)] significantly associated with increased milk citrate content. The amount of the phenotypic variation explained by the 6 SNP ranged from 10.1 to 13.7%. Also, genotype-combination analysis revealed the highest phenotypic variation was explained combining IDH1_23211, PDHB_5562, and SLC25A1_4446 genotypes. This specific genotype combination explained 21.3% of the phenotypic variation. The largest citrate associated effect was in the 3' untranslated region of the SLC25A1 gene, which is responsible for the transport of citrate across the mitochondrial inner membrane. This study provides an approach using RNA sequencing, metabolic pathway analysis, and association studies to identify genetic variation in functional target genes determining complex trait phenotypes. [ABSTRACT FROM AUTHOR]