Recent Advances and Potential Multi-Omics Approaches in the Early Phases of Inflammatory Bowel Disease.

Inflammatory bowel disease leads to debilitating gastrointestinal symptoms and reduced quality of life, resulting in a significant burden on healthcare utilization and costs. Despite substantial advancements in diagnosis and treatment, there may still be considerable delays in diagnosing some patients. To reduce disease progression before the full disease spectrum appears and improve prognostic outcomes, several strategies have concentrated on early intervention and prevention. Recent evidence shows that initial immune response changes and endoscopic lesions may exist for years before diagnosis, implying the existence of a preclinical phase of inflammatory bowel disease comparable to findings in other immune-mediated disorders. In this review, we highlight the most relevant findings regarding preclinical inflammatory bowel disease and the prospective role of novel omics techniques in this field.

microRNA-based signatures obtained from endometrial fluid identify implantative endometrium.

STUDY QUESTION: Is it possible to use free and extracellular vesicle-associated microRNAs (miRNAs) from human endometrial fluid (EF) samples as non-invasive biomarkers for implantative endometrium? SUMMARY ANSWER: The free and extracellular vesicle-associated miRNAs can be used to detect implantative endometrium in a non-invasive manner. WHAT IS KNOWN ALREADY: miRNAs and extracellular vesicles (EVs) from EF have been described as mediators of the embryo-endometrium crosstalk. Therefore, the analysis of miRNA from this fluid could become a non-invasive technique for recognizing implantative endometrium. This analysis could potentially help improve the implantation rates in ART. STUDY DESIGN, SIZE, DURATION: In this prospective study, we first optimized different protocols for EVs and miRNA analyses using the EF of a setup cohort (n = 72). Then, we examined differentially expressed miRNAs in the EF of women with successful embryo implantation (discovery cohort n = 15/validation cohort n = 30) in comparison with those for whom the implantation had failed (discovery cohort n = 15/validation cohort n = 30). Successful embryo implantation was considered when pregnancy was confirmed by vaginal ultrasound showing a gestational sac 4 weeks after embryo transfer (ET). PARTICIPANTS/MATERIALS, SETTING, METHODS: The EF of the setup cohort was obtained before starting fertility treatment during the natural cycle, 16-21 days after the beginning of menstruation. For the discovery and validation cohorts, the EF was collected from women undergoing frozen ET on Day 5, and the samples were collected immediately before ET. In this study, we compared five different methods; two of them based on direct extraction of RNA and the other three with an EV enrichment step before the RNA extraction. Small RNA sequencing was performed to determine the most efficient method and find a predictive model differentiating between implantative and non-implantative endometrium. The models were confirmed using quantitative PCR in two sets of samples (discovery and validation cohorts) with different implantation outcomes. MAIN RESULTS AND THE ROLE OF CHANCE: The protocols using EV enrichment detected more miRNAs than the methods based on direct RNA extraction. The two most efficient protocols (using polymer-based precipitation (PBP): PBP-M and PBP-N) were used to obtain two predictive models (based on three miRNAs) allowing us to distinguish between an implantative and non-implantative endometrium. The first Model 1 (PBP-M) (discovery: AUC = 0.93; P-value = 0.003; validation: AUC = 0.69; P-value = 0.019) used hsa-miR-200b-3p, hsa-miR-24-3p and hsa-miR-148b-3p. Model 2 (PBP-N) (discovery: AUC = 0.92; P-value = 0.0002; validation: AUC = 0.78; P-value = 0.0002) used hsa-miR-200b-3p, hsa-miR-24-3p and hsa-miR-99b-5p. Functional analysis of these miRNAs showed strong association with key implantation processes such as in utero embryonic development or transforming growth factor-beta signaling. LARGE SCALE DATA: The FASTQ data are available in the GEO database (access number GSE178917). LIMITATIONS, REASONS FOR CAUTION: One important factor to consider is the inherent variability among the women involved in the trial and among the transferred embryos. The embryos were pre-selected based on morphology, but neither genetic nor molecular studies were conducted, which would have improved the accuracy of our tests. In addition, a limitation in miRNA library construction is the low amount of input RNA. WIDER IMPLICATIONS OF THE FINDINGS: We describe new non-invasive protocols to analyze miRNAs from small volumes of EF. These protocols could be implemented in clinical practice to assess the status of the endometrium before attempting ET. Such evaluation could help to avoid the loss of embryos transferred to a non-implantative endometrium. STUDY FUNDING/COMPETING INTEREST(S): J.I.-P. was supported by a predoctoral grant from the Basque Government (PRE_2017_0204). This study was partially funded by the Grant for Fertility Innovation (GFI, 2011) from Merck (Darmstadt, Germany). It was also supported by the Spanish Ministry of Economy and Competitiveness MINECO within the National Plan RTI2018-094969-B-I00, the European Union's Horizon 2020 research and innovation program (860303), the Severo Ochoa Centre of Excellence Innovative Research Grant (SEV-2016-0644) and the Instituto de Salud Carlos III (PI20/01131). The funding entities did not play any role in the study design, collection, analysis and interpretation of data, writing of the report or the decision to submit the article for publication. The authors declare no competing interests.

Local genetic variation of inflammatory bowel disease in Basque population and its effect in risk prediction.

Inflammatory bowel disease (IBD) is characterised by chronic inflammation of the gastrointestinal tract. Although its aetiology remains unknown, environmental and genetic factors are involved in its development. Regarding genetics, more than 200 loci have been associated with IBD but the transferability of those signals to the Basque population living in Northern Spain, a population with distinctive genetic background, remains unknown. We have analysed 5,411,568 SNPs in 498 IBD cases and 935 controls from the Basque population. We found 33 suggestive loci (p < 5 × 10-6) in IBD and its subtypes, namely Crohn's Disease (CD) and Ulcerative Colitis (UC), detecting a genome-wide significant locus located in HLA region in patients with UC. Those loci contain previously associated genes with IBD (IL23R, JAK2 or HLA genes) and new genes that could be involved in its development (AGT, BZW2 or FSTL1). The overall genetic correlation between European populations and Basque population was high in IBD and CD, while in UC was lower. Finally, the use of genetic risk scores based on previous GWAS findings reached area under the curves > 0.68. In conclusion, we report on the genetic architecture of IBD in the Basque population, and explore the performance of European-descent genetic risk scores in this population.

Metabolic subtypes of patients with NAFLD exhibit distinctive cardiovascular risk profiles.

BACKGROUND AND AIMS: We previously identified subsets of patients with NAFLD with different metabolic phenotypes. Here we align metabolomic signatures with cardiovascular disease (CVD) and genetic risk factors. APPROACH AND RESULTS: We analyzed serum metabolome from 1154 individuals with biopsy-proven NAFLD, and from four mouse models of NAFLD with impaired VLDL-triglyceride (TG) secretion, and one with normal VLDL-TG secretion. We identified three metabolic subtypes: A (47%), B (27%), and C (26%). Subtype A phenocopied the metabolome of mice with impaired VLDL-TG secretion; subtype C phenocopied the metabolome of mice with normal VLDL-TG; and subtype B showed an intermediate signature. The percent of patients with NASH and fibrosis was comparable among subtypes, although subtypes B and C exhibited higher liver enzymes. Serum VLDL-TG levels and secretion rate were lower among subtype A compared with subtypes B and C. Subtype A VLDL-TG and VLDL-apolipoprotein B concentrations were independent of steatosis, whereas subtypes B and C showed an association with these parameters. Serum TG, cholesterol, VLDL, small dense LDL5,6 , and remnant lipoprotein cholesterol were lower among subtype A compared with subtypes B and C. The 10-year high risk of CVD, measured with the Framingham risk score, and the frequency of patatin-like phospholipase domain-containing protein 3 NAFLD risk allele were lower in subtype A. CONCLUSIONS: Metabolomic signatures identify three NAFLD subgroups, independent of histological disease severity. These signatures align with known CVD and genetic risk factors, with subtype A exhibiting a lower CVD risk profile. This may account for the variation in hepatic versus cardiovascular outcomes, offering clinically relevant risk stratification.

Results of the Seventh Scientific Workshop of ECCO: Precision Medicine in IBD-Disease Outcome and Response to Therapy.

Inflammatory bowel diseases [IBD] are a heterogeneous spectrum with two extreme phenotypes, Crohn's disease [CD] and ulcerative colitis [UC], which both represent numerous phenotypical variations. Hence, we should no longer approach all IBD patients similarly, but rather aim to rethink clinical classifications and modify treatment algorithms to usher in a new era of precision medicine in IBD. This scientific ECCO workshop aims to provide a state-of-the-art overview on prognostic and predictive markers, shed light on key questions in biomarker development, propose best practices in IBD biomarker development [including trial design], and discuss the potential for multi-omic data integration to help drive further advances to make precision medicine a reality in IBD.

Clinical and biological predictors of response to standardised paediatric colitis therapy (PROTECT): a multicentre inception cohort study.

BACKGROUND: Lack of evidence-based outcomes data leads to uncertainty in developing treatment regimens in children who are newly diagnosed with ulcerative colitis. We hypothesised that pretreatment clinical, transcriptomic, and microbial factors predict disease course. METHODS: In this inception cohort study, we recruited paediatric patients aged 4-17 years with newly diagnosed ulcerative colitis from 29 centres in the USA and Canada. Patients initially received standardised mesalazine or corticosteroids, with pre-established criteria for escalation to immunomodulators (ie, thiopurines) or anti-tumor necrosis factor-α (TNFα) therapy. We used RNA sequencing to define rectal gene expression before treatment, and 16S sequencing to characterise rectal and faecal microbiota. The primary outcome was week 52 corticosteroid-free remission with no therapy beyond mesalazine. We assessed factors associated with the primary outcome using logistic regression models of the per-protocol population. This study is registered with ClinicalTrials.gov, number NCT01536535. FINDINGS: Between July 10, 2012, and April 21, 2015, of 467 patients recruited, 428 started medical therapy, of whom 400 (93%) were evaluable at 52 weeks and 386 (90%) completed the study period with no protocol violations. 150 (38%) of 400 participants achieved week 52 corticosteroid-free remission, of whom 147 (98%) were taking mesalazine and three (2%) were taking no medication. 74 (19%) of 400 were escalated to immunomodulators alone, 123 (31%) anti-TNFα therapy, and 25 (6%) colectomy. Low baseline clinical severity, high baseline haemoglobin, and week 4 clinical remission were associated with achieving week 52 corticosteroid-free remission (n=386, logistic model area under the curve [AUC] 0·70, 95% CI 0·65-0·75; specificity 77%, 95% CI 71-82). Baseline severity and remission by week 4 were validated in an independent cohort of 274 paediatric patients with newly diagnosed ulcerative colitis. After adjusting for clinical predictors, an antimicrobial peptide gene signature (odds ratio [OR] 0·57, 95% CI 0·39-0·81; p=0·002) and abundance of Ruminococcaceae (OR 1·43, 1·02-2·00; p=0·04), and Sutterella (OR 0·81, 0·65-1·00; p=0·05) were independently associated with week 52 corticosteroid-free remission. INTERPRETATION: Our findings support the utility of initial clinical activity and treatment response by 4 weeks to predict week 52 corticosteroid-free remission with mesalazine alone in children who are newly diagnosed with ulcerative colitis. The development of personalised clinical and biological signatures holds the promise of informing ulcerative colitis therapeutic decisions. FUNDING: US National Institutes of Health.

Blood-Derived DNA Methylation Signatures of Crohn's Disease and Severity of Intestinal Inflammation.

BACKGROUND & AIMS: Crohn's disease is a relapsing and remitting inflammatory disorder with a variable clinical course. Although most patients present with an inflammatory phenotype (B1), approximately 20% of patients rapidly progress to complicated disease, which includes stricturing (B2), within 5 years. We analyzed DNA methylation patterns in blood samples of pediatric patients with Crohn's disease at diagnosis and later time points to identify changes that associate with and might contribute to disease development and progression. METHODS: We obtained blood samples from 164 pediatric patients (1-17 years old) with Crohn's disease (B1 or B2) who participated in a North American study and were followed for 5 years. Participants without intestinal inflammation or symptoms served as controls (n = 74). DNA methylation patterns were analyzed in samples collected at time of diagnosis and 1-3 years later at approximately 850,000 sites. We used genetic association and the concept of Mendelian randomization to identify changes in DNA methylation patterns that might contribute to the development of or result from Crohn's disease. RESULTS: We identified 1189 5'-cytosine-phosphate-guanosine-3' (CpG) sites that were differentially methylated between patients with Crohn's disease (at diagnosis) and controls. Methylation changes at these sites correlated with plasma levels of C-reactive protein. A comparison of methylation profiles of DNA collected at diagnosis of Crohn's disease vs during the follow-up period showed that, during treatment, alterations identified in methylation profiles at the time of diagnosis of Crohn's disease more closely resembled patterns observed in controls, irrespective of disease progression to B2. We identified methylation changes at 3 CpG sites that might contribute to the development of Crohn's disease. Most CpG methylation changes associated with Crohn's disease disappeared with treatment of inflammation and might be a result of Crohn's disease. CONCLUSIONS: Methylation patterns observed in blood samples from patients with Crohn's disease accompany acute inflammation; with treatment, these change to resemble methylation patterns observed in patients without intestinal inflammation. These findings indicate that Crohn's disease-associated patterns of DNA methylation observed in blood samples are a result of the inflammatory features of the disease and are less likely to contribute to disease development or progression.

Distinct Effector B Cells Induced by Unregulated Toll-like Receptor 7 Contribute to Pathogenic Responses in Systemic Lupus Erythematosus.

Systemic Lupus Erythematosus (SLE) is characterized by B cells lacking IgD and CD27 (double negative; DN). We show that DN cell expansions reflected a subset of CXCR5- CD11c+ cells (DN2) representing pre-plasma cells (PC). DN2 cells predominated in African-American patients with active disease and nephritis, anti-Smith and anti-RNA autoantibodies. They expressed a T-bet transcriptional network; increased Toll-like receptor-7 (TLR7); lacked the negative TLR regulator TRAF5; and were hyper-responsive to TLR7. DN2 cells shared with activated naive cells (aNAV), phenotypic and functional features, and similar transcriptomes. Their PC differentiation and autoantibody production was driven by TLR7 in an interleukin-21 (IL-21)-mediated fashion. An in vivo developmental link between aNAV, DN2 cells, and PC was demonstrated by clonal sharing. This study defines a distinct differentiation fate of autoreactive naive B cells into PC precursors with hyper-responsiveness to innate stimuli, as well as establishes prominence of extra-follicular B cell activation in SLE, and identifies therapeutic targets.

Pervasive Modulation of Obesity Risk by the Environment and Genomic Background.

The prevalence of the so-called diseases of affluence, such as type 2 diabetes or hypertension, has increased dramatically in the last two generations. Although genome-wide association studies (GWAS) have discovered hundreds of genes involved in disease etiology, the sudden increase in disease incidence suggests a major role for environmental risk factors. Obesity constitutes a case example of a modern trait shaped by contemporary environment, although with considerable debates about the extent to which gene-by-environment (G×E) interactions accentuate obesity risk in individuals following obesogenic lifestyles. Although interaction effects have been robustly confirmed at the FTO locus, accumulating evidence at the genome-wide level implicates a role for polygenic risk-by-environment interactions. Through a variety of analyses using the UK Biobank, we confirm that the genomic background plays a major role in shaping the expressivity of alleles that increase body mass index (BMI).

Disease-specific regulation of gene expression in a comparative analysis of juvenile idiopathic arthritis and inflammatory bowel disease.

BACKGROUND: The genetic and immunological factors that contribute to differences in susceptibility and progression between sub-types of inflammatory and autoimmune diseases continue to be elucidated. Inflammatory bowel disease and juvenile idiopathic arthritis are both clinically heterogeneous and known to be due in part to abnormal regulation of gene activity in diverse immune cell types. Comparative genomic analysis of these conditions is expected to reveal differences in underlying genetic mechanisms of disease. METHODS: We performed RNA-Seq on whole blood samples from 202 patients with oligoarticular, polyarticular, or systemic juvenile idiopathic arthritis, or with Crohn's disease or ulcerative colitis, as well as healthy controls, to characterize differences in gene expression. Gene ontology analysis combined with Blood Transcript Module and Blood Informative Transcript analysis was used to infer immunological differences. Comparative expression quantitative trait locus (eQTL) analysis was used to quantify disease-specific regulation of transcript abundance. RESULTS: A pattern of differentially expressed genes and pathways reveals a gradient of disease spanning from healthy controls to oligoarticular, polyarticular, and systemic juvenile idiopathic arthritis (JIA); Crohn's disease; and ulcerative colitis. Transcriptional risk scores also provide good discrimination of controls, JIA, and IBD. Most eQTL are found to have similar effects across disease sub-types, but we also identify disease-specific eQTL at loci associated with disease by GWAS. CONCLUSION: JIA and IBD are characterized by divergent peripheral blood transcriptomes, the genetic regulation of which displays limited disease specificity, implying that disease-specific genetic influences are largely independent of, or downstream of, eQTL effects.

Correction: Genetic factors affecting EBV copy number in lymphoblastoid cell lines derived from the 1000 Genome Project samples.

[This corrects the article DOI: 10.1371/journal.pone.0179446.].

Replicability and Prediction: Lessons and Challenges from GWAS.

Since the publication of the Wellcome Trust Case Control Consortium (WTCCC) landmark study a decade ago, genome-wide association studies (GWAS) have led to the discovery of thousands of risk variants involved in disease etiology. This success story has two angles that are often overlooked. First, GWAS findings are highly replicable. This is an unprecedented phenomenon in complex trait genetics, and indeed in many areas of science, which in past decades have been plagued by false positives. At a time of increasing concerns about the lack of reproducibility, we examine the biological and methodological reasons that account for the replicability of GWAS and identify the challenges ahead. In contrast to the exemplary success of disease gene discovery, at present GWAS findings are not useful for predicting phenotypes. We close with an overview of the prospects for individualized prediction of disease risk and its foreseeable impact in clinical practice.

Enhanced Contribution of HLA in Pediatric Onset Ulcerative Colitis.

Background: The genetic contributions to pediatric onset ulcerative colitis (UC), characterized by severe disease and extensive colonic involvement, are largely unknown. In adult onset UC, Genome Wide Association Study (GWAS) has identified numerous loci, most of which have a modest susceptibility risk (OR 0.84-1.14), with the exception of the human leukocyte antigen (HLA) region on Chromosome 6 (OR 3.59). Method: To study the genetic contribution to exclusive pediatric onset UC, a GWAS was performed on 466 cases with 2099 healthy controls using UK Biobank array. SNP2HLA was used to impute classical HLA alleles and their corresponding amino acids, and the results are compared with adult onset UC. Results: HLA explained the almost entire association signal, dominated with 191 single nucleotide polymorphisms (SNPs) (p = 5 x 10-8 to 5 x 10-10). Although very small effects, established SNPs in adult onset UC loci had similar direction and magnitude in pediatric onset UC. SNP2HLA imputation identified HLA-DRB1*0103 (odds ratio [OR] = 6.941, p = 1.92*10-13) as the most significant association for pediatric UC compared with adult onset UC (OR = 3.59). Further conditioning showed independent effects for HLA-DRB1*1301 (OR = 2.25, p = 7.92*10-9) and another SNP rs17188113 (OR = 0.48, p = 7.56*10-9). Two HLA-DRB1 causal alleles are shared with adult onset UC, while at least 2 signals are unique to pediatric UC. Subsequent stratified analyses indicated that HLA-DRB1*0103 has stronger association for extensive disease (E4: OR = 8.28, p = 4.66x10-10) and female gender (OR = 8.85, p = 4.82x10-13). Conclusion: In pediatric onset UC, the HLA explains almost the entire genetic associations. In addition, the HLA association is approximately twice as strong in pediatric UC compared with adults, due to a combination of novel and shared effects. We speculate the paramount importance of antigenic stimulation either by infectious or noninfectious stimuli as a causal event in pediatric UC onset.

Transcriptional risk scores link GWAS to eQTLs and predict complications in Crohn's disease.

Gene expression profiling can be used to uncover the mechanisms by which loci identified through genome-wide association studies (GWAS) contribute to pathology. Given that most GWAS hits are in putative regulatory regions and transcript abundance is physiologically closer to the phenotype of interest, we hypothesized that summation of risk-allele-associated gene expression, namely a transcriptional risk score (TRS), should provide accurate estimates of disease risk. We integrate summary-level GWAS and expression quantitative trait locus (eQTL) data with RNA-seq data from the RISK study, an inception cohort of pediatric Crohn's disease. We show that TRSs based on genes regulated by variants linked to inflammatory bowel disease (IBD) not only outperform genetic risk scores (GRSs) in distinguishing Crohn's disease from healthy samples, but also serve to identify patients who in time will progress to complicated disease. Our dissection of eQTL effects may be used to distinguish genes whose association with disease is through promotion versus protection, thereby linking statistical association to biological mechanism. The TRS approach constitutes a potential strategy for personalized medicine that enhances inference from static genotypic risk assessment.

Antagonistic pleiotropy and mutation accumulation influence human senescence and disease.

Senescence has long been a public health challenge as well as a fascinating evolutionary problem. There is neither a universally accepted theory for its ultimate causes, nor a consensus about what may be its impact on human health. Here we test the predictions of two evolutionary explanations of senescence-mutation accumulation and antagonistic pleiotropy-which postulate that genetic variants with harmful effects in old ages can be tolerated, or even favoured, by natural selection at early ages. Using data from genome-wide association studies (GWAS), we study the effects of genetic variants associated with diseases appearing at different periods in life, when they are expected to have different impacts on fitness. Data fit theoretical expectations. Namely, we observe higher risk allele frequencies combined with large effect sizes for late-onset diseases, and detect a significant excess of early-late antagonistically pleiotropic variants that, strikingly, tend to be harboured by genes related to ageing. Beyond providing systematic, genome-wide evidence for evolutionary theories of senescence in our species and contributing to the long-standing question of whether senescence is the result of adaptation, our approach reveals relationships between previously unrelated pathologies, potentially contributing to tackling the problem of an ageing population.

Constraints on eQTL Fine Mapping in the Presence of Multisite Local Regulation of Gene Expression.

Expression quantitative trait locus (eQTL) detection has emerged as an important tool for unraveling of the relationship between genetic risk factors and disease or clinical phenotypes. Most studies use single marker linear regression to discover primary signals, followed by sequential conditional modeling to detect secondary genetic variants affecting gene expression. However, this approach assumes that functional variants are sparsely distributed and that close linkage between them has little impact on estimation of their precise location and the magnitude of effects. We describe a series of simulation studies designed to evaluate the impact of linkage disequilibrium (LD) on the fine mapping of causal variants with typical eQTL effect sizes. In the presence of multisite regulation, even though between 80 and 90% of modeled eSNPs associate with normally distributed traits, up to 10% of all secondary signals could be statistical artifacts, and at least 5% but up to one-quarter of credible intervals of SNPs within r2 > 0.8 of the peak may not even include a causal site. The Bayesian methods eCAVIAR and DAP (Deterministic Approximation of Posteriors) provide only modest improvement in resolution. Given the strong empirical evidence that gene expression is commonly regulated by more than one variant, we conclude that the fine mapping of causal variants needs to be adjusted for multisite influences, as conditional estimates can be highly biased by interference among linked sites, but ultimately experimental verification of individual effects is needed. Presumably similar conclusions apply not just to eQTL mapping, but to multisite influences on fine mapping of most types of quantitative trait.