GIMS-a data warehouse for storage and analysis of genome sequence and functional data

INTRODUCTION:Genome-wide association studies (GWAS) have identified loci associated with Alzheimer’s disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.METHODS:We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases=2,184, N controls=2,383) and targeted analyses in sub-populations using WGS data from the Alzheimer’s Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.RESULTS:Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.DISCUSSION:This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

BackgroundNew "next generation" DNA sequencing technologies offer individual researchers the ability to rapidly generate large amounts of genome sequence data at dramatically reduced costs. As a result, a need has arisen for new software tools for storage, management and analysis of genome sequence data. Although bioinformatic tools are available for the analysis and management of genome sequences, limitations still remain. For example, restrictions on the submission of data and use of these tools may be imposed, thereby making them unsuitable for sequencing projects that need to remain in-house or proprietary during their initial stages. Furthermore, the availability and use of next generation sequencing in industrial, governmental and academic environments requires biologist to have access to computational support for the curation and analysis of the data generated; however, this type of support is not always immediately available.ResultsTo address these limitations, we have developed DraGnET (Draft Genome Evaluation Tool). DraGnET is an open source web application which allows researchers, with no experience in programming and database management, to setup their own in-house projects for storing, retrieving, organizing and managing annotated draft and complete genome sequence data. The software provides a web interface for the use of BLAST, allowing users to perform preliminary comparative analysis among multiple genomes. We demonstrate the utility of DraGnET for performing comparative genomics on closely related bacterial strains. Furthermore, DraGnET can be further developed to incorporate additional tools for more sophisticated analyses.ConclusionsDraGnET is designed for use either by individual researchers or as a collaborative tool available through Internet (or Intranet) deployment. For genome projects that require genome sequencing data to initially remain proprietary, DraGnET provides the means for researchers to keep their data in-house for analysis using local programs or until it is made publicly available, at which point it may be uploaded to additional analysis software applications. The DraGnET home page is available at http://www.dragnet.cvm.iastate.edu and includes example files for examining the functionalities, a link for downloading the DraGnET setup package and a link to the DraGnET source code hosted with full documentation on SourceForge.

Modification of DNA bases plays important roles in the epigenetic regulation of eukaryotic gene expression. Among the different types of DNA methylation, 5-methylcytosine (5mC) is common in higher eukaryotes. Although bisulfite sequencing is the established detection method for this modification, newer methods, such as Oxford nanopore sequencing, have been developed as quick and reliable alternatives. An earlier study using sensitive liquid chromatography tandem mass spectrometry (LC-MS/MS) indicated the presence of 5mC at very low concentration in Saccharomyces cerevisiae. More recently, a comprehensive study of the yeast genome found 40 5mC sites using the computational tool Nanopolish on nanopore sequencing output raw data. In the present study, we are trying to validate the prediction of the 5mC modifications in yeast with Nanopolish and two other nanopore software tools, Tombo and DeepSignal. Using publicly available genome sequencing data, we compared the open-access computational tools, including Tombo, Nanopolish and DeepSignal, for predicting 5mC. Our results suggest that these tools are indeed capable of predicting DNA 5mC modifications at a specific location from Oxford nanopore sequencing data. We also predicted that 5mC present in the S. cerevisiae genome might be located predominantly at the RDN locus of chromosome 12.

In this paper author presents an experiment, which shows that it is possible to form a federation of data warehouses that may simulate effectively one, “super” data warehouse. There is no need to create complete ETL tool to load data from source data warehouses into one, dedicated data warehouse. Good relations between global schema and local schemas extracted during schema integration are indispensable to create an effective federation.

Because susceptibility to diseases increases with age, studying aging gains importance. Analyses of gene expression or sequence data, which have been indispensable for investigating aging, have been limited to studying genes and their protein products in isolation, ignoring their connectivities. However, proteins function by interacting with other proteins, and this is exactly what biological networks (BNs) model. Thus, analyzing the proteins' BN topologies could contribute to the understanding of aging. Current methods for analyzing systems-level BNs deal with their static representations, even though cells are dynamic. For this reason, and because different data types can give complementary biological insights, we integrate current static BNs with aging-related gene expression data to construct dynamic age-specific BNs. Then, we apply sensitive measures of topology to the dynamic BNs to study cellular changes with age. While global BN topologies do not significantly change with age, local topologies of a number of genes do. We predict such genes to be aging-related. We demonstrate credibility of our predictions by (i) observing significant overlap between our predicted aging-related genes and 'ground truth' aging-related genes; (ii) observing significant overlap between functions and diseases that are enriched in our aging-related predictions and those that are enriched in 'ground truth' aging-related data; (iii) providing evidence that diseases which are enriched in our aging-related predictions are linked to human aging; and (iv) validating our high-scoring novel predictions in the literature. Software executables are available upon request.

Two billion people are infected with Mycobacterium tuberculosis, leading to 10 million new cases of active tuberculosis and 1.5 million deaths annually. Universal access to drug susceptibility testing (DST) has become a World Health Organization priority. We previously developed a software tool, Mykrobe predictor, which provided offline species identification and drug resistance predictions for M. tuberculosis from whole genome sequencing (WGS) data. Performance was insufficient to support the use of WGS as an alternative to conventional phenotype-based DST, due to mutation catalogue limitations. Here we present a new tool, Mykrobe, which provides the same functionality based on a new software implementation. Improvements include i) an updated mutation catalogue giving greater sensitivity to detect pyrazinamide resistance, ii) support for user-defined resistance catalogues, iii) improved identification of non-tuberculous mycobacterial species, and iv) an updated statistical model for Oxford Nanopore Technologies sequencing data. Mykrobe is released under MIT license at https://github.com/mykrobe-tools/mykrobe. We incorporate mutation catalogues from the CRyPTIC consortium et al. (2018) and from Walker et al. (2015), and make improvements based on performance on an initial set of 3206 and an independent set of 5845 M. tuberculosis Illumina sequences. To give estimates of error rates, we use a prospectively collected dataset of 4362 M. tuberculosis isolates. Using culture based DST as the reference, we estimate Mykrobe to be 100%, 95%, 82%, 99% sensitive and 99%, 100%, 99%, 99% specific for rifampicin, isoniazid, pyrazinamide and ethambutol resistance prediction respectively. We benchmark against four other tools on 10207 (=5845+4362) samples, and also show that Mykrobe gives concordant results with nanopore data. We measure the ability of Mykrobe-based DST to guide personalized therapeutic regimen design in the context of complex drug susceptibility profiles, showing 94% concordance of implied regimen with that driven by phenotypic DST, higher than all other benchmarked tools.

BackgroundKnowledge of the genetic variation underlying Primary Immune Deficiency (PID) is increasing. Reanalysis of genome-wide sequencing data from undiagnosed patients with suspected PID may improve the diagnostic rate.MethodsWe included patients monitored at the Department of Infectious Diseases or the Child and Adolescent Department, Rigshospitalet, Denmark, for a suspected PID, who had been analysed previously using a targeted PID gene panel (457 PID-related genes) on whole exome- (WES) or whole genome sequencing (WGS) data. A literature review was performed to extend the PID gene panel used for reanalysis of single nucleotide variation (SNV) and small indels. Structural variant (SV) calling was added on WGS data.ResultsGenetic data from 94 patients (86 adults) including 36 WES and 58 WGS was reanalysed a median of 23 months after the initial analysis. The extended gene panel included 208 additional PID-related genes. Genetic reanalysis led to a small increase in the proportion of patients with new suspicious PID related variants of uncertain significance (VUS). The proportion of patients with a causal genetic diagnosis was constant. In total, five patients (5%, including three WES and two WGS) had a new suspicious PID VUS identified due to reanalysis. Among these, two patients had a variant added due to the expansion of the PID gene panel, and three patients had a variant reclassified to a VUS in a gene included in the initial PID gene panel. The total proportion of patients with PID related VUS, likely pathogenic, and pathogenic variants increased from 43 (46%) to 47 (50%), as one patient had a VUS detected in both initial- and reanalysis. In addition, we detected new suspicious SNVs and SVs of uncertain significance in PID candidate genes with unknown inheritance and/or as heterozygous variants in genes with autosomal recessive inheritance in 8 patients.ConclusionThese data indicate a possible diagnostic gain of reassessing WES/WGS data from patients with suspected PID. Reasons for the possible gain included improved knowledge of genotype-phenotype correlation, expanding the gene panel, and adding SV analyses. Future studies of genotype-phenotype correlations may provide additional knowledge on the impact of the new suspicious VUSs.

In the present era, genomics represent a crucial role in the field of life sciences. Advancements in genomics and the development of high-throughput techniques facilitate the characterization of a wide range of genes according to their functions such as regulation of genes, metabolic pathways, and their reconstruction. In the era of genomics we face the challenge of storage and analysis of a huge amount of important data. Even in this early stage of the era there are many commercial techniques and tools/software available to analyze next-generation sequencing (NGS) data. All of these programs can be used for many uses such as sequence alignments, polymorphism detection, and functional and structural comparative genomics. In this chapter, we focus on advances in bioinformatics and also computational biology in genome sequencing and NGS data analysis and on the potential applications for the efficient collection, storage, and analysis of the huge amount of genomic data generated by researchers and the information retrieved from different sources and web browsers in relation to NGS analysis.

We consider modeling and analysis of functional data arising in method comparison studies. The observed data consist of repeated measurements of a continuous variable obtained using multiple methods of measurement on a sample of subjects. The data are treated as multivariate functional data that are observed with noise at a common set of discrete time points which may vary from subject to subject. The proposed methodology uses functional principal components analysis within the framework of a mixed-effects model to represent the observations in terms of a small number of method-specific principal components. Two approaches for estimating the unknowns in the model, both adaptations of general techniques developed for multivariate functional principal components analysis, are presented. Bootstrapping is employed to get estimates of bias and covariance matrix of model parameter estimates. These in turn are used to compute confidence intervals for parameters and functions thereof, such as the measures of similarity and agreement between the measurement methods, that are necessary for data analysis. The estimation approaches are evaluated using simulation. The methodology is illustrated by analyzing two datasets.

Where classical epidemiology has proven to be inadequate for surveillance and control of foodborne pathogens, molecular epidemiology, using genomic typing methods, can add value. However, the analysis of whole genome sequencing (WGS) data varies widely and is not yet fully harmonised. We used genomic data on 494 Listeria monocytogenes isolates from ready-to-eat food products and food processing environments deposited in the strain collection of the German National Reference Laboratory to compare various procedures for WGS data analysis and to evaluate compatibility of results. Two different core genome multilocus sequence typing (cgMLST) schemes, different reference genomes in single nucleotide polymorphism (SNP) analysis and commercial as well as open-source software were compared. Correlation of allele distances from the different cgMLST approaches was high, ranging from 0.97 to 1, and unified thresholds yielded higher clustering concordance than scheme-specific thresholds. The number of detected SNP differences could be increased up to a factor of 3.9 using a specific reference genome compared with a general one. Additionally, specific reference genomes improved comparability of SNP analysis results obtained using different software tools. The use of a closed or a draft specific reference genome did not make a difference. The harmonisation of WGS data analysis will finally guarantee seamless data exchange, but, in the meantime, knowledge on threshold values that lead to comparable clustering of isolates by different methods may improve communication between laboratories. We therefore established a translation code between commonly applied cgMLST and SNP methods based on optimised clustering concordances. This code can work as a first filter to identify WGS-based typing matches resulting from different methods, which opens up a new perspective for data exchange and thereby accelerates time-critical analyses, such as in outbreak investigations.

Enhanced precision of epidemiological typing in clinically suspected nosocomial outbreaks is crucial. Our aim was to investigate whether single nucleotide polymorphism (SNP) analysis and core genome (cg) multilocus sequence typing (MLST) of whole genome sequencing (WGS) data would more reliably identify a nosocomial outbreak, compared to earlier molecular typing methods. Sixteen isolates from a nosocomial outbreak of ESBL E.coli ST-131 in southeastern Sweden and three control strains were subjected to WGS. Sequences were explored by SNP analysis and cgMLST. cgMLST clearly differentiated between the outbreak isolates and the control isolates (>1400 differences). All clinically identified outbreak isolates showed close clustering (≥2 allele differences), except for two isolates (>50 allele differences). These data confirmed that the isolates with >50 differing genes did not belong to the nosocomial outbreak. The number of SNPs within the outbreak was ≤7, whereas the two discrepant isolates had >700 SNPs. Two of the ESBL E.coli ST-131 isolates did not belong to the clinically identified outbreak. Our results illustrate the power of WGS in terms of resolution, which may avoid overestimation of patients belonging to outbreaks as judged from epidemiological data and previously employed molecular methods with lower discriminatory ability.

Although many genetic factors have been successfully identified for human diseases in genome-wide association studies, genes discovered to date only account for a small proportion of overall genetic contributions to many complex traits. Association studies have difficulty in detecting the remaining true genetic variants that are either common variants with weak allelic effects, or rare variants that have strong allelic effects but are weakly associated at the population level. In this work, we applied a goodness-of-fit test for detecting sets of common and rare variants associated with quantitative or binary traits by using whole genome sequencing data. This test has been proved optimal for detecting weak and sparse signals in the literature, which fits the requirements for targeting the genetic components of missing heritability. Furthermore, this p value-combining method allows one to incorporate different data and/or research results for meta-analysis. The method was used to simultaneously analyse the whole genome sequencing and genome-wide association studies data of Genetic Analysis Workshop 18 for detecting true genetic variants. The results show that goodness-of-fit test is comparable or better than the influential sequence kernel association test in many cases.

Sign, Wilcoxon and Mann-Whitney tests are nonparametric methods in one or two-sample problems. The nonparametric methods are alternatives used for testing hypothesis when the standard methods based on the Gaussianity assumption are not suitable to be applied. Recently, the functional data analysis (FDA) has gained relevance in statistical modeling. In FDA, each observation is a curve or function which usually is a realization of a stochastic process. In the literature of FDA, several methods have been proposed for testing hypothesis with samples coming from Gaussian processes. However, when this assumption is not realistic, it is necessary to utilize other approaches. Clustering and regression methods, among others, for non-Gaussian functional data have been proposed recently. In this paper, we propose extensions of the sign, Wilcoxon and Mann-Whitney tests to the functional data context as methods for testing hypothesis when we have one or two samples of non-Gaussian functional data. We use random projections to transform the functional problem into a scalar one, and then we proceed as in the standard case. Based on a simulation study, we show that the proposed tests have a good performance. We illustrate the methodology by applying it to a real data set.

Gene expression patterns assayed across development can offer key clues about a gene’s function and regulatory role. Drosophila melanogaster is ideal for such investigations as multiple individual and high-throughput efforts have captured the spatiotemporal patterns of thousands of embryonic expressed genes in the form of in situ images. FlyExpress (www.flyexpress.net), a knowledgebase based on a massive and unique digital library of standardized images and a simple search engine to find coexpressed genes, was created to facilitate the analytical and visual mining of these patterns. Here, we introduce the next generation of FlyExpress resources to facilitate the integrative analysis of sequence data and spatiotemporal patterns of expression from images. FlyExpress 7 now includes over 100,000 standardized in situ images and implements a more efficient, user-defined search algorithm to identify coexpressed genes via Genomewide Expression Maps (GEMs). Shared motifs found in the upstream 5′ regions of any pair of coexpressed genes can be visualized in an interactive dotplot. Additional webtools and link-outs to assist in the downstream validation of candidate motifs are also provided. Together, FlyExpress 7 represents our largest effort yet to accelerate discovery via the development and dispersal of new webtools that allow researchers to perform data-driven analyses of coexpression (image) and genomic (sequence) data.

This research implements two important main fields of the study which is functional data analysis and extreme value theory. The aims of this study are to convert the extreme rainfall data into functions and to propose method development for functional data analysis in extreme value theory. Functional data analysis is one of the techniques to convert or transform discrete observations into functions or curves. The daily rainfall data in Petaling Jaya from 1974 to 2003 will be tested on extreme value theory approach which is block maxima in order to practice the functional data analysis technique. The least squares and akaike criteria information method are used to smooth the fourier basis of the functional data analysis. A descriptive statistics which represents the behavior of a rainfall data for Petaling Jaya is analyzed and illustrated by a smooth curve from functional data analysis. This study indicates that there is an increase in rates of the rainfall in Petaling Jaya for the recent period compared to the last 30 years data both in the mean and standard deviation values.