DNA-m6A calling and integrated long-read epigenetic and genetic analysis with fibertools.

Long-read DNA sequencing has recently emerged as a powerful tool for studying both genetic and epigenetic architectures at single-molecule and single-nucleotide resolution. Long-read epigenetic studies encompass both the direct identification of native cytosine methylation as well as the identification of exogenously placed DNA N6-methyladenine (DNA-m6A). However, detecting DNA-m6A modifications using single-molecule sequencing, as well as coprocessing single-molecule genetic and epigenetic architectures, is limited by computational demands and a lack of supporting tools. Here, we introduce fibertools, a state-of-the-art toolkit that features a semisupervised convolutional neural network for fast and accurate identification of m6A-marked bases using PacBio single-molecule long-read sequencing, as well as the coprocessing of long-read genetic and epigenetic data produced using either PacBio or Oxford Nanopore sequencing platforms. We demonstrate accurate DNA-m6A identification (>90% precision and recall) along >20 kilobase long DNA molecules with a ~1,000-fold improvement in speed. In addition, we demonstrate that fibertools can readily integrate genetic and epigenetic data at single-molecule resolution, including the seamless conversion between molecular and reference coordinate systems, allowing for accurate genetic and epigenetic analyses of long-read data within structurally and somatically variable genomic regions.

Evidence for Proline Catabolic Enzymes in the Metabolism of Thiazolidine Carboxylates.

Thiazolidine carboxylates such as thiazolidine-4-carboxylate (T4C) and thiazolidine-2-carboxylate (T2C) are naturally occurring sulfur analogues of proline. These compounds have been observed to have both beneficial and toxic effects in cells. Given that proline dehydrogenase has been proposed to be a key enzyme in the oxidative metabolism of thioprolines, we characterized T4C and T2C as substrates of proline catabolic enzymes using proline utilization A (PutA), which is a bifunctional enzyme with proline dehydrogenase (PRODH) and l-glutamate-γ-semialdehyde dehydrogenase (GSALDH) activities. PutA is shown here to catalyze the FAD-dependent PRODH oxidation of both T4C and T2C with catalytic efficiencies significantly higher than with proline. Stopped-flow experiments also demonstrate that l-T4C and l-T2C reduce PutA-bound FAD at rates faster than proline. Unlike proline, however, oxidation of T4C and T2C does not generate a substrate for NAD+-dependent GSALDH. Instead, PutA/PRODH oxidation of T4C leads to cysteine formation, whereas oxidation of T2C generates an apparently stable Δ4-thiazoline-2-carboxylate species. Our results provide new insights into the metabolism of T2C and T4C.

Structural analysis of prolines and hydroxyprolines binding to the l-glutamate-γ-semialdehyde dehydrogenase active site of bifunctional proline utilization A.

Proline utilization A (PutA) proteins are bifunctional proline catabolic enzymes that catalyze the 4-electron oxidation of l-proline to l-glutamate using spatially-separated proline dehydrogenase and l-glutamate-γ-semialdehyde dehydrogenase (GSALDH, a.k.a. ALDH4A1) active sites. The observation that l-proline inhibits both the GSALDH activity of PutA and monofunctional GSALDHs motivated us to study the inhibition of PutA by proline stereoisomers and analogs. Here we report five high-resolution crystal structures of PutA with the following ligands bound in the GSALDH active site: d-proline, trans-4-hydroxy-d-proline, cis-4-hydroxy-d-proline, l-proline, and trans-4-hydroxy-l-proline. Three of the structures are of ternary complexes of the enzyme with an inhibitor and either NAD+ or NADH. To our knowledge, the NADH complex is the first for any GSALDH. The structures reveal a conserved mode of recognition of the inhibitor carboxylate, which results in the pyrrolidine rings of the d- and l-isomers having different orientations and different hydrogen bonding environments. Activity assays show that the compounds are weak inhibitors with millimolar inhibition constants. Curiously, although the inhibitors occupy the aldehyde binding site, kinetic measurements show the inhibition is uncompetitive. Uncompetitive inhibition may involve proline binding to a remote site or to the enzyme-NADH complex. Together, the structural and kinetic data expand our understanding of how proline-like molecules interact with GSALDH, reveal insight into the relationship between stereochemistry and inhibitor affinity, and demonstrate the pitfalls of inferring the mechanism of inhibition from crystal structures alone.

Resolving the chromatin impact of mosaic variants with targeted Fiber-seq.

Accurately quantifying the functional consequences of non-coding mosaic variants requires the pairing of DNA sequence with both accessible and closed chromatin architectures along individual DNA molecules-a pairing that cannot be achieved using traditional fragmentation-based chromatin assays. We demonstrate that targeted single-molecule chromatin fiber sequencing (Fiber-seq) achieves this, permitting single-molecule, long-read genomic and epigenomic profiling across targeted >100 kilobase loci with ~10-fold enrichment over untargeted sequencing. Targeted Fiber-seq reveals that pathogenic expansions of the DMPK CTG repeat that underlie Myotonic Dystrophy 1 are characterized by somatic instability and disruption of multiple nearby regulatory elements, both of which are repeat length-dependent. Furthermore, we reveal that therapeutic adenine base editing of the segmentally duplicated γ-globin ( HBG1/HBG2 ) promoters in primary human hematopoietic cells induced towards an erythroblast lineage increases the accessibility of the HBG1 promoter as well as neighboring regulatory elements. Overall, we find that these non-protein coding mosaic variants can have complex impacts on chromatin architectures, including extending beyond the regulatory element harboring the variant.

DNA-m6A calling and integrated long-read epigenetic and genetic analysis with fibertools.

Long-read DNA sequencing has recently emerged as a powerful tool for studying both genetic and epigenetic architectures at single-molecule and single-nucleotide resolution. Long-read epigenetic studies encompass both the direct identification of native cytosine methylation as well as the identification of exogenously placed DNA N6-methyladenine (DNA-m6A). However, detecting DNA-m6A modifications using single-molecule sequencing, as well as co-processing single-molecule genetic and epigenetic architectures, is limited by computational demands and a lack of supporting tools. Here, we introduce fibertools, a state-of-the-art toolkit that features a semi-supervised convolutional neural network for fast and accurate identification of m6A-marked bases using PacBio single-molecule long-read sequencing, as well as the co-processing of long-read genetic and epigenetic data produced using either PacBio or Oxford Nanopore sequencing platforms. We demonstrate accurate DNA-m6A identification (>90% precision and recall) along >20 kilobase long DNA molecules with a ~1,000-fold improvement in speed. In addition, we demonstrate that fibertools can readily integrate genetic and epigenetic data at single-molecule resolution, including the seamless conversion between molecular and reference coordinate systems, allowing for accurate genetic and epigenetic analyses of long-read data within structurally and somatically variable genomic regions.

Synchronized long-read genome, methylome, epigenome, and transcriptome for resolving a Mendelian condition.

Resolving the molecular basis of a Mendelian condition (MC) remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome, and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion, and structural variant calling and diploid de novo genome assembly, and permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility, and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network (UDN) participant with a chromosome X;13 balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1, and RB1) previously associated with single-gene MCs. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included nonsense-mediated decay, fusion transcript formation, enhancer adoption, transcriptional readthrough silencing, and inappropriate X chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.