PacBio long-read amplicon sequencing enables scalable high-resolution population allele typing of the complex CYP2D6 locus

Closing the gap: Solving complex medically relevant genes at scale

Comprehending the mechanism behind human diseases with an established heritable component represents the forefront of personalized medicine. Nevertheless, numerous medically important genes are inaccurately represented in short-read sequencing data analysis due to their complexity and repetitiveness or the so-called 'dark regions' of the human genome. The advent of PacBio as a long-read platform has provided new insights, yet HiFi whole-genome sequencing (WGS) cost remains frequently prohibitive. We introduce a targeted sequencing and analysis framework, Twist Alliance Dark Genes Panel (TADGP), designed to offer phased variants across 389 medically important yet complex autosomal genes. We highlight TADGP accuracy across eleven control samples and compare it to WGS. This demonstrates that TADGP achieves variant calling accuracy comparable to HiFi-WGS data, but at a fraction of the cost. Thus, enabling scalability and broad applicability for studying rare diseases or complementing previously sequenced samples to gain insights into these complex genes. TADGP revealed several candidate variants across all cases and provided insight into LPA diversity when tested on samples from rare disease and cardiovascular disease cohorts. In both cohorts, we identified novel variants affecting individual disease-associated genes (e.g., IKZF1, KCNE1). Nevertheless, the annotation of the variants across these 389 medically important genes remains challenging due to their underrepresentation in ClinVar and gnomAD. Consequently, we also offer an annotation resource to enhance the evaluation and prioritization of these variants. Overall, we can demonstrate that TADGP offers a cost-efficient and scalable approach to routinely assess the dark regions of the human genome with clinical relevance.

PharmVar Tutorial on CYP2D6 Structural Variation Testing and Recommendations on Reporting

The Pharmacogene Variation Consortium (PharmVar) provides nomenclature for the highly polymorphic human CYP2D6 gene locus and a comprehensive summary of structural variation. CYP2D6 contributes to the metabolism of numerous drugs and, thus, genetic variation in its gene impacts drug efficacy and safety. To accurately predict a patient's CYP2D6 phenotype, testing must include structural variants including gene deletions, duplications, hybrid genes, and combinations thereof. This tutorial offers a comprehensive overview of CYP2D6 structural variation, terms, and definitions, a review of methods suitable for their detection and characterization, and practical examples to address the lack of standards to describe CYP2D6 structural variants or any other pharmacogene. This PharmVar tutorial offers practical guidance on how to detect the many, often complex, structural variants, as well as recommends terms and definitions for clinical and research reporting. Uniform reporting is not only essential for electronic health record-keeping but also for accurate translation of a patient's genotype into phenotype which is typically utilized to guide drug therapy.

Case Series of Primaquine-Induced Haemolytic Events in Controlled Trials with G6PD Screening

Primaquine for radical cure of Plasmodium vivax malaria poses a potentially life-threatening risk of haemolysis in G6PD-deficient patients. Herein, we review five events of acute haemolytic anaemia following the administration of primaquine in four malaria trials from Indonesia, the Solomon Islands, and Vietnam. Five males aged 9 to 48 years were improperly classified as G6PD-normal by various screening procedures and included as subjects in trials of anti-relapse therapy with daily primaquine. Routine safety monitoring by physical examination, urine inspection, and blood haemoglobin (Hb) assessment were performed in all those trials. Early signs of acute haemolysis, i.e., dark urine and haemoglobin drop >20%, occurred only after day 3 and as late as day 8 of primaquine dosing. All patients were hospitalized and fully recovered, all but one following blood transfusion rescue. Hb nadir was 4.7 to 7.9 g/dL. Hospitalization was for 1 to 7 days. Hb levels returned to baseline values 3 to 10 days after transfusion. Failed G6PD screening procedures in these trials led G6PD-deficient patients to suffer harmful exposures to primaquine. The safe application of primaquine anti-relapse therapy requires G6PD screening and anticipation of its failure with a means of prompt detection and rescue from the typically abrupt haemolytic crisis.

Applications of advanced technologies for detecting genomic structural variation

Chromosomal structural variation (SV) encompasses a heterogenous class of genetic variants that exerts strong influences on human health and disease. Despite their importance, many structural variants (SVs) have remained poorly characterized at even a basic level, a discrepancy predicated upon the technical limitations of prior genomic assays. However, recent advances in genomic technology can identify and localize SVs accurately, opening new questions regarding SV risk factors and their impacts in humans. Here, we first define and classify human SVs and their generative mechanisms, highlighting characteristics leveraged by various SV assays. We next examine the first-ever gapless assembly of the human genome and the technical process of assembling it, which required third-generation sequencing technologies to resolve structurally complex loci. The new portions of that "telomere-to-telomere" and subsequent pangenome assemblies highlight aspects of SV biology likely to develop in the near-term. We consider the strengths and limitations of the most promising new SV technologies and when they or longstanding approaches are best suited to meeting salient goals in the study of human SV in population-scale genomics research, clinical, and public health contexts. It is a watershed time in our understanding of human SV when new approaches are expected to fundamentally change genomic applications.

Opinion paper on the systematic application of integrated bioinformatic tools to actuate routine precision medicine in poly-treated patients

Precision Medicine is a reality in selected medical areas, as oncology, or in excellent healthcare structures, but it is still far to reach million patients who could benefit from this medical concept. Here, we sought to highlight how the time is ripe to achieve horizontal delivery to a significant larger audience of patients, represented by the poly-treated patients. Combination therapies are frequent (especially in the elderly, to treat comorbidities) and are related to decreased drug safety and efficacy, disease's exacerbation, additional treatments, hospitalization. But the recent development and validation of bioinformatic tools, aimed to automatic evaluation and optimization of poly-therapies, according to the unique individual characteristics (including genotype), is ready to change the daily approach to pharmacological prescription.

Allele mining, amplicon sequencing and computational prediction of Solanum melongena L. FT/TFL1 gene homologs uncovers putative variants associated to seed dormancy and germination

The FT/TFL1 gene homolog family plays a crucial role in the regulation of floral induction, seed dormancy and germination in angiosperms. Despite its importance, the FT/TFL1 gene homologs in eggplant (Solanum melongena L.) have not been characterized to date. In this study, we performed a genome-wide identification of FT/TFL1 genes in eggplant using in silico genome mining. The presence of these genes was validated in four economically important eggplant cultivars (Surya, EP-47 Annamalai, Pant Samrat and Arka Nidhi) through Pacbio RSII amplicon sequencing. Our results revealed the presence of 12 FT/TFL1 gene homologs in eggplant, with evidence of diversification among FT-like genes suggesting their possible adaptations towards various environmental stimuli. The amplicon sequencing also revealed the presence of two alleles for certain genes (SmCEN-1, SmCEN-2, SmMFT-1 and SmMFT-2) of which SmMFT-2 was associated with seed dormancy and germination. This association was further supported by the observation that seed dormancy is rarely reported in domesticated eggplant cultivars, but is commonly observed in wild species. A survey of the genetic regions in domesticated cultivars and a related wild species, S. incanum, showed that the alternative allele of S. incanum was present in some members of the Pant Samrat cultivar, but was absent in most other cultivars. This difference could contribute to the differences in seed traits between wild and domesticated eggplants.

DNA methylation episignatures: insight into copy number variation

In this review we discuss epigenetic disorders that result from aberrations in genes linked to epigenetic regulation. We describe current testing methods for the detection of copy number variants (CNVs) in Mendelian disorders, dosage sensitivity, reciprocal phenotypes and the challenges of test selection and overlapping clinical features in genetic diagnosis. We discuss aberrations of DNA methylation and propose a role for episignatures as a novel clinical testing method in CNV disorders. Finally, we postulate that episignature mapping in CNV disorders may provide novel insights into the molecular mechanisms of disease and unlock key findings of the genome-wide impact on disease gene networks.

Allelic diversity of the pharmacogene CYP2D6 in New Zealand Māori and Pacific peoples

The enzyme cytochrome P450 2D6 (CYP2D6) metabolises approximately 25% of commonly prescribed drugs, including analgesics, anti-hypertensives, and anti-depressants, among many others. Genetic variation in drug metabolising genes can alter how an individual responds to prescribed drugs, including predisposing to adverse drug reactions. The majority of research on the CYP2D6 gene has been carried out in European and East Asian populations, with many Indigenous and minority populations, such as those from Oceania, greatly underrepresented. However, genetic variation is often population specific and analysis of diverse ethnic groups can reveal differences in alleles that may be of clinical significance. For this reason, we set out to examine the range and frequency of CYP2D6 variants in a sample of 202 Māori and Pacific people living in Aotearoa (New Zealand). We carried out long PCR to isolate the CYP2D6 region before performing nanopore sequencing to identify all variants and alleles in these samples. We identified twelve variants which have previously not been reported in the PharmVar CYP2D6 database, three of which were exonic missense variations. Six of these occurred in single samples and one was found in 19 samples (9.4% of the cohort). The remaining five variants were identified in two samples each. Identified variants formed twelve new CYP2D6 suballeles and four new star alleles, now recorded in the PharmVar database. One striking finding was that CYP2D6*71, an allele of uncertain functional status which has been rarely observed in previous studies, occurs at a relatively high frequency (8.9%) within this cohort. These data will help to ensure that CYP2D6 genetic analysis for pharmacogenetic purposes can be carried out accurately and effectively in this population group.

Challenges and opportunities associated with rare-variant pharmacogenomics

Recent advances in next-generation sequencing (NGS) have resulted in the identification of tens of thousands of rare pharmacogenetic variations with unknown functional effects. However, although such pharmacogenetic variations have been estimated to account for a considerable amount of the heritable variability in drug response and toxicity, accurate interpretation at the level of the individual patient remains challenging. We discuss emerging strategies and concepts to close this translational gap. We illustrate how massively parallel experimental assays, artificial intelligence (AI), and machine learning can synergize with population-scale biobank projects to facilitate the interpretation of NGS data to individualize clinical decision-making and personalized medicine.