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.

Genome-Wide Structural Variation Detection by Genome Mapping on Nanochannel Arrays.

Comprehensive whole-genome structural variation detection is challenging with current approaches. With diploid cells as DNA source and the presence of numerous repetitive elements, short-read DNA sequencing cannot be used to detect structural variation efficiently. In this report, we show that genome mapping with long, fluorescently labeled DNA molecules imaged on nanochannel arrays can be used for whole-genome structural variation detection without sequencing. While whole-genome haplotyping is not achieved, local phasing (across >150-kb regions) is routine, as molecules from the parental chromosomes are examined separately. In one experiment, we generated genome maps from a trio from the 1000 Genomes Project, compared the maps against that derived from the reference human genome, and identified structural variations that are >5 kb in size. We find that these individuals have many more structural variants than those published, including some with the potential of disrupting gene function or regulation.

Assembly and diploid architecture of an individual human genome via single-molecule technologies.

We present the first comprehensive analysis of a diploid human genome that combines single-molecule sequencing with single-molecule genome maps. Our hybrid assembly markedly improves upon the contiguity observed from traditional shotgun sequencing approaches, with scaffold N50 values approaching 30 Mb, and we identified complex structural variants (SVs) missed by other high-throughput approaches. Furthermore, by combining Illumina short-read data with long reads, we phased both single-nucleotide variants and SVs, generating haplotypes with over 99% consistency with previous trio-based studies. Our work shows that it is now possible to integrate single-molecule and high-throughput sequence data to generate de novo assembled genomes that approach reference quality.

Epigenetic regulation of Kaposi's sarcoma-associated herpesvirus latency by virus-encoded microRNAs that target Rta and the cellular Rbl2-DNMT pathway.

Kaposi's sarcoma-associated herpesvirus (KSHV) encodes a cluster of 12 microRNAs (miRNAs) that are processed from a transcript that is embedded within the major latency control region. We have generated a deletion mutation that eliminates 10 of the 12 viral miRNAs from the KSHV bacmid by using recombineering methods. The KSHV miRNA deletion mutant (BAC36 DeltamiR) behaved similarly to wild-type (wt) BAC36 in viral production, latency gene transcription, and viral DNA copy number in 293 and dermal microvascular endothelial cells (DMVECs). However, BAC36 DeltamiR consistently expressed elevated levels of viral lytic genes, including the immediate-early transcriptional activator Rta (ORF50). At least one KSHV microRNA (miRK12-5) was capable of suppressing ORF50 mRNA, but poor seed sequence alignments suggest that these targets may be indirect. Comparison of epigenetic marks in DeltamiR KSHV genomes revealed decreases in histone H3 K9 methylation, increases in histone H3 acetylation, and a striking loss of DNA methylation throughout the viral and cellular genome. One viral miRNA, K12-4-5p, was found to have a sequence targeting retinoblastoma (Rb)-like protein 2 (Rbl2), which is a known repressor of DNA methyl transferase 3a and 3b mRNA transcription. We show that ectopic expression of miR-K12-4-5p reduces Rbl2 protein expression and increases DNMT1, -3a, and -3b mRNA levels relative to the levels for control cells. We conclude that KSHV miRNA targets multiple pathways to maintain the latent state of the KSHV genome, including repression of the viral immediate-early protein Rta and a cellular factor, Rbl2, that regulates global epigenetic reprogramming.

Cohesins localize with CTCF at the KSHV latency control region and at cellular c-myc and H19/Igf2 insulators.

Cohesins, which mediate sister chromatin cohesion, and CTCF, which functions at chromatin boundaries, play key roles in the structural and functional organization of chromosomes. We examined the binding of these two factors on the Kaposi's sarcoma-associated herpesvirus (KSHV) episome during latent infection and found a striking colocalization within the control region of the major latency transcript responsible for expressing LANA (ORF73), vCyclin (ORF72), vFLIP (ORF71), and vmiRNAs. Deletion of the CTCF-binding site from the viral genome disrupted cohesin binding, and crippled colony formation in 293 cells. Clonal instability correlated with elevated expression of lytic cycle gene products, notably the neighbouring promoter for K14 and vGPCR (ORF74). siRNA depletion of RAD21 from latently infected cells caused an increase in K14 and ORF74, and lytic inducers caused a rapid dissociation of RAD21 from the viral genome. RAD21 and SMC1 also associate with the cellular CTCF sites at mammalian c-myc promoter and H19/Igf2 imprinting control region. We conclude that cohesin subunits associate with viral and cellular CTCF sites involved in complex gene regulation and chromatin organization.

Epigenetic control of replication origins.

Efficient duplication of the eukaryotic genome requires the spatial and temporal coordination of numerous replication origins on each chromosome. Epigenetic factors, like chromatin environment, can have profound effects on origin site selection, utilization frequency, and cell cycle firing time. Precisely how chromatin contributes to origin site selection and timing is not completely understood. Recently, we reported on the cell cycle changes in chromatin structure at the plasmid replication origins of Epstein-Barr Virus (EBV) and Kaposi's Sarcoma-Associated Herpes virus (KSHV). These studies and others suggest that cell cycle changes in histone modification and nucleosome remodeling regulate prereplication factor assembly and initiation of DNA replication at origins. We discuss how these studies of viral origins may provide important insights into epigenetic control of cellular chromosome origins.

ORC, MCM, and histone hyperacetylation at the Kaposi's sarcoma-associated herpesvirus latent replication origin.

The viral genome of Kaposi's sarcoma-associated herpesvirus (KSHV) persists as an extrachromosomal plasmid in latently infected cells. The KSHV latency-associated nuclear antigen (LANA) stimulates plasmid maintenance and DNA replication by binding to an approximately 150-bp region within the viral terminal repeats (TR). We have used chromatin immunoprecipitation assays to demonstrate that LANA binds specifically to the replication origin sequence within the KSHV TR in latently infected cells. The latent replication origin within the TR was also bound by LANA-associated proteins CBP, double-bromodomain-containing protein 2 (BRD2), and the origin recognition complex 2 protein (ORC2) and was enriched in hyperacetylated histones H3 and H4 relative to other regions of the latent genome. Cell cycle analysis indicated that the minichromosome maintenance complex protein, MCM3, bound TR in late-G(1)/S-arrested cells, which coincided with the loss of histone H3 K4 methylation. Micrococcal nuclease studies revealed that TRs are embedded in a highly ordered nucleosome array that becomes disorganized in late G(1)/S phase. ORC binding to TR was LANA dependent when reconstituted in transfected plasmids. DNA affinity purification confirmed that LANA, CBP, BRD2, and ORC2 bound TR specifically and identified the histone acetyltransferase HBO1 (histone acetyltransferase binding to ORC1) as a potential TR binding protein. Disruption of ORC2, MCM5, and HBO1 expression by small interfering RNA reduced LANA-dependent DNA replication of TR-containing plasmids. These findings are the first demonstration that cellular replication and origin licensing factors are required for KSHV latent cycle replication. These results also suggest that the KSHV latent origin of replication is a unique chromatin environment containing histone H3 hyperacetylation within heterochromatic tandem repeats.