Canonical pathways and networks regulated by estrogen in the bovine mammary gland

Carry-over effects of maternal late-gestation heat stress on granddaughter's growth and mammary gland development

Maternal (F0) exposure to late-gestation heat stress reduces their daughter's (F1) mammary gland fat pad mass (FP), parenchyma (PAR) mass, and epithelial cell proliferation when evaluated at birth and weaning, and go on to produce less milk in their first lactation. Herein, we investigated the effect of maternal late-gestation heat stress on whole-body growth and mammary development of their granddaughters (F2). Multiparous F0 cows had access to heat abatement (n = 41, shade, and active cooling via fans and water soakers) or not (n = 41, shade only) for the last 56 d of gestation during a subtropical summer. Consequently, the F1 daughters, born to F0 cows, were heat-stressed (HTF1, n = 36) or cooled (CLF1, n = 37) in utero during the last 2 mo of gestation. All F1 heifers were raised as an identically managed cohort until first calving. The F2 granddaughters, born to HTF1 (HTF2, n = 12) or CLF1 (CLF2, n = 17), were raised as an identically managed cohort until 70 d of age. Dry matter intake (DMI), body weight, hip height, wither height, chest girth, head circumference, mammary gland teat length, and left-right and front-rear teat distances were measured. Average daily gain (ADG) was calculated for the pre-weaned period (0-49 d). Mammary ultrasounds were performed on d 21, 49, and 70 (n = 9/group) on the rear left and right quarters to quantify PAR and FP areas. Mammary biopsies were collected for histological evaluation of epithelial structures (H&E staining), and to quantify cells positive for ERα (estrogen receptor, α subunit), cell proliferation (Ki67), and apoptosis (TUNEL). Heifer growth from birth to d 49 was similar between CLF2 and HTF2 for all parameters evaluated. Distances between teats and teat length were not different between groups. On d 70, CLF2 tended to have a greater average PAR (right and left quarters) relative to HTF2. Although the left FP was smaller in HTF2 relative to CLF2, the average FP was not different. The lumenal and non-lumenal epithelial structures in the PAR of HTF2 were significantly smaller than those of CLF2. In addition, HTF2 had a reduced percentage of proliferating cells in the epithelial and stromal compartments and a greater percentage of apoptotic cells, particularly in the stroma. The percentage of ERα positive cells was significantly reduced in HTF2. In summary, although HTF2 heifer's DMI was similar and they grew at the same rate as CLF2 heifers throughout the pre-weaning phase, their mammary glands had smaller PAR areas with fewer epithelial structures characterized by reduced cell turnover and lower ERα expression. These early changes in the microstructure and cellular turnover of the mammary gland may partly explain the reduction in lactation performance relative to CLF2 counterparts at maturity.

Progress on the Regulation of Ruminant Milk Fat by Noncoding RNAs and ceRNAs

Milk fat is not only a key factor affecting the quality of fresh milk but also a major target trait forbreeding. The regulation of milk fat involves multiple genes, network regulation and signal transduction. To explore recent discoveries of pathway regulation, we reviewed the published literature with a focus on functional noncoding RNAs and epigenetic regulation in ruminants. Results indicate that miRNAs play key roles in the regulation of milk fat synthesis and catabolism in ruminants. Although few data are available, merging evidence indicates that lncRNAs and circRNAs act on milk fat related genes through indirect action with microRNAs or RNAs in the ceRNA network to elicit positive effects on transcription. Although precise regulatory mechanisms remain unclear, most studies have focused on the regulation of the function of target genes through functional noncoding RNAs. Data to help identify factors that can regulate their own expression and function or to determine whether self-regulation involves positive and/or negative feedback are needed. Despite the growing body of research on the role of functional noncoding RNA in the control of ruminant milk fat, most data are still not translatable for field applications. Overall, the understanding of mechanisms whereby miRNA, lncRNA, circRNA, and ceRNA regulate ruminant milk fat remains an exciting area of research.

TRIENNIAL LACTATION SYMPOSIUM/BOLFA: Plasticity of mammary development in the prepubertal bovine mammary gland

Although peripubertal mammary development represents only a small fraction of the total mass of mammary parenchyma present in the udder at the end of gestation and into lactation, there is increasing evidence that the tissue foundations created in early life can affect future mammary development and function. Studies on expression of estrogen and progesterone receptors seem to confirm the relevance of these steroids in prepubertal mammary development, but connections with other growth factors, hormones, and local tissue factors remain elusive. Enhanced preweaning feeding in the bovine appears to enhance the capacity of mammary tissue to response to mammogenic stimulation. This suggests the possibility that improved early nutrition might allow for creation of stem or progenitor cell populations to better support the massive ductal growth and lobulo-alveolar development during gestation. Increasing evidence that immune cells are involved in mammary development suggests there are unexpected and poorly understood connections between the immune system and mammary development. This is nearly unexplored in ruminants. Development of new tools to identify, isolate, and characterize cell populations within the developing bovine mammary gland offer the possibility of identifying and perhaps altering populations of mammary stem cells or selected progenitor cells to modulate mammary development and, possibly, mammary function.

Determination of genetic effects of ATF3 and CDKN1A genes on milk yield and compositions in Chinese Holstein population

Background

Our previous RNA-sequencing study revealed that the ATF3 and CDKN1A genes were remarkably differentially expressed between the mammary glands of lactating Holstein cows with extremely high and low milk protein and fat percentage so that both of them were considered as candidates for milk composition. Herein, we further verified whether these genes have genetic effects on milk production traits in a Chinese Holstein cow population.

Results

By re-sequencing the entire coding and regulatory regions, we identified four SNPs in 5’promoter region, two in exons, seven in 3′ un-translated region (UTR), and six in 3’flanking region of ATF3 gene, and one SNP in exon 5, two in 3’UTR, and two in 3’flanking region of CDKN1A gene. Of these, only the SNP, c.271C > T (rs442346530), in exon 5 of CDKN1A gene was predicted to result in an amino acid replacement from arginine to tryptophan. Subsequent genotype-phenotype association analysis revealed that 19 SNPs in ATF3 and 5 SNPs in CDKN1A were evidently associated with 305-days milk yield, fat yield, protein yield, or protein percentage (P = < 0.0001 ~ 0.0494). Whilst, no significant SNPs in ATF3 gene were associated with fat percentage in both first and second lactations (P > 0.05), and only two SNPs of CDKN1A gene, c.271C > T (P = 0.0377) and c.*654C > T (P = 0.0144), were markedly associated with fat percentage in the first lactation. Further, linkage disequilibrium (LD) analyses were conducted among the identified SNPs in ATF3 and/or CDKN1A genes to further confirm the association results. We also observed that the four SNPs, g.72834301C > A, g.72834229C > A, g.72833969A > G, and g.72833562G > T altered the specific transcription factor (TF) binding sites in ATF3 promoter, and one SNP, c.271C > T, changed the CDKN1A protein secondary structure, suggesting they might be the promising potential functional mutations.

Conclusion

Our findings first profiled the genetic effects of ATF3 and CDKN1A genes for milk production traits in dairy cattle and will be available for marker-assisted breeding in dairy cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0516-4) contains supplementary material, which is available to authorized users.

Mammary transcriptome analysis of lactating dairy cows following administration of bovine growth hormone

The galactopoietic effect of growth hormone (GH) in lactating ruminants is well established; however the mechanisms that mediate these effects are not well understood. The first objective of this study was to determine the effect of GH on the synthesis of the major casein and whey proteins. The second objective was to identify the genes and pathways that may be involved in mediating the effect of GH on milk synthesis. A single subcutaneous injection of a commercially available slow release formulation of GH (Lactatropin®), or physiological saline solution (control) was administered to non-pregnant dairy cows (n=4/group) in mid-late lactation. Milk samples were collected for composition analysis and mammary lobulo-alveolar tissue was collected postmortem 6 days post injection. Gene expression profiles were evaluated using either a 22 000 bovine complementary DNA microarray or quantitative PCR (qPCR), and microarrays were validated by qPCR. The yield of all the major casein and whey proteins was increased 32% to 41% in GH-treated cows, with the exception of α-lactalbumin yield which was elevated by 70% relative to controls. Treatment with GH treatment tended to increase the concentration of α-lactalbumin but had no effect on the concentration of any of the major milk proteins. Messenger RNA (mRNA) abundance of the major whey and casein genes, with the exception of α-s2-casein, was increased in response to GH compared with controls, which is consistent with the positive effect of GH on milk production. Treatment with GH treatment influenced the mRNA abundance of genes involved in cell growth and proliferation, transcriptional and translational regulation, actin cytoskeleton signalling, lipid metabolism and cell death. This study has provided new insights into the cell signalling that may be involved in mediating the effect of GH on milk production in the mammary gland of lactating dairy cows.

Expression of CD24 in Human Bone Marrow-Derived Mesenchymal Stromal Cells Is Regulated by TGFβ3 and Induces a Myofibroblast-Like Genotype

Human bone marrow-derived stromal cells (hBMSCs) derived from the adult organism hold great promise for diverse settings in regenerative medicine. Therefore a more complete understanding of hBMSC biology to fully exploit the cells' potential for clinical settings is important. The protein CD24 has been reported to be involved in a diverse range of processes such as cancer, adaptive immunity, inflammation, and autoimmune diseases in other cell types. Its expression in hBMSCs, which has not yet been analyzed, may add an important aspect in the understanding of hBMSC biology. The present study therefore analyzes the expression, regulation, and functional implication of the surface protein CD24 in hBMSCs. Methods used are stimulation studies with TGF beta as well as shRNA-mediated knockdown and overexpression of CD24 followed by microarray, immunocytochemistry, and flow cytometric analyses. To our knowledge, we demonstrate for the first time that the expression of CD24 is an inherent property of hBMSCs. Importantly, the data links the upregulation of CD24 to the adoption of a myofibroblast-like gene expression pattern in hBMSCs. We demonstrate that CD24 is an important modulator in transforming growth factor beta 3 (TGFβ3) signaling with a reciprocal regulatory relationship between these two proteins.

MicroRNA-regulated molecular mechanism underlying bovine subclinical endometritis

An impaired uterine environment triggered by the incidence of subclinical endometritis often compromises fertility in the bovine. The uterus is a dynamic organ with tight regulation of specific genes at the transcriptional and translational levels. Herein, we hypothesised that subclinical endometritis alters the expression of uterine microRNAs (miRNAs), which may result in the dysregulation of corresponding target genes and biological pathways. To test this hypothesis, we used a genome-wide RT(2) (Exiqon, Vedbaek, Denmark) miRNA PCR array consisting of 354 miRNA primers and analysed miRNA expression in uterine cytobrush samples taken from cows with and without subclinical endometritis. The results revealed aberrant expression of 23 miRNAs in cows with subclinical endometritis compared with healthy cows. Furthermore, we designed an in vitro endometrial cell culture model challenged by lipopolysaccharide (LPS) to validate the differential regulation of miRNAs in cytobrush samples. Interestingly, we observed similar expression miRNA patterns in cytobrush samples taken from cows with or without subclinical endometritis and in vitro cultured endometrial cells challenged by LPS. To trace signalling pathways and biological functions potentially controlled by the aberrantly expressed miRNAs, we filtered high-ranking target genes from miRBase and analysed them using ingenuity pathway analysis. The gene networks, canonical pathways and biological functions strikingly converged to signalling pathways that mediate inflammatory responses, cellular proliferation, cell movement, the cell cycle and apoptosis in the bovine endometrium. In addition, expression analysis of key genes from the gene networks confirmed their presence and the potential regulation of these genes by uterine miRNAs. Furthermore, luciferase assay data substantiated the primary information from bioinformatic prediction that generated potential target genes for the dysregulated miRNAs in subclinical endometritis. Together, these data suggest the potential regulatory role of uterine miRNAs in the development and progression of bovine subclinical endometritis.

Responses of the mammary transcriptome of dairy cows to altered photoperiod during late gestation

Cows exposed to short day photoperiod (SD, 8L:16D) during the 60-day nonlactating period prior to parturition produce more milk in their subsequent lactation compared with cows exposed to long day photoperiod (LD, 16L:8D). Although this response is well established in dairy cows, the underlying mechanisms are not understood. We hypothesized that differential gene expression in cows exposed to SD or LD photoperiods during the dry period could be used to identify the functional basis for the subsequent increase in milk production during lactation. Pregnant, multiparous cows were maintained on an SD or LD photoperiod for 60 days prior to parturition. Mammary biopsies were obtained on days -24 and -9 relative to parturition and Affymetrix GeneChip Bovine Genome Arrays were used to quantify gene expression. Sixty-four genes were differentially expressed (P ≤ 0.05 and fold-change ≥ |1.5|) between SD and LD treatments. Many of these genes were associated with cell growth and proliferation, or immune function. Ingenuity Pathway Analysis predicted upstream regulators to include TNF, TGF-β1, interferon-γ, and several interleukins. In addition, expression of 125 genes was significantly different between day -24 and day -9; those genes were associated with milk component metabolism and immune function. The interaction of photoperiod and time affected 32 genes associated with insulin-like growth factor I signaling. Genes differentially expressed in response to photoperiod were associated with mammary development and immune function consistent with the enhancement of milk yield in the ensuing lactation. Our results provide insight into the mechanisms by which photoperiod affects the mammary gland and subsequently lactation.

IGF-I retards proper development of acinar structures formed by bovine mammary epithelial cells via sustained activation of Akt kinase

Insulin-like growth factor-I is involved in mammary gland development, promoting proliferation and inhibiting apoptosis of mammary epithelial cells (MECs). Mitogenic actions of IGF-I are mainly mediated by the phosphatidylinositol-3 kinase (PI3K)/Akt signaling pathway. We have found that in the presence of IGF-I bovine BME-UV1 MECs cultured on reconstituted basement membrane form large spheroids with disrupted polarity and no cavity in the center. These cells showed enhanced phosphorylation of Akt, decreased level of cleaved caspase-3, and sustained proliferative activity throughout the 16-d period of 3-dimensional culture. Inhibition of the PI3K/Akt pathway by a specific inhibitor of PI3K, LY294002, resulted in the restoration of the normal acinar phenotype. However, this effect was noted only when LY294002 was added in the second week of 3-dimensional culture, which corresponded with the time of cell cycle arrest and polarity formation under control conditions. Normal development of acini was also obtained when BME-UV1 cells were treated simultaneously with IGF-I and 17β-estradiol. The addition of 17β-estradiol regulated Akt activation, enabling the subsequent initiation of polarization processes. 17β-Estradiol also increased the level of IGFBP-3 protein in MECs cultured on Matrigel in the presence of IGF-I. The presented results indicate important interactions between signaling pathways activated by estrogen and IGF-I, which regulate alveologenesis in bovine mammary gland.

Estradiol and progesterone exhibit similar patterns of hepatic gene expression regulation in the bovine model

Female sex steroid hormones, estradiol-17β (E2-17β) and progesterone (P4) regulate reproductive function and gene expression in a broad range of tissues. Given the central role of the liver in regulating homeostasis including steroid hormone metabolism, we sought to understand how E2-17β and P4 interact to affect global gene expression in liver. Ovariectomized cows (n = 8) were randomly assigned to 4 treatment groups applied in a replicated Latin Square design: 1) No hormone supplementation, 2) E2-17β treatment (ear implant), 3) P4 treatment (intravaginal inserts), and 4) E2-17β combined with P4. After 14 d of treatment, liver biopsies were collected, allowing 28 d intervals between periods. Changes in gene expression in the liver biopsies were monitored using bovine-specific arrays. Treatment with E2-17β altered expression of 479 genes, P4 472 genes, and combined treatment significantly altered expression of 468 genes. In total, 578 genes exhibited altered expression including a remarkable number (346 genes) that responded similarly to E2-17β, P4, or combined treatment. Additional evidence for similar gene expression actions of E2-17ß and/or P4 were: principal component analysis placed almost every treatment array at a substantial distance from controls; Venn diagrams indicated overall treatment effects for most regulated genes; clustering analysis indicated the two major clusters had all treatments up-regulating (172 genes) or down-regulating (173 genes) expression. Thus, unexpectedly, common biological pathways were regulated by E2-17β and/or P4 in liver. This indicates that the mechanism of action of these steroid hormones in the liver might be either indirect or might occur through non-genomic pathways. This unusual pattern of gene expression in response to steroid hormones is consistent with the idea that there are classical and non-classical tissue-specific responses to steroid hormone actions. Future studies are needed to elucidate putative mechanism(s) responsible for overlapping actions of E2-17β and P4 on the liver transcriptome.

[Genome-wide selective sweep analysis on Large White and Tongcheng pigs]

The production performance of pigs has been significantly improved due to long-term artificial selection, and the specific variation characterizations (selection signatures) emerged from the selected genome regions. Different types of breeds are subjected to different selection intensities and had different selection signatures. Selective sweep analysis is one of major methods to detect the selection signatures. In this study, based on the 60K BeadChip genotyping data of both commercial Large White (n=45) and local Tongcheng pigs (n=45), genetic differentiation coefficient Fst was applied to detect the selection signatures. Using gPLINK software to set quality control standards, a total of 34 304 SNPs were selected for statistical analysis. Fst values between two breeds were estimated with Genepop package and the average Fst value was 0.3209. Setting Fst>0.7036 (1% of total number of Fst values) as selection threshold, 344 SNPs were obtained and SNP location annotation indicated that there were 79 candidate genes (Sus scrofa Build 9). Furthermore, network analysis was performed using Ingenuity Pathway Analysis and the preliminary results suggested that most genes were involved in growth, reproduction, and immune response, such as NCOA6, ERBB4, RUNX2, and APOB genes. The findings from this study will contribute to further identification of candidate genes and causal mutations implying for meat production and disease resistance in pig.

Granule exocytosis of granulysin and granzyme B as a potential key mechanism in vaccine-induced immunity in cattle against the nematode Ostertagia ostertagi

Ostertagia ostertagi is considered one of the most economically important bovine parasites. As an alternative to anthelmintic treatment, an experimental host-protective vaccine was previously developed on the basis of ASP proteins derived from adult worms. Intramuscular injection of this vaccine, combined with QuilA as an adjuvant, significantly reduced fecal egg counts by 59%. However, the immunological mechanisms triggered by the vaccine are still unclear. Therefore, in this study, the differences in immune responses at the site of infection, i.e., the abomasal mucosa, between ASP-QuilA-vaccinated animals and QuilA-vaccinated control animals were investigated on a transcriptomic level by using a whole-genome bovine microarray combined with histological analysis. Sixty-nine genes were significantly impacted in animals protected by the vaccine, 48 of which were upregulated. A correlation study between the parasitological parameters and gene transcription levels showed that the transcription levels of two of the upregulated genes, those for granulysin (GNLY) and granzyme B (GZMB), were negatively correlated with cumulative fecal egg counts and total worm counts, respectively. Both genes were also positively correlated with each other and with another upregulated gene, that for the IgE receptor subunit (FCER1A). Surprisingly, these three genes were also correlated significantly with CMA1, which encodes a mast cell marker, and with counts of mast cells and cells previously described as globule leukocytes. Furthermore, immunohistochemical data showed that GNLY was present in the granules of globule leukocytes and that it was secreted in mucus. Overall, the results suggest a potential role for granule exocytosis by globule leukocytes, potentially IgE mediated, in vaccine-induced protection against O. ostertagi.

β-Estradiol-dependent activation of the JAK/STAT pathway requires p/CIP and CARM1

The steroid receptor coactivator p/CIP, also known as SRC-3, is an oncogene commonly amplified in breast and ovarian cancers. p/CIP is known to associate with coactivator arginine methyltransferase 1 (CARM1) on select estrogen responsive genes. We have shown, using a ChIP-on-chip approach, that in response to stimulation with 17β-estradiol (E2), the p/CIP/CARM1 complex is recruited to 204 proximal promoters in MCF-7 cells. Many of the complex target genes have been previously implicated in signaling pathways related to oncogenesis. Jak2, a member of the Jak/Stat signaling cascade, is one of the direct E2-dependent targets of the p/CIP/CARM1 complex. Following E2-treatment, histone modifications at the Jak2 promoter are reflective of a transcriptionally permissive gene, and modest changes in RNA and protein expression lead us to suggest that an additional factor(s) may be required for a more notable transcriptional and functional response. Bioinformatic examination of the 204 proximal promoter sequences of p/CIP/CARM1 targets supports the idea that transcription factor crosstalk is likely the favored mechanism of E2-dependent p/CIP/CARM1 complex recruitment. This data may have implications towards understanding the oncogenic role of p/CIP in breast cancer and ultimately allow for the identification of new prognostic indicators and/or viable therapeutic targets.

Level of nutrient intake affects mammary gland gene expression profiles in preweaned Holstein heifers

Bovine mammary parenchyma (PAR) and fat pad (MFP) development are responsive to preweaning level of nutrient intake. We studied transcriptome alterations in PAR and MFP from Holstein heifer calves (n=6/treatment) fed different nutrient intakes from birth to ca. 65 d age. Conventional nutrient intake received 441 g of dry matter (DM)/d of a control milk replacer (MR) [CON; 20% crude protein (CP), 20% fat, DM basis]. Calves in the accelerated nutrition groups received 951 g/d of high-protein/low-fat MR (HPLF; 28% CP, 20% fat, DM basis), 951 g/d of high-protein/high-fat MR (HPHF; 28% CP, 28% fat, DM basis), or 1,431 g/d of HPHF (HPHF+) MR. Out of 13,000 genes evaluated, over 1,500 differentially expressed genes (DEG) were affected (false discovery rate <0.10) by level of nutrient intake in PAR or MFP. Feeding HPLF versus CON resulted in the most dramatic changes in gene expression, with 278 and 588 DEG having ≥1.5-fold change in PAR and MFP. In PAR, the most-altered molecular functions were associated with metabolism of the cell (molecular transport and lipid metabolism) with most of the genes downregulated in HPLF versus CON. In MFP, DEG also were primarily associated with metabolism but changes also occurred in genes linked to cell morphology, cell-to-cell signaling, and immune response. Compared with CON, feeding HPHF or HPHF+ did not result in substantial additional effects on DEG beyond those observed with HPLF. The pentose phosphate, mitochondrial dysfunction, and ubiquinone biosynthesis pathways were among the most enriched due to HPLF versus CON in PAR and were inhibited, whereas glycosphingolipid biosynthesis, arachidonic acid metabolism, and eicosanoid synthesis pathways were among the most enriched due to HPLF versus CON in MFP and were inhibited. These responses suggest that, in PAR, doubling nutrient intake from standard feeding rates inhibited energy metabolism and activity of oxidative pathways that partly serve to protect cells against oxidative stress. The MFP in those heifers appeared to decrease production of lipid-derived metabolites that may play roles in signaling pathways within the adipocyte. Overall, results indicated that prepubertal/preweaned mammary transcriptome is responsive to long-term enhanced nutrient supply to achieve greater growth rates before weaning. The biological significance of these results to future milk production remains to be elucidated.

Neospora caninum: comparative gene expression profiling of Neospora caninum wild type and a temperature sensitive clone

To understand the genetic basis of virulence, gene expression profiles of a temperature-sensitive clone (NCts-8, relatively avirulent) and its wild type (NC-1) of Neospora caninum were characterized and compared using a high-density microarray with approximately 63,000 distinct oligonucleotides. This microarray consists of 5692 unique N. caninum sequences, including 1980 Tentative Consensus sequences and 3712 singleton ESTs from the TIGR N. caninum Gene Index (NCGI, release 5.0). Each sequence was represented by 11 distinct 60mer oligonucleotides synthesized in situ on the microarray. The results showed that 111 genes were significantly repressed and no up-regulated genes were identified in the NCts-8 clone. The level of 10 randomly selected genes from the repressed genes was confirmed using real-time RT-PCR. Of the 111 repressed genes, 58 were hypothetical protein products and 53 were annotated genes. Over 70% of the repressed genes identified in this study are clustered on five chromosomes (I, VII, VIII, X and XII). These results suggest that the down-regulated genes may be in part responsible for the reduced pathogenesis of NCts-8; further characterization of the regulated genes may aid in understanding of molecular basis of virulence and development of countermeasures against neosporosis.

The vitamin D receptor and inducible nitric oxide synthase associated pathways in acquired resistance to Cooperia oncophora infection in cattle

Cooperia oncophora is an economically important gastrointestinal nematode in ruminants. Acquired resistance to Cooperia oncophora infection in cattle develops rapidly as a result of prior infections. Naïve cattle, when given a primary infection of high-dose infective L3 larvae, develop a strong immunity to subsequent reinfection. Compared to primary infection, reinfection resulted in a marked reduction in worm establishment. In order to understand molecular mechanisms underlying the development of acquired resistance, we characterized the transcriptomic responses of the bovine small intestine to a primary infection and reinfection. A total of 23 pathways were significantly impacted during infection. The vitamin D receptor activation was strongly induced only during reinfection, suggesting that this pathway may play an important role in the development of acquired resistance via its potential roles in immune regulation and intestinal mucosal integrity maintenance. The expression of inducible nitric oxide synthase (NOS2) was strongly induced during reinfection but not during primary infection. As a result, several canonical pathways associated with NOS2 were impacted. The genes involved in eicosanoid synthesis, including prostaglandin synthase 2 (PTGS2 or COX2), remained largely unchanged during infection. The rapid development of acquired resistance may help explain the lack of relative pathogenicity by Cooperia oncophora infection in cattle. Our findings facilitate the understanding of molecular mechanisms underlying the development of acquired resistance, which could have an important implication in vaccine design.

Mastitis associated transcriptomic disruptions in cattle

Mastitis is ranked as the top disease for dairy cattle based on traditional cost analysis. Greater than 100 organisms from a broad phylogenetic spectrum are able to cause bovine mastitis. Transcriptomic characterization facilitates our understanding of host-pathogen relations and provides mechanistic insight into host resistance to mastitis. In this review, we discuss effector mechanisms and transcriptomic changes within the mammary gland in response to experimental infections. We compare temporal, spatial and pathogen-specific local transcriptomic disruptions in the mammary gland as well as pathogen-induced systemic responses and transcriptional changes in distant organs. We attempt to explain why studies on transcriptomic changes during critical physiological periods and in response to non-mastitic pathogens may have important implications for mastitis studies. Future perspectives on revealing bidirectional molecular cross-talk between mastitis pathogens and host cells using cutting-edge genomic technologies are also discussed.

Localized complement activation in the development of protective immunity against Ostertagia ostertagi infections in cattle

The abomasal nematode Ostertagia ostertagi is a major causal agent contributing to production inefficiencies in the cattle industry in temperate regions of the world. Protective immunity to infections develops very slowly and resistance to reinfection manifests only after prolonged exposure. Mechanisms underlying the development of protective immunity remain largely unexplored. Immune animals, which have significantly reduced worm burdens, were developed after multiple drug-attenuated experimental infections and were compared to a primary infected group and their respective uninfected controls. In this study, transcriptomic analysis identified three signaling pathways significantly impacted during both primary and repeat infections, the complement system, leukocyte extravasation and acute phase responses. Increased mRNA levels of complement components C3, factor B (CFB) and factor I (CFI) in the abomasal mucosa of the infected cattle were confirmed using quantitative PCR while Western blot analysis established the presence of elevated levels of activated C3 proteins in the mucosa. One of the initiators of local complement activation could be related to secretory IgA and IgM because infections significantly up-regulated expression of J chain (IGJ), as well as polymeric Ig receptor (PIGR) and an IgM-specific receptor (FAIM3), suggesting sustained increases in both synthesis and transepithelial transport of IgA and IgM during the infection. The elevated levels of pro-inflammatory cytokines, such as IL-4 and IL-1β, during infection may be involved in gene regulation of complement components. Our results suggest enhanced tissue repair and mucin secretion in immune animals may also contribute to protective immunity. These results are the first evidence that local complement activation may be involved in the development of long-term protective immunity and provide a novel mechanistic insight into resistance against O. ostertagi in cattle.

Microarray analysis of equine endometrium at days 8 and 12 of pregnancy

Establishment and maintenance of pregnancy in equids is only partially understood. To provide new insights into early events of this process, we performed a systematic analysis of transcriptome changes in the endometrium at Days 8 and 12 of pregnancy. Endometrial biopsy samples from pregnant and nonpregnant stages were taken from the same mares. Composition of the collected biopsy samples was analyzed using quantitative stereological techniques to determine proportions of surface and glandular epithelium and blood vessels. Microarray analysis did not reveal detectable changes in gene expression at Day 8, whereas at Day 12 of pregnancy 374 differentially expressed genes were identified, 332 with higher and 42 with lower transcript levels in pregnant endometrium. Expression of selected genes was validated by quantitative real-time RT-PCR. Gene set enrichment analysis, functional annotation clustering, and cocitation analysis were performed to characterize the genes differentially expressed in Day 12 pregnant endometrium. Many known estrogen-induced genes and genes involved in regulation of estrogen signaling were found, but also genes known to be regulated by progesterone and prostaglandin E2. Additionally, differential expression of a number of genes related to angiogenesis and vascular remodeling suggests an important role of this process. Furthermore, genes that probably have conserved functions across species, such as CRYAB, ERRFI1, FGF9, IGFBP2, NR2F2, STC1, and TNFSF10, were identified. This study revealed the potential target genes and pathways of conceptus-derived estrogens, progesterone, and prostaglandin E2 in the equine endometrium probably involved in the early events of establishment and maintenance of pregnancy in the mare.

Functional and gene network analyses of transcriptional signatures characterizing pre-weaned bovine mammary parenchyma or fat pad uncovered novel inter-tissue signaling networks during development

BACKGROUND

The neonatal bovine mammary fat pad (MFP) surrounding the mammary parenchyma (PAR) is thought to exert proliferative effects on the PAR through secretion of local modulators of growth induced by systemic hormones. We used bioinformatics to characterize transcriptomics differences between PAR and MFP from approximately 65 d old Holstein heifers. Data were mined to uncover potential crosstalk through the analyses of signaling molecules preferentially expressed in one tissue relative to the other.

RESULTS

Over 9,000 differentially expressed genes (DEG; False discovery rate <or= 0.05) were found of which 1,478 had a >or=1.5-fold difference between PAR and MFP. Within the DEG highly-expressed in PAR vs. MFP (n = 736) we noted significant enrichment of functions related to cell cycle, structural organization, signaling, and DNA/RNA metabolism. Only actin cytoskeletal signaling was significant among canonical pathways. DEG more highly-expressed in MFP vs. PAR (n = 742) belong to lipid metabolism, signaling, cell movement, and immune-related functions. Canonical pathways associated with metabolism and signaling, particularly immune- and metabolism-related were significantly-enriched. Network analysis uncovered a central role of MYC, TP53, and CTNNB1 in controlling expression of DEG highly-expressed in PAR vs. MFP. Similar analysis suggested a central role for PPARG, KLF2, EGR2, and EPAS1 in regulating expression of more highly-expressed DEG in MFP vs. PAR. Gene network analyses revealed putative inter-tissue crosstalk between cytokines and growth factors preferentially expressed in one tissue (e.g., ANGPTL1, SPP1, IL1B in PAR vs. MFP; ADIPOQ, IL13, FGF2, LEP in MFP vs. PAR) with DEG preferentially expressed in the other tissue, particularly transcription factors or pathways (e.g., MYC, TP53, and actin cytoskeletal signaling in PAR vs. MFP; PPARG and LXR/RXR Signaling in MFP vs. PAR).

CONCLUSIONS

Functional analyses underscored a reciprocal influence in determining the biological features of MFP and PAR during neonatal development. This was exemplified by the potential effect that the signaling molecules (cytokines, growth factors) released preferentially (i.e., more highly-expressed) by PAR or MFP could have on molecular functions or signaling pathways enriched in the MFP or PAR. These bidirectional interactions might be required to coordinate mammary tissue development under normal circumstances or in response to nutrition.

Profile of estrogen-responsive genes in an estrogen-specific mammary gland outgrowth model

Both ovarian and pituitary hormones are required for the pubertal development of the mouse mammary gland. Estradiol directs ductal elongation and branching, while progesterone leads to tertiary branching and alveolar development. The purpose of this investigation was to identify estrogen-responsive genes associated with pubertal ductal growth in the mouse mammary gland in the absence of other ovarian hormones and at different stages of development. We hypothesized that the estrogen-induced genes and their associated functions at early stages of ductal elongation would be distinct from those induced after significant ductal elongation had occurred. Therefore, ovariectomized prepubertal mice were exposed to 17beta-estradiol from two to 28 days, and mammary gland global gene expression analyzed by microarray analysis at various times during this period. We found that: (a) gene expression changes in our estrogen-only model mimic those changes that occur in normal pubertal development in intact mice, (b) both distinct and overlapping gene profiles were observed at varying extents of ductal elongation, and (c) cell proliferation, the immune response, and metabolism/catabolism were the most common functional categories associated with mammary ductal growth. Particularly striking was the novel observation that genes active during carbohydrate metabolism were rapidly and robustly decreased in response to estradiol. Lastly, we identified mammary estradiol-responsive genes that are also co-expressed with estrogen receptor alpha in human breast cancer. In conclusion, our genomic data support the physiological observation that estradiol is one of the primary hormonal signals driving ductal elongation during pubertal mammary development.

The genome response to artificial selection: a case study in dairy cattle

Dairy cattle breeds have been subjected over the last fifty years to intense artificial selection towards improvement of milk production traits. In this study, we performed a whole genome scan for differentiation using 42,486 SNPs in the three major French dairy cattle breeds (Holstein, Normande and Montbéliarde) to identify the main physiological pathways and regions which were affected by this selection. After analyzing the population structure, we estimated F(ST) within and across the three breeds for each SNP under a pure drift model. We further considered two different strategies to evaluate the effect of selection at the genome level. First, smoothing F(ST) values over each chromosome with a local variable bandwidth kernel estimator allowed identifying 13 highly significant regions subjected to strong and/or recent positive selection. Some of them contained genes within which causal variants with strong effect on milk production traits (GHR) or coloration (MC1R) have already been reported. To go further in the interpretation of the observed signatures of selection we subsequently concentrated on the annotation of differentiated genes defined according to the F(ST) value of SNPs localized close or within them. To that end we performed a comprehensive network analysis which suggested a central role of somatotropic and gonadotropic axes in the response to selection. Altogether, these observations shed light on the antagonism, at the genome level, between milk production and reproduction traits in highly producing dairy cows.

ASAS centennial paper: Lactation biology for the twenty-first century

Knowledge of general aspects of mammary gland function, including metabolic pathways and hormonal regulation of mammary gland development and lactation, in livestock species was obtained several decades ago. As basic biological information of growth factor action, apoptotic mechanisms, and signal transduction events has exploded, the mouse became the model of choice for studying fundamental mechanisms regulating mammary function. A complete sequenced genome also has made the mouse amenable for studies of mammary gene network expression. Advances in molecular biology techniques currently allow researchers to genetically modify mice to either overexpress or completely lack specific genes, thereby studying their function in an in vivo setting. Furthermore, the use of mammary-specific promoters has allowed genes related to mammary gland function to be eliminated from the mammary gland in a developmental and tissue-specific manner. These studies have demonstrated the complexity that underlies mammary gland development and function in rodents and may provide insight into the mechanisms that ultimately allow the ruminant or swine mammary gland to function in a coordinated fashion throughout puberty, pregnancy, lactation, and involution. The challenge facing animal scientists is how to obtain similar information in much larger and expensive livestock. One possible approach is to manipulate gene expression in vitro using mammary cell culture models derived from domestic animals (e.g., genes can be "knocked down" using small interfering RNA approaches). Ultimately, major advances in understanding lactation biology may come from coupling basic mechanistic information with functional genomics, proteomics, and metabolomics approaches. A strong foundation in bioinformatics will also be required to optimize use of these new technologies. Stem cell biology also represents an exciting area in the next decade that holds promise for improving lactation efficiency. Strong training of our future scientists in these areas should facilitate livestock-focused mammary gland research that will allow basic information to be gained at unprecedented amounts, ultimately leading to optimization of livestock production.

Motif-specific sampling of phosphoproteomes

Phosphoproteomics, the targeted study of a subfraction of the proteome which is modified by phosphorylation, has become an indispensable tool to study cell signaling dynamics. We described a methodology that linked phosphoproteome and proteome analysis based on Ba2+ binding properties of amino acids. This technology selected motif-specific phosphopeptides independent of the system under analysis. MudPIT (Multidimensional Identification Technology) identified 1037 precipitated phosphopeptides from as little as 250 microg of proteins. To extend coverage of the phosphoproteome, we sampled the nuclear extract of HeLa cells with three values of Ba2+ ions molarity. The presence of more than 70% of identified phosphoproteins was further substantiated by their nonmodified peptides. Upon isoproterenol stimulation of HEK cells, we identified an increasing number of phosphoproteins from MAPK cascades and AKAP signaling hubs. We quantified changes in both protein and phosphorylation levels of 197 phosphoproteins including a critical kinase, MAPK1. Integration of differential phosphorylation of MAPK1 with knowledge bases constructed modules that correlated well with its role as node in cross-talk of canonical pathways.

Pathway-based association analysis of genome-wide screening data suggest that genes associated with the gamma-aminobutyric acid receptor signaling pathway are involved in neuroleptic-induced, treatment-resistant tardive dyskinesia

OBJECTIVE

Neuroleptic-induced tardive dyskinesia (TD) is an involuntary movement disorder that develops in patients who have undergone long-term treatment with antipsychotic medications, and its etiology is unclear. In this study, a genome-wide association screening was done to identify the pathway(s) in which genetic variations influence susceptibility to neuroleptic-induced TD.

METHODS

Screening with Sentrix Human-1 Genotyping BeadChip (Illumina, San Diego, California, USA) was done for 50 Japanese schizophrenia patients with treatment-resistant TD and 50 Japanese schizophrenia patients without TD. A total of 40 573 single nucleotide polymorphisms that were not in linkage disequilibrium with each other and were located in the exonic and intronic regions of 13 307 genes were analyzed. After gene-based corrections, P values for allelic associations were subjected to canonical pathway-based analyses with Ingenuity Pathway Analysis software (Ingenuity Systems, Inc., Redwood City, California, USA).

RESULTS

Eight genes (ABAT, ALDH9A1, GABRA3, GABRA4, GABRB2, GABRAG3, GPHN, and SLC6A11) contained polymorphisms with gene-based corrected allelic P values of less than 0.05. They were aggregated significantly in 33 genes belonging to the gamma-aminobutyric acid (GABA) receptor signaling pathway (P=0.00007, corrected P=0.01). Associations were replicated in an independent sample of 36 patients with TD and 136 patients without TD for polymorphisms in SLC6A11 (GABA transporter 3) (P=0.0004 in the total sample), GABRB2 (beta-2 subunit of GABA-A receptor) (P=0.00007 in the total sample), and GABRG3 (gamma-3 subunit of GABA-A receptor) (P=0.0006 in the total sample).

CONCLUSION

The results suggest that the GABA receptor signaling pathway may be involved in genetic susceptibility to treatment-resistant TD, at least in a subgroup of Japanese patients with schizophrenia. The present results suggest that benzodiazepines may be considered as possible treatment option for TD.