Dataset Information

Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

ABSTRACT: Telomere length is a risk factor in disease and the dynamics of telomere length are crucial to our understanding of cell replication and vitality. The proliferation of whole genome sequencing represents an unprecedented opportunity to glean new insights into telomere biology on a previously unimaginable scale. To this end, a number of approaches for estimating telomere length from whole-genome sequencing data have been proposed. Here we present Telomerecat, a novel approach to the estimation of telomere length. Previous methods have been dependent on the number of telomeres present in a cell being known, which may be problematic when analysing aneuploid cancer data and non-human samples. Telomerecat is designed to be agnostic to the number of telomeres present, making it suited for the purpose of estimating telomere length in cancer studies. Telomerecat also accounts for interstitial telomeric reads and presents a novel approach to dealing with sequencing errors. We show that Telomerecat performs well at telomere length estimation when compared to leading experimental and computational methods. Furthermore, we show that it detects expected patterns in longitudinal data, repeated measurements, and cross-species comparisons. We also apply the method to a cancer cell data, uncovering an interesting relationship with the underlying telomerase genotype.

SUBMITTER: Farmery JHR

PROVIDER: S-EPMC5778012 | biostudies-literature | 2018 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

Farmery James H R JHR Smith Mike L ML Lynch Andy G AG

Scientific reports 20180122 1

Telomere length is a risk factor in disease and the dynamics of telomere length are crucial to our understanding of cell replication and vitality. The proliferation of whole genome sequencing represents an unprecedented opportunity to glean new insights into telomere biology on a previously unimaginable scale. To this end, a number of approaches for estimating telomere length from whole-genome sequencing data have been proposed. Here we present Telomerecat, a novel approach to the estimation of ...[more]

PMID: 29358629

Similar Datasets

Project description:BackgroundAmyotrophic lateral sclerosis (ALS) is a neurodegenerative disease characterized by the loss of upper and lower motor neurons, leading to progressive weakness of voluntary muscles, with death following from neuromuscular respiratory failure, typically within 3 to 5 years. There is a strong genetic contribution to ALS risk. In 10% or more, a family history of ALS or frontotemporal dementia is obtained, and the Mendelian genes responsible for ALS in such families have now been identified in about 50% of cases. Only about 14% of apparently sporadic ALS is explained by known genetic variation, suggesting that other forms of genetic variation are important. Telomeres maintain DNA integrity during cellular replication, differ between sexes, and shorten naturally with age. Sex and age are risk factors for ALS and we therefore investigated telomere length in ALS.MethodsSamples were from Project MinE, an international ALS whole genome sequencing consortium that includes phenotype data. For validation we used donated brain samples from motor cortex from people with ALS and controls. Ancestry and relatedness were evaluated by principal components analysis and relationship matrices of DNA microarray data. Whole genome sequence data were from Illumina HiSeq platforms and aligned using the Isaac pipeline. TelSeq was used to quantify telomere length using whole genome sequence data. We tested the association of telomere length with ALS and ALS survival using Cox regression.ResultsThere were 6,580 whole genome sequences, reducing to 6,195 samples (4,315 from people with ALS and 1,880 controls) after quality control, and 159 brain samples (106 ALS, 53 controls). Accounting for age and sex, there was a 20% (95% CI 14%, 25%) increase of telomere length in people with ALS compared to controls (p = 1.1 × 10-12), validated in the brain samples (p = 0.03). Those with shorter telomeres had a 10% increase in median survival (p = 5.0×10-7). Although there was no difference in telomere length between sporadic ALS and familial ALS (p=0.64), telomere length in 334 people with ALS due to expanded C9orf72 repeats was shorter than in those without expanded C9orf72 repeats (p = 5.0×10-4).DiscussionAlthough telomeres shorten with age, longer telomeres are a risk factor for ALS and worsen prognosis. Longer telomeres are associated with ALS.

Project description:Epistasis plays an essential rule in understanding the regulation mechanisms and is an essential component of the genetic architecture of the gene expressions. However, interaction analysis of gene expressions remains fundamentally unexplored due to great computational challenges and data availability. Due to variation in splicing, transcription start sites, polyadenylation sites, post-transcriptional RNA editing across the entire gene, and transcription rates of the cells, RNA-seq measurements generate large expression variability and collectively create the observed position level read count curves. A single number for measuring gene expression which is widely used for microarray measured gene expression analysis is highly unlikely to sufficiently account for large expression variation across the gene. Simultaneously analyzing epistatic architecture using the RNA-seq and whole genome sequencing (WGS) data poses enormous challenges.We develop a nonlinear functional regression model (FRGM) with functional responses where the position-level read counts within a gene are taken as a function of genomic position, and functional predictors where genotype profiles are viewed as a function of genomic position, for epistasis analysis with RNA-seq data. Instead of testing the interaction of all possible pair-wises SNPs, the FRGM takes a gene as a basic unit for epistasis analysis, which tests for the interaction of all possible pairs of genes and use all the information that can be accessed to collectively test interaction between all possible pairs of SNPs within two genome regions.By large-scale simulations, we demonstrate that the proposed FRGM for epistasis analysis can achieve the correct type 1 error and has higher power to detect the interactions between genes than the existing methods. The proposed methods are applied to the RNA-seq and WGS data from the 1000 Genome Project. The numbers of pairs of significantly interacting genes after Bonferroni correction identified using FRGM, RPKM and DESeq were 16,2361, 260 and 51, respectively, from the 350 European samples.The proposed FRGM for epistasis analysis of RNA-seq can capture isoform and position-level information and will have a broad application. Both simulations and real data analysis highlight the potential for the FRGM to be a good choice of the epistatic analysis with sequencing data.

Project description:BACKGROUND:Mixed infections of Mycobacterium tuberculosis and antibiotic heteroresistance continue to complicate tuberculosis (TB) diagnosis and treatment. Detection of mixed infections has been limited to molecular genotyping techniques, which lack the sensitivity and resolution to accurately estimate the multiplicity of TB infections. In contrast, whole genome sequencing offers sensitive views of the genetic differences between strains of M. tuberculosis within a sample. Although metagenomic tools exist to classify strains in a metagenomic sample, most tools have been developed for more divergent species, and therefore cannot provide the sensitivity required to disentangle strains within closely related bacterial species such as M. tuberculosis. Here we present QuantTB, a method to identify and quantify individual M. tuberculosis strains in whole genome sequencing data. QuantTB uses SNP markers to determine the combination of strains that best explain the allelic variation observed in a sample. QuantTB outputs a list of identified strains, their corresponding relative abundances, and a list of drugs for which resistance-conferring mutations (or heteroresistance) have been predicted within the sample. RESULTS:We show that QuantTB has a high degree of resolution and is capable of differentiating communities differing by less than 25 SNPs and identifying strains down to 1× coverage. Using simulated data, we found QuantTB outperformed other metagenomic strain identification tools at detecting strains and quantifying strain multiplicity. In a real-world scenario, using a dataset of 50 paired clinical isolates from a study of patients with either reinfections or relapses, we found that QuantTB could detect mixed infections and reinfections at rates concordant with a manually curated approach. CONCLUSION:QuantTB can determine infection multiplicity, identify hetero-resistance patterns, enable differentiation between relapse and re-infection, and clarify transmission events across seemingly unrelated patients - even in low-coverage (1×) samples. QuantTB outperforms existing tools and promises to serve as a valuable resource for both clinicians and researchers working with clinical TB samples.

Dataset Information

Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

Publications

Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets