@article {705, title = {SMN1 copy-number and sequence variant analysis from next-generation sequencing data.}, journal = {Hum Mutat}, volume = {41}, year = {2020}, month = {2020 12}, pages = {2073-2077}, abstract = {

Spinal muscular atrophy (SMA) is a severe neuromuscular autosomal recessive disorder affecting 1/10,000 live births. Most SMA patients present homozygous deletion of SMN1, while the vast majority of SMA carriers present only a single SMN1 copy. The sequence similarity between SMN1 and SMN2, and the complexity of the SMN locus makes the estimation of the SMN1 copy-number by next-generation sequencing (NGS) very difficult. Here, we present SMAca, the first python tool to detect SMA carriers and estimate the absolute SMN1 copy-number using NGS data. Moreover, SMAca takes advantage of the knowledge of certain variants specific to SMN1 duplication to also identify silent carriers. This tool has been validated with a cohort of 326 samples from the Navarra 1000 Genomes Project (NAGEN1000). SMAca was developed with a focus on execution speed and easy installation. This combination makes it especially suitable to be integrated into production NGS pipelines. Source code and documentation are available at https://www.github.com/babelomics/SMAca.

}, keywords = {Base Sequence, DNA Copy Number Variations, High-Throughput Nucleotide Sequencing, Humans, Reproducibility of Results, Software, Survival of Motor Neuron 1 Protein}, issn = {1098-1004}, doi = {10.1002/humu.24120}, author = {L{\'o}pez-L{\'o}pez, Daniel and Loucera, Carlos and Carmona, Rosario and Aquino, Virginia and Salgado, Josefa and Pasalodos, Sara and Miranda, Mar{\'\i}a and Alonso, {\'A}ngel and Dopazo, Joaquin} } @article {388, title = {Reference genome assessment from a population scale perspective: an accurate profile of variability and noise.}, journal = {Bioinformatics}, volume = {33}, year = {2017}, month = {2017 Nov 15}, pages = {3511-3517}, abstract = {

MOTIVATION: Current plant and animal genomic studies are often based on newly assembled genomes that have not been properly consolidated. In this scenario, misassembled regions can easily lead to false-positive findings. Despite quality control scores are included within genotyping protocols, they are usually employed to evaluate individual sample quality rather than reference sequence reliability. We propose a statistical model that combines quality control scores across samples in order to detect incongruent patterns at every genomic region. Our model is inherently robust since common artifact signals are expected to be shared between independent samples over misassembled regions of the genome.

RESULTS: The reliability of our protocol has been extensively tested through different experiments and organisms with accurate results, improving state-of-the-art methods. Our analysis demonstrates synergistic relations between quality control scores and allelic variability estimators, that improve the detection of misassembled regions, and is able to find strong artifact signals even within the human reference assembly. Furthermore, we demonstrated how our model can be trained to properly rank the confidence of a set of candidate variants obtained from new independent samples.

AVAILABILITY AND IMPLEMENTATION: This tool is freely available at http://gitlab.com/carbonell/ces.

CONTACT: jcarbonell.cipf@gmail.com or joaquin.dopazo@juntadeandalucia.es.

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

}, keywords = {Animals, Genetic Variation, Genome, Genomics, Genotype, Humans, Models, Statistical, Quality Control, Reproducibility of Results, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btx482}, url = {https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btx482}, author = {Carbonell-Caballero, Jos{\'e} and Amadoz, Alicia and Alonso, Roberto and Hidalgo, Marta R and Cubuk, Cankut and Conesa, David and L{\'o}pez-Qu{\'\i}lez, Antonio and Dopazo, Joaquin} } @article {439, title = {Integrated gene set analysis for microRNA studies.}, journal = {Bioinformatics}, volume = {32}, year = {2016}, month = {2016 Sep 15}, pages = {2809-16}, abstract = {

MOTIVATION: Functional interpretation of miRNA expression data is currently done in a three step procedure: select differentially expressed miRNAs, find their target genes, and carry out gene set overrepresentation analysis Nevertheless, major limitations of this approach have already been described at the gene level, while some newer arise in the miRNA scenario.Here, we propose an enhanced methodology that builds on the well-established gene set analysis paradigm. Evidence for differential expression at the miRNA level is transferred to a gene differential inhibition score which is easily interpretable in terms of gene sets or pathways. Such transferred indexes account for the additive effect of several miRNAs targeting the same gene, and also incorporate cancellation effects between cases and controls. Together, these two desirable characteristics allow for more accurate modeling of regulatory processes.

RESULTS: We analyze high-throughput sequencing data from 20 different cancer types and provide exhaustive reports of gene and Gene Ontology-term deregulation by miRNA action.

AVAILABILITY AND IMPLEMENTATION: The proposed methodology was implemented in the Bioconductor library mdgsa http://bioconductor.org/packages/mdgsa For the purpose of reproducibility all of the scripts are available at https://github.com/dmontaner-papers/gsa4mirna

CONTACT: : david.montaner@gmail.com

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

}, keywords = {Computational Biology, Gene Expression Profiling, Gene ontology, Gene Regulatory Networks, High-Throughput Nucleotide Sequencing, Humans, MicroRNAs, Neoplasms, Reproducibility of Results}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btw334}, author = {Garcia-Garcia, Francisco and Panadero, Joaquin and Dopazo, Joaquin and Montaner, David} } @article {453, title = {The Mutational Landscape of Acute Promyelocytic Leukemia Reveals an Interacting Network of Co-Occurrences and Recurrent Mutations.}, journal = {PLoS One}, volume = {11}, year = {2016}, month = {2016}, pages = {e0148346}, abstract = {

Preliminary Acute Promyelocytic Leukemia (APL) whole exome sequencing (WES) studies have identified a huge number of somatic mutations affecting more than a hundred different genes mainly in a non-recurrent manner, suggesting that APL is a heterogeneous disease with secondary relevant changes not yet defined. To extend our knowledge of subtle genetic alterations involved in APL that might cooperate with PML/RARA in the leukemogenic process, we performed a comprehensive analysis of somatic mutations in APL combining WES with sequencing of a custom panel of targeted genes by next-generation sequencing. To select a reduced subset of high confidence candidate driver genes, further in silico analysis were carried out. After prioritization and network analysis we found recurrent deleterious mutations in 8 individual genes (STAG2, U2AF1, SMC1A, USP9X, IKZF1, LYN, MYCBP2 and PTPN11) with a strong potential of being involved in APL pathogenesis. Our network analysis of multiple mutations provides a reliable approach to prioritize genes for additional analysis, improving our knowledge of the leukemogenesis interactome. Additionally, we have defined a functional module in the interactome of APL. The hypothesis is that the number, or the specific combinations, of mutations harbored in each patient might not be as important as the disturbance caused in biological key functions, triggered by several not necessarily recurrent mutations.

}, keywords = {Exome, Gene Regulatory Networks, Genome, Human, Humans, INDEL Mutation, Leukemia, Promyelocytic, Acute, mutation, Mutation Rate, Polymorphism, Single Nucleotide, Reproducibility of Results}, issn = {1932-6203}, doi = {10.1371/journal.pone.0148346}, author = {Ib{\'a}{\~n}ez, Mariam and Carbonell-Caballero, Jos{\'e} and Garc{\'\i}a-Alonso, Luz and Such, Esperanza and Jim{\'e}nez-Almaz{\'a}n, Jorge and Vidal, Enrique and Barrag{\'a}n, Eva and L{\'o}pez-Pav{\'\i}a, Mar{\'\i}a and LLop, Marta and Mart{\'\i}n, Iv{\'a}n and G{\'o}mez-Segu{\'\i}, In{\'e}s and Montesinos, Pau and Sanz, Miguel A and Dopazo, Joaquin and Cervera, Jos{\'e}} } @article {532, title = {Early transcriptional defense responses in Arabidopsis cell suspension culture under high-light conditions.}, journal = {Plant Physiol}, volume = {156}, year = {2011}, month = {2011 Jul}, pages = {1439-56}, abstract = {

The early transcriptional defense responses and reactive oxygen species (ROS) production in Arabidopsis (Arabidopsis thaliana) cell suspension culture (ACSC), containing functional chloroplasts, were examined at high light (HL). The transcriptional analysis revealed that most of the ROS markers identified among the 449 transcripts with significant differential expression were transcripts specifically up-regulated by singlet oxygen ((1)O(2)). On the contrary, minimal correlation was established with transcripts specifically up-regulated by superoxide radical or hydrogen peroxide. The transcriptional analysis was supported by fluorescence microscopy experiments. The incubation of ACSC with the (1)O(2) sensor green reagent and 2{\textquoteright},7{\textquoteright}-dichlorofluorescein diacetate showed that the 30-min-HL-treated cultures emitted fluorescence that corresponded with the production of (1)O(2) but not of hydrogen peroxide. Furthermore, the in vivo photodamage of the D1 protein of photosystem II indicated that the photogeneration of (1)O(2) took place within the photosystem II reaction center. Functional enrichment analyses identified transcripts that are key components of the ROS signaling transduction pathway in plants as well as others encoding transcription factors that regulate both ROS scavenging and water deficit stress. A meta-analysis examining the transcriptional profiles of mutants and hormone treatments in Arabidopsis showed a high correlation between ACSC at HL and the fluorescent mutant family of Arabidopsis, a producer of (1)O(2) in plastids. Intriguingly, a high correlation was also observed with ABA deficient1 and more axillary growth4, two mutants with defects in the biosynthesis pathways of two key (apo)carotenoid-derived plant hormones (i.e. abscisic acid and strigolactones, respectively). ACSC has proven to be a valuable system for studying early transcriptional responses to HL stress.

}, keywords = {Arabidopsis, Blotting, Western, Cell Culture Techniques, Cells, Cultured, Chloroplasts, Cluster Analysis, Gene Expression Profiling, Gene Expression Regulation, Plant, Hydrogen Peroxide, Light, mutation, Oligonucleotide Array Sequence Analysis, Photosystem II Protein Complex, Plant Growth Regulators, Reproducibility of Results, Reverse Transcriptase Polymerase Chain Reaction, RNA, Messenger, Signal Transduction, Stress, Physiological, Transcription, Genetic}, issn = {1532-2548}, doi = {10.1104/pp.111.177766}, author = {Gonz{\'a}lez-P{\'e}rez, Sergio and Guti{\'e}rrez, Jorge and Garcia-Garcia, Francisco and Osuna, Daniel and Dopazo, Joaquin and Lorenzo, Oscar and Revuelta, Jos{\'e} L and Arellano, Juan B} } @article {536, title = {Mutation screening of multiple genes in Spanish patients with autosomal recessive retinitis pigmentosa by targeted resequencing.}, journal = {PLoS One}, volume = {6}, year = {2011}, month = {2011}, pages = {e27894}, abstract = {

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.

}, keywords = {Alleles, DNA Mutational Analysis, Exons, Genetic Variation, Genome, Hispanic or Latino, Humans, Introns, Language, mutation, Mutation, Missense, Oligonucleotide Array Sequence Analysis, Polymerase Chain Reaction, Reproducibility of Results, Retinitis pigmentosa, United States}, issn = {1932-6203}, doi = {10.1371/journal.pone.0027894}, author = {Gonz{\'a}lez-del Pozo, Mar{\'\i}a and Borrego, Salud and Barrag{\'a}n, Isabel and Pieras, Juan I and Santoyo, Javier and Matamala, Nerea and Naranjo, Bel{\'e}n and Dopazo, Joaquin and Anti{\v n}olo, Guillermo} } @article {582, title = {Gene set internal coherence in the context of functional profiling.}, journal = {BMC Genomics}, volume = {10}, year = {2009}, month = {2009 Apr 27}, pages = {197}, abstract = {

BACKGROUND: Functional profiling methods have been extensively used in the context of high-throughput experiments and, in particular, in microarray data analysis. Such methods use available biological information to define different types of functional gene modules (e.g. gene ontology -GO-, KEGG pathways, etc.) whose representation in a pre-defined list of genes is further studied. In the most popular type of microarray experimental designs (e.g. up- or down-regulated genes, clusters of co-expressing genes, etc.) or in other genomic experiments (e.g. Chip-on-chip, epigenomics, etc.) these lists are composed by genes with a high degree of co-expression. Therefore, an implicit assumption in the application of functional profiling methods within this context is that the genes corresponding to the modules tested are effectively defining sets of co-expressing genes. Nevertheless not all the functional modules are biologically coherent entities in terms of co-expression, which will eventually hinder its detection with conventional methods of functional enrichment.

RESULTS: Using a large collection of microarray data we have carried out a detailed survey of internal correlation in GO terms and KEGG pathways, providing a coherence index to be used for measuring functional module co-regulation. An unexpected low level of internal correlation was found among the modules studied. Only around 30\% of the modules defined by GO terms and 57\% of the modules defined by KEGG pathways display an internal correlation higher than the expected by chance.This information on the internal correlation of the genes within the functional modules can be used in the context of a logistic regression model in a simple way to improve their detection in gene expression experiments.

CONCLUSION: For the first time, an exhaustive study on the internal co-expression of the most popular functional categories has been carried out. Interestingly, the real level of coexpression within many of them is lower than expected (or even inexistent), which will preclude its detection by means of most conventional functional profiling methods. If the gene-to-function correlation information is used in functional profiling methods, the results obtained improve the ones obtained by conventional enrichment methods.

}, keywords = {Algorithms, Breast Neoplasms, Carcinoma, Intraductal, Noninfiltrating, Computational Biology, Databases, Nucleic Acid, Female, Gene Expression Profiling, Genomics, Humans, Oligonucleotide Array Sequence Analysis, Papillomavirus Infections, Reproducibility of Results}, issn = {1471-2164}, doi = {10.1186/1471-2164-10-197}, author = {Montaner, David and Minguez, Pablo and Al-Shahrour, F{\'a}tima and Dopazo, Joaquin} }