UVPAR: fast detection of functional shifts in duplicate genes

BACKGROUND: The imprint of natural selection on gene sequences is often difficult to detect. A plethora of methods have been devised to detect genetic changes due to selective processes. However, many of those methods depend heavily on underlying assumptions regarding the mode of change of DNA sequences and often require sophisticated mathematical treatments that made them computationally slow. The development of fast and effective methods to detect modifications in the selective constraints of genes is therefore of great interest. RESULTS: We describe UVPAR, a program designed to quickly test for changes in the functional constraints of duplicate genes. Starting with alignments of the proteins encoded by couples of duplicate genes in two different species, UVPAR detects the regions in which modifications of the functional constraints in the paralogs occurred since both species diverged. Sequences can be analyzed with UVPAR in just a few minutes on a standard PC computer. To demonstrate the power of the program, we first show how the results obtained with UVPAR compare to those based on other approaches, using data for vertebrate Hox genes. We then describe a comprehensive study of the RBR family of ubiquitin ligases in which we have performed 529 analyses involving 14 duplicate genes in seven model species. A significant increase in the number of functional shifts was observed for the species Danio rerio and for the gene Ariadne-2. CONCLUSION: These results show that UVPAR can be used to generate sensitive analyses to detect changes in the selection constraints acting on paralogs. The high speed of the program allows its application to genome-scale analyses

Non-bee insects are important contributors to global crop pollination

Wild andmanaged bees arewell documented as effective pollinators of global crops of economic importance. However, the contributions by pollinators other than bees have been little explored despite their potential to contribute to crop production and stability in the face of environmental change. Non-bee pollinators include flies, beetles, moths, butterflies, wasps, ants, birds, and bats, among others. Here we focus on non-bee insects and synthesize 39 field studies from five continents that directly measured the crop pollination services provided by non-bees, honey bees, and other bees to compare the relative contributions of these taxa. Non-bees performed 25-50% of the total number of flower visits. Although non-bees were less effective pollinators than bees per flower visit, they made more visits; thus these two factors compensated for each other, resulting in pollination services rendered by non-bees that were similar to those provided by bees. In the subset of studies that measured fruit set, fruit set increased with non-bee insect visits independently of bee visitation rates, indicating that non-bee insects provide a unique benefit that is not provided by bees. We also show that non-bee insects are not as reliant as bees on the presence of remnant natural or seminatural habitat in the surrounding landscape. These results strongly suggest that non-bee insect pollinators play a significant role in global crop production and respond differently than bees to landscape structure, probably making their crop pollination services more robust to changes in land use. Non-bee insects provide a valuable service and provide potential insurance against bee population declines.Peer Reviewe

CropPol: a dynamic, open and global database on crop pollination

Seventy five percent of the world's food crops benefit from insect pollination. Hence, there has been increased interest in how global change drivers impact this critical ecosystem service. Because standardized data on crop pollination are rarely available, we are limited in our capacity to understand the variation in pollination benefits to crop yield, as well as to anticipate changes in this service, develop predictions, and inform management actions. Here, we present CropPol, a dynamic, open and global database on crop pollination. It contains measurements recorded from 202 crop studies, covering 3,394 field observations, 2,552 yield measurements (i.e. berry weight, number of fruits and kg per hectare, among others), and 47,752 insect records from 48 commercial crops distributed around the globe. CropPol comprises 32 of the 87 leading global crops and commodities that are pollinator dependent. Malus domestica is the most represented crop (32 studies), followed by Brassica napus (22 studies), Vaccinium corymbosum (13 studies), and Citrullus lanatus (12 studies). The most abundant pollinator guilds recorded are honey bees (34.22% counts), bumblebees (19.19%), flies other than Syrphidae and Bombyliidae (13.18%), other wild bees (13.13%), beetles (10.97%), Syrphidae (4.87%), and Bombyliidae (0.05%). Locations comprise 34 countries distributed among Europe (76 studies), Northern America (60), Latin America and the Caribbean (29), Asia (20), Oceania (10), and Africa (7). Sampling spans three decades and is concentrated on 2001-05 (21 studies), 2006-10 (40), 2011-15 (88), and 2016-20 (50). This is the most comprehensive open global data set on measurements of crop flower visitors, crop pollinators and pollination to date, and we encourage researchers to add more datasets to this database in the future. This data set is released for non-commercial use only. Credits should be given to this paper (i.e., proper citation), and the products generated with this database should be shared under the same license terms (CC BY-NC-SA). This article is protected by copyright. All rights reserved

Pollination supply models from a local to global scale

Ecological intensification has been embraced with great interest by the academic sector but is still rarely taken up by farmers because monitoring the state of different ecological functions is not straightforward. Modelling tools can represent a more accessible alternative of measuring ecological functions, which could help promote their use amongst farmers and other decision-makers. In the case of crop pollination, modelling has traditionally followed either a mechanistic or a data-driven approach. Mechanistic models simulate the habitat preferences and foraging behaviour of pollinators, while data-driven models associate georeferenced variables with real observations. Here, we test these two approaches to predict pollination supply and validate these predictions using data from a newly released global dataset on pollinator visitation rates to different crops. We use one of the most extensively used models for the mechanistic approach, while for the data-driven approach, we select from among a comprehensive set of state-of-the-art machine-learning models. Moreover, we explore a mixed approach, where data-derived inputs, rather than expert assessment, inform the mechanistic model. We find that, at a global scale, machine-learning models work best, offering a rank correlation coefficient between predictions and observations of pollinator visitation rates of 0.56. In turn, the mechanistic model works moderately well at a global scale for wild bees other than bumblebees. Biomes characterized by temperate or Mediterranean forests show a better agreement between mechanistic model predictions and observations, probably due to more comprehensive ecological knowledge and therefore better parameterization of input variables for these biomes. This study highlights the challenges of transferring input variables across multiple biomes, as expected given the different composition of species in different biomes. Our results provide clear guidance on which pollination supply models perform best at different spatial scales – the first step towards bridging the stakeholder–academia gap in modelling ecosystem service delivery under ecological intensification

UVPAR: fast detection of functional shifts in duplicate genes

Abstract Background The imprint of natural selection on gene sequences is often difficult to detect. A plethora of methods have been devised to detect genetic changes due to selective processes. However, many of those methods depend heavily on underlying assumptions regarding the mode of change of DNA sequences and often require sophisticated mathematical treatments that made them computationally slow. The development of fast and effective methods to detect modifications in the selective constraints of genes is therefore of great interest. Results We describe UVPAR, a program designed to quickly test for changes in the functional constraints of duplicate genes. Starting with alignments of the proteins encoded by couples of duplicate genes in two different species, UVPAR detects the regions in which modifications of the functional constraints in the paralogs occurred since both species diverged. Sequences can be analyzed with UVPAR in just a few minutes on a standard PC computer. To demonstrate the power of the program, we first show how the results obtained with UVPAR compare to those based on other approaches, using data for vertebrate Hox genes. We then describe a comprehensive study of the RBR family of ubiquitin ligases in which we have performed 529 analyses involving 14 duplicate genes in seven model species. A significant increase in the number of functional shifts was observed for the species Danio rerio and for the gene Ariadne-2. Conclusion These results show that UVPAR can be used to generate sensitive analyses to detect changes in the selection constraints acting on paralogs. The high speed of the program allows its application to genome-scale analyses.</p

Impaired development of neocortical circuits contributes to the neurological alterations in DYRK1A haploinsufficiency syndrome

Autism spectrum disorders are early onset neurodevelopmental disorders characterized by deficits in social communication and restricted repetitive behaviors, yet they are quite heterogeneous in terms of their genetic basis and phenotypic manifestations. Recently, de novo pathogenic mutations in DYRK1A, a chromosome 21 gene associated to neuropathological traits of Down syndrome, have been identified in patients presenting a recognizable syndrome included in the autism spectrum. These mutations produce DYRK1A kinases with partial or complete absence of the catalytic domain, or they represent missense mutations located within this domain. Here, we undertook an extensive biochemical characterization of the DYRK1A missense mutations reported to date and show that most of them, but not all, result in enzymatically dead DYRK1A proteins. We also show that haploinsufficient Dyrk1a+/- mutant mice mirror the neurological traits associated with the human pathology, such as defective social interactions, stereotypic behaviors and epileptic activity. These mutant mice present altered proportions of excitatory and inhibitory neocortical neurons and synapses. Moreover, we provide evidence that alterations in the production of cortical excitatory neurons are contributing to these defects. Indeed, by the end of the neurogenic period, the expression of developmental regulated genes involved in neuron differentiation and/or activity is altered. Therefore, our data indicate that altered neocortical neurogenesis could critically affect the formation of cortical circuits, thereby contributing to the neuropathological changes in DYRK1A haploinsufficiency syndrome.This work was supported by grants from the Spanish Ministry of Economia, Industria y Competitividad (MINECO) (BFU2016-81887-REDT, SAF2013-46676-P and SAF2016-77971-R to M.L.A. and BFU2013-44513 and BFU2016-76141 to S.L.), the Secretariat of Universities and Research-Generalitat de Catalunya (2014SGR674), and the Fundación Alicia Koplowitz (Spain). The group of S.L. acknowledges the support of the MINECO Centro de Excelencia Severo Ochoa Programme and of the CERCA Programme (Generalitat de Catalunya). G.S-E. and M.P.S. acknowledge the support of the National Institute of Neurological Disorders and Stroke of the National Institutes of Health(USA) (P01NS097197) and the Spanish Instituto de Salud Carlos III (ISCIII), (PI13/00865, Fondo Europeo de Desarrollo Regional -FEDER “A way of making Europe”, Spain). M.J.B. is supported by the CIBERER, an initiative of the ISCIII. J.A., E.B. and S.N. were supported by MINECO predoctoral fellowships (AP2012-3064 and BES2011-047472) and G.S-E. by a predoctoral fellowship from the Fundación Conchita Rábago, Spain

A faecal microbiota signature with high specificity for pancreatic cancer

Recent evidence suggests a role for the microbiome in pancreatic ductal adenocarcinoma (PDAC) aetiology and progression. To explore the faecal and salivary microbiota as potential diagnostic biomarkers. We applied shotgun metagenomic and 16S rRNA amplicon sequencing to samples from a Spanish case-control study (n=136), including 57 cases, 50 controls, and 29 patients with chronic pancreatitis in the discovery phase, and from a German case-control study (n=76), in the validation phase. Faecal metagenomic classifiers performed much better than saliva-based classifiers and identified patients with PDAC with an accuracy of up to 0.84 area under the receiver operating characteristic curve (AUROC) based on a set of 27 microbial species, with consistent accuracy across early and late disease stages. Performance further improved to up to 0.94 AUROC when we combined our microbiome-based predictions with serum levels of carbohydrate antigen (CA) 19-9, the only current non-invasive, Food and Drug Administration approved, low specificity PDAC diagnostic biomarker. Furthermore, a microbiota-based classification model confined to PDAC-enriched species was highly disease-specific when validated against 25 publicly available metagenomic study populations for various health conditions (n=5792). Both microbiome-based models had a high prediction accuracy on a German validation population (n=76). Several faecal PDAC marker species were detectable in pancreatic tumour and non-tumour tissue using 16S rRNA sequencing and fluorescence in situ hybridisation. Taken together, our results indicate that non-invasive, robust and specific faecal microbiota-based screening for the early detection of PDAC is feasible

Impaired development of neocortical circuits contributes to the neurological alterations in DYRK1A haploinsufficiency syndrome

Autism spectrum disorders are early onset neurodevelopmental disorders characterized by deficits in social communication and restricted repetitive behaviors, yet they are quite heterogeneous in terms of their genetic basis and phenotypic manifestations. Recently, de novo pathogenic mutations in DYRK1A, a chromosome 21 gene associated to neuropathological traits of Down syndrome, have been identified in patients presenting a recognizable syndrome included in the autism spectrum. These mutations produce DYRK1A kinases with partial or complete absence of the catalytic domain, or they represent missense mutations located within this domain. Here, we undertook an extensive biochemical characterization of the DYRK1A missense mutations reported to date and show that most of them, but not all, result in enzymatically dead DYRK1A proteins. We also show that haploinsufficient Dyrk1a+/- mutant mice mirror the neurological traits associated with the human pathology, such as defective social interactions, stereotypic behaviors and epileptic activity. These mutant mice present altered proportions of excitatory and inhibitory neocortical neurons and synapses. Moreover, we provide evidence that alterations in the production of cortical excitatory neurons are contributing to these defects. Indeed, by the end of the neurogenic period, the expression of developmental regulated genes involved in neuron differentiation and/or activity is altered. Therefore, our data indicate that altered neocortical neurogenesis could critically affect the formation of cortical circuits, thereby contributing to the neuropathological changes in DYRK1A haploinsufficiency syndrome.This work was supported by grants from the Spanish Ministry of Economia, Industria y Competitividad (MINECO) (BFU2016-81887-REDT, SAF2013-46676-P and SAF2016-77971-R to M.L.A. and BFU2013-44513 and BFU2016-76141 to S.L.), the Secretariat of Universities and Research-Generalitat de Catalunya (2014SGR674), and the Fundación Alicia Koplowitz (Spain). The group of S.L. acknowledges the support of the MINECO Centro de Excelencia Severo Ochoa Programme and of the CERCA Programme (Generalitat de Catalunya). G.S-E. and M.P.S. acknowledge the support of the National Institute of Neurological Disorders and Stroke of the National Institutes of Health(USA) (P01NS097197) and the Spanish Instituto de Salud Carlos III (ISCIII), (PI13/00865, Fondo Europeo de Desarrollo Regional -FEDER “A way of making Europe”, Spain). M.J.B. is supported by the CIBERER, an initiative of the ISCIII. J.A., E.B. and S.N. were supported by MINECO predoctoral fellowships (AP2012-3064 and BES2011-047472) and G.S-E. by a predoctoral fellowship from the Fundación Conchita Rábago, Spain

Deciphering the complex interplay between pancreatic cancer, diabetes mellitus subtypes and obesity/BMI through causal inference and mediation analyses

[EN]Objectives: To characterise the association between type 2 diabetes mellitus (T2DM) subtypes (new-onset T2DM (NODM) or long-standing T2DM (LSDM)) and pancreatic cancer (PC) risk, to explore the direction of causation through Mendelian randomisation (MR) analysis and to assess the mediation role of body mass index (BMI). Design: Information about T2DM and related factors was collected from 2018 PC cases and 1540 controls from the PanGenEU (European Study into Digestive Illnesses and Genetics) study. A subset of PC cases and controls had glycated haemoglobin, C-peptide and genotype data. Multivariate logistic regression models were applied to derive ORs and 95% CIs. T2DM and PC-related single nucleotide polymorphism (SNP) were used as instrumental variables (IVs) in bidirectional MR analysis to test for two-way causal associations between PC, NODM and LSDM. Indirect and direct effects of the BMI-T2DM-PC association were further explored using mediation analysis. Results: T2DM was associated with an increased PC risk when compared with non-T2DM (OR=2.50; 95% CI: 2.05 to 3.05), the risk being greater for NODM (OR=6.39; 95% CI: 4.18 to 9.78) and insulin users (OR=3.69; 95% CI: 2.80 to 4.86). The causal association between T2DM (57-SNP IV) and PC was not statistically significant (ORLSDM=1.08, 95% CI: 0.86 to 1.29, ORNODM=1.06, 95% CI: 0.95 to 1.17). In contrast, there was a causal association between PC (40-SNP IV) and NODM (OR=2.85; 95% CI: 2.04 to 3.98), although genetic pleiotropy was present (MR-Egger: p value=0.03). Potential mediating effects of BMI (125-SNPs as IV), particularly in terms of weight loss, were evidenced on the NODM-PC association (indirect effect for BMI in previous years=0.55). Conclusion: Findings of this study do not support a causal effect of LSDM on PC, but suggest that PC causes NODM. The interplay between obesity, PC and T2DM is complex.The work was partially supported by Fondo de Investigaciones Sanitarias (FIS), Instituto de Salud Carlos III, Spain (PI11/01542, PI0902102, PI12/01635, PI12/00815, PI15/01573); Red Temática de Investigación Cooperativa en Cáncer, Spain (RD12/0036/0034, RD12/0036/0050, RD12/0036/0073); WCR (15-0391); European Cooperation in Science and Technology - COST Action BM1204: EUPancreas. EU-6FP Integrated Project (018771-MOLDIAG-PACA), EU-FP7-HEALTH (259737- CANCERALIA, 256974-EPC-TM-Net); Associazione Italiana Ricerca sul Cancro (12182); Cancer Focus Northern Ireland and Department for Employment and Learning; and ALF (Agreement between the Swedish state and some county councils on cooperation on basic education of doctors, medical research, and the development of health care,SLL20130022), Sweden