Splice-Break: exploiting an RNA-seq splice junction algorithm to discover mitochondrial DNA deletion breakpoints and analyses of psychiatric disorders.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

2023 White Paper on Recent Issues in Bioanalysis: Deuterated Drugs; LNP; Tumor/FFPE Biopsy; Targeted Proteomics; Small Molecule Covalent Inhibitors; Chiral Bioanalysis; Remote Regulatory Assessments; Sample Reconciliation/Chain of Custody (PART 1A - Recommendations on Mass Spectrometry, Chromatography, Sample Preparation Latest Developments, Challenges, and Solutions and BMV/Regulated Bioanalysis PART 1B - Regulatory Agencies' Inputs on Regulated Bioanalysis/BMV, Biomarkers/IVD/CDx/BAV, Immunogenicity, Gene & Cell Therapy and Vaccine).

The 17th Workshop on Recent Issues in Bioanalysis (17th WRIB) took place in Orlando, FL, USA on June 19-23, 2023. Over 1000 professionals representing pharma/biotech companies, CROs, and multiple regulatory agencies convened to actively discuss the most current topics of interest in bioanalysis. The 17th WRIB included 3 Main Workshops and 7 Specialized Workshops that together spanned 1 week to allow an exhaustive and thorough coverage of all major issues in bioanalysis of biomarkers, immunogenicity, gene therapy, cell therapy and vaccines.Moreover, in-depth workshops on "EU IVDR 2017/746 Implementation and impact for the Global Biomarker Community: How to Comply with this NEW Regulation" and on "US FDA/OSIS Remote Regulatory Assessments (RRAs)" were the special features of the 17th edition.As in previous years, WRIB continued to gather a wide diversity of international, industry opinion leaders and regulatory authority experts working on both small and large molecules as well as gene, cell therapies and vaccines to facilitate sharing and discussions focused on improving quality, increasing regulatory compliance, and achieving scientific excellence on bioanalytical issues.This 2023 White Paper encompasses recommendations emerging from the extensive discussions held during the workshop and is aimed to provide the bioanalytical community with key information and practical solutions on topics and issues addressed, in an effort to enable advances in scientific excellence, improved quality and better regulatory compliance. Due to its length, the 2023 edition of this comprehensive White Paper has been divided into three parts for editorial reasons.This publication covers the recommendations on Mass Spectrometry Assays, Regulated Bioanalysis/BMV (Part 1A) and Regulatory Inputs (Part 1B). Part 2 (Biomarkers, IVD/CDx, LBA and Cell-Based Assays) and Part 3 (Gene Therapy, Cell therapy, Vaccines and Biotherapeutics Immunogenicity) are published in volume 16 of Bioanalysis, issues 7 and 8 (2024), respectively.

Large Common Mitochondrial DNA Deletions Are Associated with a Mitochondrial SNP T14798C Near the 3' Breakpoints.

Introduction: Large somatic deletions of mitochondrial DNA (mtDNA) accumulate with aging in metabolically active tissues such as the brain. We have cataloged the breakpoints and frequencies of large mtDNA deletions in the human brain. Methods: We quantified 112 high-frequency mtDNA somatic deletions across four human brain regions with the Splice-Break2 pipeline. In addition, we utilized PLINK/Seq to test the association of mitochondrial genotypes with the abundance of these high-frequency mtDNA deletions. A conservative p value threshold of 5E-08 was used to find the significant loci. Results: One mtDNA SNP (T14798C) was significantly associated with mtDNA deletions in two brain regions, the dorsolateral prefrontal cortex (DLPFC) and the superior temporal gyrus. Since the DLPFC showed the most robust association between T14798C and two deletion breakpoints (7816-14807 and 5462-14807), this association was tested in the DLPFC of a replication sample and validated the first results. Incorporating the C allele at 14,798 bp increased the perfect/imperfect length of the repeat at the 3' breakpoint of the two associated deletions. Conclusion: This is the first study to identify the association of mtDNA SNP with large mtDNA deletions in the human brain. The T14798C allele located in the MT-CYB gene is a common polymorphism that occurs in several mitochondrial haplogroups. We hypothesize that the T14798C association with two deletions occurs by extending the repeat length around the 3' deletion breakpoints. This simple mechanism suggests that mtDNA SNPs can affect the mitochondrial genome structure, especially in brain where high levels of reactive oxygen species lead to deletion accumulation with aging.

Circular extrachromosomal DNA promotes tumor heterogeneity in high-risk medulloblastoma.

Circular extrachromosomal DNA (ecDNA) in patient tumors is an important driver of oncogenic gene expression, evolution of drug resistance and poor patient outcomes. Applying computational methods for the detection and reconstruction of ecDNA across a retrospective cohort of 481 medulloblastoma tumors from 465 patients, we identify circular ecDNA in 82 patients (18%). Patients with ecDNA-positive medulloblastoma were more than twice as likely to relapse and three times as likely to die within 5 years of diagnosis. A subset of tumors harbored multiple ecDNA lineages, each containing distinct amplified oncogenes. Multimodal sequencing, imaging and CRISPR inhibition experiments in medulloblastoma models reveal intratumoral heterogeneity of ecDNA copy number per cell and frequent putative 'enhancer rewiring' events on ecDNA. This study reveals the frequency and diversity of ecDNA in medulloblastoma, stratified into molecular subgroups, and suggests copy number heterogeneity and enhancer rewiring as oncogenic features of ecDNA.

3D genome mapping identifies subgroup-specific chromosome conformations and tumor-dependency genes in ependymoma.

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

Novel sequential immunocapture microflow LC/MS/MS approach to measuring PTH-Fc protein in human serum.

The quantitation of PTH-Fc in circulation by ligand binding assay presented a significant challenge due to the extremely low doses of administration, interference from the endogenous. A robust LC-MS/MS method to quantify the extremely low concentration of PTH-Fc in human serum utilized sequential immunoaffinity enrichment at PTH and Fc domains in conjunction with microflow LC-MS/MS technology significantly improved the sensitivity and selectivity. The assay displayed a quantitation range of 0.025-5.0 ng/ml and acceptable intraday and interday precision (%CV ≤ 15%) and accuracy (%bias ≤ ±15%) and can be routinely used for pharmacokinetic measurement of the drug. The novel sequential immunocapture workflow described herein can be applied to the quantitation of other recombinant therapeutic proteins to support clinical studies.

Mitochondria DNA copy number, mitochondria DNA total somatic deletions, Complex I activity, synapse number, and synaptic mitochondria number are altered in schizophrenia and bipolar disorder.

Mitochondrial dysfunction is a neurobiological phenomenon implicated in the pathophysiology of schizophrenia and bipolar disorder that can synergistically affect synaptic neurotransmission. We hypothesized that schizophrenia and bipolar disorder share molecular alterations at the mitochondrial and synaptic levels. Mitochondria DNA (mtDNA) copy number (CN), mtDNA common deletion (CD), mtDNA total deletion, complex I activity, synapse number, and synaptic mitochondria number were studied in the postmortem human dorsolateral prefrontal cortex (DLPFC), superior temporal gyrus (STG), primary visual cortex (V1), and nucleus accumbens (NAc) of controls (CON), and subjects with schizophrenia (SZ), and bipolar disorder (BD). The results showed (i) the mtDNA CN is significantly higher in DLPFC of both SZ and BD, decreased in the STG of BD, and unaltered in V1 and NAc of both SZ and BD; (ii) the mtDNA CD is significantly higher in DLPFC of BD while unaltered in STG, V1, and NAc of both SZ and BD; (iii) The total deletion burden is significantly higher in DLPFC in both SZ and BD while unaltered in STG, V1, and NAc of SZ and BD; (iv) Complex I activity is significantly lower in DLPFC of both SZ and BD, which is driven by the presence of medications, with no alteration in STG, V1, and NAc. In addition, complex I protein concentration, by ELISA, was decreased across three cortical regions of SZ and BD subjects; (v) The number of synapses is decreased in DLPFC of both SZ and BD, while the synaptic mitochondria number was significantly lower in female SZ and female BD compared to female controls. Overall, these findings will pave the way to understand better the pathophysiology of schizophrenia and bipolar disorder for therapeutic interventions.

Analysis of whole genome biomarker expression in blood and brain.

The consistency of peripheral gene expression data and the overlap with brain expression has not been evaluated in biomarker discovery, nor has it been reported in multiple tissues from the same subjects on a genome wide transcript level. The effects of processing whole blood, transformation, and passaged cell lines on gene expression profiling was studied in healthy subjects using Affymetrix arrays. Ficoll extracted peripheral blood mononuclear cells (PBMCs), Epstein-Barr virus (EBV) transformed lymphocytes, passaged lymphoblastic cell lines (LCLs), and whole blood from Tempus tubes were compared. There were 6,813 transcripts differentially expressed between different methods of blood preparation. Principal component analysis resolved two partitions involving pre- and post-transformation EBV effects. Combining results from Affymetrix arrays, postmortem subjects' brain and PBMC profiles showed co-expression levels of summarized transcripts for 4,103 of 17,859 (22.9%) RefSeq transcripts. In a control experiment, rat hemi-brain and blood showed similar expression levels for 19% of RefSeq transcripts. After filtering transcripts that were not significantly different in abundance between human cerebellum and PBMCs from the Affymetrix exon array the correlation in mean transcript abundance was high as expected (r = 0.98). Differences in the alternative splicing index in brain and blood were found for about 90% of all transcripts examined. This study demonstrates over 4,100 brain transcripts co-expressed in blood samples can be further examined by in vitro and in vivo experimental studies of blood and cell lines from patients with psychiatric disorders.

Differential time- and NADPH-dependent inhibition of CYP2C19 by enantiomers of fluoxetine.

Fluoxetine [+/--N-methyl-3-phenyl-3-[(alpha, alpha, (-trifluoro-p-tolyl)oxy]-propylamine)] a selective serotonin reuptake inhibitor, is widely used in treating depression and other serotonin-dependent disease conditions. Racemic, (R)- and (S)-fluoxetine are potent reversible inhibitors of CYP2D6, and the racemate has been shown to be a mechanism-based inhibitor of CYP3A4. Racemic fluoxetine also demonstrates time- and concentration-dependent inhibition of CYP2C19 catalytic activity in vitro. In this study, we compared fluoxetine, its (R)- and (S)-enantiomers, ticlopidine, and S-benzylnirvanol as potential time-dependent inhibitors of human liver microsomal CYP2C19. In a reversible inhibition protocol (30 min preincubation with liver microsomes without NADPH), we found (R)-, (S)- and racemic fluoxetine to be moderate inhibitors with IC(50) values of 21, 93, and 27 microM, respectively. However, when the preincubation was supplemented with NADPH, IC(50) values shifted to 4.0, 3.4, and 3.0 microM, respectively resulting in IC(50) shifts of 5.2-, 28-, and 9.3-fold. Ticlopidine showed a 1.8-fold shift in IC(50) value, and S-benzylnirvanol shifted right (0.41-fold shift). Follow-up K(I) and k(inact) determinations with fluoxetine confirmed time-dependent inhibition [K(I) values of 6.5, 47, and 14 microM; k(inact) values of 0.023, 0.085, 0.030 min(-1) for (R)-, (S)-, and racemate, respectively]. Although the (S)-isomer exhibits a much lower affinity for CYP2C19 inactivation relative to the (R)-enantiomer, it exhibits a more rapid rate of inactivation. Racemic norfluoxetine exhibited an 11-fold shift (18-1.5 microM) in IC(50) value, suggesting that conversion of fluoxetine to this metabolite represents a metabolic pathway leading to time-dependent inhibition. These data provide an improved understanding of the drug-interaction potential of fluoxetine.

Mitochondrial Complex I Deficiency in Schizophrenia and Bipolar Disorder and Medication Influence.

Subjects with schizophrenia (SZ) and bipolar disorder (BD) show decreased protein and transcript levels for mitochondrial complex I. In vitro results suggest antipsychotic and antidepressant drugs may be responsible. We measured complex I activity in BD, SZ, and controls and presence of antipsychotic and antidepressant medications, mitochondrial DNA (mtDNA) copy number, and the mtDNA "common deletion" in the brain. Complex I activity in the prefrontal cortex was decreased by 45% in SZ compared to controls (p = 0.02), while no significant difference was found in BD. Complex I activity was significantly decreased (p = 0.01) in pooled cases (SZ and BD) that had detectable psychotropic medications and drugs compared to pooled cases with no detectable levels. Subjects with age at onset in their teens and psychotropic medications showed decreased (p < 0.05) complex I activity compared to subjects with an adult age at onset. Both SZ and BD groups displayed significant increases (p < 0.05) in mtDNA copy number compared to controls; however, common deletion burden was not altered. Complex I deficiency is found in SZ brain tissue, and psychotropic medications may play a role in mitochondrial dysfunction. Studies of medication-free first-episode psychosis patients are needed to elucidate whether mitochondrial pathophysiology occurs independent of medication effects.

Quantitative Trait Locus and Brain Expression of HLA-DPA1 Offers Evidence of Shared Immune Alterations in Psychiatric Disorders.

Genome-wide association studies of schizophrenia encompassing the major histocompatibility locus (MHC) were highly significant following genome-wide correction. This broad region implicates many genes including the MHC complex class II. Within this interval we examined the expression of two MHC II genes (HLA-DPA1 and HLA-DRB1) in brain from individual subjects with schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and controls by differential gene expression methods. A third MHC II mRNA, CD74, was studied outside of the MHC II locus, as it interacts within the same immune complex. Exon microarrays were performed in anterior cingulate cortex (ACC) in BD compared to controls, and both HLA-DPA1 and CD74 were decreased in expression in BD. The expression of HLA-DPA1 and CD74 were both reduced in hippocampus, amygdala, and dorsolateral prefrontal cortex regions in SZ and BD compared to controls by specific qPCR assay. We found several novel HLA-DPA1 mRNA variants spanning HLA-DPA1 exons 2-3-4 as suggested by exon microarrays. The intronic rs9277341 SNP was a significant cis expression quantitative trait locus (eQTL) that was associated with the total expression of HLA-DPA1 in five brain regions. A biomarker study of MHC II mRNAs was conducted in SZ, BD, MDD, and control lymphoblastic cell lines (LCL) by qPCR assay of 87 subjects. There was significantly decreased expression of HLA-DPA1 and CD74 in BD, and trends for reductions in SZ in LCLs. The discovery of multiple splicing variants in brain for HLA-DPA1 is important as the HLA-DPA1 gene is highly conserved, there are no reported splicing variants, and the functions in brain are unknown. Future work on the function and localization of MHC Class II proteins in brain will help to understand the role of alterations in neuropsychiatric disorders. The HLA-DPA1 eQTL is located within a large linkage disequilibrium block that has an irrefutable association with schizophrenia. Future tests in a larger cohort are needed to determine the significance of this eQTL association with schizophrenia. Our findings support the long-held hypothesis that alterations in immune function are associated with the pathophysiology of psychiatric disorders.

The somatic common deletion in mitochondrial DNA is decreased in schizophrenia.

Large deletions in mitochondrial DNA (mtDNA) can occur during or result from oxidative stress leading to a vicious cycle that increases reactive oxygen species (ROS) damage and decreases mitochondrial function, thereby causing further oxidative stress. The objective of this study was to determine if disease specific brain differences of the somatic mtDNA common deletion (4977 bp) could be observed in major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SZ) compared to a control group. The accumulation of the mtDNA common deletion was measured using a quantitative assay across 10 brain regions (anterior cingulate cortex, amygdala, caudate nucleus, dorsolateral prefrontal cortex, hippocampus, nucleus accumbens, orbitofrontal cortex, putamen, substantia nigra, and thalamus). The correlation with age of the mtDNA deletion was highly significant across brain regions as previously shown. A significant decrease in the global accumulation of common deletion in subjects with SZ compared to MDD, BD, and controls was observed after correcting for age, pH, PMI, and gender. The decreases in SZ were largest in dopaminergic regions. One potential side effect of antipsychotic drugs on mitochondria is the impairment of mitochondria function, which might explain these findings. The decreased global brain mtDNA common deletion levels suggests that mitochondrial function is impaired and might be part of an overall mitochondria dysfunction signature in subjects with schizophrenia.

Evidence of allelic imbalance in the schizophrenia susceptibility gene ZNF804A in human dorsolateral prefrontal cortex.

The rs1344706, an intronic SNP within the zinc-finger protein 804A gene (ZNF804A), was identified as one of the most compelling risk SNPs for schizophrenia (SZ) and bipolar disorder (BD). It is however not clear by which molecular mechanisms ZNF804A increases disease risk. We evaluated the role of ZNF804A in SZ and BD by genotyping the originally associated rs1344706 SNP and an exonic SNP (rs12476147) located in exon four of ZNF804A in a sample of 422 SZ, 382 BD, and 507 controls from the isolated population of the Costa Rica Central Valley. We also investigated the rs1344706 SNP for allelic specific expression (ASE) imbalance in the dorsolateral prefrontal cortex (DLPFC) of 46 heterozygous postmortem brains. While no significant association between rs1344706 and SZ or BD was observed in the Costa Rica sample, we observed an increased risk of SZ for the minor allele (A) of the exonic rs12476147 SNP (p=0.026). Our ASE assay detected a significant over-expression of the rs12476147 A allele in DLPFC of rs1344706 heterozygous subjects. Interestingly, cDNA allele ratios were significantly different according to the intronic rs1344706 genotypes (p-value=0.03), with the rs1344706 A allele associated with increased ZNF804A rs12476147 A allele expression (average 1.06, p-value=0.02, for heterozygous subjects vs. genomic DNA). In conclusion, we have demonstrated a significant association of rs12476147 with SZ, and using a powerful within-subject design, an allelic expression imbalance of ZNF804A exonic SNP rs12476147 in the DLPFC. Although this data does not preclude the possibility of other functional variants in ZNF804A, it provides evidence that the rs1344706 SZ risk allele is the cis-regulatory variant directly responsible for this allelic expression imbalance in adult cortex.

Gene expression changes in the prefrontal cortex, anterior cingulate cortex and nucleus accumbens of mood disorders subjects that committed suicide.

Suicidal behaviors are frequent in mood disorders patients but only a subset of them ever complete suicide. Understanding predisposing factors for suicidal behaviors in high risk populations is of major importance for the prevention and treatment of suicidal behaviors. The objective of this project was to investigate gene expression changes associated with suicide in brains of mood disorder patients by microarrays (Affymetrix HG-U133 Plus2.0) in the dorsolateral prefrontal cortex (DLPFC: 6 Non-suicides, 15 suicides), the anterior cingulate cortex (ACC: 6NS, 9S) and the nucleus accumbens (NAcc: 8NS, 13S). ANCOVA was used to control for age, gender, pH and RNA degradation, with P ≤ 0.01 and fold change ± 1.25 as criteria for significance. Pathway analysis revealed serotonergic signaling alterations in the DLPFC and glucocorticoid signaling alterations in the ACC and NAcc. The gene with the lowest p-value in the DLPFC was the 5-HT2A gene, previously associated both with suicide and mood disorders. In the ACC 6 metallothionein genes were down-regulated in suicide (MT1E, MT1F, MT1G, MT1H, MT1X, MT2A) and three were down-regulated in the NAcc (MT1F, MT1G, MT1H). Differential expression of selected genes was confirmed by qPCR, we confirmed the 5-HT2A alterations and the global down-regulation of members of the metallothionein subfamilies MT 1 and 2 in suicide completers. MTs 1 and 2 are neuro-protective following stress and glucocorticoid stimulations, suggesting that in suicide victims neuroprotective response to stress and cortisol may be diminished. Our results thus suggest that suicide-specific expression changes in mood disorders involve both glucocorticoids regulated metallothioneins and serotonergic signaling in different regions of the brain.

Validated LC-MS/MS methods for the determination of dapagliflozin, a sodium-glucose co-transporter 2 inhibitor in normal and ZDF rat plasma.

BACKGROUND: Dapagliflozin is an inhibitor of sodium-glucose co-transporter 2 (SGLT-2) in development for the treatment of Type 2 diabetes. To support toxicology studies, LC-MS/MS methods were developed and validated for the quantitation of dapagliflozin in rat plasma. RESULTS: The assay uses solid phase extraction and LC-MS/MS analysis in negative ion electrospray ionization mode. Because dapagliflozin readily forms adducts in the presence of formic acid, the mobile phases were simple mixtures of water and acetonitrile. The assay was validated in the concentration range of 5-2000 ng/ml with good intra- and inter-day precisions and acceptable sample stability. CONCLUSION: The validated assay was successfully applied to the quantitation of dapagliflozin in plasma in support of preclinical studies in both normal and diabetic rats.

Mitochondrial variants in schizophrenia, bipolar disorder, and major depressive disorder.

BACKGROUND: Mitochondria provide most of the energy for brain cells by the process of oxidative phosphorylation. Mitochondrial abnormalities and deficiencies in oxidative phosphorylation have been reported in individuals with schizophrenia (SZ), bipolar disorder (BD), and major depressive disorder (MDD) in transcriptomic, proteomic, and metabolomic studies. Several mutations in mitochondrial DNA (mtDNA) sequence have been reported in SZ and BD patients. METHODOLOGY/PRINCIPAL FINDINGS: Dorsolateral prefrontal cortex (DLPFC) from a cohort of 77 SZ, BD, and MDD subjects and age-matched controls (C) was studied for mtDNA sequence variations and heteroplasmy levels using Affymetrix mtDNA resequencing arrays. Heteroplasmy levels by microarray were compared to levels obtained with SNaPshot and allele specific real-time PCR. This study examined the association between brain pH and mtDNA alleles. The microarray resequencing of mtDNA was 100% concordant with conventional sequencing results for 103 mtDNA variants. The rate of synonymous base pair substitutions in the coding regions of the mtDNA genome was 22% higher (p = 0.0017) in DLPFC of individuals with SZ compared to controls. The association of brain pH and super haplogroup (U, K, UK) was significant (p = 0.004) and independent of postmortem interval time. CONCLUSIONS: Focusing on haplogroup and individual susceptibility factors in psychiatric disorders by considering mtDNA variants may lead to innovative treatments to improve mitochondrial health and brain function.