RNA sequencing differential gene expression analysis of isolated perfused bovine udders experimentally inoculated with Streptococcus agalactiae.

The aim of this study was to elucidate the differential gene expression in the RNA sequencing transcriptome of isolated perfused udders collected from 4 slaughtered Holstein × Zebu crossbred dairy cows experimentally inoculated with Streptococcus agalactiae. We studied 3 different statistical tools (edgeR, baySeq, and Cuffdiff 2). In summary, 2 quarters of each udder were experimentally inoculated with Strep. agalactiae and the other 2 were used as a control. Mammary tissue biopsies were collected at times 0 and 3 h after infection. The total RNA was extracted and sequenced on an Illumina HiSeq 2000 (Illumina Inc., San Diego, CA). Transcripts were assembled from the reads aligned to the bovine UMD 3.1 reference genome, and the statistical analyses were performed using the previously mentioned tools (edgeR, baySeq, and Cuffdiff 2). Finally, the identified genes were submitted to pathway enrichment analysis. A total of 1,756, 1,161, and 3,389 genes with differential gene expression were identified when using edgeR, baySeq, and Cuffdiff 2, respectively. A total of 122 genes were identified by the overlapping of the 3 methods; however, only the platelet activation presented a significantly enriched pathway. From the results, we suggest the FCER1G, GNAI2, ORAI1, and VASP genes shared among the 3 methods in this pathway for posterior biological validation.

Pre- and post-puberty expression of genes and proteins in the uterus of Bos indicus heifers: the luteal phase effect post-puberty.

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.

Performance strategies affect mammary gland development in prepubertal heifers.

In Brazil, the majority of dairy cattle are Holstein × Gyr (H×G). It is unknown whether excessive energy intake negatively affects their mammary development to the same extent as in purebred Holsteins. We hypothesized that mammary development of H×G heifers can be affected by dietary energy supply. We evaluated the effect of different average daily gains (ADG) achieved by feeding different amounts of a standard diet during the growing period on biometric measurements, development of mammary parenchyma (PAR) and mammary fat pad (MFP), and blood hormones. At the outset of this 84-d experiment, H×G heifers (n = 18) weighed 102.2 ± 3.4 kg and were 3 to 4 mo of age. Heifers were randomly assigned to 1 of 3 ADG programs using a completely randomized design. Treatments were high gain (HG; n = 6), where heifers were fed to gain 1 kg/d; low gain (LG; n = 6), where heifers were fed to gain 0.5 kg/d; and maintenance (MA; n = 6), where heifers were fed to gain a minimal amount of weight per day. Heifers were fed varying amounts of a single TMR to support desired BW gains. Over the 84 d, periodic biometric and blood hormone measurements were obtained. On d 84, all heifers were slaughtered and carcass and mammary samples were collected. At the end, HG heifers weighed the most (181 ± 7.5 kg), followed by LG (146 ± 7.5 kg) and MA (107 ± 7.5 kg) heifers. The ADG were near expected values and averaged 0.907, 0.500, and 0.105 ± 0.03 kg/d for HG, LG, and MA, respectively. In addition, body lengths, heart girths, and withers heights were affected by dietary treatment, with MA heifers generally being the smallest and HG heifers generally being the largest. Body condition scores differed by treatment and were highest in HG and lowest in MA heifers; in vivo subcutaneous fat thickness measurement and direct analysis of carcass composition supported this. The HG heifers had the heaviest MFP, followed by LG and then MA heifers. Amount of PAR was highest in LG heifers and was the same for HG and MA heifers. The percentage of udder mass occupied by PAR was lowest in HG heifers, differing from LG and MA heifers. Composition of MFP was not evaluated. Regarding PAR composition, no differences in ash or DM were found. On the other hand, CP concentration of PAR for HG heifers was lower than that for LG heifers, which was lower than that for MA heifers. Regarding the fat content, HG treatment was higher than LG and MA treatment, which did not differ from each other. In PAR, differences in relative abundance of genes related to both stimulation and inhibition of mammary growth were observed to depend on dietary treatment, sampling day, or both. The same can be said for most of the blood hormones that were measured in this experiment. In this experiment, high ADG achieved by feeding different amounts of a standard diet during the growing period negatively affected mammary development.

Transcriptome analyses identify five transcription factors differentially expressed in the hypothalamus of post- versus prepubertal Brahman heifers.

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 (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: , , , , and . 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 ( < 0.05). Two of these 5 TF, and , were zinc fingers, belonging to a gene family previously reported to have a central regulatory role in mammalian puberty. The gene belongs to the family of homologues of Drosophila sine oculis () genes implicated in transcriptional regulation of gonadotrope gene expression. Tumor-related genes such as and are known to affect basic cellular processes that are relevant in both cancer and developmental processes. Mutations in 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.

Effects of nutrient intake level on mammary parenchyma growth and gene expression in crossbred (Holstein × Gyr) prepubertal heifers.

This study investigated the effects of increased nutrient intake levels on prepubertal mammary parenchyma development in crossbreed (Holstein × Gyr) dairy heifers. Eighteen heifers age 3 to 4 mo were fed 1 of 3 nutrient intake levels (n=6 per treatment) designed to sustain an average daily gain of 0.0kg/d (maintenance, MA), 0.5kg/d (low gain, LG), or 1.0kg/d (high gain, HG). Serum blood samples collected on d 42 and 84 after a 12-h fast were analyzed for triglycerides, leptin, insulin, and insulin-like growth factor 1 (IGF-1). Liver and mammary parenchyma were biopsied on d 42 and harvested on d 84 for gene expression analysis. Parenchyma samples were also used for biochemical and histological analysis. Mammary parenchyma weight was lower in HG than in MA or LG heifers, but mammary extraparenchymal fat was greater in HG heifers than in other groups. Heifers fed the HG diet had a greater fraction of ether extract in their parenchyma than the others and a smaller fraction of crude protein in their parenchyma than MA heifers. Moreover, the HG and LG heifers had greater body fat mass than MA heifers. Nutrient intake level had no effect on the number of intraparenchymal adipocytes. Heifers fed the HG diet had greater serum IGF-1 than the others, and serum insulin was lower in the MA than the HG or LG heifers. Liver GHR, IGF1, and IGFBP3 mRNA expression was higher, but IGFBP2 mRNA was lower in HG heifers than in others. The parenchyma mRNA expression of lipogenic markers, such as CD36, ACCA, FASN, and ADIPOR1, was upregulated by nutrient intake level. Significant nutrient intake × time interactions for lipogenic genes during the experimental period indicated variable gene expression depending on the time point of prepubertal mammary gland development. Overall, our data suggest that enhancing nutrient intake increased body fat accumulation and lipogenesis in the mammary gland to the detriment of parenchyma growth. Moreover, increased lipogenesis in the parenchyma of HG heifers may indicate that fat accumulation occurred because of adipocyte hypertrophy and not differences in adipogenesis. The implications of these results for milk yield needs to be elucidated.

Effect of maternal nutrition and days of gestation on pituitary gland and gonadal gene expression in cattle.

This study investigated effects of maternal overnutrition on gonadal development and pituitary-gonadal gene expression in cattle fetuses at mid- and late-gestation. Twenty-seven multiparous dry cows were fed either high (ad libitum, H) or moderate (M) intake of the same diet. Twelve cows from H (n=6) and M (n=6) intake carrying females fetuses were euthanized at 199 and 268d of gestation (DG; n=3 for H or M on each DG). Fifteen cows from H (n=6) and M intake (n=9) carrying male fetuses were euthanized at 139, 199, and 241 DG (n=2 for H and n=3 for M on each DG). Fetal gonads and pituitary gland were sampled for gene expression and histological analyses. Sex-specific responses to maternal intake were observed. Primordial and total follicle numbers were lower in fetal ovaries from H than in M intake cows. These results were the reverse for preantral and antral follicles. Volumetric proportion and diameter of seminiferous cord were lower in fetal testis of H than M intake cows. The expression level of FSHB was greater in pituitary gland of the female fetus from H compared with M intake cows, irrespective of DG, whereas LHB gene expression did not differ. In males, FSHB and LHB gene expression levels were similar between maternal intake groups. Fetal ovarian expression of P450 aromatase, StAR, BMPR2, TGFBR1, GDF9, FSHR, Bax, and CASP3 genes were higher in H than in M intake cows, irrespective of DG. Fetal testicular expression of StAR, HSD17B3, IGF1, IGF2, and IGF1R genes was higher in M than in H intake cows. The differences in gene expression for steroidogenesis, folliculogenesis, and apoptosis may explain the distinct pattern of follicular growth between offspring of M and H intake cows. By contrast, the lower volumetric proportion, diameter, and length of seminiferous cord may relate to decreased gene expression in fetal testis from H intake cows. In conclusion, maternal H intake seems to affect fetal ovarian follicular growth and number of follicles, which may affect the size of ovarian reserve in their offspring. In male fetus, maternal H intake seems to disturb testicular development and may have implications on sperm production. The underlying mechanism of differential gene expression and the effect on offspring reproductive potential should be the focus of further research, especially considering larger sample size, reducing the chance for type I errors.

Effect of heat stress and feeding phosphorus levels on pig electron transport chain gene expression.

The purpose of this study was to evaluate the effect of temperature and different levels of available phosphorus (aP) on the expression of nine genes encoding electron transport chain proteins in the Longissimus dorsi (LD) muscle of pigs. Two trials were carried out using 48 high-lean growth pigs from two different growth phases: from 15 to 30 kg (phase 1) and from 30 to 60 kg (phase 2). Pigs from growth phase 1 were fed with three different levels of dietary aP (0.107%, 0.321% or 0.535%) and submitted either to a thermoneutral (24°C and RH at 76%) or to a heat stress (34°C and RH at 70%) environment. Pigs from growth phase 2 were fed with three different levels of dietary aP (0.116%, 0.306% or 0.496%) and submitted either to a thermoneutral (22ºC and RH at 77%) or to a heat stress (32ºC and RH at 73%) environment. Heat stress decreased (P<0.001) average daily feed intake at both growth phases. At 24°C, pigs in phase 1 fed the 0.321% aP diet had greater average daily gain and feed conversion (P<0.05) than those fed the 0.107% or 0.535% while, at 34°C pigs fed the 0.535% aP had the best performance (P<0.05). Pigs from phase 2 fed the 0.306% aP had best performance in both thermal environments. Gene expression profile was analyzed by quantitative real-time polymerase chain reaction. Irrespective of growing phase, the expression of six genes was lower (P<0.05) at high temperature than at thermoneutrality. The lower expression of these genes under high temperatures evidences the effects of heat stress by decreasing oxidative metabolism, through adaptive physiological mechanisms in order to reduce heat production. In pigs from phase 1, six genes were differentially expressed across aP levels (P<0.05) in the thermoneutral and one gene in the heat stress. In pigs from phase 2, two genes were differentially expressed across aP levels (P<0.05) in both thermal environments. These data revealed strong evidence that phosphorus and thermal environments are key factors to regulate oxidative phosphorylation with direct implications on animal performance.