Abstract:Phasing of heterozygous alleles is critical for interpretation of cis -effects of disease-relevant variation. For population studies, phase is often inferred from external data but read-based phasing approaches that span long genomic distances would be more accurate because they enable both genotype and phase to be obtained from a single dataset. To demonstrate how read-based phasing can provide functional insights, we sequenced 477 individuals with Cystic Fibrosis (CF) using linked-read sequencing. We benchmark read-based phasing with different short- and long-read sequencing technologies, prioritize linked-read technology as the most informative and produce a benchmark phase call set from reference sample HG002 for the community. The 477 samples display an average phase block N50 of 4.39 Mb. We use these samples to construct a graph representation of CFTR haplotypes, which facilitates understanding of complex CF alleles. Fine-mapping and phasing of the chr7q35 trypsinogen locus associated with CF meconium ileus demonstrates a 20 kb deletion and a PRSS2 missense variant p.Thr8Ile (rs62473563) independently contribute to meconium ileus risk (p=0.0028, p=0.011, respectively) and are PRSS2 pancreas eQTLs (p=9.5e-7 and p=1.4e-4, respectively), explaining the mechanism by which these polymorphisms contribute to CF. Phase enables access to haplotypes that can be used for genome graph or reference panel construction, identification of cis -effects, and for understanding disease associated loci. The phase information from linked-reads provides a causal explanation for variation at a CF-relevant locus which also has implications for the genetic basis of non-CF pancreatitis to which this locus has been reported to contribute. ### Competing Interest Statement DMC received an honorarium for teaching module development for Vertex Pharmaceuticals. NM is doing contract research trials for Vertex Phaemaceuticals and Abbvie. ALS has received speaking fees for educational programs sponsored by Vertex Pharmaceuticals. BSQ has received speaker fees from Vertex Pharmaceuticals and has served as site PI for several Vertex-sponsored clinical trials. WML is a study investigator for Vertex Pharmaceuticals. ET and FR act as consultants for Vertex Pharmaceuticals. MS participated in Vertex clinical trials and received payment for education modules. SM, AC, JG, FL, BT, WWLS, JW, ZW, RVP, KK, AH, NP, JA, CW, GCM, SB, DA, EB, CB, MC, AP, MP, RVW, DH, MJS, ET, PW, LS, FR, and LJS have no conflicts of interest. * CF : cystic fibrosis CFTR : cystic fibrosis transmembrane conductance WGS : whole genome sequencing GWAS : genome-wide association studies MI : meconium ileus LD : linkage disequilibrium PacBio : Pacific Biosciences 10XG : 10x Genomics CGMS : Canadian CF Gene Modifier Study Consortium CLR : PacBio continuous long-reads CCS : PacBio circular consensus sequence VCF : variant call format GIAB : Genome in a Bottle HMW : high molecular weight TCAG : The Centre for Applied Genomics GTEx : Genotype-Tissue Expression ER : endoplasmic reticulum SRP : signal recognition particle QC : quality control

StarPhase: Comprehensive Phase-Aware Pharmacogenomic Diplotyper for Long-Read Sequencing Data

ClinPharmSeq: A targeted sequencing panel for clinical pharmacogenetics implementation

HiPhase: jointly phasing small, structural, and tandem repeat variants from HiFi sequencing

Universal Targeted Haplotyping by Droplet Digital PCR Sequencing and Its Applications in Noninvasive Prenatal Testing and Pharmacogenetics Analysis

Pharmacogenetic analysis of structural variation in the 1000 genomes project using whole genome sequences

A systematic comparison of pharmacogene star allele calling bioinformatics algorithms: a focus on CYP2D6 genotyping

Closing the gap: Solving complex medically relevant genes at scale

High Quality Phasing Using Linked-Read Whole Genome Sequencing of Patient Cohorts Informs Genetic Understanding of Complex Traits

Haplotype phasing of CYP2D6 : an allelic ratio method using Agena MassARRAY data

Barthélémy et al. Direct Comparative Analysis of a Pharmacogenomics Panel with PacBio Hifi® Long-Read and Illumina Short-Read Sequencing. J. Pers. Med. 2023, 13, 1655

High-throughput method for the hybridisation-based targeted enrichment of long genomic fragments for PacBio third-generation sequencing

Approach for Phased Sequence-Based Genotyping of the Critical Pharmacogene Dihydropyrimidine Dehydrogenase (DPYD)

Star allele search: a pharmacogenetic annotation database and user-friendly search tool of publicly available 1000 Genomes Project biospecimens

Characterization of complex structural variation in the CYP2D6-CYP2D7-CYP2D8 gene loci using single-molecule long-read sequencing

Efficient whole genome haplotyping and high-throughput single molecule phasing with barcode-linked reads

Genome-wide profiling of highly similar paralogous genes using HiFi sequencing

Contradiction in Star-Allele Nomenclature of Pharmacogenes between Common Haplotypes and Rare Variants

Whole-genome Haplotyping Using Long Reads and Statistical Methods.

Targeted phasing of 2–200 kilobase DNA fragments with a short-read sequencer and a single-tube linked-read library method

Phasing millions of samples achieves near perfect accuracy, enabling parent-of-origin classification of variants

Characterization of 137 Genomic DNA Reference Materials for 28 Pharmacogenetic Genes: A GeT-RM Collaborative Project