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.

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.

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.

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.

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.

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.

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.

Decoding Neuromuscular Disorders Using Phenotypic Clusters Obtained From Co-Occurrence Networks

Neuromuscular disorders (NMDs) represent an important subset of rare diseases associated with elevated morbidity and mortality whose diagnosis can take years. Here we present a novel approach using systems biology to produce functionally-coherent phenotype clusters that provide insight into the cellular functions and phenotypic patterns underlying NMDs, using the Human Phenotype Ontology as a common framework. Gene and phenotype information was obtained for 424 NMDs in OMIM and 126 NMDs in Orphanet, and 335 and 216 phenotypes were identified as typical for NMDs, respectively. ‘Elevated serum creatine kinase’ was the most specific to NMDs, in agreement with the clinical test of elevated serum creatinine kinase that is conducted on NMD patients. The approach to obtain co-occurring NMD phenotypes was validated based on co-mention in PubMed abstracts. A total of 231 (OMIM) and 150 (Orphanet) clusters of highly connected co-occurrent NMD phenotypes were obtained. In parallel, a tripartite network based on phenotypes, diseases and genes was used to associate NMD phenotypes with functions, an approach also validated by literature co-mention, with KEGG pathways showing proportionally higher overlap than Gene Ontology and Reactome. Phenotype-function pairs were crossed with the co-occurrent NMD phenotype clusters to obtain 40 (OMIM) and 72 (Orphanet) functionally coherent phenotype clusters. As expected, many of these overlapped with known diseases and confirmed existing knowledge. Other clusters revealed interesting new findings, indicating informative phenotypes for differential diagnosis, providing deeper knowledge of NMDs, and pointing towards specific cell dysfunction caused by pleiotropic genes. This work is an example of reproducible research that i) can help better understand NMDs and support their diagnosis by providing a new tool that exploits existing information to obtain novel clusters of functionally-related phenotypes, and ii) takes us another step towards personalised medicine for NMDs.

Gene Expression Profile of Mexican Lime (Citrus aurantifolia) Trees in Response to Huanglongbing Disease caused by Candidatus Liberibacter asiaticus

Nowadays, Huanglongbing (HLB) disease, associated with Candidatus Liberibacter asiaticus (CLas), seriously affects citriculture worldwide, and no cure is currently available. Transcriptomic analysis of host-pathogen interaction is the first step to understand the molecular landscape of a disease. Previous works have reported the transcriptome profiling in response to HLB in different susceptible citrus species; however, similar studies in tolerant citrus species, including Mexican lime, are limited. In this work, we have obtained an RNA-seq-based differential expression profile of Mexican lime plants challenged against CLas infection, at both asymptomatic and symptomatic stages. Typical HLB-responsive differentially expressed genes (DEGs) are involved in photosynthesis, secondary metabolism, and phytohormone homeostasis. Enrichment of DEGs associated with biotic response showed that genes related to cell wall, secondary metabolism, transcription factors, signaling, and redox reactions could play a role in the tolerance of Mexican lime against CLas infection. Interestingly, despite some concordance observed between transcriptional responses of different tolerant citrus species, a subset of DEGs appeared to be species-specific. Our data highlights the importance of studying the host response during HLB disease using as model tolerant citrus species, in order to design new and opportune diagnostic and management methods.

Expression-based, consistent biomarkers for prognosis and diagnosis in lung cancer

OVERVIEW

Lung cancer is one of the deadliest cancers in the world. Its histological classification depends on early diagnosis and successful treatment. Therefore, having specific biomarkers for a quick sorting widens the successful output of lung cancer treatment.

MATERIAL AND METHODS

High-throughput sequencing (RNA-seq) was performed of small cohorts of BioBanco samples from healthy and tumour cells from lung adenocarcinoma (LUAD) and squamous cell carcinoma of the lung (lSCC). RNA-seq samples from small cell lung cancer (SCLC) were downloaded from databases. A bioinformatic workflow has been programmed for the identification of differentially expressed genes (DEGs).

RESULTS

A total of 4777 DEGs were differentially expressed in SCLC, 3676 DEGs were in lSCC, while the lowest number of DEGs, 2819, appeared in LUAD. Among them, 945 DEGs were common to the three histological types. Once validated their expression profile and their survival predictive capacity in large, public cohorts, three DEGs can be exclusively considered as diagnostic biomarkers, three as prognosis biomarkers, and other three exhibit both diagnosis and prognosis capabilities.

CONCLUSIONS

This prospective study presents evidences for the diagnostic and prognostic capabilities of expression changes in CAPN8-2, TMC5 and MUC1 in LUAD, while they are non-significant in SCLC and lSCC. Their translation to clinical practice is proposed.

De novo assembly and functional annotation of Citrus aurantifolia transcriptome from Candidatus Liberibacter asiaticus infected and non-infected trees

Mexican lime (Citrus aurantifolia) belongs to the Rutaceae family and nowadays is one of the major commercial citrus crops in different countries. In Mexico, Mexican lime production is impaired by Huanglongbing (HLB) disease associated to Candidatus Liberibacter asiaticus (CLas) bacteria. To date, transcriptomic studies of CLas-Citrus interaction, have been performed mainly in sweet citrus models at symptomatic (early) stage where pleiotropic responses could mask important, pathogen-driven host modulation as well as, host antibacterial responses. Additionally, well-assembled reference transcriptomes for acid limes including C. aurantifolia are not available. The development of improved transcriptomic resources for CLas-citrus pathosystem, including both asymptomatic (early) and symptomatic (late) stages, could accelerate the understanding of the disease. Here, we provide the first transcriptomic analysis from healthy and HLB-infected C. aurantifolia leaves at both asymptomatic and symptomatic stages, using a RNA-seq approach in the Illumina NexSeq500 platform. The construction of the assembled transcriptome was conducted using the predesigned workflow Transflow and a total of 41,522 tentative transcripts (TTs) obtained. These C. aurantifolia TTs were functionally annotated using TAIR10 and UniProtKB databases. All raw reads were deposited in the NCBI SRA with accession numbers SRR10353556, SRR10353558, SRR10353560 and SRR10353562. Overall, this dataset adds new transcriptomic valuable tools for future breeding programs, will allow the design of novel diagnostic molecular markers, and will be an essential tool for studying the HLB disease.

Biomarker potential of repetitive-element transcriptome in lung cancer

Since repetitive elements (REs) account for nearly 53% of the human genome, profiling its transcription after an oncogenic change might help in the search for new biomarkers. Lung cancer was selected as target since it is the most frequent cause of cancer death. A bioinformatic workflow based on well-established bioinformatic tools (such as RepEnrich, RepBase, SAMTools, edgeR and DESeq2) has been developed to identify differentially expressed RNAs from REs. It was trained and tested with public RNA-seq data from matched sequencing of tumour and healthy lung tissues from the same patient to reveal differential expression within the RE transcriptome. Healthy lung tissues express a specific set of REs whose expression, after an oncogenic process, is strictly and specifically changed. Discrete sets of differentially expressed REs were found for lung adenocarcinoma, for small-cell lung cancer, and for both cancers. Differential expression affects more HERV-than LINE-derived REs and seems biased towards down-regulation in cancer cells. REs behaving consistently in all patients were tested in a different patient cohort to validate the proposed biomarkers. Down-regulation of AluYg6 and LTR18B was confirmed as potential lung cancer biomarkers, while up-regulation of HERVK11D-Int is specific for lung adenocarcinoma and up-regulation of UCON88 is specific for small cell lung cancer. Hence, the study of RE transcriptome might be considered another research target in cancer, making REs a promising source of lung cancer biomarkers.

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.

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5938-0) contains supplementary material, which is available to authorized users.

TarSynFlow, a workflow for bacterial genome comparisons that revealed genes putatively involved in the probiotic character of Shewanella putrefaciens strain Pdp11

Probiotic microorganisms are of great interest in clinical, livestock and aquaculture. Knowledge of the genomic basis of probiotic characteristics can be a useful tool to understand why some strains can be pathogenic while others are probiotic in the same species. An automatized workflow called TarSynFlow (Targeted Synteny Workflow) has been then developed to compare finished or draft bacterial genomes based on a set of proteins. When used to analyze the finished genome of the probiotic strain Pdp11 of Shewanella putrefaciens and genome drafts from seven known non-probiotic strains of the same species obtained in this work, 15 genes were found exclusive of Pdp11. Their presence was confirmed by PCR using Pdp11-specific primers. Functional inspection of the 15 genes allowed us to hypothesize that Pdp11 underwent genome rearrangements spurred by plasmids and mobile elements. As a result, Pdp11 presents specific proteins for gut colonization, bile salt resistance and gut pathogen adhesion inhibition, which can explain some probiotic features of Pdp11.

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.

First draft genome assembly of the Argane tree ( Argania spinosa)

Background: The Argane tree ( Argania spinosa L. Skeels) is an endemic tree of mid-western Morocco that plays an important socioeconomic and ecologic role for a dense human population in an arid zone. Several studies confirmed the importance of this species as a food and feed source and as a resource for both pharmaceutical and cosmetic compounds. Unfortunately, the argane tree ecosystem is facing significant threats from environmental changes (global warming, over-population) and over-exploitation. Limited research has been conducted, however, on argane tree genetics and genomics, which hinders its conservation and genetic improvement. Methods: Here, we present a draft genome assembly of A. spinosa. A reliable reference genome of  A. spinosa was created using a hybrid  de novo assembly approach combining short and long sequencing reads. Results: In total, 144 Gb Illumina HiSeq reads and 7.6 Gb PacBio reads were produced and assembled. The final draft genome comprises 75 327 scaffolds totaling 671 Mb with an N50 of 49 916 kb. The draft assembly is close to the genome size estimated by k-mers distribution and covers 89% of complete and 4.3 % of partial Arabidopsis orthologous groups in BUSCO. Conclusion: The A. spinosa genome will be useful for assessing biodiversity leading to efficient conservation of this endangered endemic tree. Furthermore, the genome may enable genome-assisted cultivar breeding, and provide a better understanding of important metabolic pathways and their underlying genes for both cosmetic and pharmacological.

First draft genome assembly of the Argane tree ( Argania spinosa)

Background: The Argane tree ( Argania spinosa L. Skeels) is an endemic tree of mid-western Morocco that plays an important socioeconomic and ecologic role for a dense human population in an arid zone. Several studies confirmed the importance of this species as a food and feed source and as a resource for both pharmaceutical and cosmetic compounds. Unfortunately, the argane tree ecosystem is facing significant threats from environmental changes (global warming, over-population) and over-exploitation. Limited research has been conducted, however, on argane tree genetics and genomics, which hinders its conservation and genetic improvement. Methods: Here, we present a draft genome assembly of A. spinosa. A reliable reference genome of  A. spinosa was created using a hybrid  de novo assembly approach combining short and long sequencing reads. Results: In total, 144 Gb Illumina HiSeq reads and 7.6 Gb PacBio reads were produced and assembled. The final draft genome comprises 75 327 scaffolds totaling 671 Mb with an N50 of 49 916 kb. The draft assembly is close to the genome size estimated by k-mers distribution and covers 89% of complete and 4.3 % of partial Arabidopsis orthologous groups in BUSCO. Conclusion: The A. spinosa genome will be useful for assessing biodiversity leading to efficient conservation of this endangered endemic tree. Furthermore, the genome may enable genome-assisted cultivar breeding, and provide a better understanding of important metabolic pathways and their underlying genes for both cosmetic and pharmacological.

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.

S-nitroso- and nitro- proteomes in the olive (Olea europaea L.) pollen. Predictive versus experimental data by nano-LC-MS

The data presented here are related to the research article entitled "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" doi:10.1016/j.niox.2017.06.005 (Jimenez-Quesada et al., 2017) [1]. Predicted cysteine S-nitrosylation and Tyr-nitration sites in proteins derived from a de novo assembled and annotated pollen transcriptome from olive tree (Olea europaea L.) were obtained after using well-established predictive tools in silico. Predictions were performed using both default and highly restrictive thresholds. Numerous gene products identified with these characteristics are listed here. An experimental validation of the data, consisting in nano-LC-MS (Liquid Chromatography-Mass Spectrometry) determination of olive pollen proteins after immunoprecipitation with antibodies to anti-S-nitrosoCys and anti-3-NT (NitroTyrosine) allowed identification of numerous proteins subjected to these two post-translational modifications, which are listed here together with information regarding their cross-presence among the predictions.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12938-017-0356-5) contains supplementary material, which is available to authorized users.

Genomic analysis of the marine fish pathogen Photobacterium damselae subsp. piscicida: Insertion sequences proliferation is associated with chromosomal reorganisations and rampant gene decay

Photobacterium damselae subsp. piscicida (Pdp) is an intracellular fish pathogen that causes photobacteriosis, a disease proven deadly in farmed fish worldwide. This work focuses on the analysis of genome sequences, chromosomes structure and gene contents of two strains from Sparus aurata (DI21) and Solea senegalensis (L091106-03H), isolated on the Spanish Atlantic coast. The comparative genomic analysis revealed that DI21 and L091106-03H share 98% of their genomes, including two virulence plasmids: pPHDP70 encoding siderophore piscibactin synthesis and pPHDP10 encoding the apoptotic toxin AIP56. Both genomes harbour a surprisingly large number of IS elements accounting for 12-17% of the total genome, representing an IS density of 0.15 elements per kb, one of the highest IS density values in a bacterial pathogen. This massive proliferation of ISs is responsible for the generation of a high number of pseudogenes that caused extensive loss of biological functions. Pseudogene formation is one of the main features of Pdp genome that explains most of the ecological and phenotypic differences with respect to its sibling subspecies P. damselae subsp. damselae and to other Vibrionaceae. Evidence was also found proving the existence of two chromosomal configurations depending on the origin of the strains: an European and an Asian/American types of genome organisation, reinforcing the idea of the existence of two geographically-linked clonal lineages in Pdp. In short, our study suggests that the host-dependent lifestyle of Pdp allowed massive IS proliferation and gene decay processes, which are major evolutionary forces in the shaping of the Pdp genome.

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.

Identification and Functional Annotation of Genes Differentially Expressed in the Reproductive Tissues of the Olive Tree (Olea europaea L.) through the Generation of Subtractive Libraries

The olive tree is a crop of high socio-economical importance in the Mediterranean area. Sexual reproduction in this plant is an essential process, which determines the yield. Successful fertilization is mainly favored and sometimes needed of the presence of pollen grains from a different cultivar as the olive seizes a self-incompatibility system allegedly determined of the sporophytic type. The purpose of the present study was to identify key gene products involved in the function of olive pollen and pistil, in order to help elucidate the events and signaling processes, which happen during the courtship, pollen grain germination, and fertilization in olive. The use of subtractive SSH libraries constructed using, on the one hand one specific stage of the pistil development with germinating pollen grains, and on the other hand mature pollen grains may help to reveal the specific transcripts involved in the cited events. Such libraries have also been created by subtracting vegetative mRNAs (from leaves), in order to identify reproductive sequences only. A variety of transcripts have been identified in the mature pollen grains and in the pistil at the receptive stage. Among them, those related to defense, transport and oxidative metabolism are highlighted mainly in the pistil libraries where transcripts related to stress, and response to biotic and abiotic stimulus have a prominent position. Extensive lists containing information as regard to the specific transcripts determined for each stage and tissue are provided, as well as functional classifications of these gene products. Such lists were faced up to two recent datasets obtained in olive after transcriptomic and genomic approaches. The sequences and the differential expression level of the SSH-transcripts identified here, highly matched the transcriptomic information. Moreover, the unique presence of a representative number of these transcripts has been validated by means of qPCR approaches. The construction of SSH libraries using pistil and pollen, considering the high interaction between male-female counterparts, allowed the identification of transcripts with important roles in stigma physiology. The functions of many of the transcripts obtained are intimately related, and most of them are of pivotal importance in defense, pollen-stigma interaction and signaling.

Development of genomic tools in a widespread tropical tree, Symphonia globulifera L.f.: a new low-coverage draft genome, SNP and SSR markers

Population genetic studies in tropical plants are often challenging because of limited information on taxonomy, phylogenetic relationships and distribution ranges, scarce genomic information and logistic challenges in sampling. We describe a strategy to develop robust and widely applicable genetic markers based on a modest development of genomic resources in the ancient tropical tree species Symphonia globulifera L.f. (Clusiaceae), a keystone species in African and Neotropical rainforests. We provide the first low-coverage (11X) fragmented draft genome sequenced on an individual from Cameroon, covering 1.027 Gbp or 67.5% of the estimated genome size. Annotation of 565 scaffolds (7.57 Mbp) resulted in the prediction of 1046 putative genes (231 of them containing a complete open reading frame) and 1523 exact simple sequence repeats (SSRs, microsatellites). Aligning a published transcriptome of a French Guiana population against this draft genome produced 923 high-quality single nucleotide polymorphisms. We also preselected genic SSRs in silico that were conserved and polymorphic across a wide geographical range, thus reducing marker development tests on rare DNA samples. Of 23 SSRs tested, 19 amplified and 18 were successfully genotyped in four S. globulifera populations from South America (Brazil and French Guiana) and Africa (Cameroon and São Tomé island, F_ST  = 0.34). Most loci showed only population-specific deviations from Hardy-Weinberg proportions, pointing to local population effects (e.g. null alleles). The described genomic resources are valuable for evolutionary studies in Symphonia and for comparative studies in plants. The methods are especially interesting for widespread tropical or endangered taxa with limited DNA availability.

Whole-exome sequencing analysis in twin sibling males with an anterior cruciate ligament rupture

Familial predisposition is among the major genetic risk factors for non-contact musculoskeletal tissue injuries. Personal genome sequence shows that different polymorphism profiles may account for the number and the degree of injuries and the recovery time. Genotyping studies allow investigation into genome factors with potential impact on pathogenesis of non-contact ligament injuries. We have studied a family with twin sibling males surgically diagnosed of an anterior cruciate ligament non-contact rupture and non-affected progenitors (father and mother) were subjected to whole exome sequencing (WES) analysis. WES analysis previously carried out on 16 individuals, without ACL injury medical records, were also included in this study for single nucleotide variants (SNVs) and small insertions and deletions detection (indels), variant filtering and to prioritize variants relative to the disease. WES analysis to identify SNVs and indels was performed using open web-based bioinformatics tools. A set of 11 new variants shared by family members can be associated to ACL non-contact injury, including SerpinA11, ARSI, NOCHT4, EPB41, FDFT1, POMC, KIF26A, OLFML2B, ATG7, FAH and WDR6. All of them, except ATG7 and WDR6, have shown a damaging predictive pattern by combinatorial standard predictive scores. In combination to the identified SNVs of EPB41 and SerpinA11 genes, ACTL7A gene showed a predicted deleterious variant reinforcing the idea these variants impact on of fibroblast-like cells deformability and ECM misbalance, Differential gene expression and RNA sequencing analysis will help to understand the combined participation of these protein coding genes in ACL non-contact injuries.

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.

Genomic characterization and expression analysis of four apolipoprotein A-IV paralogs in Senegalese sole (Solea senegalensis Kaup)

The apolipoprotein A-IV (ApoA-IV) plays a key role in lipid transport and feed intake regulation. In this work, four cDNA sequences encoding ApoA-IV paralogs were identified. Sequence analysis revealed conserved structural features including the common 33-codon block and nine repeated motifs. Gene structure analysis identified four exons and three introns except for apoA-IVAa1 (with only 3 exons). Synteny analysis showed that the four paralogs were structured into two clusters (cluster A containing apoA-IVAa1 and apoA-IVAa2 and cluster B with apoA-IVBa3 and apoA-IVBa4) linked to an apolipoprotein E. Phylogenetic analysis clearly separated the paralogs according to their cluster organization as well as revealed four subclades highly conserved in Acanthopterygii. Whole-mount analyses (WISH) in early larvae (0 and 1day post-hatch (dph)) showed that the four paralogs were mainly expressed in yolk syncytial layer surrounding the oil globules. Later, at 3 and 5dph, the four paralogs were mainly expressed in liver and intestine although with differences in their relative abundance and temporal expression patterns. Diet supply triggered the intensity of WISH signals in the intestine of the four paralogs. Quantification of mRNA abundance by qPCR using whole larvae only detected the induction by diet at 5dph. Moreover, transcript levels increased progressively with age except for apoA-IVAa2, which appeared as a low-expressed isoform. Expression analysis in juvenile tissues confirmed that the four paralogs were mainly expressed in liver and intestine and secondary in other tissues. The role of these ApoA-IV genes in lipid transport and the possible role of apoA-IVAa2 as a regulatory form are discussed.

ReprOlive: a database with linked data for the olive tree (Olea europaea L.) reproductive transcriptome

Plant reproductive transcriptomes have been analyzed in different species due to the agronomical and biotechnological importance of plant reproduction. Here we presented an olive tree reproductive transcriptome database with samples from pollen and pistil at different developmental stages, and leaf and root as control vegetative tissues http://reprolive.eez.csic.es). It was developed from 2,077,309 raw reads to 1,549 Sanger sequences. Using a pre-defined workflow based on open-source tools, sequences were pre-processed, assembled, mapped, and annotated with expression data, descriptions, GO terms, InterPro signatures, EC numbers, KEGG pathways, ORFs, and SSRs. Tentative transcripts (TTs) were also annotated with the corresponding orthologs in Arabidopsis thaliana from TAIR and RefSeq databases to enable Linked Data integration. It results in a reproductive transcriptome comprising 72,846 contigs with average length of 686 bp, of which 63,965 (87.8%) included at least one functional annotation, and 55,356 (75.9%) had an ortholog. A minimum of 23,568 different TTs was identified and 5,835 of them contain a complete ORF. The representative reproductive transcriptome can be reduced to 28,972 TTs for further gene expression studies. Partial transcriptomes from pollen, pistil, and vegetative tissues as control were also constructed. ReprOlive provides free access and download capability to these results. Retrieval mechanisms for sequences and transcript annotations are provided. Graphical localization of annotated enzymes into KEGG pathways is also possible. Finally, ReprOlive has included a semantic conceptualisation by means of a Resource Description Framework (RDF) allowing a Linked Data search for extracting the most updated information related to enzymes, interactions, allergens, structures, and reactive oxygen species.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-952) contains supplementary material, which is available to authorized users.

De novo assembly of maritime pine transcriptome: implications for forest breeding and biotechnology

Maritime pine (Pinus pinasterAit.) is a widely distributed conifer species in Southwestern Europe and one of the most advanced models for conifer research. In the current work, comprehensive characterization of the maritime pine transcriptome was performed using a combination of two different next-generation sequencing platforms, 454 and Illumina. De novo assembly of the transcriptome provided a catalogue of 26 020 unique transcripts in maritime pine trees and a collection of 9641 full-length cDNAs. Quality of the transcriptome assembly was validated by RT-PCR amplification of selected transcripts for structural and regulatory genes. Transcription factors and enzyme-encoding transcripts were annotated. Furthermore, the available sequencing data permitted the identification of polymorphisms and the establishment of robust single nucleotide polymorphism (SNP) and simple-sequence repeat (SSR) databases for genotyping applications and integration of translational genomics in maritime pine breeding programmes. All our data are freely available at SustainpineDB, the P. pinaster expressional database. Results reported here on the maritime pine transcriptome represent a valuable resource for future basic and applied studies on this ecological and economically important pine species.

A microRNA signature associated with early recurrence in breast cancer

Recurrent breast cancer occurring after the initial treatment is associated with poor outcome. A bimodal relapse pattern after surgery for primary tumor has been described with peaks of early and late recurrence occurring at about 2 and 5 years, respectively. Although several clinical and pathological features have been used to discriminate between low- and high-risk patients, the identification of molecular biomarkers with prognostic value remains an unmet need in the current management of breast cancer. Using microarray-based technology, we have performed a microRNA expression analysis in 71 primary breast tumors from patients that either remained disease-free at 5 years post-surgery (group A) or developed early (group B) or late (group C) recurrence. Unsupervised hierarchical clustering of microRNA expression data segregated tumors in two groups, mainly corresponding to patients with early recurrence and those with no recurrence. Microarray data analysis and RT-qPCR validation led to the identification of a set of 5 microRNAs (the 5-miRNA signature) differentially expressed between these two groups: miR-149, miR-10a, miR-20b, miR-30a-3p and miR-342-5p. All five microRNAs were down-regulated in tumors from patients with early recurrence. We show here that the 5-miRNA signature defines a high-risk group of patients with shorter relapse-free survival and has predictive value to discriminate non-relapsing versus early-relapsing patients (AUC = 0.993, p-value<0.05). Network analysis based on miRNA-target interactions curated by public databases suggests that down-regulation of the 5-miRNA signature in the subset of early-relapsing tumors would result in an overall increased proliferative and angiogenic capacity. In summary, we have identified a set of recurrence-related microRNAs with potential prognostic value to identify patients who will likely develop metastasis early after primary breast surgery.

Sma3s: a three-step modular annotator for large sequence datasets

Automatic sequence annotation is an essential component of modern 'omics' studies, which aim to extract information from large collections of sequence data. Most existing tools use sequence homology to establish evolutionary relationships and assign putative functions to sequences. However, it can be difficult to define a similarity threshold that achieves sufficient coverage without sacrificing annotation quality. Defining the correct configuration is critical and can be challenging for non-specialist users. Thus, the development of robust automatic annotation techniques that generate high-quality annotations without needing expert knowledge would be very valuable for the research community. We present Sma3s, a tool for automatically annotating very large collections of biological sequences from any kind of gene library or genome. Sma3s is composed of three modules that progressively annotate query sequences using either: (i) very similar homologues, (ii) orthologous sequences or (iii) terms enriched in groups of homologous sequences. We trained the system using several random sets of known sequences, demonstrating average sensitivity and specificity values of ~85%. In conclusion, Sma3s is a versatile tool for high-throughput annotation of a wide variety of sequence datasets that outperforms the accuracy of other well-established annotation algorithms, and it can enrich existing database annotations and uncover previously hidden features. Importantly, Sma3s has already been used in the functional annotation of two published transcriptomes.

Robust gene signatures from microarray data using genetic algorithms enriched with biological pathway keywords

Genetic algorithms are widely used in the estimation of expression profiles from microarrays data. However, these techniques are unable to produce stable and robust solutions suitable to use in clinical and biomedical studies. This paper presents a novel two-stage evolutionary strategy for gene feature selection combining the genetic algorithm with biological information extracted from the KEGG database. A comparative study is carried out over public data from three different types of cancer (leukemia, lung cancer and prostate cancer). Even though the analyses only use features having KEGG information, the results demonstrate that this two-stage evolutionary strategy increased the consistency, robustness and accuracy of a blind discrimination among relapsed and healthy individuals. Therefore, this approach could facilitate the definition of gene signatures for the clinical prognosis and diagnostic of cancer diseases in a near future. Additionally, it could also be used for biological knowledge discovery about the studied disease.

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.

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.

Evidence for an operative glutamine translocator in chloroplasts from maritime pine (Pinus pinaster Ait.) cotyledons

In higher plants, ammonium is assimilated into amino acids through the glutamine synthetase (GS)/glutamate synthase (GOGAT) cycle. This metabolic cycle is distributed in different cellular compartments in conifer seedlings: glutamine synthesis occurs in the cytosol and glutamate synthesis within the chloroplast. A method for preparing intact chloroplasts of pine cotyledons is presented with the aim of identifying a glutamine-glutamate translocator. Glutamine-glutamate exchange has been studied using the double silicone layer system, suggesting the existence of a translocator that imports glutamine into the chloroplast and exports glutamate to the cytoplasm. The translocator identified is specific for glutamine and glutamate, and the kinetic constants for both substrates indicate that it is unsaturated at intracellular concentrations. Thus, the experimental evidence obtained supports the model of the GS/GOGAT cycle in developing pine seedlings that accounts for the stoichiometric balance of metabolites. As a result, the efficient assimilation of free ammonia produced by photorespiration, nitrate reduction, storage protein mobilisation, phenylpropanoid pathway or S-adenosylmethionine synthesis is guaranteed.

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/.

AlignMiner: a Web-based tool for detection of divergent regions in multiple sequence alignments of conserved sequences

Background

Multiple sequence alignments are used to study gene or protein function, phylogenetic relations, genome evolution hypotheses and even gene polymorphisms. Virtually without exception, all available tools focus on conserved segments or residues. Small divergent regions, however, are biologically important for specific quantitative polymerase chain reaction, genotyping, molecular markers and preparation of specific antibodies, and yet have received little attention. As a consequence, they must be selected empirically by the researcher. AlignMiner has been developed to fill this gap in bioinformatic analyses.

Results

AlignMiner is a Web-based application for detection of conserved and divergent regions in alignments of conserved sequences, focusing particularly on divergence. It accepts alignments (protein or nucleic acid) obtained using any of a variety of algorithms, which does not appear to have a significant impact on the final results. AlignMiner uses different scoring methods for assessing conserved/divergent regions, Entropy being the method that provides the highest number of regions with the greatest length, and Weighted being the most restrictive. Conserved/divergent regions can be generated either with respect to the consensus sequence or to one master sequence. The resulting data are presented in a graphical interface developed in AJAX, which provides remarkable user interaction capabilities. Users do not need to wait until execution is complete and can.even inspect their results on a different computer. Data can be downloaded onto a user disk, in standard formats. In silico and experimental proof-of-concept cases have shown that AlignMiner can be successfully used to designing specific polymerase chain reaction primers as well as potential epitopes for antibodies. Primer design is assisted by a module that deploys several oligonucleotide parameters for designing primers "on the fly".

Conclusions

AlignMiner can be used to reliably detect divergent regions via several scoring methods that provide different levels of selectivity. Its predictions have been verified by experimental means. Hence, it is expected that its usage will save researchers' time and ensure an objective selection of the best-possible divergent region when closely related sequences are analysed. AlignMiner is freely available at http://www.scbi.uma.es/alignminer.

PreP+07: improvements of a user friendly tool to preprocess and analyse microarray data

Background

Nowadays, microarray gene expression analysis is a widely used technology that scientists handle but whose final interpretation usually requires the participation of a specialist. The need for this participation is due to the requirement of some background in statistics that most users lack or have a very vague notion of. Moreover, programming skills could also be essential to analyse these data. An interactive, easy to use application seems therefore necessary to help researchers to extract full information from data and analyse them in a simple, powerful and confident way.

Results

PreP+07 is a standalone Windows XP application that presents a friendly interface for spot filtration, inter- and intra-slide normalization, duplicate resolution, dye-swapping, error removal and statistical analyses. Additionally, it contains two unique implementation of the procedures – double scan and Supervised Lowess-, a complete set of graphical representations – MA plot, RG plot, QQ plot, PP plot, PN plot – and can deal with many data formats, such as tabulated text, GenePix GPR and ArrayPRO. PreP+07 performance has been compared with the equivalent functions in Bioconductor using a tomato chip with 13056 spots. The number of differentially expressed genes considering p-values coming from the PreP+07 and Bioconductor Limma packages were statistically identical when the data set was only normalized; however, a slight variability was appreciated when the data was both normalized and scaled.

Conclusion

PreP+07 implementation provides a high degree of freedom in selecting and organizing a small set of widely used data processing protocols, and can handle many data formats. Its reliability has been proven so that a laboratory researcher can afford a statistical pre-processing of his/her microarray results and obtain a list of differentially expressed genes using PreP+07 without any programming skills. All of this gives support to scientists that have been using previous PreP releases since its first version in 2003.

Interoperability with Moby 1.0--it's better than sharing your toothbrush!

The BioMoby project was initiated in 2001 from within the model organism database community. It aimed to standardize methodologies to facilitate information exchange and access to analytical resources, using a consensus driven approach. Six years later, the BioMoby development community is pleased to announce the release of the 1.0 version of the interoperability framework, registry Application Programming Interface and supporting Perl and Java code-bases. Together, these provide interoperable access to over 1400 bioinformatics resources worldwide through the BioMoby platform, and this number continues to grow. Here we highlight and discuss the features of BioMoby that make it distinct from other Semantic Web Service and interoperability initiatives, and that have been instrumental to its deployment and use by a wide community of bioinformatics service providers. The standard, client software, and supporting code libraries are all freely available at http://www.biomoby.org/.

Genetic Correction of Mitochondrial Diseases: Using the Natural Migration of Mitochondrial Genes to the Nucleus in Chlorophyte Algae as a Model System

Mitochondrial diseases display great diversity in clinical symptoms and biochemical characteristics. Although mtDNA mutations have been identified in many patients, there are currently no effective treatments. A number of human diseases result from mutations in mtDNA-encoded proteins, a group of proteins that are hydrophobic and have multiple membrane-spanning regions. One method that has great potential for overcoming the pathogenic consequences of these mutations is to place a wild-type copy of the affected gene in the nucleus, and target the expressed protein to the mitochondrion to function in place of the defective protein. Several respiratory chain subunit genes, which are typically mtDNA encoded, are nucleus encoded in the chlorophyte algae Chlamydomonas reinhardtii and Polytomella sp. Analysis of these genes has revealed adaptations that facilitated their expression from the nucleus. The nucleus-encoded proteins exhibited diminished physical constraints for import as compared to their mtDNA-encoded homologues. The hydrophobicity of the nucleus-encoded proteins is diminished in those regions that are not involved in subunit-subunit interactions or that contain amino acids critical for enzymatic reactions of the proteins. In addition, these proteins have unusually large mitochondrial targeting sequences. Information derived from these studies should be applicable toward the development of genetic therapies for human diseases resulting from mutations in mtDNA-encoded polypeptides.

Cytochrome f and subunit IV, two essential components of the photosynthetic bf complex typically encoded in the chloroplast genome, are nucleus-encoded in Euglena gracilis

The photosynthetic protist Euglena gracilis contains chloroplasts surrounded by three membranes which arise from secondary endosymbiosis. The genes petA and petD, encoding cytochrome f and subunit IV of the cytochrome bf complex, normally present in chloroplast genomes, are lacking from the chloroplast DNA (cpDNA) of E. gracilis. The bf complex of E. gracilis was isolated, and the identities of cytochrome f and subunit IV were established immunochemically, by heme-specific staining, and by Edman degradation. Based on N-terminal and conserved internal protein sequences, primers were designed and used for PCR gene amplification and cDNA sequencing. The complete sequence of the petA cDNA and the partial sequence of the petD cDNA from E. gracilis are described. Evidence is provided that in this protist, the petA and petD genes have migrated from the chloroplast to the nucleus. Both genes exhibit a typical nuclear codon usage, clearly distinct from the usage of chloroplast genes. The petA gene encodes an atypical cytochrome f, with a unique insertion of 62 residues not present in other f-type cytochromes. The petA gene also acquired a region that encodes a large tripartite chloroplast transit peptide (CTP), which is thought to allow the import of apocytochrome f through the three-membrane envelope of E. gracilis chloroplasts. This is the first description of petA and petD genes that are nucleus-localized.

Molecular and enzymatic analysis of ammonium assimilation in woody plants

Ammonium is assimilated into amino acids through the sequential action of glutamine synthetase (GS) and glutamate synthase (GOGAT) enzymes. This metabolic pathway is driven by energy, reducing power and requires the net supply of 2-oxoglutarate that can be provided by the reaction catalysed by isocitrate dehydrogenase (IDH). Most studies on the biochemistry and molecular biology of N-assimilating enzymes have been carried out on annual plant species and the available information on woody models is far more limited. This is in spite of their economic and ecological importance and the fact that nitrogen is a common limiting factor for tree growth. GS, GOGAT and IDH enzymes have been purified from several woody species and their kinetic and molecular properties determined. A number of cDNA clones have also been isolated and characterized. Although the enzymes are remarkably well conserved along the evolutionary scale, major differences have been found in their compartmentation within the cell between angiosperms and conifers, suggesting possible adaptations to specific functional roles. The analysis of the gene expression patterns in a variety of biological situations such as changes in N nutrition, development, biotic or abiotic stresses and senescence, suggest that cytosolic GS plays a central and pivotal role in ammonium assimilation and metabolism in woody plants. The modification of N assimilation efficiency has been recently approached in trees by overexpression of a cytosolic pine GS in poplar. The results obtained, suggest that an increase in cytosolic GS might lead to a global effect on the synthesis of nitrogenous compounds in the leaves, with enhanced vegetative growth of transgenic trees. All these data suggest that manipulation of cytosolic GS may have consequences for plant growth and biomass production.