BACKGROUND: Cell free DNA (cfDNA)-based assays hold great potential in detecting early cancer signals yet determining the tissue-of-origin (TOO) for cancer signals remains a challenging task. Here, we investigated the contribution of a methylation atlas to TOO detection in low depth cfDNA samples.
METHODS: We constructed a tumor-specific methylation atlas (TSMA) using whole-genome bisulfite sequencing (WGBS) data from five types of tumor tissues (breast, colorectal, gastric, liver and lung cancer) and paired white blood cells (WBC). TSMA was used with a non-negative least square matrix factorization (NNLS) deconvolution algorithm to identify the abundance of tumor tissue types in a WGBS sample. We showed that TSMA worked well with tumor tissue but struggled with cfDNA samples due to the overwhelming amount of WBC-derived DNA. To construct a model for TOO, we adopted the multi-modal strategy and used as inputs the combination of deconvolution scores from TSMA with other features of cfDNA.
RESULTS: Our final model comprised of a graph convolutional neural network using deconvolution scores and genome-wide methylation density features, which achieved an accuracy of 69% in a held-out validation dataset of 239 low-depth cfDNA samples.
CONCLUSIONS: In conclusion, we have demonstrated that our TSMA in combination with other cfDNA features can improve TOO detection in low-depth cfDNA samples.

DNA polymerase κ (Polκ) is a specialized polymerase that has multiple cellular roles such as translesion DNA synthesis, replication of repetitive sequences, and nucleotide excision repair. We have developed a method for capturing DNA synthesized by Polκ utilizing a Polκ-specific substrate, N2-(4-ethynylbenzyl)-2\'-deoxyguanosine (EBndG). After shearing of the DNA into 200 to 500 bp lengths, the EBndG-containing DNA was covalently bound to biotin using the Cu(I)-catalyzed alkyne-azide cycloaddition reaction and isolated with streptavidin beads. Isolated DNA was then ligated to adaptors, followed by PCR amplification and next-generation sequencing to generate genome-wide repair maps. We have termed this method polymerase κ sequencing. Here, we present the human genome maps for Polκ activity in an undamaged cell line. We found that Polκ activity was enhanced in GC-rich regions, euchromatin regions, the promoter of genes, and in DNA that is replicated early in the S phase.

Cycling cells replicate their DNA during the S phase through a defined temporal program known as replication timing. Mutation frequencies, epigenetic chromatin states, and transcriptional activities are different for genomic regions that are replicated early and late in the S phase. Here, we found from ChIP-Seq analysis that DNA polymerase (Pol) κ is enriched in early-replicating genomic regions in HEK293T cells. In addition, by feeding cells with N 2-heptynyl-2\'-deoxyguanosine followed by click chemistry-based enrichment and high-throughput sequencing, we observed elevated Pol κ activities in genomic regions that are replicated early in the S phase. On the basis of the established functions of Pol κ in accurate and efficient nucleotide insertion opposite endogenously induced N 2-modified dG lesions, our work suggests that active engagement of Pol κ may contribute to diminished mutation rates observed in early-replicating regions of the human genome, including cancer genomes. Together, our work expands the functions of Pol κ and offered a plausible mechanism underlying replication timing-dependent mutation accrual in the human genome.

The clinical success of CRISPR therapies hinges on the safety and efficacy of Cas proteins. The Cas9 from Francisella novicida (FnCas9) is highly precise, with a negligible affinity for mismatched substrates, but its low cellular targeting efficiency limits therapeutic use. Here, we rationally engineer the protein to develop enhanced FnCas9 (enFnCas9) variants and broaden their accessibility across human genomic sites by ~3.5-fold. The enFnCas9 proteins with single mismatch specificity expanded the target range of FnCas9-based CRISPR diagnostics to detect the pathogenic DNA signatures. They outperform Streptococcus pyogenes Cas9 (SpCas9) and its engineered derivatives in on-target editing efficiency, knock-in rates, and off-target specificity. enFnCas9 can be combined with extended gRNAs for robust base editing at sites which are inaccessible to PAM-constrained canonical base editors. Finally, we demonstrate an RPE65 mutation correction in a Leber congenital amaurosis 2 (LCA2) patient-specific iPSC line using enFnCas9 adenine base editor, highlighting its therapeutic utility.

Gene regulation is essential to placental function and fetal development. We built a genome-scale transcriptional regulatory network (TRN) of the human placenta using digital genomic footprinting and transcriptomic data. We integrated 475 transcriptomes and 12 DNase hypersensitivity datasets from placental samples to globally and quantitatively map transcription factor (TF)-target gene interactions. In an independent dataset, the TRN model predicted target gene expression with an out-of-sample R2 greater than 0.25 for 73% of target genes. We performed siRNA knockdowns of four TFs and achieved concordance between the predicted gene targets in our TRN and differences in expression of knockdowns with an accuracy of >0.7 for three of the four TFs. Our final model contained 113,158 interactions across 391 TFs and 7712 target genes and is publicly available. We identified 29 TFs which were significantly enriched as regulators for genes previously associated with preterm birth, and eight of these TFs were decreased in preterm placentas.

CONCLUSIONS: Improvements in nanopore sequencing necessitate efficient classification methods, including pre-filtering and adaptive sampling algorithms that enrich for reads of interest. Signal-based approaches circumvent the computational bottleneck of basecalling. But past methods for signal-based classification do not scale efficiently to large, repetitive references like pangenomes, limiting their utility to partial references or individual genomes. We introduce Sigmoni: a rapid, multiclass classification method based on the r-index that scales to references of hundreds of Gbps. Sigmoni quantizes nanopore signal into a discrete alphabet of picoamp ranges. It performs rapid, approximate matching using matching statistics, classifying reads based on distributions of picoamp matching statistics and co-linearity statistics, all in linear query time without the need for seed-chain-extend. Sigmoni is 10-100× faster than previous methods for adaptive sampling in host depletion experiments with improved accuracy, and can query reads against large microbial or human pangenomes. Sigmoni is the first signal-based tool to scale to a complete human genome and pangenome while remaining fast enough for adaptive sampling applications.
METHODS: Sigmoni is implemented in Python, and is available open-source at https://github.com/vshiv18/sigmoni.

Endogenous retroviruses (ERVs) are related to long terminal repeat (LTR) retrotransposons, comprising gene sequences of exogenous retroviruses integrated into the host genome and inherited according to Mendelian law. They are considered to have contributed greatly to the evolution of host genome structure and function. We previously characterized HERV-K HML-9 in the human genome. However, the biological function of this type of element in the genome of the chimpanzee, which is the closest living relative of humans, largely remains elusive. Therefore, the current study aims to characterize HML-9 in the chimpanzee genome and to compare the results with those in the human genome. Firstly, we report the distribution and genetic structural characterization of the 26 proviral elements and 38 solo LTR elements of HML-9 in the chimpanzee genome. The results showed that the distribution of these elements displayed a non-random integration pattern, and only six elements maintained a relatively complete structure. Then, we analyze their phylogeny and reveal that the identified elements all cluster together with HML-9 references and with those identified in the human genome. The HML-9 integration time was estimated based on the 2-LTR approach, and the results showed that HML-9 elements were integrated into the chimpanzee genome between 14 and 36 million years ago and into the human genome between 18 and 49 mya. In addition, conserved motifs, cis-regulatory regions, and enriched PBS sequence features in the chimpanzee genome were predicted based on bioinformatics. The results show that pathways significantly enriched for ERV LTR-regulated genes found in the chimpanzee genome are closely associated with disease development, including neurological and neurodevelopmental psychiatric disorders. In summary, the identification, characterization, and genomics of HML-9 presented here not only contribute to our understanding of the role of ERVs in primate evolution but also to our understanding of their biofunctional significance.

Deficiencies in DNA mismatch repair (MMRd) leave characteristic footprints of microsatellite instability (MSI) in cancer genomes. We used data from the Cancer Genome Atlas and International Cancer Genome Consortium to conduct a comprehensive analysis of MSI-associated cancers, focusing on indel mutational signatures. We classified MSI-high genomes into two subtypes based on their indel profiles: deletion-dominant (MMRd-del) and insertion-dominant (MMRd-ins). Compared with MMRd-del genomes, MMRd-ins genomes exhibit distinct mutational and transcriptomic features, including a higher prevalence of T>C substitutions and related mutation signatures. Short insertions and deletions in MMRd-ins and MMRd-del genomes target different sets of genes, resulting in distinct indel profiles between the two subtypes. In addition, indels in the MMRd-ins genomes are enriched with subclonal alterations that provide clues about a distinct evolutionary relationship between the MMRd-ins and MMRd-del genomes. Notably, the transcriptome analysis indicated that MMRd-ins cancers upregulate immune-related genes, show a high level of immune cell infiltration, and display an elevated neoantigen burden. The genomic and transcriptomic distinctions between the two types of MMRd genomes highlight the heterogeneity of genetic mechanisms and resulting genomic footprints and transcriptomic changes in cancers, which has potential clinical implications.

The Merovingian period (5th to 8th cc AD) was a time of demographic, socioeconomic, cultural, and political realignment in Western Europe. Here, we report the whole-genome shotgun sequence data of 30 human skeletal remains from a coastal Late Merovingian site of Koksijde (675 to 750 AD), alongside 18 remains from two Early to Late Medieval sites in present-day Flanders, Belgium. We find two distinct ancestries, one shared with Early Medieval England and the Netherlands, while the other, minor component, reflecting likely continental Gaulish ancestry. Kinship analyses identified no large pedigrees characteristic to elite burials revealing instead a high modularity of distant relationships among individuals of the main ancestry group. In contrast, individuals with >90% Gaulish ancestry had no kinship links among sampled individuals. Evidence for population structure and major differences in the extent of Gaulish ancestry in the main group, including in a mother-daughter pair, suggests ongoing admixture in the community at the time of their burial. The isotopic and genetic evidence combined supports a model by which the burials, representing an established coastal nonelite community, had incorporated migrants from inland populations. The main group of burials at Koksijde shows an abundance of >5 cM long shared allelic intervals with the High Medieval site nearby, implying long-term continuity and suggesting that similarly to Britain, the Early Medieval ancestry shifts left a significant and long-lasting impact on the genetic makeup of the Flemish population. We find substantial allele frequency differences between the two ancestry groups in pigmentation and diet-associated variants, including those linked with lactase persistence, likely reflecting ancestry change rather than local adaptation.

Hodgkin lymphoma (HL) is one of the most common lymphomas, with an incidence of 3 per 100,000 persons. Current treatment uses a cocktail of genotoxic agents, including adriamycin, bleomycin, vinblastine, and dacarbazine (ABVD), along with or without radiotherapy. This treatment regimen has proved to be efficient in killing cancer cells, resulting in HL patients having a survival rate of >90% cancer-free survival at five years. However, this therapy does not have a specific cell target, and it can induce damage in the genome of non-cancerous cells. Previous studies have shown that HL survivors often exhibit karyotypes characterized by complex chromosomal abnormalities that are difficult to analyze by conventional banding. Multicolor fluorescence in situ hybridization (M-FISH) is a powerful tool to analyze complex karyotypes; we used M-FISH to investigate the presence of chromosomal damage in peripheral blood lymphocytes from five healthy individuals and five HL patients before, during, and one year after anti-cancer treatment. Our results show that this anti-cancer treatment-induced genomic chaos that persists in the hematopoietic stem cells from HL patients one year after finishing therapy. This chromosomal instability may play a role in the occurrence of second primary cancers that are observed in 10% of HL survivors. This chapter will describe a protocol for utilizing M-FISH to study treatment-induced genome chaos in Hodgkin\'s lymphoma (HL) patients, following a brief discussion.