Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Identification and functional relevance of de novo DNA methylation in cancerous B-cell populations.

Epigenetic remodeling is a hallmark of cancer, with the frequent acquisition of de novo DNA methylation in CpG islands. However, the functional relevance of de novo DNA methylation in cancer is less well defined. To begin to address this issue in B-cells, we used BeadArray assays to survey the methylation status of over 1,500 cancer-related CpG loci in two molecular subtypes of diffuse large B-cell lymphoma (ABC-DLBCL and GCB-DLBCL) and cognate normal B-cell populations. We identified 81 loci that showed frequent de novo DNA methylation in GCB-DLBCL and 67 loci that showed frequent de novo DNA methylation in ABC-DLBCL. These de novo methylated CpG loci included reported targets of polycomb repressive complexes (PRC) in stem cells. All candidate loci in GCB-DLBCL are proximal to genes that are poorly expressed or silent in purified normal germinal center (GC) B-cells. This is consistent with the hypothesis that de novo DNA methylation in cancer is more frequently involved in the maintenance rather than the initiation of gene silencing (de novo repression). This suggests that epigenetic switching occurs during tumorigenesis with de novo DNA methylation locking in gene silencing normally mediated by transcriptional repressors. Furthermore, we propose that similar to de novo genetic mutations, the majority of de novo DNA methylation events observed in tumors are passengers not causally involved in tumorigenesis.

Large-scale profiling of archival lymph nodes reveals pervasive remodeling of the follicular lymphoma methylome.

Emerging technologies allow broad profiling of the cancer genome for differential DNA methylation relative to benign cells. Herein, bisulfite-modified DNA from lymph nodes with either reactive hyperplasia or follicular lymphoma (FL) were analyzed using a commercial C/UpG genotyping assay. Two hundred fifty-nine differentially methylated targets (DMT) distributed among 183 unique genes were identified in FL. Comparison of matched formalin-fixed, paraffin-embedded and frozen surgical pathology replicates showed the complete preservation of the cancer methylome among differently archived tissue specimens. Analysis of the DMT profile is consistent with a pervasive epigenomic remodeling process in FL that affects predominantly nonlymphoid genes.

Inhibition of prostate cancer metastatic colonization by approximately 4.2 Mb of human chromosome 12.

Our previous studies demonstrate that introduction of a approximately 70 cM region (now estimated at 63.75 Mb by the Human Genome Project) of human chromosome 12 into the highly metastatic Dunning rat prostate cancer cell line AT6.1 results in >30-fold (>/=90%) reduction in the number of overt metastases in spontaneous metastasis assays. We report the further localization and biological characterization of the metastasis-suppressor activity encoded by a reduced region of chromosome 12. To localize this metastasis-suppressor activity, a panel of AT6.1 microcell hybrids that retain varying portions of human chromosome 12 was constructed and subjected to sequence-tagged site (STS)-based PCR analysis and assessment of in vivo metastatic ability. Data from these complementary approaches localized the metastasis-suppressor activity to a approximately 4.2 Mb portion of human chromosome 12q24.3 comprised of 3 separate regions. Reverse transcriptase-polymerase chain reaction (RT-PCR) and immunoblotting were used for differential expression analyses to identify which characterized genes, predicted gene sequences and expressed sequence tags (ESTs) within this region could be responsible for the observed metastasis suppression. Comprehensive in vivo studies showed that suppressed AT6.1-12 hybrids that retain the metastasis-suppressor region on 12q24.3 are capable of arriving at the secondary site, but are not able to persist there. Thus, unlike other metastasis-suppressor genes characterized to date, the metastasis-suppressor gene encoded by this region appears to utilize a different biologic mechanism to suppress the growth of overt metastases at the secondary site.