The genomic landscape of 2,023 colorectal cancers.

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

Role of human plasma metabolites in prediabetes and type 2 diabetes from the IMI-DIRECT study.

AIMS/HYPOTHESIS: Type 2 diabetes is a chronic condition that is caused by hyperglycaemia. Our aim was to characterise the metabolomics to find their association with the glycaemic spectrum and find a causal relationship between metabolites and type 2 diabetes. METHODS: As part of the Innovative Medicines Initiative - Diabetes Research on Patient Stratification (IMI-DIRECT) consortium, 3000 plasma samples were measured with the Biocrates AbsoluteIDQ p150 Kit and Metabolon analytics. A total of 911 metabolites (132 targeted metabolomics, 779 untargeted metabolomics) passed the quality control. Multivariable linear and logistic regression analysis estimates were calculated from the concentration/peak areas of each metabolite as an explanatory variable and the glycaemic status as a dependent variable. This analysis was adjusted for age, sex, BMI, study centre in the basic model, and additionally for alcohol, smoking, BP, fasting HDL-cholesterol and fasting triacylglycerol in the full model. Statistical significance was Bonferroni corrected throughout. Beyond associations, we investigated the mediation effect and causal effects for which causal mediation test and two-sample Mendelian randomisation (2SMR) methods were used, respectively. RESULTS: In the targeted metabolomics, we observed four (15), 34 (99) and 50 (108) metabolites (number of metabolites observed in untargeted metabolomics appear in parentheses) that were significantly different when comparing normal glucose regulation vs impaired glucose regulation/prediabetes, normal glucose regulation vs type 2 diabetes, and impaired glucose regulation vs type 2 diabetes, respectively. Significant metabolites were mainly branched-chain amino acids (BCAAs), with some derivatised BCAAs, lipids, xenobiotics and a few unknowns. Metabolites such as lysophosphatidylcholine a C17:0, sum of hexoses, amino acids from BCAA metabolism (including leucine, isoleucine, valine, N-lactoylvaline, N-lactoylleucine and formiminoglutamate) and lactate, as well as an unknown metabolite (X-24295), were associated with HbA1c progression rate and were significant mediators of type 2 diabetes from baseline to 18 and 48 months of follow-up. 2SMR was used to estimate the causal effect of an exposure on an outcome using summary statistics from UK Biobank genome-wide association studies. We found that type 2 diabetes had a causal effect on the levels of three metabolites (hexose, glutamate and caproate [fatty acid (FA) 6:0]), whereas lipids such as specific phosphatidylcholines (PCs) (namely PC aa C36:2, PC aa C36:5, PC ae C36:3 and PC ae C34:3) as well as the two n-3 fatty acids stearidonate (18:4n3) and docosapentaenoate (22:5n3) potentially had a causal role in the development of type 2 diabetes. CONCLUSIONS/INTERPRETATION: Our findings identify known BCAAs and lipids, along with novel N-lactoyl-amino acid metabolites, significantly associated with prediabetes and diabetes, that mediate the effect of diabetes from baseline to follow-up (18 and 48 months). Causal inference using genetic variants shows the role of lipid metabolism and n-3 fatty acids as being causal for metabolite-to-type 2 diabetes whereas the sum of hexoses is causal for type 2 diabetes-to-metabolite. Identified metabolite markers are useful for stratifying individuals based on their risk progression and should enable targeted interventions.

A Multi-omic Integrative Scheme Characterizes Tissues of Action at Loci Associated with Type 2 Diabetes.

Resolving the molecular processes that mediate genetic risk remains a challenge because most disease-associated variants are non-coding and functional characterization of these signals requires knowledge of the specific tissues and cell-types in which they operate. To address this challenge, we developed a framework for integrating tissue-specific gene expression and epigenomic maps to obtain "tissue-of-action" (TOA) scores for each association signal by systematically partitioning posterior probabilities from Bayesian fine-mapping. We applied this scheme to credible set variants for 380 association signals from a recent GWAS meta-analysis of type 2 diabetes (T2D) in Europeans. The resulting tissue profiles underscored a predominant role for pancreatic islets and, to a lesser extent, adipose and liver, particularly among signals with greater fine-mapping resolution. We incorporated resulting TOA scores into a rule-based classifier and validated the tissue assignments through comparison with data from cis-eQTL enrichment, functional fine-mapping, RNA co-expression, and patterns of physiological association. In addition to implicating signals with a single TOA, we found evidence for signals with shared effects in multiple tissues as well as distinct tissue profiles between independent signals within heterogeneous loci. Lastly, we demonstrated that TOA scores can be directly coupled with eQTL colocalization to further resolve effector transcripts at T2D signals. This framework guides mechanistic inference by directing functional validation studies to the most relevant tissues and can gain power as fine-mapping resolution and cell-specific annotations become richer. This method is generalizable to all complex traits with relevant annotation data and is made available as an R package.

A meta-analysis and a functional study support the influence of mtDNA variant m.16519C on the risk of rapid progression of knee osteoarthritis.

OBJECTIVES: To identify mitochondrial DNA (mtDNA) genetic variants associated with the risk of rapid progression of knee osteoarthritis (OA) and to characterise their functional significance using a cellular model of transmitochondrial cybrids. METHODS: Three prospective cohorts contributed participants. The osteoarthritis initiative (OAI) included 1095 subjects, the Cohort Hip and Cohort Knee included 373 and 326 came from the PROspective Cohort of Osteoarthritis from A Coruña. mtDNA variants were screened in an initial subset of 450 subjects from the OAI by in-depth sequencing of mtDNA. A meta-analysis of the three cohorts was performed. A model of cybrids was constructed to study the functional consequences of harbouring the risk mtDNA variant by assessing: mtDNA copy number, mitochondrial biosynthesis, mitochondrial fission and fusion, mitochondrial reactive oxygen species (ROS), oxidative stress, autophagy and a whole transcriptome analysis by RNA-sequencing. RESULTS: mtDNA variant m.16519C is over-represented in rapid progressors (combined OR 1.546; 95% CI 1.163 to 2.054; p=0.0027). Cybrids with this variant show increased mtDNA copy number and decreased mitochondrial biosynthesis; they produce higher amounts of mitochondrial ROS, are less resistant to oxidative stress, show a lower expression of the mitochondrial fission-related gene fission mitochondrial 1 and an impairment of autophagic flux. In addition, its presence modulates the transcriptome of cybrids, especially in terms of inflammation, where interleukin 6 emerges as one of the most differentially expressed genes. CONCLUSIONS: The presence of the mtDNA variant m.16519C increases the risk of rapid progression of knee OA. Among the most modulated biological processes associated with this variant, inflammation and negative regulation of cellular process stand out. The design of therapies based on the maintenance of mitochondrial function is recommended.

Impact of insufficient sleep on dysregulated blood glucose control under standardised meal conditions.

AIMS/HYPOTHESIS: Sleep, diet and exercise are fundamental to metabolic homeostasis. In this secondary analysis of a repeated measures, nutritional intervention study, we tested whether an individual's sleep quality, duration and timing impact glycaemic response to a breakfast meal the following morning. METHODS: Healthy adults' data (N = 953 [41% twins]) were analysed from the PREDICT dietary intervention trial. Participants consumed isoenergetic standardised meals over 2 weeks in the clinic and at home. Actigraphy was used to assess sleep variables (duration, efficiency, timing) and continuous glucose monitors were used to measure glycaemic variation (>8000 meals). RESULTS: Sleep variables were significantly associated with postprandial glycaemic control (2 h incremental AUC), at both between- and within-person levels. Sleep period time interacted with meal type, with a smaller effect of poor sleep on postprandial blood glucose levels when high-carbohydrate (low fat/protein) (pinteraction = 0.02) and high-fat (pinteraction = 0.03) breakfasts were consumed compared with a reference 75 g OGTT. Within-person sleep period time had a similar interaction (high carbohydrate: pinteraction = 0.001, high fat: pinteraction = 0.02). Within- and between-person sleep efficiency were significantly associated with lower postprandial blood glucose levels irrespective of meal type (both p < 0.03). Later sleep midpoint (time deviation from midnight) was found to be significantly associated with higher postprandial glucose, in both between-person and within-person comparisons (p = 0.035 and p = 0.051, respectively). CONCLUSIONS/INTERPRETATION: Poor sleep efficiency and later bedtime routines are associated with more pronounced postprandial glycaemic responses to breakfast the following morning. A person's deviation from their usual sleep pattern was also associated with poorer postprandial glycaemic control. These findings underscore sleep as a modifiable, non-pharmacological therapeutic target for the optimal regulation of human metabolic health. Trial registration ClinicalTrials.gov NCT03479866.

Genome-wide associations for birth weight and correlations with adult disease.

Birth weight (BW) has been shown to be influenced by both fetal and maternal factors and in observational studies is reproducibly associated with future risk of adult metabolic diseases including type 2 diabetes (T2D) and cardiovascular disease. These life-course associations have often been attributed to the impact of an adverse early life environment. Here, we performed a multi-ancestry genome-wide association study (GWAS) meta-analysis of BW in 153,781 individuals, identifying 60 loci where fetal genotype was associated with BW (P < 5 × 10-8). Overall, approximately 15% of variance in BW was captured by assays of fetal genetic variation. Using genetic association alone, we found strong inverse genetic correlations between BW and systolic blood pressure (Rg = -0.22, P = 5.5 × 10-13), T2D (Rg = -0.27, P = 1.1 × 10-6) and coronary artery disease (Rg = -0.30, P = 6.5 × 10-9). In addition, using large -cohort datasets, we demonstrated that genetic factors were the major contributor to the negative covariance between BW and future cardiometabolic risk. Pathway analyses indicated that the protein products of genes within BW-associated regions were enriched for diverse processes including insulin signalling, glucose homeostasis, glycogen biosynthesis and chromatin remodelling. There was also enrichment of associations with BW in known imprinted regions (P = 1.9 × 10-4). We demonstrate that life-course associations between early growth phenotypes and adult cardiometabolic disease are in part the result of shared genetic effects and identify some of the pathways through which these causal genetic effects are mediated.

Correction: Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

[This corrects the article DOI: 10.1371/journal.pgen.1007813.].

Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

Polycystic ovary syndrome (PCOS) is a disorder characterized by hyperandrogenism, ovulatory dysfunction and polycystic ovarian morphology. Affected women frequently have metabolic disturbances including insulin resistance and dysregulation of glucose homeostasis. PCOS is diagnosed with two different sets of diagnostic criteria, resulting in a phenotypic spectrum of PCOS cases. The genetic similarities between cases diagnosed based on the two criteria have been largely unknown. Previous studies in Chinese and European subjects have identified 16 loci associated with risk of PCOS. We report a fixed-effect, inverse-weighted-variance meta-analysis from 10,074 PCOS cases and 103,164 controls of European ancestry and characterisation of PCOS related traits. We identified 3 novel loci (near PLGRKT, ZBTB16 and MAPRE1), and provide replication of 11 previously reported loci. Only one locus differed significantly in its association by diagnostic criteria; otherwise the genetic architecture was similar between PCOS diagnosed by self-report and PCOS diagnosed by NIH or non-NIH Rotterdam criteria across common variants at 13 loci. Identified variants were associated with hyperandrogenism, gonadotropin regulation and testosterone levels in affected women. Linkage disequilibrium score regression analysis revealed genetic correlations with obesity, fasting insulin, type 2 diabetes, lipid levels and coronary artery disease, indicating shared genetic architecture between metabolic traits and PCOS. Mendelian randomization analyses suggested variants associated with body mass index, fasting insulin, menopause timing, depression and male-pattern balding play a causal role in PCOS. The data thus demonstrate 3 novel loci associated with PCOS and similar genetic architecture for all diagnostic criteria. The data also provide the first genetic evidence for a male phenotype for PCOS and a causal link to depression, a previously hypothesized comorbid disease. Thus, the genetics provide a comprehensive view of PCOS that encompasses multiple diagnostic criteria, gender, reproductive potential and mental health.

Quantitative proteomic profiling of human articular cartilage degradation in osteoarthritis.

Osteoarthritis (OA) is the most common rheumatic pathology and is characterized primarily by articular cartilage degradation. Despite its high prevalence, there is no effective therapy to slow disease progression or regenerate the damaged tissue. Therefore, new diagnostic and monitoring tests for OA are urgently needed, which would also promote the development of alternative therapeutic strategies. In the present study, we have performed an iTRAQ-based quantitative proteomic analysis of secretomes from healthy human articular cartilage explants, comparing their protein profile to those from unwounded (early disease) and wounded (advanced disease) zones of osteoarthritic tissue. This strategy allowed us to identify a panel of 76 proteins that are distinctively released by the diseased tissue. Clustering analysis allowed the classification of proteins according to their different profile of release from cartilage. Among these proteins, the altered release of osteoprotegerin (decreased in OA) and periostin (increased in OA), both involved in bone remodelling processes, was verified in further analyses. Moreover, periostin was also increased in the synovial fluid of OA patients. Altogether, the present work provides a novel insight into the mechanisms of human cartilage degradation and a number of new cartilage-characteristic proteins with possible biomarker value for early diagnosis and prognosis of OA.

Analysis of autoantibody profiles in osteoarthritis using comprehensive protein array concepts.

Osteoarthritis (OA) is the most common rheumatic disease and one of the most disabling pathologies worldwide. To date, the diagnostic methods of OA are very limited, and there are no available medications capable of halting its characteristic cartilage degeneration. Therefore, there is a significant interest in new biomarkers useful for the early diagnosis, prognosis, and therapeutic monitoring. In the recent years, protein microarrays have emerged as a powerful proteomic tool to search for new biomarkers. In this study, we have used two concepts for generating protein arrays, antigen microarrays, and NAPPA (nucleic acid programmable protein arrays), to characterize differential autoantibody profiles in a set of 62 samples from OA, rheumatoid arthritis (RA), and healthy controls. An untargeted screen was performed on 3840 protein fragments spotted on planar antigen arrays, and 373 antigens were selected for validation on bead-based arrays. In the NAPPA approach, a targeted screening was performed on 80 preselected proteins. The autoantibody targeting CHST14 was validated by ELISA in the same set of patients. Altogether, nine and seven disease related autoantibody target candidates were identified, and this work demonstrates a combination of these two array concepts for biomarker discovery and their usefulness for characterizing disease-specific autoantibody profiles.

Genome-wide DNA methylation analysis of articular chondrocytes reveals a cluster of osteoarthritic patients.

OBJECTIVE: Alterations in DNA methylation patterns have been found to correlate with several diseases including osteoarthritis (OA). The aim of this study was to identify, for the first time, the genome-wide DNA methylation profiles of human articular chondrocytes from OA cartilage and healthy control cartilage samples. METHODS: DNA methylation profiling was performed using Illumina Infinium HumanMethylation27 in 25 patients with OA and 20 healthy controls. Subsequent validation was performed by genome-wide expression analysis using the Affymetrix Human Gene 1.1 ST array in an independent cohort of 24 patients with OA. Finally, the most consistent genes in both assays were amplified by quantitative reverse transcriptase PCR in a validation cohort of 48 patients using microfluidic real-time quantitative PCR. Appropriate bioinformatics analyses were carried out using R bioconductor software packages and qBase plus software from Biogazelle. RESULTS: We found 91 differentially methylated (DM) probes, which permitted us to separate patients with OA from healthy controls. Among the patients with OA, we detected 1357 DM probes that identified a tight cluster of seven patients who were different from the rest. This cluster was also identified by genome-wide expression in which 450 genes were differentially expressed. Further validation of the most consistent genes in an independent cohort of patients with OA permitted us to identify this cluster, which was characterised by increased inflammatory processes. CONCLUSIONS: We were able to identify a tight subgroup of patients with OA, characterised by an increased inflammatory response that could be regulated by epigenetics. The identification and isolation of this subgroup may be critical for the development of effective treatment and disease prevention.

An exploratory analysis of associations of genetic variation with the efficacy of tocilizumab in severe COVID-19 patients. A pharmacogenetic study based on next-generation sequencing.

Background: In the context of the cytokine storm the takes place in severe COVID-19 patients, the Interleukin 6 (IL6) pathway emerges as one of the key pathways involved in the pathogenesis of this hyperinflammatory state. The strategy of blocking the inflammatory storm by targeting the IL6 is a promising therapy to mitigate mortality. The use of Tocilizumab was recommended by the World Health Organization (WHO) to treat severe COVID-19 patients. However, the efficacy of Tocilizumab is variable. We hypothesize that the genetic background could be behind the efficacy of Tocilizumab in terms of mortality. Methods: We performed a targeted-next generation sequencing of 287 genes, of which 264 belong to a community panel of ThermoFisher for the study of genetic causes of primary immunodeficiency disorders, and 23 additional genes mostly related to inflammation, not included in the original community panel. This panel was sequenced in an initial cohort of 425 COVID-19 patients, of which 232 were treated with Tocilizumab and standard therapy, and 193 with standard therapy only. Selected genetic variants were genotyped by single base extension in additional 245 patients (95 treated with Tocilizumab and 150 non-treated with Tocilizumab). Appropriate statistical analyses and internal validation, including logistic regression models, with the interaction between Tocilizumab and genetic variants, were applied to assess the impact of these genetic variants in the efficacy of Tocilizumab in terms of mortality. Results: Age (p < 0.001) and cardiovascular disease (p < 0.001) are risk factors for mortality in COVID-19 patients. The presence of GG and TT genotypes at IL10Rß (rs2834167) and IL1ß (rs1143633) genes significantly associates with a reduced risk of mortality in patients treated with Tocilizumab (OR = 0.111; 95%CI = 0.015-0.829; p = 0.010 and OR = 0.378; 95%CI = 0.154-0.924; p = 0.028 respectively). The presence of CC genotype at IL1RN (rs2234679) significantly associates with an increased risk of mortality, but only in patients not treated with Tocilizumab (OR = 3.200; 95%CI = 1.512-6.771; p = 0.002). Exhaustive internal validation using a bootstrap method (B = 500 replicates) validated the accuracy of the predictive models. Conclusion: We developed a series of predictive models based on three genotypes in genes with a strong implication in the etiopathogenesis of COVID-19 disease capable of predicting the risk of mortality in patients treated with Tocilizumab.

A set of circulating microRNAs belonging to the 14q32 chromosome locus identifies two subgroups of individuals with recent-onset type 1 diabetes.

Circulating microRNAs (miRNAs) are linked to the onset and progression of type 1 diabetes mellitus (T1DM), thus representing potential disease biomarkers. In this study, we employed a multiplatform sequencing approach to analyze circulating miRNAs in an extended cohort of prospectively evaluated recent-onset T1DM individuals from the INNODIA consortium. Our findings reveal that a set of miRNAs located within T1DM susceptibility chromosomal locus 14q32 distinguishes two subgroups of individuals. To validate our results, we conducted additional analyses on a second cohort of T1DM individuals, confirming the identification of these subgroups, which we have named cluster A and cluster B. Remarkably, cluster B T1DM individuals, who exhibit increased expression of a set of 14q32 miRNAs, show better glycemic control and display a different blood immunomics profile. Our findings suggest that this set of circulating miRNAs can identify two different T1DM subgroups with distinct blood immunomics at baseline and clinical outcomes during follow-up.

The Association of Cardiometabolic, Diet and Lifestyle Parameters With Plasma Glucagon-like Peptide-1: An IMI DIRECT Study.

CONTEXT: The role of glucagon-like peptide-1 (GLP-1) in type 2 diabetes (T2D) and obesity is not fully understood. OBJECTIVE: We investigate the association of cardiometabolic, diet, and lifestyle parameters on fasting and postprandial GLP-1 in people at risk of, or living with, T2D. METHODS: We analyzed cross-sectional data from the two Innovative Medicines Initiative (IMI) Diabetes Research on Patient Stratification (DIRECT) cohorts, cohort 1 (n = 2127) individuals at risk of diabetes; cohort 2 (n = 789) individuals with new-onset T2D. RESULTS: Our multiple regression analysis reveals that fasting total GLP-1 is associated with an insulin-resistant phenotype and observe a strong independent relationship with male sex, increased adiposity, and liver fat, particularly in the prediabetes population. In contrast, we showed that incremental GLP-1 decreases with worsening glycemia, higher adiposity, liver fat, male sex, and reduced insulin sensitivity in the prediabetes cohort. Higher fasting total GLP-1 was associated with a low intake of wholegrain, fruit, and vegetables in people with prediabetes, and with a high intake of red meat and alcohol in people with diabetes. CONCLUSION: These studies provide novel insights into the association between fasting and incremental GLP-1, metabolic traits of diabetes and obesity, and dietary intake, and raise intriguing questions regarding the relevance of fasting GLP-1 in the pathophysiology T2D.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Effect of okadaic acid on carpet shell clam (Ruditapes decussatus) haemocytes by in vitro exposure and harmful algal bloom simulation assays.

Okadaic acid (OA), produced by dinoflagellates during harmful algal blooms (HAB), belongs to the Diarrheic Shellfish Poisoning toxins that cause gastrointestinal symptoms in humans after consumption. In the present work, Ruditapes decussatus haemocytes were selected to evaluate the effect of OA on cell viability, enzymatic status and immune capacity through the measure by flow cytometry of apoptosis-cell death, non-specific esterase activity and phagocytosis. In order to compare different exposure conditions, two experiments were developed: in vitro exposure to OA and HAB simulation by feeding clams with the OA producer, Prorocentrum lima. Apoptosis was not OA dose-dependent and cell death increased in both assays. Phagocytosis of latex beads and esterase activity decreased in haemocytes incubated with OA. In contrast, esterases increased during the feeding with P. lima. Our results showed that OA and the simulated HAB caused damages on haemocyte functions and viability.

Bivalve omics: state of the art and potential applications for the biomonitoring of harmful marine compounds.

The extraordinary progress experienced by sequencing technologies and bioinformatics has made the development of omic studies virtually ubiquitous in all fields of life sciences nowadays. However, scientific attention has been quite unevenly distributed throughout the different branches of the tree of life, leaving molluscs, one of the most diverse animal groups, relatively unexplored and without representation within the narrow collection of well established model organisms. Within this Phylum, bivalve molluscs play a fundamental role in the functioning of the marine ecosystem, constitute very valuable commercial resources in aquaculture, and have been widely used as sentinel organisms in the biomonitoring of marine pollution. Yet, it has only been very recently that this complex group of organisms became a preferential subject for omic studies, posing new challenges for their integrative characterization. The present contribution aims to give a detailed insight into the state of the art of the omic studies and functional information analysis of bivalve molluscs, providing a timely perspective on the available data resources and on the current and prospective applications for the biomonitoring of harmful marine compounds.

The CHROMEVALOA database: a resource for the evaluation of Okadaic Acid contamination in the marine environment based on the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis.

Okadaic Acid (OA) constitutes the main active principle in Diarrhetic Shellfish Poisoning (DSP) toxins produced during Harmful Algal Blooms (HABs), representing a serious threat for human consumers of edible shellfish. Furthermore, OA conveys critical deleterious effects for marine organisms due to its genotoxic potential. Many efforts have been dedicated to OA biomonitoring during the last three decades. However, it is only now with the current availability of detailed molecular information on DNA organization and the mechanisms involved in the maintenance of genome integrity, that a new arena starts opening up for the study of OA contamination. In the present work we address the links between OA genotoxicity and chromatin by combining Next Generation Sequencing (NGS) technologies and bioinformatics. To this end, we introduce CHROMEVALOAdb, a public database containing the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis (a sentinel model organism) in response to OA exposure. This resource constitutes a leap forward for the development of chromatin-based biomarkers, paving the road towards the generation of powerful and sensitive tests for the detection and evaluation of the genotoxic effects of OA in coastal areas.

The marine toxin okadaic acid induces alterations in the expression level of cancer-related genes in human neuronal cells.

Okadaic acid (OA) is one of the most common and highly distributed marine toxins. It can be accumulated in several molluscs and other marine organisms and cause acute gastrointestinal symptoms after oral consumption by humans, called diarrheic shellfish poisoning. However other toxic effects beyond these gastrointestinal symptoms were also reported. Thus, OA was found to induce important chromosomal abnormalities and other genetic injuries that can lead to severe pathologies, including cancer. Furthermore, the relationship between OA and carcinogenic processes has been previously demonstrated in in vivo studies with rodents, and also suggested in human epidemiological studies. In this context, further research is required to better understand the underlying mechanisms of OA-related tumourigenesis. In a previous study, we identified 247 genes differentially expressed in SHSY5Y neuroblastoma cells exposed to 100nM OA at different times (3, 24 and 48h) by means of suppression subtractive hybridization. These genes were involved in relevant cell functions such as signal transduction, cell cycle, metabolism, and transcription and translation processes. However, due to the high potential percentage of false positives that may be obtained by this approach, results from SSH are recommended to be analyzed by an independent method. In the present study, we selected ten genes related to cancer initiation or progression, directly or indirectly, for further quantitative PCR analysis (ANAPC13, PTTG1, CALM2, CLU, HN1, MALAT1, MAPRE2, MLLT11, SGA-81M and TAX1BP1). Results obtained showed important alterations in the expression patterns of all the genes evaluated at one or more treatment times, providing, for the first time, a possible explanation at the molecular level of the potential relationship between the consumption of OA-contaminated shellfish and the incidence of different cancers in humans. Nevertheless, given the complexity of this process, more exhaustive studies are required before drawing any final conclusion.

A phenome-wide comparative analysis of genetic discordance between obesity and type 2 diabetes.

Obesity and type 2 diabetes are causally related, yet there is considerable heterogeneity in the consequences of both conditions and the mechanisms of action are poorly defined. Here we show a genetic-driven approach defining two obesity profiles that convey highly concordant and discordant diabetogenic effects. We annotate and then compare association signals for these profiles across clinical and molecular phenotypic layers. Key differences are identified in a wide range of traits, including cardiovascular mortality, fat distribution, liver metabolism, blood pressure, specific lipid fractions and blood levels of proteins involved in extracellular matrix remodelling. We find marginal differences in abundance of Bacteroidetes and Firmicutes bacteria in the gut. Instrumental analyses reveal prominent causal roles for waist-to-hip ratio, blood pressure and cholesterol content of high-density lipoprotein particles in the development of diabetes in obesity. We prioritize 17 genes from the discordant signature that convey protection against type 2 diabetes in obesity, which may represent logical targets for precision medicine approaches.