Mitochondrial genomes and exceptional longevity in a Chinese population: the Rugao longevity study

Mother's Curse effects on lifespan and aging

The Mother's Curse hypothesis posits that mothers curse their sons with harmful mitochondria, because maternal mitochondrial inheritance makes selection blind to mitochondrial mutations that harm only males. As a result, mitochondrial function may be evolutionarily optimized for females. This is an attractive explanation for ubiquitous sex differences in lifespan and aging, given the prevalence of maternal mitochondrial inheritance and the established relationship between mitochondria and aging. This review outlines patterns expected under the hypothesis, and traits most likely to be affected, chiefly those that are sexually dimorphic and energy intensive. A survey of the literature shows that evidence for Mother's Curse is limited to a few taxonomic groups, with the strongest support coming from experimental crosses in Drosophila. Much of the evidence comes from studies of fertility, which is expected to be particularly vulnerable to male-harming mitochondrial mutations, but studies of lifespan and aging also show evidence of Mother's Curse effects. Despite some very compelling studies supporting the hypothesis, the evidence is quite patchy overall, with contradictory results even found for the same traits in the same taxa. Reasons for this scarcity of evidence are discussed, including nuclear compensation, factors opposing male-specific mutation load, effects of interspecific hybridization, context dependency and demographic effects. Mother's Curse effects may indeed contribute to sex differences, but the complexity of other contributing factors make Mother's Curse a poor general predictor of sex-specific lifespan and aging.

Common mtDNA variations at C5178a and A249d/T6392C/G10310A decrease the risk of severe COVID-19 in a Han Chinese population from Central China

BACKGROUND

Mitochondria have been shown to play vital roles during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and coronavirus disease 2019 (COVID-19) development. Currently, it is unclear whether mitochondrial DNA (mtDNA) variants, which define mtDNA haplogroups and determine oxidative phosphorylation performance and reactive oxygen species production, are associated with COVID-19 risk.

METHODS

A population-based case-control study was conducted to compare the distribution of mtDNA variations defining mtDNA haplogroups between healthy controls (n = 615) and COVID-19 patients (n = 536). COVID-19 patients were diagnosed based on molecular diagnostics of the viral genome by qPCR and chest X-ray or computed tomography scanning. The exclusion criteria for the healthy controls were any history of disease in the month preceding the study assessment. MtDNA variants defining mtDNA haplogroups were identified by PCR-RFLPs and HVS-I sequencing and determined based on mtDNA phylogenetic analysis using Mitomap Phylogeny. Student's t-test was used for continuous variables, and Pearson's chi-squared test or Fisher's exact test was used for categorical variables. To assess the independent effect of each mtDNA variant defining mtDNA haplogroups, multivariate logistic regression analyses were performed to calculate the odds ratios (ORs) and 95% confidence intervals (CIs) with adjustments for possible confounding factors of age, sex, smoking and diseases (including cardiopulmonary diseases, diabetes, obesity and hypertension) as determined through clinical and radiographic examinations.

RESULTS

Multivariate logistic regression analyses revealed that the most common investigated mtDNA variations (> 10% in the control population) at C5178a (in NADH dehydrogenase subunit 2 gene, ND2) and A249d (in the displacement loop region, D-loop)/T6392C (in cytochrome c oxidase I gene, CO1)/G10310A (in ND3) were associated with a reduced risk of severe COVID-19 (OR = 0.590, 95% CI 0.428-0.814, P = 0.001; and OR = 0.654, 95% CI 0.457-0.936, P = 0.020, respectively), while A4833G (ND2), A4715G (ND2), T3394C (ND1) and G5417A (ND2)/C16257a (D-loop)/C16261T (D-loop) were related to an increased risk of severe COVID-19 (OR = 2.336, 95% CI 1.179-4.608, P = 0.015; OR = 2.033, 95% CI 1.242-3.322, P = 0.005; OR = 3.040, 95% CI 1.522-6.061, P = 0.002; and OR = 2.890, 95% CI 1.199-6.993, P = 0.018, respectively).

CONCLUSIONS

This is the first study to explore the association of mtDNA variants with individual's risk of developing severe COVID-19. Based on the case-control study, we concluded that the common mtDNA variants at C5178a and A249d/T6392C/G10310A might contribute to an individual's resistance to developing severe COVID-19, whereas A4833G, A4715G, T3394C and G5417A/C16257a/C16261T might increase an individual's risk of developing severe COVID-19.

Genetic and epigenetic regulation of human aging and longevity

Here we summarize the latest data on genetic and epigenetic contributions to human aging and longevity. Whereas environmental and lifestyle factors are important at younger ages, the contribution of genetics appears more important in reaching extreme old age. Genome-wide studies have implicated ~57 gene loci in lifespan. Epigenomic changes during aging profoundly affect cellular function and stress resistance. Dysregulation of transcriptional and chromatin networks is likely a crucial component of aging. Large-scale bioinformatic analyses have revealed involvement of numerous interaction networks. As the young well-differentiated cell replicates into eventual senescence there is drift in the highly regulated chromatin marks towards an entropic middle-ground between repressed and active, such that genes that were previously inactive "leak". There is a breakdown in chromatin connectivity such that topologically associated domains and their insulators weaken, and well-defined blocks of constitutive heterochromatin give way to generalized, senescence-associated heterochromatin, foci. Together, these phenomena contribute to aging.

Genetic cartography of longevity in humans and mice: Current landscape and horizons

Aging is a complex and highly variable process. Heritability of longevity among humans and other species is low, and this finding has given rise to the idea that it may be futile to search for DNA variants that modulate aging. We argue that the problem in mapping longevity genes is mainly one of low power and the genetic and environmental complexity of aging. In this review we highlight progress made in mapping genes and molecular networks associated with longevity, paying special attention to work in mice and humans. We summarize 40 years of linkage studies using murine cohorts and 15 years of studies in human populations that have exploited candidate gene and genome-wide association methods. A small but growing number of gene variants contribute to known longevity mechanisms, but a much larger set have unknown functions. We outline these and other challenges and suggest some possible solutions, including more intense collaboration between research communities that use model organisms and human cohorts. Once hundreds of gene variants have been linked to differences in longevity in mammals, it will become feasible to systematically explore gene-by-environmental interactions, dissect mechanisms with more assurance, and evaluate the roles of epistasis and epigenetics in aging. A deeper understanding of complex networks-genetic, cellular, physiological, and social-should position us well to improve healthspan.

Bagging Nearest-Neighbor Prediction independence Test: an efficient method for nonlinear dependence of two continuous variables

Testing dependence/correlation of two variables is one of the fundamental tasks in statistics. In this work, we proposed an efficient method for nonlinear dependence of two continuous variables (X and Y). We addressed this research question by using BNNPT (Bagging Nearest-Neighbor Prediction independence Test, software available at https://sourceforge.net/projects/bnnpt/). In the BNNPT framework, we first used the value of X to construct a bagging neighborhood structure. We then obtained the out of bag estimator of Y based on the bagging neighborhood structure. The square error was calculated to measure how well Y is predicted by X. Finally, a permutation test was applied to determine the significance of the observed square error. To evaluate the strength of BNNPT compared to seven other methods, we performed extensive simulations to explore the relationship between various methods and compared the false positive rates and statistical power using both simulated and real datasets (Rugao longevity cohort mitochondrial DNA haplogroups and kidney cancer RNA-seq datasets). We concluded that BNNPT is an efficient computational approach to test nonlinear correlation in real world applications.

MtDNA analysis reveals enriched pathogenic mutations in Tibetan highlanders

Tibetan highlanders, including Tibetans, Monpas, Lhobas, Dengs and Sherpas, are considered highly adaptive to severe hypoxic environments. Mitochondrial DNA (mtDNA) might be important in hypoxia adaptation given its role in coding core subunits of oxidative phosphorylation. In this study, we employed 549 complete highlander mtDNA sequences (including 432 random samples) to obtain a comprehensive view of highlander mtDNA profile. In the phylogeny of a total of 36,914 sequences, we identified 21 major haplogroups representing founding events of highlanders, most of which were coalesced in 10 kya. Through founder analysis, we proposed a three-phase model of colonizing the plateau, i.e., pre-LGM Time (30 kya, 4.68%), post-LGM Paleolithic Time (16.8 kya, 29.31%) and Neolithic Time (after 8 kya, 66.01% in total). We observed that pathogenic mutations occurred far more frequently in 22 highlander-specific lineages (five lineages carrying two pathogenic mutations and six carrying one) than in the 6,857 haplogroups of all the 36,914 sequences (P = 4.87 × 10⁻⁸). Furthermore, the number of possible pathogenic mutations carried by highlanders (in average 3.18 ± 1.27) were significantly higher than that in controls (2.82 ± 1.40) (P = 1.89 × 10⁻⁴). Considering that function-altering and pathogenic mutations are enriched in highlanders, we therefore hypothesize that they may have played a role in hypoxia adaptation.

Compatibility between mitochondrial and nuclear genomes correlates with the quantitative trait of lifespan in Caenorhabditis elegans

Mutations in mitochondrial genome have epistatic effects on organisms depending on the nuclear background, but a role for the compatibility of mitochondrial-nuclear genomes (mit-n) in the quantitative nature of a complex trait remains unexplored. We studied a panel of recombinant inbred advanced intercrossed lines (RIAILs) of C. elegans that were established from a cross between the N2 and HW strains. We determined the HW nuclear genome content and the mitochondrial type (HW or N2) of each RIAIL strain. We found that the degree of mit-n compatibility was correlated with the lifespans but not the foraging behaviors of RIAILs. Several known aging-associated QTLs individually showed no relationship with mitotypes but collectively a weak trend consistent with a role in mit-n compatibility. By association mapping, we identified 293 SNPs that showed linkage with lifespan and a relationship with mitotypes consistent with a role in mit-n compatibility. We further found an association between mit-n compatibility and several functional characteristics of mitochondria as well as the expressions of genes involved in the respiratory oxidation pathway. The results provide the first evidence implicating mit-n compatibility in the quantitative nature of a complex trait, and may be informative to certain evolutionary puzzles on hybrids.