Significance: There are hundreds to thousands of copies of mitochondrial DNA (mtDNA) in each human cell in contrast to only two copies of nuclear DNA. High-frequency pathogenic mtDNA mutations have been found in patients with classic mitochondrial diseases, premature aging, cancers, and neurodegenerative diseases. In this study we investigated the distribution of heteroplasmic mutations, their pathogenic potential, and their underlying evolutionary forces using genome sequence data from the 1000 Genomes Project. Our results demonstrated the prevalence of low-frequency high-pathogenic-potential mtDNA mutations in healthy human individuals. These deleterious mtDNA mutations, when reaching high frequency, could provide a likely source of mitochondrial dysfunction. Managing the expansion of deleterious mtDNA mutations could be a promising means of preventing disease progression.

The accumulation of mitochondrial DNA (mtDNA) mutations, and the reduction of mtDNA copy number, both disrupt mitochondrial energetics, and may contribute to aging and age-associated phenotypes. However, there are few genetic and epidemiological studies on the spectra of blood mtDNA heteroplasmies, and the distribution of mtDNA copy numbers in different age groups and their impact on age-related phenotypes. In this work, we used whole-genome sequencing data of isolated peripheral blood mononuclear cells (PBMCs) from the UK10K project to investigate in parallel mtDNA heteroplasmy and copy number in 1511 women, between 17 and 85 years old, recruited in the TwinsUK cohorts. We report a high prevalence of pathogenic mtDNA heteroplasmies in this population. We also find an increase in mtDNA heteroplasmies with age (β = 0.011, P = 5.77e-6), and showed that, on average, individuals aged 70-years or older had 58.5% more mtDNA heteroplasmies than those under 40-years old. Conversely, mtDNA copy number decreased by an average of 0.4 copies per year (β = −0.395, P = 0.0097). Multiple regression analyses also showed that age had independent effects on mtDNA copy number decrease and heteroplasmy accumulation. Finally, mtDNA copy number was positively associated with serum bicarbonate level (P = 4.46e-5), and inversely correlated with white blood cell count (P = 0.0006). Moreover, the aggregated heteroplasmy load was associated with blood apolipoprotein B level (P = 1.33e-5), linking the accumulation of mtDNA mutations to age-related physiological markers. Our population-based study indicates that both mtDNA quality and quantity are influenced by age. An open question for the future is whether interventions that would contribute to maintain optimal mtDNA copy number and prevent the expansion of heteroplasmy could promote healthy aging.

Granzyme A from killer lymphocytes cleaves gasdermin B (GSDMB) and triggers pyroptosis in targeted human tumor cells, eliciting antitumor immunity. However, GSDMB has a controversial role in pyroptosis and has been linked to both anti- and protumor functions. Here, we found that GSDMB splicing variants are functionally distinct. Cleaved N-terminal (NT) fragments of GSDMB isoforms 3 and 4 caused pyroptosis, but isoforms 1, 2, and 5 did not. The nonfunctional isoforms have a deleted or modified exon 6 and therefore lack a stable belt motif. The belt likely contributes to the insertion of oligomeric GSDMB-NTs into the membrane. Consistently, noncytotoxic GSDMB-NTs blocked pyroptosis caused by cytotoxic GSDMB-NTs in a dominant-negative manner. Upon natural killer (NK) cell attack, GSDMB3-expressing cells died by pyroptosis, whereas GSDMB4-expressing cells died by mixed pyroptosis and apoptosis, and GSDMB1/2-expressing cells died only by apoptosis. GSDMB4 partially resisted NK cell-triggered cleavage, suggesting that only GSDMB3 is fully functional. GSDMB1-3 were the most abundant isoforms in the tested tumor cell lines and were similarly induced by interferon-γ and the chemotherapy drug methotrexate. Expression of cytotoxic GSDMB3/4 isoforms, but not GSDMB1/2 isoforms that are frequently up-regulated in tumors, was associated with better outcomes in bladder and cervical cancers, suggesting that GSDMB3/4-mediated pyroptosis was protective in those tumors. Our study indicates that tumors may block and evade killer cell-triggered pyroptosis by generating noncytotoxic GSDMB isoforms. Therefore, therapeutics that favor the production of cytotoxic GSDMB isoforms by alternative splicing may improve antitumor immunity.

Fatty acid desaturase (FADS) genes encode rate-limiting enzymes for the biosynthesis of omega-6 and omega-3 long-chain polyunsaturated fatty acids (LCPUFAs). This biosynthesis is essential for individuals subsisting on LCPUFA-poor diets (for example, plant-based). Positive selection on FADS genes has been reported in multiple populations, but its cause and pattern in Europeans remain unknown. Here we demonstrate, using ancient and modern DNA, that positive selection acted on the same FADS variants both before and after the advent of farming in Europe, but on opposite (that is, alternative) alleles. Recent selection in farmers also varied geographically, with the strongest signal in southern Europe. These varying selection patterns concur with anthropological evidence of varying diets, and with the association of farming-adaptive alleles with higher FADS1 expression and thus enhanced LCPUFA biosynthesis. Genome-wide association studies reveal that farming-adaptive alleles not only increase LCPUFAs, but also affect other lipid levels and protect against several inflammatory diseases.