Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation

Long-read sequencing technologies substantially overcome the limitations of short-reads but have not been considered as a feasible replacement for population-scale projects, being a combination of too expensive, not scalable enough or too error-prone. Here we develop an efficient and scalable wet lab and computational protocol, Napu, for Oxford Nanopore Technologies long-read sequencing that seeks to address those limitations. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the National Institutes of Health Center for Alzheimer's and Related Dementias. Using a single PromethION flow cell, we can detect single nucleotide polymorphisms with F1-score comparable to Illumina short-read sequencing. Small indel calling remains difficult within homopolymers and tandem repeats, but achieves good concordance to Illumina indel calls elsewhere. Further, we can discover structural variants with F1-score on par with state-of-the-art de novo assembly methods. Our protocol phases small and structural variants at megabase scales and produces highly accurate, haplotype-specific methylation calls.

Scalable, accessible, and reproducible reference genome assembly and evaluation in Galaxy

Improvements in genome sequencing and assembly are enabling high-quality reference genomes for all species. However, the assembly process is still laborious, computationally and technically demanding, lacks standards for reproducibility, and is not readily scalable. Here we present the latest Vertebrate Genomes Project assembly pipeline and demonstrate that it delivers high-quality reference genomes at scale across a set of vertebrate species arising over the last ~500 million years. The pipeline is versatile and combines PacBio HiFi long-reads and Hi-C-based haplotype phasing in a new graph-based paradigm. Standardized quality control is performed automatically to troubleshoot assembly issues and assess biological complexities. We make the pipeline freely accessible through Galaxy, accommodating researchers even without local computational resources and enhanced reproducibility by democratizing the training and assembly process. We demonstrate the flexibility and reliability of the pipeline by assembling reference genomes for 51 vertebrate species from major taxonomic groups (fish, amphibians, reptiles, birds, and mammals).

Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation

Long-read sequencing technologies substantially overcome the limitations of short-reads but to date have not been considered as feasible replacement at scale due to a combination of being too expensive, not scalable enough, or too error-prone. Here, we develop an efficient and scalable wet lab and computational protocol for Oxford Nanopore Technologies (ONT) long-read sequencing that seeks to provide a genuine alternative to short-reads for large-scale genomics projects. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the NIH Center for Alzheimer's and Related Dementias (CARD). Using a single PromethION flow cell, we can detect SNPs with F1-score better than Illumina short-read sequencing. Small indel calling remains to be difficult inside homopolymers and tandem repeats, but is comparable to Illumina calls elsewhere. Further, we can discover structural variants with F1-score comparable to state-of the-art methods involving Pacific Biosciences HiFi sequencing and trio information (but at a lower cost and greater throughput). Using ONT based phasing, we can then combine and phase small and structural variants at megabase scales. Our protocol also produces highly accurate, haplotype-specific methylation calls. Overall, this makes large-scale long-read sequencing projects feasible; the protocol is currently being used to sequence thousands of brain-based genomes as a part of the NIH CARD initiative. We provide the protocol and software as open-source integrated pipelines for generating phased variant calls and assemblies.

Towards complete and error-free genome assemblies of all vertebrate species

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

An international virtual hackathon to build tools for the analysis of structural variants within species ranging from coronaviruses to vertebrates

In October 2020, 62 scientists from nine nations worked together remotely in the Second Baylor College of Medicine & DNAnexus hackathon, focusing on different related topics on Structural Variation, Pan-genomes, and SARS-CoV-2 related research.   The overarching focus was to assess the current status of the field and identify the remaining challenges. Furthermore, how to combine the strengths of the different interests to drive research and method development forward. Over the four days, eight groups each designed and developed new open-source methods to improve the identification and analysis of variations among species, including humans and SARS-CoV-2. These included improvements in SV calling, genotyping, annotations and filtering. Together with advancements in benchmarking existing methods. Furthermore, groups focused on the diversity of SARS-CoV-2. Daily discussion summary and methods are available publicly at  https://github.com/collaborativebioinformatics provides valuable insights for both participants and the research community.

An international virtual hackathon to build tools for the analysis of structural variants within species ranging from coronaviruses to vertebrates

In October 2020, 62 scientists from nine nations worked together remotely in the Second Baylor College of Medicine & DNAnexus hackathon, focusing on different related topics on Structural Variation, Pan-genomes, and SARS-CoV-2 related research.   The overarching focus was to assess the current status of the field and identify the remaining challenges. Furthermore, how to combine the strengths of the different interests to drive research and method development forward. Over the four days, eight groups each designed and developed new open-source methods to improve the identification and analysis of variations among species, including humans and SARS-CoV-2. These included improvements in SV calling, genotyping, annotations and filtering. Together with advancements in benchmarking existing methods. Furthermore, groups focused on the diversity of SARS-CoV-2. Daily discussion summary and methods are available publicly at  https://github.com/collaborativebioinformatics provides valuable insights for both participants and the research community.

HySA: a Hybrid Structural variant Assembly approach using next-generation and single-molecule sequencing technologies

Achieving complete, accurate, and cost-effective assembly of human genomes is of great importance for realizing the promise of precision medicine. The abundance of repeats and genetic variations in human genomes and the limitations of existing sequencing technologies call for the development of novel assembly methods that can leverage the complementary strengths of multiple technologies. We propose a Hybrid Structural variant Assembly (HySA) approach that integrates sequencing reads from next-generation sequencing and single-molecule sequencing technologies to accurately assemble and detect structural variants (SVs) in human genomes. By identifying homologous SV-containing reads from different technologies through a bipartite-graph-based clustering algorithm, our approach turns a whole genome assembly problem into a set of independent SV assembly problems, each of which can be effectively solved to enhance the assembly of structurally altered regions in human genomes. We used data generated from a haploid hydatidiform mole genome (CHM1) and a diploid human genome (NA12878) to test our approach. The result showed that, compared with existing methods, our approach had a low false discovery rate and substantially improved the detection of many types of SVs, particularly novel large insertions, small indels (10–50 bp), and short tandem repeat expansions and contractions. Our work highlights the strengths and limitations of current approaches and provides an effective solution for extending the power of existing sequencing technologies for SV discovery.

Abstract 136: Refining the molecular profile of colorectal tumors to expand prevention and treatment opportunities

The completion of The Cancer Genome Atlas (TCGA) project for colorectal cancer (CRC) is ushering in a new phase of identifying treatment strategies tailored to the molecular profile of each person's tumor. Precision medicine approaches to cancer treatment rely on the identification of molecular profiles that can be used to identify effective therapies and can be used in a targeted sequencing setting to make treatment decisions.

The initial TCGA colorectal effort included 276 samples and focused on integrating data from exome sequencing with genome-wide DNA copy number alterations (CNAs), DNA methylation, and mRNA and microRNA expression. Since then a total of 626 samples have been completed with the potential to refine CRC subtypes, identify novel mutated pathways, and further functional understanding. Such a large data set presents opportunities to identify new recurrent drug targets and to stratify patients into groups that are predictive of treatment response. However, large data sets also present substantial challenges, since hand-curation becomes intractable, while computational tools can be overwhelmed by hypermutation and copy number changes.

Here we present a comprehensive molecular analysis of all 626 TCGA colorectal cancer samples, including exome sequencing, CNAs, DNA methylation, and mRNA expression. For each data type, we identified recurrently altered genes. Using MutSigCV on 525 samples yielded 27 and 87 significantly mutated genes in non-hypermutated and hypermutated samples, respectively, a substantial increase over the 15 and 17 somatically recurrently mutated genes identified using MutSig in non-hypermutated and hypermutated samples, respectively, in the previously published TCGA colorectal study. For example, PTEN, a known tumor suppressor, was not reported as significantly recurrently mutated in the initial TCGA non-hypermutated set; however, it was in the larger non-hypermutated set, demonstrating the power of a larger data set for assessing the significance and relative frequency of mutations in the context of known subtypes.

In addition, we integrated the somatic mutation data, copy number data, LOH data, and hyper-methylation data to identify genes, like MLH1, that are recurrently disrupted by different mechanisms. We also considered somatic mutations that are likely gain-of-function mutations based on nonrandom clustering; and we used recurrent indels to identify loss-of-function drivers in samples positive for microsatellite instability (MSI). We further classified each sample using the previously identified subtypes – BRAF+, KRAS+, APC+, CTNNB1+ (beta-catenin+), TGFBR2/SMAD4+, PTEN+ and PIK3CA+, and R-spondin fusion positive, as well as CpG Island Methylator Phenotype (CIMP) and MSI - in order to refine the relevant molecular signatures driving CRC etiology and thereby prevention and treatment paradigms.

Citation Format: Catherine S. Grasso, Eve Shinbrot, Ming Yu, Max Liesersen, Mark Chaisson, Andrew Chan, Charles Connolly, James Dai, Margaret Du, Charles Fuchs, Levi Garraway, Marios Giannakis, Tabitha Harrison, Li Hsu, Jeroen Huyghe, Jasmine Mu, Shuji Ogino, Colin Pritchard, Stephen Salipante, Wei Sun, Syed H. Zaidi, Ni Zhao, William Grady, Ben Raphael, Thomas Hudson, David Wheeler, Ulrike Peters. Refining the molecular profile of colorectal tumors to expand prevention and treatment opportunities. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 136.

Extensive sequencing of seven human genomes to characterize benchmark reference materials

The Genome in a Bottle Consortium, hosted by the National Institute of Standards and Technology (NIST) is creating reference materials and data for human genome sequencing, as well as methods for genome comparison and benchmarking. Here, we describe a large, diverse set of sequencing data for seven human genomes; five are current or candidate NIST Reference Materials. The pilot genome, NA12878, has been released as NIST RM 8398. We also describe data from two Personal Genome Project trios, one of Ashkenazim Jewish ancestry and one of Chinese ancestry. The data come from 12 technologies: BioNano Genomics, Complete Genomics paired-end and LFR, Ion Proton exome, Oxford Nanopore, Pacific Biosciences, SOLiD, 10X Genomics GemCode WGS, and Illumina exome and WGS paired-end, mate-pair, and synthetic long reads. Cell lines, DNA, and data from these individuals are publicly available. Therefore, we expect these data to be useful for revealing novel information about the human genome and improving sequencing technologies, SNP, indel, and structural variant calling, and de novo assembly.

STAR: ultrafast universal RNA-seq aligner

MOTIVATION

Accurate alignment of high-throughput RNA-seq data is a challenging and yet unsolved problem because of the non-contiguous transcript structure, relatively short read lengths and constantly increasing throughput of the sequencing technologies. Currently available RNA-seq aligners suffer from high mapping error rates, low mapping speed, read length limitation and mapping biases.

RESULTS

To align our large (>80 billon reads) ENCODE Transcriptome RNA-seq dataset, we developed the Spliced Transcripts Alignment to a Reference (STAR) software based on a previously undescribed RNA-seq alignment algorithm that uses sequential maximum mappable seed search in uncompressed suffix arrays followed by seed clustering and stitching procedure. STAR outperforms other aligners by a factor of >50 in mapping speed, aligning to the human genome 550 million 2 × 76 bp paired-end reads per hour on a modest 12-core server, while at the same time improving alignment sensitivity and precision. In addition to unbiased de novo detection of canonical junctions, STAR can discover non-canonical splices and chimeric (fusion) transcripts, and is also capable of mapping full-length RNA sequences. Using Roche 454 sequencing of reverse transcription polymerase chain reaction amplicons, we experimentally validated 1960 novel intergenic splice junctions with an 80-90% success rate, corroborating the high precision of the STAR mapping strategy.

AVAILABILITY AND IMPLEMENTATION

STAR is implemented as a standalone C++ code. STAR is free open source software distributed under GPLv3 license and can be downloaded from http://code.google.com/p/rna-star/.

Paired de bruijn graphs: a novel approach for incorporating mate pair information into genome assemblers

The recent proliferation of next generation sequencing with short reads has enabled many new experimental opportunities but, at the same time, has raised formidable computational challenges in genome assembly. One of the key advances that has led to an improvement in contig lengths has been mate pairs, which facilitate the assembly of repeating regions. Mate pairs have been algorithmically incorporated into most next generation assemblers as various heuristic post-processing steps to correct the assembly graph or to link contigs into scaffolds. Such methods have allowed the identification of longer contigs than would be possible with single reads; however, they can still fail to resolve complex repeats. Thus, improved methods for incorporating mate pairs will have a strong effect on contig length in the future. Here, we introduce the paired de Bruijn graph, a generalization of the de Bruijn graph that incorporates mate pair information into the graph structure itself instead of analyzing mate pairs at a post-processing step. This graph has the potential to be used in place of the de Bruijn graph in any de Bruijn graph based assembler, maintaining all other assembly steps such as error-correction and repeat resolution. Through assembly results on simulated perfect data, we argue that this can effectively improve the contig sizes in assembly.

Fragment assembly with short reads

MOTIVATION

Current DNA sequencing technology produces reads of about 500-750 bp, with typical coverage under 10x. New sequencing technologies are emerging that produce shorter reads (length 80-200 bp) but allow one to generate significantly higher coverage (30x and higher) at low cost. Modern assembly programs and error correction routines have been tuned to work well with current read technology but were not designed for assembly of short reads.

RESULTS

We analyze the limitations of assembling reads generated by these new technologies and present a routine for base-calling in reads prior to their assembly. We demonstrate that while it is feasible to assemble such short reads, the resulting contigs will require significant (if not prohibitive) finishing efforts.

AVAILABILITY

Available from the web at http://www.cse.ucsd.edu/groups/bioinformatics/software.html

Disruption of redox homeostasis in tumor necrosis factor-induced apoptosis in a murine hepatocyte cell line

Tumor necrosis factor (TNF) is a mediator of the acute phase response in the liver and can initiate proliferation and cause cell death in hepatocytes. We investigated the mechanisms by which TNF causes apoptosis in hepatocytes focusing on the role of oxidative stress, antioxidant defenses, and mitochondrial damage. The studies were conducted in cultured AML12 cells, a line of differentiated murine hepatocytes. As is the case for hepatocytes in vivo, AML12 cells were not sensitive to cell death by TNF alone, but died by apoptosis when exposed to TNF and a small dose of actinomycin D (Act D). Morphological signs of apoptosis were not detected until 6 hours after the treatment and by 18 hours approximately 50% of the cells had died. Exposure of the cells to TNF+Act D did not block NFkappaB nuclear translocation, DNA binding, or its overall transactivation capacity. Induction of apoptosis was characterized by oxidative stress indicated by the loss of NAD(P)H and glutathione followed by mitochondrial damage that included loss of mitochondrial membrane potential, inner membrane structural damage, and mitochondrial condensation. These changes coincided with cytochrome C release and the activation of caspases-8, -9, and -3. TNF-induced apoptosis was dependent on glutathione levels. In cells with decreased levels of glutathione, TNF by itself in the absence of transcriptional blocking acted as an apoptotic agent. Conversely, the antioxidant alpha-lipoic acid, that protected against the loss of glutathione in cells exposed to TNF+Act D completely prevented mitochondrial damage, caspase activation, cytochrome C release, and apoptosis. The results demonstrate that apoptosis induced by TNF+Act D in AML12 cells involves oxidative injury and mitochondrial damage. As injury was regulated to a larger extent by the glutathione content of the cells, we suggest that the combination of TNF+Act D causes apoptosis because Act D blocks the transcription of genes required for antioxidant defenses.

Tumor necrosis factor induces DNA replication in hepatic cells through nuclear factor kappaB activation

Tumor necrosis factor (TNF) signaling through TNF receptor 1 (TNFR1) with downstream participation of nuclear factor kappaB (NFkappaB), interleukin 6 (IL-6), and signal transducers and activators of transcription 3 (STAT3) is required for initiation of liver regeneration. It is not known whether the proliferative effect of TNF on hepatocytes is direct or requires the participation of Kupffer cells, the liver resident macrophages. Moreover, it has not been determined whether NFkappaB activation is an essential step in TNF-induced proliferation. To answer these questions, we conducted studies in LE6 cells, a rat liver epithelial cell line with hepatocyte progenitor capacity. We report that TNF induces DNA replication in growth-arrested LE6 cells and that its effect involves the activation of NFkappaB and STAT3 and an increase in c-myc and IL-6 mRNAs. All of these effects, which mimic the events that initiate liver regeneration in vivo, are blocked if NFKB activation is inhibited by expression of a dominant-inhibitor IkappaBalpha mutant (deltaN-IkappaBalpha). Although NFkappaB blockage by deltaN-IkappaBalpha causes caspase activation and massive death of cells stimulated by TNF, inhibition of NFkappaB and STAT3 binding by the serine protease inhibitor N-tosyl-L-phenylalanine chloromethyl ketone results in G0-G1 cell cycle arrest without death. We conclude that NFkappaB is an essential component of the TNF proliferative pathway and that TNF-induced changes in IL-6 mRNA, STAT3, and c-myc mRNA are dependent on NFkappaB activation. Blockage of NFkappaB inhibits TNF-induced proliferation but does not necessarily cause cell death.