@article {406, title = {The first complete genomic structure of Butyrivibrio fibrisolvens and its chromid.}, journal = {Microb Genom}, volume = {4}, year = {2018}, month = {2018 Oct}, abstract = {

Butyrivibrio fibrisolvens forms part of the gastrointestinal microbiome of ruminants and other mammals, including humans. Indeed, it is one of the most common bacteria found in the rumen and plays an important role in ruminal fermentation of polysaccharides, yet, to date, there is no closed reference genome published for this species in any ruminant animal. We successfully assembled the nearly complete genome sequence of B. fibrisolvens strain INBov1 isolated from cow rumen using Illumina paired-end reads, 454 Roche single-end and mate pair sequencing technology. Additionally, we constructed an optical restriction map of this strain to aid in scaffold ordering and positioning, and completed the first genomic structure of this species. Moreover, we identified and assembled the first chromid of this species (pINBov266). The INBov1 genome encodes a large set of genes involved in the cellulolytic process but lacks key genes. This seems to indicate that B. fibrisolvens plays an important role in ruminal cellulolytic processes, but does not have autonomous cellulolytic capacity. When searching for genes involved in the biohydrogenation of unsaturated fatty acids, no linoleate isomerase gene was found in this strain. INBov1 does encode oleate hydratase genes known to participate in the hydrogenation of oleic acids. Furthermore, INBov1 contains an enolase gene, which has been recently determined to participate in the synthesis of conjugated linoleic acids. This work confirms the presence of a novel chromid in B. fibrisolvens and provides a new potential reference genome sequence for this species, providing new insight into its role in biohydrogenation and carbohydrate degradation.

}, keywords = {Animals, Butyrivibrio fibrisolvens, Cattle, Genome, Bacterial, Genomics, Humans, Milk, Rumen, Sequence Analysis, DNA}, issn = {2057-5858}, doi = {10.1099/mgen.0.000216}, author = {Rodr{\'\i}guez Hern{\'a}ez, Javier and Cer{\'o}n Cucchi, Maria Esperanza and Cravero, Silvio and Martinez, Maria Carolina and Gonzalez, Sergio and Puebla, Andrea and Dopazo, Joaquin and Farber, Marisa and Paniego, Norma and Rivarola, M{\'a}ximo} } @article {472, title = {Fast inexact mapping using advanced tree exploration on backward search methods.}, journal = {BMC Bioinformatics}, volume = {16}, year = {2015}, month = {2015 Jan 28}, pages = {18}, abstract = {

BACKGROUND: Short sequence mapping methods for Next Generation Sequencing consist on a combination of seeding techniques followed by local alignment based on dynamic programming approaches. Most seeding algorithms are based on backward search alignment, using the Burrows Wheeler Transform, the Ferragina and Manzini Index or Suffix Arrays. All these backward search algorithms have excellent performance, but their computational cost highly increases when allowing errors. In this paper, we discuss an inexact mapping algorithm based on pruning strategies for search tree exploration over genomic data.

RESULTS: The proposed algorithm achieves a 13x speed-up over similar algorithms when allowing 6 base errors, including insertions, deletions and mismatches. This algorithm can deal with 400 bps reads with up to 9 errors in a high quality Illumina dataset. In this example, the algorithm works as a preprocessor that reduces by 55\% the number of reads to be aligned. Depending on the aligner the overall execution time is reduced between 20-40\%.

CONCLUSIONS: Although not intended as a complete sequence mapping tool, the proposed algorithm could be used as a preprocessing step to modern sequence mappers. This step significantly reduces the number reads to be aligned, accelerating overall alignment time. Furthermore, this algorithm could be used for accelerating the seeding step of already available sequence mappers. In addition, an out-of-core index has been implemented for working with large genomes on systems without expensive memory configurations.

}, keywords = {Algorithms, Genome, Human, Genomics, High-Throughput Nucleotide Sequencing, Humans, Sequence Alignment, Sequence Analysis, DNA, Software}, issn = {1471-2105}, doi = {10.1186/s12859-014-0438-3}, author = {Salavert, Jos{\'e} and Tom{\'a}s, Andr{\'e}s and T{\'a}rraga, Joaqu{\'\i}n and Medina, Ignacio and Dopazo, Joaquin and Blanquer, Ignacio} } @article {467, title = {Re-evaluation casts doubt on the pathogenicity of homozygous USH2A p.C759F.}, journal = {Am J Med Genet A}, volume = {167}, year = {2015}, month = {2015 Jul}, pages = {1597-600}, abstract = {

Mutations in USH2A are a common cause of Retinitis Pigmentosa (RP). Among the most frequently reported USH2A variants, c.2276G>T (p.C759F) has been found in both affected and healthy individuals. The pathogenicity of this variant remains controversial since it was detected in homozygosity in two healthy siblings of a Spanish family (S23), eleven years ago. The fact that these individuals remain asymptomatic today, prompted us to study the presence of other pathogenic variants in this family using targeted resequencing of 26 retinal genes in one of the affected individuals. This approach allowed us to identify one novel pathogenic homozygous mutation in exon 13 of PDE6B (c.1678C>T; p.R560C). This variant cosegregated with the disease and was absent in 200 control individuals. Remarkably, the identified variant in PDE6B corresponds to the mutation responsible of the retinal degeneration in the naturally occurring rd10 mutant mice. To our knowledge, this is the first report of the identification of the rd10 mice mutation in a RP family. These findings, together with a review of the literature, support the hypothesis that homozygous p.C759F mutations are not pathogenic and led us to exclude the implication of p.C759F in the RP of family S23. Our results indicate the need of re-evaluating all families genetically diagnosed with this mutation.

}, keywords = {Base Sequence, Cyclic Nucleotide Phosphodiesterases, Type 6, Extracellular Matrix Proteins, Gene Library, Humans, Molecular Sequence Data, Mutation, Missense, Pedigree, Retinitis pigmentosa, Sequence Analysis, DNA, Spain}, issn = {1552-4833}, doi = {10.1002/ajmg.a.37003}, author = {Pozo, Mar{\'\i}a Gonz{\'a}lez-Del and Bravo-Gil, Nereida and M{\'e}ndez-Vidal, Cristina and Montero-de-Espinosa, Ignacio and Mill{\'a}n, Jos{\'e} M and Dopazo, Joaquin and Borrego, Salud and Anti{\v n}olo, Guillermo} } @article {482, title = {ngsCAT: a tool to assess the efficiency of targeted enrichment sequencing.}, journal = {Bioinformatics}, volume = {30}, year = {2014}, month = {2014 Jun 15}, pages = {1767-8}, abstract = {

MOTIVATION: Targeted enrichment sequencing by next-generation sequencing is a common approach to interrogate specific loci or the whole exome in the human genome. The efficiency and the lack of bias in the enrichment process need to be assessed as a quality control step before performing downstream analysis of the sequence data. Tools that can report on the sensitivity, specificity, uniformity and other enrichment-specific features are needed.

RESULTS: We have implemented the next-generation sequencing data Capture Assessment Tool (ngsCAT), a tool that takes the information of the mapped reads and the coordinates of the targeted regions as input files, and generates a report with metrics and figures that allows the evaluation of the efficiency of the enrichment process. The tool can also take as input the information of two samples allowing the comparison of two different experiments.

AVAILABILITY AND IMPLEMENTATION: Documentation and downloads for ngsCAT can be found at http://www.bioinfomgp.org/ngscat.

}, keywords = {Exome, Genome, Human, High-Throughput Nucleotide Sequencing, Humans, Sequence Analysis, DNA, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btu108}, author = {L{\'o}pez-Domingo, Francisco J and Florido, Javier P and Rueda, Antonio and Dopazo, Joaquin and Santoyo-L{\'o}pez, Javier} } @article {494, title = {The role of the interactome in the maintenance of deleterious variability in human populations.}, journal = {Mol Syst Biol}, volume = {10}, year = {2014}, month = {2014 Sep 26}, pages = {752}, abstract = {

Recent genomic projects have revealed the existence of an unexpectedly large amount of deleterious variability in the human genome. Several hypotheses have been proposed to explain such an apparently high mutational load. However, the mechanisms by which deleterious mutations in some genes cause a pathological effect but are apparently innocuous in other genes remain largely unknown. This study searched for deleterious variants in the 1,000 genomes populations, as well as in a newly sequenced population of 252 healthy Spanish individuals. In addition, variants causative of monogenic diseases and somatic variants from 41 chronic lymphocytic leukaemia patients were analysed. The deleterious variants found were analysed in the context of the interactome to understand the role of network topology in the maintenance of the observed mutational load. Our results suggest that one of the mechanisms whereby the effect of these deleterious variants on the phenotype is suppressed could be related to the configuration of the protein interaction network. Most of the deleterious variants observed in healthy individuals are concentrated in peripheral regions of the interactome, in combinations that preserve their connectivity, and have a marginal effect on interactome integrity. On the contrary, likely pathogenic cancer somatic deleterious variants tend to occur in internal regions of the interactome, often with associated structural consequences. Finally, variants causative of monogenic diseases seem to occupy an intermediate position. Our observations suggest that the real pathological potential of a variant might be more a systems property rather than an intrinsic property of individual proteins.

}, keywords = {Alleles, Exome, Gene Library, Genetic Variation, Genetics, Population, Genome, Human, Genomics, Humans, Models, Genetic, mutation, Phenotype, Protein Conformation, Protein Interaction Maps, Sequence Analysis, DNA, Whites}, issn = {1744-4292}, doi = {10.15252/msb.20145222}, author = {Garc{\'\i}a-Alonso, Luz and Jim{\'e}nez-Almaz{\'a}n, Jorge and Carbonell-Caballero, Jos{\'e} and Vela-Boza, Alicia and Santoyo-L{\'o}pez, Javier and Anti{\v n}olo, Guillermo and Dopazo, Joaquin} } @article {483, title = {Sequencing and functional analysis of the genome of a nematode egg-parasitic fungus, Pochonia chlamydosporia.}, journal = {Fungal Genet Biol}, volume = {65}, year = {2014}, month = {2014 Apr}, pages = {69-80}, abstract = {

Pochonia chlamydosporia is a worldwide-distributed soil fungus with a great capacity to infect and destroy the eggs and kill females of plant-parasitic nematodes. Additionally, it has the ability to colonize endophytically roots of economically-important crop plants, thereby promoting their growth and eliciting plant defenses. This multitrophic behavior makes P. chlamydosporia a potentially useful tool for sustainable agriculture approaches. We sequenced and assembled \~{}41 Mb of P. chlamydosporia genomic DNA and predicted 12,122 gene models, of which many were homologous to genes of fungal pathogens of invertebrates and fungal plant pathogens. Predicted genes (65\%) were functionally annotated according to Gene Ontology, and 16\% of them found to share homology with genes in the Pathogen Host Interactions (PHI) database. The genome of this fungus is highly enriched in genes encoding hydrolytic enzymes, such as proteases, glycoside hydrolases and carbohydrate esterases. We used RNA-Seq technology in order to identify the genes expressed during endophytic behavior of P. chlamydosporia when colonizing barley roots. Functional annotation of these genes showed that hydrolytic enzymes and transporters are expressed during endophytism. This structural and functional analysis of the P. chlamydosporia genome provides a starting point for understanding the molecular mechanisms involved in the multitrophic lifestyle of this fungus. The genomic information provided here should also prove useful for enhancing the capabilities of this fungus as a biocontrol agent of plant-parasitic nematodes and as a plant growth-promoting organism.

}, keywords = {Animals, Ascomycota, Female, Gene Expression Regulation, Fungal, Gene ontology, Genome, Fungal, Hordeum, Host-Pathogen Interactions, Nematoda, Ovum, Phylogeny, Plant Roots, Sequence Analysis, DNA, Signal Transduction, Transcriptome}, issn = {1096-0937}, doi = {10.1016/j.fgb.2014.02.002}, author = {Larriba, Eduardo and Jaime, Mar{\'\i}a D L A and Carbonell-Caballero, Jos{\'e} and Conesa, Ana and Dopazo, Joaquin and Nislow, Corey and Mart{\'\i}n-Nieto, Jos{\'e} and Lopez-Llorca, Luis Vicente} } @article {513, title = {Diversification of the expanded teleost-specific toll-like receptor family in Atlantic cod, Gadus morhua.}, journal = {BMC Evol Biol}, volume = {12}, year = {2012}, month = {2012 Dec 29}, pages = {256}, abstract = {

BACKGROUND: Toll-like receptors (Tlrs) are major molecular pattern recognition receptors of the innate immune system. Atlantic cod (Gadus morhua) is the first vertebrate known to have lost most of the mammalian Tlr orthologues, particularly all bacterial recognising and other cell surface Tlrs. On the other hand, its genome encodes a unique repertoire of teleost-specific Tlrs. The aim of this study was to investigate if these duplicate Tlrs have been retained through adaptive evolution to compensate for the lack of other cell surface Tlrs in the cod genome.

RESULTS: In this study, one tlr21, 12 tlr22 and two tlr23 genes representing the teleost-specific Tlr family have been cloned and characterised in cod. Phylogenetic analysis grouped all tlr22 genes under a single clade, indicating that the multiple cod paralogues have arisen through lineage-specific duplications. All tlrs examined were transcribed in immune-related tissues as well as in stomach, gut and gonads of adult cod and were differentially expressed during early development. These tlrs were also differentially regulated following immune challenge by immersion with Vibrio anguillarum, indicating their role in the immune response. An increase in water temperature from 4 to 12{\textdegree}C was associated with a 5.5-fold down-regulation of tlr22d transcript levels in spleen. Maximum likelihood analysis with different evolution models revealed that tlr22 genes are under positive selection. A total of 24 codons were found to be positively selected, of which 19 are in the ligand binding region of ectodomain.

CONCLUSION: Positive selection pressure coupled with experimental evidence of differential expression strongly support the hypothesis that teleost-specific tlr paralogues in cod are undergoing neofunctionalisation and can recognise bacterial pathogen-associated molecular patterns to compensate for the lack of other cell surface Tlrs.

}, keywords = {Amino Acid Sequence, Animals, Binding Sites, Evolution, Molecular, Fish Diseases, Fish Proteins, Gadus morhua, Gene Expression Profiling, Genetic Variation, Gills, Head Kidney, Host-Pathogen Interactions, Models, Molecular, Molecular Sequence Data, Multigene Family, Phylogeny, Protein Structure, Tertiary, Reverse Transcriptase Polymerase Chain Reaction, Selection, Genetic, Sequence Analysis, DNA, Sequence Homology, Amino Acid, Temperature, Toll-Like Receptors, Vibrio}, issn = {1471-2148}, doi = {10.1186/1471-2148-12-256}, author = {Sundaram, Arvind Y M and Kiron, Viswanath and Dopazo, Joaquin and Fernandes, Jorge M O} } @article {524, title = {Whole-genome bisulfite DNA sequencing of a DNMT3B mutant patient.}, journal = {Epigenetics}, volume = {7}, year = {2012}, month = {2012 Jun 01}, pages = {542-50}, abstract = {

The immunodeficiency, centromere instability and facial anomalies (ICF) syndrome is associated to mutations of the DNA methyl-transferase DNMT3B, resulting in a reduction of enzyme activity. Aberrant expression of immune system genes and hypomethylation of pericentromeric regions accompanied by chromosomal instability were determined as alterations driving the disease phenotype. However, so far only technologies capable to analyze single loci were applied to determine epigenetic alterations in ICF patients. In the current study, we performed whole-genome bisulphite sequencing to assess alteration in DNA methylation at base pair resolution. Genome-wide we detected a decrease of methylation level of 42\%, with the most profound changes occurring in inactive heterochromatic regions, satellite repeats and transposons. Interestingly, transcriptional active loci and ribosomal RNA repeats escaped global hypomethylation. Despite a genome-wide loss of DNA methylation the epigenetic landscape and crucial regulatory structures were conserved. Remarkably, we revealed a mislocated activity of mutant DNMT3B to H3K4me1 loci resulting in hypermethylation of active promoters. Functionally, we could associate alterations in promoter methylation with the ICF syndrome immunodeficient phenotype by detecting changes in genes related to the B-cell receptor mediated maturation pathway.

}, keywords = {B-Lymphocytes, Cell Line, Transformed, Child, Preschool, DNA (Cytosine-5-)-Methyltransferases, DNA Methylation, Epigenesis, Genetic, Face, Female, Genome, Human, High-Throughput Nucleotide Sequencing, Humans, Immunologic Deficiency Syndromes, mutation, Primary Immunodeficiency Diseases, Sequence Analysis, DNA, Sulfites}, issn = {1559-2308}, doi = {10.4161/epi.20523}, author = {Heyn, Holger and Vidal, Enrique and Sayols, Sergi and Sanchez-Mut, Jose V and Moran, Sebastian and Medina, Ignacio and Sandoval, Juan and Sim{\'o}-Riudalbas, Laia and Szczesna, Karolina and Huertas, Dori and Gatto, Sole and Matarazzo, Maria R and Dopazo, Joaquin and Esteller, Manel} } @article {537, title = {Natural selection on functional modules, a genome-wide analysis.}, journal = {PLoS Comput Biol}, volume = {7}, year = {2011}, month = {2011 Mar}, pages = {e1001093}, abstract = {

Classically, the functional consequences of natural selection over genomes have been analyzed as the compound effects of individual genes. The current paradigm for large-scale analysis of adaptation is based on the observed significant deviations of rates of individual genes from neutral evolutionary expectation. This approach, which assumed independence among genes, has not been able to identify biological functions significantly enriched in positively selected genes in individual species. Alternatively, pooling related species has enhanced the search for signatures of selection. However, grouping signatures does not allow testing for adaptive differences between species. Here we introduce the Gene-Set Selection Analysis (GSSA), a new genome-wide approach to test for evidences of natural selection on functional modules. GSSA is able to detect lineage specific evolutionary rate changes in a notable number of functional modules. For example, in nine mammal and Drosophilae genomes GSSA identifies hundreds of functional modules with significant associations to high and low rates of evolution. Many of the detected functional modules with high evolutionary rates have been previously identified as biological functions under positive selection. Notably, GSSA identifies conserved functional modules with many positively selected genes, which questions whether they are exclusively selected for fitting genomes to environmental changes. Our results agree with previous studies suggesting that adaptation requires positive selection, but not every mutation under positive selection contributes to the adaptive dynamical process of the evolution of species.

}, keywords = {Animals, Databases, Genetic, Drosophila, Genome, Insect, Genome-Wide Association Study, Genomics, Mammals, Phylogeny, Selection, Genetic, Sequence Analysis, DNA}, issn = {1553-7358}, doi = {10.1371/journal.pcbi.1001093}, author = {Serra, Fran{\c c}ois and Arbiza, Leonardo and Dopazo, Joaquin and Dopazo, Hern{\'a}n} } @article {594, title = {High-throughput functional annotation and data mining with the Blast2GO suite.}, journal = {Nucleic Acids Res}, volume = {36}, year = {2008}, month = {2008 Jun}, pages = {3420-35}, abstract = {

Functional genomics technologies have been widely adopted in the biological research of both model and non-model species. An efficient functional annotation of DNA or protein sequences is a major requirement for the successful application of these approaches as functional information on gene products is often the key to the interpretation of experimental results. Therefore, there is an increasing need for bioinformatics resources which are able to cope with large amount of sequence data, produce valuable annotation results and are easily accessible to laboratories where functional genomics projects are being undertaken. We present the Blast2GO suite as an integrated and biologist-oriented solution for the high-throughput and automatic functional annotation of DNA or protein sequences based on the Gene Ontology vocabulary. The most outstanding Blast2GO features are: (i) the combination of various annotation strategies and tools controlling type and intensity of annotation, (ii) the numerous graphical features such as the interactive GO-graph visualization for gene-set function profiling or descriptive charts, (iii) the general sequence management features and (iv) high-throughput capabilities. We used the Blast2GO framework to carry out a detailed analysis of annotation behaviour through homology transfer and its impact in functional genomics research. Our aim is to offer biologists useful information to take into account when addressing the task of functionally characterizing their sequence data.

}, keywords = {Animals, Computational Biology, Computer Graphics, Databases, Genetic, Expressed Sequence Tags, Genes, Genomics, Sequence Analysis, DNA, Sequence Analysis, Protein, Software, Vocabulary, Controlled}, issn = {1362-4962}, doi = {10.1093/nar/gkn176}, author = {G{\"o}tz, Stefan and Garc{\'\i}a-G{\'o}mez, Juan Miguel and Terol, Javier and Williams, Tim D and Nagaraj, Shivashankar H and Nueda, Maria Jos{\'e} and Robles, Montserrat and Talon, Manuel and Dopazo, Joaquin and Conesa, Ana} }