73 research outputs found

    Evaluating Maximum Likelihood Estimation Methods to Determine the Hurst Coefficient

    Get PDF
    A maximum likelihood estimation method implemented in S-PLUS (S-MLE) to estimate the Hurst coefficient (H) is evaluated. The Hurst coefficient, with 0.5\u3cHS-MLE was developed to estimate H for fractionally differenced (fd) processes. However, in practice it is difficult to distinguish between fd processes and fractional Gaussian noise (fGn) processes. Thus, the method is evaluated for estimating H for both fd and fGn processes. S-MLE gave biased results of H for fGn processes of any length and for fd processes of lengths less than 210. A modified method is proposed to correct for this bias. It gives reliable estimates of H for both fd and fGn processes of length greater than or equal to 211

    Cost-Effective Cloud Computing: A Case Study Using the Comparative Genomics Tool, Roundup

    Get PDF
    Background Comparative genomics resources, such as ortholog detection tools and repositories are rapidly increasing in scale and complexity. Cloud computing is an emerging technological paradigm that enables researchers to dynamically build a dedicated virtual cluster and may represent a valuable alternative for large computational tools in bioinformatics. In the present manuscript, we optimize the computation of a large-scale comparative genomics resource—Roundup—using cloud computing, describe the proper operating principles required to achieve computational efficiency on the cloud, and detail important procedures for improving cost-effectiveness to ensure maximal computation at minimal costs. Methods Utilizing the comparative genomics tool, Roundup, as a case study, we computed orthologs among 902 fully sequenced genomes on Amazon's Elastic Compute Cloud. For managing the ortholog processes, we designed a strategy to deploy the web service, Elastic MapReduce, and maximize the use of the cloud while simultaneously minimizing costs. Specifically, we created a model to estimate cloud runtime based on the size and complexity of the genomes being compared that determines in advance the optimal order of the jobs to be submitted. Results We computed orthologous relationships for 245,323 genome-to-genome comparisons on Amazon's computing cloud, a computation that required just over 200 hours and cost $8,000 USD, at least 40% less than expected under a strategy in which genome comparisons were submitted to the cloud randomly with respect to runtime. Our cost savings projections were based on a model that not only demonstrates the optimal strategy for deploying RSD to the cloud, but also finds the optimal cluster size to minimize waste and maximize usage. Our cost-reduction model is readily adaptable for other comparative genomics tools and potentially of significant benefit to labs seeking to take advantage of the cloud as an alternative to local computing infrastructure

    Streaming Support for Data Intensive Cloud-Based Sequence Analysis

    Get PDF
    Cloud computing provides a promising solution to the genomics data deluge problem resulting from the advent of next-generation sequencing (NGS) technology. Based on the concepts of “resources-on-demand” and “pay-as-you-go”, scientists with no or limited infrastructure can have access to scalable and cost-effective computational resources. However, the large size of NGS data causes a significant data transfer latency from the client's site to the cloud, which presents a bottleneck for using cloud computing services. In this paper, we provide a streaming-based scheme to overcome this problem, where the NGS data is processed while being transferred to the cloud. Our scheme targets the wide class of NGS data analysis tasks, where the NGS sequences can be processed independently from one another. We also provide the elastream package that supports the use of this scheme with individual analysis programs or with workflow systems. Experiments presented in this paper show that our solution mitigates the effect of data transfer latency and saves both time and cost of computation

    Variant Characterization of a Representative Large Pedigree Suggests “Variant Risk Clusters” Convey Varying Predisposition of Risk to Lynch Syndrome

    Get PDF
    Recently, worldwide incidences of young adult aggressive colorectal cancer (CRC) have rapidly increased. Of these incidences diagnosed as familial Lynch syndrome (LS) CRC, outcomes are extremely poor. In this study, we seek novel familial germline variants from a large pedigree Tunisian family with 12 LS-affected individuals to identify putative germline variants associated with varying risk of LS. Whole-genome sequencing analysis was performed to identify known and novel germline variants shared between affected and non-affected pedigree members. SNPs, indels, and structural variants (SVs) were computationally identified, and their oncological influence was predicted using the Genetic Association of Complex Diseases and Disorders, OncoKB, and My Cancer Genome databases. Of 94 germline familial variants identified with predicted functional impact, 37 SNPs/indels were detected in 28 genes, 2 of which (MLH1 and PRH1-TAS2R14) have known association with CRC and 4 others (PPP1R13B, LAMA5, FTO, and NLRP14) have known association with non-CRC cancers. In addition, 48 of 57 identified SVs overlap with 43 genes. Three of these genes (RELN, IRS2, and FOXP1) have a known association with non-CRC digestive cancers and one (RRAS2) has a known association with non-CRC cancer. Our study identified 83 novel, predicted functionally impactful germline variants grouped in three “variant risk clusters” shared in three familiarly associated LS groups (high, intermediate and low risk). This variant characterization study demonstrates that large pedigree investigations provide important evidence supporting the hypothesis that different “variant risk clusters” can convey different mechanisms of risk and oncogenesis of LS-CRC even within the same pedigree

    Genotator: A disease-agnostic tool for genetic annotation of disease

    Get PDF
    <p>Abstract</p> <p>Background</p> <p>Disease-specific genetic information has been increasing at rapid rates as a consequence of recent improvements and massive cost reductions in sequencing technologies. Numerous systems designed to capture and organize this mounting sea of genetic data have emerged, but these resources differ dramatically in their disease coverage and genetic depth. With few exceptions, researchers must manually search a variety of sites to assemble a complete set of genetic evidence for a particular disease of interest, a process that is both time-consuming and error-prone.</p> <p>Methods</p> <p>We designed a real-time aggregation tool that provides both comprehensive coverage and reliable gene-to-disease rankings for any disease. Our tool, called Genotator, automatically integrates data from 11 externally accessible clinical genetics resources and uses these data in a straightforward formula to rank genes in order of disease relevance. We tested the accuracy of coverage of Genotator in three separate diseases for which there exist specialty curated databases, Autism Spectrum Disorder, Parkinson's Disease, and Alzheimer Disease. Genotator is freely available at <url>http://genotator.hms.harvard.edu</url>.</p> <p>Results</p> <p>Genotator demonstrated that most of the 11 selected databases contain unique information about the genetic composition of disease, with 2514 genes found in only one of the 11 databases. These findings confirm that the integration of these databases provides a more complete picture than would be possible from any one database alone. Genotator successfully identified at least 75% of the top ranked genes for all three of our use cases, including a 90% concordance with the top 40 ranked candidates for Alzheimer Disease.</p> <p>Conclusions</p> <p>As a meta-query engine, Genotator provides high coverage of both historical genetic research as well as recent advances in the genetic understanding of specific diseases. As such, Genotator provides a real-time aggregation of ranked data that remains current with the pace of research in the disease fields. Genotator's algorithm appropriately transforms query terms to match the input requirements of each targeted databases and accurately resolves named synonyms to ensure full coverage of the genetic results with official nomenclature. Genotator generates an excel-style output that is consistent across disease queries and readily importable to other applications.</p
    corecore