Transcriptomic analysis of the stress response to weaning at housing in bovine leukocytes using RNA-seq technology

Can a gradual weaning and separation process reduce weaning distress in dam-reared dairy calves? A comparison with the 2-step method

The weaning and separation phase remains one of the biggest challenges for cow-calf contact systems, but a gradual process that better mimics the naturally occurring reduction in milk intake has not yet been scientifically investigated. Therefore, the aim of our study was to compare behavioral and physiological indicators of distress in 3-mo-old dam-reared dairy calves (with previous full-time cow-calf contact) weaned and separated either via gradual reduction of contact time with the dam (GR; 1 wk of half-day contact, 1 wk of morning contact, and 1 wk of fence-line contact before complete separation, n = 18) or via 2-step weaning using a nose flap (NF, 2 wk of access to the dam with a nose flap, 1 wk of fence-line contact before complete separation, n = 18). Behavior was recorded 1 wk before (or for lying 3 wk before) weaning start and during the 3 wk weaning and separation period with direct observations on 4 d/wk or via accelerometers (locomotor play, lying behavior). Blood and fecal samples were taken twice per week from weaning start until 3 wk after weaning start. Calves were weighed weekly. Statistical analysis was conducted using (generalized) linear mixed models. Over the whole weaning and separation phase, NF calves showed a stronger decrease in the number of lying bouts, amount of locomotor play, and ADG, as well as a higher increase in TMR feeding time compared with GR calves, whereas GR calves vocalized more often and showed more searching behavior than NF calves. Also, the neutrophil:lymphocyte ratio of NF calves was elevated on d 3 after insertion of the nose flaps compared with baseline, but showed no change for GR calves on any sampling day. Overall, results point toward a favorable effect of a gradual weaning strategy on reduction of weaning and separation distress in dam-reared dairy calves, but the method requires further improvement from the protocol used in our study.

Resting and activated bovine neutrophils and eosinophils differ in their responses to adrenergic agonists

Polymorphonuclear cells (PMN) provide a rapid response to infection and tissue damage and stress can modify these critical innate immune defences. The study of adrenergic receptor (AR) expression and function in bovine PMNs is limited but both neutrophils and eosinophils express numerous AR genes but differ significantly in their expression of individual AR genes. A flow cytometric technique was developed to differentiate between bovine neutrophils and eosinophils so both neutrophil and eosinophil responses to adrenergic agonists could be analysed. Neutrophils and eosinophils displayed significantly different changes in CD11b, L-selectin, and CD44 expression when activated by bovine serum opsonized zymosan and recombinant bovine interferon gamma. The responses of activated and resting neutrophils and eosinophils were then compared following stimulation with endogenous adrenergic agonists, epinephrine (E) norepinephrine (NE), and synthetic agonists targeting α1-, α2-, or β-ARs. Both resting and activated neutrophils and eosinophils displayed differences in iROS, CD44, and L-selectin expression following stimulation with E and NE. Resting neutrophils displayed pro-inflammatory responses to both E and NE, while resting eosinophils displayed a pro-inflammatory response to only NE. No single synthetic adrenergic agonist fully recapitulated responses observed with either E or NE and responses to adrenergic agonists were dose-dependent. In conclusion, bovine eosinophils and neutrophils responded to multiple adrenergic agonists by altering expression of proteins involved in immune surveillance and pro-inflammatory responses. Significant differences in neutrophil and eosinophil responses to adrenergic agonists are consistent with their differences in AR gene expression. This highlights the importance of analysing separately these two PMN subpopulations when investigating the effects of either endogenous or synthetic AR agonists.

Effects of dietary supplementation of a blend of Saccharomyces cerevisiae, multiple live probiotic bacteria, and their fermentation products on performance, health, and rumen bacterial community of newly weaned beef steers during a 56-d receiving period

We examined the effects of a blend of Saccharomyces cerevisiae, multiple live probiotic bacteria, and their fermentation products on performance, health, and the ruminal bacterial community of newly weaned beef steers during a 56-d receiving period. Forty newly weaned Angus crossbred steers (221 ± 25.6 kg BW; 180 ± 17 d of age) were stratified by body weight (BW) into four pens (10 steers per pen) such that each pen had a similar average BW at the beginning of the experiment. The pens were randomly assigned to receive a corn silage basal diet (CON; n = 20) or the basal diet supplemented with 9 g/steer/d of PRO feed additive (PRO; n = 20). The PRO additive is a blend of S. cerevisiae and the fermentation products of Enterococcus faecium, Bacillus licheniformis, B. subtilis, Lactobacillus animalis, and Propionibacterium freudenreichii. The DMI and water consumed were monitored using the GrowSafe intake nodes and custom flow meters, respectively. BWs were recorded weekly to calculate average daily gain (ADG). Before morning feeding, 10 mL of blood was taken from each steer on days 0-7, and thereafter weekly for analyses of immune cells, plasma glucose, and NEFAs. On day 56, rumen fluid samples (200 mL each) were collected from all the steers for microbiome analysis. Over the 56-d receiving period, the supplemental PRO had no effects on DMI, water intake, or ADG. However, compared to CON, beef steers fed supplemental PRO tended to have greater ADG (P = 0.08) and BW (P = 0.07) during the first 14 d of the study. There was a treatment × day interaction (P ≤ 0.05) for WBC, neutrophils and monocytes over the 56 d such that beef steers fed supplemental PRO had lower blood concentrations on certain days during the first 7 d after weaning, indicating reduced inflammation or stress response. The results of the rumen microbiome analysis revealed that the relative abundance of complex fiber degrading or obligate proton-reducing bacterial genera such as Bacteroides, Ruminococcus gauvreauii group, Desulfovibrio, Syntrophococcus, and Acetitomaculum were greater (P ≤ 0.05) in beef steers fed supplemental PRO compared to CON. This study demonstrated that dietary supplementation of PRO improved the growth performance, reduced stress or inflammatory response during the initial days after weaning, and altered the ruminal bacterial community toward increased relative abundance of bacterial genera associated with improved rumen function.

SWATH-MS Analysis of Blood Plasma and Circulating Small Extracellular Vesicles Enables Detection of Putative Protein Biomarkers of Fertility in Young and Aged Dairy Cows

The development of biomarkers of fertility could provide benefits for the genetic improvement of dairy cows. Circulating small extracellular vesicles (sEVs) show promise as diagnostic or prognostic markers since their cargo reflects the metabolic state of the cell of origin; thus, they mirror the physiological status of the host. Here, we employed data-independent acquisition mass spectrometry to survey the plasma and plasma sEV proteomes of two different cohorts of Young (Peripubertal; n = 30) and Aged (Primiparous; n = 20) dairy cows (Bos taurus) of high- and low-genetic merit of fertility and known pregnancy outcomes (ProteomeXchange data set identifier PXD042891). We established predictive models of fertility status with an area under the curve of 0.97 (sEV; p value = 3.302e-07) and 0.95 (plasma; p value = 6.405e-08). Biomarker candidates unique to high-fertility Young cattle had a sensitivity of 0.77 and specificity of 0.67 (*p = 0.0287). Low-fertility biomarker candidates uniquely identified in sEVs from Young and Aged cattle had a sensitivity and specificity of 0.69 and 1.0, respectively (***p = 0.0005). Our bioinformatics pipeline enabled quantification of plasma and circulating sEV proteins associated with fertility phenotype. Further investigations are warranted to validate this research in a larger population, which may lead to improved classification of fertility status in cattle.

Differential gene expression analysis using RNA-seq in the blood of goats exposed to transportation stress

Transportation stress causes significant changes in physiological responses in goats; however, studies exploring the transcriptome of stress are very limited. The objective of this study was to determine the differential gene expressions and related pathways in the blood samples using RNA-seq procedure in Spanish goats subjected to different durations of transportation stress. Fifty-four male Spanish goats (8-mo old; BW = 29.7 ± 2.03 kg) were randomly subjected to one of three treatments (TRT; n = 18 goats/treatment): (1) transported for 180 min, (2) transported for 30 min, or (3) held in pens (control). Blood samples were collected before and after treatment for stress hormone, metabolite, and transcriptomic analysis. RNA-seq technology was used to obtain the transcriptome profiles of blood. Analysis of physiological data using SAS showed that plasma cortisol concentrations were higher (P < 0.01) in 180 min and 30 min groups compared to the control group. Enrichment analysis of DEGs related to transportation stress through Gene Ontology and KEGG databases revealed that the differentially expressed genes related to inflammatory pathways, caspases, and apoptosis such as IL1R2, CASP14, CD14, TLR4, and MAPK14 were highly enriched in the transported group of goats compared to non-transported goats. Stress in goats leads to a sequence of events at cellular and molecular levels that causes inflammation and apoptosis.

Asinine milk mitigates stress-mediated immune, cortisol and behavioral responses of piglets to weaning: A study to foster future interventions in humans

Introduction

The present study assessed whether asinine milk supplementation improved the immune and behavioral responses of piglets during an early life weaning stress event as a model for its future use in humans.

Methods

For this, 48 piglets from 4 different litters were used. At 20 days of age, piglets were weighed and allocated with their litter and dam into group pens until 28 days of age. Four piglets from each litter were then randomly assigned to either (1) asinine milk supplementation (n = 16) (2), skimmed cow milk supplementation (n = 16) or (3) no supplementation (n = 16; control group). The supplementations were voluntarily administered for 3 days preweaning and 3 days postweaning using a baby bottle. The effects on the weaning stress response were assessed through salivary cortisol measurements; behavioral tests such as the open field, novel object end elevated plus maze tests; and gene expression of HSD11B1, NR3C1 and IL1B in PBMCs, which was determined by RT-qPCR and normalized to GAPDH and UBB. To test the effect of the supplementations on weight, milk intake, gene expression, and behavior, a randomized block design was used with repeated measurements over time by the PROC MIXED procedure.

Results and discussion

The effects on salivary cortisol were determined using the ratio between the morning and afternoon concentrations, considering the time before and after the weaning event. Principal component analysis (PCA) and Fisher's test were performed to evaluate the behavior test data. When comparing salivary cortisol concentrations between the pre- and postweaning periods, there was a difference (p < 0.05) between the supplementation groups in the afternoon period, suggesting that piglets fed asinine milk had lower afternoon cortisol concentrations postweaning than their counterparts. For the behavioral tests, the supplementations had no measurable effects. No difference was between groups pre- and postweaning for the expression of HSD11B2, which codes for an enzyme that breaks down cortisol. However, the expression of NR3C1, which encodes the glucocorticoid receptor, was significantly upregulated in piglets supplemented with cow milk (mean 1.245; p < 0.05).

Conclusion

Asinine milk downregulated 1L1B gene expression, which codes for an inflammatory cytokine. In conclusion, these results suggest that supplementation with asinine milk may represent a strategy to diminish the damage associated with an early life event by modulating IL1B expression and reducing salivary cortisol levels in piglets undergoing weaning stress. Further transcriptomic and metabolomic studies may improve our understanding of the molecular pathways that mediate this systemic immune-mediated response.

Systems biology under heat stress in Indian cattle

Transcriptome profiling of Vrindavani and Tharparkar cattle (n = 5 each) revealed that more numbers of genes were dysregulated in Vrindavani than in Tharparkar. A contrast in gene expression was observed with 18.9 % of upregulated genes in Vrindavani downregulated in Tharparkar and 17.8% upregulated genes in Tharparkar downregulated in Vrindavani. Functional annotation of genes differentially expressed in Tharparkar and Vrindavani revealed that the systems biology in Tharparkar is moving towards counteracting the effects due to heat stress. Unlike Vrindavani, Tharparkar is not only endowed with higher expression of the scavengers (UBE2G1, UBE2S, and UBE2H) of misfolded proteins but also with protectors (VCP, Serp1, and CALR) of naïve unfolded proteins. Further, higher expression of the antioxidants in Tharparkar enables it to cope up with higher levels of free radicals generated as a result of heat stress. In this study, we found relevant genes dysregulated in Tharparkar in the direction that can counter heat stress.

RNA-Seq Bayesian Network Exploration of Immune System in Bovine

Background:

The stress is one of main factors effects on production system. Several factors (both genetic and environmental elements) regulate immune response to stress.

Objectives:

In order to determine the major immune system regulatory genes underlying stress responses, a learning Bayesian network approach for those regulatory genes was applied to RNA-Seq data from a bovine leukocyte model system.

Material and Methods:

The transcriptome dataset GSE37447 was used from GEO and a Bayesian network on differentially expressed genes was learned to investigate the gene regulatory network.

Results:

Applying the method produced a strongly interconnected network with four genes (TERF2IP, PDCD10, DDX10 and CENPE) acting as nodes, suggesting these genes may be important in the transcriptome regulation program of stress response. Of these genes TERF2IP has been shown previously to regulate gene expression, act as a regulator of the nuclear factor-kappa B (NF-κB) signalling, and to activate expression of NF-κB target genes; PDCD10 encodes a conserved protein associated with cell apoptosis; DDX10 encodes a DEAD box protein and is believed to be associated with cellular growth and division; and CENPE involves unstable spindle microtubule capture at kinetochores. Together these genes are involved in DNA damage of apoptosis, RNA splicing, DNA repairing, and regulating cell division in the bovine genome. The topology of the learned Bayesian gene network indicated that the genes had a minimal interrelationship with each other. This type of structure, using the publically available computational tool, was also observed on human orthologous genes of the differentially expressed genes.

Conclusions:

Overall, the results might be used in transcriptomic-assisted selection and design of new drug targets to treat stress-related problems in bovines.

Transcriptome analysis of the effect of pyrroloquinoline quinone disodium (PQQ·Na2) on reproductive performance in sows during gestation and lactation

Background

Pyrroloquinoline quinone (PQQ), which is a water soluble, thermo-stable triglyceride-quinone, is widely distributed in nature and characterized as a mammalian vitamin-like redox cofactor. The objective of this study was to investigate the effects of pyrroloquinoline quinone disodium (PQQ·Na₂) on reproductive performance in sows.

Results

Dietary supplementation with PQQ·Na₂ significantly increased the total number of piglets born, the number of piglets born alive and the born alive litter weight. It also increased the antioxidant status in the placenta, plasma and milk. The concentration of NO was significantly increased in the plasma and placenta. RNA-seq analysis showed that 462 unigenes were differentially expressed between the control (Con) treatment and PQQ treatment groups. Among these unigenes, 199 were upregulated, while 263 unigenes were downregulated. The assigned functions of the unigenes covered a broad range of GO categories. Reproduction (27, 7.03%) and the reproduction process (27, 7.03%) were assigned to the biological process category. By matching DEGs to the KEGG database, we identified 29 pathways.

Conclusions

In conclusion, dietary supplementation with PQQ·Na₂ in gestating and lactating sows had positive effects on their reproductive performance.

Electronic supplementary material

The online version of this article (10.1186/s40104-019-0369-y) contains supplementary material, which is available to authorized users.

Effect of supplementation with n-3 polyunsaturated fatty acids and/or β-glucans on performance, feeding behaviour and immune status of Holstein Friesian bull calves during the pre- and post-weaning periods

Background

Previous research in both calves and other species has suggested n-3 polyunsaturated fatty acids (PUFA) and β-glucans may have positive effects on immune function. This experiment measured performance, behaviour, metabolite and immunological responses to pre-weaning supplementation of dairy bull calves with n-3 PUFA in the form of fish oil and β-glucans derived from seaweed extract. 44 Holstein Friesian bull calves, aged 13.7 ± 2.5 d and weighing 48.0 ± 5.8 kg were artificially reared using an electronic feeding system. Each calf was offered 5 L (120 g/L) per day of milk replacer (MR) and assigned to one of four treatments included in the MR, (1) Control (CON); (2) 40 g n-3 PUFA per day (FO); (3) 1 g β-glucans per day (GL) and (4) 40 g n-3 PUFA per day & 1 g/d β-glucans (FOGL) in a 2 × 2 factorial design. Milk replacer and concentrate was offered from d 0–62 (pre-weaning), while concentrate provision continued for a further 31 d post-weaning period. Individual daily feed intake and feeding behaviour was recorded throughout, while bodyweight and blood analyte data were collected at regular intervals.

Results

Overall mean concentrate DMI from d 0–93 was 1.39, 1.27, 1.00 and 0.72 kg/d for CON, FO, GL and FOGL calves, respectively (SEM = 0.037; P < 0.0001). Calves supplemented with GL were significantly lighter (P < 0.0001) at both weaning (d 62) and turnout to pasture (d 93) than un-supplemented calves, with a similar effect (P < 0.0001) evident for calves receiving FO compared to un-supplemented contemporaries. Supplementation with GL reduced the number of unrewarded visits where milk was not consumed (P < 0.0001) while supplementation with FO increased mean drinking speed (P < 0.0001). Supplementation with GL resulted in greater concentrations of haptoglobin (P = 0.034), greater serum osmolality (P = 0.021) and lower lymphocyte levels (P = 0.027). In addition, cells from GL supplemented calves exhibited a lower response than un-supplemented contemporaries to both Phytohaemagglutinin A stimulated IFN-γ (P = 0.019) and Concanavalin A stimulated IFN-γ (P = 0.012) following in vitro challenges.

Conclusions

Pre-weaning supplementation of bull calves with either n-3 PUFA or β-glucan resulted in reduced voluntary feed intake of concentrate and consequently poorer pre-weaning calf performance. There was no evidence for any beneficial effect of either supplementation strategy on calves’ immune responses.

An intensive milk replacer feeding program benefits immune response and intestinal microbiota of lambs during weaning

Background

Pre-weaning milk replacer (MR) feeding program is a key factor affecting the health and welfare of lambs during their weaning. Weaning stress is well known as an inducement that negatively impacts the immune system of young ruminants, whose physiological and immune state is closely linked to the community of microbiota in their intestines. This study had two objectives: 1) To evaluate the innate immune response to weaning stress at both the physiological and molecular level; 2) To investigate changes to the jejunal chyme and mucosal adhesive microbiota between the control and high plane of MR groups.

Results

In this experiment, the plasma concentrations of cortisol, norepinephrine (NE) and tumor necrosis factor-α (TNFα) were higher in the C than the H group (P < 0.05), as was the expression of pro-inflammatory cytokines such as TNFα and CXCL8 (P < 0.05) in plasma. In jejunal tissue, the expression of TLR4 and TNFα were also higher in the C group (P < 0.01); histopathology showed the H group had lower lymphocyte infiltration. In the C group, however, major pathological changes were associated with extensive infiltration of lymphocytes, eosinophils, and neutrophils. Principal component analysis indicated the lamb immune response was influenced by weaning stress and modulated by the MR treatments. 16S-rRNA sequencing was used to evaluate jejunal mucosa and chyme bacterial diversity and composition. The C group’s chyme had a greater alpha index (ACE: P = 0.095; Chao1: P = 0.085) than H group. In jejunal mucosa, the relative abundance of Plesiomonas was 4-fold higher (P = 0.017) in the C than the H group.

Conclusions

This study’s results revealed that weaning stress induced alterations to the lambs’ immune system that lasted beyond the 21 d measured, and that a long-term inflammatory response effect was evidenced by changes in their hematological and expressed pro-inflammatory cytokines. Pre-weaning with a differing MR allowance resulted in complicated biological responses and compositional changes to the lambs’ jejunal microbiota. Clearly, an intensive MR feeding program induced a milder immunity response and lower relative abundance of pathogenic bacteria when compared with the traditional feeding program.

Impact of energy restriction during late gestation on the muscle and blood transcriptome of beef calves after preconditioning

Background

Maternal nutrition has been highlighted as one of the main factors affecting intra-uterine environment. The increase in nutritional requirements by beef cows during late gestation can cause nutritional deficiency in the fetus and impact the fetal regulation of genes associated with myogenesis and immune response.

Methods

Forty days before the expected calving date, cows were assigned to one of two diets: 100% (control) or 70% (restricted group) of the daily energy requirement. Muscle samples were collected from 12 heifers and 12 steers, and blood samples were collected from 12 steers. The objective of this work was to identify and to assess the biological relevance of differentially expressed genes (DEG) in the skeletal muscle and blood of beef calves born from cows that experienced [or not] a 30% energy restriction during the last 40 days of gestation.

Results

A total of 160, 164, and 346 DEG (q-value< 0.05) were identified in the skeletal muscle for the effects of diet, sex, and diet-by-sex interaction, respectively. For blood, 452, 1392, and 155 DEG were identified for the effects of diet, time, and diet-by-time interaction, respectively. For skeletal muscle, results based on diet identified genes involved in muscle metabolism. In muscle, from the 10 most DEG down-regulated in the energy-restricted group (REST), we identified 5 genes associated with muscle metabolism and development: SLCO3A1, ATP6V0D1, SLC2A1, GPC4, and RASD2. In blood, among the 10 most DEG, we found genes related to response to stress up-regulated in the REST after weaning, such as SOD3 and INO80D, and to immune response down-regulated in the REST after vaccination, such as OASL, KLRF1, and LOC104968634.

Conclusion

In conclusion, maternal energy restriction during late gestation may limit the expression of genes in the muscle and increase expression in the blood of calves. In addition, enrichment analysis showed that a short-term maternal energy restriction during pregnancy affects the expression of genes related to energy metabolism and muscle contraction, and immunity and stress response in the blood. Therefore, alterations in the intra-uterine environment can modify prenatal development with lasting consequences to adult life.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5089-8) contains supplementary material, which is available to authorized users.

RNA Sequencing (RNA-Seq) Reveals Extremely Low Levels of Reticulocyte-Derived Globin Gene Transcripts in Peripheral Blood From Horses (Equus caballus) and Cattle (Bos taurus)

RNA-seq has emerged as an important technology for measuring gene expression in peripheral blood samples collected from humans and other vertebrate species. In particular, transcriptomics analyses of whole blood can be used to study immunobiology and develop novel biomarkers of infectious disease. However, an obstacle to these methods in many mammalian species is the presence of reticulocyte-derived globin mRNAs in large quantities, which can complicate RNA-seq library sequencing and impede detection of other mRNA transcripts. A range of supplementary procedures for targeted depletion of globin transcripts have, therefore, been developed to alleviate this problem. Here, we use comparative analyses of RNA-seq data sets generated from human, porcine, equine, and bovine peripheral blood to systematically assess the impact of globin mRNA on routine transcriptome profiling of whole blood in cattle and horses. The results of these analyses demonstrate that total RNA isolated from equine and bovine peripheral blood contains very low levels of globin mRNA transcripts, thereby negating the need for globin depletion and greatly simplifying blood-based transcriptomic studies in these two domestic species.

Construction of an integrated gene regulatory network link to stress-related immune system in cattle

The immune system is an important biological system that is negatively impacted by stress. This study constructed an integrated regulatory network to enhance our understanding of the regulatory gene network used in the stress-related immune system. Module inference was used to construct modules of co-expressed genes with bovine leukocyte RNA-Seq data. Transcription factors (TFs) were then assigned to these modules using Lemon-Tree algorithms. In addition, the TFs assigned to each module were confirmed using the promoter analysis and protein-protein interactions data. Therefore, our integrated method identified three TFs which include one TF that is previously known to be involved in immune response (MYBL2) and two TFs (E2F8 and FOXS1) that had not been recognized previously and were identified for the first time in this study as novel regulatory candidates in immune response. This study provides valuable insights on the regulatory programs of genes involved in the stress-related immune system.

Breed and adaptive response modulate bovine peripheral blood cells' transcriptome

BACKGROUND

Adaptive response includes a variety of physiological modifications to face changes in external or internal conditions and adapt to a new situation. The acute phase proteins (APPs) are reactants synthesized against environmental stimuli like stress, infection, inflammation.

METHODS

To delineate the differences in molecular constituents of adaptive response to the environment we performed the whole-blood transcriptome analysis in Italian Holstein (IH) and Italian Simmental (IS) breeds. For this, 663 IH and IS cows from six commercial farms were clustered according to the blood level of APPs. Ten extreme individuals (five APP+ and APP- variants) from each farm were selected for the RNA-seq using the Illumina sequencing technology. Differentially expressed (DE) genes were analyzed using dynamic impact approach (DIA) and DAVID annotation clustering. Milk production data were statistically elaborated to assess the association of APP+ and APP- gene expression patterns with variations in milk parameters.

RESULTS

The overall de novo assembly of cDNA sequence data generated 13,665 genes expressed in bovine blood cells. Comparative genomic analysis revealed 1,152 DE genes in the comparison of all APP+ vs. all APP- variants; 531 and 217 DE genes specific for IH and IS comparison respectively. In all comparisons overexpressed genes were more represented than underexpressed ones. DAVID analysis revealed 369 DE genes across breeds, 173 and 73 DE genes in IH and IS comparison respectively. Among the most impacted pathways for both breeds were vitamin B6 metabolism, folate biosynthesis, nitrogen metabolism and linoleic acid metabolism.

CONCLUSIONS

Both DIA and DAVID approaches produced a high number of significantly impacted genes and pathways with a narrow connection to adaptive response in cows with high level of blood APPs. A similar variation in gene expression and impacted pathways between APP+ and APP- variants was found between two studied breeds. Such similarity was also confirmed by annotation clustering of the DE genes. However, IH breed showed higher and more differentiated impacts compared to IS breed and such particular features in the IH adaptive response could be explained by its higher metabolic activity. Variations of milk production data were significantly associated with APP+ and APP- gene expression patterns.

Osteocalcin and Sex Hormone Binding Globulin Compete on a Specific Binding Site of GPRC6A

The undercarboxylated form of osteocalcin (ucOC) regulates male fertility and energy metabolism, acting through the G protein-coupled receptor (GPRC)6A, thus forming a new pancreas-bone-testis axis. Recently, GPRC6A has also been suggested to mediate the nongenomic responses of free testosterone (T). However, these data did not consider the physiological scenario, where circulating T is mainly bound to sex hormone-binding globulin (SHBG) and only a small percentage circulates freely in the blood. Here, by the use of computational modelling, we document the existence of similar structural moieties between ucOC and SHBG that are predicted to bind to GPRC6A at docking analysis. This hypothesis of competition was assessed by binding experiments on human embryonic kidney-293 cells transfected with human GPRC6A gene. Unliganded SHBG specifically bound the membrane of human embryonic kidney-293 cells transfected with GPRC6A and was displaced by ucOC when coincubated at 100-fold molar excess. Furthermore, specific downstream Erk1/2 phosphorylation after stimulation of GPRC6A with ucOC was significantly blunted by 100-fold molar excess of unliganded SHBG. Intriguingly previous incubation with unliganded SHBG, followed by incubation with T, induced Erk1/2 phosphorylation in a dose-dependent manner. Neither binding nor stimulating activities were shown for SHBG saturated with T. Experiments on mutation constructs of GPRC6A strengthened the hypothesis of a common binding site of ucOC and SHBG. Given the role of GPRC6A on energy metabolism, these data agree with epidemiological association between SHBG levels and insulin sensitivity, suggest GPRC6A as a likely SHBG receptor, and add bases for the possible regulation of androgen activity in a nonsteroidal manner.

Adrenocortical Expression Profiling of Cattle with Distinct Juvenile Temperament Types

Temperament affects ease of handling, animal welfare, and economically important production traits in cattle. The use of gene expression profiles as molecular traits provides a novel means of gaining insight into behavioural genetics. In this study, differences in adrenocortical expression profiles between 60 F₂ cows (Charolais × German Holstein) of distinct temperament types were analysed. The cows were assessed in a novel-human test at an age of 90 days. Most of the adrenal cortex transcripts which were differentially expressed (FDR <0.05) were found between temperament types of 'fearful/neophobic-alert' and all other temperament types. These transcripts belong to several biological functions like NRF2-mediated oxidative stress response, Glucocorticoid Receptor Signalling and Complement System. Overall, the present study provides new insight into transcriptional differences in the adrenal cortex between cows of distinct temperament types. Genetic regulations of such molecular traits facilitate uncovering positional and functional gene candidates for temperament type in cattle.

The developmental transcriptome landscape of bovine skeletal muscle defined by Ribo-Zero ribonucleic acid sequencing

Ribonucleic acid sequencing (RNA-Seq) libraries are normally prepared with oligo(dT) selection of poly(A)+ mRNA, but it depends on intact total RNA samples. Recent studies have described Ribo-Zero technology, a novel method that can capture both poly(A)+ and poly(A)- transcripts from intact or fragmented RNA samples. We report here the first application of Ribo-Zero RNA-Seq for the analysis of the bovine embryonic, neonatal, and adult skeletal muscle whole transcriptome at an unprecedented depth. Overall, 19,893 genes were found to be expressed, with a high correlation of expression levels between the calf and the adult. Hundreds of genes were found to be highly expressed in the embryo and decreased at least 10-fold after birth, indicating their potential roles in embryonic muscle development. In addition, we present for the first time the analysis of global transcript isoform discovery in bovine skeletal muscle and identified 36,694 transcript isoforms. Transcriptomic data were also analyzed to unravel sequence variations; 185,036 putative SNP and 12,428 putative short insertions-deletions (InDel) were detected. Specifically, many stop-gain, stop-loss, and frameshift mutations were identified that probably change the relative protein production and sequentially affect the gene function. Notably, the numbers of stage-specific transcripts, alternative splicing events, SNP, and InDel were greater in the embryo than in the calf and the adult, suggesting that gene expression is most active in the embryo. The resulting view of the transcriptome at a single-base resolution greatly enhances the comprehensive transcript catalog and uncovers the global trends in gene expression during bovine skeletal muscle development.

Decoding the Regulatory Landscape of Ageing in Musculoskeletal Engineered Tissues Using Genome-Wide DNA Methylation and RNASeq

Mesenchymal stem cells (MSC) are capable of multipotent differentiation into connective tissues and as such are an attractive source for autologous cell-based regenerative medicine and tissue engineering. Epigenetic mechanisms, like DNA methylation, contribute to the changes in gene expression in ageing. However there was a lack of sufficient knowledge of the role that differential methylation plays during chondrogenic, osteogenic and tenogenic differentiation from ageing MSCs. This study undertook genome level determination of the effects of DNA methylation on expression in engineered tissues from chronologically aged MSCs. We compiled unique DNA methylation signatures from chondrogenic, osteogenic, and tenogenic engineered tissues derived from young; n = 4 (21.8 years ± 2.4 SD) and old; n = 4 (65.5 years±8.3SD) human MSCs donors using the Illumina HumanMethylation 450 Beadchip arrays and compared these to gene expression by RNA sequencing. Unique and common signatures of global DNA methylation were identified. There were 201, 67 and 32 chondrogenic, osteogenic and tenogenic age-related DE protein-coding genes respectively. Findings inferred the nature of the transcript networks was predominantly for 'cell death and survival', 'cell morphology', and 'cell growth and proliferation'. Further studies are required to validate if this gene expression effect translates to cell events. Alternative splicing (AS) was dysregulated in ageing with 119, 21 and 9 differential splicing events identified in chondrogenic, osteogenic and tenogenic respectively, and enrichment in genes associated principally with metabolic processes. Gene ontology analysis of differentially methylated loci indicated age-related enrichment for all engineered tissue types in 'skeletal system morphogenesis', 'regulation of cell proliferation' and 'regulation of transcription' suggesting that dynamic epigenetic modifications may occur in genes associated with shared and distinct pathways dependent upon engineered tissue type. An altered phenotype in engineered tissues was observed with ageing at numerous levels. These changes represent novel insights into the ageing process, with implications for stem cell therapies in older patients. In addition we have identified a number of tissue-dependant pathways, which warrant further studies.

Characterisation of the Whole Blood mRNA Transcriptome in Holstein-Friesian and Jersey Calves in Response to Gradual Weaning

Weaning of dairy calves is an early life husbandry management practice which involves the changeover from a liquid to a solid feed based diet. The objectives of the study were to use RNA-seq technology to examine the effect of (i) breed and (ii) gradual weaning, on the whole blood mRNA transcriptome of artificially reared Holstein-Friesian and Jersey calves. The calves were gradually weaned over 14 days (day (d) -13 to d 0) and mRNA transcription was examined one day before gradual weaning was initiated (d -14), one day after weaning (d 1), and 8 days after weaning (d 8). On d -14, 550 genes were differentially expressed between Holstein-Friesian and Jersey calves, while there were 490 differentially expressed genes (DEG) identified on d 1, and 411 DEG detected eight days after weaning (P < 0.05; FDR < 0.1). No genes were differentially expressed within breed, in response to gradual weaning (P > 0.05). The pathways, gene ontology terms, and biological functions consistently over-represented among the DEG between Holstein-Friesian and Jersey were associated with the immune response and immune cell signalling, specifically chemotaxis. Decreased transcription of several cytokines, chemokines, immunoglobulin-like genes, phagocytosis-promoting receptors and g-protein coupled receptors suggests decreased monocyte, natural killer cell, and T lymphocyte, chemotaxis and activation in Jersey compared to Holstein-Friesian calves. Knowledge of breed-specific immune responses could facilitate health management practices better tailored towards specific disease sensitivities of Holstein-Friesian and Jersey calves. Gradual weaning did not compromise the welfare of artificially-reared dairy calves, evidenced by the lack of alterations in the expression of any genes in response to gradual weaning.

Developmental expression of STATs, nuclear factor-κB and inflammatory genes in the jejunum of piglets during weaning

The signal transducer and activator of transcription (STAT) proteins play essential roles in apoptosis, proliferation and survival. However, the role of STATs in intestinal inflammation during weaning is unclear. This study aimed to investigate developmental expression of STATs, nuclear factor-κB (NF-κB) and inflammatory genes in the jejunum of piglets during weaning. Thirty-two piglets were weaned at 21d and sacrificed at 0, 1, 7, or 14d (n=8) after weaning. Villus height and the villus height/crypt depth ratio were decreased, whereas crypt depth was increased in the jejunum at 7 and 14d after weaning. In addition, the mRNA levels of interferon-γ (IFN-γ), inducible nitric oxide synthase (iNOS), IL-6, IL-8, IL-12 and IL-22 were increased in the jejunum at 7 and 14d after weaning, whereas transforming growth factor-β (TGF-β), suppressor of cytokine signaling 3 (SCOS3) and arginase-1 was decreased. Neutrophil infiltration was increased in the mucosa of the jejunum after weaning. Moreover, phosphorylation of IκB-α, NF-κB, AKT and STAT-3 was increased. However, the phosphorylation of STAT-1 (at 7 and 14d) and STAT-6 (at 1 and 7d) was suppressed in the jejunum after weaning. Treatment of porcine jejunal epithelial (IPEC-J2) cells with the STAT inhibitors fludarabine, niclosamide and teriflunomide, which inhibit the phosphorylation of STAT-1, STAT-3 and STAT-6, respectively, weakened the defense capacity of these cells against bacterial infection. In conclusion, weaning caused severe inflammation associated with activation of the NF-κB and STAT-3 pathways and suppression of STAT-1 and STAT-6 in the jejunum of piglets.

Fertility and genomics: comparison of gene expression in contrasting reproductive tissues of female cattle

To compare gene expression among bovine tissues, large bovine RNA-seq datasets were used, comprising 280 samples from 10 different bovine tissues (uterine endometrium, granulosa cells, theca cells, cervix, embryos, leucocytes, liver, hypothalamus, pituitary, muscle) and generating 260 Gbases of data. Twin approaches were used: an information–theoretic analysis of the existing annotated transcriptome to identify the most tissue-specific genes and a de-novo transcriptome annotation to evaluate general features of the transcription landscape. Expression was detected for 97% of the Ensembl transcriptome with at least one read in one sample and between 28% and 66% at a level of 10 tags per million (TPM) or greater in individual tissues. Over 95% of genes exhibited some level of tissue-specific gene expression. This was mostly due to different levels of expression in different tissues rather than exclusive expression in a single tissue. Less than 1% of annotated genes exhibited a highly restricted tissue-specific expression profile and approximately 2% exhibited classic housekeeping profiles. In conclusion, it is the combined effects of the variable expression of large numbers of genes (73%–93% of the genome) and the specific expression of a small number of genes (<1% of the transcriptome) that contribute to determining the outcome of the function of individual tissues.

Evaluating the Stability of RNA-Seq Transcriptome Profiles and Drug-Induced Immune-Related Expression Changes in Whole Blood

Methods were developed to evaluate the stability of rat whole blood expression obtained from RNA sequencing (RNA-seq) and assess changes in whole blood transcriptome profiles in experiments replicated over time. Expression was measured in globin-depleted RNA extracted from the whole blood of Sprague-Dawley rats, given either saline (control) or neurotoxic doses of amphetamine (AMPH). The experiment was repeated four times (paired control and AMPH groups) over a 2-year span. The transcriptome of the control and AMPH-treated groups was evaluated on: 1) transcript levels for ribosomal protein subunits; 2) relative expression of immune-related genes; 3) stability of the control transcriptome over 2 years; and 4) stability of the effects of AMPH on immune-related genes over 2 years. All, except one, of the 70 genes that encode the 80s ribosome had levels that ranked in the top 5% of all mean expression levels. Deviations in sequencing performance led to significant changes in the ribosomal transcripts. The overall expression profile of immune-related genes and genes specific to monocytes, T-cells or B-cells were well represented and consistent within treatment groups. There were no differences between the levels of ribosomal transcripts in time-matched control and AMPH groups but significant differences in the expression of immune-related genes between control and AMPH groups. AMPH significantly increased expression of some genes related to monocytes but down-regulated those specific to T-cells. These changes were partially due to changes in the two types of leukocytes present in blood, which indicate an activation of the innate immune system by AMPH. Thus, the stability of RNA-seq whole blood transcriptome can be verified by assessing ribosomal protein subunits and immune-related gene expression. Such stability enables the pooling of samples from replicate experiments to carry out differential expression analysis with acceptable power.

Acute phase cytokines, TAC1, and toll-like receptor4 mRNA expression and health associated with group size in veal calves

Chronic stressors are a major health and well-being issue in animals. Immune status of animals under chronic stress is compromised, thus reducing disease resistance and compromising well-being of the animal. The objective of this study was to determine the influence of group size of veal calves on immune status and leukocyte mRNA expression of acute phase cytokines, toll-like receptor 4 (TLR4) and tachykinin 1 (TAC1) over a five-month finishing period. Holstein bull calves (n=168), 44±3 days of age were assigned to one of three treatments; 2, 4, or 8 calves/pen (pen space allowance of 1.82m(2)/calf). Jugular blood samples were collected at the day of grouping and then monthly for 4 months. The differential leukocyte counts were determined and mRNA was extracted from the leukocytes. Reverse transcription-qPCR was used to measure the gene expression of interleukin-1 (IL-1β), IL-1 receptor antagonist (IL-1Ra), tumor necrosis factor (TNF-α), TLR4, and TAC1 in leukocytes. Health was evaluated before grouping and monthly for 4 months. On the 1st month after grouping, veal calves that were housed in groups of 8 have greater expression of IL-1β mRNA than calves housed in groups of 4 or 2 (treatment×month, P=0.04). Also at 1 month, groups of 8 had greater TAC1 expression (P<0.05) than calves housed in groups of 4 or 2. However, the expression of IL-1Ra, TNF-α, and TLR4 were not influenced by group size. In the first month of the trial, calves in groups of 8 coughed more (P<0.05) than calves in groups of 2 and coughed more than calves in groups of 4 and 2 during the 2nd month (treatment×month, P=0.03). Calves housed in groups of 8 tended to have greater neutrophil percentage (P=0.09), neutrophil to lymphocyte ratio (P=0.06), and had lower lymphocyte percentage (P=0.06) than those housed in groups of 4 or 2. In conclusion, the number of veal calves in a group, given the same space during the finishing period did not alter IL-1Ra, TNF-α, and TLR4 mRNA expression. However, housing of calves in groups of 8 was associated with greater expression of IL-1β and TAC1 mRNA in peripheral blood leukocytes, and coughing during the first 2 months after grouping. Therefore, housing of veal calves in larger groups may lead to greater susceptibility to respiratory disease and stress.

Variation in RNA-Seq transcriptome profiles of peripheral whole blood from healthy individuals with and without globin depletion

Background

The molecular profile of circulating blood can reflect physiological and pathological events occurring in other tissues and organs of the body and delivers a comprehensive view of the status of the immune system. Blood has been useful in studying the pathobiology of many diseases. It is accessible and easily collected making it ideally suited to the development of diagnostic biomarker tests. The blood transcriptome has a high complement of globin RNA that could potentially saturate next-generation sequencing platforms, masking lower abundance transcripts. Methods to deplete globin mRNA are available, but their effect has not been comprehensively studied in peripheral whole blood RNA-Seq data. In this study we aimed to assess technical variability associated with globin depletion in addition to assessing general technical variability in RNA-Seq from whole blood derived samples.

Results

We compared technical and biological replicates having undergone globin depletion or not and found that the experimental globin depletion protocol employed removed approximately 80% of globin transcripts, improved the correlation of technical replicates, allowed for reliable detection of thousands of additional transcripts and generally increased transcript abundance measures. Differential expression analysis revealed thousands of genes significantly up-regulated as a result of globin depletion. In addition, globin depletion resulted in the down-regulation of genes involved in both iron and zinc metal ion bonding.

Conclusions

Globin depletion appears to meaningfully improve the quality of peripheral whole blood RNA-Seq data, and may improve our ability to detect true biological variation. Some concerns remain, however. Key amongst them the significant reduction in RNA yields following globin depletion. More generally, our investigation of technical and biological variation with and without globin depletion finds that high-throughput sequencing by RNA-Seq is highly reproducible within a large dynamic range of detection and provides an accurate estimation of RNA concentration in peripheral whole blood. High-throughput sequencing is thus a promising technology for whole blood transcriptomics and biomarker discovery.

Effect of dietary n-3 polyunsaturated fatty acids on transcription factor regulation in the bovine endometrium

Dietary n-3 polyunsaturated fatty acid (n-3 PUFA) supplementation is postulated to have positive effects on fertility. The impact of dietary n-3 PUFA supplementation on physiological and biochemical processes involved in reproduction is likely to be associated with significant alterations in gene expression in key reproductive tissues which is in turn regulated by transcription factors. Beef heifers were supplemented with a rumen protected source of either a saturated fatty acid or high n-3 PUFA diet per animal per day for 45 days and uterine endometrial tissue was harvested post slaughter. A microarray analysis was conducted and bioinformatic tools were employed to evaluate the effect of n-3 PUFA supplementation on gene expression in the bovine endometrium. Clustering of microarray gene expression data was performed to identify co-expressed genes. Functional annotation of each cluster of genes was carried out using Ingenuity Pathway Analysis. Furthermore, oPOSSUM was employed to identify transcription factors involved in gene expression changes due to supplementary PUFA. Gene functions which showed a significant response to n-3 PUFA supplementation included tissue development, immune function and reproductive function. Numerous transcription factors such as FOXD1, FOXD3, NFKB1, ESR1, PGR, FOXA2, NKX3-1 and PPARα were identified as potential regulators of gene expression in the endometrium of cattle supplemented with n-3 PUFA. This study demonstrates the complex nature of the alterations in the transcriptional regulation process in the uterine endometrium of cattle following dietary supplementation which may positively influence the uterine environment.

Identification of novel transcripts and noncoding RNAs in bovine skin by deep next generation sequencing

BACKGROUND

Deep RNA sequencing (RNAseq) has opened a new horizon for understanding global gene expression. The functional annotation of non-model mammalian genomes including bovines is still poor compared to that of human and mouse. This particularly applies to tissues without direct significance for milk and meat production, like skin, in spite of its multifunctional relevance for the individual. Thus, applying an RNAseq approach, we performed a whole transcriptome analysis of pigmented and nonpigmented bovine skin to describe the comprehensive transcript catalogue of this tissue.

RESULTS

A total of 39,577 unique primary skin transcripts were mapped to the bovine reference genome assembly. The majority of the transcripts were mapped to known transcriptional units (65%). In addition to the reannotation of known genes, a substantial number (10,884) of unknown transcripts (UTs) were discovered, which had not previously been annotated. The classification of UTs was based on the prediction of their coding potential and comparative sequence analysis, subsequently followed by meticulous manual curation. The classification analysis and experimental validation of selected UTs confirmed that RNAseq data can be used to amend the annotation of known genes by providing evidence for additional exons, untranslated regions or splice variants, by approving genes predicted in silico and by identifying novel bovine loci. A large group of UTs (4,848) was predicted to potentially represent long noncoding RNA (lncRNA). Predominantly, potential lncRNAs mapped in intergenic chromosome regions (4,365) and therefore, were classified as potential intergenic lncRNA. Our analysis revealed that only about 6% of all UTs displayed interspecies conservation and discovered a variety of unknown transcripts without interspecies homology but specific expression in bovine skin.

CONCLUSIONS

The results of our study demonstrate a complex transcript pattern for bovine skin and suggest a possible functional relevance of novel transcripts, including lncRNA, in the modulation of pigmentation processes. The results also indicate that the comprehensive identification and annotation of unknown transcripts from whole transcriptome analysis using RNAseq data remains a tremendous future challenge.

Monitoring the immune response to vaccination with an inactivated vaccine associated to bovine neonatal pancytopenia by deep sequencing transcriptome analysis in cattle

Bovine neonatal pancytopenia (BNP) is a new fatal, alloimmune/alloantibody mediated disease of new-born calves induced by ingestion of colostrum from cows, which had been vaccinated with a specific vaccine against the Bovine Virus Diarrhoea Virus (BVDV). The hypothesis of pathogenic MHC class I molecules in the vaccine had been put up, but no formal proof of specific causal MHC class I alleles has been provided yet. However, the unique features of the vaccine obviously result in extremely high specific antibody titres in the vaccinated animals, but apparently also in further molecules inducing BNP. Thus, a comprehensive picture of the immune response to the vaccine is essential. Applying the novel approach of next generation RNA sequencing (RNAseq), our study provides a new holistic, comprehensive analysis of the blood transcriptome regulation after vaccination with the specific BVDV vaccine. Our RNAseq approach identified a novel cytokine-like gene in the bovine genome that is highly upregulated after vaccination. This gene has never been described before in any other species and might be specific to ruminant immune response. Furthermore, our data revealed a very coordinated immune response to double-stranded (ds) RNA or a dsRNA analogue after vaccination with the inactivated single-stranded (ss) RNA vaccine. This would suggest either a substantial contamination of the vaccine with dsRNA from host cells after virus culture or a dsRNA analogue applied to the vaccine. The first option would highlight the potential risks associated with virus culture on homologous cells during vaccine production; the latter option would emphasise the potential risks associated with immune stimulating adjuvants used in vaccine production.

Failure of respiratory defenses in the pathogenesis of bacterial pneumonia of cattle

The respiratory system is well defended against inhaled bacteria by a dynamic system of interacting layers, including mucociliary clearance, host defense factors including antimicrobial peptides in the epithelial lining fluid, proinflammatory responses of the respiratory epithelium, resident alveolar macrophages, and recruited neutrophils and monocytes. Nevertheless, these manifold defenses are susceptible to failure as a result of stress, glucocorticoids, viral infections, abrupt exposure to cold air, and poor air quality. When some of these defenses fail, the lung can be colonized by bacterial pathogens that are equipped to evade the remaining defenses, resulting in the development of pneumonia. This review considers the mechanisms by which these predisposing factors compromise the defenses of the lung, with a focus on the development of bacterial pneumonia in cattle and supplemented with advances based on mouse models and the study of human disease. Deepening our understanding of how the respiratory defenses fail is expected to lead to interventions that restore these dynamic immune responses and prevent disease.

Transcriptomic signatures in cartilage ageing

Introduction

Age is an important factor in the development of osteoarthritis. Microarray studies provide insight into cartilage aging but do not reveal the full transcriptomic phenotype of chondrocytes such as small noncoding RNAs, pseudogenes, and microRNAs. RNA-Seq is a powerful technique for the interrogation of large numbers of transcripts including nonprotein coding RNAs. The aim of the study was to characterise molecular mechanisms associated with age-related changes in gene signatures.

Methods

RNA for gene expression analysis using RNA-Seq and real-time PCR analysis was isolated from macroscopically normal cartilage of the metacarpophalangeal joints of eight horses; four young donors (4 years old) and four old donors (>15 years old). RNA sequence libraries were prepared following ribosomal RNA depletion and sequencing was undertaken using the Illumina HiSeq 2000 platform. Differentially expressed genes were defined using Benjamini-Hochberg false discovery rate correction with a generalised linear model likelihood ratio test (P < 0.05, expression ratios ± 1.4 log₂ fold-change). Ingenuity pathway analysis enabled networks, functional analyses and canonical pathways from differentially expressed genes to be determined.

Results

In total, the expression of 396 transcribed elements including mRNAs, small noncoding RNAs, pseudogenes, and a single microRNA was significantly different in old compared with young cartilage (± 1.4 log₂ fold-change, P < 0.05). Of these, 93 were at higher levels in the older cartilage and 303 were at lower levels in the older cartilage. There was an over-representation of genes with reduced expression relating to extracellular matrix, degradative proteases, matrix synthetic enzymes, cytokines and growth factors in cartilage derived from older donors compared with young donors. In addition, there was a reduction in Wnt signalling in ageing cartilage.

Conclusion

There was an age-related dysregulation of matrix, anabolic and catabolic cartilage factors. This study has increased our knowledge of transcriptional networks in cartilage ageing by providing a global view of the transcriptome.

Challenges in drug target discovery in bipolar disorder

INTRODUCTION

Misdiagnosis and subsequent inappropriate treatment of patients with bipolar disorder (BD) can worsen their clinical condition and outcome.

AREAS COVERED

This review focuses on the therapeutic targets which have been implicated in BD, including the glycogen synthase kinase 3 (GSK-3) and phosphoinositide signaling pathways. In addition, evidence is presented for potential new molecular strategies which involve targeting neuropeptide-converting endopeptidases, glutamatergic excitotoxicity, insulin signaling and dysfunctions in mitochondrial metabolism. Current limitations in study design, molecular platforms, preclinical and cellular models in the context of BD drug target discovery, suggest that there are many areas for improvement.

EXPERT OPINION

For the future outlook, this review outlines the importance of developments such as the use of BD patient-derived cellular models for providing better understanding of the BD etiology and robust translational drug screening tools in combination with developments in the fields of bioinformatics and systems biology.