SEARCH · PubMed Health

Results for “high-throughput functional genomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

PubMed2026

A de novo algorithm for allele reconstruction from Oxford nanopore amplicon reads, with application to CYP2D6.

MOTIVATION: The Oxford Nanopore Technologies' sequencing platform offers a path towards bedside genomics, producing long reads that can completely cover a gene of interest, and detect any known or novel variant the gene contains. However, the analysis of these long reads to identify actionable genotypes remains challenging and typically requires customization depending on the target gene. RESULTS: Here, we describe a generic algorithm to accurately reconstruct allele sequences derived from long-reads of amplicon-based data. Rather than calling variants directly from these long-reads, our method takes a "sequence-first" approach, performing an unbiased reconstruction of the underlying amplicon sequences to generate high-confidence reconstructed allele sequences. This is done without user input of the target gene, allowing for any source amplicon to be reconstructed. These high-confidence reconstructed allele sequences are then compared to the genomic reference sequence of the gene to infer the specific diplotype present in the sample. This approach is agnostic towards the number of genes and alleles present and readily detects novel variants. We demonstrate our approach using three independent data sets for CYP2D6, a diverse and complex gene with over 175 known alleles of clinical significance. We show how our approach can accurately recover validated CYP2D6 diplotypes from 20 Coriell samples covering 14 distinct alleles, using different amplicons, flow cell versions, and depths. This includes inferring occurrences of allele duplication events from relative abundances of each allele, a critical factor for ascribing functional effects to a diplotype. Further, we demonstrate our approach's utility for other genomic regions, including HLA. AVAILABILITY: Custom code is available at the following GitHub repository, along with instructions for use and test data: https://github.com/scottdbrown/allele-reconstruction-long-read-amplicon-data. A snapshot of the code at the time of publication is available on Zenodo.org; doi 10.5281/zenodo.19716004. Raw .fastq sequence data for our three sequencing runs is available at the SRA under Bioproject PRJNA1357883 (https://www.ncbi.nlm.nih.gov/bioproject/1357883).

Alleles↗

PubMed2026

Functional Prediction of Epitranscriptome.

N6-methyladenosine (m6A) is one of the most prevalent and well-studied RNA modifications, playing a pivotal role in many biological processes. With the recent advances in high-throughput sequencing technologies, tens of thousands of m6A sites have been reported. However, not all m6A sites are important or functionally significant, highlighting the need to distinguish biologically relevant m6As from non-functional or technically artefactual ones. Here, we describe ConsRM, which is a web-based resource that was designed to evaluate the importance of m6As from an evolutionary perspective. It introduced a novel scoring framework for quantifying the conservation degree of m6As in humans. Its web interface includes a database of 177998 distinct human m6A sites along with their calculated conservation score, and allows users to analyze their own data via the web server. ConsRM is freely accessible at: http://180.208.58.19/conservation/browser.html .

Humans↗

PubMed2026

A novel peptide encoded by circTLL1 drives osimertinib resistance in lung cancer by modulating the NT5C2/Ras/PI3K axis.

BACKGROUND: Acquired resistance to osimertinib, a third-generation EGFR tyrosine kinase inhibitor, remains a major clinical challenge in the treatment of non-small cell lung cancer (NSCLC). Although circular RNAs (circRNAs) have been increasingly implicated in drug resistance, most studies have focused on their canonical role as microRNA sponges, while their capacity to encode functional micropeptides remains largely unexplored. This study aimed to identify novel circRNAs involved in osimertinib resistance and to characterize their regulatory functions at the protein level. METHODS: Osimertinib-resistant (OR) NSCLC cell lines were established and validated. High-throughput RNA sequencing was performed to compare the circRNA expression profiles between parental and OR cells. The function of the candidate circRNA was assessed through a series of in vitro and in vivo experiments, including cell viability assays, apoptosis analysis, and xenograft mouse models. Mechanistic investigations involved mass spectrometry, co-immunoprecipitation and western blotting to explore its protein-coding potential and downstream signaling pathways. RESULTS: We identified a novel circRNA, termed circTLL1, that was stably and significantly upregulated in OR-NSCLC cells. Functionally, overexpression of circTLL1 promoted osimertinib resistance, whereas its knockdown restored drug sensitivity both in vitro and in vivo. Mechanistically, we discovered that circTLL1 harbors an open reading frame (ORF) that is translated into a novel 90-amino-acid protein, which we designated circTLL1-90aa. Further investigation revealed that circTLL1-90aa directly interacts with and promotes the degradation of 5'-nucleotidase, cytosolic II (NT5C2), thereby uncoupling nucleotide metabolism from its normal regulatory constraints. The consequent downregulation of NT5C2 leads to elevated GTP levels and leading to the sustained activation of the downstream Ras/PI3K/AKT signaling pathway. CONCLUSION: Our findings unveil a previously unrecognized circRNA/micropeptide/metabolism cascade underlying osimertinib resistance. The identification of the circTLL1-90aa/NT5C2/Ras/PI3K axis not only expands the functional repertoire of the non-coding genome but also provides new insights into the complexity of drug resistance. Given its selective upregulation in resistant cells, circTLL1-90aa holds promise both as a predictive biomarker for treatment stratification and as an actionable therapeutic target, offering a novel strategy to overcome osimertinib resistance in NSCLC patients.

Pyrimidines↗

PubMed2026

esloco: simulation-based estimation of local coverage in long-read DNA sequencing.

SUMMARY: Long-read DNA sequencing is increasingly applied for whole-genome studies, yet experimental planning often lacks reliable estimates of target region coverage, leading to costly and time-consuming pilot studies and replicates. We present esloco, a Monte Carlo-based simulation framework for estimating local coverage in long-read sequencing experiments, including scenarios with unknown target regions (e.g. viral integration, CRISPR-Cas9) or PCR-free designs (e.g. base modifications). By modeling coverage as a function of sequencing depth and read length distribution, esloco enables informed predictions of local sequencing outcomes. Benchmarking across a 45-gene panel demonstrated close agreement with empirical data, underscoring the framework's reliability. AVAILABILITY AND IMPLEMENTATION: esloco is a Python package available on PyPI (https://pypi.org/project/esloco/), GitHub (https://github.com/aweich/esloco), and Zenodo (https://doi.org/10.5281/zenodo.17776161).

Sequence Analysis, DNA↗

PubMed2025

Decipher RNA isoform combinations from minigene splicing assays and massive parallel sequencing with MAGIC.

SUMMARY: Functional testing of RNA using minigene splicing assays is increasingly being realized to demonstrate the effects of variants on splicing. In complex cases, variant pathogenicity is assessed by Sanger sequencing, which can be time consuming and may be replaced by short read sequencing. Moreover, strategies based on long read sequencing of the amplified minigene construct are promising and allow the isoforms to be fully characterized. We introduce MAGIC, a user-friendly tool that first generates the artificial construction genome files required to then perform alignment, assembly and annotation of the isoforms obtained by either short or long read minigene splicing assay sequencing. AVAILABILITY AND IMPLEMENTATION: MAGIC is available at https://github.com/LBGC-CFB/MAGIC. Zenodo DOI: 10.5281/zenodo.17052752.

High-Throughput Nucleotide Sequencing↗

PubMed2025

DNA Methylation Analysis by Bisulfite Pyrosequencing of Mouse Embryonic Fibroblasts with Reprogramming Enhanced by Thyroid Hormones.

DNA methylation is a widely studied epigenetic mark which in mammals involves the incorporation of a methyl group to the fifth carbon of cytosines, mainly those belonging to CpG dinucleotides. It has been linked to context-dependent regulatory functions ranging from gene and repetitive DNA silencing to gene body transcriptional activity. Because of its important roles during embryonic development and cell differentiation, DNA methylation can be used to track cell reprogramming by measuring the methylation levels of pluripotency-associated factors. In this scenario, bisulfite pyrosequencing is a simple, robust, and widely used technique which allows for the quantification of DNA methylation levels at small, specific regions of the genome. It involves the amplification and biotin tagging of bisulfite-converted DNA. Single amplified strands are then purified using streptavidin and finally pyrosequenced using a sequencing primer. Thus, it is an ideal method for the quantitative profiling of specific genomic regions, with applications ranging from biomarker discovery and epigenetic clock tracking to omic validation studies.

Animals↗

PubMed2025

Evaluation of sequencing reads at scale using rdeval.

MOTIVATION: Large sequencing datasets are being produced and deposited into public archives at unprecedented rates. The availability of tools that can reliably and efficiently generate and store sequencing read summary statistics has become critical. RESULTS: As part of the effort by the Vertebrate Genomes Project (VGP) to generate high-quality reference genomes at scale, we sought to address the community's need for efficient sequence data evaluation by developing rdeval, a standalone tool to quickly compute and interactively display sequencing read metrics. Rdeval can either run on the fly or store key sequence data metrics in tiny read 'snapshot' files. Statistics can then be efficiently recalled from snapshots for additional processing. Rdeval can convert fa*[.gz] files to and from other popular formats including BAM and CRAM for better compression. Overall, while CRAM achieves the best compression, the gain compared to BAM is marginal, and BAM achieves the best compromise between data compression and access speed. Rdeval also generates a detailed visual report with multiple data analytics that can be exported in various formats. We showcase rdeval's functionalities using long-read data from different sequencing platforms and species, including human. For PacBio long-read sequencing, our analysis shows dramatic improvements in both read length and quality over time, as well as the benefit of increased coverage for genome assembly, though the magnitude varies by taxa. AVAILABILITY AND IMPLEMENTATION: Rdeval is implemented in C++ for data processing and in R for data visualization. Precompiled releases (Linux, MacOS, Windows) and commented source code for rdeval are available under MIT license at https://github.com/vgl-hub/rdeval. Documentation is available on ReadTheDocs (https://rdeval-documentation.readthedocs.io). Rdeval is also available in Bioconda and in Galaxy (https://usegalaxy.org). An automated test workflow ensures the consistency of software updates.

Software↗

PubMed2025

Genome-Wide Profiling of Histone Modifications in Fission Yeast Using CUT&Tag.

Eukaryotic DNA is organized in the nucleus in the form of chromatin. Nucleosomes, the fundamental unit of chromatin, are subject to many posttranslational modifications (PTMs) as well as compositional variations through incorporation of histone variants. These alterations play important roles in regulation of genome structure and activity. Genome-wide profiling of these regulatory features is essential for understanding of genome function. Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-Seq) is a widely used method to assay genome-wide localization in fission yeast but suffers from the requirement for a large amount of input chromatin, antibodies, and a cumbersome experimental pipeline. New methods such as Cleavage Under Targets and Tagmentation (CUT&Tag), which combine the specificity of targeted cleavage and adapter insertion with the sensitivity of next-generation sequencing, enable identification and characterization of various epigenetic marks affording low input requirement as well as more streamlined protocols. However, these approaches have not been adapted for use in fission yeast, Schizosaccharomyces pombe. Here, we describe an adapted CUT&Tag protocol for epigenomic profiling in fission yeast using the heterochromatin-associated histone H3K9 methylation PTM for benchmarking.

Schizosaccharomyces↗

PubMed2024

Multiome Perturb-seq unlocks scalable discovery of integrated perturbation effects on the transcriptome and epigenome.

Single-cell CRISPR screens link genetic perturbations to transcriptional states, but high-throughput methods connecting these induced changes to their regulatory foundations are limited. Here, we introduce Multiome Perturb-seq, extending single-cell CRISPR screens to simultaneously measure perturbation-induced changes in gene expression and chromatin accessibility. We apply Multiome Perturb-seq in a CRISPRi screen of 13 chromatin remodelers in human RPE-1 cells, achieving efficient assignment of sgRNA identities to single nuclei via an improved method for capturing barcode transcripts from nuclear RNA. We organize expression and accessibility measurements into coherent programs describing the integrated effects of perturbations on cell state, finding that ARID1A and SUZ12 knockdowns induce programs enriched for developmental features. Modeling of perturbation-induced heterogeneity connects accessibility changes to changes in gene expression, highlighting the value of multimodal profiling. Overall, our method provides a scalable and simply implemented system to dissect the regulatory logic underpinning cell state. A record of this paper's transparent peer review process is included in the supplemental information.

Humans↗

PubMed2026

Utilization of long-read sequencing for the detection of structural rearrangements with AgileStructure.

MOTIVATION: Changes in genome organisation contribute to genetic disease when they disrupt gene function or regulation. Structural rearrangements may interrupt coding sequence or alter expression through promoter loss or gain, chromatin changes, copy-number variation, or disruption of short-range regulatory elements. Although short-read sequencing excels at detecting small variants, it performs poorly at resolving breakpoints of large rearrangements, especially in repetitive or low-complexity regions. Long-read sequencing overcomes these limitations, but analytical tools have not kept pace, making accurate identification and annotation of large structural variants challenging. RESULTS: We developed AgileStructure, a desktop application for locating and annotating large‑scale genomic rearrangements using aligned long‑read data. The software enables user‑guided exploration of breakpoint‑spanning reads, supporting accurate interpretation of complex events and filling a key gap in current structural variant analysis workflows. AVAILABILITY AND IMPLEMENTATION: Source code, binaries, user guide, and example aligned read data, are available on GitHub: https://github.com/msjimc/AgileStructure. An archived version is also available on Zenodo at https://doi.org/10.5281/zenodo.18610110.

Software↗

PubMed2026

Comparison of paralog identification methods and their impact on species tree topologies in target capture phylogenomics within the Sindora clade (Detarioideae: Leguminosae).

Target capture is a common method of generating high throughput DNA sequencing data for phylogenetic reconstruction of species relationships, for which single copy genes are usually most informative. However, a pervasive problem with target capture is that putatively single copy genes may in fact be paralogs resulting from gene duplication, which are problematic for phylogenetic inference because their evolutionary history may differ from the divergence history of species. Here, we use as a case study a target enrichment dataset of 88 species of Detarioideae (Leguminosae) with a focus on the Sindora clade to examine approaches for handling paralogs, including the built-in paralog handling functions in HybPiper and CAPTUS, plus subsequent steps using Putative Paralog Detection and the tree-based Yang & Smith orthology inference approach. We compare the paralogs flagged using these methods and verify their performance with BLAST mapping against a reference genome sequence of Sindora glabra, and then subsequently compare the species tree topologies produced across these methods. Our comparisons of paralogs flagged across the Sindora clade show that the Putative Paralog Detection pipeline was the most accurate in identifying paralogs in terms of its similarity to the BLAST mapping, followed by the built-in paralog identification function of CAPTUS. However, the results we recovered for the Detarioideae subfamily suggest that the largest differences in species tree topology resulted from the use of paralog-filtered alignments (such as with the Putative Paralog Detection pipeline and the Yang & Smith orthology inference approaches) rather than just by removing the sequences of identified paralogous genes. This was the true for HybPiper-assembled datasets but was not seen in CAPTUS-assembled datasets. In all comparisons, the topological differences caused by different paralog handling methods tended to be confined to clades where processes such as hybridisation and introgression are prevalent. Our study provides a roadmap to establish the best approach to identify, eliminate or separate paralogs in the absence of a chromosomally contiguous reference genome for a study group, and highlights the importance of careful data inspection and processing in addition to understanding the extent of paralogy and paralog characteristics (e.g. sequence divergence between copies) for their study group.

Phylogeny↗

PubMed2026

Obtaining a Diagnostic Yield via Scan findings prior to the introduction of SEquencing retrospectivelY (ODYSSEY): a cohort study.

OBJECTIVE: To determine the retrospective yield of prenatal exome sequencing (PES) by establishing the proportion of children with a postnatal monogenic diagnosis that could have been diagnosed prenatally if PES had been available. METHODS: The study cohort comprised a sample of children in Northern Ireland, born between January 2010 and January 2018 (predating routine availability of PES), who received a monogenic diagnosis postnatally via next generation sequencing as part of either of two UK-wide studies (the 100 000 Genomes Project (2015-2018) or the Deciphering Developmental Disorders study (2011-2015)). Clinical data were collected retrospectively and correlated with the current UK National Health Service PES protocol, including the phenotypic eligibility criteria for PES and the associated fetal anomalies gene panel. Cases were considered retrospective diagnoses if the fetal phenotype would have been eligible for PES and the diagnostic gene was included on the test panel, meaning prenatal diagnosis in this current era could have been feasible. RESULTS: Of 101 children, 17.8% (95% CI, 10.3-25.3%) had both an eligible fetal structural anomaly (FSA) (i.e. high-risk FSA) and a diagnostic gene on the associated test panel, meaning that they could have been diagnosed prenatally in the current clinical landscape. The median length of the diagnostic odyssey for this subgroup of children was 3.7 years (1354 (range, 822-2450) days). Moreover, 58.4% (n = 59) of cases had no anomalies detected prenatally and 19.8% (n = 20) had a FSA that would not meet the eligibility criteria for PES (low-risk FSA). Although these cases would have been ineligible for PES under the current clinical pathway, 89.9% (n = 71/79) were affected by severe or profound syndromes. Postnatally, the most common functional anomalies were neurodevelopmental delay/intellectual disability and/or behavioral abnormality, which were observed in 80.2% (n = 81) of the included children. However, 80.2% (n = 65/81) of these affected children did not present with fetal anomalies eligible for PES. CONCLUSIONS: Almost one-fifth of children with a monogenic condition included in this study could have received a diagnosis via modern PES, avoiding a diagnostic odyssey lasting almost 4 years. However, despite having a monogenic condition, over half of the children did not present with any structural anomalies in utero. This demonstrates the degree to which fetal imaging is limited in its ability to reassure parents of the absence of a fetal genetic syndrome. © 2026 The Author(s). Ultrasound in Obstetrics & Gynecology published by John Wiley & Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology.

Humans↗

PubMed2026

Bilateral Conversion Risk in Unilateral Retinoblastoma Using Age and Genetic Testing.

IMPORTANCE: Metachronous bilateral conversion in initially unilateral retinoblastoma is uncommon but clinically consequential, potentially requiring intensified treatment and carrying worse prognosis. Clarifying how age at diagnosis refines genetic-risk stratification could enable safer, more efficient surveillance protocols. OBJECTIVE: To estimate the incidence and timing of metachronous bilateral conversion in unilateral retinoblastoma and assess whether age at diagnosis and RB1 testing are associated with bilateral conversion risk. DESIGN, SETTING AND PARTICIPANTS: This was a retrospective cohort study at a tertiary center in Shanghai, China, including 1108 consecutive children with initially unilateral retinoblastoma diagnosed from July 2010 to October 2024 (after exclusions for short follow-up [n = 139], missing data [n = 53], or synchronous bilateral disease [n = 10]). The median (IQR) follow-up was 43.4 (24.2-67.6) months. EXPOSURES: Age at diagnosis and RB1 genetic status/subtypes assessed by next-generation sequencing and multiplex ligation-dependent probe amplification, including penetrance class (high vs low) and mosaic vs germline categorization. MAIN OUTCOMES AND MEASURES: Time to metachronous bilateral conversion; cumulative incidence functions with death as a competing risk; spatial distribution of fellow-eye tumors. RESULTS: Among 1108 patients (median [IQR] age at diagnosis, 22.2 [12.0-31.4] months; 591 [53.3%] male), 24 (2.2%) developed metachronous bilateral disease. At 24 months, cumulative incidence was 2.2% (95% CI, 1.3-3.1) overall. By genetic status, the 24-month cumulative incidence was 24.8% (95% CI, 13.8-35.9) in RB1 variant-positive vs 1.6% (95% CI, 0.0-3.1) in RB1 variant-negative patients. Among RB1 variant-positive patients, risk clustered among those diagnosed before 9 months, whereas no conversions were observed among those diagnosed at older than 9 months. Four RB1 variant-negative patients who were initially diagnosed at notably late ages (20.9, 42.7, 79.6, and 118 months) subsequently converted; these cases likely represent undetected low-level mosaicism, somatic variants below detection thresholds, or rare genomic events not captured by standard sequencing panels. Fellow-eye tumors did not involve macula and showed a nasal-predominant distribution. CONCLUSIONS AND RELEVANCE: The findings in this study suggest that age at diagnosis may refine genetic risk stratification for metachronous bilateral conversion. RB1 variant-positive patients diagnosed at 9 months or later represent a very low-risk subgroup that may warrant surveillance deescalation, while rare late conversions in RB1 variant-negative patients necessitate continued long-term monitoring.

Humans↗

PubMed2026

Discovery and Engineering of a Rat Endogenous Retrovirus Reverse Transcriptase for Efficient Prime Editing.

CRISPR-based prime editors (PEs) install precise edits into genomic DNA without generating double-strand breaks. Their editing efficiency is highly dependent on reverse transcriptases (RTs), but efficient RT candidates remain limited. Here, we identified 19 novel active RTs by screening 558 candidates. Among them, RERV-RT, derived from Rattus norvegicus, exhibited the highest activity. Through structure-guided engineering and deep mutational scanning, we developed an optimized variant, enRERV-RT, which outperforms conventional M-MLV-RT-based PE systems by 1.20-fold in mammalian and plant cells, and by 1.88-fold at hard-to-edit loci, while enabling precise multiplex editing of functionally relevant genes. Additionally, we developed a high-throughput platform, TRAP-seq-PE, to systematically evaluate prime editor performance. Across diverse mutation types, we found that PE systems based on enRERV-RT exhibited higher editing efficiencies than those based on M-MLV-RT. Collectively, our work establishes a versatile, high-efficiency PE system, thereby facilitating advances in clinical gene therapy and precise crop breeding.

Animals↗

PubMed2026

Human Systems Immunology in the Omics Era: Challenges, Methods, and Emerging Directions.

The human immune system is a highly complex, dynamic, and heterogeneous network shaped by genetic, environmental, and temporal influences. Advances in high-throughput omics technologies have transformed our ability to study this complexity directly and comprehensively in human cohorts. These developments have positioned systems immunology as a powerful framework for investigating coordinated immune responses, identifying regulatory mechanisms, and linking molecular patterns to clinical phenotypes. However, the analytical challenges inherent to large-scale, multimodal datasets-including batch effects, small sample sizes, high dimensionality, and substantial interindividual heterogeneity-require rigorous study design, robust statistical modeling, and thoughtful data analysis strategies. In this review, we summarize key technological foundations enabling modern human systems immunology, outline common analytical pitfalls and effective mitigation approaches, discuss data integration concepts, and highlight emerging opportunities in the field. Together, these technological and analytical advances are redefining how immune function is measured and interpreted in real-world human biology and hold significant promise for enhancing mechanistic insight, biomarker discovery, and precision medicine across immunological diseases and interventions.

Humans↗

PubMed2025

Mapping Active RNA Polymerases in Proliferating and Quiescent Fission Yeast Cells Using Precision Run-On Sequencing.

The development of next-generation sequencing (NGS) approaches to investigate the functioning of RNA polymerases has led to groundbreaking advances in the field of transcriptional regulation. One powerful method, Precision nuclear Run-On sequencing (PRO-seq), maps the locations of RNA polymerase active sites genome-wide at high resolution. PRO-seq provides a snapshot of strand-specific transcriptional activity and does not rely on immunoprecipitation of the polymerase of interest. Notably, this technique has been utilized to investigate the control of the RNA polymerase II transcription cycle in a variety of model systems. However, the initially published PRO-seq method required significant amounts of starting sample and was technically challenging, both of which were deterrents for its broader use. Recently, an improved and simplified version called qPRO-seq that reduced the length of the experiment and the quantity of necessary input sample was developed for human and Drosophila cell lines. Here we provide an updated, step-by-step protocol in which we have validated and optimized qPRO-seq for the fission yeast Schizosaccharomyces pombe. Importantly, we have implemented this method for assessing RNA polymerase activity in nutrient-limiting conditions, for both proliferating and nitrogen-depleted quiescent cells.

Schizosaccharomyces↗

PubMed2025

Sawfish: improving long-read structural variant discovery and genotyping with local haplotype modeling.

MOTIVATION: Structural variants (SVs) play an important role in evolutionary and functional genomics but are challenging to characterize. High-accuracy, long-read sequencing can substantially improve SV characterization when coupled with effective calling methods. While state-of-the-art long-read SV callers are highly accurate, further improvements are achievable by systematically modeling local haplotypes during SV discovery and genotyping. RESULTS: We describe sawfish, an SV caller for mapped high-quality long reads incorporating systematic SV haplotype modeling to improve accuracy and resolution. Assessment against the draft Genome in a Bottle (GIAB) SV benchmark from the T2T-HG002-Q100 diploid assembly shows that sawfish has the highest accuracy among state-of-the-art long-read SV callers across every tested SV size group. Additionally, sawfish maintains the highest accuracy at every tested depth level from 10- to 32-fold coverage, such that other callers required at least 30-fold coverage to match sawfish accuracy at 15-fold coverage. Sawfish also shows the highest accuracy in the GIAB challenging medically relevant genes benchmark, demonstrating improvements in both comprehensive and medically relevant contexts.When joint-genotyping seven samples from CEPH-1463, sawfish has over 9000 more pedigree-concordant calls than other state-of-the-art SV callers, with the highest proportion of concordant SVs (81%). Sawfish's quality model enables selection for an even higher proportion of concordant SVs (88%), while still calling nearly 5000 more pedigree-concordant SVs than other callers. These results demonstrate that sawfish improves on the state-of-the-art for long-read SV calling accuracy across both individual and joint-sample analyses. AVAILABILITY AND IMPLEMENTATION: Sawfish source code, pre-compiled Linux binaries, and documentation are released on GitHub: https://github.com/PacificBiosciences/sawfish.

Haplotypes↗

PubMed2025

Genome-wide profiling the integration patterns with T7-PCR.

Integration of exogenous gene fragments into the host genomes is a widely used and powerful method for studying gene functions, advancing molecular breeding, and conducting gene therapy. Accurately identifying the integration sites is essential for ensuring both the safety and efficacy of genome engineering efforts. However, current mapping techniques are constrained by high costs and a low signal-to-noise ratio. In this study, we developed an innovative tool for mapping integration sites, leveraging T7 polymerase-mediated in vitro transcription (T7-IVT) to capture the junction fragments surrounding integration loci. This approach converts genomic flanking sequences into RNA, enabling the simultaneous enrichment of junction fragments and the elimination of background genomic DNA, thereby significantly enhancing the signal-to-noise ratio. We have validated the efficiency of this method, named T7-PCR, across yeast, plant, and human cells under diverse integration scenarios. T7-PCR outperforms current next-generation sequencing (NGS)-based mapping strategies in terms of efficiency and accuracy, with minimal positional effects. This method is highly applicable for high-throughput transgene screening and also supports the development of next-generation tools for targeted integration of large fragments.

Humans↗