Transcriptional analysis of sodium valproate in a serotonergic cell line reveals gene regulation through both HDAC inhibition-dependent and independent mechanisms

Sodium valproate (VPA) is a histone deacetylase (HDAC) inhibitor, widely prescribed in the treatment of bipolar disorder, and yet the precise modes of therapeutic action for this drug are not fully understood. After exposure of the rat serotonergic cell line RN46A to VPA, RNA-sequencing (RNA-Seq) analysis showed widespread changes in gene expression. Analysis by four bioinformatic pipelines revealed as many as 230 genes were significantly upregulated and 72 genes were significantly downregulated. A subset of 23 differentially expressed genes was selected for validation using the nCounter® platform, and of these we obtained robust validation for ADAM23, LSP1, MAOB, MMP13, PAK3, SERPINB2, SNAP91, WNT6, and ZCCHC12. We investigated the effect of lithium on this subset and found four genes, CDKN1C, LSP1, SERPINB2, and WNT6 co-regulated by lithium and VPA. We also explored the effects of other HDAC inhibitors and the VPA analogue valpromide on the subset of 23 selected genes. Expression of eight of these genes, CDKN1C, MAOB, MMP13, NGFR, SHANK3, VGF, WNT6 and ZCCHC12, was modified by HDAC inhibition, whereas others did not appear to respond to several HDAC inhibitors tested. These results suggest VPA may regulate genes through both HDAC-dependent and independent mechanisms. Understanding the broader gene regulatory effects of VPA in this serotonergic cell model should provide insights into how this drug works and whether other HDAC inhibitor compounds may have similar gene regulatory effects, as well as highlighting molecular processes that may underlie regulation of mood.

G-quadruplex structures bind to EZ-Tn5 transposase

Next generation DNA sequencing and analysis of amplicons spanning the pharmacogene CYP2D6 suggested that the Nextera transposase used for fragmenting and providing sequencing priming sites displayed a targeting bias. This manifested as dramatically lower sequencing coverage at sites in the amplicon that appeared likely to form G-quadruplex structures. Since secondary DNA structures such as G-quadruplexes are abundant in the human genome, and are known to interact with many other proteins, we further investigated these sites of low coverage. Our investigation revealed that G-quadruplex structures are formed in vitro within the CYP2D6 pharmacogene at these sites, and G-quadruplexes can interact with the hyperactive Tn5 transposase (EZ-Tn5) with high affinity. These findings indicate that secondary DNA structures such as G-quadruplexes may represent preferential transposon integration sites and provide additional evidence for the role of G-quadruplex structures in transposition or viral integration processes.

Nanopore sequencing of the pharmacogene CYP2D6 allows simultaneous haplotyping and detection of duplications

Aim: Long read sequencing offers the promise of overcoming some of the challenges in accurate genotyping of complex genes, along with the advantage of straightforward variant phasing. We have established methods for sequencing and haplotyping of the whole CYP2D6 gene using nanopore sequencing. Materials and methods: 32 samples covering various haplotypes including gene duplication were sequenced on the GridION platform. Results: Haplotypes of 52 alleles matched accurately to known star (*) allele subvariants, with the remaining 12 being assigned as new alleles, or new subvariants of known alleles. Duplicated alleles could be detected by analyzing the allelic balance. Conclusion: Nanopore sequencing of CYP2D6 offers a high throughput method for accurate haplotyping, detection of new variants and determination of duplicated alleles.

Comprehensive Assessment of BARD1 Messenger Ribonucleic Acid Splicing With Implications for Variant Classification

Introduction: Case-control analyses have shown BARD1 variants to be associated with up to >2-fold increase in risk of breast cancer, and potentially greater risk of triple negative breast cancer. BARD1 is included in several gene sequencing panels currently marketed for the prediction of risk of cancer, however there are no gene-specific guidelines for the classification of BARD1 variants. We present the most comprehensive assessment of BARD1 messenger RNA splicing, and demonstrate the application of these data for the classification of truncating and splice site variants according to American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines. Methods: Nanopore sequencing, short-read RNA-seq (whole transcriptome and targeted), and capillary electrophoresis analysis were performed by four laboratories to investigate alternative BARD1 splicing in blood, breast, and fimbriae/ovary related specimens from non-cancer affected tissues. Splicing data were also collated from published studies of nine different tissues. The impact of the findings for PVS1 annotation was assessed for truncating and splice site variants. Results: We identified 62 naturally occurring alternative spliced BARD1 splicing events, including 19 novel events found by next generation sequencing and/or reverse transcription PCR analysis performed for this study. Quantitative analysis showed that naturally occurring splicing events causing loss of clinically relevant domains or nonsense mediated decay can constitute up to 11.9% of overlapping natural junctions, suggesting that aberrant splicing can be tolerated up to this level. Nanopore sequencing of whole BARD1 transcripts characterized 16 alternative isoforms from healthy controls, revealing that the most complex transcripts combined only two alternative splicing events. Bioinformatic analysis of ClinVar submitted variants at or near BARD1 splice sites suggest that all consensus splice site variants in BARD1 should be considered likely pathogenic, with the possible exception of variants at the donor site of exon 5. Conclusions: No BARD1 candidate rescue transcripts were identified in this study, indicating that all premature translation-termination codons variants can be annotated as PVS1. Furthermore, our analysis suggests that all donor and acceptor (IVS+/-1,2) variants can be considered PVS1 or PVS1_strong, with the exception of variants targeting the exon 5 donor site, that we recommend considering as PVS1_moderate.

Multiplexed Nanopore Sequencing of HLA-B Locus in Māori and Pacific Island Samples

The human leukocyte antigen (HLA) system encodes the human major histocompatibility complex (MHC). HLA-B is the most polymorphic gene in the MHC class I region and many HLA-B alleles have been associated with adverse drug reactions (ADRs) and disease susceptibility. The frequency of such HLA-B alleles varies by ethnicity, and therefore it is important to understand the prevalence of such alleles in different population groups. Research into HLA involvement in ADRs would be facilitated by improved methods for genotyping key HLA-B alleles. Here, we describe an approach to HLA-B typing using next generation sequencing (NGS) on the MinION™ nanopore sequencer, combined with data analysis with the SeqNext-HLA software package. The nanopore sequencer offers the advantages of long-read capability and single molecule reads, which can facilitate effective haplotyping. We developed this method using reference samples as well as individuals of New Zealand Māori or Pacific Island descent, because HLA-B diversity in these populations is not well understood. We demonstrate here that nanopore sequencing of barcoded, pooled, 943 bp polymerase chain reaction (PCR) amplicons of 49 DNA samples generated ample read depth for all samples. HLA-B alleles were assigned to all samples at high-resolution with very little ambiguity. Our method is a scaleable and efficient approach for genotyping HLA-B and potentially any other HLA locus. Finally, we report our findings on HLA-B genotypes of this cohort, which adds to our understanding of HLA-B allele frequencies among Māori and Pacific Island people.

DNA G-quadruplexes show strong interaction with DNA methyltransferases in vitro

The DNA methyltransferase enzymes (DNMTs) catalyzing cytosine methylation do so at specific locations of the genome, although with some level of redundancy. The de novo methyltransferases DNMT3A and 3B play a vital role in methylating the genome of the developing embryo in regions devoid of methylation marks. The ability of DNMTs to colocalize at sites of DNA damage is suggestive that recognition of mispaired bases and unusual structures is inherent to the function of these proteins. We provide evidence for G-quadruplex formation within imprinted gene promoters, and report high-affinity binding of recombinant human DNMTs to such DNA G-quadruplexes in vitro. These observations suggest a potential interaction of G-quadruplexes with the DNA methylation machinery, which may be of epigenetic and biological significance.

Cross-Comparison of Exome Analysis, Next-Generation Sequencing of Amplicons, and the iPLEX(®) ADME PGx Panel for Pharmacogenomic Profiling

Whole-exome sequencing (WES) has been widely used for analysis of human genetic diseases, but its value for the pharmacogenomic profiling of individuals is not well studied. Initially, we performed an in-depth evaluation of the accuracy of WES variant calling in the pharmacogenes CYP2D6 and CYP2C19 by comparison with MiSeq(®) amplicon sequencing data (n = 36). This analysis revealed that the concordance rate between WES and MiSeq(®) was high, achieving 99.60% for variants that were called without exceeding the truth-sensitivity threshold (99%), defined during variant quality score recalibration (VQSR). Beyond this threshold, the proportion of discordant calls increased markedly. Subsequently, we expanded our findings beyond CYP2D6 and CYP2C19 to include more genes genotyped by the iPLEX(®) ADME PGx Panel in the subset of twelve samples. WES performed well, agreeing with the genotyping panel in approximately 99% of the selected pass-filter variant calls. Overall, our results have demonstrated WES to be a promising approach for pharmacogenomic profiling, with an estimated error rate of lower than 1%. Quality filters, particularly VQSR, are important for reducing the number of false variants. Future studies may benefit from examining the role of WES in the clinical setting for guiding drug therapy.

G-quadruplex structures and CpG methylation cause drop-out of the maternal allele in polymerase chain reaction amplification of the imprinted MEST gene promoter

We observed apparent non-Mendelian behaviour of alleles when genotyping a region in a CpG island at the 5′ end of the maternally imprinted human MEST isoform. This region contains three single nucleotide polymorphisms (SNPs) in total linkage disequilibrium, such that only two haplotypes occur in the human population. Only one haplotype was detectable in each subject, never both, despite the use of multiple primers and several genotyping methods. We observed that this region contains motifs capable of forming several G-quadruplex structures. Circular dichroism spectroscopy and native polyacrylamide gel electrophoresis confirmed that at least three G-quadruplexes form in vitro in the presence of potassium ions, and one of these structures has a T_m of greater than 99°C in polymerase chain reaction (PCR) buffer. We demonstrate that it is the methylated maternal allele that is always lost during PCR amplification, and that formation of G-quadruplexes and presence of methylated cytosines both contributed to this phenomenon. This observed parent-of-origin specific allelic drop-out has important implications for analysis of imprinted genes in research and diagnostic settings.

Relevance of G-quadruplex structures to pharmacogenetics

G-quadruplexes are non-canonical secondary structures formed within nucleic acids that are involved in modulating cellular processes such as replication, gene regulation, recombination and epigenetics. Within genes, there is mounting evidence of G-quadruplex involvement in transcriptional and post transcriptional regulation. We report the presence of potential G-quadruplex motifs within relevant sites of some important pharmacogenes and discuss the possible implications of this on the function and expression of these genes. Appreciating the location and potential functions of these motifs may be of value when considering the impacts of some pharmacogenetic variants. G-quadruplexes are also the focus of drug development efforts in oncology and we highlight the broader pharmacological implications of treatment strategies that may target G-quadruplexes.