Riza Daza

High-quality draft assemblies of mammalian genomes from massively parallel sequence data

Massively parallel DNA sequencing technologies are revolutionizing genomics by making it possible to generate billions of relatively short (~100-base) sequence reads at very low cost. Whereas such data can be readily used for a wide range of biomedical applications, it has proven difficult to use them to generate high-quality de novo genome assemblies of large, repeat-rich vertebrate genomes. To date, the genome assemblies generated from such data have fallen far short of those obtained with the older (but much more expensive) capillary-based sequencing approach. Here, we report the development of an algorithm for genome assembly, ALLPATHS-LG, and its application to massively parallel DNA sequence data from the human and mouse genomes, generated on the Illumina platform. The resulting draft genome assemblies have good accuracy, short-range contiguity, long-range connectivity, and coverage of the genome. In particular, the base accuracy is high (≥99.95%) and the scaffold sizes (N50 size = 11.5 Mb for human and 7.2 Mb for mouse) approach those obtained with capillary-based sequencing. The combination of improved sequencing technology and improved computational methods should now make it possible to increase dramatically the de novo sequencing of large genomes. The ALLPATHS-LG program is available at http://www.broadinstitute.org/science/programs/genome-biology/crd.

Multiplex single-cell profiling of chromatin accessibility by combinatorial cellular indexing

Chromatin state and the single cell Identifying the chromatin state of any single cell, which may or may not have a different function or represent different stages relative to others collected within any single culture, experiment, or tissue, has been challenging. Cusanovitch et al. skirted previously identified technological limitations to identify regions of accessible chromatin at single-cell resolution. Combinatorial cellular indexing, a strategy for multiplex barcoding of thousands of single cells per experiment, was successfully used to investigate the genome-wide chromatin accessibility landscape in each of over 15,000 single cells. Science, this issue p. 910 Combinatorial indexing can identify chromatin states at single-cell resolution. Technical advances have enabled the collection of genome and transcriptome data sets with single-cell resolution. However, single-cell characterization of the epigenome has remained challenging. Furthermore, because cells must be physically separated before biochemical processing, conventional single-cell preparatory methods scale linearly. We applied combinatorial cellular indexing to measure chromatin accessibility in thousands of single cells per assay, circumventing the need for compartmentalization of individual cells. We report chromatin accessibility profiles from more than 15,000 single cells and use these data to cluster cells on the basis of chromatin accessibility landscapes. We identify modules of coordinately regulated chromatin accessibility at the level of single cells both between and within cell types, with a scalable method that may accelerate progress toward a human cell atlas.

Decoding long nanopore sequencing reads of natural DNA

Nanopore sequencing of DNA is a single-molecule technique that may achieve long reads, low cost and high speed with minimal sample preparation and instrumentation. Here, we build on recent progress with respect to nanopore resolution and DNA control to interpret the procession of ion current levels observed during the translocation of DNA through the pore MspA. As approximately four nucleotides affect the ion current of each level, we measured the ion current corresponding to all 256 four-nucleotide combinations (quadromers). This quadromer map is highly predictive of ion current levels of previously unmeasured sequences derived from the bacteriophage phi X 174 genome. Furthermore, we show nanopore sequencing reads of phi X 174 up to 4,500 bases in length, which can be unambiguously aligned to the phi X 174 reference genome, and demonstrate proof-of-concept utility with respect to hybrid genome assembly and polymorphism detection. This work provides a foundation for nanopore sequencing of long, natural DNA strands.

A scalable, fully automated process for construction of sequence-ready barcoded libraries for 454

We present an automated, high throughput library construction process for 454 technology. Sample handling errors and cross-contamination are minimized via end-to-end barcoding of plasticware, along with molecular DNA barcoding of constructs. Automation-friendly magnetic bead-based size selection and cleanup steps have been devised, eliminating major bottlenecks and significant sources of error. Using this methodology, one technician can create 96 sequence-ready 454 libraries in 2 days, a dramatic improvement over the standard method.

Mutations causing medullary cystic kidney disease type 1 lie in a large VNTR in MUC1 missed by massively parallel sequencing

Although genetic lesions responsible for some mendelian disorders can be rapidly discovered through massively parallel sequencing of whole genomes or exomes, not all diseases readily yield to such efforts. We describe the illustrative case of the simple mendelian disorder medullary cystic kidney disease type 1 (MCKD1), mapped more than a decade ago to a 2-Mb region on chromosome 1. Ultimately, only by cloning, capillary sequencing and de novo assembly did we find that each of six families with MCKD1 harbors an equivalent but apparently independently arising mutation in sequence markedly under-represented in massively parallel sequencing data: the insertion of a single cytosine in one copy (but a different copy in each family) of the repeat unit comprising the extremely long (∼1.5–5 kb), GC-rich (>80%) coding variable-number tandem repeat (VNTR) sequence in the MUC1 gene encoding mucin 1. These results provide a cautionary tale about the challenges in identifying the genes responsible for mendelian, let alone more complex, disorders through massively parallel sequencing.

Fragment Length of Circulating Tumor DNA

Malignant tumors shed DNA into the circulation. The transient half-life of circulating tumor DNA (ctDNA) may afford the opportunity to diagnose, monitor recurrence, and evaluate response to therapy solely through a non-invasive blood draw. However, detecting ctDNA against the normally occurring background of cell-free DNA derived from healthy cells has proven challenging, particularly in non-metastatic solid tumors. In this study, distinct differences in fragment length size between ctDNAs and normal cell-free DNA are defined. Human ctDNA in rat plasma derived from human glioblastoma multiforme stem-like cells in the rat brain and human hepatocellular carcinoma in the rat flank were found to have a shorter principal fragment length than the background rat cell-free DNA (134–144 bp vs. 167 bp, respectively). Subsequently, a similar shift in the fragment length of ctDNA in humans with melanoma and lung cancer was identified compared to healthy controls. Comparison of fragment lengths from cell-free DNA between a melanoma patient and healthy controls found that the BRAF V600E mutant allele occurred more commonly at a shorter fragment length than the fragment length of the wild-type allele (132–145 bp vs. 165 bp, respectively). Moreover, size-selecting for shorter cell-free DNA fragment lengths substantially increased the EGFR T790M mutant allele frequency in human lung cancer. These findings provide compelling evidence that experimental or bioinformatic isolation of a specific subset of fragment lengths from cell-free DNA may improve detection of ctDNA.

Comprehensive single-cell transcriptional profiling of a multicellular organism

Sequencing each cell of the nematode Single-cell sequencing is challenging owing to the limited biological material available in an individual cell and the high cost of sequencing across multiple cells. Cao et al. developed a two-step combinatorial barcoding method to profile both single-cell and single-nucleus transcriptomes without requiring physical isolation of each cell. The authors profiled almost 50,000 single cells from an individual Caenorhabditis elegans larval stage and were able to identify and recover information from different, even rare, cell types. Science, this issue p. 661 Single-cell combinatorial indexing RNA sequencing achieves more than 50-fold cellular coverage of a developing nematode worm. To resolve cellular heterogeneity, we developed a combinatorial indexing strategy to profile the transcriptomes of single cells or nuclei, termed sci-RNA-seq (single-cell combinatorial indexing RNA sequencing). We applied sci-RNA-seq to profile nearly 50,000 cells from the nematode Caenorhabditis elegans at the L2 larval stage, which provided >50-fold “shotgun” cellular coverage of its somatic cell composition. From these data, we defined consensus expression profiles for 27 cell types and recovered rare neuronal cell types corresponding to as few as one or two cells in the L2 worm. We integrated these profiles with whole-animal chromatin immunoprecipitation sequencing data to deconvolve the cell type–specific effects of transcription factors. The data generated by sci-RNA-seq constitute a powerful resource for nematode biology and foreshadow similar atlases for other organisms.

Copy-Number Variation and False Positive Prenatal Aneuploidy Screening Results

Investigations of noninvasive prenatal screening for aneuploidy by analysis of circulating cell-free DNA (cfDNA) have shown high sensitivity and specificity in both high-risk and low-risk cohorts. However, the overall low incidence of aneuploidy limits the positive predictive value of these tests. Currently, the causes of false positive results are poorly understood. We investigated four pregnancies with discordant prenatal test results and found in two cases that maternal duplications on chromosome 18 were the likely cause of the discordant results. Modeling based on population-level copy-number variation supports the possibility that some false positive results of noninvasive prenatal screening may be attributable to large maternal copy-number variants. (Funded by the National Institutes of Health and others.).

Tsukamurella strandjordae sp. nov., a Proposed New Species Causing Sepsis

ABSTRACT We have isolated a gram-positive, weakly acid-alcohol-fast, irregular rod-shaped bacterium from cultures of blood from a 5-year-old girl with acute myelogenous leukemia. This isolate was compared with 14 other strains including reference strains of Tsukamurellaspecies by a polyphasic approach based on physiological and biochemical properties, whole-cell short-chain fatty acid and mycolic acid analyses, DNA-DNA hybridization, and sequencing of the 16S rRNA gene. This isolate represents a new taxon within the genusTsukamurella for which we propose the nameTsukamurella strandjordae sp. nov. Our study also revealed that Tsukamurella paurometabola ATCC 25938 represents a misnamed Tsukamurella inchonensis isolate and confirms thatTsukamurella wratislaviensis belongs to the genusRhodococcus.

Yeast Sterol Regulatory Element-binding Protein (SREBP) Cleavage Requires Cdc48 and Dsc5, a Ubiquitin Regulatory X Domain-containing Subunit of the Golgi Dsc E3 Ligase

Background: Yeast sterol regulatory element-binding protein (SREBP) proteolytic activation requires the Golgi Dsc E3 ligase and the proteasome. Results: Genetic selection identified additional genes required for SREBP activation. Conclusion: UBX domain protein Dsc5 and AAA ATPase Cdc48 are Dsc E3 ligase subunits required for SREBP proteolysis. Significance: Dsc5 and Cdc48 provide a mechanistic link between the Dsc E3 ligase and proteasome in SREBP proteolysis. Schizosaccharomyces pombe Sre1 is a membrane-bound transcription factor that controls adaptation to hypoxia. Like its mammalian homolog, sterol regulatory element-binding protein (SREBP), Sre1 activation requires release from the membrane. However, in fission yeast, this release occurs through a strikingly different mechanism that requires the Golgi Dsc E3 ubiquitin ligase complex and the proteasome. The mechanistic details of Sre1 cleavage, including the link between the Dsc E3 ligase complex and proteasome, are not well understood. Here, we present results of a genetic selection designed to identify additional components required for Sre1 cleavage. From the selection, we identified two new components of the fission yeast SREBP pathway: Dsc5 and Cdc48. The AAA (ATPase associated with diverse cellular activities) ATPase Cdc48 and Dsc5, a ubiquitin regulatory X domain-containing protein, interact with known Dsc complex components and are required for SREBP cleavage. These findings provide a mechanistic link between the Dsc E3 ligase complex and the proteasome in SREBP cleavage and add to a growing list of similarities between the Dsc E3 ligase and membrane E3 ligases involved in endoplasmic reticulum-associated degradation.