Combining Genetic and Transcriptomic Approaches to Identify Transporter-Coding Genes as Likely Responsible for a Repeatable Salt Tolerance QTL in Citrus.

The excessive accumulation of chloride (Cl-) in leaves due to salinity is frequently related to decreased yield in citrus. Two salt tolerance experiments to detect quantitative trait loci (QTLs) for leaf concentrations of Cl-, Na+, and other traits using the same reference progeny derived from the salt-tolerant Cleopatra mandarin (Citrus reshni) and the disease-resistant donor Poncirus trifoliata were performed with the aim to identify repeatable QTLs that regulate leaf Cl- (and/or Na+) exclusion across independent experiments in citrus, as well as potential candidate genes involved. A repeatable QTL controlling leaf Cl- was detected in chromosome 6 (LCl-6), where 23 potential candidate genes coding for transporters were identified using the C. clementina genome as reference. Transcriptomic analysis revealed two important candidate genes coding for a member of the nitrate transporter 1/peptide transporter family (NPF5.9) and a major facilitator superfamily (MFS) protein. Cell wall biosynthesis- and secondary metabolism-related processes appeared to play a significant role in differential gene expression in LCl-6. Six likely gene candidates were mapped in LCl-6, showing conserved synteny in C. reshni. In conclusion, markers to select beneficial Cleopatra mandarin alleles of likely candidate genes in LCl-6 to improve salt tolerance in citrus rootstock breeding programs are provided.

Insights into ROS-dependent signalling underlying transcriptomic plant responses to the herbicide 2,4-D.

The synthetic auxin 2,4-dichlorophenoxyacetic acid (2,4-D) functions as an agronomic weed control herbicide. High concentrations of 2,4-D induce plant growth defects, particularly leaf epinasty and stem curvature. Although the 2,4-D triggered reactive oxygen species (ROS) production, little is known about its signalling. In this study, by using a null mutant in peroxisomal acyl CoA oxidase 1 (acx1-2), we identified acyl-coenzyme A oxidase 1 (ACX1) as one of the main sources of ROS production and, in part, also causing the epinastic phenotype following 2,4-D application. Transcriptomic analyses of wild type (WT) plants after treatment with 2,4-D revealed a ROS-related peroxisomal footprint in early plant responses, while other organelles, such as mitochondria and chloroplasts, are involved in later responses. Interestingly, a group of 2,4-D-responsive ACX1-dependent transcripts previously associated with epinasty is related to auxin biosynthesis, metabolism, and signalling. We found that the auxin receptor auxin signalling F-box 3 (AFB3), a component of Skp, Cullin, F-box containing complex (SCF) (ASK-cullin-F-box) E3 ubiquitin ligase complexes, which mediates auxin/indole acetic acid (AUX/IAA) degradation by the 26S proteasome, acts downstream of ACX1 and is involved in the epinastic phenotype induced by 2,4-D. We also found that protein degradation associated with ubiquitin E3-RING and E3-SCF-FBOX in ACX1-dependent signalling in plant responses to 2,4-D is significantly regulated over longer treatment periods.

Study of microRNA expression in Salmonella Typhimurium-infected porcine ileum reveals miR-194a-5p as an important regulator of the TLR4-mediated inflammatory response.

Infection with Salmonella Typhimurium (S. Typhimurium) is a common cause of food-borne zoonosis leading to acute gastroenteritis in humans and pigs, causing economic losses to producers and farmers, and generating a food security risk. In a previous study, we demonstrated that S. Typhimurium infection produces a severe transcriptional activation of inflammatory processes in ileum. However, little is known regarding how microRNAs regulate this response during infection. Here, small RNA sequencing was used to identify 28 miRNAs differentially expressed (DE) in ileum of S. Typhimurium-infected pigs, which potentially regulate 14 target genes involved in immune system processes such as regulation of cytokine production, monocyte chemotaxis, or cellular response to interferon gamma. Using in vitro functional and gain/loss of function (mimics/CRISPR-Cas system) approaches, we show that porcine miR-194a-5p (homologous to human miR-194-5p) regulates TLR4 gene expression, an important molecule involved in pathogen virulence, recognition and activation of innate immunity in Salmonella infection.

The haustorial transcriptome of the cucurbit pathogen Podosphaera xanthii reveals new insights into the biotrophy and pathogenesis of powdery mildew fungi.

BACKGROUND: Podosphaera xanthii is the main causal agent of powdery mildew disease in cucurbits and is responsible for important yield losses in these crops worldwide. Powdery mildew fungi are obligate biotrophs. In these parasites, biotrophy is determined by the presence of haustoria, which are specialized structures of parasitism developed by these fungi for the acquisition of nutrients and the delivery of effectors. Detailed molecular studies of powdery mildew haustoria are scarce due mainly to difficulties in their isolation. Therefore, their analysis is considered an important challenge for powdery mildew research. The aim of this work was to gain insights into powdery mildew biology by analysing the haustorial transcriptome of P. xanthii. RESULTS: Prior to RNA isolation and massive-scale mRNA sequencing, a flow cytometric approach was developed to isolate P. xanthii haustoria free of visible contaminants. Next, several commercial kits were used to isolate total RNA and to construct the cDNA and Illumina libraries that were finally sequenced by the Illumina NextSeq system. Using this approach, the maximum amount of information from low-quality RNA that could be obtained was used to accomplish the de novo assembly of the P. xanthii haustorial transcriptome. The subsequent analysis of this transcriptome and comparison with the epiphytic transcriptome allowed us to identify the importance of several biological processes for haustorial cells such as protection against reactive oxygen species, the acquisition of different nutrients and genetic regulation mediated by non-coding RNAs. In addition, we could also identify several secreted proteins expressed exclusively in haustoria such as cell adhesion proteins that have not been related to powdery mildew biology to date. CONCLUSIONS: This work provides a novel approach to study the molecular aspects of powdery mildew haustoria. In addition, the results of this study have also allowed us to identify certain previously unknown processes and proteins involved in the biology of powdery mildews that could be essential for their biotrophy and pathogenesis.

Nitric oxide-dependent regulation of sweet pepper fruit ripening.

Ripening is a complex physiological process that involves changes in reactive nitrogen and oxygen species that govern the shelf-life and quality of fruits. Nitric oxide (NO)-dependent changes in the sweet pepper fruit transcriptome were determined by treating fruits at the initial breaking point stage with NO gas. Fruits were also harvested at the immature (green) and ripe (red) stages. Fruit ripening in the absence of NO resulted in changes in the abundance of 8805 transcripts whose function could be identified. Among these, functional clusters associated with reactive oxygen/nitrogen species and lipid metabolism were significantly modified. NO treatment resulted in the differential expression of 498 genes framed within these functional categories. Biochemical analysis revealed that NO treatment resulted in changes in fatty acid profiling, glutathione and proline contents, and the extent of lipid peroxidation, as well as increases in the activity of ascorbate peroxidase and lipoxygenase. These data provide supporting evidence for the crucial role of NO in the ripening of pepper fruit.

TransFlow: a modular framework for assembling and assessing accurate de novo transcriptomes in non-model organisms.

BACKGROUND: The advances in high-throughput sequencing technologies are allowing more and more de novo assembling of transcriptomes from many new organisms. Some degree of automation and evaluation is required to warrant reproducibility, repetitivity and the selection of the best possible transcriptome. Workflows and pipelines are becoming an absolute requirement for such a purpose, but the issue of assembling evaluation for de novo transcriptomes in organisms lacking a sequenced genome remains unsolved. An automated, reproducible and flexible framework called TransFlow to accomplish this task is described. RESULTS: TransFlow with its five independent modules was designed to build different workflows depending on the nature of the original reads. This architecture enables different combinations of Illumina and Roche/454 sequencing data, and can be extended to other sequencing platforms. Its capabilities are illustrated with the selection of reliable plant reference transcriptomes and the assembling six transcriptomes (three case studies for grapevine leaves, olive tree pollen, and chestnut stem, and other three for haustorium, epiphytic structures and their combination for the phytopathogenic fungus Podosphaera xanthii). Arabidopsis and poplar transcriptomes revealed to be the best references. A common result regarding de novo assemblies is that Illumina paired-end reads of 100 nt in length assembled with OASES can provide reliable transcriptomes, while the contribution of longer reads is noticeable only when they complement a set of short, single-reads. CONCLUSIONS: TransFlow can handle up to 181 different assembling strategies. Evaluation based on principal component analyses allows its self-adaptation to different sets of reads to provide a suitable transcriptome for each combination of reads and assemblers. As a result, each case study has its own behaviour, prioritises evaluation parameters, and gives an objective and automated way for detecting the best transcriptome within a pool of them. Sequencing data type and quantity (preferably several hundred millions of 2×100 nt or longer), assemblers (OASES for Illumina, MIRA4 and EULER-SR reconciled with CAP3 for Roche/454) and strategy (preferably scaffolding with OASES, and probably merging with Roche/454 when available) arise as the most impacting factors.

Regulatory role of microRNA in mesenteric lymph nodes after Salmonella Typhimurium infection.

Salmonellosis is a gastrointestinal disease caused by non-typhoidal Salmonella serovars such as Salmonella Typhimurium. This pathology is a zoonosis, and food animals with subclinical infection constitute a vast reservoir for disease. After intestinal colonization, Salmonella Typhimurium reaches mesenteric lymph nodes (MLN), where infection is controlled avoiding systemic spread. Although the molecular basis of this infection has been extensively studied, little is known about how microRNA (miRNA) regulate the expression of proteins involved in the Salmonella-host interaction. Using small RNA-seq, we examined expression profiles of MLN 2 days after infection with Salmonella Typhimurium, and we found 110 dysregulated miRNA. Among them, we found upregulated miR-21, miR-155, miR-150, and miR-221, as well as downregulated miR-143 and miR-125, all of them previously linked to other bacterial infections. Integration with proteomic data revealed 30 miRNA potentially regulating the expression of 15 proteins involved in biological functions such as cell death and survival, inflammatory response and antigenic presentation. The inflammatory response was found increased via upregulation of miRNA such as miR-21 and miR-155. Downregulation of miR-125a/b, miR-148 and miR-1 were identified as potential regulators of MHC-class I components PSMB8, HSP90B1 and PDIA3, respectively. Furthermore, we confirmed that miR-125a is a direct target of immunoproteasome component PSMB8. Since we also found miR-130 downregulation, which is associated with upregulation of HSPA8, we suggest induction of both MHC-I and MHC-II antigen presentation pathways. In conclusion, our study identifies miRNA that could regulate critical networks for antigenic presentation, inflammatory response and cytoskeletal rearrangements.

Generation of nitric oxide by olive (Olea europaea L.) pollen during in vitro germination and assessment of the S-nitroso- and nitro-proteomes by computational predictive methods.

Nitric oxide is recognized as a signaling molecule involved in a broad range of physiological processes in plants including sexual reproduction. NO has been detected in the pollen grain at high levels and regulates pollen tube growth. Previous studies demonstrated that NO as well as ROS are produced in the olive reproductive tissues in a stage- and tissue-specific manner. The aim of this study was to assess the production of NO throughout the germination of olive (Olea europaea L.) pollen in vitro. The NO fluorescent probe DAF-2DA was used to image NO production in situ, which was correlated to pollen viability. Moreover, by means of a fluorimetric assay we showed that growing pollen tubes release NO. GSNO -a mobile reservoir of NO, formed by the S-nitrosylation of NO with reduced glutathione (GSH) - was for the first time detected and quantified at different stages of pollen tube growth using a LC-ES/MS analysis. Exogenous NO donors inhibited both pollen germination and pollen tube growth and these effects were partially reverted by the specific NO-scavenger c-PTIO. However, little is known about how NO affects the germination process. With the aim of elucidating the putative relevance of protein S-nitrosylation and Tyr-nitration as important post-translational modifications in the development and physiology of the olive pollen, a de novo assembled and annotated reproductive transcriptome from olive was challenged in silico for the putative capability of transcripts to become potentially modified by S-nitrosylation/Tyr-nitration according to well-established criteria. Numerous gene products with these characteristics were identified, and a broad discussion as regards to their potential role in plant reproduction was built after their functional classification. Moreover, the importance of both S-nitrosylation/Tyr-nitrations was experimentally assessed and validated by using Western blotting, immunoprecipitation and proteomic approaches.

Automated identification of reference genes based on RNA-seq data.

BACKGROUND: Gene expression analyses demand appropriate reference genes (RGs) for normalization, in order to obtain reliable assessments. Ideally, RG expression levels should remain constant in all cells, tissues or experimental conditions under study. Housekeeping genes traditionally fulfilled this requirement, but they have been reported to be less invariant than expected; therefore, RGs should be tested and validated for every particular situation. Microarray data have been used to propose new RGs, but only a limited set of model species and conditions are available; on the contrary, RNA-seq experiments are more and more frequent and constitute a new source of candidate RGs. RESULTS: An automated workflow based on mapped NGS reads has been constructed to obtain highly and invariantly expressed RGs based on a normalized expression in reads per mapped million and the coefficient of variation. This workflow has been tested with Roche/454 reads from reproductive tissues of olive tree (Olea europaea L.), as well as with Illumina paired-end reads from two different accessions of Arabidopsis thaliana and three different human cancers (prostate, small-cell cancer lung and lung adenocarcinoma). Candidate RGs have been proposed for each species and many of them have been previously reported as RGs in literature. Experimental validation of significant RGs in olive tree is provided to support the algorithm. CONCLUSION: Regardless sequencing technology, number of replicates, and library sizes, when RNA-seq experiments are designed and performed, the same datasets can be analyzed with our workflow to extract suitable RGs for subsequent PCR validation. Moreover, different subset of experimental conditions can provide different suitable RGs.

Establishing gene models from the Pinus pinaster genome using gene capture and BAC sequencing.

BACKGROUND: In the era of DNA throughput sequencing, assembling and understanding gymnosperm mega-genomes remains a challenge. Although drafts of three conifer genomes have recently been published, this number is too low to understand the full complexity of conifer genomes. Using techniques focused on specific genes, gene models can be established that can aid in the assembly of gene-rich regions, and this information can be used to compare genomes and understand functional evolution. RESULTS: In this study, gene capture technology combined with BAC isolation and sequencing was used as an experimental approach to establish de novo gene structures without a reference genome. Probes were designed for 866 maritime pine transcripts to sequence genes captured from genomic DNA. The gene models were constructed using GeneAssembler, a new bioinformatic pipeline, which reconstructed over 82% of the gene structures, and a high proportion (85%) of the captured gene models contained sequences from the promoter regulatory region. In a parallel experiment, the P. pinaster BAC library was screened to isolate clones containing genes whose cDNA sequence were already available. BAC clones containing the asparagine synthetase, sucrose synthase and xyloglucan endotransglycosylase gene sequences were isolated and used in this study. The gene models derived from the gene capture approach were compared with the genomic sequences derived from the BAC clones. This combined approach is a particularly efficient way to capture the genomic structures of gene families with a small number of members. CONCLUSIONS: The experimental approach used in this study is a valuable combined technique to study genomic gene structures in species for which a reference genome is unavailable. It can be used to establish exon/intron boundaries in unknown gene structures, to reconstruct incomplete genes and to obtain promoter sequences that can be used for transcriptional studies. A bioinformatics algorithm (GeneAssembler) is also provided as a Ruby gem for this class of analyses.

Transcriptional analysis of porcine intestinal mucosa infected with Salmonella Typhimurium revealed a massive inflammatory response and disruption of bile acid absorption in ileum.

Infected pork meat is an important source of non-typhoidal human salmonellosis. Understanding of molecular mechanisms involved in disease pathogenesis is important for the development of therapeutic and preventive strategies. Thus, hereby we study the transcriptional profiles along the porcine intestine during infection with Salmonella Typhimurium, as well as post-transcriptional gene modulation by microRNAs (miRNA). Sixteen piglets were orally challenged with S. Typhimurium. Samples from jejunum, ileum and colon, collected 1, 2 and 6 days post infection (dpi) were hybridized to mRNA and miRNA expression microarrays and analyzed. Jejunum showed a reduced transcriptional response indicating mild inflammation only at 2 dpi. In ileum inflammatory genes were overexpressed (e.g., IL-1B, IL-6, IL-8, IL1RAP, TNFα), indicating a strong immune response at all times of infection. Infection also down-regulated genes of the FXR pathway (e.g., NR1H4, FABP6, APOA1, SLC10A2), indicating disruption of the bile acid absorption in ileum. This result was confirmed by decreased high-density lipoprotein cholesterol in serum of infected pigs. Ileal inflammatory gene expression changes peaked at 2 dpi and tended to resolve at 6 dpi. Furthermore, miRNA analysis of ileum at 2 dpi revealed 62 miRNAs potentially regulating target genes involved in this inflammatory process (e.g., miR-374 and miR-451). In colon, genes involved in epithelial adherence, proliferation and cellular reorganization were down-regulated at 2 and 6 dpi. In summary, here we show the transcriptional changes occurring at the intestine at different time points of the infection, which are mainly related to inflammation and disruption of the bile acid metabolism.

De novo assembly, characterization and functional annotation of Senegalese sole (Solea senegalensis) and common sole (Solea solea) transcriptomes: integration in a database and design of a microarray.

BACKGROUND: Senegalese sole (Solea senegalensis) and common sole (S. solea) are two economically and evolutionary important flatfish species both in fisheries and aquaculture. Although some genomic resources and tools were recently described in these species, further sequencing efforts are required to establish a complete transcriptome, and to identify new molecular markers. Moreover, the comparative analysis of transcriptomes will be useful to understand flatfish evolution. RESULTS: A comprehensive characterization of the transcriptome for each species was carried out using a large set of Illumina data (more than 1,800 millions reads for each sole species) and 454 reads (more than 5 millions reads only in S. senegalensis), providing coverages ranging from 1,384x to 2,543x. After a de novo assembly, 45,063 and 38,402 different transcripts were obtained, comprising 18,738 and 22,683 full-length cDNAs in S. senegalensis and S. solea, respectively. A reference transcriptome with the longest unique transcripts and putative non-redundant new transcripts was established for each species. A subset of 11,953 reference transcripts was qualified as highly reliable orthologs (>97% identity) between both species. A small subset of putative species-specific, lineage-specific and flatfish-specific transcripts were also identified. Furthermore, transcriptome data permitted the identification of single nucleotide polymorphisms and simple-sequence repeats confirmed by FISH to be used in further genetic and expression studies. Moreover, evidences on the retention of crystallins crybb1, crybb1-like and crybb3 in the two species of soles are also presented. Transcriptome information was applied to the design of a microarray tool in S. senegalensis that was successfully tested and validated by qPCR. Finally, transcriptomic data were hosted and structured at SoleaDB. CONCLUSIONS: Transcriptomes and molecular markers identified in this study represent a valuable source for future genomic studies in these economically important species. Orthology analysis provided new clues regarding sole genome evolution indicating a divergent evolution of crystallins in flatfish. The design of a microarray and establishment of a reference transcriptome will be useful for large-scale gene expression studies. Moreover, the integration of transcriptomic data in the SoleaDB will facilitate the management of genomic information in these important species.

Pyroptosis and adaptive immunity mechanisms are promptly engendered in mesenteric lymph-nodes during pig infections with Salmonella enterica serovar Typhimurium.

In this study, we explored the transcriptional response and the morphological changes occurring in porcine mesenteric lymph-nodes (MLN) along a time course of 1, 2 and 6 days post infection (dpi) with Salmonella Typhimurium. Additionally, we analysed the expression of some Salmonella effectors in tissue to complete our view of the processes triggered in these organs upon infection. The results indicate that besides dampening apoptosis, swine take advantage of the flagellin and prgJ expression by Salmonella Typhimuriun to induce pyroptosis in MLN, preventing bacterial dissemination. Furthermore, cross-presentation of Salmonella antigens was inferred as a mechanism that results in a rapid clearance of pathogen by cytotoxic T cells. In summary, although the Salmonella Typhimurium strain employed in this study was able to express some of its major virulence effectors in porcine MLN, a combination of early innate and adaptive immunity mechanisms might overcome virulence strategies employed by the pathogen, enabling the host to protect itself against bacterial spread beyond gut-associated lymph-nodes. Interestingly, we deduced that clathrin-mediated endocytosis could contribute to mechanisms of pathogen virulence and/or host defence in MLN of Salmonella infected swine. Taken together, our results are useful for a better understanding of the critical protective mechanisms against Salmonella that occur in porcine MLN to prevent the spread of infection beyond the intestine.

Salmonella Typhimurium induces genome-wide expression and phosphorylation changes that modulate immune response, intracellular survival and vesicle transport in infected neutrophils.

Salmonella Typhimurium is a food-borne pathogen that causes salmonellosis. When in contact with the host, neutrophils are rapidly recruited to act as first line of defense. To better understand the pathogenesis of this infection, we used an in vitro model of neutrophil infection to perform dual RNA-sequencing (both host and pathogen). In addition, and given that many pathogens interfere with kinase-mediated phosphorylation in host signaling, we performed a phosphoproteomic analysis. The immune response was overall diminished in infected neutrophils, mainly JAK/STAT and toll-like receptor signaling pathways. We found decreased expression of proinflammatory cytokine receptor genes and predicted downregulation of the mitogen-activated protein (MAPK) signaling pathway. Also, Salmonella infection inhibited interferons I and II signaling pathways by upregulation of SOCS3 and subsequent downregulation of STAT1 and STAT2. Additionally, phosphorylation of PSMC2 and PSMC4, proteasome regulatory proteins, was decreased in infected neutrophils. Cell viability and survival was increased by p53 signaling, cell cycle arrest and NFkB-proteasome pathways activation. Combined analysis of RNA-seq and phosphoproteomics also revealed inhibited vesicle transport mechanisms mediated by dynein/dynactin and exocyst complexes, involved in ER-to-Golgi transport and centripetal movement of lysosomes and endosomes. Among the overexpressed virulence genes from Salmonella we found potential effectors responsible of these dysregulations, such as spiC, sopD2, sifA or pipB2, all of them involved in intracellular replication. Our results suggest that Salmonella induces (through overexpression of virulence factors) transcriptional and phosphorylation changes that increases neutrophil survival and shuts down immune response to minimize host response, and impairing intracellular vesicle transport likely to keep nutrients for replication and Salmonella-containing vacuole formation and maintenance.

Construction of miRNA-mRNA networks for the identification of lung cancer biomarkers in liquid biopsies.

Lung cancer (LC) is the most common cause of cancer death worldwide mostly due to the low survival rate: 75% of cases are identified in advanced stages. In this study, the list of useful biomarkers to make an early diagnosis using liquid biopsies was expanded. A total of 30 samples of LC were analyzed to define potential miRNA biomarkers in liquid biopsies for LC. The biomarkers have been identified in interaction networks miRNA-mRNA. The potential biomarkers have been then validated in large cohorts. A total of 15 candidate miRNAs, that regulate the repression of 30 mRNAs, have been identified as a specific functional interaction network for squamous carcinoma, while the specific functional interaction network of adenocarcinoma consists of four candidate miRNAs that seem to handle the repression of five mRNA. Inspection of expression levels in larger cohorts validates the usefulness of the 11 candidates as biomarkers in liquid biopsies. The 11 candidate miRNAs found could be utilized to form diagnostic predictive biomarkers for LC in liquid biopsies.

An innovative framework to determine the implementation level of personalized medicine: A systematic review.

Background: Personalized medicine (PM) is now the new frontier in patient care. The application of this new paradigm extends to various pathologies and different patient care phases, such as diagnosis and treatment. Translating biotechnological advances to clinical routine means adapting health services at all levels is necessary. Purpose: This article aims to identify the elements for devising a framework that will allow the level of PM implementation in the country under study to be quantitatively and qualitatively assessed and that can be used as a guideline for future implementation plans. Methods: A systematic review was conducted per the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement. The research question was: What are the domains for determining the level of implementation of PM at the national level? The domains for assessing the degree of PM implementation, which would form the framework, were established. Results: 19 full-text studies that met the inclusion criteria were peer-selected in the systematic review. From all the studies that were included, 37 elements-encompassed in 11 domains-were extracted for determining the degree of PM implementation. These domains and their constituent elements comprise the qualitative and quantitative assessment framework presented herein. Each of the elements can be assessed individually. On the other hand, the domains were standardized to all have the same weight in an overall assessment. Conclusions: A framework has been developed that takes a multi-factorial approach to determine the degree of implementation of PM at the national level. This framework could also be used to rank countries and their implementation strategies according to the score they receive in the application of the latter. It could also be used as a guide for developing future national PM implementation strategies. Systematic review registration: https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42022338611, Identifier: CRD42022338611.

MOWServ: a web client for integration of bioinformatic resources.

The productivity of any scientist is affected by cumbersome, tedious and time-consuming tasks that try to make the heterogeneous web services compatible so that they can be useful in their research. MOWServ, the bioinformatic platform offered by the Spanish National Institute of Bioinformatics, was released to provide integrated access to databases and analytical tools. Since its release, the number of available services has grown dramatically, and it has become one of the main contributors of registered services in the EMBRACE Biocatalogue. The ontology that enables most of the web-service compatibility has been curated, improved and extended. The service discovery has been greatly enhanced by Magallanes software and biodataSF. User data are securely stored on the main server by an authentication protocol that enables the monitoring of current or already-finished user's tasks, as well as the pipelining of successive data processing services. The BioMoby standard has been greatly extended with the new features included in the MOWServ, such as management of additional information (metadata such as extended descriptions, keywords and datafile examples), a qualified registry, error handling, asynchronous services and service replication. All of them have increased the MOWServ service quality, usability and robustness. MOWServ is available at http://www.inab.org/MOWServ/ and has a mirror at http://www.bitlab-es.com/MOWServ/.

EuroPineDB: a high-coverage web database for maritime pine transcriptome.

BACKGROUND: Pinus pinaster is an economically and ecologically important species that is becoming a woody gymnosperm model. Its enormous genome size makes whole-genome sequencing approaches are hard to apply. Therefore, the expressed portion of the genome has to be characterised and the results and annotations have to be stored in dedicated databases. DESCRIPTION: EuroPineDB is the largest sequence collection available for a single pine species, Pinus pinaster (maritime pine), since it comprises 951 641 raw sequence reads obtained from non-normalised cDNA libraries and high-throughput sequencing from adult (xylem, phloem, roots, stem, needles, cones, strobili) and embryonic (germinated embryos, buds, callus) maritime pine tissues. Using open-source tools, sequences were optimally pre-processed, assembled, and extensively annotated (GO, EC and KEGG terms, descriptions, SNPs, SSRs, ORFs and InterPro codes). As a result, a 10.5× P. pinaster genome was covered and assembled in 55 322 UniGenes. A total of 32 919 (59.5%) of P. pinaster UniGenes were annotated with at least one description, revealing at least 18 466 different genes. The complete database, which is designed to be scalable, maintainable, and expandable, is freely available at: http://www.scbi.uma.es/pindb/. It can be retrieved by gene libraries, pine species, annotations, UniGenes and microarrays (i.e., the sequences are distributed in two-colour microarrays; this is the only conifer database that provides this information) and will be periodically updated. Small assemblies can be viewed using a dedicated visualisation tool that connects them with SNPs. Any sequence or annotation set shown on-screen can be downloaded. Retrieval mechanisms for sequences and gene annotations are provided. CONCLUSIONS: The EuroPineDB with its integrated information can be used to reveal new knowledge, offers an easy-to-use collection of information to directly support experimental work (including microarray hybridisation), and provides deeper knowledge on the maritime pine transcriptome.

Chromosome anchoring in Senegalese sole (Solea senegalensis) reveals sex-associated markers and genome rearrangements in flatfish.

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.