Dataset Information

GRACy: A tool for analysing human cytomegalovirus sequence data.

ABSTRACT: Modern DNA sequencing has instituted a new era in human cytomegalovirus (HCMV) genomics. A key development has been the ability to determine the genome sequences of HCMV strains directly from clinical material. This involves the application of complex and often non-standardized bioinformatics approaches to analysing data of variable quality in a process that requires substantial manual intervention. To relieve this bottleneck, we have developed GRACy (Genome Reconstruction and Annotation of Cytomegalovirus), an easy-to-use toolkit for analysing HCMV sequence data. GRACy automates and integrates modules for read filtering, genotyping, genome assembly, genome annotation, variant analysis, and data submission. These modules were tested extensively on simulated and experimental data and outperformed generic approaches. GRACy is written in Python and is embedded in a graphical user interface with all required dependencies installed by a single command. It runs on the Linux operating system and is designed to allow the future implementation of a cross-platform version. GRACy is distributed under a GPL 3.0 license and is freely available at https://bioinformatics.cvr.ac.uk/software/ with the manual and a test dataset.

SUBMITTER: Camiolo S

PROVIDER: S-EPMC7816668 | biostudies-literature | 2021 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

GRACy: A tool for analysing human cytomegalovirus sequence data.

Camiolo Salvatore S Suárez Nicolás M NM Chalka Antonia A Venturini Cristina C Breuer Judith J Davison Andrew J AJ

Virus evolution 20201230 1

Modern DNA sequencing has instituted a new era in human cytomegalovirus (HCMV) genomics. A key development has been the ability to determine the genome sequences of HCMV strains directly from clinical material. This involves the application of complex and often non-standardized bioinformatics approaches to analysing data of variable quality in a process that requires substantial manual intervention. To relieve this bottleneck, we have developed GRACy (Genome Reconstruction and Annotation of Cyto ...[more]

PMID: 33505707

Similar Datasets

Project description:BackgroundShort-sequence repeats (SSRs) occur in both prokaryotic and eukaryotic DNA, inter- and intragenically, and may be exact or inexact copies. When heterogeneous SSRs are present in a given locus, we can take advantage of the pattern of different repeats to genotype strains based on the SSRs. Cataloguing and tracking these repeats can be difficult as diverse groups of researchers are involved in the identification of the repeats. Additionally, the task is error-prone when done manually.ResultsWe developed RepeatAnalyzer, a new software tool capable of tracking, managing, analysing and cataloguing SSRs and genotypes using Anaplasma marginale as a model species. RepeatAnalyzer's analysis capability includes novel metrics for measuring regional genetic diversity (corresponding to variety and regularity of SSR occurrence). As a part of its visualization capabilities, RepeatAnalyzer produces high quality maps of the geographic distribution of genotypes or SSRs over a region of interest. RepeatAnalyzer's repeat identification functionality was validated for all SSRs and genotypes reported in 21 publications, using 380 A. marginale isolates gathered from the five publications within that list that provided access to their isolates. The tool produced accurate genotyping results in every case. In addition, it uncovered a number of errors in the published literature: 11 cases where SSRs were misreported, 5 cases where two different SSRs had been given the same name, and 16 cases where two or more names had been given to a single SSR. The analysis and visualization functionalities of the tool are demonstrated using several examples.ConclusionsRepeatAnalyzer is a robust software tool that can be used for storing, managing, and analysing short-sequence repeats for the purpose of strain identification. The tool can be used for any set of SSRs regardless of species. When applied to A. marginale, our test case, we show that genotype lengths for a given region follow a normal distribution, while SSR frequencies follow a power-law-like distribution. Further, we find that over 90 % of repeats are 28 to 29 amino acids long, which is in agreement with conventional wisdom. Lastly, our analysis reveals that the most common edit distance is five or six, which is counter-intuitive since we expected that result to be closer to one, resulting from the simplest change from one repeat to another.

Project description:The mechanisms of action and regulation of the human cytomegalovirus (HCMV) lytic-phase DNA replicator, oriLyt, which spans more than 2 kbp in a structurally complex region near the middle of the unique long region (UL), are not understood. Because oriLyt is thought to be essential for promoting initiation of lytic DNA synthesis and may participate in regulating the switch between lytic and latent phases, we undertook a mutational study to better define its sequence requirements. Kanr gene cassette insertions located an oriLyt core region between nucleotides (nt) 91751 and 93299 that is necessary but not sufficient for replicator activity in transient assays. In contrast, insertions into auxiliary regions flanking either side of this core-also required for significant replicator activity-had little effect. To search for essential components within the core region, we made a series of overlapping, roughly 200-bp deletions, and qualitatively and quantitatively assessed the abilities of the resulting constructs to mediate replication. All but one of these deletions produced a significant (i.e., greater than twofold) loss of activity, arguing that sequences across this entire region contribute to replicator function. However, two particularly critical segments separated by a dispensable region, here called essential regions I and II, were identified. Within essential region I, which overlaps the previously identified early transcript SRT, two adjacent but nonoverlapping, roughly 200-bp deletions abolished detectable replication. No single element or motif from the left half of essential region I was found to be essential. Thus, essential region I probably promotes replication through the cooperation of multiple elements. However, four small deletions in the right half of essential region I, which included or lay adjacent to the conserved 31-nt oligopyrimidine tract (referred to as the Y block), abolished or virtually abolished oriLyt activity. Together, these results identify candidate oriLyt sequences within which molecular interactions essential for initiation of oriLyt-mediated DNA synthesis are likely to occur.

Dataset Information

GRACy: A tool for analysing human cytomegalovirus sequence data.

Publications

GRACy: A tool for analysing human cytomegalovirus sequence data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure