We present a novel approach for the analysis of transcriptomics data that integrates functional annotation of gene sets with expression values in a multivariate fashion, and directly assesses the relation of functional features to a multivariate space of response phenotypical variables. Multivariate projection methods are used to obtain new correlated variables for a set of genes that share a given function. These new functional variables are then related to the response variables of interest. The analysis of the principal directions of the multivariate regression allows for the identification of gene function features correlated with the phenotype. Two different transcriptomics studies are used to illustrate the statistical and interpretative aspects of the methodology. We demonstrate the superiority of the proposed method over equivalent approaches.

Blastomere fate and embryonic genome activation (EGA) during human embryonic development are unsolved areas of high scientific and clinical interest. Forty-nine blastomeres from 5- to 8-cell human embryos have been investigated following an efficient single-cell cDNA amplification protocol to provide a template for high-density microarray analysis. The previously described markers, characteristic of Inner Cell Mass (ICM) (n = 120), stemness (n = 190) and Trophectoderm (TE) (n = 45), were analyzed, and a housekeeping pattern of 46 genes was established. All the human blastomeres from the 5- to 8-cell stage embryo displayed a common gene expression pattern corresponding to ICM markers (e.g., DDX3, FOXD3, LEFTY1, MYC, NANOG, POU5F1), stemness (e.g., POU5F1, DNMT3B, GABRB3, SOX2, ZFP42, TERT), and TE markers (e.g., GATA6, EOMES, CDX2, LHCGR). The EGA profile was also investigated between the 5-6- and 8-cell stage embryos, and compared to the blastocyst stage. Known genes (n = 92) such as depleted maternal transcripts (e.g., CCNA1, CCNB1, DPPA2) and embryo-specific activation (e.g., POU5F1, CDH1, DPPA4), as well as novel genes, were confirmed. In summary, the global single-cell cDNA amplification microarray analysis of the 5- to 8-cell stage human embryos reveals that blastomere fate is not committed to ICM or TE. Finally, new EGA features in human embryogenesis are presented.

Recent observations point towards the existence of a large number of neighborhoods composed of functionally-related gene modules that lie together in the genome. This local component in the distribution of the functionality across chromosomes is probably affecting the own chromosomal architecture by limiting the possibilities in which genes can be arranged and distributed across the genome. As a direct consequence of this fact it is therefore presumable that diseases such as cancer, harboring DNA copy number alterations (CNAs), will have a symptomatology strongly dependent on modules of functionally-related genes rather than on a unique "important" gene.

Phylemon 2.0 is a new release of the suite of web tools for molecular evolution, phylogenetics, phylogenomics and hypotheses testing. It has been designed as a response to the increasing demand of molecular sequence analyses for experts and non-expert users. Phylemon 2.0 has several unique features that differentiates it from other similar web resources: (i) it offers an integrated environment that enables evolutionary analyses, format conversion, file storage and edition of results; (ii) it suggests further analyses, thereby guiding the users through the web server; and (iii) it allows users to design and save phylogenetic pipelines to be used over multiple genes (phylogenomics). Altogether, Phylemon 2.0 integrates a suite of 30 tools covering sequence alignment reconstruction and trimming; tree reconstruction, visualization and manipulation; and evolutionary hypotheses testing.

Recent advances in whole-genome association studies (WGASs) for human cancer risk are beginning to provide the part lists of low-penetrance susceptibility genes. However, statistical analysis in these studies is complicated by the vast number of genetic variants examined and the weak effects observed, as a result of which constraints must be incorporated into the study design and analytical approach. In this scenario, biological attributes beyond the adjusted statistics generally receive little attention and, more importantly, the fundamental biological characteristics of low-penetrance susceptibility genes have yet to be determined.

BRCA1 is a key protein in cell network, involved in DNA repair pathways and cell cycle. Recently, the ENIGMA consortium has reported a high number of alternative splicing (AS) events at this locus in blood-derived samples. However, BRCA1 splicing pattern in breast tissue samples is unknown. Here, we provide an accurate description of BRCA1 splicing events distribution in breast tissue samples.

Molecular diagnosis of Inherited Retinal Dystrophies (IRD) has long been challenging due to the extensive clinical and genetic heterogeneity present in this group of disorders. Here, we describe the clinical application of an integrated next-generation sequencing approach to determine the underlying genetic defects in a Spanish family with a provisional clinical diagnosis of autosomal recessive Retinitis Pigmentosa (arRP).

Disease targeted sequencing is gaining importance as a powerful and cost-effective application of high throughput sequencing technologies to the diagnosis. However, the lack of proper tools to process the data hinders its extensive adoption. Here we present TEAM, an intuitive and easy-to-use web tool that fills the gap between the predicted mutations and the final diagnostic in targeted enrichment sequencing analysis. The tool searches for known diagnostic mutations, corresponding to a disease panel, among the predicted patient's variants. Diagnostic variants for the disease are taken from four databases of disease-related variants (HGMD-public, HUMSAVAR, ClinVar and COSMIC.) If no primary diagnostic variant is found, then a list of secondary findings that can help to establish a diagnostic is produced. TEAM also provides with an interface for the definition of and customization of panels, by means of which, genes and mutations can be added or discarded to adjust panel definitions. TEAM is freely available at: http://team.babelomics.org.

Most diagnostic laboratories are confronted with the increasing demand for molecular diagnosis from patients and families and the ever-increasing genetic heterogeneity of visual disorders. Concerning Retinal Dystrophies (RD), almost 200 causative genes have been reported to date, and most families carry private mutations. We aimed to approach RD genetic diagnosis using all the available genetic information to prioritize candidates for mutational screening, and then restrict the number of cases to be analyzed by massive sequencing. We constructed and optimized a comprehensive cosegregation RD-chip based on SNP genotyping and haplotype analysis. The RD-chip allows to genotype 768 selected SNPs (closely linked to 100 RD causative genes) in a single cost-, time-effective step. Full diagnosis was attained in 17/36 Spanish pedigrees, yielding 12 new and 12 previously reported mutations in 9 RD genes. The most frequently mutated genes were USH2A and CRB1. Notably, RD3-up to now only associated to Leber Congenital Amaurosis- was identified as causative of Retinitis Pigmentosa. The main assets of the RD-chip are: i) the robustness of the genetic information that underscores the most probable candidates, ii) the invaluable clues in cases of shared haplotypes, which are indicative of a common founder effect, and iii) the detection of extended haplotypes over closely mapping genes, which substantiates cosegregation, although the assumptions in which the genetic analysis is based could exceptionally lead astray. The combination of the genetic approach with whole exome sequencing (WES) greatly increases the diagnosis efficiency, and revealed novel mutations in USH2A and GUCY2D. Overall, the RD-chip diagnosis efficiency ranges from 16% in dominant, to 80% in consanguineous recessive pedigrees, with an average of 47%, well within the upper range of massive sequencing approaches, highlighting the validity of this time- and cost-effective approach whilst high-throughput methodologies become amenable for routine diagnosis in medium sized labs.

The auditory and vestibular organs form the inner ear and have a common developmental origin. Insulin like growth factor 1 (IGF-1) has a central role in the development of the cochlea and maintenance of hearing. Its deficiency causes sensorineural hearing loss in man and mice. During chicken early development, IGF-1 modulates neurogenesis of the cochleovestibular ganglion but no further studies have been conducted to explore the potential role of IGF-1 in the vestibular system. In this study we have compared the whole transcriptome of the vestibular organ from wild type and Igf1(-/-) mice at different developmental and postnatal times. RNA was prepared from E18.5, P15 and P90 vestibular organs of Igf1(-/-) and Igf1(+/+) mice and the transcriptome analysed in triplicates using Affymetrix(®) Mouse Gene 1.1 ST Array Plates. These plates are whole-transcript arrays that include probes to measure both messenger (mRNA) and long intergenic non-coding RNA transcripts (lincRNA), with a coverage of over 28 thousand coding transcripts and over 7 thousands non-coding transcripts. Given the complexity of the data we used two different methods VSN-RMA and mmBGX to analyse and compare the data. This is to better evaluate the number of false positives and to quantify uncertainty of low signals. We identified a number of differentially expressed genes that we described using functional analysis and validated using RT-qPCR. The morphology of the vestibular organ did not show differences between genotypes and no evident alterations were observed in the vestibular sensory areas of the null mice. However, well-defined cellular alterations were found in the vestibular neurons with respect their number and size. Although these mice did not show a dramatic vestibular phenotype, we conducted a functional analysis on differentially expressed genes between genotypes and across time. This was with the aim to identify new pathways that are involved in the development of the vestibular organ as well as pathways that maybe affected by the lack of IGF-1 and be associated to the morphological changes of the vestibular neurons that we observed in the Igf1(-/-) mice.

Models of plant-virus interaction assume that the ability of a virus to infect a host genotype depends on the matching between virulence and resistance genes. Recently, we evolved tobacco etch potyvirus (TEV) lineages on different ecotypes of Arabidopsis thaliana, and found that some ecotypes selected for specialist viruses whereas others selected for generalists. Here we sought to evaluate the transcriptomic basis of such relationships. We have characterized the transcriptomic responses of five ecotypes infected with the ancestral and evolved viruses. Genes and functional categories differentially expressed by plants infected with local TEV isolates were identified, showing heterogeneous responses among ecotypes, although significant parallelism existed among lineages evolved in the same ecotype. Although genes involved in immune responses were altered upon infection, other functional groups were also pervasively over-represented, suggesting that plant resistance genes were not the only drivers of viral adaptation. Finally, the transcriptomic consequences of infection with the generalist and specialist lineages were compared. Whilst the generalist induced very similar perturbations in the transcriptomes of the different ecotypes, the perturbations induced by the specialist were divergent. Plant defense mechanisms were activated when the infecting virus was specialist but they were down-regulated when infecting with generalist.

Mesenchymal stem cells are often transplanted into inflammatory environments where they are able to survive and modulate host immune responses through a poorly understood mechanism. In this paper we analyzed the responses of MSC to IL-1β: a representative inflammatory mediator. Microarray analysis of MSC treated with IL-1β revealed that this cytokine activateds a set of genes related to biological processes such as cell survival, cell migration, cell adhesion, chemokine production, induction of angiogenesis and modulation of the immune response. Further more detailed analysis by real-time PCR and functional assays revealed that IL-1β mainly increaseds the production of chemokines such as CCL5, CCL20, CXCL1, CXCL3, CXCL5, CXCL6, CXCL10, CXCL11 and CX(3)CL1, interleukins IL-6, IL-8, IL23A, IL32, Toll-like receptors TLR2, TLR4, CLDN1, metalloproteins MMP1 and MMP3, growth factors CSF2 and TNF-α, together with adhesion molecules ICAM1 and ICAM4. Functional analysis of MSC proliferation, migration and adhesion to extracellular matrix components revealed that IL-1β did not affect proliferation but also served to induce the secretion of trophic factors and adhesion to ECM components such as collagen and laminin. IL-1β treatment enhanced the ability of MSC to recruit monocytes and granulocytes in vitro. Blockade of NF-κβ transcription factor activation with IκB kinase beta (IKKβ) shRNA impaired MSC migration, adhesion and leucocyte recruitment, induced by IL-1β demonstrating that NF-κB pathway is an important downstream regulator of these responses. These findings are relevant to understanding the biological responses of MSC to inflammatory environments.

Finding gene associations in rare diseases is frequently hampered by the reduced numbers of patients accessible. Conventional gene-based association tests rely on the availability of large cohorts, which constitutes a serious limitation for its application in this scenario. To overcome this problem we have used here a combined strategy in which a pathway-based analysis (PBA) has been initially conducted to prioritize candidate genes in a Spanish cohort of 53 trios of short-segment Hirschsprung's disease. Candidate genes have been further validated in an independent population of 106 trios. The study revealed a strong association of 11 gene ontology (GO) modules related to signal transduction and its regulation, enteric nervous system (ENS) formation and other HSCR-related processes. Among the preselected candidates, a total of 4 loci, RASGEF1A, IQGAP2, DLC1 and CHRNA7, related to signal transduction and migration processes, were found to be significantly associated to HSCR. Network analysis also confirms their involvement in the network of already known disease genes. This approach, based on the study of functionally-related gene sets, requires of lower sample sizes and opens new opportunities for the study of rare diseases.

Retinitis Pigmentosa (RP) is a heterogeneous group of inherited retinal dystrophies characterised ultimately by the loss of photoreceptor cells. RP is the leading cause of visual loss in individuals younger than 60 years, with a prevalence of about 1 in 4000. The molecular genetic diagnosis of autosomal recessive RP (arRP) is challenging due to the large genetic and clinical heterogeneity. Traditional methods for sequencing arRP genes are often laborious and not easily available and a screening technique that enables the rapid detection of the genetic cause would be very helpful in the clinical practice. The goal of this study was to develop and apply microarray-based resequencing technology capable of detecting both known and novel mutations on a single high-throughput platform. Hence, the coding regions and exon/intron boundaries of 16 arRP genes were resequenced using microarrays in 102 Spanish patients with clinical diagnosis of arRP. All the detected variations were confirmed by direct sequencing and potential pathogenicity was assessed by functional predictions and frequency in controls. For validation purposes 4 positive controls for variants consisting of previously identified changes were hybridized on the array. As a result of the screening, we detected 44 variants, of which 15 are very likely pathogenic detected in 14 arRP families (14%). Finally, the design of this array can easily be transformed in an equivalent diagnostic system based on targeted enrichment followed by next generation sequencing.

The issue of whether coelomates form a single clade, the Coelomata, or whether all animals that moult an exoskeleton (such as the coelomate arthropods and the pseudocoelomate nematodes) form a distinct clade, the Ecdysozoa, is the most puzzling issue in animal systematics and a major open-ended subject in evolutionary biology. Previous single-gene and genome-scale analyses designed to resolve the issue have produced contradictory results. Here we present the first genome-scale phylogenetic evidence that strongly supports the Ecdysozoa hypothesis.

We have developed a web tool, PupasView, for the selection of single nucleotide polymorphisms (SNPs) with potential phenotypic effect. PupasView constitutes an interactive environment in which functional information and population frequency data can be used as sequential filters over linkage disequilibrium parameters to obtain a final list of SNPs optimal for genotyping purposes. PupasView is the first resource that integrates phenotypic effects caused by SNPs at both the translational and the transcriptional level. PupasView retrieves SNPs that could affect conserved regions that the cellular machinery uses for the correct processing of genes (intron/exon boundaries or exonic splicing enhancers), predicted transcription factor binding sites and changes in amino acids in the proteins for which a putative pathological effect is calculated. The program uses the mapping of SNPs in the genome provided by Ensembl. PupasView will be of much help in studies of multifactorial disorders, where the use of functional SNPs will increase the sensitivity of the identification of the genes responsible for the disease. The PupasView web interface is accessible through http://pupasview.ochoa.fib.es and through http://www.pupasnp.org.

Phylemon is an online platform for phylogenetic and evolutionary analyses of molecular sequence data. It has been developed as a web server that integrates a suite of different tools selected among the most popular stand-alone programs in phylogenetic and evolutionary analysis. It has been conceived as a natural response to the increasing demand of data analysis of many experimental scientists wishing to add a molecular evolution and phylogenetics insight into their research. Tools included in Phylemon cover a wide yet selected range of programs: from the most basic for multiple sequence alignment to elaborate statistical methods of phylogenetic reconstruction including methods for evolutionary rates analyses and molecular adaptation. Phylemon has several features that differentiates it from other resources: (i) It offers an integrated environment that enables the direct concatenation of evolutionary analyses, the storage of results and handles required data format conversions, (ii) Once an outfile is produced, Phylemon suggests the next possible analyses, thus guiding the user and facilitating the integration of multi-step analyses, and (iii) users can define and save complete pipelines for specific phylogenetic analysis to be automatically used on many genes in subsequent sessions or multiple genes in a single session (phylogenomics). The Phylemon web server is available at http://phylemon.bioinfo.cipf.es.

In silico drug prescription tools for precision cancer medicine can match molecular alterations with tailored candidate treatments. These methodologies require large and well-annotated datasets to systematically evaluate their performance, but this is currently constrained by the lack of complete patient clinicopathological data. Moreover, in silico drug prescription performance could be improved by integrating additional tumour information layers like intra-tumour heterogeneity (ITH) which has been related to drug response and tumour progression. PanDrugs is an in silico drug prescription method which prioritizes anticancer drugs combining both biological and clinical evidence. We have systematically evaluated PanDrugs in the Genomic Data Commons repository (GDC). Our results showed that PanDrugs is able to establish an a priori stratification of cancer patients treated with Epidermal Growth Factor Receptor (EGFR) inhibitors. Patients labelled as responders according to PanDrugs predictions showed a significantly increased overall survival (OS) compared to non-responders. PanDrugs was also able to suggest alternative tailored treatments for non-responder patients. Additionally, PanDrugs usefulness was assessed considering spatial and temporal ITH in cancer patients and showed that ITH can be approached therapeutically proposing drugs or combinations potentially capable of targeting the clonal diversity. In summary, this study is a proof of concept where PanDrugs predictions have been correlated to OS and can be useful to manage ITH in patients while increasing therapeutic options and demonstrating its clinical utility.

Cyclooxygenase 2 (COX-2) is a key enzyme of the tumorigenesis-inflammation interface and can be induced by hypoxia. A pseudohypoxic transcriptional signature characterizes pheochromocytomas and paragangliomas (PPGLs) of the cluster I, mainly represented by tumors with mutations in von Hippel-Lindau (VHL), endothelial PAS domain-containing protein 1 (EPAS1), or succinate dehydrogenase (SDH) subunit genes. The aim of this study was to investigate a possible association between underlying tumor driver mutations and COX-2 in PPGLs. COX-2 gene expression and immunoreactivity were examined in clinical specimens with documented mutations, as well as in spheroids and allografts derived from mouse pheochromocytoma (MPC) cells. COX-2 in vivo imaging was performed in allograft mice. We observed significantly higher COX-2 expression in cluster I, especially in VHL-mutant PPGLs, however, no specific association between COX-2 mRNA levels and a hypoxia-related transcriptional signature was found. COX-2 immunoreactivity was present in about 60% of clinical specimens as well as in MPC spheroids and allografts. A selective COX-2 tracer specifically accumulated in MPC allografts. This study demonstrates that, although pseudohypoxia is not the major determinant for high COX-2 levels in PPGLs, COX-2 is a relevant molecular target. This potentially allows for employing selective COX-2 inhibitors as targeted chemotherapeutic agents and radiosensitizers. Moreover, available models are suitable for preclinical testing of these treatments.

Over the past few years, next generation technologies have been applied to unravel the genetics of rare inherited diseases, facilitating the discovery of new susceptibility genes. We recently found germline DNMT3A gain-of-function variants in two patients with head and neck paragangliomas causing a characteristic hypermethylated DNA profile. Here, whole-exome sequencing identifies a novel germline DNMT3A variant (p.Gly332Arg) in a patient with bilateral carotid paragangliomas, papillary thyroid carcinoma and idiopathic intellectual disability. The variant, located in the Pro-Trp-Trp-Pro (PWWP) domain of the protein involved in chromatin targeting, affects a residue mutated in papillary thyroid tumors and located between the two residues found mutated in microcephalic dwarfism patients. Structural modelling of the variant in the DNMT3A PWWP domain predicts that the interaction with H3K36me3 will be altered. An increased methylation of DNMT3A target genes, compatible with a gain-of-function effect of the alteration, was observed in saliva DNA from the proband and in one independent acute myeloid leukemia sample carrying the same p.Gly332Arg variant. Although further studies are needed to support a causal role of DNMT3A variants in paraganglioma, the description of a new DNMT3A alteration in a patient with multiple clinical features suggests a heterogeneous phenotypic spectrum related to DNMT3A germline variants.