Desarrollo de técnicas bioinformáticas para el análisis de datos de secuenciación masiva en sistemática y genómica evolutiva: Aplicación en el análisis del sistema quimiosensorial en artrópodos

[spa] Las tecnologías de secuenciación de próxima generación (NGS) proporcionan datos potentes para investigar cuestiones biológicas y evolutivas fundamentales, como estudios relacionados con la genómica evolutiva de la adaptación y la filogenética. Actualmente, es posible llevar a cabo proyectos genómicos complejos analizando genomas completos y / o transcriptomas, incluso de organismos no modelo. En esta tesis, hemos realizado dos estudios complementarios utilizando datos NGS. En primer lugar, hemos analizado el transcriptoma (RNAseq) de los principales órganos quimiosensoriales del quelicerado Macrothele calpeiana, Walckenaer, 1805, la única araña protegida en Europa, para investigar el origen y la evolución del sistema quimiosensorial (SQ) en los artrópodos. El SQ es un proceso fisiológico esencial para la supervivencia de los organismos, y está involucrado en procesos biológicos vitales, como la detección de alimentos, parejas o depredadores y sitios de ovoposición. Este sistema, está relativamente bien caracterizado en hexápodos, pero existen pocos estudios en otros linajes de artrópodos. El análisis de nuestro transcriptoma permitió detectar algunos genes expresados en los supuestos órganos quimiosensoriales de los quelicerados, como cinco NPC2 y dos IR. Además, también detectamos 29 tránscritos adicionales después de incluir en los perfiles de HMM nuevos miembros del SQ de genomas de artrópodos recientemente disponibles, como algunos genes de las familias de los SNMP, ENaC, TRP, GR y una OBP-like. Desafortunadamente, muchos de ellos eran fragmentos parciales. En segundo lugar, también hemos desarrollado algunas herramientas bioinformáticas para analizar datos de RNAseq y desarrollar marcadores moleculares. Los investigadores interesados en la aplicación biológica de datos NGS pueden carecer de la experiencia bioinformática requerida para el tratamiento de la gran cantidad de datos generados. En este contexto, principalmente, es necesario el desarrollo de herramientas fáciles de usar para realizar todos los procesos relacionados con el procesamiento básico de datos NGS y la integración de utilidades para realizar análisis posteriores. En esta tesis, hemos desarrollado dos herramientas bioinformáticas con interfaz gráfica, que permite realizar todos los procesos comunes del procesamiento de datos NGS y algunos de los principales análisis posteriores: i) TRUFA (TRanscriptome User-Friendly Analysis), que permite analizar datos RNAseq de organismos que no modelos, incluyendo la anotación funcional y el análisis de expresión génica diferencial; y ii) DOMINO (Development Of Molecular markers In Non-model Organisms), que permite identificar y seleccionar marcadores moleculares apropiados para análisis de biología evolutiva. Estas herramientas han sido validadas utilizando simulaciones por ordenador y datos experimentales, principalmente de arañas.[eng] The Next Generation Sequencing (NGS) technologies are providing powerful data to investigate fundamental biological and evolutionary questions including phylogenetic and adaptive genomic topics. Currently, it is possible to carry out complex genomic projects analyzing the complete genomes and/or transcriptomes even in non-model organisms. In this thesis, we have performed two complementary studies using NGS data. Firstly, we have analyzed the transcriptome (RNAseq) of the main chemosensory organs of the chelicerate Macrothele calpeiana, Walckenaer, 1805, the only spider protected in Europe, to investigate the origin and evolution of the Chemosensory System (CS) in arthropods. The CS is an essential physiological process for the survival of organisms, and it is involved in vital biological processes, such as the detection of food, partners or predators and oviposition sites. This system, which has it relatively well characterized in hexapods, is completely unknown in other arthropod lineages. Our transcriptome analysis allowed to detect some genes expressed in the putative chemosensory organs of chelicerates, such as five NPC2s and two IRs. Furthermore, we detected 29 additional transcripts after including new CS members from recently available genomes in the HMM profiles, such as the SNMPs, ENaCs, TRPs, GRs and one OBP-like. Unfortunately, many of them were partial fragments. Secondly, we have also developed some bioinformatics tools to analyze RNAseq data, and to develop molecular markers. Researchers interested in the biological application of NGS data may lack the bioinformatic expertise required for the treatment of the large amount of data generated. In this context, the development of user-friendly tools for common data processing and the integration of utilities to perform downstream analysis is mostly needed. In this thesis, we have developed two bioinformatics tools with an easy to use graphical interface to perform all the basics processes of the NGS data processing: i) TRUFA (TRanscriptome User-Friendly Analysis), that allows analyzing RNAseq data from non-model organisms, including the functional annotation and differential gene expression analysis; and ii) DOMINO (Development of Molecular markers in Non-model Organisms), which allows identifying and selecting molecular markers appropriated for evolutionary biology analysis. These tools have been validated using computer simulations and experimental data, mainly from spiders

TRUFA: A user-friendly web server for de novo RNA-seq analysis using cluster computing.

Application of next-generation sequencing (NGS) methods for transcriptome analysis (RNA-seq) has become increasingly accessible in recent years and are of great interest to many biological disciplines including, eg, evolutionary biology, ecology, biomedicine, and computational biology. Although virtually any research group can now obtain RNA-seq data, only a few have the bioinformatics knowledge and computation facilities required for transcriptome analysis. Here, we present TRUFA (TRanscriptome User-Friendly Analysis), an open informatics platform offering a web-based interface that generates the outputs commonly used in de novo RNA-seq analysis and comparative transcriptomics. TRUFA provides a comprehensive service that allows performing dynamically raw read cleaning, transcript assembly, annotation, and expression quantification. Due to the computationally intensive nature of such analyses, TRUFA is highly parallelized and benefits from accessing high-performance computing resources. The complete TRUFA pipeline was validated using four previously published transcriptomic data sets. TRUFA's results for the example datasets showed globally similar results when comparing with the original studies, and performed particularly better when analyzing the green tea dataset. The platform permits analyzing RNA-seq data in a fast, robust, and user-friendly manner. Accounts on TRUFA are provided freely upon request at https://trufa.ifca.es

Comparative analysis of tissue-specific transcriptomes in the funnel-web spider Macrothele calpeiana (Araneae, Hexathelidae)

Altres ajuts: Beatriu de Pinós (Generalitat de Catalunya, 2010BP-A 00438 i 2010BP-B 00175)The funnel-web spider Macrothele calpeiana is a charismatic Mygalomorph with a great interest in basic, applied and translational research. Nevertheless, current scarcity of genomic and transcriptomic data of this species clearly limits the research in this non-model organism. To overcome this limitation, we launched the first tissue-specific enriched RNA-seq analysis in this species using a subtractive hybridization approach, with two main objectives, to characterize the specific transcriptome of the putative chemosensory appendages (palps and first pair of legs), and to provide a new set of DNA markers for further phylogenetic studies. We have characterized the set of transcripts specifically expressed in putative chemosensory tissues of this species, much of them showing features shared by chemosensory system genes. Among specific candidates, we have identified some members of the iGluR and NPC2 families. Moreover, we have demonstrated the utility of these newly generated data as molecular markers by inferring the phylogenetic position M. calpeina in the phylogenetic tree of Mygalomorphs. Our results provide novel resources for researchers interested in spider molecular biology and systematics, which can help to expand our knowledge on the evolutionary processes underlying fundamental biological questions, as species invasion or biodiversity origin and maintenance

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome

Data processing; DNA sequencing; Genomic analysisTratamiento de datos; Secuenciación de ADN; Análisis genómicoTractament de dades; Seqüenciació d'ADN; Anàlisi genòmicaMethods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.This research has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 824110 – EASI-Genomics (I.G.G.) and the ERC Synergy project BCLL@las under grant agreement No 810287 (I.G.G.). Institutional support was from the Spanish Instituto de Salud Carlos III, Fondo de Investigaciones Sanitarias and cofunded with ERDF funds (PI19/01772). We acknowledge the institutional support of the Spanish Ministry of Science and Innovation through the Instituto de Salud Carlos III and the 2014–2020 Smart Growth Operating Program, to the EMBL partnership and institutional co-financing with the European Regional Development Fund (MINECO/FEDER, BIO2015-71792-P). We also acknowledge the support of the Centro de Excelencia Severo Ochoa, and the Generalitat de Catalunya through the Departament de Salut, Departament d’Empresa i Coneixement and the CERCA Programme to the institute

The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates

Background We present the draft genome sequence of Dysdera silvatica, a nocturnal ground-dwelling spider from a genus that has undergone a remarkable adaptive radiation in the Canary Islands. Results The draft assembly was obtained using short (Illumina) and long (PaciBio and Nanopore) sequencing reads. Our de novo assembly (1.36 Gb), which represents 80% of the genome size estimated by flow cytometry (1.7 Gb), is constituted by a high fraction of interspersed repetitive elements (53.8%). The assembly completeness, using BUSCO and core eukaryotic genes, ranges from 90% to 96%. Functional annotations based on both ab initio and evidence-based information (including D. silvatica RNA sequencing) yielded a total of 48,619 protein-coding sequences, of which 36,398 (74.9%) have the molecular hallmark of known protein domains, or sequence similarity with Swiss-Prot sequences. The D. silvatica assembly is the first representative of the superfamily Dysderoidea, and just the second available genome of Synspermiata, one of the major evolutionary lineages of the 'true spiders' (Araneomorphae). Conclusions Dysderoids, which are known for their numerous instances of adaptation to underground environments, include some of the few examples of trophic specialization within spiders and are excellent models for the study of cryptic female choice. This resource will be therefore useful as a starting point to study fundamental evolutionary and functional questions, including the molecular bases of the adaptation to extreme environments and ecological shifts, as well of the origin and evolution of relevant spider traits, such as the venom and silk

Effect of the methacrylate-based endodontic sealer Epiphany on rat peritoneal macrophages viability

Objective: To evaluate the effect of the endodontic sealer Epiphany on rat peritoneal macrophages viability. Materials and methods: Peritoneal macrophages were obtained from Wistar rats and resuspended in RPMI- 1640 medium. Undiluted (crude extract) and diluted extracts to 10%, 1%, 0.1%, 0.01%, 0.001% and 0.0001% of Epiphany, AH 26 and AH Plus sealers on RPMI-1640 medium were tested for cytotoxicity to rat peritoneal macrophages using the trypan blue dye exclusion assay. Data were analyzed statistically by the Kruskal-Wallis and Mann-Whitney tests at 5% significance level. Results: Crude extract of Epiphany killed 51% of cells, but was less cytotoxic that crude extracts of AH Plus and AH 26, which killed 81% and 86% of cells, respectively. Ten-fold dilutions of Epiphany, AH Plus and AH 26 killed 44%, 56%, 62% of macrophages, respectively. A hundred dilution of Epiphany only killed 7% of macrophages, but the same dilution of AH Plus and AH 26 killed 10% and 31% of macrophages, respectively. Lower dilutions of sealer extracts caused minimal cell death as compared to the control groups (p>0.05). Conclusions: The methacrylate-based endodontic sealer Epiphany showed lower cytotoxicity on macrophages than resin-based sealers AH Plus and AH 26. Dilution of elutes of the three materials by tenfold markedly reduced their effect

Randomised multicentre clinical trial to evaluate voriconazole pre-emptive genotyping strategy in patients with risk of aspergillosis: vorigenipharm study protocol.

Introduction Invasive aspergillosis is the most important cause of morbidity and mortality in patients with haematological diseases. At present, voriconazole is the first-line treatment for invasive fungal disease. The pharmacokinetic interindividual variability of voriconazole depends on genetic factors. CYP450 is involved in 70%–75% of total metabolism of voriconazole, mainly CYP3A4 and CYP2C19, with the remaining 25%–30% of metabolism conducted by monooxygenase flavins. CYP2C19 single nucleotide polymorphisms could explain 50%–55% of variability in voriconazole metabolism. Materials and methods The main objective is to compare efficiency of pre-emptive voriconazole genotyping with routine practice. The primary outcome is serum voriconazole on the fifth day within the therapeutic range. The secondary outcome is the combined variables of therapeutic failure and adverse events within 90 days of first administration, associated with voriconazole. A total of 146 patients at risk of invasive aspergillosis who will potentially receive voriconazole will be recruited, and CYP2C19 will be genotyped. If the patient ultimately receives voriconazole, they will be randomised (1:1 experimental/control). In the experimental arm, patients will receive a dose according to a pharmacogenetic algorithm, including CYP2C19 genotype and clinical and demographic information. In the control arm, patients will receive a dose according to clinical practice guidelines. In addition, a Spanish National Healthcare System (NHS) point-of-view cost-effectiveness evaluation will be performed. Direct cost calculations for each arm will be performed. Conclusion This trial will provide information about the viability and cost-effectiveness of the mplementation of a pre-emptive voriconazole genotyping strategy in the Spanish NHS. Ethics and dissemination A Spanish version of this protocol has been evaluated and approved by the La Paz University Hospital Ethics Committee and the Spanish Agency of Medicines and Medical Devices. Trial results will be submitted for publication in an open peer-reviewed medical speciality-specific publication. Trial registration number Eudra-CT: 2019-000376-41 and NCT04238884; Pre-results.post-print441 K

Local wind speed forecasting based on WRF-HDWind coupling

[EN] Wind speed forecasts obtained by Numerical Weather Prediction models are limited for fine interpretation in heterogeneous terrain, in which different roughnesses and orographies occur. This limitation is derived from the use of low-resolution and grid-box averaged data. In this paper a dynamical downscaling method is presented to increase the local accuracy of wind speed forecasts. The proposed method divides the wind speed forecasting into two steps. In the first one, the mesoscale model WRF (Weather Research and Forecasting) is used for getting wind speed forecasts at specific points of the study domain. On a second stage, these values are used for feeding the HDWind microscale model. HDWind is a local model that provides both a high-resolution wind field that covers the entire study domain and values of wind speed and direction at very located points. As an example of use of the proposed method, we calculate a high-resolution wind field in an urban-interface area from Badajoz, a South-West Spanish city located near the Portugal border. The results obtained are compared with the values read by a weathervane tower of the Spanish State Meteorological Agency (AEMET) in order to prove that the microscale model improves the forecasts obtained by the mesoscale model

Deep Phenotyping and Genetic Characterization of a Cohort of 70 Individuals With 5p Minus Syndrome.

Chromosome-5p minus syndrome (5p-Sd, OMIM #123450) formerly known as Cri du Chat syndrome results from the loss of genetic material at the distal region of the short arm of chromosome 5. It is a neurodevelopmental disorder of genetic cause. So far, about 400 patients have been reported worldwide. Individuals affected by this syndrome have large phenotypic heterogeneity. However, a specific phenotype has emerged including global developmental delay, microcephaly, delayed speech, some dysmorphic features, and a characteristic and monochromatic high-pitch voice, resembling a cat's cry. We here describe a cohort of 70 patients with clinical features of 5p- Sd characterized by means of deep phenotyping, SNP arrays, and other genetic approaches. Individuals have a great clinical and molecular heterogeneity, which can be partially explained by the existence of additional significant genomic rearrangements in around 39% of cases. Thus, our data showed significant statistical differences between subpopulations (simple 5p deletions versus 5p deletions plus additional rearrangements) of the cohort. We also determined significant "functional" differences between male and female individuals.S