Resolving the Full Spectrum of Human Genome Variation using Linked-Reads
Ontology highlight
ABSTRACT:
INSTRUMENT(S): Illumina HiSeq X Ten, Illumina HiSeq 4000
ORGANISM(S): Homo Sapiens
SUBMITTER: 10X Genomics
PROVIDER: PRJEB28297 | EVA | 2018-08-21
REPOSITORIES: EVA
Similar Datasets
Project description:Large-scale population based analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. Standard short-read approaches, used primarily due to accuracy, throughput and costs, fail to give a complete picture of a genome. They struggle to identify large, balanced structural events, cannot access repetitive regions of the genome and fail to resolve the human genome into its two haplotypes. Here we describe an approach that retains long range information while harnessing the power of short reads. Starting from only 1ng of DNA, we produce barcoded short read libraries. The use of novel informatic approaches allows for the barcoded short reads to be associated with the long molecules of origin producing a novel datatype known as 'Linked-Reads'. This approach allows for simultaneous detection of small and large variants from a single Linked-Read library. We have previously demonstrated the utility of whole genome Linked-Reads (lrWGS) for performing diploid, de novo assembly of individual genomes. In this manuscript, we show the utility of reference based analysis using a single Linked-Read library for full spectrum genome analysis. We demonstrate the ability of Linked-Reads to reconstruct megabase scale haplotypes and to recover parts of the genome that are typically inaccessible to short reads, including phenotypically important genes such as STRC, SMN1 and SMN2. We demonstrate the ability of both lrWGS and Linked-Read Whole Exome Sequencing (lrWES) to identify complex structural variations, including balanced events, single exon deletions, and single exon duplications. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.
| EGAS00001003121 | EGA
Project description:<p>Recently developed methods that utilize partitioning of long genomic DNA fragments, and barcoding of shorter fragments derived from them, have succeeded in retaining long-range information in short sequencing reads. These so-called read cloud approaches represent a powerful, accurate, and cost-effective alternative to single-molecule long-read sequencing. We developed software, GROC-SVs, that takes advantage of read clouds for structural variant detection and assembly. We apply the method to two 10x Genomics data sets, one chromothriptic sarcoma with several spatially separated samples, and one breast cancer cell line, all Illumina-sequenced to high coverage. Comparison to short-fragment data from the same samples, and validation by mate-pair data from a subset of the sarcoma samples, demonstrate substantial improvement in specificity of breakpoint detection compared to short-fragment sequencing, at comparable sensitivity, and vice versa. The embedded long-range information also facilitates sequence assembly of a large fraction of the breakpoints; importantly, consecutive breakpoints that are closer than the average length of the input DNA molecules can be assembled together and their order and arrangement reconstructed, with some events exhibiting remarkable complexity. These features facilitated an analysis of the structural evolution of the sarcoma. In the chromothripsis, rearrangements occurred before copy number amplifications, and using the phylogenetic tree built from point mutation data, we show that single nucleotide variants and structural variants are not correlated. We predict significant future advances in structural variant science using 10x data analyzed with GROC-SVs and other read cloud-specific methods.</p>
| phs001255 | dbGaP
Project description:Pioneering studies (PXD014844) have identified many interesting molecules in tick saliva by LC-MS/MS proteomics, but the protein databases used to assign mass spectra were based on short Illumina reads of the Amblyomma americanum transcriptome and may not have captured the diversity and complexity of longer transcripts. Here we apply long-read Pacific Bioscience technologies to complement the previously reported short-read Illumina transcriptome-based proteome in an effort to increase spectrum assignments. Our dataset reveals a small increase in assignable spectra to supplement the previously released short-read transcriptome-based proteome.
2022-05-31 | PXD033870 | Pride
Project description:Pioneering studies (PXD014844) have identified many interesting molecules by LC-MS/MS proteomics, but the protein databases used to assign mass spectra were based on short Illumina reads of the Amblyomma americanum transcriptome and may not have captured the diversity and complexity of longer transcripts. Here we apply long-read Pacific Bioscience technologies to complement the previously reported short-read Illumina transcriptome-based proteome in an effort to increase spectrum assignments. Our dataset reveals a small increase in assignable spectra to supplement previously released short-read transcriptome-based proteome.
2020-01-23 | PXD014844 | Pride
Project description:Stable ETV5 knockdown SH-SY5Y were used to investigate the role of ETV5 in vitro and in vivo. Libraries for mRNA sequencing were prepared using the TruSeq stranded mRNA sample prep kit (Illumina), involving polyA-selection, fragmentation, adapter ligation, reverse transcription and PCR amplification. Libraries were quantified using the Qubit 2.0 Fluorometer prior to sequencing with read length of 75 bp on a NextSeq 500 sequencer (Illumina). On average 24.7 M high-quality reads were generated per sample, with a minimal read count of 20.0 M reads for sample 2015_027_009 (sample L105). For each sample, gene-level read counts and KPKM-values were generated with Sailfish (Patro et al. 2014). Read counts are the number of reads overlapping a given feature such as a gene. KPKM-values (K-mers Per Kilobase of exon per Million reads) represent read counts normalized for sequencing depth and transcript length. KPKM-values allow you to compare genes across and within samples.
2020-05-01 | E-MTAB-6713 | biostudies-arrayexpress
Project description:In this study we generated 5'P libraries in wt and RNAse III mutant strains, grown to exponential and stationary phases. Libraries that retain short RNA fragments were also generated in both growth phases. After sequencing by Illumina NextSeq 500 system, reads were mapped to E. coli genome NC_000913.3. By comparing the read start counts per position in the wt and mutant strain libraries we identified the cleavage sites of RNase III.
2018-09-30 | E-MTAB-6704 | biostudies-arrayexpress
Project description:Purpose: We aimed to identify the targets of the RNA binding protein ZFP36L1 in thymoctes. Methods: Total naïve thymocytes from control or DCKO mice were treated with UV light to crosslink RNA and proteins, then RNA-protein complexes were pulled down with anti-ZFP36L1. RNA was extracted and used to make cDNS libraries that were then sequenced by MiSeq 150bp single-end read (Sample 1) and HiSeq2500 RapidRun 50bp single-end read (sample 2, 3). Results: Sample demultiplexing was performed by identification of the 3 known bases of the 7 bases barcode introduced in the 5’ end of the read by the RCLIP primer. The remaining four random bases were used to remove PCR duplicate reads. Reads were trimmed to remove any adaptor sequence and barcodes before mapping reads to genome GRCm38 using Bowtie. After read mapping, the single-nucleotide at position -1 was annotated as unique ZFP36L1 crosslink site. Identification of highly significant ZFP36L1 binding sites was performed using iCount to assign a FDR to each crosslink site. Conclusions: We identified ZFP36L1 binding sites in 8675 thymocyte mRNAs.
2016-09-05 | GSE79178 | GEO
Project description:Purpose: We aimed to identify the targets of the RNA binding protein ZFP36L1 in thymoctes. Methods: Total naïve thymocytes from control or DCKO mice were treated with UV light to crosslink RNA and proteins, then RNA-protein complexes were pulled down with anti-ZFP36L1. RNA was extracted and used to make cDNS libraries that were then sequenced by MiSeq 150bp single-end read (Sample 1) and HiSeq2500 RapidRun 50bp single-end read (sample 2, 3). Results: Sample demultiplexing was performed by identification of the 3 known bases of the 7 bases barcode introduced in the 5â end of the read by the RCLIP primer. The remaining four random bases were used to remove PCR duplicate reads. Reads were trimmed to remove any adaptor sequence and barcodes before mapping reads to genome GRCm38 using Bowtie. After read mapping, the single-nucleotide at position -1 was annotated as unique ZFP36L1 crosslink site. Identification of highly significant ZFP36L1 binding sites was performed using iCount to assign a FDR to each crosslink site. Conclusions: We identified ZFP36L1 binding sites in 8675 thymocyte mRNAs. Individual nucleotide resolution cross linking immunoprecipitation (iCLIP) of ZFP36L1-bound RNAs in total naive thymocytes from C57BL/6 mice.
2016-09-05 | E-GEOD-79178 | biostudies-arrayexpress
Project description:New tools for improved long-read transcript assembly and coalescence with its short-read counterpart are required. Using our short- and long-read measurements from different cell lines with spiked-in standards, we systematically compared key parameters and biases in the read alignment and assembly of transcripts. We report a cDNA synthesis artifact in long-read datasets that impacts the identity and quantitation of assembled transcripts. We developed a computational pipeline to strand long-read cDNA libraries that markedly improves assembly of transcripts from long-reads. Incorporating stranded long-reads in a new hybrid assembly approach, we demonstrate its efficacy for improved characterization of challenging lncRNA transcripts. Our workflow can be applied to a wide range of transcriptomics datasets for superior demarcation of transcript ends and refined isoform structure, which can enable better differential gene expression analyses and molecular manipulations of transcripts.
2023-10-14 | GSE215357 | GEO
Project description:New tools for improved long-read transcript assembly and coalescence with its short-read counterpart are required. Using our short- and long-read measurements from different cell lines with spiked-in standards, we systematically compared key parameters and biases in the read alignment and assembly of transcripts. We report a cDNA synthesis artifact in long-read datasets that impacts the identity and quantitation of assembled transcripts. We developed a computational pipeline to strand long-read cDNA libraries that markedly improves assembly of transcripts from long-reads. Incorporating stranded long-reads in a new hybrid assembly approach, we demonstrate its efficacy for improved characterization of challenging lncRNA transcripts. Our workflow can be applied to a wide range of transcriptomics datasets for superior demarcation of transcript ends and refined isoform structure, which can enable better differential gene expression analyses and molecular manipulations of transcripts.
2023-10-14 | GSE215355 | GEO