Evaluation of short-term hair follicle storage conditions for maintenance of RNA integrity.

Hair follicles provide an easily accessible tissue for interrogating gene expression for multiple purposes in mammals. RNAlater® is a liquid storage solution that stabilises and preserves cellular RNA, eliminating the need to immediately process or freeze tissue specimens. The manufacturer advises storage of samples at 2-8°C overnight before transfer to -20°C. This study aimed to evaluate RNA integrity in hair follicle samples collected from horses, stabilized in RNAlater®, and stored under three short-term storage conditions. Mane hair samples complete with follicles were collected from four horses at a single time point. Approximately 15 hairs were placed in each of three 2 mL tubes containing 0.75ml RNAlater® solution. Test group A was stored at 4°C for 24-h, then decanted and stored at -20°C. Test groups B and C were stored at 4°C and 19°C (room temperature) respectively for 7 days, then decanted and stored at -20°C. RNA was isolated from all samples and RNA quantity and quality were measured. One-way ANOVA revealed no difference in RNA concentration (A:516 +/-125 ng/ml, B:273+/-93 ng/ml, C:476+/-176 ng/ml;P = 0.2) or quality (A:9.5 +/-0.19, B:9.8+/-0.09, C:9.2+/-0.35 RIN; P = 0.46) between the test groups. There were no group differences in mean Cycle Threshold values from qPCR validation assays confirming high-quality template cDNA. The results suggest that storage of hair follicles for one week in RNAlater® at cool or room temperature conditions will not compromise RNA integrity and will permit extended transport times from remote sampling locations without the need for freezing.

Genomic insights into the population history and adaptive traits of Latin American Criollo cattle.

Criollo cattle, the descendants of animals brought by Iberian colonists to the Americas, have been the subject of natural and human-mediated selection in novel tropical agroecological zones for centuries. Consequently, these breeds have evolved distinct characteristics such as resistance to diseases and exceptional heat tolerance. In addition to European taurine (Bos taurus) ancestry, it has been proposed that gene flow from African taurine and Asian indicine (Bos indicus) cattle has shaped the ancestry of Criollo cattle. In this study, we analysed Criollo breeds from Colombia and Venezuela using whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) array data to examine population structure and admixture at high resolution. Analysis of genetic structure and ancestry components provided evidence for African taurine and Asian indicine admixture in Criollo cattle. In addition, using WGS data, we detected selection signatures associated with a myriad of adaptive traits, revealing genes linked to thermotolerance, reproduction, fertility, immunity and distinct coat and skin coloration traits. This study underscores the remarkable adaptability of Criollo cattle and highlights the genetic richness and potential of these breeds in the face of climate change, habitat flux and disease challenges. Further research is warranted to leverage these findings for more effective and sustainable cattle breeding programmes.

Transcriptomics analysis of the bovine endometrium during the perioestrus period.

During the oestrous cycle, the bovine endometrium undergoes morphological and functional changes, which are regulated by alterations in the levels of oestrogen and progesterone and consequent changes in gene expression. To clarify these changes before and after oestrus, RNA-seq was used to profile the transcriptome of oestrus-synchronized beef heifers. Endometrial samples were collected from 29 animals, which were slaughtered in six groups beginning 12 h after the withdrawal of intravaginal progesterone releasing devices until seven days post-oestrus onset (luteal phase). The groups represented proestrus, early oestrus, metoestrus and early dioestrus (luteal phase). Changes in gene expression were estimated relative to gene expression at oestrus. Ingenuity Pathway Analysis (IPA) was used to identify canonical pathways and functional processes of biological importance. A total of 5,845 differentially expressed genes (DEGs) were identified. The lowest number of DEGs was observed at the 12 h post-oestrus time point, whereas the greatest number was observed at Day 7 post-oestrus onset (luteal phase). A total of 2,748 DEGs at this time point did not overlap with any other time points. Prior to oestrus, Neurological disease and Organismal injury and abnormalities appeared among the top IPA diseases and functions categories, with upregulation of genes involved in neurogenesis. Lipid metabolism was upregulated before oestrus and downregulated at 48h post-oestrus, at which point an upregulation of immune-related pathways was observed. In contrast, in the luteal phase the Lipid metabolism and Small molecule biochemistry pathways were upregulated.

Integrative and comparative genomic analyses of mammalian macrophage responses to intracellular mycobacterial pathogens.

Mycobacterium tuberculosis, the causative agent of human tuberculosis (hTB), is a close evolutionary relative of Mycobacterium bovis, which causes bovine tuberculosis (bTB), one of the most damaging infectious diseases to livestock agriculture. Previous studies have shown that the pathogenesis of bTB disease is comparable to hTB disease, and that the bovine and human alveolar macrophage (bAM and hAM, respectively) transcriptomes are extensively reprogrammed in response to infection with these intracellular mycobacterial pathogens. In this study, a multi-omics integrative approach was applied with functional genomics and GWAS data sets across the two primary hosts (Bos taurus and Homo sapiens) and both pathogens (M. bovis and M. tuberculosis). Four different experimental infection groups were used: 1) bAM infected with M. bovis, 2) bAM infected with M. tuberculosis, 3) hAM infected with M. tuberculosis, and 4) human monocyte-derived macrophages (hMDM) infected with M. tuberculosis. RNA-seq data from these experiments 24 h post-infection (24 hpi) was analysed using three computational pipelines: 1) differentially expressed genes, 2) differential gene expression interaction networks, and 3) combined pathway analysis. The results were integrated with high-resolution bovine and human GWAS data sets to detect novel quantitative trait loci (QTLs) for resistance to mycobacterial infection and resilience to disease. This revealed common and unique response macrophage pathways for both pathogens and identified 32 genes (12 bovine and 20 human) significantly enriched for SNPs associated with disease resistance, the majority of which encode key components of the NF-κB signalling pathway and that also drive formation of the granuloma.

Optimised Stable Lighting Strengthens Circadian Clock Gene Rhythmicity in Equine Hair Follicles.

Hair follicles (HF) represent a useful tissue for monitoring the circadian clock in mammals. Irregular light exposure causes circadian disruption and represents a welfare concern for stabled horses. We aimed to evaluate the impact of two stable lighting regimes on circadian clock gene rhythmicity in HF from racehorses. Two groups of five Thoroughbred racehorses in training at a commercial racehorse yard were exposed to standard incandescent light or a customized LED lighting system. The control group received light from incandescent bulbs used according to standard yard practice. The treatment group received timed, blue-enriched white LED light by day and dim red LED light at night. On weeks 0 and 20, mane hairs were collected at 4 h intervals for 24 h. Samples were stored in RNAlater at -20 °C. RNA was isolated and samples interrogated by quantitative PCR for the core clock genes: ARNTL, CRY1, PER1, PER2, NR1D2, and the clock-controlled gene DBP. Cosinor analyses revealed 24 h rhythmicity for NR1D2 and PER2 and approached significance for CRY1 (p = 0.013, p = 0.013, and p = 0.051, respectively) in week 20 in the treatment group only. No rhythmicity was detected in week 0 or in week 20 in the HF of control horses. Results suggest that lighting practices in racehorse stables may be improved to better stimulate optimum functioning of the circadian system.

The comparative performance of a custom Canine NanoString® panel on FFPE and snap frozen liver biopsies.

Formalin-Fixed Paraffin Embedded (FFPE) biopsies would provide a critical mass of cases to allow investigation of canine liver disease, however their use is often limited by challenges typically associated with transcriptomic analysis. This study evaluates the capability of NanoString® to measure the expression of a broad panel of genes in FFPE liver samples. RNA was isolated from matched histopathologically normal liver samples using FFPE (n = 6) and snap frozen in liquid nitrogen (n = 6) and measured using a custom NanoString® panel. Out of the 40 targets on the panel, 27 and 23 targets were above threshold for non-diseased snap frozen and FFPE tissue respectively. The binding density and total counts were significantly reduced in the FFPE samples relative to the snap frozen samples (p = 0.005, p = 0.01, respectively), confirming a reduction in sensitivity. The concordance between the snap frozen and FFPE samples was high, with correlations (R) ranging between 0.88 and 0.99 between the paired samples. An additional 14 immune-related targets, undetectable the non-diseased FFPE liver, were above threshold when the technique was applied to a series of diseased samples, further supporting their inclusion on this panel. This use of NanoString® based analysis opens up huge opportunity for retrospective evaluation of gene signatures in larger caseloads through harnessing the capacity of archived FFPE samples This information used alongside clinical and histological data will not only afford a way to explore disease etiopathogenesis, it may also offer insight into sub-types of liver disease in dogs, which cannot be discerned using more traditional diagnostic methods.

Variation in salivary cortisol responses in yearling Thoroughbred racehorses during their first year of training.

Thoroughbred horses are bred for competitive racing and undergo intense training regimes. The maintenance of physical soundness and desirable behavioural characteristics are critical to the longevity of a racing career. Horses intended for Flat racing generally enter training as yearlings and undergo introductory training prior to exercise conditioning for racing. This period requires rapid adjustment to a novel environment. As a prey animal, a horse's 'fight-or-flight' response is highly adapted, in which a well-understood component of this response, the hypothalamic-pituitary-axis, is activated in response to a stress stimulus, releasing cortisol. In the Thoroughbred, a significant difference in salivary cortisol concentrations between pre- and post-first time ridden (i.e., first backing) by a jockey have previously been identified. Here, to test the hypothesis that salivary cortisol concentrations may be used to objectively detect individual variations in the acute physiological stress response we investigate individual variation in cortisol response to training milestones. Saliva samples were collected from a cohort of n = 96 yearling Flat racehorses, at the same training yard, across three timepoints at rest: before entering the training yard (n = 66), within three days of entry to the training yard (n = 67) and following 2-3 weeks in the training yard (n = 50). Salivary cortisol concentration was measured using an ELISA. There was no significant difference in cortisol concentration (ANOVA, P > 0.05) across the samples collected at timepoints at rest. Samples were also collected before and 30 minutes after exposure to three novel training events: first time long-reined (n = 6), first time backed by a jockey (n = 34), and first time ridden on the gallops (n = 10). Mean salivary cortisol concentration after all three novel training events was significantly higher than prior to the training event (Paired t-test, P <0.005). The ranges of post-event salivary cortisol concentration across all timepoints suggest individual variation in the measured stress response, reflecting individual differences in stress response to the early training environment. This measure may be used as an objective assessment of the stress response of Thoroughbred racehorses during training.

High-resolution transcriptomics of bovine purified protein derivative-stimulated peripheral blood from cattle infected with Mycobacterium bovis across an experimental time course.

OBJECTIVES: Improved bovine tuberculosis (bTB) diagnostics with higher sensitivity and specificity are urgently required. A better understanding of the peripheral blood transcriptional response of Mycobacterium bovis-infected animals after bovine purified protein derivative (PPD-b) stimulation of whole blood-an important component of current bTB diagnostics-will provide new information for development of better diagnostics. METHODS: RNA sequencing (RNA-seq) was used to study the peripheral blood transcriptome after stimulation with PPD-b across four time points (-1 wk pre-infection, and +1 wk, +2 wk, and +10 wk post-infection) from a 14-week M. bovis infection time course experiment with ten age-matched Holstein-Friesian cattle. RESULTS: In vitro PPD-b stimulation of peripheral blood from M. bovis-infected and non-infected cattle elicited a strong transcriptional response. Comparison of PPD-b stimulated, and unstimulated samples revealed higher expression of genes encoding cytokine receptors, transcription factors, and interferon-inducible proteins. Lower expression was seen for genes encoding proteins involved in antimicrobial activity, C-type lectin receptors, inhibition of signal transduction, and genes encoding metal ion transporters. CONCLUSIONS: A transcriptional signature associated with the peripheral blood response to PPD-b stimulation consisting of 170 genes was identified exclusively in the post-infection time points. Therefore, this represents a panel of potential biomarkers of M. bovis infection.

Corrigendum: The TbD1 Locus Mediates a Hypoxia-Induced Copper Response in Mycobacterium bovis.

[This corrects the article DOI: 10.3389/fmicb.2022.817952.].

The TbD1 Locus Mediates a Hypoxia-Induced Copper Response in Mycobacterium bovis.

The Mycobacterium tuberculosis complex (MTBC) contains the causative agents of tuberculosis (TB) in mammals. The archetypal members of the MTBC, Mycobacterium tuberculosis and Mycobacterium bovis, cause human tuberculosis and bovine tuberculosis, respectively. Although M. tuberculosis and M. bovis share over 99.9% genome identity, they show distinct host adaptation for humans and animals; hence, while the molecular basis of host adaptation is encoded in their genomes, the mechanistic basis of host tropism is still unclear. Exploration of the in vitro phenotypic consequences of known genetic difference between M. bovis and M. tuberculosis offers one route to explore genotype-phenotype links that may play a role in host adaptation. The TbD1 ("Mycobacterium tuberculosis deletion 1 region") locus encompasses the mmpS6 and mmpL6 genes. TbD1 is absent in M. tuberculosis "modern" lineages (Lineages 2, 3, and 4) but present in "ancestral" M. tuberculosis (Lineages 1 and 7), Mycobacterium africanum lineages (Lineages 5 and 6), newly identified M. tuberculosis lineages (Lineages 8 and 9), and animal adapted strains, such as M. bovis. The function of TbD1 has previously been investigated in M. tuberculosis, where conflicting data has emerged on the role of TbD1 in sensitivity to oxidative stress, while the underlying mechanistic basis of such a phenotype is unclear. In this study, we aimed to shed further light on the role of the TbD1 locus by exploring its function in M. bovis. Toward this, we constructed an M. bovis TbD1 knockout (ΔTbD1) strain and conducted comparative transcriptomics to define global gene expression profiles of M. bovis wild-type (WT) and the ΔTbD1 strains under in vitro culture conditions (rolling and standing cultures). This analysis revealed differential induction of a hypoxia-driven copper response in WT and ΔTbD1 strains. In vitro phenotypic assays demonstrated that the deletion of TbD1 sensitized M. bovis to H2O2 and hypoxia-specific copper toxicity. Our study provides new information on the function of the TbD1 locus in M. bovis and its role in stress responses in the MTBC.

Evidence for local and international spread of Mycobacterium avium subspecies paratuberculosis through whole genome sequencing of isolates from the island of Ireland.

We describe application of whole genome sequencing (WGS) to a collection of 197 Mycobacterium avium subsp paratuberculosis (MAP) isolates gathered from 122 cattle herds across 27 counties of the island of Ireland. We compare WGS to MAP diversity quantified using mycobacterial interspersed random unit - variable number tandem repeats (MIRU-VNTR). While MIRU-VNTR showed only two major types, WGS could split the 197 isolates into eight major groups. We also found six isolates corresponding to INMV 13, a novel MIRU-VNTR type for Ireland. Evidence for dispersal of MAP across Ireland via cattle movement could be discerned from the data, with mixed infections present in several herds. Furthermore, comparisons of MAP WGS data from Ireland to data from Great Britain and continental Europe revealed many instances of close genetic similarity and hence evidence for international transmission of infection. BEAST MASCOT structured coalescent analyses, with relaxed and strict molecular clocks, estimated the substitution rate to be 0.10-0.13 SNPs/site/year and disclosed greater transitions per lineage per year from Europe to Ireland, indicating transmission into Ireland. Our work therefore reveals new insight into the seeding of MAP infection across Ireland, highlighting how WGS can inform policy formulation to ultimately control MAP transmission at local, national and international scales.

Timing of Transcriptomic Peripheral Blood Mononuclear Cell Responses of Sheep to Fasciola hepatica Infection Differs From Those of Cattle, Reflecting Different Disease Phenotypes.

Infection with the zoonotic trematode Fasciola hepatica, common in many regions with a temperate climate, leads to delayed growth and loss of productivity in cattle, while infection in sheep can have more severe effects, potentially leading to death. Previous transcriptomic analyses revealed upregulation of TGFB1, cell death and Toll-like receptor signalling, T-cell activation, and inhibition of nitric oxide production in macrophages in response to infection. However, the differences between ovine and bovine responses have not yet been explored. The objective of this study was to further investigate the transcriptomic response of ovine peripheral blood mononuclear cells (PBMC) to F. hepatica infection, and to elucidate the differences between ovine and bovine PBMC responses. Sixteen male Merino sheep were randomly assigned to infected or control groups (n = 8 per group) and orally infected with 120 F. hepatica metacercariae. Transcriptomic data was generated from PBMC at 0, 2 and 16 weeks post-infection (wpi), and analysed for differentially expressed (DE) genes between infected and control animals at each time point (analysis 1), and for each group relative to time 0 (analysis 2). Analysis 2 was then compared to a similar study performed previously on bovine PBMC. A total of 453 DE genes were found at 2 wpi, and 2 DE genes at 16 wpi (FDR < 0.1, analysis 1). Significantly overrepresented biological pathways at 2 wpi included role of PKR in interferon induction and anti-viral response, death receptor signalling and RIG-I-like receptor signalling, which suggested that an activation of innate response to intracellular nucleic acids and inhibition of cellular apoptosis were taking place. Comparison of analysis 2 with the previous bovine transcriptomic study revealed that anti-inflammatory response pathways which were significantly overrepresented in the acute phase in cattle, including IL-10 signalling, Th2 pathway, and Th1 and Th2 activation were upregulated only in the chronic phase in sheep. We propose that the earlier activation of anti-inflammatory responses in cattle, as compared with sheep, may be related to the general absence of acute clinical signs in cattle. These findings offer scope for "smart vaccination" strategies for this important livestock parasite.

Insights into the ancestry evolution of the Mycobacterium tuberculosis complex from analysis of Mycobacterium riyadhense.

Current evolutionary scenarios posit the emergence of Mycobacterium tuberculosis from an environmental saprophyte through a cumulative process of genome adaptation. Mycobacterium riyadhense, a related bacillus, is being increasingly isolated from human clinical cases with tuberculosis-like symptoms in various parts of the world. To elucidate the evolutionary relationship between M. riyadhense and other mycobacterial species, including members of the M. tuberculosis complex (MTBC), eight clinical isolates of M. riyadhense were sequenced and analyzed. We show, among other features, that M. riyadhense shares a large number of conserved orthologs with M. tuberculosis and shows the expansion of toxin/antitoxin pairs, PE/PPE family proteins compared with other non-tuberculous mycobacteria. We observed M. riyadhense lacks wecE gene which may result in the absence of lipooligosaccharides (LOS) IV. Comparative transcriptomic analysis of infected macrophages reveals genes encoding inducers of Type I IFN responses, such as cytosolic DNA sensors, were relatively less expressed by macrophages infected with M. riyadhense or M. kansasii, compared to BCG or M. tuberculosis. Overall, our work sheds new light on the evolution of M. riyadhense, its relationship to the MTBC, and its potential as a system for the study of mycobacterial virulence and pathogenesis.

Immunological aspects of ovarian follicle ovulation and corpus luteum formation in cattle.

Ovulation has been described as an inflammatory event, characterized by an influx of leukocytes into the ovulatory follicle and changes in the expression profile of immune factors in both the theca and granulosa tissue layers. Since information on this process is limited in cattle, our objective was to elucidate the contribution of the immune system to dominant follicle luteinization, ovulation and corpus luteum (CL) formation in cattle. Beef heifers (n = 50) were oestrous synchronized, slaughtered and ovarian follicular or luteal tissue collected during a 96 h window around ovulation. Follicular fluid cytokine concentration, temporal immune cell infiltration and inflammatory status were determined by Luminex multiplex analysis, immunohistochemistry and quantitative real-time PCR-analysis, respectively, in pre- and peri-ovulatory follicular tissues. The concentrations of IL10 and VEGF-A were highest in pre-ovulatory and the concentration of CXCL10 was highest in peri-ovulatory follicular fluid samples. The pre and peri-ovulatory follicles play host to a broad repertoire of immune cells, including T-cells, granulocytes and monocytes. Dendritic cells were the most abundant cells in ovulatory follicular and luteal-tissue at all times. The mRNA expression of candidate genes associated with inflammation was highest in pre- and peri-ovulatory tissue, whereas tissue growth and modelling factors were highest in the post-ovulatory follicular and early luteal tissue. In conclusion, ovulation in cattle is characterized by the presence of neutrophils, macrophages and dendritic cells in the ovulatory follicle, reflected in compartmentalized cytokine and growth factor expression. These findings indicate a tightly regulated sterile inflammatory response to the LH surge in the ovulatory follicle which is rapidly resolved in advance of CL formation.

Mycobacterial Infection of Precision-Cut Lung Slices Reveals Type 1 Interferon Pathway Is Locally Induced by Mycobacterium bovis but Not M. tuberculosis in a Cattle Breed.

Tuberculosis exacts a terrible toll on human and animal health. While Mycobacterium tuberculosis (Mtb) is restricted to humans, Mycobacterium bovis (Mb) is present in a large range of mammalian hosts. In cattle, bovine TB (bTB) is a noticeable disease responsible for important economic losses in developed countries and underestimated zoonosis in the developing world. Early interactions that take place between mycobacteria and the lung tissue early after aerosol infection govern the outcome of the disease. In cattle, these early steps remain poorly characterized. The precision-cut lung slice (PCLS) model preserves the structure and cell diversity of the lung. We developed this model in cattle in order to study the early lung response to mycobacterial infection. In situ imaging of PCLS infected with fluorescent Mb revealed bacilli in the alveolar compartment, in adjacent or inside alveolar macrophages, and in close contact with pneumocytes. We analyzed the global transcriptional lung inflammation signature following infection of PCLS with Mb and Mtb in two French beef breeds: Blonde d'Aquitaine and Charolaise. Whereas, lungs from the Blonde d'Aquitaine produced high levels of mediators of neutrophil and monocyte recruitment in response to infection, such signatures were not observed in the Charolaise in our study. In the Blonde d'Aquitaine lung, whereas the inflammatory response was highly induced by two Mb strains, AF2122 isolated from cattle in the UK and Mb3601 circulating in France, the response against two Mtb strains, H37Rv, the reference laboratory strain, and BTB1558, isolated from zebu in Ethiopia, was very low. Strikingly, the type I interferon pathway was only induced by Mb but not Mtb strains, indicating that this pathway may be involved in mycobacterial virulence and host tropism. Hence, the PCLS model in cattle is a valuable tool to deepen our understanding of early interactions between lung host cells and mycobacteria. It revealed striking differences between cattle breeds and mycobacterial strains. This model could help in deciphering biomarkers of resistance vs. susceptibility to bTB in cattle as such information is still critically needed for bovine genetic selection programs and would greatly help the global effort to eradicate bTB.

Transcriptomic Analysis of Ovine Hepatic Lymph Node Following Fasciola hepatica Infection - Inhibition of NK Cell and IgE-Mediated Signaling.

Fasciola hepatica is a trematode parasite responsible for major economic losses in livestock production, and is also a food-borne zoonotic agent in developing rural regions. For years, the immunoregulatory mechanisms employed by the parasite have hampered efforts to develop a successful vaccine candidate. Given that a comprehensive understanding of the immune response to infection is needed, we investigated the gene expression changes in ovine hepatic lymph nodes after experimental infection with F. hepatica. Lymph nodes from uninfected and infected animals were processed for RNA sequencing (RNA-seq) at 16 weeks post-infection. Comparison of groups revealed 5,132 differentially-expressed genes (DEGs). An inhibition of pro-inflammatory pathways, which has previously been described during fasciolosis, was evident in our data. However, other signals previously identified in ruminant peripheral blood mononuclear cells (PBMC) or liver tissue, such as activation of TGF-ß or apoptosis-related pathways were not detected. We found inhibition of some key immunological pathways, including natural killer (NK) cell activity and IgE-mediated signaling. These may point to additional some as yet unrecognized mechanisms employed by the parasite to evade the host immune response. Understanding these, and leveraging information from this and other omics studies, will be important for the development of future vaccine prototypes against this parasite.

RNA-Seq Transcriptome Analysis of Peripheral Blood From Cattle Infected With Mycobacterium bovis Across an Experimental Time Course.

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including RNA sequencing, has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analysed the transcriptome of bovine whole peripheral blood samples collected at -1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the -1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-infection.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria.

BACKGROUND: Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection. RESULTS: The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool-gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively. CONCLUSIONS: The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

MicroRNAs in amniotic fluid and maternal blood plasma associated with sex determination and early gonad differentiation in cattle†.

We hypothesized that sexually dimorphic differences exist in the expression of miRNAs in amniotic fluid (AF) and maternal blood plasma (MP) in association with the process of sex determination and gonad differentiation in cattle. Amniotic fluid and MP were collected from six pregnant heifers (three carrying a single male and three a single female embryo) following slaughter on Day 39 postinsemination, coinciding with the peak of SRY expression. Samples (six AF and six MP) were profiled using an miRNA Serum/Plasma Focus PCR Panel. Differentially expressed (DE) miRNAs were identified in AF (n = 5) and associated MP (n = 56) of male vs. female embryos (P < 0.05). Functional analysis showed that inflammatory and immune response were among the 13 biological processes enriched by miRNAs DE in MP in the male group (FDR < 0.05), suggesting that these sex-dependent DE miRNAs may be implicated in modulating the receptivity of the dam to a male embryo. Further, we compared the downstream targets of the sex-dependent DE miRNAs detected in MP with genes previously identified as DE in male vs. female genital ridges. The analyses revealed potential targets that might be important during this developmental stage such as SHROOM2, DDX3Y, SOX9, SRY, PPP1CB, JARID2, USP9X, KDM6A, and EIF2S3. Results from this study highlight novel aspects of sex determination and embryo-maternal communication in cattle such as the potential role of miRNAs in gonad development as well as in the modulation of the receptivity of the dam to a male embryo.

Plasma extracellular vesicle miRNAs as potential biomarkers of superstimulatory response in cattle.

The ability to predict superstimulatory response would be a beneficial tool in assisted reproduction. Using small RNAseq technology, we profiled extracellular vesicle microRNA (EV-miRNA) abundance in the blood plasma of heifers exhibiting variable responses to superstimulation. Estrous synchronized crossbred beef heifers (n = 25) were superstimulated and blood samples were collected from each heifer on Day 7 of consecutive unstimulated (U) and superstimulated (S) cycles. A subset of high (H) and low (L) responders was selected depending on their response to superstimulation and EV-miRNA profiles were analysed at both time-points in each heifer. Approximately 200 known miRNAs were detected in each sample with 144 commonly detected in all samples. A total of 12 and 14 miRNAs were dysregulated in UH vs. UL and in SH vs. SL heifers, respectively. Interestingly, miR-206 and miR-6517 exhibited the same differential expression pattern in H compared to L heifers both before and after superstimulation. Pathway analysis indicated that circadian rhythm and signaling pathways were among the top pathways enriched with genes targeted by dysregulated miRNAs in H vs. L responding heifers. In conclusion, heifers with divergent ovarian responses exhibited differential expression of plasma EV-miRNAs which may be used as a potential biomarker to predict superstimulation response.