Adam Yao

FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization

Single nucleotide polymorphism (SNP) prioritization based on the phenotypic risk is essential for association studies. Assessment of the risk requires access to a variety of heterogeneous biological databases and analytical tools. FASTSNP (function analysis and selection tool for single nucleotide polymorphisms) is a web server that allows users to efficiently identify and prioritize high-risk SNPs according to their phenotypic risks and putative functional effects. A unique feature of FASTSNP is that the functional effect information used for SNP prioritization is always up-to-date, because FASTSNP extracts the information from 11 external web servers at query time using a team of web wrapper agents. Moreover, FASTSNP is extendable by simply deploying more Web wrapper agents. To validate the results of our prioritization, we analyzed 1569 SNPs from the SNP500Cancer database. The results show that SNPs with a high predicted risk exhibit low allele frequencies for the minor alleles, consistent with a well-known finding that a strong selective pressure exists for functional polymorphisms. We have been using FASTSNP for 2 years and FASTSNP enables us to discover a novel promoter polymorphism. FASTSNP is available at .

PrimerZ: streamlined primer design for promoters, exons and human SNPs

PrimerZ (http://genepipe.ngc.sinica.edu.tw/primerz/) is a web application dedicated primarily to primer design for genes and human SNPs. PrimerZ accepts genes by gene name or Ensembl accession code, and SNPs by dbSNP rs or AFFY_Probe IDs. The promoter and exon sequence information of all gene transcripts fetched from the Ensembl database (http://www.ensembl.org) are processed before being passed on to Primer3 (http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) for individual primer design. All results returned from Primer 3 are organized and integrated in a specially designed web page for easy browsing. Besides the web page presentation, csv text file export is also provided for enhanced user convenience. PrimerZ automates highly standard but tedious gene primer design to improve the success rate of PCR experiments. More than 2000 primers have been designed with PrimerZ at our institute since 2004 and the success rate is over 70%. The addition of several new features has made PrimerZ even more useful to the research community in facilitating primer design for promoters, exons and SNPs.

An optimal algorithm for maximum-sum segment and its application in bioinformatics

We study a fundamental sequence algorithm arising from bioinformatics. Given two integers L and U and a sequence A of n numbers, the maximum-sum segment problem is to find a segment A[i, j] of A with L < j-i+1 ≤ U that maximizes A[i]+A[i+1]+...+A[j]. The problem finds applications in finding repeats, designing low complexity filter, and locating segments with rich C+G content for biomolecular sequences. The best known algorithm, due to Lin, Jiang, and Chao, runs in O(n) time, based upon a clever technique called left-negative decomposition for A. In the present paper, we present a new O(n)-time algorithm that bypasses the left-negative decomposition. As a result, our algorithm has the capability to handle the input sequence in an online manner, which is clearly an important feature to cope with genome-scale sequences. We also show how to exploit the sparsity in the input sequence: If A is representable in O(k) space in some format, then our algorithm runs in O(k) time. Moreover, practical implementation of our algorithm running on the rice genome helps us to identify a very long repeat structure in rice chromosome 1 that is previously unknown.

VarioWatch: providing large-scale and comprehensive annotations on human genomic variants in the next generation sequencing era

VarioWatch (http://genepipe.ncgm.sinica.edu.tw/variowatch/) has been vastly improved since its former publication GenoWatch in the 2008 Web Server Issue. It is now at least 10 000-times faster in annotating a variant. Drastic speed increase, through complete re-design of its working mechanism, makes VarioWatch capable of annotating millions of human genomic variants generated from next generation sequencing in minutes, if not seconds. While using MegaQuery of VarioWatch to quickly annotate variants, users can apply various filters to retrieve a subgroup of variants according to the risk levels, interested regions, etc. that satisfy users’ requirements. In addition to performance leap, many new features have also been added, such as annotation on novel variants, functional analyses on splice sites and in/dels, detailed variant information in tabulated form, plus a risk level decision tree regarding the analyzed variant. Up to 1000 target variants can be visualized with our carefully designed Genome View, Gene View, Transcript View and Variation View. Two commonly used reference versions, NCBI build 36.3 and NCBI build 37.2, are supported. VarioWatch is unique in its ability to annotate comprehensively and efficiently millions of variants online, immediately delivering the results in real time, plus visualizes up to 1000 annotated variants.

Functional analysis of novel SNPs and mutations in human and mouse genomes

BackgroundWith the flood of information generated by the new generation of sequencing technologies, more efficient bioinformatics tools are needed for in-depth impact analysis of novel genomic variations. FANS (Functional Analysis of Novel SNPs) was developed to streamline comprehensive but tedious functional analysis steps into a few clicks and to offer a carefully designed presentation of results so researchers can focus more on thinking instead of typing and calculating.ResultsFANS http://fans.ngc.sinica.edu.tw/ harnesses the power of public information databases and powerful tools from six well established websites to enhance the efficiency of analysis of novel variations. FANS can process any point change in any coding region or GT-AG splice site to provide a clear picture of the disease risk of a prioritized variation by classifying splicing and functional alterations into one of nine risk subtypes with five risk levels.ConclusionFANS significantly simplifies the analysis operations to a four-step procedure while still covering all major areas of interest to researchers. FANS offers a convenient way to prioritize the variations and select the ones with most functional impact for validation. Additionally, the program offers a distinct improvement in efficiency over manual operations in our benchmark test.

Identification of differentially expressed microRNAs in human hepatocellular adenoma associated with type I glycogen storage disease: a potential utility as biomarkers

BackgroundIt is known that malignant transformation to hepatocellular carcinoma (HCC) occurs at a higher frequency in hepatocellular adenoma (HCA) from type I glycogen storage disease (GSD I) compared to HCA from other etiologies. In this study, we aimed to identify differentially expressed miRNAs in GSD Ia HCA as candidates that could serve as putative biomarkers for detection of GSD Ia HCA and/or risk assessment of malignant transformation.MethodsUtilizing massively parallel sequencing, the miRNA profiling was performed for paired adenomas and normal liver tissues from seven GSD Ia patients. Differentially expressed miRNAs were validated in liver tumor tissues, HCC cell lines and serum using quantitative RT-PCR.ResultsmiR-34a, miR-34a*, miR-224, miR-224*, miR-424, miR-452 and miR-455-5p were found to be commonly deregulated in GSD Ia HCA, general population HCA, and HCC cell lines at compatible levels. In comparison with GSD Ia HCA, the upregulation of miR-130b and downregulation of miR-199a-5p, miR-199b-5p, and miR-214 were more significant in HCC cell lines. Furthermore, serum level of miR-130b in GSD Ia patients with HCA was moderately higher than that in either GSD Ia patients without HCA or healthy individuals.ConclusionWe make the first observation of distinct miRNA deregulation in HCA associated with GSD Ia. We also provide evidence that miR-130b could serve as a circulating biomarker for detection of GSD Ia HCA. This work provides prominent candidate miRNAs worth evaluating as biomarkers for monitoring the development and progress of liver tumors in GSD Ia patients in the future.

GenoWatch: a disease gene mining browser for association study

A human gene association study often involves several genomic markers such as single nucleotide polymorphisms (SNPs) or short tandem repeat polymorphisms, and many statistically significant markers may be identified during the study. GenoWatch can efficiently extract up-to-date information about multiple markers and their associated genes in batch mode from many relevant biological databases in real-time. The comprehensive gene information retrieved includes gene ontology, function, pathway, disease, related articles in PubMed and so on. Subsequent SNP functional impact analysis and primer design of a target gene for re-sequencing can also be done in a few clicks. The presentation of results has been carefully designed to be as intuitive as possible to all users. The GenoWatch is available at the website http://genepipe.ngc.sinica.edu.tw/genowatch

LISE: a server using ligand-interacting and site-enriched protein triangles for prediction of ligand-binding sites.

LISE is a web server for a novel method for predicting small molecule binding sites on proteins. It differs from a number of servers currently available for such predictions in two aspects. First, rather than relying on knowledge of similar protein structures, identification of surface cavities or estimation of binding energy, LISE computes a score by counting geometric motifs extracted from sub-structures of interaction networks connecting protein and ligand atoms. These network motifs take into account spatial and physicochemical properties of ligand-interacting protein surface atoms. Second, LISE has now been more thoroughly tested, as, in addition to the evaluation we previously reported using two commonly used small benchmark test sets and targets of two community-based experiments on ligand-binding site predictions, we now report an evaluation using a large non-redundant data set containing >2000 protein–ligand complexes. This unprecedented test, the largest ever reported to our knowledge, demonstrates LISE’s overall accuracy and robustness. Furthermore, we have identified some hard to predict protein classes and provided an estimate of the performance that can be expected from a state-of-the-art binding site prediction server, such as LISE, on a proteome scale. The server is freely available at http://lise.ibms.sinica.edu.tw.

Abstract 518: Identification of differentially expressed microRNAs in human hepatocellular adenoma associated with type I glycogen storage disease: a potential utility as biomarkers

Proceedings: AACR Annual Meeting 2014; April 5-9, 2014; San Diego, CA Background: It is known that malignant transformation to hepatocellular carcinoma (HCC) occurs at a higher frequency in hepatocellular adenoma (HCA) from type I glycogen storage disease (GSD I) compared to HCA from other etiologies. In this study, we aimed to identify differentially expressed miRNAs in GSD Ia HCA as candidates that could serve as putative biomarkers for detection of GSD Ia HCA and/or risk assessment of malignant transformation. Methods: Utilizing massively parallel sequencing, the miRNA profiling was performed for paired adenomas and normal liver tissues from seven GSD Ia patients. Differentially expressed miRNAs were validated in liver tumor tissues, HCC cell lines and serum using quantitative RT-PCR. Results: miR-34a, miR-34a*, miR-224, miR-224*, miR-424, miR-452 and miR-455-5p were found to be commonly deregulated in GSD Ia HCA, general population HCA, and HCC cell lines at compatible levels. In comparison with GSD Ia HCA, the upregulation of miR-130b and downregulation of miR-199a-5p, miR-199b-5p, and miR-214 were more significant in HCC cell lines. Furthermore, serum level of miR-130b in GSD Ia patients with HCA was moderately higher than that in either GSD Ia patients without HCA or healthy individuals. Conclusion: We make the first observation of the distinct miRNA deregulation in HCA associated with GSD Ia. We also provide evidence that miR-130b could serve as a circulating biomarker for detection of GSD Ia HCA. This work provides prominent candidate miRNAs worth evaluating as biomarkers for monitoring the development and progress of liver tumors in GSD Ia patients in the future. Citation Format: Ling-Hui Li, Li-Ya Chiu, Priya S. Kishnani, Tzu-Po Chuang, Cheng-Yang Tang, Cheng-Yuan Liu, Deeksha Bali, Dwight Koeberl, Stephanie Austin, Keri Boyette, David A. Weinstein, Elaine Murphy, Adam Yao, Yuan-Tsong Chen. Identification of differentially expressed microRNAs in human hepatocellular adenoma associated with type I glycogen storage disease: a potential utility as biomarkers. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 518. doi:10.1158/1538-7445.AM2014-518

Qualiseq: Quality Genomic Sequence Retrieval

Genomic sequence analysis starts from a sequence with poor quality such as mis-orientation, incorrect location or contig order can hardly produce fruitful results in final assays. Unfortunately, more than often the circumstance is not found early enough until a lot of time and effort has been invested. To alleviate this situation, we have developed QualiSeq, a web application with a friendly interface providing information on sequence quality of a genomic region so that researchers can easily download much more reliable sequences for further analysis and assay design.