BISQUE: locus- and variant-specific conversion of genomic, transcriptomic and proteomic database identifiers.

Acetylcholinesterase in man is encoded by a single gene, ACHE, located on chromosome 7q22. In this study, the transcription start sites and major DNA promoter elements controlling the expression of this gene have been characterized by structural and functional studies. Immediately upstream of the first untranslated exon of the gene are GC-rich sequences containing consensus binding sites for several transcription factors, including Sp1, EGR-1 and AP2. In vitro transcription studies and RNase protection analyses of mRNA isolated from human NT2/D1 teratocarcinoma cells reveal that two closely spaced transcription cap sites are located at a consensus initiator (Inr) element similar to that found in the terminal transferase gene. Transient transfection of mutant genes shows that removal of three bases of this initiator sequence reduces promoter activity by 98% in NT2/D1 cells. In vitro transcription studies and transient transfection of a series of 5' deletion mutants of the ACHE promoter linked to a luciferase reporter show an Sp1 site at -71 to be essential for promoter activity. Purified Sp1 protein protects this site from DNase cleavage during in vitro footprinting experiments. A conserved AP2 consensus binding site, located between the GC box elements and the Inr, is protected by recombinant AP2 protein in DNase footprinting experiments, induces a mobility shift with AP2 protein and AP2-containing cell extracts, and fosters inhibition of transcription by AP2 as measured by transient transfection in mouse and human cell lines and in in vitro transcription reactions. These results indicate that AP2 functions as a repressor of human ACHE and mouse Ache transcription.

Acyl-coenzyme A synthetases (ACSs) catalyze the fundamental, initial reaction in fatty acid metabolism. "Activation" of fatty acids by thioesterification to CoA allows their participation in both anabolic and catabolic pathways. The availability of the sequenced human genome has facilitated the investigation of the number of ACS genes present. Using two conserved amino acid sequence motifs to probe human DNA databases, 26 ACS family genes/proteins were identified. ACS activity in either humans or rodents was demonstrated previously for 20 proteins, but 6 remain candidate ACSs. For two candidates, cDNA was cloned, protein was expressed in COS-1 cells, and ACS activity was detected. Amino acid sequence similarities were used to assign enzymes into subfamilies, and subfamily assignments were consistent with acyl chain length preference. Four of the 26 proteins did not fit into a subfamily, and bootstrap analysis of phylograms was consistent with evolutionary divergence. Three additional conserved amino acid sequence motifs were identified that likely have functional or structural roles. The existence of many ACSs suggests that each plays a unique role, directing the acyl-CoA product to a specific metabolic fate. Knowing the full complement of ACS genes in the human genome will facilitate future studies to characterize their specific biological functions.

Duplication in biological sequence databases has persisted for 20 years. Duplicate records introduce redundancies to databases, delay biocuration processes, and undermine the accuracy of studies based on sequence analysis such as GC content and melting temperature. Rapid growth of data makes purely manual de-duplication nearly impossible, and existing automatic systems cannot detect duplicates as precisely as experts. Supervised learning has the potential to address such problems by building automatic systems that learn from expert curation to detect duplicates precisely and efficiently. While a mature approach in other duplicate detection contexts, machine learning has seen only preliminary application in the large biological sequence databases. We developed a supervised duplicate detection method, employing an over one million-pair expert curated dataset of duplicates across five organisms derived from genomic sequence databases. We selected 22 features to represent distinct attributes of the database records, and developed both binary and multi-class models. Both models achieve promising performance; the binary model had over 90% accuracy in all the 5 organisms while the multi-class model maintains high accuracy and is more robust in generalisation. We performed an ablation study to quantify the impact of different sequence record features, finding that features derived from meta-data, sequence identity, and alignment quality impact performance most strongly. In particular, better measurement on sequences drives the performance.

The evolutionary history of a protein reflects the functional history of its ancestors. Recent phylogenetic studies identified distinct evolutionary signatures that characterize proteins involved in cancer, Mendelian disease, and different ontogenic stages. Despite the potential to yield insight into the cellular functions and interactions of proteins, such comparative phylogenetic analyses are rarely performed, because they require custom algorithms. We developed ProteinHistorian to make tools for performing analyses of protein origins widely available. Given a list of proteins of interest, ProteinHistorian estimates the phylogenetic age of each protein, quantifies enrichment for proteins of specific ages, and compares variation in protein age with other protein attributes. ProteinHistorian allows flexibility in the definition of protein age by including several algorithms for estimating ages from different databases of evolutionary relationships. We illustrate the use of ProteinHistorian with three example analyses. First, we demonstrate that proteins with high expression in human, compared to chimpanzee and rhesus macaque, are significantly younger than those with human-specific low expression. Next, we show that human proteins with annotated regulatory functions are significantly younger than proteins with catalytic functions. Finally, we compare protein length and age in many eukaryotic species and, as expected from previous studies, find a positive, though often weak, correlation between protein age and length. ProteinHistorian is available through a web server with an intuitive interface and as a set of command line tools; this allows biologists and bioinformaticians alike to integrate these approaches into their analysis pipelines. ProteinHistorian's modular, extensible design facilitates the integration of new datasets and algorithms. The ProteinHistorian web server, source code, and pre-computed ages for 32 eukaryotic genomes are freely available under the GNU public license at http://lighthouse.ucsf.edu/ProteinHistorian/.

BackgroundAlternative splicing can increase the diversity of gene functions by generating multiple isoforms with different sequences and functions. However, the extent to which splicing events have functional consequences remains unclear and predicting the impact of splicing events on protein activity is limited to gene-specific analysis.ResultsTo accelerate the identification of functionally relevant alternative splicing events we created SAPFIR, a predictor of protein features associated with alternative splicing events. This webserver tool uses InterProScan to predict protein features such as functional domains, motifs and sites in the human and mouse genomes and link them to alternative splicing events. Alternative protein features are displayed as functions of the transcripts and splice sites. SAPFIR could be used to analyze proteins generated from a single gene or a group of genes and can directly identify alternative protein features in large sequence data sets. The accuracy and utility of SAPFIR was validated by its ability to rediscover previously validated alternative protein domains. In addition, our de novo analysis of public datasets using SAPFIR indicated that only a small portion of alternative protein domains was conserved between human and mouse, and that in human, genes involved in nervous system process, regulation of DNA-templated transcription and aging are more likely to produce isoforms missing functional domains due to alternative splicing.ConclusionOverall SAPFIR represents a new tool for the rapid identification of functional alternative splicing events and enables the identification of cellular functions affected by a defined splicing program. SAPFIR is freely available at https://bioinfo-scottgroup.med.usherbrooke.ca/sapfir/, a website implemented in Python, with all major browsers supported. The source code is available at https://github.com/DelongZHOU/SAPFIR.

BackgroundA variety of high-throughput analyses, such as transcriptome, proteome, and metabolome analysis, have been developed, producing unprecedented amounts of omics data. These studies generate large gene lists, of which the biological significance shall be deeply understood. However, manually interpreting these lists is difficult, especially for non-bioinformatics-savvy scientists.ResultsWe developed an R package and a corresponding web server—Genekitr, to assist biologists in exploring large gene sets. Genekitr comprises four modules: gene information retrieval, ID (identifier) conversion, enrichment analysis and publication-ready plotting. Currently, the information retrieval module can retrieve information on up to 23 attributes for genes of 317 organisms. The ID conversion module assists in ID-mapping of genes, probes, proteins, and aliases. The enrichment analysis module organizes 315 gene set libraries in different biological contexts by over-representation analysis and gene set enrichment analysis. The plotting module performs customizable and high-quality illustrations that can be used directly in presentations or publications.ConclusionsThis web server tool will make bioinformatics more accessible to scientists who might not have programming expertise, allowing them to perform bioinformatics tasks without coding.

Structural genomics provides an important approach for characterizing and understanding systems biology. As a step toward better integrating protein three-dimensional (3D) structural information in cancer systems biology, we have constructed a Human Cancer Pathway Protein Interaction Network (HCPIN) by analysis of several classical cancer-associated signaling pathways and their physical protein-protein interactions. Many well known cancer-associated proteins play central roles as "hubs" or "bottlenecks" in the HCPIN. At least half of HCPIN proteins are either directly associated with or interact with multiple signaling pathways. Although some 45% of residues in these proteins are in sequence segments that meet criteria sufficient for approximate homology modeling (Basic Local Alignment Search Tool (BLAST) E-value <10(-6)), only approximately 20% of residues in these proteins are structurally covered using high accuracy homology modeling criteria (i.e. BLAST E-value <10(-6) and at least 80% sequence identity) or by actual experimental structures. The HCPIN Website provides a comprehensive description of this biomedically important multipathway network together with experimental and homology models of HCPIN proteins useful for cancer biology research. To complement and enrich cancer systems biology, the Northeast Structural Genomics Consortium is targeting >1000 human proteins and protein domains from the HCPIN for sample production and 3D structure determination. The long range goal of this effort is to provide a comprehensive 3D structure-function database for human cancer-associated proteins and protein complexes in the context of their interaction networks. The network-based target selection (BioNet) approach described here is an example of a general strategy for targeting co-functioning proteins by structural genomics projects.

Over the last decade, Haplogrep has become a standard tool for haplogroup classification in the field of human mitochondrial DNA and is widely used by medical, forensic, and evolutionary researchers. Haplogrep scales well for thousands of samples, supports many file formats and provides an intuitive graphical web interface. Nevertheless, the currently available version has limitations when applying it to large biobank-scale data. In this paper, we present a major upgrade to the software by adding (a) haplogroup summary statistics and variant annotations from various publicly available genome databases, (b) an interface to connect new phylogenetic trees, (c) a new state-of-the-art web framework managing large scale data, (d) algorithmic adaptions to improve FASTA classification using BWA-specific alignment rules and (e) a pre-classification quality control step for VCF samples. These improvements will give researchers the opportunity to classify thousands of samples as usual but providing additional ways to investigate the dataset directly in the browser. The web service and its documentation can be accessed freely without any registration at https://haplogrep.i-med.ac.at.

Chromosome 21 is the smallest human chromosome and represents a model for physical mapping of the entire human genome. Three copies of this autosome cause Down syndrome, the most frequent genetic disorder associated with significant mental retardation. Five years ago, in the framework of the international human genome project, and as an integrated part of the chromosome 21 community, a consortium, of academic groups from Japan and Germany was formed to map and sequence this chromosome. Several whole-genome and chromosome-specific DNA libraries were constructed for mapping purposes. Using a combination of nested-deletion [1] and shotgun sequencing approaches, our center determined and analyzed over 17 million bases of high-quality data. In total, over 33.5 million basepairs of DNA, distributed over four contigs, were sequenced with very high accuracy [2]. The largest of these contigs is nearly 25.5Mb long. Only three small clone gaps remain, which together comprise approximately 100kb. Thus, we achieved a coverage of 99.7% of 21q. In addition, 211,116 bp from the short arm were also sequenced. Analysis of the chromosome revealed far fewer genes than previously estimated. Instead of 800-1000 expected genes, we found only 225 genes, of which 127 are known genes and 98 are predicted genes. In addition, 59 pseudogenes were identified. The completion of the 21q sequence provides a unique resource for understanding the molecular pathophysiology of Down syndrome, as well as all other monogenic and complex disorders that map to this chromosome, including Alzheimer’s disease, leukemia, autoimmune disease, epilepsy and manic-depressive psychosis. It also stands as a structural framework from which the complete molecular architecture of the chromosome can be determined. With the announcement of the first “working draft” of the human genome by the international human genome sequencing consortium this past June, there has been a great interest in the total number of genes in the human genome[3]. Some of the current estimates range from as low as 34,000 genes to as high as 140,000 genes. The low gene-density of chromosome 21 was quite unexpected. Most striking is a 7-Mb region near the region near the centromere that contains only one gene. This region is much larger than the whole genome of some species, such as Escherichia coli, yet evolutionary processes permitted the existence of such a gene-poor SNA segment. This finding leads us to propose that there are additional large gene-less regions in other human/mammalian chromosomes. By combining the gene numbers of chromosomes 21 and 22 [4] (770genes; 2-3% of human genome) and assuming that together they represent the average gene content of the human genome, we estimate that the total number of human genes may be close to 40,000. As part of the international human genome sequencing project, we are also sequencing parts of chromosomes 11 and 18. Our high-throughput sequence production line generates about 10 Mb per day. In order to handle such a high volume of output, we have developed a system for automated data assembly, annotation and release. All the data is released through our website and through the DNA Database of Japan immediately according to the policies established by the international consortium. In addition to chromosomes 11, 18 and 21, we are also interested in the structure of the entire human genome. In collaboration with the University of Tokyo, we have developed a database called HGREP [5] (Human Genome Reconstruction Project) that contains working draft and finished sequences which cover more than 85% of the human genomic sequences. Sequence entries are aligned along the chromosomes based on sequence similarities to STS markers, BAC-end and other entry sequences. Further, biological features, such as genes, gene functions, repeats and CpG islands, are fully annotated on these entries. To identify evolutionally conserved and biologically important information in the genome, we are taking the strategy of comparative genomic sequence analysis between human and other related genomes such as those of rodents and primates. We are currently sequencing several regions of the mouse, including some regions that are counterpart to human chromosome 21, and are preparing genomic resources for sequencing of the chimpanzee. Through comparative methods, we are also trying to understand telomeric regions, structures at the end of chromosomes that play an important role in several biological functions and have been associated with human diseases like an oncogenesis and a cellular aging. The structure and function of all genes and their regulatory regions can only be fully understood by looking at how they interact under different circumstances and in several different backgrounds. Determining which genes and elements are human-specific will lead to a deeper understanding of who we are and will advance medical sciences drastically.

N-Glycolylneuraminic acid (NeuGc) is abundantly expressed in most mammals, but it is not detectable in humans. The expression of NeuGc is controlled by cytidine monophospho-N-acetylneuraminic acid (CMP-NeuAc) hydroxylase activity. We previously cloned a cDNA for mouse CMP-NeuAc hydroxylase and found that the human genome contains a homologue. We report here the molecular basis for the absence of NeuGc in humans. We cloned a cDNA for human CMP-NeuAc hydroxylase from a HeLa cell cDNA library. The cDNA encodes a 486-amino acid protein, and its deduced amino acid sequence lacks a domain corresponding to the N-terminal 104 amino acids of the mouse CMP-NeuAc hydroxylase protein, although the human protein is highly identical (93%) to the rest of the mouse hydroxylase protein. The N-terminal truncation of the human hydroxylase is caused by deletion of a 92-base pair-long exon in human genomic DNA. The human hydroxylase expressed in COS-7 cells exhibited no enzymatic activity, and a mouse hydroxylase mutant, which lacks the N-terminal domain, was also inactive. A chimera composed of the human hydroxylase and the N-terminal domain of the mouse hydroxylase displayed the enzyme activity. These results indicate that the human homologue of CMP-NeuAc hydroxylase is inactive because it lacks an N-terminal domain that is essential for enzyme activity. The absence of NeuGc in human glycoconjugates is due to a partial deletion in the gene that encodes CMP-NeuAc hydroxylase.

Microbial secondary metabolism constitutes a rich source of antibiotics, chemotherapeutics, insecticides and other high-value chemicals. Genome mining of gene clusters that encode the biosynthetic pathways for these metabolites has become a key methodology for novel compound discovery. In 2011, we introduced antiSMASH, a web server and stand-alone tool for the automatic genomic identification and analysis of biosynthetic gene clusters, available at http://antismash.secondarymetabolites.org. Here, we present version 3.0 of antiSMASH, which has undergone major improvements. A full integration of the recently published ClusterFinder algorithm now allows using this probabilistic algorithm to detect putative gene clusters of unknown types. Also, a new dereplication variant of the ClusterBlast module now identifies similarities of identified clusters to any of 1172 clusters with known end products. At the enzyme level, active sites of key biosynthetic enzymes are now pinpointed through a curated pattern-matching procedure and Enzyme Commission numbers are assigned to functionally classify all enzyme-coding genes. Additionally, chemical structure prediction has been improved by incorporating polyketide reduction states. Finally, in order for users to be able to organize and analyze multiple antiSMASH outputs in a private setting, a new XML output module allows offline editing of antiSMASH annotations within the Geneious software.

Multiple matching algorithms are used to locate the occurrences of patterns from a finite pattern set in a large input string. Aho–Corasick and Wu–Manber, two of the most well known algorithms for multiple matching require an increased computing power, particularly in cases where large-size datasets must be processed, as is common in computational biology applications. Over the past years, Graphics Processing Units (GPUs) have evolved to powerful parallel processors outperforming Central Processing Units (CPUs) in scientific calculations. Moreover, multiple GPUs can be used in parallel, forming hybrid computer cluster configurations to achieve an even higher processing throughput. This paper evaluates the speedup of the parallel implementation of the Aho–Corasick and Wu–Manber algorithms on a hybrid GPU cluster, when used to process a snapshot of the Expressed Sequence Tags of the human genome and for different problem parameters.

The recent release of the first draft of the human genome provides an unprecedented opportunity to integrate human genes and their functions in a complete positional context. However, at least three significant technical hurdles remain: first, to assemble a complete and nonredundant human transcript index; second, to accurately place the individual transcript indices on the human genome; and third, to functionally annotate all human genes. Here, we report the extension of the UNIGENE database through the assembly of its sequence clusters into nonredundant sequence contigs. Each resulting consensus was aligned to the human genome draft. A unique location for each transcript within the human genome was determined by the integration of the restriction fingerprint, assembled genomic contig, and radiation hybrid (RH) maps. A total of 59,500 UNIGENE clusters were mapped on the basis of at least three independent criteria as compared with the 30,000 human genes/ESTs currently mapped in Genemap'99. Finally, the extension of the human transcript consensus in this study enabled a greater number of putative functional assignments than the 11,000 annotated entries in UNIGENE. This study reports a draft physical map with annotations for a majority of the human transcripts, called the Human Index of Nonredundant Transcripts (HINT). Such information can be immediately applied to the discovery of new genes and the identification of candidate genes for positional cloning.

Recent technological advances allowed an exponential expansion of biological sequence data, and the extraction of meaningful information through Machine Learning (ML) algorithms. This knowledge improved the understanding of the mechanisms related to several fatal diseases, e.g., Cancer and COVID-19, helping to develop innovative solutions, such as CRISPR-based gene editing, coronavirus vaccine, and precision medicine. These advances benefit our society and economy, directly impacting peoples lives in various areas, such as health care, drug discovery, forensic analysis, and food analysis. Nevertheless, ML approaches applied to biological data require representative, quantitative, and informative features. Necessarily, as many ML algorithms can handle only numerical data, sequences need to be translated into a feature vector. This process, known as feature extraction, is a fundamental step for the elaboration of high-quality ML-based models in bioinformatics, by allowing the feature engineering stage, with the design and selection of suitable features. Feature engineering, ML algorithm selection, and hyperparameter tuning are often manual and time-consuming processes, requiring extensive domain knowledge, and performed manually by a human expert. To deal with this problem, we developed a new package, BioAutoML, which automatically runs an end-to-end ML pipeline. BioAutoML extracts numerical and informative features from biological sequence databases, automating feature selection, recommendation of ML algorithm(s), and tuning of hyperparameters, using Automated ML (AutoML). BioAutoML has two components, divided into four modules, (1) automated feature engineering (feature extraction and selection modules) and (2) Metalearning (algorithm recommendation and hyperparameter tuning modules). Our experimental results, assessing the relevance of our proposal, indicate robust results for different problem domains, such as SARS-CoV-2, anticancer peptides, HIV sequences, and non-coding RNAs. According to our systematic review, our proposal is innovative compared to available studies in the literature, being the first study to propose automated feature engineering and metalearning for biological sequences. BioAutoML has a high potential to significantly reduce the expertise required to use ML pipelines, aiding researchers in combating diseases, particularly in low- and middle-income countries. This initiative can provide biologists, physicians, epidemiologists, and other stakeholders with an opportunity for widespread use of these techniques to enhance the health and well-being of their communities.

It is easy for today's students and researchers to believe that modern bioinformatics emerged recently to assist next-generation sequencing data analysis. However, the very beginnings of bioinformatics occurred more than 50 years ago, when desktop computers were still a hypothesis and DNA could not yet be sequenced. The foundations of bioinformatics were laid in the early 1960s with the application of computational methods to protein sequence analysis (notably, de novo sequence assembly, biological sequence databases and substitution models). Later on, DNA analysis also emerged due to parallel advances in (i) molecular biology methods, which allowed easier manipulation of DNA, as well as its sequencing, and (ii) computer science, which saw the rise of increasingly miniaturized and more powerful computers, as well as novel software better suited to handle bioinformatics tasks. In the 1990s through the 2000s, major improvements in sequencing technology, along with reduced costs, gave rise to an exponential increase of data. The arrival of 'Big Data' has laid out new challenges in terms of data mining and management, calling for more expertise from computer science into the field. Coupled with an ever-increasing amount of bioinformatics tools, biological Big Data had (and continues to have) profound implications on the predictive power and reproducibility of bioinformatics results. To overcome this issue, universities are now fully integrating this discipline into the curriculum of biology students. Recent subdisciplines such as synthetic biology, systems biology and whole-cell modeling have emerged from the ever-increasing complementarity between computer science and biology.