Sphingosine 1-phosphate modulation and immune cell trafficking in inflammatory bowel disease.

Immune cell trafficking is a critical element of the intestinal immune response, both in homeostasis and in pathological conditions associated with inflammatory bowel disease (IBD). This process involves adhesion molecules, chemoattractants and receptors expressed on immune cell surfaces, blood vessels and stromal intestinal tissue as well as signalling pathways, including those modulated by sphingosine 1-phosphate (S1P). The complex biological processes of leukocyte recruitment, activation, adhesion and migration have been targeted by various monoclonal antibodies (vedolizumab, etrolizumab, ontamalimab). Promising preclinical and clinical data with several oral S1P modulators suggest that inhibition of lymphocyte egress from the lymph nodes to the bloodstream might be a safe and efficacious alternative mechanism for reducing inflammation in immune-mediated disorders, including Crohn's disease and ulcerative colitis. Although various questions remain, including the potential positioning of S1P modulators in treatment algorithms and their long-term safety, this novel class of compounds holds great promise. This Review summarizes the critical mediators and mechanisms involved in immune cell trafficking in IBD and the available evidence for efficacy, safety and pharmacokinetics of S1P receptor modulators in IBD and other immune-mediated disorders. Further, it discusses potential future approaches to incorporate S1P modulators into the treatment of IBD.

Combined Histo-endoscopic Remission but not Endoscopic Healing Alone in Ulcerative Colitis is Associated with a Mucosal Transcriptional Profile Resembling Healthy Mucosa.

BACKGROUND AND AIMS: A composite endpoint of histological and endoscopic remission is proposed to be the most complete measure of mucosal healing in ulcerative colitis [UC]. We aim to establish the prognosis, and transcriptional and microbial features of histo-endoscopic remission and activity. METHODS: A cross-sectional endoscopic rectosigmoid colon sample collection from UC patients and healthy controls [HC] was performed for histopathology and host genome-wide RNA-sequencing. Histo-endoscopic remission and histo-endoscopic activity were defined as Mayo endoscopic subscore [MES] 0-1 with and without histological activity, respectively. UC relapse, defined as symptomatic and endoscopic worsening, was retrospectively recorded for survival analysis. Unsupervised and differential gene expression analyses were performed, and the interaction between transcriptomics and mucosal gut microbiota was analysed based on the 16S rRNA gene sequencing profile. RESULTS: UC patients with histo-endoscopic remission showed a significantly lower risk of relapse compared to histo-endoscopic activity. Unsupervised analysis of the transcriptomic profile showed that histo-endoscopic remission and histo-endoscopic activity samples clustered with HC and MES 2-3 samples, respectively. A total of 452 host genes enriched for humoral immune response, antimicrobial defence, chemokine and TH17 signalling pathway were upregulated in histo-endoscopic activity compared to histo-endoscopic remission. A set of host genes with antimicrobial properties showed significant associations with mucosal microbiota. CONCLUSIONS: The rectosigmoid mucosa transcriptional profile of UC patients in histo-endoscopic remission resembles that of HC mucosa and confers a lower risk of relapse. These data support that the combination of histo-endoscopic remission could be the most appropriate definition of mucosal healing in UC.

A Role for CXCR3 Ligands as Biomarkers of Post-Operative Crohn's Disease Recurrence.

BACKGROUND AND AIMS: Crohn's disease [CD] recurrence following ileocolic resection [ICR] is common. We sought to identify blood-based biomarkers associated with CD recurrence. METHODS: CD patients undergoing ICR were recruited across six centres. Serum samples were obtained at post-operative colonoscopy. A multiplex immunoassay was used to analyse 92 inflammation-related proteins [Olink Proteomics]. Bayesian analysis was used to identify proteins associated with increasing Rutgeerts score. Identified proteins were used in receiver operating characteristic [ROC] analysis to examine the ability to identify CD recurrence [Rutgeerts score ≥i2]. Existing single cell data were interrogated to further elucidate the role of the identified proteins. RESULTS: Data from 276 colonoscopies in 213 patients were available. Median time from surgery to first and second colonoscopy was 7 (interquartile range [IQR] 6-9) and 19 [IQR 16-23] months, respectively. Disease recurrence was evident at 60 [30%] first and 36 [49%] second colonoscopies. Of 14 proteins significantly associated with Rutgeerts score, the strongest signal was seen for CXCL9 and MMP1. Among patients on anti-tumour necrosis factor drugs, CXCL9 and CXCL11 were most strongly associated with Rutgeerts score. Both are CXCR3 ligands. Incorporation of identified proteins into ROC analysis improved the ability to identify disease recurrence as compared to C-reactive protein alone: area under the curve [AUC] 0.75 (95% confidence interval [CI]: 0.66-0.82] vs 0.64 [95% CI 0.56-0.72], p = 0.012. Single cell transcriptomic data provide evidence that innate immune cells are the primary source of the identified proteins. CONCLUSIONS: CXCR3 ligands are associated with CD recurrence following ICR. Incorporation of novel blood-based candidate biomarkers may aid in identification of CD recurrence.

Meta-analysis of gene expression disease signatures in colonic biopsy tissue from patients with ulcerative colitis.

Publicly available ulcerative colitis (UC) gene expression datasets from observational studies and clinical trials include inherently heterogeneous disease characteristics and methodology. We used meta-analysis to identify a robust UC gene signature from inflamed biopsies. Eight gene expression datasets derived from biopsy tissue samples from noninflammatory bowel disease (IBD) controls and areas of active inflammation from patients with UC were publicly available. Expression- and meta-data were downloaded with GEOquery. Differentially expressed genes (DEG) in individual datasets were defined as those with fold change > 1.5 and a Benjamini-Hochberg adjusted P value < .05. Meta-analysis of all DEG used a random effects model. Reactome pathway enrichment analysis was conducted. Meta-analysis identified 946 up- and 543 down-regulated genes in patients with UC compared to non-IBD controls (1.2 and 1.7 times fewer up- and down-regulated genes than the median of the individual datasets). Top-ranked up- and down-regulated DEG were LCN2 and AQP8. Multiple immune-related pathways (e.g., 'Chemokine receptors bind chemokine' and 'Interleukin-10 signaling') were significantly up-regulated in UC, while 'Biological oxidations' and 'Fatty acid metabolism' were downregulated. A web-based data-mining tool with the meta-analysis results was made available ( https://premedibd.com/genes.html ). A UC inflamed biopsy disease gene signature was derived. This signature may be an unbiased reference for comparison and improve the efficiency of UC biomarker studies by increasing confidence for identification of disease-related genes and pathways.

Integrative Analysis of Colonic Biopsies from Inflammatory Bowel Disease Patients Identifies an Interaction Between Microbial Bile Acid-inducible Gene Abundance and Human Angiopoietin-like 4 Gene Expression.

BACKGROUND AND AIMS: Microbial-derived bile acids can modulate host gene expression, and their faecal abundance is decreased in active inflammatory bowel disease [IBD]. We analysed the impact of endoscopic inflammation on microbial genes involved in bile acid biotransformation, and their interaction with host transcriptome in the intestinal mucosa of IBD patients. METHODS: Endoscopic mucosal biopsies were collected from non-inflamed and inflamed terminal ileum, ascending and sigmoid colon of IBD patients. Prediction of imputed metagenome functional content from 16S rRNA profile and real-time quantitative polymerase chain reaction [qPCR] were utsed to assess microbial bile acid biotransformation gene abundance, and RNA-seq was used for host transcriptome analysis. Linear regression and partial Spearman correlation accounting for age, sex, and IBD type were used to assess the association between microbial genes, inflammation, and host transcriptomics in each biopsy location. A Bayesian network [BN] analysis was fitted to infer the direction of interactions between IBD traits and microbial and host genes. RESULTS: The inferred microbial gene pathway involved in secondary bile acid biosynthesis [ko00121 pathway] was depleted in inflamed terminal ileum of IBD patients compared with non-inflamed tissue. In non-inflamed sigmoid colon, the relative abundance of bile acid-inducible [baiCD] microbial genes was positively correlated with the host Angiopoietin-like 4 [Angptl4] gene expression. The BN analysis suggests that the microbial baiCD gene abundance could affect Angptl4 expression, and this interaction appears to be lost in the presence of inflammation. CONCLUSIONS: Endoscopic inflammation affects the abundance of crucial microbial bile acid-metabolising genes and their interaction with Angptl4 in intestinal mucosa of IBD patients.

Conducting Translational Gastrointestinal Research in the Era of COVID-19.

Spread of the novel coronavirus SARS-CoV-2 has resulted in a global pandemic that is affecting the health and economy of all World Health Organization [WHO] regions. Clinical and translational research activities have been affected drastically by this global catastrophe. In this document we provide a suggested roadmap for resuming gastrointestinal translational research activities, emphasising physical distancing and use of personal protective equipment. We discuss modes of virus transmission in enclosed environments [including clinical workplaces and laboratories] and potential risks of exposure in the endoscopy environment for research staff. The proposed guidelines should be considered in conjunction with local institutional and government guidelines so that translational research can be resumed as safely as possible.

Estimation of genetic parameters for mid-infrared-predicted lactoferrin and milk fat globule size in Holstein cattle.

Lactoferrin (LF) and milk fat globule (MFG) are 2 biologically active components of milk with great economical and nutritional value in the dairy industry. The objectives of this study were to estimate (1) the heritability of mid-infrared (MIR)-predicted LF and MFG size (MFGS) and (2) the genetic correlations between predicted LF and MFGS with milk, fat, and protein yields, fat and protein percentages, and somatic cell score in first-parity Canadian Holstein cattle. A total of 109,029 test-day records from 22,432 cows and 1,572 farms for MIR-predicted LF and 109,212 test-day records from 22,424 cows and 1,559 farms for MIR-predicted MFGS were used in the analyses. Four separate 5-trait random regression test-day models were used. The models included days in milk, herd test date, and a polynomial regression on DIM nested in age-season of calving classes as fixed effects, random polynomial regressions on DIM nested in herd-year of calving, animal additive genetic and permanent environment classes, and a residual effect. Regression curves were modeled using orthogonal Legendre polynomials of order 4 for the fixed age-season of calving effect and of order 5 for the random effects. Moderate overall heritability estimates of 0.34 and 0.46 were estimated for the MIR-predicted LF and MIR-predicted MFGS, respectively. These heritability estimates were similar to the ones estimated for the direct measure of MFGS in a previous study. The genetic correlations between predicted MFGS and fat percentage (0.53) and between predicted LF and protein percentage (0.41) were both moderate and positive. Predicted LF and somatic cell score showed a weaker correlation (0.06) compared with other studies. The moderate genetic correlation between MIR-predicted MFGS and fat percentage and between MIR-predicted LF and protein percentage suggests that MIR predictions of MFGS and LF are not simply a function of the amount of fat and protein percentage, respectively, in the milk (i.e., the prediction equations are not simply predicting fat or protein percentages). Thus, these MIR-predicted values may provide additional information for selecting for fine milk components in Holstein cattle.

A review of traditional and machine learning methods applied to animal breeding.

The current livestock management landscape is transitioning to a high-throughput digital era where large amounts of information captured by systems of electro-optical, acoustical, mechanical, and biosensors is stored and analyzed on a daily and hourly basis, and actionable decisions are made based on quantitative and qualitative analytic results. While traditional animal breeding prediction methods have been used with great success until recently, the deluge of information starts to create a computational and storage bottleneck that could lead to negative long-term impacts on herd management strategies if not handled properly. A plethora of machine learning approaches, successfully used in various industrial and scientific applications, made their way in the mainstream approaches for livestock breeding techniques, and current results show that such methods have the potential to match or surpass the traditional approaches, while most of the time they are more scalable from a computational and storage perspective. This article provides a succinct view on what traditional and novel prediction methods are currently used in the livestock breeding field, how successful they are, and how the future of the field looks in the new digital agriculture era.

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle.

BACKGROUND: Our aim was to identify genomic regions via genome-wide association studies (GWAS) to improve the predictability of genetic merit in Holsteins for 10 calving and 28 body conformation traits. Animals were genotyped using the Illumina Bovine 50 K BeadChip and imputed to the Illumina BovineHD BeadChip (HD). GWAS were performed on 601,717 real and imputed single nucleotide polymorphism (SNP) genotypes using a single-SNP mixed linear model on 4841 Holstein bulls with breeding value predictions and followed by gene identification and in silico functional analyses. The association results were further validated using five scenarios with different numbers of SNPs. RESULTS: Seven hundred and eighty-two SNPs were significantly associated with calving performance at a genome-wise false discovery rate (FDR) of 5%. Most of these significant SNPs were on chromosomes 18 (71.9%), 17 (7.4%), 5 (6.8%) and 7 (2.4%) and mapped to 675 genes, among which 142 included at least one significant SNP and 532 were nearby one (100 kbp). For body conformation traits, 607 SNPs were significant at a genome-wise FDR of 5% and most of them were located on chromosomes 5 (30%), 18 (27%), 20 (13%), 6 (6%), 7 (5%), 14 (5%) and 13 (3%). SNP enrichment functional analyses for calving traits at a FDR of 1% suggested potential biological processes including musculoskeletal movement, meiotic cell cycle, oocyte maturation and skeletal muscle contraction. Furthermore, pathway analyses suggested potential pathways associated with calving performance traits including tight junction, oxytocin signaling, and MAPK signaling (P < 0.10). The prediction ability of the 1206 significant SNPs was between 78 and 83% of the prediction ability of the BovineSNP50 SNPs for calving performance traits and between 35 and 79% for body conformation traits. CONCLUSIONS: Various SNPs that are significantly associated with calving performance are located within or nearby genes with potential roles in tight junction, oxytocin signaling, and MAPK signaling. Combining the significant SNPs or SNPs within or nearby gene(s) from the HD panel with the BovineSNP50 panel yielded a marginal increase in the accuracy of prediction of genomic estimated breeding values for all traits compared to the use of the BovineSNP50 panel alone.

Genome-wide association for milk production and female fertility traits in Canadian dairy Holstein cattle.

BACKGROUND: Genome-wide association studies (GWAS) are a powerful tool for detecting genomic regions explaining variation in phenotype. The objectives of the present study were to identify or refine the positions of genomic regions affecting milk production, milk components and fertility traits in Canadian Holstein cattle, and to use these positions to identify genes and pathways that may influence these traits. RESULT: Several QTL regions were detected for milk production (MILK), fat production (FAT), protein production (PROT) and fat and protein deviation (FATD, PROTD respectively). The identified QTL regions for production traits (including milk production) support previous findings and some overlap with genes with known relevant biological functions identified in earlier studies such as DGAT1 and CPSF1. A significant region on chromosome 21 overlapping with the gene FAM181A and not previous linked to fertility in dairy cattle was identified for the calving to first service interval and days open. A functional enrichment analysis of the GWAS results yielded GO terms consistent with the specific phenotypes tested, for example GO terms GO:0007595 (lactation) and GO:0043627 (response to estrogen) for milk production (MILK), GO:0051057 (positive regulation of small GTPase mediated signal transduction) for fat production (FAT), GO:0040019 (positive regulation of embryonic development) for first service to calving interval (CTFS) and GO:0043268 (positive regulation of potassium ion transport) for days open (DO). In other cases the connection between the enriched GO terms and the traits were less clear, for example GO:0003279 (cardiac septum development) for FAT and GO:0030903 (notochord development) for DO trait. CONCLUSION: The chromosomal regions and enriched pathways identified in this study confirm several previous findings and highlight new regions and pathways that may contribute to variation in production or fertility traits in dairy cattle.