Search

Climate change challenges the adaptive capacity of several forest tree species in the face of increasing drought and rising temperatures. Therefore, understanding the mechanistic connections between genetic diversity and drought resilience is highly valuable for conserving drought-sensitive forests. Nonetheless, the post-drought recovery in trees from a transcriptomic perspective has not yet been studied by comparing contrasting phenotypes. Here, experimental drought treatments, gas-exchange dynamics and transcriptomic analysis (RNA-seq) were performed in the relict and drought-sensitive fir Abies pinsapo Boiss. to identify gene expression differences over immediate (24 h) and extended drought (20 days). Post-drought responses were investigated to define resilient and sensitive phenotypes. Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of A. pinsapo drought resilience. Weighted gene co-expression network analysis showed an activation of stomatal closing and an inhibition of plant growth-related genes during the immediate drought, consistent with an isohydric dynamic. During the extended drought, transcription factors, as well as cellular damage and homeostasis protection-related genes prevailed. Resilient individuals activate photosynthesis-related genes and inhibit aerial growth-related genes, suggesting a shifting shoot/root biomass allocation to improve water uptake and whole-plant carbon balance. About, 152 fixed SNPs were found between resilient and sensitive seedlings, which were mostly located in RNA-activity-related genes, including epigenetic regulation. Contrasting gene expression and SNPs were found between different post-drought resilience phenotypes for the first time in a forest tree, suggesting a transcriptomic and genomic basis for drought resilience. The obtained drought-related transcriptomic profile and drought-resilience candidate genes may guide conservation programs for this threatened tree species.

Data de publicació electrònica: 01-01-2023 Includes supplementary materials for the online appendix. Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis. This study was supported by the Spanish Ministry of Economy and Competitiveness, FEDER Grants (RTI2018-096847-B-C21, RTI2018-096847-B-C22 and RTI2018-097110-B-C21), Junta de Andalucía-FEDER Grant (P20-00938) and the European Union's Horizon 2020 research and innovation programme under grant agreement No. 817923 (AQUA-FAANG). We thank Geneaqua SL for their participation and financial support of sequencing. We acknowledge the bioinformatic support of the Centro de Supercomputación de Galicia (CESGA).

The resolution of cryptic diversity is essential for understanding the evolutionary diversification of lineages and establishing conservation priorities. We examine relationships in Eriogoneae (Polygonaceae), a diverse lineage in western North America. We ask whether Eriogonum umbellatum, a morphologically and ecologically diverse species, is monophyletic and whether its varieties represent evolutionary lineages. We use genotyping-by-sequencing to assemble a SNP dataset for 51 species in the genera Chorizanthe, Eriogonum and Sidotheca. We report a hierarchical phylogenetic analysis using maximum likelihood to estimate the evolutionary history of Eriogoneae. We illustrate admixture components for 21 populations of E. umbellatum, representing four varieties, and test for lineage structure using TreeMix. We identify strongly supported clades within Eriogoneae. Many relationships in the Eucycla + Oregonium and Latifolia clades are supported, while most relationships within the Eriogonum subg. Oligogonum clade and a clade with most Chorizanthe remain unresolved. Eriogonum congdonii resolves within the main E. umbellatum clade, while populations of three varieties of E. umbellatum are closely related to E. ursinum and are associated with serpentine soils. ADmixture and TreeMix analyses suggest E. umbellatum varieties represent evolutionary lineages. These results from SNP data are largely consistent with previous phylogenetic studies of Eriogoneae based on sequence variation. Structure within Oligogonum suggests consistent environmental association and radiation after initial colonization of serpentine. Morphology is unreliable for the infraspecific taxonomy of E. umbellatum. Additional molecular studies are needed to resolve the evolutionary relationships and ecological diversification within this species, in Oligogonum, and in Eriogoneae. We thank SGIker research support services at the University of the Basque Country, Leioa, Spain for DNA extraction and quality control, and acknowledge Centro Nacional de d-Analisi Genomica in Barcelona, Spain for GBS sequencing. We thank J. Andre, N.J. Jensen and J. Steele for contributing samples and D. Marino for collaboration in nucleic acid extractions. We acknowledge the contributions of O. Lao Grueso to the ADmixture and TreeMix analyses. We thank U.S. Forest Service employees D. Austin, J. Fedorchuk, M. Friend, J. Haas, D. Ikeda, L. Janeway, J. Nelson, D. Netz, A. Sanger, and S. Weis among others for facilitating permitting and collecting. This work was supported by funds from the Basque Government in support of the Terrestrial Plant Diversity group of the Department of Plant Biology and Ecology, University of the Basque Country, and an ERC Advanced Grant, FP7-IDEAS-ERC, `ADAPT', project 339941 awarded to T. Brown.

Genomic resources for amphibians are still hugely under-represented in vertebrate genomic research, despite being a group of major interest for ecology, evolution and conservation. Amphibians constitute a highly threatened group of vertebrates, present a vast diversity in reproductive modes, are extremely diverse in morphology, occupy most ecoregions of the world, and present the widest range in genome sizes of any major group of vertebrates. We combined Illumina, Nanopore and Hi-C sequencing technologies to assemble a chromosome-level genome sequence for an anuran with a moderate genome size (assembly span 3.09 Gb); Pelobates cultripes, the western spadefoot toad. The genome has an N50 length of 330 Mb with 98.6% of the total sequence length assembled into 14 super scaffolds, and 87.7% complete BUSCO genes. We use published transcriptomic data to provide annotations, identifying 32,684 protein-coding genes. We also reconstruct the P. cultripes phylome and identify 2,527 gene expansions. We contribute the first draft of the genome of the western spadefoot toad, P. cultripes. This species represents a relatively basal lineage in the anuran tree with an interesting ecology and a high degree of developmental plasticity, and thus is an important resource for amphibian genomic research. This research was funded by Plan Nacional I + D grant #CGL2017-83407-P awarded to I.G.-M. and a Ministerio de Ciencia e Innovación fellowship #IJC2018-035780-I awarded to H.C.L. TG’s group acknowledges support from the Spanish Ministry of Science and Innovation for grant PGC2018-099921-B-I00, cofounded by European Regional Development Fund (ERDF); and from the European Union’s Horizon 2020 research and innovation programme (ERC-2016-724173).

The integration of physical and high-density genetic maps is a very useful approach to achieve chromosome-level genome assemblies. Here, the genome of a male Senegalese sole (Solea senegalensis) was de novo assembled and the contigs were anchored to a high-quality genetic map for chromosome-level scaffolding. Hybrid assembled genome was 609.3 Mb long and contained 3403 contigs with a N50 of 513 kb. The linkage map was constructed using 16,287 informative SNPs derived from ddRAD sequencing in 327 sole individuals from five families. Markers were assigned to 21 linkage groups with an average number of 21.9 markers per megabase. The anchoring of the physical to the genetic map positioned 1563 contigs into 21 pseudo-chromosomes covering 548.6 Mb. Comparison of genetic and physical distances indicated that the average genome-wide recombination rate was 0.23 cM/Mb and the female-to-male ratio 1.49 (female map length: 2,698.4 cM, male: 2,036.6 cM). Genomic recombination landscapes were different between sexes with crossovers mainly concentrated toward the telomeres in males while they were more uniformly distributed in females. A GWAS analysis using seven families identified 30 significant sex-associated SNP markers located in linkage group 18. The follicle-stimulating hormone receptor appeared as the most promising locus associated with sex within a region with very low recombination rates. An incomplete penetrance of sex markers with males as the heterogametic sex was determined. An interspecific comparison with other Pleuronectiformes genomes identified a high sequence similarity between homologous chromosomes, and several chromosomal rearrangements including a lineage-specific Robertsonian fusion in S. senegalensis. This study was funded by project RTA2017-00054-C03-01 and RTA2017-00054-C03-funded from MCIU/AEI/FEDER, UE and cofunded 80% by Programa Operativo FEDER de Andalucía 2014-2020, project PP.AVA.AVA201601.9 SOLEALGAE. Moreover, the study has received funding from EU H2020 research and innovation program under grant agreement 817992 ERANET-BLUEBIO COFUND project PCI2020-111994 BestBrood/AEI/10.13039/501100011033. IGC is funded by a predoctoral fellowship from INIA. This work would not have been possible without the computer resources and the technical support provided by the Plataforma Andaluza de Bioinformática of the University of Málaga and CNAG.

Transcriptional quiescence of post-meiotic spermatozoa associated with protamine-mediated chromatin condensation is widely recognized in animals. How sperm acquire the extratesticular maturational competence to move and fertilize the egg is therefore thought to occur via non-transcriptional mechanisms. Here, using transcriptional profiling during spermatozoon differentiation in a fish that does not condense chromatin with protamines, we uncover spatially distinct roles of the GnRH receptor and PDGF signaling pathways between the somatic epithelia of the extratesticular ducts and the maturing spermatozoa. In vitro induction and inhibition experiments demonstrate that the endocrine signaling pathways are conserved in different lineages of fish and activate de novo transcription of spermatozoon genes required for the acquisition of full motility. These experiments further confirmed that mitochondrial translation is important for sperm maturation in anamniotes as in amniotes, but that transcriptional quiescence of post-meiotic spermatozoa is not a pan vertebrate phenomenon. On the contrary, the data show that the identified signal transduction pathways between the soma and the sperm upregulate effector genes essential for maturational competence and male fertility.

Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) has recently been used for the direct detection of KPC-producing isolates by analysis of the 11,109 Da mass peak representing the P019 protein. In this study, we evaluate the presence of the 11,109 Da mass peak in a collection of 435 unduplicated Klebsiella pneumoniae clinical isolates. The prevalence of the P019 peak in the blaKPC K. pneumoniae isolates was 49.2% (32/65). The 11,109 Da mass peak was not observed in any of the other carbapenemase (319) or noncarbapenemase producers (116). Computational analysis of the presence of the p019 gene was performed in the aforementioned carbapenemase-producing K. pneumoniae isolates fully characterized by whole-genome sequencing (WGS) and in a further collection of 1,649 K. pneumoniae genomes included in EuSCAPE. Herein, we have demonstrated that the p019 gene is not exclusively linked to the pKpQil plasmid but that it is present in the following plasmids: IncFIB(K)/IncFII(K)/ColRNAI, IncFIB(pQil), IncFIB(pQil)/ColRNAI, IncFIB(pQil)/IncFII(K), IncFIB(K)/IncFII(K), and IncX3. In addition, we have proven the independent movement of the Tn4401 and the ISKpn31, of which the p019 gene is a component. The absence of the p019 gene was obvious in Col440I, Col(pHAD28), IncFIB(K)/IncX3/IncFII(K), and IncFIB(K)/IncFII(K) plasmids. In addition, we also observed another plasmid in which neither Tn4401 nor ISKpn31 was found, IncP6. In the EuSCAPE, the occurrence of p019 varied from 0% to 100% among the different geographical locations. The adverse clinical impact of the diminished prevalence of the p019 gene within the plasmid encoding KPC-producing Klebsiella pneumoniae puts forward the need for reconsideration when applying this technique in a clinical setting.

The Eastern woodchuck (Marmota monax) has been extensively used in research of chronic hepatitis B and liver cancer because its infection with the woodchuck hepatitis virus closely resembles a human hepatitis B virus infection. Development of novel immunotherapeutic approaches requires genetic information on immune pathway genes in this animal model. The woodchuck genome was assembled with a combination of high-coverage whole-genome shotgun sequencing of Illumina paired-end, mate-pair libraries and fosmid pool sequencing. The result is a 2.63 Gigabase (Gb) assembly with a contig N50 of 74.5 kilobases (kb), scaffold N50 of 892 kb, and genome completeness of 99.2%. RNA sequencing (RNA-seq) from seven different tissues aided in the annotation of 30,873 protein-coding genes, which in turn encode 41,826 unique protein products. More than 90% of the genes have been functionally annotated, with 82% of them containing open reading frames. This genome sequence and its annotation will enable further research in chronic hepatitis B and hepatocellular carcinoma and contribute to the understanding of immunological responses in the woodchuck.

Aphids (Aphidoidea) are a diverse group of hemipteran insects that feed on plant phloem sap. A common finding in studies of aphid genomes is the presence of a large number of duplicated genes. However, when these duplications occurred remains unclear, partly due to the high relatedness of sequenced species. To better understand the origin of aphid duplications we sequenced and assembled the genome of Cinara cedri, an early branching lineage (Lachninae) of the Aphididae family. We performed a phylogenomic comparison of this genome with 20 other sequenced genomes, including the available genomes of five other aphids, along with the transcriptomes of two species belonging to Adelgidae (a closely related clade to the aphids) and Coccoidea. We found that gene duplication has been pervasive throughout the evolution of aphids, including many parallel waves of recent, species-specific duplications. Most notably, we identified a consistent set of very ancestral duplications, originating from a large-scale gene duplication predating the diversification of Aphidomorpha (comprising aphids, phylloxerids, and adelgids). Genes duplicated in this ancestral wave are enriched in functions related to traits shared by Aphidomorpha, such as association with endosymbionts, and adaptation to plant defenses and phloem-sap-based diet. The ancestral nature of this duplication wave (106–227 Ma) and the lack of sufficiently conserved synteny make it difficult to conclude whether it originated from a whole-genome duplication event or, alternatively, from a burst of large-scale segmental duplications. Genome sequencing of other aphid species belonging to different Aphidomorpha and related lineages may clarify these findings., This research was funded by European Regional Development Fund (ERDF) and Ministerio de Economía y Competitividad (Spain) (Grant Nos. PGC2018-099344-B-100 and BFU2015-67107). T.G. group also acknowledges support from the Catalan Research Agency (AGAUR) SGR857, and grants from the European Union’s Horizon 2020 research and innovation program under the grant agreements ERC-2016-724173 and MSC-747607. T.G. also receives support from an INB (Grant No. PT17/0009/0023—ISCIII-SGEFI/ERDF).

Amphibians are a highly vulnerable and highly diverse group of vertebrates for which we still have modest genomic resources. In particular, amphibian larvae are a key component of continental wetlands, where they have strong influences on energy fluxes, nutrient cycling, and community structure. Amphibian larvae are highly sensitive to environmental conditions and can often alter their physiology, behaviour and even morphology in response to the local conditions experienced, although we still know relatively little about the transcriptomic changes that enable such plasticity. Here we contribute the transcriptomes of two spadefoot toad species with divergent degrees of developmental plasticity in response to one of the most common hazards to amphibians breeding in temporary ponds, risk of pond drying.

The increasing emergence of carbapenemase-producing Klebsiella pneumoniae (CPK) isolates is a global health alarm. Rapid methods that require minimum sample preparation and rapid data analysis are urgently required. Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) has recently been used by clinical laboratories for identification of antibiotic-resistant bacteria; however, discrepancies have arisen regarding biological and technical issues. The aim of this study was to standardize an operating procedure and data analysis for identification of CPK by MALDI-TOF MS. To evaluate this approach, a series of 162 K. pneumoniae isolates (112 CPK and 50 non-CPK) were processed in the MALDI BioTyper system (Bruker Daltonik, Germany) following a standard operating procedure. The study was conducted in two stages; the first is denominated the "reproducibility stage" and the second "CPK identification." The first stage was designed to evaluate the biological and technical variation associated with the entire analysis of CPK and the second stage to assess the final accuracy of MALDI-TOF MS for the identification of CPK. Therefore, we present an improved MALDI-TOF MS data analysis pipeline using neural network analysis implemented in Clover MS Data Analysis Software (Clover Biosoft, Spain) that is designed to reduce variability, guarantee interlaboratory reproducibility, and maximize the information selected from the bacterial proteome. Using the random forest (RF) algorithm, 100% of CPK isolates were correctly identified when all the peaks in the spectra were selected as input features and total ion current (TIC) normalization was applied. Thus, we have demonstrated that real-time direct tracking of CPK is possible using MALDI-TOF MS.

Matrix-assisted laser desorption ionization–time of flight mass spectrometry (MALDI-TOF MS) has recently been used for the direct detection of KPC-producing isolates by analysis of the 11,109 Da mass peak representing the P019 protein. In this study, we evaluate the presence of the 11,109 Da mass peak in a collection of 435 unduplicated Klebsiella pneumoniae clinical isolates. The prevalence of the P019 peak in the bla(KPC) K. pneumoniae isolates was 49.2% (32/65). The 11,109 Da mass peak was not observed in any of the other carbapenemase (319) or noncarbapenemase producers (116). Computational analysis of the presence of the p019 gene was performed in the aforementioned carbapenemase-producing K. pneumoniae isolates fully characterized by whole-genome sequencing (WGS) and in a further collection of 1,649 K. pneumoniae genomes included in EuSCAPE. Herein, we have demonstrated that the p019 gene is not exclusively linked to the pKpQil plasmid but that it is present in the following plasmids: IncFIB(K)/IncFII(K)/ColRNAI, IncFIB(pQil), IncFIB(pQil)/ColRNAI, IncFIB(pQil)/IncFII(K), IncFIB(K)/IncFII(K), and IncX3. In addition, we have proven the independent movement of the Tn4401 and the ISKpn31, of which the p019 gene is a component. The absence of the p019 gene was obvious in Col440I, Col(pHAD28), IncFIB(K)/IncX3/IncFII(K), and IncFIB(K)/IncFII(K) plasmids. In addition, we also observed another plasmid in which neither Tn4401 nor ISKpn31 was found, IncP6. In the EuSCAPE, the occurrence of p019 varied from 0% to 100% among the different geographical locations. The adverse clinical impact of the diminished prevalence of the p019 gene within the plasmid encoding KPC-producing Klebsiella pneumoniae puts forward the need for reconsideration when applying this technique in a clinical setting.

Ostreococcus tauri is a simple unicellular green alga representing an ecologically important group of phytoplankton in oceans worldwide. Modern molecular techniques must be developed in order to understand the mechanisms that permit adaptation of microalgae to their environment. We present for the first time in O. tauri a detailed characterization of individual genomic integration events of foreign DNA of plasmid origin after PEG-mediated transformation. Vector integration occurred randomly at a single locus in the genome and mainly as a single copy. Thus, we confirmed the utility of this technique for insertional mutagenesis. While the mechanism of double-stranded DNA repair in the O. tauri model remains to be elucidated, we clearly demonstrate by genome resequencing that the integration of the vector leads to frequent structural variations (deletions/insertions and duplications) and some chromosomal rearrangements in the genome at the insertion loci. Furthermore, we often observed variations in the vector sequence itself. From these observations, we speculate that a nonhomologous end-joining-like mechanism is employed during random insertion events, as described in plants and other freshwater algal models. PEG-mediated transformation is therefore a promising molecular biology tool, not only for functional genomic studies, but also for biotechnological research in this ecologically important marine alga. This research received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 824110—EASI-Genomics and the French National Research Agency project grant (ANR Algalvirus ANR-17-CE02-0012). We acknowledge the 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 cofinancing 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

Rift Valley fever phlebovirus (RVFV) causes an emerging zoonotic disease and is mainly transmitted by Culex and Aedes mosquitoes. While Aedes aegypti-dengue virus (DENV) is the most studied model, less is known about the genes involved in infection-responses in other mosquito-arboviruses pairing. The main objective was to investigate the molecular responses of Cx. pipiens to RVFV exposure focusing mainly on genes implicated in innate immune responses. Mosquitoes were fed with blood spiked with RVFV. The fully-engorged females were pooled at 3 different time points: 2 hours post-exposure (hpe), 3- and 14-days post-exposure (dpe). Pools of mosquitoes fed with non-infected blood were also collected for comparisons. Total RNA from each mosquito pool was subjected to RNA-seq analysis and a de novo transcriptome was constructed. A total of 451 differentially expressed genes (DEG) were identified. Most of the transcriptomic alterations were found at an early infection stage after RVFV exposure. Forty-eight DEG related to immune infection-response were characterized. Most of them were related with the RNAi system, Toll and IMD pathways, ubiquitination pathway and apoptosis. Our findings provide for the first time a comprehensive view on Cx. pipiens-RVFV interactions at the molecular level. The early depletion of RNAi pathway genes at the onset of the RVFV infection would allow viral replication in mosquitoes. While genes from the Toll and IMD immune pathways were altered in response to RVFV none of the DEG were related to the JAK/STAT pathway. The fact that most of the DEG involved in the Ubiquitin-proteasome pathway (UPP) or apoptosis were found at an early stage of infection would suggest that apoptosis plays a regulatory role in infected Cx. pipiens midguts. This study provides a number of target genes that could be used to identify new molecular targets for vector control., This study was funded by CERCA Programme/Generalitat de Catalunya (http://cerca. cat/) and NB received an award to perform this work from Ministerio de Economı´a, Industria y Competitividad, Gobierno de España (grant no. MINECO AGL2013-47257-P (http://www.mineco. gob.es/portal/site/mineco/).

Background Olive tree (Olea europaea L. subsp. europaea, Oleaceae) has been the most emblematic perennial crop for Mediterranean countries since its domestication around 6000 years ago in the Levant. Two taxonomic varieties are currently recognized: cultivated (var. europaea) and wild (var. sylvestris) trees. However, it remains unclear whether olive cultivars derive from a single initial domestication event followed by secondary diversification, or whether cultivated lineages are the result of more than a single, independent primary domestication event. To shed light into the recent evolution and domestication of the olive tree, here we analyze a group of newly sequenced and available genomes using a phylogenomics and population genomics framework. Results We improved the assembly and annotation of the reference genome, newly sequenced the genomes of twelve individuals: ten var. europaea, one var. sylvestris, and one outgroup taxon (subsp. cuspidata)—and assembled a dataset comprising whole genome data from 46 var. europaea and 10 var. sylvestris. Phylogenomic and population structure analyses support a continuous process of olive tree domestication, involving a major domestication event, followed by recurrent independent genetic admixture events with wild populations across the Mediterranean Basin. Cultivated olives exhibit only slightly lower levels of genetic diversity than wild forms, which can be partially explained by the occurrence of a mild population bottleneck 3000–14,000 years ago during the primary domestication period, followed by recurrent introgression from wild populations. Genes associated with stress response and developmental processes were positively selected in cultivars, but we did not find evidence that genes involved in fruit size or oil content were under positive selection. This suggests that complex selective processes other than directional selection of a few genes are in place. Conclusions Altogether, our results suggest that a primary domestication area in the eastern Mediterranean basin was followed by numerous secondary events across most countries of southern Europe and northern Africa, often involving genetic admixture with genetically rich wild populations, particularly from the western Mediterranean Basin.

BackgroundThe olive tree (Olea europaea L. subsp. europaea, Oleaceae) has been the most economic perennial crop for Mediterranean countries since its domestication around 6,000 years ago. Two taxonomic varieties are currently recognized: cultivated (var. europaea) and wild (var. sylvestris) trees. To shed light into the recent evolution and domestication of the olive tree, we sequenced the genomes of twelve individuals: ten var. europaea, one var. sylvestris, and one outgroup taxon (subsp. cuspidata). All of them were analysed together with an improved assembly of var. europaea reference genome and the available assembly of var. sylvestris.ResultsOur analyses show that cultivated olives exhibit slightly lower levels of overall genetic diversity than wild forms, and that this can be partially explained by the occurrence of a mild population bottleneck 5000-7000 years ago during the primary domestication period. We also provide the first phylogenetic analysis of genome-wide sequences, which supports a continuous process of domestication of the olive tree. This, together with population structure and introgression analyses highlights genetic admixture with wild populations across the Mediterranean Basin in the course of domestication.ConclusionsAltogether, our results suggest that a primary domestication area in the eastern Mediterranean basin was followed by numerous secondary events across most countries of southern Europe and northern Africa, often involving genetic admixture with genetically rich wild populations, particularly from the western Mediterranean Basin. Based on selection tests and a search for selective sweeps, we found that genes associated with stress response and developmental processes were positively selected in cultivars. However, we did not find evidence that genes involved in fruit size or oil content were under positive selection.

21 pages, 4 figures, Background: The Mediterranean mussel Mytilus galloprovincialis is an ecologically and economically relevant edible marine bivalve, highly invasive and resilient to biotic and abiotic stressors causing recurrent massive mortalities in other bivalves. Although these traits have been recently linked with the maintenance of a high genetic variation within natural populations, the factors underlying the evolutionary success of this species remain unclear. Results: Here, after the assembly of a 1.28-Gb reference genome and the resequencing of 14 individuals from two independent populations, we reveal a complex pan-genomic architecture in M. galloprovincialis, with a core set of 45,000 genes plus a strikingly high number of dispensable genes (20,000) subject to presence-absence variation, which may be entirely missing in several individuals. We show that dispensable genes are associated with hemizygous genomic regions affected by structural variants, which overall account for nearly 580 Mb of DNA sequence not included in the reference genome assembly. As such, this is the first study to report the widespread occurrence of gene presence-absence variation at a whole-genome scale in the animal kingdom. Conclusions: Dispensable genes usually belong to young and recently expanded gene families enriched in survival functions, which might be the key to explain the resilience and invasiveness of this species. This unique pan-genome architecture is characterized by dispensable genes in accessory genomic regions that exceed by orders of magnitude those observed in other metazoans, including humans, and closely mirror the open pan-genomes found in prokaryotes and in a few nonmetazoan eukaryotes., This work was conducted with the support of the projects AGL2011-14507-E, AGL2015-65705-R, RTI2018-095997-B-I00 (Ministerio de Ciencia, Innovación y Universidades, Spain) and INCITE 10PXIB402096PR, IN607B 2016/12 (Consellería de Economía, Emprego e Industria - GAIN, Xunta de Galicia). Antonio Figueras, Beatriz Novoa, Rebeca Moreira, Alberto Pallavicini, Marco Gerdol, Paola Venier, and Umberto Rosani are supported by the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 678589. David Posada is supported by the European Research Council, the Spanish Ministry of Economy and Competitiveness, and Xunta de Galicia. We acknowledge the support of the Spanish Ministry of Science and Innovation to the EMBL partnership, the Centro de Excelencia Severo Ochoa, the CERCA Programme/Generalitat de Catalunya, the Spanish Ministry of Science and Innovation through the Instituto de Salud Carlos III, the Generalitat de Catalunya through Departament de Salut and Departament d’Empresa i Coneixement, and the Co-financing by the Spanish Ministry of Science and Innovation with funds from the European Regional Development Fund (ERDF) corresponding to the 2014-2020 Smart Growth Operating Program.

The evolution of winged insects revolutionized terrestrial ecosystems and led to the largest animal radiation on Earth. However, we still have an incomplete picture of the genomic changes that underlay this diversification. Mayflies, as one of the sister groups of all other winged insects, are key to understanding this radiation. Here, we describe the genome of the mayfly Cloeon dipterum and its gene expression throughout its aquatic and aerial life cycle and specific organs. We discover an expansion of odorant-binding-protein genes, some expressed specifically in breathing gills of aquatic nymphs, suggesting a novel sensory role for this organ. In contrast, flying adults use an enlarged opsin set in a sexually dimorphic manner, with some expressed only in males. Finally, we identify a set of wing-associated genes deeply conserved in the pterygote insects and find transcriptomic similarities between gills and wings, suggesting a common genetic program. Globally, this comprehensive genomic and transcriptomic study uncovers the genetic basis of key evolutionary adaptations in mayflies and winged insects., This project was mainly funded by the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie Grant Agreement 657732 to I.A., Grant BFU2015-66040-P to F.Ca., institutional Grant MDM-2016-0687 (MINECO, Spain). Additional funding was provided by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (ERC-StG-LS2-637591 to M.I.), the Spanish Ministerio de Ciencia (BFU2017-89201-P to M.I., RYC-2016-20089 and PGC2018-099392-A-I00 to I.M.).

Phytoplasmas are bacterial plant pathogens that are detrimental to many plants and cause devastating effects on crops. They are not viable outside their host plants and depend on specific insect vectors for their transmission. So far, research has largely focused on plant-pathogen interactions, while the complex interactions between phytoplasmas and insect vectors are far less understood. Here, we used next-generation sequencing to investigate how transcriptional profiles of the vector psyllid Cacopsylla melanoneura (Hemiptera, Psyllidae) are altered during infection by the bacterium Candidatus Phytoplasma mali (P. mali), which causes the economically important apple proliferation disease. This first de novo transcriptome assembly of an apple proliferation vector revealed that mainly genes involved in small GTPase mediated signal transduction, nervous system development, adhesion, reproduction, actin-filament based and rhythmic processes are significantly altered upon P. mali infection. Furthermore, the presence of P. mali is accompanied by significant changes in carbohydrate and polyol levels, as revealed by metabolomics analysis. Taken together, our results suggest that infection with P. mali impacts on the insect vector physiology, which in turn likely affects the ability of the vector to transmit phytoplasma.

12 pages, 5 figures, 6 tables, supplemental material https://doi.org/10.25387/g3.12640691, The octocoral, Paramuricea clavata, is a habitat-forming anthozoan with a key ecological role in rocky benthic and biodiversity-rich communities in the Mediterranean and Eastern Atlantic. Shallow populations of P. clavata in the North-Western Mediterranean are severely affected by warming-induced mass mortality events (MMEs). These MMEs have differentially impacted individuals and populations of P. clavata (i.e., varied levels of tissue necrosis and mortality rates) over thousands of kilometers of coastal areas. The eco-evolutionary processes, including genetic factors, contributing to these differential responses remain to be characterized. Here, we sequenced a P. clavata individual with short and long read technologies, producing 169.98 Gb of Illumina paired-end and 3.55 Gb of Oxford Nanopore Technologies (ONT) reads. We obtained a de novo genome assembly accounting for 607 Mb in 64,145 scaffolds. The contig and scaffold N50s are 19.15 Kb and 23.92 Kb, respectively. Despite of the low contiguity of the assembly, its gene completeness is relatively high, including 75.8% complete and 9.4% fragmented genes out of the 978 metazoan genes contained in the metazoa_odb9 database. A total of 62,652 protein-coding genes have been annotated. This assembly is one of the few octocoral genomes currently available. This is undoubtedly a valuable resource for characterizing the genetic bases of the differential responses to thermal stress and for the identification of thermo-resistant individuals and populations. Overall, having the genome of P. clavata will facilitate studies of various aspects of its evolutionary ecology and elaboration of effective conservation plans such as active restoration to overcome the threats of global change, We acknowledge the funding support of the European Union’s Horizon 2020 research and innovation program under grant agreement No 689518 (MERCES) and the Strategic Funding UID/Multi/04423/2013 through national funds provided by FCT – Foundation for Science and Technology and European Regional Development Fund (ERDF), in the framework of the programme PT2020. JBL is funded by an assistant researcher contract framework of the RD Unit - UID/Multi/04423/2019 - Interdisciplinary Centre of Marine and Environmental Research – financed by the European Regional Development Fund (ERDF) through COMPETE2020 - Operational Program for Competitiveness and Internationalisation (POCI) and national funds through FCT/MCTES (PIDDAC). AA was partially supported by the FCT project PTDC/CTA-AMB/31774/2017 (POCI-01-0145-FEDER/031774/2017). DGG was supported by an FPU grant (FPU15/05457) from the Spanish Ministry of Education. We acknowledge support of the Spanish Ministry of Science and Innovation to the EMBL partnership, the Centro de Excelencia Severo Ochoa, the CERCA Programme / Generalitat de Catalunya, the Spanish Ministry of Science and Innovation through the Instituto de Salud Carlos III, the Generalitat de Catalunya through Departament de Salut and Departament d'Empresa i Coneixement, and the Co-financing by the Spanish Ministry of Science and Innovation with funds from the European Regional Development Fund (ERDF) corresponding to the 2014-2020 Smart Growth Operating Program. JBL and JG planned and granted the funding to start the project, With the funding support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S), of the Spanish Research Agency (AEI)

The first winged insects underwent profound morphological and functional transformations leading to the most successful animal radiations in the history of earth. Despite this, we still have a very incomplete picture of the changes in their genomes that underlay this radiation. Mayflies (Ephemeroptera) are one of the extant sister groups of all other winged insects and therefore are at a key phylogenetic position to understand this radiation. Here, we describe the genome of the cosmopolitan mayfly Cloeon dipterum and study its expression along development and in specific organs. We discover an expansion of odorant-binding proteins, some expressed specifically in the breathing gills of aquatic nymphs, suggesting a novel sensory role for gills. In contrast, as flying adults, mayflies make use of an enlarged set of opsins and utilise these visual genes in a sexually dimorphic manner, with some opsins expressed only in males. Finally, to illuminate the origin of wings, we identify a core set of deeply conserved wing-specific genes at the root of the pterygote insects. Globally, this is the first comprehensive study of the structure and expression of the genome of a paleopteran insect and shows how its genome has kept a record of its functional adaptations.

The octocoral, Paramuricea clavata, is a habitat-forming anthozoan with a key ecological role in rocky benthic and biodiversity-rich communities in the Mediterranean and Eastern Atlantic. Shallow populations of P. clavata in the North-Western Mediterranean are severely affected by warming-induced mass mortality events (MMEs). These MMEs have differentially impacted individuals and populations of P. clavata (i.e. varied levels of tissue necrosis and mortality rates) over thousands of kilometers of coastal areas. The eco-evolutionary processes and genetic factors contributing to these differential responses remain to be characterized. Here, we sequenced a P. clavata individual with short and long read technologies, producing 169.98 Gb of Illumina paired-end and 3.55 Gb of Oxford Nanopore Technologies (ONT) reads. We obtained a de novo hybrid assembly accounting for 712.4 Mb and 107,682 scaffolds. The contig and scaffold N50 are 15.85 Kb and 17.01 Kb, respectively. Despite of the low contiguity of the assembly, the gene completeness was relatively high, including 86% of the 978 metazoan genes contained in the metazoa_odb9 database. A total of 76,508 protein-coding genes and 85,763 transcripts have been annotated. This assembly is one of the few octocoral genomes currently available. This is undoubtedly a valuable resource for characterizing the genetic bases of the differential responses to thermal stress and for the identification of thermo-resistant individuals and populations. Overall, the genome of P. clavata will help to understand various aspects of its evolutionary ecology and to elaborate effective conservation plans such as active restoration actions to overcome the threats of global change.

Campylobacter is the causal agent of campylobacteriosis in humans, a self-limiting gastroenteritis. Campylobacteriosis is a zoonosis, commonly transmitted from contaminated chicken meat by either direct consumption or cross contamination during food manipulation. Presence of plasmids encoding for resistance to antibiotics such as tetracycline is common among Campylobacter isolates. In this report, we studied the effect of the temperature in the conjugation frequency of several tet(O) carrying plasmids, providing tetracycline resistance to the recipient cells. The conjugation frequency from donor cells carrying three previously characterized plasmids (pCjA13, pCjA9 and pTet) and from two clinical isolates was determined. Two temperatures, 37 and 42 degrees C, mimicking the conditions encountered by C. jejuni in the human and broiler chicken gastrointestinal tracts, respectively, were assessed. Our results clearly indicate that the conjugation process is promoted at high temperature. Accordingly, the transcriptional expression of some putative conjugative apparatus genes is thermoregulated, being induced at 42 degrees C. The two plasmids present in the clinical isolates were sequenced and assembled. Both plasmids are highly related among them and to the pTet plasmid. The high identity of the genes putatively involved in the conjugation process among the plasmids is in agreement with the similar behavior regarding the temperature dependency of the conjugative process. This report suggest that conjugation of plasmids carrying antibiotic resistance genes occurs preferentially at temperatures that resemble the gastrointestinal tract of birds, the main reservoir of C. jejuni.

Mussels are ecologically and economically relevant edible marine bivalves, highly invasive and resilient to biotic and abiotic stressors causing recurrent massive mortalities in other species. Here we show that the Mediterranean mussel Mytilus galloprovincialis has a complex pan-genomic architecture, which includes a core set of 45,000 genes shared by all individuals plus a surprisingly high number of dispensable genes (∼15,000). The latter are subject to presence/absence variation (PAV), i.e., they may be entirely missing in a given individual and, when present, they are frequently found as a single copy. The enrichment of dispensable genes in survival functions suggests an adaptive value for PAV, which might be the key to explain the extraordinary capabilities of adaptation and invasiveness of this species. Our study underpins a unique metazoan pan-genome architecture only previously described in prokaryotes and in a few non-metazoan eukaryotes, but that might also characterize other marine invertebrates.Significance statementIn animals, intraspecific genomic diversity is generally thought to derive from relatively small-scale variants, such as single nucleotide polymorphisms, small indels, duplications, inversions and translocations. On the other hand, large-scale structural variations which involve the loss of genomic regions encoding protein-coding genes in some individuals (i.e. presence/absence variation, PAV) have been so far only described in bacteria and, occasionally, in plants and fungi. Here we report the first evidence of a pan-genome in the animal kingdom, revealing that 25% of the genes of the Mediterranean mussel are subject to PAV. We show that this unique feature might have an adaptive value, due to the involvement of dispensable genes in functions related with defense and survival.

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated 238 Mb almond genome size, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27,969 protein-coding genes and 6,747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (P. persica) diverged around 5.88 Mya. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions/kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). TEs have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. TEs may also be at the origin of important phenotypic differences between both species, and in particular, for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach., This research was supported in part by grants from: the Ministry of Economy and Competitiveness (MINECO/FEDER projects AGL2012-40228-C02-01, AGL2015-68329-R and RTA2015-00050-00-00, Severo Ochoa Program for Centres of Excellence in R&D 201-2019 SEV-2015-0533 and CERCA Programme-Generalitat de Catalunya) from Spain. MW acknowledges grant support from the Waite Research Institute of the University of Adelaide.

Silver fir (Abies alba Mill.) is a keystone conifer of European montane forest ecosystems that has experienced large fluctuations in population size during during the Quaternary and, more recently, due to land-use change. To forecast the species’ future distribution and survival, it is important to investigate the genetic basis of adaptation to environmental change, notably to extreme events. For this purpose, we here provide a first draft genome assembly and annotation of the silver fir genome, established through a community-based initiative. DNA obtained from haploid megagametophyte and diploid needle tissue was used to construct and sequence Illumina paired-end and mate-pair libraries, respectively, to high depth. The assembled A. alba genome sequence accounted for over 37 million scaffolds corresponding to 18.16 Gb, with a scaffold N50 of 14,051 bp. Despite the fragmented nature of the assembly, a total of 50,757 full-length genes were functionally annotated in the nuclear genome. The chloroplast genome was also assembled into a single scaffold (120,908 bp) that shows a high collinearity with both the A. koreana and A. sibirica complete chloroplast genomes. This first genome assembly of silver fir is an important genomic resource that is now publicly available in support of a new generation of research. By genome-enabling this important conifer, this resource will open the gate for new research and more precise genetic monitoring of European silver fir forests.

Combining both short and long-read sequencing, we have estimated the almondPrunus dulciscv. Texas genome size in 235 Mbp and assembled 227.6 Mb of its sequence. The highly heterozygous compact genome of Texas comprises eight chromosomes, to which we have anchored over 91% of the assembly. We annotated 27,042 protein-coding genes and 6,800 non-coding transcripts. High levels of genetic variability were characterized after resequencing a collection of ten almond accessions. Phylogenomic comparison with the genomes of 16 other close and distant species allowed estimating that almond and peach diverged around 5.88 Mya. Comparison between peach and almond genomes confirmed the high synteny between these close relatives, but also revealed high numbers of presence-absence variants, many attributable to the movement of transposable elements (TEs). The number and distribution of TEs between peach and almond was similar, but the history of TE movement was distinct, with peach having a larger proportion of recent transpositions and almond preserving a higher level of polymorphism in the older TEs. When focusing on specific genes involved in key characters such as the bitter vs. sweet kernel taste and the formation of a fleshy mesocarp, we found that for one gene associated with the biosynthesis of amygdalin that confers the bitter kernel taste, several TEs were inserted in its vicinity only in sweet almond cultivars but not in bitter cultivars andPrunusbitter kernel relatives, includingP. webbii,P. mume, and other species like peach and cherry. TE insertions likely to produce affects in the expression of six more genes involved in the formation of the fleshy mesocarp were also identified. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond with respect to peach.

12 páginas, 5 figuras.-- Antonio Figueras ... et al.-- This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited, The turbot is a flatfish (Pleuronectiformes) with increasing commercial value, which has prompted active genomic research aimed at more efficient selection. Here we present the sequence and annotation of the turbot genome, which represents a milestone for both boosting breeding programmes and ascertaining the origin and diversification of flatfish. We compare the turbot genome with model fish genomes to investigate teleost chromosome evolution. We observe a conserved macrosyntenic pattern within Percomorpha and identify large syntenic blocks within the turbot genome related to the teleost genome duplication. We identify gene family expansions and positive selection of genes associated with vision and metabolism of membrane lipids, which suggests adaptation to demersal lifestyle and to cold temperatures, respectively. Our data indicate a quick evolution and diversification of flatfish to adapt to benthic life and provide clues for understanding their controversial origin. Moreover, we investigate the genomic architecture of growth, sex determination and disease resistance, key traits for understanding local adaptation and boosting turbot production, by mapping candidate genes and previously reported quantitative trait loci. The genomic architecture of these productive traits has allowed the identification of candidate genes and enriched pathways that may represent useful information for future marker-assisted selection in turbot, This work was funded by the Spanish Government: projects Consolider Ingenio: Aquagenomics (CSD2007-00002) and Metagenoma de la Península Ibérica (CSD2007-00005), Ministerio de Economía y Competitividad and European Regional Development Funds (AGL2012-35904), and Ministerio de Economía y Competitividad (AGL2014-51773 and AGL2014-57065-R); and Local Government Xunta de Galicia (GRC2014/010). P.P. and D.R. gratefully acknowledge the Spanish Ministerio de Educación for their FPU fellowships (AP2010-2408, AP2012-0254). Funding to pay the Open Access publication charges for this article was provided by the Ministerio de Economía y Competitividad (AGL2014-51773) and Xunta de Galicia (GRC2014/010)

This unit describes the usage of geneid, an efficient gene-finding program that allows for the analysis of large genomic sequences, including whole mammalian chromosomes. These sequences can be partially annotated, and geneid can be used to refine this initial annotation. Training geneid is relatively easy, and parameter configurations exist for a number of eukaryotic species. geneid produces output in a variety of standard formats. The results, thus, can be processed by a variety of software tools, including visualization programs. geneid software is in the public domain, and is undergoing constant development. It is easy to install and use. Exhaustive benchmark evaluations show that geneid compares favorably with other existing gene-finding tools. © 2018 by John WileySons, Inc.

The existence of dynamic cellular phenotypes in changing environmental conditions is of major interest for cell biologists who aim to understand the mechanism and sequence of regulation of gene expression. In the context of therapeutic protein production by Chinese Hamster Ovary (CHO) cells, a detailed temporal understanding of cell-line behavior and control is necessary to achieve a more predictable and reliable process performance. Of particular interest are data on dynamic, temporally resolved transcriptional regulation of genes in response to altered substrate availability and culture conditions. In this study, the gene transcription dynamics throughout a 9-day batch culture of CHO cells was examined by analyzing histone modifications and gene expression profiles in regular 12- and 24-hr intervals, respectively. Three levels of regulation were observed: (a) the presence or absence of DNA methylation in the promoter region provides an ON/OFF switch; (b) a temporally resolved correlation is observed between the presence of active transcription- and promoter-specific histone marks and the expression level of the respective genes; and (c) a major mechanism of gene regulation is identified by interaction of coding genes with long non-coding RNA (lncRNA), as observed in the regulation of the expression level of both neighboring coding/lnc gene pairs and of gene pairs where the lncRNA is able to form RNA-DNA-DNA triplexes. Such triplex-forming regions were predominantly found in the promoter or enhancer region of the targeted coding gene. Significantly, the coding genes with the highest degree of variation in expression during the batch culture are characterized by a larger number of possible triplex-forming interactions with differentially expressed lncRNAs. This indicates a specific role of lncRNA-triplexes in enabling rapid and large changes in transcription. A more comprehensive understanding of these regulatory mechanisms will provide an opportunity for new tools to control cellular behavior and to engineer enhanced phenotypes. Funding: IHL received support from FWF, the Austrian Science Fund, within the Doctoral Program “BioTechnology of Proteins” W1224, and HD from the “eCHO Systems” ITN PhD program funded by the Marie Sklodowska-Curie grant no. 642663. Additional grants were received from the National Bioinformatics Institute (INB), PRB2-ISCIII (PT13/0001/0044 to JG and AEC) and MINECO (PTA2014-09515 to MD)

Mammalian Y chromosomes are often neglected from genomic analysis. Due to their inherent assembly difficulties, high repeat content, and large ampliconic regions1, only a handful of species have their Y chromosome properly characterized. To date, just a single human reference quality Y chromosome, of European ancestry, is available due to a lack of accessible methodology2–5. To facilitate the assembly of such complicated genomic territory, we developed a novel strategy to sequence native, unamplified flow sorted DNA on a MinION nanopore sequencing device. Our approach yields a highly continuous and complete assembly of the first human Y chromosome of African origin. It constitutes a significant improvement over comparable previous methods, increasing continuity by more than 800%6, thus allowing a chromosome scale analysis of human Y chromosomes. Sequencing native DNA also allows to take advantage of the nanopore signal data to detect epigenetic modifications in situ7. This approach is in theory generalizable to any species simplifying the assembly of extremely large and repetitive genomes.

Viral diseases are responsible for high rates of mortality and subsequent economic losses in modern aquaculture. The nervous necrosis virus (NNV) produces viral encephalopathy and retinopathy (VER), which affects the fish central nervous system. It is considered one of the most serious viral diseases in marine aquaculture, the European sea bass (Dicentrarchus labrax) being amongst the most susceptible. We have evaluated the European sea bass brain derived cell line (DLB-1) susceptibility to NNV genotypes and evaluated its transcriptomic profile. DLB-1 cells supported NNV gene transcription and replication since strains belonging to the four NNV genotypes produce cytopathic effects. Afterwards, DLB-1 cells were infected with an RGNNV strain, the one which showed the highest replication, for 12 and 72 h and an RNA-seq analysis was performed to identify potential genes involved in the host-NNV interactions. Differential expression analysis showed the up-regulation of many genes related to immunity, heat-shock proteins or apoptosis but not to proteasome or autophagy processes. These data suggest that the immune response, mainly the interferon (IFN) pathway, is not powerful enough to abrogate the infection, and cells finally suffer stress and die by apoptosis liberating infective particles. GO enrichment also revealed, for the first time, the down-regulation of terms related to brain/neuron biology indicating molecular mechanisms causing the pathogenic effect of NNV. This study opens the way to understand key elements in sea bass brain and NNV interactions., NODAMED, Detección de nodavirus en el mar Mediterráneo. Implicación y riesgo para la acuicultura, SI

Nervous necrosis virus (NNV) causes high mortalities in several marine species. We aimed to evaluate the innate cell-mediated cytotoxic (CMC) activity of head-kidney leucocytes (HKLs) isolated from naïve European sea bass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata), a very susceptible and resistant fish species to NNV, respectively, against fish cell lines infected with NNV. Seabream HKLs showed significantly increased innate CMC activity against NNV-infected cells, compared to those uninfected, while sea bass HKLs failed to do so. Thus, we performed a RNA-seq study to identify genes related to the CMC activity of sea bass leucocytes. Thus, we found that sea bass HKLs incubated with DLB-1 cells alone (CMC_DLB1) or with NNV-infected DLB-1 cells (CMC_DLB1-NNV) showed very similar transcriptomic profiles and the GO analysis revealed that most of the up-regulated genes were related to immunity. Strikingly, when the CMC samples with and without NNV were compared, GO analysis revealed that most of the up-regulated genes in CMC_DLB1-NNV samples were related to metabolism and very few to immunity. This is also in agreement with the functional data. These data point to the escape of CMC activity by NNV infection as an important factor involved in the high susceptibility to nodavirus infections of European sea bass., SI

Background: Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes. Results: We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools. Conclusions: The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop., La lista completa de autores que integran el documento puede consultarse en el archivo., Facultad de Ciencias Exactas, Instituto de Biotecnologia y Biologia Molecular

12 p.-6 fig.-5 tab. Fernando Cruz et al., Background: The Mediterranean olive tree (Olea europaea subsp. europaea) was one of the first trees to be domesticated and is currently of major agricultural importance in the Mediterranean region as the source of olive oil. The molecular bases underlying the phenotypic differences among domesticated cultivars, or between domesticated olive trees and their wild relatives, remain poorly understood. Both wild and cultivated olive trees have 46 chromosomes (2n)., Findings: A total of 543 Gb of raw DNA sequence from whole genome shotgun sequencing, and a fosmid library containing 155,000 clones from a 1,000+ year-old olive tree (cv. Farga) were generated by Illumina sequencing using different combinations of mate-pair and pair-end libraries. Assembly gave a final genome with a scaffold N50 of 443 kb, and a total length of 1.31 Gb, which represents 95 % of the estimated genome length (1.38 Gb). In addition, the associated fungus Aureobasidium pullulans was partially sequenced. Genome annotation, assisted by RNA sequencing from leaf, root, and fruit tissues at various stages, resulted in 56,349 unique protein coding genes, suggesting recent genomic expansion. Genome completeness, as estimated using the CEGMA pipeline, reached 98.79 %., Conclusions: The assembled draft genome of O. europaea will provide a valuable resource for the study of the evolution and domestication processes of this important tree, and allow determination of the genetic bases of key phenotypic traits. Moreover, it will enhance breeding programs and the formation of new varieties., This project was funded by Banco Santander, which provided plant material and financed whole genome sequencing of the olive tree. The authors especially want to thank the late Mr. Emilio Botín for his support in driving this project. IJ was supported by a grant from the Peruvian Ministry of Education: ‘Beca Presidente de la República’ (2013-III).

Recent technological developments have revolutionized the way we perform genetic analyses. In particular whole-genome sequencing provides access to the entire genetic makeup of an individual, and it is now an affordable approach for many research groups. As a consequence genome sequencing is pervading many fields of biological research. Sequencing technologies are evolving rapidly and so do their applications. Here we provide a first primer on whole-genome sequencing, focusing on two of the most popular applications: (1) de novo genome sequencing, in which the objective is obtaining a high-quality genome assembly that can serve as a reference for a species or variety, and (2) genome resequencing, when there is an available reference genome and the objective is to map sequence variation of an individual or a set of individuals. It is not our intention to provide a comprehensive overview of current methodologies that will likely soon become obsolete, but rather focus on general principles that will have a more general applicability.

Alternative splicing (AS) has the potential to greatly expand the functional repertoire of mammalian transcriptomes. However, few variant transcripts have been characterized functionally, making it difficult to assess the contribution of AS to the generation of phenotypic complexity and to study the evolution of splicing patterns. We have compared the AS of 309 protein-coding genes in the human ENCODE pilot regions against their mouse orthologs in unprecedented detail, utilizing traditional transcriptomic and RNAseq data. The conservation status of every transcript has been investigated, and each functionally categorized as coding (separated into coding sequence [CDS] or nonsense-mediated decay [NMD] linked) or noncoding. In total, 36.7% of human and 19.3% of mouse coding transcripts are species specific, and we observe a 3.6 times excess of human NMD transcripts compared with mouse; in contrast to previous studies, the majority of species-specific AS is unlinked to transposable elements. We observe one conserved CDS variant and one conserved NMD variant per 2.3 and 11.4 genes, respectively. Subsequently, we identify and characterize equivalent AS patterns for 22.9% of these CDS or NMD-linked events in nonmammalian vertebrate genomes, and our data indicate that functional NMD-linked AS is more widespread and ancient than previously thought. Furthermore, although we observe an association between conserved AS and elevated sequence conservation, as previously reported, we emphasize that 30% of conserved AS exons display sequence conservation below the average score for constitutive exons. In conclusion, we demonstrate the value of detailed comparative annotation in generating a comprehensive set of AS transcripts, increasing our understanding of AS evolution in vertebrates. Our data supports a model whereby the acquisition of functional AS has occurred throughout vertebrate evolution and is considered alongside amino acid change as a key mechanism in gene evolution.

Septin 3 is a novel member of the septin subfamily of GTPase domain proteins that was recently identified in human neuronal cells. These proteins are involved in vesicle trafficking, neurite outgrowth, and neurofibrillary tangle formation; however, the expression and functional role of septin 3 in normal neuronal tissues and as an etiological agent in neurological disorders is currently unclear. To further characterize these parameters, the present study analyzed the expression of three isoforms of septin 3 (A, B, and C) in fetal and adult human brains and polymorphism of the septin 3 exon 11 microsatellite in control, pure Alzheimer's disease (AD), Lewy body variant (LBV) of AD, and Parkinson's disease. Septin 3 mRNAs for isoforms A and B, but not C, were detected in the frontal cortex of fetus and adult human samples, as measured by reverse transcription-coupled polymerase chain reaction. Genotype analyses indicated that polymorphic septin 3 alleles were distributed in two peaks of frequency in both control and disease groups. Categorization of the alleles into short (S) and long (L) types revealed a significant difference between AD patients and controls (p = 0.034 by chi-square test). Furthermore, the S-allele homozygosity was significantly underrepresented in AD compared with control (p = 0.015 by chi-square test). These results suggest that polymorphism in exon 11 of septin 3 may have a determinative role in the pathogenesis of AD.

The emergence of next generation DNA sequencing technology is enabling high-resolution cancer genome analysis. Large-scale projects like the International Cancer Genome Consortium (ICGC) are systematically scanning cancer genomes to identify recurrent somatic mutations. Second generation DNA sequencing, however, is still an evolving technology and procedures, both experimental and analytical, are constantly changing. Thus the research community is still defining a set of best practices for cancer genome data analysis, with no single protocol emerging to fulfil this role. Here we describe an extensive benchmark exercise to identify and resolve issues of somatic mutation calling. Whole genome sequence datasets comprising tumor-normal pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, were shared within the ICGC and submissions of somatic mutation calls were compared to verified mutations and to each other. Varying strategies to call mutations, incomplete awareness of sources of artefacts, and even lack of agreement on what constitutes an artefact or real mutation manifested in widely varying mutation call rates and somewhat low concordance among submissions. We conclude that somatic mutation calling remains an unsolved problem. However, we have identified many issues that are easy to remedy that are presented here. Our study highlights critical issues that need to be addressed before this valuable technology can be routinely used to inform clinical decision-making.