Massively multiplex single-molecule oligonucleosome footprinting
Bibliographic record
Abstract
These are the intermediate data used in "Massively multiplex single-molecule oligonucleosome footprinting", where the nonspecific adenine methyltransferase EcoGII was used to footprint accessible chromatin, and the methylation was then read using the Pacific Biosciences sequencing platform. The files here are intermediate outputs that capture metrics about the inter-pulse distance values as well as predictions of methylation status. The .npy, .feather, and .pickle files are the output of extractIPD.py, and callNucPeaks.py. The .csv is the output of a cell in SAMOSA_analyses.ipynb. All of this code can be found at https://github.com/RamaniLab/SAMOSA, including SAMOSA_analyses.ipynb which contains all downstream analyses that were performed on this data. meanIPDinfoChrControls.csv: This contains the data used to generate Supplementary Figures 2 and 3. It contains various summary measurements of the IPD values in each molecule of the in vitro samples Files ending in _bingmm.npy: These contain the posterior probability of adenines being methylated for the in vitro data. The files beginning with pbrun3 were sequenced on the Sequel I, and the files beginning with pbrun4 or pbrun5 were sequenced on the Sequel II. Other than Supplementary Figures 2 and 3, the analysis in the paper was based on the Sequel II data. naked_neg and DNA_minusM are both negative controls. naked_methyl and DNA_plusM are positive controls. chromatin samples are in vitro assembled chromatin. pbrun4_gold_nuc47_chromatin_peaks.feather: The estimated nucleosome centers from the in vitro assembled chromatin, in a data frame. Each row is an individual nucleosome dyad prediction. Files ending in _onlyT_zmwinfo.pickle: These files each contain a pandas data frame containing information about each molecule sequenced in the in vivo samples. These should be read in in python using the pandas read_pickle function, and require the same namespace as was used when saving them, so pandas must be imported as pd, and numpy as np. The neg samples are deproteinated unmethylated molecules, the pos samples are deproteinated methylated molecules, and the chromatin samples are methylated chromatin. Files ending in _bingmm.pickle: These files contain the posterior probability of being methylated for each adenine in each DNA molecule in the sample. They each have a corresponding zmwinfo file described above, and similarly must be read in with numpy imported as np. Each file is a dictionary with the zero-mode waveguide (ZMW) hole number as a key, and the value a numpy array with length equal to the unaligned CCS of that molecule, with methylation posterior probabilities at each A/T base. The zmwinfo dataframe has a 'zmw' column that can be used to match the information in that file with the methylation information in this one.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.002 | 0.003 |
| Meta-epidemiology (narrow) | 0.002 | 0.002 |
| Meta-epidemiology (broad) | 0.001 | 0.002 |
| Bibliometrics | 0.002 | 0.002 |
| Science and technology studies | 0.002 | 0.001 |
| Scholarly communication | 0.003 | 0.001 |
| Open science | 0.002 | 0.002 |
| Research integrity | 0.001 | 0.004 |
| Insufficient payload (model declined to judge) | 0.113 | 0.080 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".