iWhale: a computational pipeline based on Docker and SCons for detection and annotation of somatic variants in cancer WES data.

Whole exome sequencing (WES) is a powerful approach for discovering sequence variants in cancer cells but its time effectiveness is limited by the complexity and issues of WES data analysis. Here we present iWhale, a customizable pipeline based on Docker and SCons, reliably detecting somatic variants by three complementary callers (MuTect2, Strelka2 and VarScan2). The results are combined to obtain a single variant call format file for each sample and variants are annotated by integrating a wide range of information extracted from several reference databases, ultimately allowing variant and gene prioritization according to different criteria. iWhale allows users to conduct a complex series of WES analyses with a powerful yet customizable and easy-to-use tool, running on most operating systems (macOs, GNU/Linux and Windows). iWhale code is freely available at https://github.com/alexcoppe/iWhale and the docker image is downloadable from https://hub.docker.com/r/alexcoppe/iwhale.

Somatic mutations activating Wiskott-Aldrich syndrome protein concomitant with RAS pathway mutations in juvenile myelomonocytic leukemia patients.

The WAS gene product is expressed exclusively in the cytoplasm of hematopoietic cells and constitutional genetic abrogation of WASP leads to Wiskott-Aldrich syndrome (WAS). Moreover, mutational activation of WASP has been associated with X-linked neutropenia. Although studies reported that patients with constitutional WAS mutations affecting functional WASP expression may present juvenile myelomonocytic leukemia (JMML)-like features, confounding differential diagnosis above all in the copresence of mutated RAS, an activating somatic mutation of WASP has not been previously described in JMML patients. In our ongoing studies on JMML genomics, we at first detected a somatic WAS mutation in a major clone found at two consecutive relapses in one of two twins with JMML. Both twins were treated with hematopoietic stem cell transplantation after diagnosis of JMML. The somatic WAS mutation detected here displayed an activating WASP phenotype. Screening of 46 sporadic JMML patients at disease onset for mutations in the same PBD domain of WAS revealed two additional singleton patients carrying minor mutated clones. This is the first study to associate somatically acquired WASP mutations with a hematopoietic malignancy and increases insight in the complexity of the genomic landscape of JMML that shows low recurrent mutations concomitant with general hyperactivation of RAS pathway signaling.

Somatic mutations in specific and connected subpathways are associated with short neuroblastoma patients' survival and indicate proteins targetable at onset of disease.

Neuroblastoma (NB) is an embryonic malignancy of the sympathetic nervous system with heterogeneous biological, morphological, genetic and clinical characteristics. Although genomic studies revealed the specific biological features of NB pathogenesis useful for new therapeutic approaches, the improvement of high-risk (HR)-NB patients overall survival remains unsatisfactory. To further clarify the biological basis of disease aggressiveness, we used whole-exome sequencing to examine the genomic landscape of HR-NB patients at stage M with short survival (SS) and long survival (LS). Only a few genes, including SMARCA4, SMO, ZNF44 and CHD2, were recurrently and specifically mutated in the SS group, confirming the low recurrence of common mutations in this tumor. A systems biology approach revealed that in the two patient groups, mutations occurred in different pathways. Mutated genes (ARHGEF11, CACNA1G, FGF4, PTPRA, PTK2, ANK3, SMO, NTNG2, VCL and NID2) regulate the MAPK pathway associated with the organization of the extracellular matrix, cell motility through PTK2 signaling and matrix metalloproteinase activity. Moreover, we detected mutations in LAMA2, PTK2, LAMA4, and MMP14 genes, impairing MET signaling, in SFI1 and CHD2 involved in centrosome maturation and chromosome remodeling, in AK7 and SPTLC2, which regulate the metabolism of nucleotides and lipoproteins, and in NALCN, SLC12A1, SLC9A9, which are involved in the transport of small molecules. Notably, connected networks of somatically mutated genes specific for SS patients were identified. The detection of mutated genes present at the onset of disease may help to address an early treatment of HR-NB patients using FDA-approved compounds targeting the deregulated pathways.

MAGIA²: from miRNA and genes expression data integrative analysis to microRNA-transcription factor mixed regulatory circuits (2012 update).

MAGIA(2) (http://gencomp.bio.unipd.it/magia2) is an update, extension and evolution of the MAGIA web tool. It is dedicated to the integrated analysis of in silico target prediction, microRNA (miRNA) and gene expression data for the reconstruction of post-transcriptional regulatory networks. miRNAs are fundamental post-transcriptional regulators of several key biological and pathological processes. As miRNAs act prevalently through target degradation, their expression profiles are expected to be inversely correlated to those of the target genes. Low specificity of target prediction algorithms makes integration approaches an interesting solution for target prediction refinement. MAGIA(2) performs this integrative approach supporting different association measures, multiple organisms and almost all target predictions algorithms. Nevertheless, miRNAs activity should be viewed as part of a more complex scenario where regulatory elements and their interactors generate a highly connected network and where gene expression profiles are the result of different levels of regulation. The updated MAGIA(2) tries to dissect this complexity by reconstructing mixed regulatory circuits involving either miRNA or transcription factor (TF) as regulators. Two types of circuits are identified: (i) a TF that regulates both a miRNA and its target and (ii) a miRNA that regulates both a TF and its target.

Transcriptome sequencing and de novo annotation of the critically endangered Adriatic sturgeon.

BACKGROUND: Sturgeons are a group of Condrostean fish with very high evolutionary, economical and conservation interest. The eggs of these living fossils represent one of the most high prized foods of animal origin. The intense fishing pressure on wild stocks to harvest caviar has caused in the last decades a dramatic decline of their distribution and abundance leading the International Union for Conservation of Nature to list them as the more endangered group of species. As a direct consequence, world-wide efforts have been made to develop sturgeon aquaculture programmes for caviar production. In this context, the characterization of the genes involved in sex determination could provide relevant information for the selective farming of the more profitable females. RESULTS: The 454 sequencing of two cDNA libraries from the gonads and brain of one male and one female full-sib A. naccarii, yielded 182,066 and 167,776 reads respectively, which, after strict quality control, were iterative assembled into more than 55,000 high quality ESTs. The average per-base coverage reached by assembling the two libraries was 4X. The multi-step annotation process resulted in 16% successfully annotated sequences with GO terms. We screened the transcriptome for 32 sex-related genes and highlighted 7 genes that are potentially specifically expressed, 5 in male and 2 in females, at the first life stage at which sex is histologically identifiable. In addition we identified 21,791 putative EST-linked SNPs and 5,295 SSRs. CONCLUSIONS: This study represents the first large massive release of sturgeon transcriptome information that we organized into the public database AnaccariiBase, which is freely available at http://compgen.bio.unipd.it/anaccariibase/. This transcriptomic data represents an important source of information for further studies on sturgeon species. The hundreds of putative EST-linked molecular makers discovered in this study will be invaluable for sturgeon reintroduction and breeding programs.

mRNA-Seq and microarray development for the Grooved Carpet shell clam, Ruditapes decussatus: a functional approach to unravel host-parasite interaction.

BACKGROUND: The Grooved Carpet shell clam Ruditapes decussatus is the autochthonous European clam and the most appreciated from a gastronomic and economic point of view. The production is in decline due to several factors such as Perkinsiosis and habitat invasion and competition by the introduced exotic species, the manila clam Ruditapes philippinarum. After we sequenced R. decussatus transcriptome we have designed an oligo microarray capable of contributing to provide some clues on molecular response of the clam to Perkinsiosis. RESULTS: A database consisting of 41,119 unique transcripts was constructed, of which 12,479 (30.3%) were annotated by similarity. An oligo-DNA microarray platform was then designed and applied to profile gene expression in R. decussatus heavily infected by Perkinsus olseni. Functional annotation of differentially expressed genes between those two conditionswas performed by gene set enrichment analysis. As expected, microarrays unveil genes related with stress/infectious agents such as hydrolases, proteases and others. The extensive role of innate immune system was also analyzed and effect of parasitosis upon expression of important molecules such as lectins reviewed. CONCLUSIONS: This study represents a first attempt to characterize Ruditapes decussatus transcriptome, an important marine resource for the European aquaculture. The trancriptome sequencing and consequent annotation will increase the available tools and resources for this specie, introducing the possibility of high throughput experiments such as microarrays analysis. In this specific case microarray approach was used to unveil some important aspects of host-parasite interaction between the Carpet shell clam and Perkinsus, two non-model species, highlighting some genes associated with this interaction. Ample information was obtained to identify biological processes significantly enriched among differentially expressed genes in Perkinsus infected versus non-infected gills. An overview on the genes related with the immune system on R. decussatus transcriptome is also reported.

Multiomic analysis of HER2-enriched and AR-positive breast carcinoma with apocrine differentiation and an oligometastatic course: a case report.

Breast carcinoma is the most prevalent cancer among women globally. It has variable clinical courses depending on the stage and clinical-biological features. This case report describes a 56-year-old female with invasive breast cancer without estrogen or progesterone receptor expression, with apocrine differentiation, and with no germline variants in the BRCA1 and BRCA2 genes. Throughout the clinical course, the patient exhibited discordant results for HER2 in immunohistochemistry and in situ hybridization. During the second relapse, the disease displayed apocrine microscopic features. The tumor underwent analysis for the androgen receptor, GCDFP-15, RNA-seq, and whole-genome sequencing (WGS) to identify the breast cancer subtype and to characterize the cancer genome. Our bioinformatic analysis revealed 20,323 somatic SNV/Indels, including five mutations in cancer-related genes that are believed to be responsible for the tumor's development. Two of these mutations were found in the PIK3CA and TP53 genes. Furthermore, the tumor tissue exhibited large copy number alterations to the chromosomes, which could impact gene expression through complex mechanisms and contribute to the tumor phenotype. Clustering algorithms applied on RNA-sequencing data categorized this cancer as a HER2+ subtype. The second-line capecitabine chemotherapy treatment is ongoing, and the patient is responding well. Bioinformatic results support the current treatment decision and open the way to further treatments.

Surviving in a toxic world: transcriptomics and gene expression profiling in response to environmental pollution in the critically endangered European eel.

BACKGROUND: Genomic and transcriptomic approaches have the potential for unveiling the genome-wide response to environmental perturbations. The abundance of the catadromous European eel (Anguilla anguilla) stock has been declining since the 1980s probably due to a combination of anthropogenic and climatic factors. In this paper, we explore the transcriptomic dynamics between individuals from high (river Tiber, Italy) and low pollution (lake Bolsena, Italy) environments, which were measured for 36 PCBs, several organochlorine pesticides and brominated flame retardants and nine metals. RESULTS: To this end, we first (i) updated the European eel transcriptome using deep sequencing data with a total of 640,040 reads assembled into 44,896 contigs (Eeelbase release 2.0), and (ii) developed a transcriptomic platform for global gene expression profiling in the critically endangered European eel of about 15,000 annotated contigs, which was applied to detect differentially expressed genes between polluted sites. Several detoxification genes related to metabolism of pollutants were upregulated in the highly polluted site, including genes that take part in phase I of the xenobiotic metabolism (CYP3A), phase II (glutathione-S-transferase) and oxidative stress (glutathione peroxidase). In addition, key genes in the mitochondrial respiratory chain and oxidative phosphorylation were down-regulated at the Tiber site relative to the Bolsena site. CONCLUSIONS: Together with the induced high expression of detoxification genes, the suggested lowered expression of genes supposedly involved in metabolism suggests that pollution may also be associated with decreased respiratory and energy production.

Impact of probe annotation on the integration of miRNA-mRNA expression profiles for miRNA target detection.

MicroRNAs (miRNAs) are small non-coding RNAs that mediate gene expression at the post-transcriptional and translational levels by an imperfect binding to target mRNA 3'UTR regions. While the ab-initio computational prediction of miRNA-mRNA interactions still poses significant challenges, it is possible to overcome some of its limitations by carefully integrating into the analysis the paired expression profiles of miRNAs and mRNAs. In this work, we show how the choice of a proper probe annotation for microarray platforms is an essential requirement to achieve good sensitivity in the identification of miRNA-mRNA interactions. We compare the results obtained from the analysis of the same expression profiles using both gene and transcript based custom CDFs that we have developed for a number of different annotations (ENSEMBL, RefSeq, AceView). In all cases, transcript-based annotations clearly improve the effectiveness of data integration and thus provide a more reliable confirmation of computationally predicted miRNA-mRNA interactions.

MAGIA, a web-based tool for miRNA and Genes Integrated Analysis.

MAGIA (miRNA and genes integrated analysis) is a novel web tool for the integrative analysis of target predictions, miRNA and gene expression data. MAGIA is divided into two parts: the query section allows the user to retrieve and browse updated miRNA target predictions computed with a number of different algorithms (PITA, miRanda and Target Scan) and Boolean combinations thereof. The analysis section comprises a multistep procedure for (i) direct integration through different functional measures (parametric and non-parametric correlation indexes, a variational Bayesian model, mutual information and a meta-analysis approach based on P-value combination) of mRNA and miRNA expression data, (ii) construction of bipartite regulatory network of the best miRNA and mRNA putative interactions and (iii) retrieval of information available in several public databases of genes, miRNAs and diseases and via scientific literature text-mining. MAGIA is freely available for Academic users at http://gencomp.bio.unipd.it/magia.

TRAM (Transcriptome Mapper): database-driven creation and analysis of transcriptome maps from multiple sources.

BACKGROUND: Several tools have been developed to perform global gene expression profile data analysis, to search for specific chromosomal regions whose features meet defined criteria as well as to study neighbouring gene expression. However, most of these tools are tailored for a specific use in a particular context (e.g. they are species-specific, or limited to a particular data format) and they typically accept only gene lists as input. RESULTS: TRAM (Transcriptome Mapper) is a new general tool that allows the simple generation and analysis of quantitative transcriptome maps, starting from any source listing gene expression values for a given gene set (e.g. expression microarrays), implemented as a relational database. It includes a parser able to assign univocal and updated gene symbols to gene identifiers from different data sources. Moreover, TRAM is able to perform intra-sample and inter-sample data normalization, including an original variant of quantile normalization (scaled quantile), useful to normalize data from platforms with highly different numbers of investigated genes. When in 'Map' mode, the software generates a quantitative representation of the transcriptome of a sample (or of a pool of samples) and identifies if segments of defined lengths are over/under-expressed compared to the desired threshold. When in 'Cluster' mode, the software searches for a set of over/under-expressed consecutive genes. Statistical significance for all results is calculated with respect to genes localized on the same chromosome or to all genome genes. Transcriptome maps, showing differential expression between two sample groups, relative to two different biological conditions, may be easily generated. We present the results of a biological model test, based on a meta-analysis comparison between a sample pool of human CD34+ hematopoietic progenitor cells and a sample pool of megakaryocytic cells. Biologically relevant chromosomal segments and gene clusters with differential expression during the differentiation toward megakaryocyte were identified. CONCLUSIONS: TRAM is designed to create, and statistically analyze, quantitative transcriptome maps, based on gene expression data from multiple sources. The release includes FileMaker Pro database management runtime application and it is freely available at http://apollo11.isto.unibo.it/software/, along with preconfigured implementations for mapping of human, mouse and zebrafish transcriptomes.

Transcriptome sequencing and microarray development for the Manila clam, Ruditapes philippinarum: genomic tools for environmental monitoring.

BACKGROUND: The Manila clam, Ruditapes philippinarum, is one of the major aquaculture species in the world and a potential sentinel organism for monitoring the status of marine ecosystems. However, genomic resources for R. philippinarum are still extremely limited. Global analysis of gene expression profiles is increasingly used to evaluate the biological effects of various environmental stressors on aquatic animals under either artificial conditions or in the wild. Here, we report on the development of a transcriptomic platform for global gene expression profiling in the Manila clam. RESULTS: A normalized cDNA library representing a mixture of adult tissues was sequenced using a ultra high-throughput sequencing technology (Roche 454). A database consisting of 32,606 unique transcripts was constructed, 9,747 (30%) of which could be annotated by similarity. An oligo-DNA microarray platform was designed and applied to profile gene expression of digestive gland and gills. Functional annotation of differentially expressed genes between different tissues was performed by enrichment analysis. Expression of Natural Antisense Transcripts (NAT) analysis was also performed and bi-directional transcription appears a common phenomenon in the R. philippinarum transcriptome. A preliminary study on clam samples collected in a highly polluted area of the Venice Lagoon demonstrated the applicability of genomic tools to environmental monitoring. CONCLUSIONS: The transcriptomic platform developed for the Manila clam confirmed the high level of reproducibility of current microarray technology. Next-generation sequencing provided a good representation of the clam transcriptome. Despite the known limitations in transcript annotation and sequence coverage for non model species, sufficient information was obtained to identify a large set of genes potentially involved in cellular response to environmental stress.

Motif discovery in promoters of genes co-localized and co-expressed during myeloid cells differentiation.

Genes co-expressed may be under similar promoter-based and/or position-based regulation. Although data on expression, position and function of human genes are available, their true integration still represents a challenge for computational biology, hampering the identification of regulatory mechanisms. We carried out an integrative analysis of genomic position, functional annotation and promoters of genes expressed in myeloid cells. Promoter analysis was conducted by a novel multi-step method for discovering putative regulatory elements, i.e. over-represented motifs, in a selected set of promoters, as compared with a background model. The combination of transcriptional, structural and functional data allowed the identification of sets of promoters pertaining to groups of genes co-expressed and co-localized in regions of the human genome. The application of motif discovery to 26 groups of genes co-expressed in myeloid cells differentiation and co-localized in the genome showed that there are more over-represented motifs in promoters of co-expressed and co-localized genes than in promoters of simply co-expressed genes (CEG). Motifs, which are similar to the binding sequences of known transcription factors, non-uniformly distributed along promoter sequences and/or occurring in highly co-expressed subset of genes were identified. Co-expressed and co-localized gene sets were grouped in two co-expressed genomic meta-regions, putatively representing functional domains of a high-level expression regulation.

Sequencing, de novo annotation and analysis of the first Anguilla anguilla transcriptome: EeelBase opens new perspectives for the study of the critically endangered European eel.

BACKGROUND: Once highly abundant, the European eel (Anguilla anguilla L.; Anguillidae; Teleostei) is considered to be critically endangered and on the verge of extinction, as the stock has declined by 90-99% since the 1980s. Yet, the species is poorly characterized at molecular level with little sequence information available in public databases. RESULTS: The first European eel transcriptome was obtained by 454 FLX Titanium sequencing of a normalized cDNA library, produced from a pool of 18 glass eels (juveniles) from the French Atlantic coast and two sites in the Mediterranean coast. Over 310,000 reads were assembled in a total of 19,631 transcribed contigs, with an average length of 531 nucleotides. Overall 36% of the contigs were annotated to known protein/nucleotide sequences and 35 putative miRNA identified. CONCLUSIONS: This study represents the first transcriptome analysis for a critically endangered species. EeelBase, a dedicated database of annotated transcriptome sequences of the European eel is freely available at http://compgen.bio.unipd.it/eeelbase. Considering the multiple factors potentially involved in the decline of the European eel, including anthropogenic factors such as pollution and human-introduced diseases, our results will provide a rich source of data to discover and identify new genes, characterize gene expression, as well as for identification of genetic markers scattered across the genome to be used in various applications.

Lack of Viral Load Within Chronic Lymphoproliferative Disorder of Natural Killer Cells: What Is Outside the Leukemic Clone?

Large granular lymphocyte leukemias (LGLL) are sustained by proliferating cytotoxic T cells or NK cells, as happens in Chronic Lymphoproliferative Disorder of Natural Killer cells (CLPD-NK), whose etiology is only partly understood. Different hypotheses have been proposed on the original events triggering NK cell hyperactivation and transformation, including a role of viral agents. In this perspective, we revise the lines of evidence that suggested a pathogenetic role in LGLL of the exposure to retroviruses and that identified Epstein Barr Virus (EBV) in other NK cell leukemias and lymphomas and focus on the contrasting data about the importance of viral agents in CLPD-NK. EBV was detected in aggressive NK leukemias but not in the indolent CLPD-NK, where seroreactivity against HTLV-1 retrovirus envelope BA21 protein antigens has been reported in patients, although lacking clear evidence of HTLV infection. We next present original results of whole exome sequencing data analysis that failed to identify viral sequences in CLPD-NK. We recently demonstrated that proliferating NK cells of patients harbor several somatic lesions likely contributing to sustain NK cell proliferation. Thus, we explore whether "neoantigens" similar to the BA21 antigen could be generated by aberrancies present in the leukemic clone. In light of the literature and new data, we evaluated the intriguing hypothesis that NK cell activation can be caused by retroviral agents located outside the hematopoietic compartment and on the possible mechanisms involved with the prospects of immunotherapy-based approaches to limit the growth of NK cells in CLPD-NK disease.

A high definition picture of somatic mutations in chronic lymphoproliferative disorder of natural killer cells.

The molecular pathogenesis of chronic lymphoproliferative disorder of natural killer (NK) cells (CLPD-NK) is poorly understood. Following the screening of 57 CLPD-NK patients, only five presented STAT3 mutations. WES profiling of 13 cases negative for STAT3/STAT5B mutations uncovered an average of 18 clonal, population rare and deleterious somatic variants per patient. The mutational landscape of CLPD-NK showed that most patients carry a heavy mutational burden, with major and subclonal deleterious mutations co-existing in the leukemic clone. Somatic mutations hit genes wired to cancer proliferation, survival, and migration pathways, in the first place Ras/MAPK, PI3K-AKT, in addition to JAK/STAT (PIK3R1 and PTK2). We confirmed variants with putative driver role of MAP10, MPZL1, RPS6KA1, SETD1B, TAOK2, TMEM127, and TNFRSF1A genes, and of genes linked to viral infections (DDX3X and RSF1) and DNA repair (PAXIP1). A truncating mutation of the epigenetic regulator TET2 and a variant likely abrogating PIK3R1-negative regulatory activity were validated. This study significantly furthered the view of the genes and pathways involved in CLPD-NK, indicated similarities with aggressive diseases of NK cells and detected mutated genes targetable by approved drugs, being a step forward to personalized precision medicine for CLPD-NK patients.

A-MADMAN: annotation-based microarray data meta-analysis tool.

BACKGROUND: Publicly available datasets of microarray gene expression signals represent an unprecedented opportunity for extracting genomic relevant information and validating biological hypotheses. However, the exploitation of this exceptionally rich mine of information is still hampered by the lack of appropriate computational tools, able to overcome the critical issues raised by meta-analysis. RESULTS: This work presents A-MADMAN, an open source web application which allows the retrieval, annotation, organization and meta-analysis of gene expression datasets obtained from Gene Expression Omnibus. A-MADMAN addresses and resolves several open issues in the meta-analysis of gene expression data. CONCLUSION: A-MADMAN allows i) the batch retrieval from Gene Expression Omnibus and the local organization of raw data files and of any related meta-information, ii) the re-annotation of samples to fix incomplete, or otherwise inadequate, metadata and to create user-defined batches of data, iii) the integrative analysis of data obtained from different Affymetrix platforms through custom chip definition files and meta-normalization. Software and documentation are available on-line at http://compgen.bio.unipd.it/bioinfo/amadman/.