High mitochondrial DNA stability in B-cell chronic lymphocytic leukemia

Reactive Oxygen Species and Metabolism in Leukemia: A Dangerous Liaison

Reactive oxygen species (ROS), previously considered toxic by-products of aerobic metabolism, are increasingly recognized as regulators of cellular signaling. Keeping ROS levels low is essential to safeguard the self-renewal capacity of hematopoietic stem cells (HSC). HSC reside in a hypoxic environment and have been shown to be highly dependent on the glycolytic pathway to meet their energy requirements. However, when the differentiation machinery is activated, there is an essential enhancement of ROS together with a metabolic shift toward oxidative metabolism. Initiating and sustaining leukemia depend on the activity of leukemic stem cells (LSC). LSC also show low ROS levels, but unlike HSC, LSC rely on oxygen to meet their metabolic energetic requirements through mitochondrial respiration. In contrast, leukemic blasts show high ROS levels and great metabolic plasticity, both of which seem to sustain their invasiveness. Oxidative stress and metabolism rewiring are recognized as hallmarks of cancer that are intimately intermingled. Here we present a detailed overview of these two features, sustained at different levels, that support a two-way relationship in leukemia. Modifying ROS levels and targeting metabolism are interesting therapeutic approaches. Therefore, we provide the most recent evidence on the modulation of oxidative stress and metabolism as a suitable anti-leukemic approach.

Reactive oxygen species in haematopoiesis: leukaemic cells take a walk on the wild side

Oxidative stress is related to ageing and degenerative diseases, including cancer. However, a moderate amount of reactive oxygen species (ROS) is required for the regulation of cellular signalling and gene expression. A low level of ROS is important for maintaining quiescence and the differentiation potential of haematopoietic stem cells (HSCs), whereas the level of ROS increases during haematopoietic differentiation; thus, suggesting the importance of redox signalling in haematopoiesis. Here, we will analyse the importance of ROS for haematopoiesis and include evidence showing that cells from leukaemia patients live under oxidative stress. The potential sources of ROS will be described. Finally, the level of oxidative stress in leukaemic cells can also be harnessed for therapeutic purposes. In this regard, the reliance of front-line anti-leukaemia chemotherapeutics on increased levels of ROS for their mechanism of action, as well as the active search for novel compounds that modulate the redox state of leukaemic cells, will be analysed.

A Method for Next-Generation Sequencing of Paired Diagnostic and Remission Samples to Detect Mitochondrial DNA Mutations Associated with Leukemia

Somatic mitochondrial DNA (mtDNA) mutations have been identified in many human cancers, including leukemia. To identify somatic mutations, it is necessary to have a control tissue from the same individual for comparison. When patients with leukemia achieve remission, the remission peripheral blood may be a suitable and easily accessible control tissue, but this approach has not previously been applied to the study of mtDNA mutations. We have developed and validated a next-generation sequencing approach for the identification of leukemia-associated mtDNA mutations in 26 chronic myeloid leukemia patients at diagnosis using either nonhematopoietic or remission blood samples as the control. The entire mt genome was amplified by long-range PCR and sequenced using Illumina technology. Variant caller software was used to detect mtDNA somatic mutations, and an empirically determined threshold of 2% was applied to minimize false-positive results because of sequencing errors. Mutations were called against both nonhematopoietic and remission controls: the overall concordance between the two approaches was 81% (73/90 mutations). Some discordant results were because of the presence of somatic mutations in remission samples, because of either minimal residual disease or nonleukemic hematopoietic clones. This method could be applied to study somatic mtDNA mutations in leukemia patients who achieve minimal residual disease, and in patients with nonhematopoietic cancers who have a matched uninvolved tissue available.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

Metabolism pathways in chronic lymphocytic leukemia

Alterations in chronic lymphocytic leukemia (CLL) cell metabolism have been studied by several investigators. Unlike normal B lymphocytes or other leukemia cells, CLL cells, like adipocytes, store lipids and utilize free fatty acids (FFA) to produce chemical energy. None of the recently identified mutations in CLL directly affects metabolic pathways, suggesting that genetic alterations do not directly contribute to CLL cells' metabolic reprogramming. Conversely, recent data suggest that activation of STAT3 or downregulation of microRNA-125 levels plays a crucial role in the utilization of FFA to meet the CLL cells' metabolic needs. STAT3, known to be constitutively activated in CLL, increases the levels of lipoprotein lipase (LPL) that mediates lipoprotein uptake and shifts the CLL cells' metabolism towards utilization of FFA. Herein, we review the evidence for altered lipid metabolism, increased mitochondrial activity and formation of reactive oxygen species (ROS) in CLL cells, and discuss the possible therapeutic strategies to inhibit lipid metabolism pathways in patient with CLL.

The complete mitogenome of a 500-year-old Inca child mummy

In 1985, a frozen mummy was found in Cerro Aconcagua (Argentina). Archaeological studies identified the mummy as a seven-year-old Inca sacrifice victim who lived >500 years ago, at the time of the expansion of the Inca Empire towards the southern cone. The sequence of its entire mitogenome was obtained. After querying a large worldwide database of mitogenomes (>28,000) we found that the Inca haplotype belonged to a branch of haplogroup C1b (C1bi) that has not yet been identified in modern Native Americans. The expansion of C1b into the Americas, as estimated using 203 C1b mitogenomes, dates to the initial Paleoindian settlements (~18.3 thousand years ago [kya]); however, its internal variation differs between Mesoamerica and South America. By querying large databases of control region haplotypes (>150,000), we found only a few C1bi members in Peru and Bolivia (e.g. Aymaras), including one haplotype retrieved from ancient DNA of an individual belonging to the Wari Empire (Peruvian Andes). Overall, the results suggest that the profile of the mummy represents a very rare sub-clade that arose 14.3 (5–23.6) kya and could have been more frequent in the past. A Peruvian Inca origin for present-day C1bi haplotypes would satisfy both the genetic and paleo-anthropological findings.

Comprehensive Analysis of Pan-African Mitochondrial DNA Variation Provides New Insights into Continental Variation and Demography

Africa is the cradle of all human beings, and although it has been the focus of a number of genetic studies, there are many questions that remain unresolved. We have performed one of the largest and most comprehensive meta-analyses of mitochondrial DNA (mtDNA) lineages carried out in the African continent to date. We generated high-throughput mtDNA single nucleotide polymorphism (SNP) data (230 SNPs) from 2024 Africans, where more than 500 of them were additionally genotyped for the control region. These data were analyzed together with over 12,700 control region profiles collected from the literature, representing more than 300 population samples from Africa. Insights into the African homeland of humans are discussed. Phylogeographic patterns for the African continent are shown at a high phylogeographic resolution as well as at the population and regional levels. The deepest branch of the mtDNA tree, haplogroup L0, shows the highest sub-haplogroup diversity in Southeast and East Africa, suggesting this region as the homeland for modern humans. Several demographic estimates point to the coast as a facilitator of human migration in Africa, but the data indicate complex patterns, perhaps mirroring the effect of recent continental-scaled demographic events in re-shaping African mtDNA variability.

Mitochondrial DNA aberrations and pathophysiological implications in hematopoietic diseases, chronic inflammatory diseases, and cancers

Mitochondria are important intracellular organelles that produce energy for cellular development, differentiation, and growth. Mitochondrial DNA (mtDNA) presents a 10- to 20-fold higher susceptibility to genetic mutations owing to the lack of introns and histone proteins. The mtDNA repair system is relatively inefficient, rendering it vulnerable to reactive oxygen species (ROS) produced during ATP synthesis within the mitochondria, which can then target the mtDNA. Under conditions of chronic inflammation and excess stress, increased ROS production can overwhelm the antioxidant system, resulting in mtDNA damage. This paper reviews recent literature describing the pathophysiological implications of oxidative stress, mitochondrial dysfunction, and mitochondrial genome aberrations in aging hematopoietic stem cells, bone marrow failure syndromes, hematological malignancies, solid organ cancers, chronic inflammatory diseases, and other diseases caused by exposure to environmental hazards.

Mitochondrial DNA copy number and chronic lymphocytic leukemia/small lymphocytic lymphoma risk in two prospective studies

BACKGROUND

Mitochondrial DNA copy number (mtDNA CN) may be modified by mitochondria in response to oxidative stress. Previously, mtDNA CN was associated with non-Hodgkin lymphoma (NHL) risk, particularly chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL). We conducted a replication study in the Prostate, Lung, Colorectal, and Ovarian (PLCO) study and pooled with published ATBC (Alpha-Tocopherol, Beta-Carotene) data.

METHODS

In PLCO, 292 NHL cases (95 CLL/SLL cases) and 301 controls were pooled with 142 NHL cases (47 CLL/SLL cases) and 142 controls from ATBC. Subjects answered a questionnaire and provided blood. DNA was extracted from prediagnostic peripheral white blood, and mtDNA CN assayed by quantitative polymerase chain reaction. Unconditional logistic regression estimated mtDNA CN and NHL risk by odds ratios (OR) and 95% confidence intervals (95% CI).

RESULTS

Greater mtDNA CN was associated with increased risk of CLL/SLL among males in PLCO (3rd vs. 1st tertile: OR, 2.21; 95% CI, 1.03-4.72; Ptrend: 0.049) and pooled (T3 vs. T1: OR, 3.12; 95% CI, 1.72-5.68; Ptrend: 0.0002). Association was stronger among male smokers (Ptrend: <0.0001) and essentially identical for cases diagnosed <6, >6-8, and >8 years from blood draw (pooled: Pinteraction: 0.65). mtDNA CN and risk of other NHL subtypes and multiple myeloma showed no association.

CONCLUSIONS AND IMPACT

Mitochondrial DNA CN was associated with risk of CLL/SLL in males/male smokers. The risk was observed among cases diagnosed as long as 8 years after blood draw. These results suggest that higher mtDNA CN may reflect a process involved in CLL/SLL development.

Static and moving frontiers: the genetic landscape of Southern African Bantu-speaking populations

A consensus on Bantu-speaking populations being genetically similar has emerged in the last few years, but the demographic scenarios associated with their dispersal are still a matter of debate. The frontier model proposed by archeologists postulates different degrees of interaction among incoming agropastoralist and resident foraging groups in the presence of "static" and "moving" frontiers. By combining mitochondrial DNA and Y chromosome data collected from several southern African populations, we show that Bantu-speaking populations from regions characterized by a moving frontier developing after a long-term static frontier have larger hunter-gatherer contributions than groups from areas where a static frontier was not followed by further spatial expansion. Differences in the female and male components suggest that the process of assimilation of the long-term resident groups into agropastoralist societies was gender biased. Our results show that the diffusion of Bantu languages and culture in Southern Africa was a process more complex than previously described and suggest that the admixture dynamics between farmers and foragers played an important role in shaping the current patterns of genetic diversity.

The saga of the many studies wrongly associating mitochondrial DNA with breast cancer

Background

A large body of genetic research has focused on the potential role that mitochondrial DNA (mtDNA) variants might play on the predisposition to common and complex (multi-factorial) diseases. It has been argued however that many of these studies could be inconclusive due to artifacts related to genotyping errors or inadequate design.

Methods

Analyses of the data published in case–control breast cancer association studies have been performed using a phylogenetic-based approach. Variation observed in these studies has been interpreted in the light of data available on public resources, which now include over >27,000 complete mitochondrial sequences and the worldwide phylogeny determined by these mitogenomes. Complementary analyses were carried out using public datasets of partial mtDNA sequences, mainly corresponding to control-region segments.

Results

By way of example, we show here another kind of fallacy in these medical studies, namely, the phenomenon of SNP-SNP interaction wrongly applied to haploid data in a breast cancer study. We also reassessed the mutually conflicting studies suggesting some functional role of the non-synonymous polymorphism m.10398A > G (ND3 subunit of mitochondrial complex I) in breast cancer. In some studies, control groups were employed that showed an extremely odd haplogroup frequency spectrum compared to comparable information from much larger databases. Moreover, the use of inappropriate statistics signaled spurious “significance” in several instances.

Conclusions

Every case–control study should come under scrutiny in regard to the plausibility of the control-group data presented and appropriateness of the statistical methods employed; and this is best done before potential publication.

Indian signatures in the westernmost edge of the European Romani diaspora: new insight from mitogenomes

In agreement with historical documentation, several genetic studies have revealed ancestral links between the European Romani and India. The entire mitochondrial DNA (mtDNA) of 27 Spanish Romani was sequenced in order to shed further light on the origins of this population. The data were analyzed together with a large published dataset (mainly hypervariable region I [HVS-I] haplotypes) of Romani (N=1,353) and non-Romani worldwide populations (N>150,000). Analysis of mitogenomes allowed the characterization of various Romani-specific clades. M5a1b1a1 is the most distinctive European Romani haplogroup; it is present in all Romani groups at variable frequencies (with only sporadic findings in non-Romani) and represents 18% of their mtDNA pool. Its phylogeographic features indicate that M5a1b1a1 originated 1.5 thousand years ago (kya; 95% CI: 1.3-1.8) in a proto-Romani population living in Northwest India. U3 represents the most characteristic Romani haplogroup of European/Near Eastern origin (12.4%); it appears at dissimilar frequencies across the continent (Iberia: ≈ 31%; Eastern/Central Europe: ≈ 13%). All U3 mitogenomes of our Iberian Romani sample fall within a new sub-clade, U3b1c, which can be dated to 0.5 kya (95% CI: 0.3-0.7); therefore, signaling a lower bound for the founder event that followed admixture in Europe/Near East. Other minor European/Near Eastern haplogroups (e.g. H24, H88a) were also assimilated into the Romani by introgression with neighboring populations during their diaspora into Europe; yet some show a differentiation from the phylogenetically closest non-Romani counterpart. The phylogeny of Romani mitogenomes shows clear signatures of low effective population sizes and founder effects. Overall, these results are in good agreement with historical documentation, suggesting that cultural identity and relative isolation have allowed the Romani to preserve a distinctive mtDNA heritage, with some features linking them unequivocally to their ancestral Indian homeland.

The genetic legacy of the pre-colonial period in contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.

The maintenance of mitochondrial genetic stability is crucial during the oncogenic process

The main energetic resources of the cell are the mitochondria. As such, these organelles control a number of processes related to the life and death of the cell and also have a prominent function in the maintenance of tumor cells. In the last years, several authors have proposed an active role for mitochondria in tumorigenesis, more specifically concerning somatic mutations in mitochondrial DNA (mtDNA). Here, we wanted to evaluate this hypothesis based on the conclusions obtained in a model of gliomagenesis with elevated levels of ROS (reactive oxygen species), a toxic by-product of tumor metabolism. According to our findings, none of the mtDNA variants were found relevant to the tumoral process or suggest the involvement of mitochondria in tumorigenesis beyond the metabolic requirements of the tumoral cell. We conclude that there is not enough evidence to support the claim that mitochondrial instability holds any relevant role in the tumoral process.

Reconstructing ancient mitochondrial DNA links between Africa and Europe

Mitochondrial DNA (mtDNA) lineages of macro-haplogroup L (excluding the derived L3 branches M and N) represent the majority of the typical sub-Saharan mtDNA variability. In Europe, these mtDNAs account for <1% of the total but, when analyzed at the level of control region, they show no signals of having evolved within the European continent, an observation that is compatible with a recent arrival from the African continent. To further evaluate this issue, we analyzed 69 mitochondrial genomes belonging to various L sublineages from a wide range of European populations. Phylogeographic analyses showed that ~65% of the European L lineages most likely arrived in rather recent historical times, including the Romanization period, the Arab conquest of the Iberian Peninsula and Sicily, and during the period of the Atlantic slave trade. However, the remaining 35% of L mtDNAs form European-specific subclades, revealing that there was gene flow from sub-Saharan Africa toward Europe as early as 11,000 yr ago.

Genetic continuity in the Franco-Cantabrian region: new clues from autochthonous mitogenomes

Background

The Late Glacial Maximum (LGM), ∼20 thousand years ago (kya), is thought to have forced the people inhabiting vast areas of northern and central Europe to retreat to southern regions characterized by milder climatic conditions. Archaeological records indicate that Franco-Cantabria might have been the major source for the re-peopling of Europe at the beginning of the Holocene (11.5 kya). However, genetic evidence is still scarce and has been the focus of an intense debate.

Methods/Principal Findings

Based on a survey of more than 345,000 partial control region sequences and the analysis of 53 mitochondrial DNA (mtDNA) genomes, we identified an mtDNA lineage, HV4a1a, which most likely arose in the Franco-Cantabrian area about 5.4 kya and remained confined to northern Iberia.

Conclusions/Significance

The HV4a1a lineage and several of its younger branches reveal for the first time genetic continuity in this region and long-term episodes of isolation. This, in turn, could at least in part explain the unique linguistic and cultural features of the Basque region.

The mitochondrial genome is a "genetic sanctuary" during the oncogenic process

Since Otto Warburg linked mitochondrial physiology and oncogenesis in the 1930s, a number of studies have focused on the analysis of the genetic basis for the presence of aerobic glycolysis in cancer cells. However, little or no evidence exists today to indicate that mtDNA mutations are directly responsible for the initiation of tumor onset. Based on a model of gliomagenesis in the mouse, we aimed to explore whether or not mtDNA mutations are associated with the initiation of tumor formation, maintenance and aggressiveness. We reproduced the different molecular events that lead from tumor initiation to progression in the mouse glioma. In human gliomas, most of the genetic alterations that have been previously identified result in the aberrant activation of different signaling pathways and deregulation of the cell cycle. Our data indicates that mitochondrial dysfunction is associated with reactive oxygen species (ROS) generation, leading to increased nuclear DNA (nDNA) mutagenesis, but maintaining the integrity of the mitochondrial genome. In addition, mutational stability has been observed in entire mtDNA of human gliomas; this is in full agreement with the results obtained in the cancer mouse model. We use this model as a paradigm of oncogenic transformation due to the fact that mutations commonly found in gliomas appear to be the most common molecular alterations leading to tumor development in most types of human cancer. Our results indicate that the mtDNA genome is kept by the cell as a "genetic sanctuary" during tumor development in the mouse and humans. This is compatible with the hypothesis that the mtDNA molecule plays an essential role in the control of the cellular adaptive survival response to tumor-induced oxidative stress. The integrity of mtDNA seems to be a necessary element for responding to the increased ROS production associated with the oncogenic process.

A reduced number of mtSNPs saturates mitochondrial DNA haplotype diversity of worldwide population groups

Background

The high levels of variation characterising the mitochondrial DNA (mtDNA) molecule are due ultimately to its high average mutation rate; moreover, mtDNA variation is deeply structured in different populations and ethnic groups. There is growing interest in selecting a reduced number of mtDNA single nucleotide polymorphisms (mtSNPs) that account for the maximum level of discrimination power in a given population. Applications of the selected mtSNP panel range from anthropologic and medical studies to forensic genetic casework.

Methodology/Principal Findings

This study proposes a new simulation-based method that explores the ability of different mtSNP panels to yield the maximum levels of discrimination power. The method explores subsets of mtSNPs of different sizes randomly chosen from a preselected panel of mtSNPs based on frequency. More than 2,000 complete genomes representing three main continental human population groups (Africa, Europe, and Asia) and two admixed populations (“African-Americans” and “Hispanics”) were collected from GenBank and the literature, and were used as training sets. Haplotype diversity was measured for each combination of mtSNP and compared with existing mtSNP panels available in the literature. The data indicates that only a reduced number of mtSNPs ranging from six to 22 are needed to account for 95% of the maximum haplotype diversity of a given population sample. However, only a small proportion of the best mtSNPs are shared between populations, indicating that there is not a perfect set of “universal” mtSNPs suitable for all population contexts. The discrimination power provided by these mtSNPs is much higher than the power of the mtSNP panels proposed in the literature to date. Some mtSNP combinations also yield high diversity values in admixed populations.

Conclusions/Significance

The proposed computational approach for exploring combinations of mtSNPs that optimise the discrimination power of a given set of mtSNPs is more efficient than previous empirical approaches. In contrast to precedent findings, the results seem to indicate that only few mtSNPs are needed to reach high levels of discrimination power in a population, independently of its ancestral background.