Long-term outcomes by lobular versus ductal histology in four NSABP adjuvant breast cancer trials.

We evaluated differences in long-term outcomes of invasive lobular carcinoma (ILC) vs breast cancers of no special type (NST) treated with anthracycline-based adjuvant chemotherapy using 4 National Surgical Adjuvant Breast and Bowel Project (NSABP) randomized phase III trials (B-22, B-25, B-28, B-30). Our cohort included 11,251 patients with NST and 1,231 with ILC. Patients with ILC were older, had larger and more often estrogen receptor-positive tumors, and more positive lymph nodes. During early follow-up (0-5 years), patients with ILC had fewer recurrences (HR: 0.797; 95% confidence interval [CI] 0.685-0.929) and deaths (HR: 0.756; 95% CI 0.623-0.917). After 5 years patients with ILC had more recurrences (HR: 1.30; 95% CI 1.085-1.558) and deaths (HR: 1.044; 95% CI 0.898-1.214). Conditional probability analysis showed significant interactions between time-period and histologic type for recurrences (p < .001) and deaths (p < .001). Patients with ILC have elevated risk of late recurrence and death compared to patients with NST.

Hypermetabolic state is associated with circadian rhythm disruption in mouse and human cancer cells.

Crosstalk between metabolism and circadian rhythms is a fundamental building block of multicellular life, and disruption of this reciprocal communication could be relevant to disease. Here, we investigated whether maintenance of circadian rhythms depends on specific metabolic pathways, particularly in the context of cancer. We found that in adult mouse fibroblasts, ATP levels were a major contributor to signal from a clock gene luciferase reporter, although not necessarily to the strength of circadian cycling. In contrast, we identified significant metabolic control of circadian function across a series of pancreatic adenocarcinoma cell lines. Metabolic profiling of congenic tumor cell clones revealed substantial diversity among these lines that we used to identify clones to generate circadian reporter lines. We observed diverse circadian profiles among these lines that varied with their metabolic phenotype: The most hypometabolic line [exhibiting low levels of oxidative phosphorylation (OxPhos) and glycolysis] had the strongest rhythms, while the most hypermetabolic line had the weakest rhythms. Pharmacological enhancement of OxPhos decreased the amplitude of circadian oscillation in a subset of tumor cell lines. Strikingly, inhibition of OxPhos enhanced circadian rhythms only in the tumor cell line in which glycolysis was also low, thereby establishing a hypometabolic state. We further analyzed metabolic and circadian phenotypes across a panel of human patient-derived melanoma cell lines and observed a significant negative association between metabolic activity and circadian cycling strength. Together, these findings suggest that metabolic heterogeneity in cancer directly contributes to circadian function and that high levels of glycolysis or OxPhos independently disrupt circadian rhythms in these cells.

Shared and unique 3D genomic features of substance use disorders across multiple cell types.

Recent genome-wide association studies (GWAS) have revealed shared genetic components among alcohol, opioid, tobacco and cannabis use disorders. However, the extent of the underlying shared causal variants and effector genes, along with their cellular context, remain unclear. We leveraged our existing 3D genomic datasets comprising high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq and RNA-seq across >50 diverse human cell types to focus on genomic regions that coincide with GWAS loci. Using stratified LD regression, we determined the proportion of genomewide SNP heritability attributable to the features assayed across our cell types by integrating recent GWAS summary statistics for the relevant traits: alcohol use disorder (AUD), tobacco use disorder (TUD), opioid use disorder (OUD) and cannabis use disorder (CanUD). Statistically significant enrichments (P<0.05) were observed in 14 specific cell types, with heritability reaching 9.2-fold for iPSC-derived cortical neurons and neural progenitors, confirming that they are crucial cell types for further functional exploration. Additionally, several pancreatic cell types, notably pancreatic beta cells, showed enrichment for TUD, with heritability enrichments up to 4.8-fold, suggesting genomic overlap with metabolic processes. Further investigation revealed significant positive genetic correlations between T2D with both TUD and CanUD (FDR<0.05) and a significant negative genetic correlation with AUD. Interestingly, after partitioning the heritability for each cell type's cis-regulatory elements, the correlation between T2D and TUD for pancreatic beta cells was greater (r=0.2) than the global genetic correlation value. Our study provides new genomic insights into substance use disorders and implicates cell types where functional follow-up studies could reveal causal variant-gene mechanisms underpinning these disorders.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2.

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

22q11.2 Deletion-Associated Blood-Brain Barrier Permeability Potentiates Systemic Capillary Leak Syndrome Neurologic Features.

We present a case study of a young male with a history of 22q11.2 deletion syndrome (22qDS), diagnosed with systemic capillary leak syndrome (SCLS) who presented with acute onset of diffuse anasarca and sub-comatose obtundation. We hypothesized that his co-presentation of neurological sequelae might be due to blood-brain barrier (BBB) susceptibility conferred by the 22q11.2 deletion, a phenotype that we have previously identified in 22qDS. Using pre- and post-intravenous immunoglobulins (IVIG) patient serum, we studied circulating biomarkers of inflammation and assessed the potential susceptibility of the 22qDS BBB. We employed in vitro cultures of differentiated BBB-like endothelial cells derived from a 22qDS patient and a healthy control. We found evidence of peripheral inflammation and increased serum lipopolysaccharide (LPS) alongside endothelial cells in circulation. We report that the patient's serum significantly impairs barrier function of the 22qDS BBB compared to control. Only two other cases of pediatric SCLS with neurologic symptoms have been reported, and genetic risk factors have been suggested in both instances. As the third case to be reported, our findings are consistent with the hypothesis that genetic susceptibility of the BBB conferred by genes such as claudin-5 deleted in the 22q11.2 region promoted neurologic involvement during SCLS in this patient.

Targeted blockade of aberrant sodium current in a stem cell-derived neuron model of SCN3A encephalopathy.

Missense variants in SCN3A encoding the voltage-gated sodium (Na+) channel α subunit Nav1.3 are associated with SCN3A-related neurodevelopmental disorder (SCN3A-NDD), a spectrum of disease that includes epilepsy and malformation of cortical development. How genetic variation in SCN3A leads to pathology remains unclear, as prior electrophysiological work on disease-associated variants has been performed exclusively in heterologous cell systems. To further investigate the mechanisms of SCN3A-NDD pathogenesis, we used CRISPR/Cas9 gene editing to modify a control human induced pluripotent stem cell (iPSC) line to express the recurrent de novo missense variant SCN3A c.2624T>C (p.Ile875Thr). With the established Ngn2 rapid induction protocol, we generated glutamatergic forebrain-like neurons (iNeurons), which we showed to express SCN3A mRNA and Nav1.3-mediated Na+ currents. We performed detailed whole-cell patch clamp recordings to determine the effect of the SCN3A-p.Ile875Thr variant on endogenous Na+ currents in, and intrinsic excitability of, human neurons. Compared to control iNeurons, variant-expressing iNeurons exhibit markedly increased slowly-inactivating/persistent Na+ current, abnormal firing patterns with paroxysmal bursting and plateau-like potentials with action potential failure, and a hyperpolarized voltage threshold for action potential generation. We then validated these findings using a separate iPSC line generated from a patient harbouring the SCN3A-p.Ile875Thr variant compared to a corresponding CRISPR-corrected isogenic control line. Finally, we found that application of the Nav1.3-selective blocker ICA-121431 normalizes action potential threshold and aberrant firing patterns in SCN3A-p.Ile1875Thr iNeurons; in contrast, consistent with action as a Na+ channel blocker, ICA-121431 decreases excitability of control iNeurons. Our findings demonstrate that iNeurons can model the effects of genetic variation in SCN3A yet reveal a complex relationship between gain-of-function at the level of the ion channel versus impact on neuronal excitability. Given the transient expression of SCN3A in the developing human nervous system, selective blockade or suppression of Nav1.3-containing Na+ channels could represent a therapeutic approach towards SCN3A-NDD.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity.

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

CHEX-seq detects single-cell genomic single-stranded DNA with catalytical potential.

Genomic DNA (gDNA) undergoes structural interconversion between single- and double-stranded states during transcription, DNA repair and replication, which is critical for cellular homeostasis. We describe "CHEX-seq" which identifies the single-stranded DNA (ssDNA) in situ in individual cells. CHEX-seq uses 3'-terminal blocked, light-activatable probes to prime the copying of ssDNA into complementary DNA that is sequenced, thereby reporting the genome-wide single-stranded chromatin landscape. CHEX-seq is benchmarked in human K562 cells, and its utilities are demonstrated in cultures of mouse and human brain cells as well as immunostained spatially localized neurons in brain sections. The amount of ssDNA is dynamically regulated in response to perturbation. CHEX-seq also identifies single-stranded regions of mitochondrial DNA in single cells. Surprisingly, CHEX-seq identifies single-stranded loci in mouse and human gDNA that catalyze porphyrin metalation in vitro, suggesting a catalytic activity for genomic ssDNA. We posit that endogenous DNA enzymatic activity is a function of genomic ssDNA.

Hypermetabolic state is associated with circadian rhythm disruption in mouse and human cancer cells.

Crosstalk between cellular metabolism and circadian rhythms is a fundamental building block of multicellular life, and disruption of this reciprocal communication could be relevant to degenerative disease, including cancer. Here, we investigated whether maintenance of circadian rhythms depends upon specific metabolic pathways, particularly in the context of cancer. We found that in adult mouse fibroblasts, ATP levels were a major contributor to overall levels of a clock gene luciferase reporter, although not necessarily to the strength of circadian cycling. In contrast, we identified significant metabolic control of circadian function in an in vitro mouse model of pancreatic adenocarcinoma. Metabolic profiling of a library of congenic tumor cell clones revealed significant differences in levels of lactate, pyruvate, ATP, and other crucial metabolites that we used to identify candidate clones with which to generate circadian reporter lines. Despite the shared genetic background of the clones, we observed diverse circadian profiles among these lines that varied with their metabolic phenotype: the most hypometabolic line had the strongest circadian rhythms while the most hypermetabolic line had the weakest rhythms. Treatment of these tumor cell lines with bezafibrate, a peroxisome proliferator-activated receptor (PPAR) agonist shown to increase OxPhos, decreased the amplitude of circadian oscillation in a subset of tumor cell lines. Strikingly, treatment with the Complex I antagonist rotenone enhanced circadian rhythms only in the tumor cell line in which glycolysis was also low, thereby establishing a hypometabolic state. We further analyzed metabolic and circadian phenotypes across a panel of human patient-derived melanoma cell lines and observed a significant negative association between metabolic activity and circadian cycling strength. Together, these findings suggest that metabolic heterogeneity in cancer directly contributes to circadian function, and that high levels of glycolysis or OxPhos independently disrupt circadian rhythms in these cells.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2.

The ch12q13 obesity locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via an influence on cis-regulation within the genomic region. We implicated rs7132908 as a putative causal variant at this locus leveraging a combination of our inhouse 3D genomic data, public domain datasets, and several computational approaches. Using a luciferase reporter assay in human primary astrocytes, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. Motivated by this finding, we went on to generate isogenic human embryonic stem cell lines homozygous for either rs7132908 allele with CRISPR-Cas9 homology-directed repair to assess changes in gene expression due to genotype and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. We observed that the rs7132908 obesity risk allele influenced the expression of FAIM2 along with other genes, decreased the proportion of neurons produced during differentiation, up-regulated cell death gene sets, and conversely down-regulated neuron differentiation gene sets. We have therefore functionally validated rs7132908 as a causal obesity variant which temporally regulates nearby effector genes at the ch12q13 locus and influences neurodevelopment and survival.

Immune status of the murine 22q11.2 deletion syndrome model.

Mice modeling the hemizygous deletion of chromosome 22q11.2 (22qMc) have been utilized to address various clinical phenotypes associated with the disease, including cardiac malformations, altered neural circuitry, and behavioral deficits. Yet, the status of the T cell compartment, an important clinical concern among 22q11.2 deletion syndrome (22qDS) patients, has not been addressed. While infancy and early childhood in 22qDS are associated with deficient T cell numbers and thymic hypoplasia, which can be severe in a small subset of patients, studies suggest normalization of the T cell counts by adulthood. We found that adult 22qMc do not exhibit thymic hypoplasia or altered thymic T cell development. Our findings that immune cell counts and inflammatory T cell activation are unaffected in 22qMc lend support to the hypothesis that human 22qDS immunodeficiencies are secondary to thymic hypoplasia, rather than intrinsic effects due to the deletion. Furthermore, the 22q11.2 deletion does not impact the differentiation capacity of T cells, nor their activity and response during inflammatory activation. Thus, 22qMc reflects the T cell compartment in adult 22qDS patients, and our findings suggest that 22qMc may serve as a novel model to address experimental and translational aspects of immunity in 22qDS.

Low-Dose Augmentation With Buprenorphine for Treatment-Resistant Depression: A Multisite Randomized Controlled Trial With Multimodal Assessment of Target Engagement.

Background: The experimental therapeutics approach that combines a placebo-controlled clinical trial with translational neuroscience methods can provide a better understanding of both the clinical and physiological effects of pharmacotherapy. We aimed to test the efficacy and tolerability of low-dose augmentation with buprenorphine (BPN) for treatment-resistant depression, combined with multimodal assessment of target engagement. Methods: In this multisite randomized clinical trial, 85 participants ≥50 years of age with a major depressive episode that had not responded to venlafaxine extended release were randomized to augmentation with BPN or placebo for 8 weeks. The primary outcome measure was the Montgomery-Åsberg Depression Rating Scale. In addition, three linked experiments were conducted to test target engagement: 1) functional magnetic resonance imaging using the monetary incentive delay task, 2) brain positron emission tomography of healthy participants using a novel kappa opioid receptor antagonist tracer [11C]LY2795050, and 3) transcranial magnetic stimulation measure of cortical transmission after daily BPN administration. Results: The mean ± SD dosage of BPN was 0.59 ± 0.33 mg/day. There were no significant differences between the BPN and placebo groups in Montgomery-Åsberg Depression Rating Scale changes over time or adverse effects. BPN administration had minimal effects on functional magnetic resonance imaging blood oxygen level-dependent responses in regions involved in reward anticipation and response, no significant displacement of kappa opioid receptor radioligand in positron emission tomography imaging, and no significant changes in transcranial magnetic stimulation measures of inhibitory and excitatory cortical transmission. Conclusions: Our findings suggest a lack of clinical effect of low-dose BPN augmentation and lack of target engagement with this dosage and physiological probes.

Time-varying coefficient cumulative gap time models for intensive longitudinal ecological momentary assessment data with missingness.

Ecological momentary assessment (EMA) studies investigate intensive repeated observations of the current behavior and experiences of subjects in real time. In particular, such studies aim to minimize recall bias and maximize ecological validity, thereby strengthening the investigation and inference of microprocesses that influence behavior in real-world contexts by gathering intensive information on the temporal patterning of behavior of study subjects. Throughout this paper, we focus on the data analysis of an EMA study that examined behavior of intermittent smokers (ITS). Specifically, we sought to explore the pattern of clustered smoking behavior of ITS, or smoking 'bouts', as well as the covariates that predict such smoking behavior. To do this, in this paper we introduce a framework for characterizing the temporal behavior of ITS via the functions of event gap time to distinguish the smoking bouts. We used the time-varying coefficient models for the cumulative log gap time and to characterize the temporal patterns of smoking behavior, while simultaneously adjusting for behavioral covariates, and incorporated the inverse probability weighting into the models to accommodate missing data. Simulation studies showed that irrespective of whether missing by design or missing at random, the model was able to reliably determine prespecified time-varying functional forms of a given covariate coefficient, provided the the within-subject level was small.

Racial disparity in distant recurrence-free survival in patients with localized breast cancer: A pooled analysis of National Surgical Adjuvant Breast and Bowel Project trials.

BACKGROUND: Black race is associated with worse outcome in patients with breast cancer. The distant relapse-free survival (DRFS) between Black and White women with localized breast cancer who participated in National Cancer Institute-sponsored clinical trial was evaluated. METHODS: Pooled data were analyzed from 8 National Surgical Adjuvant Breast and Bowel Project (NSABP) trials including 9702 women with localized breast cancer treated with adjuvant chemotherapy (AC, n = 7485) or neoadjuvant chemotherapy (NAC, n = 2217), who self-reported as Black (n = 1070) or White (n = 8632) race. The association between race and DRFS was analyzed using log-rank tests and multivariate Cox regression. RESULTS: After adjustment for covariates including age, tumor size, nodal status, body mass index and taxane use, and treatment (AC vs NAC), Black race was associated with an inferior DRFS in estrogen receptor-positive (ER+; hazard ratio [HR], 1.24; 95% CI, 1.05-1.46; P = .01), but not in ER- disease (HR, 0.97; 95% CI, 0.83-1.14; P = .73), and significant interaction between race and ER status was observed (P = .03). There was no racial disparity in DRFS among patients with pathologic complete response (pCR) (log-rank P = .8). For patients without pCR, Black race was associated with worse DRFS in ER+ (HR, 1.67; 95% CI, 1.14-2.45; P = .01), but not in ER- disease (HR, 0.91; 95% CI, 0.65-1.28; P = .59). CONCLUSIONS: Black race was associated with significantly inferior DRFS in ER+ localized breast cancer treated with AC or NAC, but not in ER- disease. In the NAC group, racial disparity was also observed in patients with residual ER+ breast cancer at surgery, but not in those who had pCR. LAY SUMMARY: Black women with breast cancer have worse outcomes compared with White women. We investigated if this held true in the context of clinical trials that provide controlled treatment setting. Black women with cancer expressing estrogen receptors (ERs) had worse outcome than White women. If breast cancers did not express ERs, there was no racial disparity in outcome. We also observed racial disparity in women who received chemotherapy before their cancer was removed, but only if they had cancer expressing ERs and residual disease on completion of treatment. If the cancer disappeared with presurgical chemotherapy, there was no racial disparity.

Deficits in Seizure Threshold and Other Behaviors in Adult Mice without Gross Neuroanatomic Injury after Late Gestation Transient Prenatal Hypoxia.

Intrauterine hypoxia is a common cause of brain injury in children resulting in a broad spectrum of long-term neurodevelopmental sequela, including life-long disabilities that can occur even in the absence of severe neuroanatomic damage. Postnatal hypoxia-ischemia rodent models are commonly used to understand the effects of ischemia and transient hypoxia on the developing brain. Postnatal models, however, have some limitations. First, they do not test the impact of placental pathologies on outcomes from hypoxia. Second, they primarily recapitulate severe injury because they provoke substantial cell death, which is not seen in children with mild hypoxic injury. Lastly, they do not model preterm hypoxic injury. Prenatal models of hypoxia in mice may allow us to address some of these limitations to expand our understanding of developmental brain injury. The published rodent models of prenatal hypoxia employ multiple days of hypoxic exposure or complicated surgical procedures, making these models challenging to perform consistently in mice. Furthermore, large animal models suggest that transient prenatal hypoxia without ischemia is sufficient to lead to significant functional impairment to the developing brain. However, these large animal studies are resource-intensive and not readily amenable to mechanistic molecular studies. Therefore, here we characterized the effect of late gestation (embryonic day 17.5) transient prenatal hypoxia (5% inspired oxygen) on long-term anatomical and neurodevelopmental outcomes in mice. Late gestation transient prenatal hypoxia increased hypoxia-inducible factor 1 alpha protein levels (a marker of hypoxic exposure) in the fetal brain. Hypoxia exposure predisposed animals to decreased weight at postnatal day 2, which normalized by day 8. However, hypoxia did not affect gestational age at birth, litter size at birth, or pup survival. No differences in fetal brain cell death or long-term gray or white matter changes resulted from hypoxia. Animals exposed to prenatal hypoxia did have several long-term functional consequences, including sex-dichotomous changes. Hypoxia exposure was associated with a decreased seizure threshold and abnormalities in hindlimb strength and repetitive behaviors in males and females. Males exposed to hypoxia had increased anxiety-related deficits, whereas females had deficits in social interaction. Neither sex developed any motor or visual learning deficits. This study demonstrates that late gestation transient prenatal hypoxia in mice is a simple, clinically relevant paradigm for studying putative environmental and genetic modulators of the long-term effects of hypoxia on the developing brain.

A Short History of Bernard Fisher's Contributions to Randomized Clinical Trials.

Dr. Bernard Fisher (1918-2019) was an early proponent of evidence-based medicine using the mechanism of prospective, multicenter, randomized clinical trials to test biological and clinical hypotheses. In this article, I trace how his early scientific work in striving to understand the nature of cancer metastasis through animal experiments led to a new, testable, clinical hypothesis: that surgery to remove only the tumor and a small amount of tissue around it was as effective as the more disfiguring operations that were then the standard treatment. Fisher's work with the National Surgical Adjuvant Breast and Bowel Project (NSABP) using large, randomized clinical trials to demonstrate the veracity of this hypothesis led to a new paradigm in which the emphasis was placed on how systemic therapies used at an early stage of disease could effectively eradicate breast cancer for many patients. This new therapeutic approach led to the successful development of new treatments, many of which are widely used today. Ultimately, the new paradigm led to successfully preventing breast cancer in women who were at high risk for the disease but who had not yet been diagnosed with the disease. Throughout his entire career, Fisher championed the use of large prospective, randomized clinical trials despite criticism from many in the medical community who strongly criticized his use of randomization as a mechanism for testing clinical hypotheses. The approach he and the NSABP employed is still considered to be the highest standard of evidence in conducting clinical studies.

MitoScape: A big-data, machine-learning platform for obtaining mitochondrial DNA from next-generation sequencing data.

The growing number of next-generation sequencing (NGS) data presents a unique opportunity to study the combined impact of mitochondrial and nuclear-encoded genetic variation in complex disease. Mitochondrial DNA variants and in particular, heteroplasmic variants, are critical for determining human disease severity. While there are approaches for obtaining mitochondrial DNA variants from NGS data, these software do not account for the unique characteristics of mitochondrial genetics and can be inaccurate even for homoplasmic variants. We introduce MitoScape, a novel, big-data, software for extracting mitochondrial DNA sequences from NGS. MitoScape adopts a novel departure from other algorithms by using machine learning to model the unique characteristics of mitochondrial genetics. We also employ a novel approach of using rho-zero (mitochondrial DNA-depleted) data to model nuclear-encoded mitochondrial sequences. We showed that MitoScape produces accurate heteroplasmy estimates using gold-standard mitochondrial DNA data. We provide a comprehensive comparison of the most common tools for obtaining mtDNA variants from NGS and showed that MitoScape had superior performance to compared tools in every statistically category we compared, including false positives and false negatives. By applying MitoScape to common disease examples, we illustrate how MitoScape facilitates important heteroplasmy-disease association discoveries by expanding upon a reported association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men (adjusted p-value = 0.003). The improved accuracy of mitochondrial DNA variants produced by MitoScape will be instrumental in diagnosing disease in the context of personalized medicine and clinical diagnostics.