Differences in mtDNA whole sequence between Tibetan and Han populations suggesting adaptive selection to high altitude

Mitochondrial Genomic Evidence of Selective Constraints in Small-Bodied Terrestrial Cetartiodactyla

Body size may drive the molecular evolution of mitochondrial genes in response to changes in energy requirements across species of different sizes. In this study, we perform selection pressure analysis and phylogenetic independent contrasts (PIC) to investigate the association between molecular evolution of mitochondrial genome protein-coding genes (mtDNA PCGs) and body size in terrestrial Cetartiodactyla. Employing selection pressure analysis, we observe that the average non-synonymous/synonymous substitution rate ratio (ω) of mtDNA PCGs is significantly reduced in small-bodied species relative to their medium and large counterparts. PIC analysis further confirms that ω values are positively correlated with body size (R2 = 0.162, p = 0.0016). Our results suggest that mtDNA PCGs of small-bodied species experience much stronger purifying selection as they need to maintain a heightened metabolic rate. On the other hand, larger-bodied species may face less stringent selective pressures on their mtDNA PCGs, potentially due to reduced relative energy expenditure per unit mass. Furthermore, we identify several genes that undergo positive selection, possibly linked to species adaptation to specific environments. Therefore, despite purifying selection being the predominant force in the evolution of mtDNA PCGs, positive selection can also occur during the process of adaptive evolution.

Whole mitochondrial genome analysis in highland Tibetans: further matrilineal genetic structure exploration

Introduction: The Qinghai-Tibet Plateau is one of the last terrestrial environments conquered by modern humans. Tibetans are among the few high-altitude settlers in the world, and understanding the genetic profile of Tibetans plays a pivotal role in studies of anthropology, genetics, and archaeology. Methods: In this study, we investigated the maternal genetic landscape of Tibetans based on the whole mitochondrial genome collected from 145 unrelated native Lhasa Tibetans. Molecular diversity indices, haplotype diversity (HD), Tajima's D and Fu's Fs were calculated and the Bayesian Skyline Plot was obtained to determining the genetic profile and population fluctuation of Lhasa Tibetans. To further explore the genetic structure of Lhasa Tibetans, we collected 107 East Asian reference populations to perform principal component analysis (PCA), multidimensional scaling (MDS), calculated Fst values and constructed phylogenetic tree. Results: The maternal genetic landscape of Tibetans showed obvious East Asian characteristics, M9a (28.28%), R (11.03%), F1 (12.41%), D4 (9.66%), N (6.21%), and M62 (4.14%) were the dominant haplogroups. The results of PCA, MDS, Fst and phylogenetic tree were consistent: Lhasa Tibetans clustered with other highland Tibeto-Burman speakers, there was obvious genetic homogeneity of Tibetans in Xizang, and genetic similarity between Tibetans and northern Han people and geographically adjacent populations was found. In addition, specific maternal lineages of Tibetans also be determined in this study. Discussion: In general, this study further shed light on long-time matrilineal continuity on the Tibetan Plateau and the genetic connection between Tibetans and millet famers in the Yellow River Basin, and further revealed that multiple waves of population interaction and admixture during different historical periods between lowland and highland populations shaped the maternal genetic profile of Tibetans.

The First Complete Mitochondrial Genome of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) and Phylogenetic Relationships within Infraorder Brachyura

Characterizing the complete mitochondrial genome (mitogenome) of an organism is useful for genomic studies in taxonomy and evolution. The mitogenomic characteristics of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) have never been studied. The present study decodes the first mitogenome of E. crenata by high-throughput sequencing (HTS). The length of the mitogenome is 15,597 bp, and it contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL), and 22 transfer RNA genes. There are 14 and 23 genes observed on the heavy and light strands, respectively. E. crenata possesses a trnH-cac translocation, with the trnH-cac shifted between trnE-gaa and trnF-ttc instead of the usual location between nad5 and nad4 in decapods. Phylogenetic analyses based on the current dataset of 33 Brachyuran mitogenomes indicate that E. crenata. is closely related to Ashtoret lunaris of Matutidae. The similar codon usage and rearrangements in the two species provide evidence for their close phylogenetic relationship. Positive selection analysis showed that one residue located in cox1 was identified as a positively selected site with high BEB value (>95%), indicating that this gene was under positive selection pressure. This study is the first complete mitogenome record for the family Goneplacidae, and the results obtained may improve the understanding of the phylogeny of Goneplacidae in Brachyura.

New Record of Hydrothermal Vent Squat Lobster (Munidopsis lauensis) Provides Evidence of a Dispersal Corridor between the Pacific and Indian Oceans

Hydrothermal vents are chemosynthetically driven ecosystems and one of the most extreme environments on Earth. Vent communities exhibit remarkable taxonomic novelty at the species and supra-species levels, and over 80% of vent species are endemic. Here, we used mitochondrial DNA to identify the biogeographic distribution of Munidopsis lauensis and the heme-binding regions of A1-type COX1 from six species (including M. lauensis) to investigate whether genetic variation in the protein structure affects oxygen-binding ability. We verified the identity of Indian Ocean specimens by comparing sequences from the barcoding gene mitochondrial cytochrome oxidase subunit 1 (COI) with known M. lauensis sequences from the NCBI database. The data show that these are the first recorded specimens of M. lauensis in the Indian Ocean; previously, this species had been reported only in the southwest Pacific. Our findings support the hypothesis that vent fauna in the Pacific and Indian Oceans can interact via active ridges. In the case of the mitochondrial DNA-binding site, the arrangement of heme-binding ligands and type A1 motif of M. lauensis was identical to that in other species. Moreover, our findings suggest that the mechanism of oxygen binding is well conserved among species from terrestrial organisms to hydrothermal extremophiles. Overall, dispersal of the same species to geologically separated hydrothermal vents and conserved heme-binding regions in mitochondrial proteins suggest that hydrothermal species might have evolved from shallow sea organisms and became distributed geographically using a dispersion corridor.

Insights into the Divergence of Chinese Ips Bark Beetles during Evolutionary Adaptation

Simple Summary

Bark beetle species of the genus Ips are among the major pests of Chinese conifer forests. Based on mitochondrial genome and SNP, we investigated the phylogenetic relationships and evolutionary trends of 19 populations of six Ips species that had serious outbreaks in recent years. Our results demonstrated the relationships between Ips evolution and host plants, pheromones, and altitudinal differences, and provided new insights into the mechanism of adaptive evolution of Ips bark beetles.

Abstract

Many bark beetles of the genus Ips are economically important insect pests that cause severe damage to conifer forests worldwide. In this study, sequencing the mitochondrial genome and restriction site-associated DNA of Ips bark beetles helps us understand their phylogenetic relationships, biogeographic history, and evolution of ecological traits (e.g., pheromones and host plants). Our results show that the same topology in phylogenetic trees constructed in different ways (ML/MP/BI) and with different data (mtDNA/SNP) helps us to clarify the phylogenetic relationships between Chinese Ips bark beetle populations and Euramerican species and their higher order clades; Ips bark beetles are polyphyletic. The structure of the mitochondrial genome of Ips bark beetles is similar and conserved to some extent, especially in the sibling species Ips typographus and Ips nitidus. Genetic differences among Ips species are mainly related to their geographic distribution and different hosts. The evolutionary pattern of aggregation pheromones of Ips species reflects their adaptations to the environment and differences among hosts in their evolutionary process. The evolution of Ips species is closely related to the uplift of the Qinghai-Tibet Plateau and host switching. Our study addresses the evolutionary trend and phylogenetic relationships of Ips bark beetles in China, and also provides a new perspective on the evolution of bark beetles and their relationships with host plants and pheromones.

Association of variants m.T16172C and m.T16519C in whole mtDNA sequences with high altitude pulmonary edema in Han Chinese lowlanders

Background

High altitude pulmonary edema (HAPE) is a hypoxia-induced non-cardiogenic pulmonary edema that typically occurred in un-acclimatized lowlanders, which inevitably leads to life-threatening consequences. Apart from multiple factors involved, the genetic factors also play an important role in the pathogenesis of HAPE. So far, researchers have put more energy into the nuclear genome and HAPE, and ignored the relationship between the mitochondrion DNA (mtDNA) variants and HAPE susceptibility.

Methods

We recruited a total of 366 individuals including 181 HAPE patients and 185 non-HAPE populations through two times. The first time, 49 HAPE patients and 58 non-HAPE individuals were performed through whole mtDNA sequences to search the mutations and haplogroups. The second time, 132 HAPE patients and 127 non-HAPE subjects were collected to apply verifying these mutations and haplogroups of mtDNA with the routine PCR method.

Results

We analyzed and summarized the clinical characteristics and sequence data for the 49 HAPE patients and 58 non-HAPE individuals. We found that a series of routine blood indexes including systolic arterial blood pressure (SBP), heart rate (HR), white blood cell (WBC), and C-reactive protein (CRP) in the HAPE group presented higher and displayed significant differences compared with those in the non-HAPE group. Although the average numbers of variants in different region and group samples were not statistically significant (P > 0.05), the mutation densities of different regions in the internal group showed significant differences. Then we found two mutations (T16172C and T16519C) associated with the HAPE susceptibility, the T16172C mutation increased the risk of HAPE, and the T16519C mutation decreased the HAPE rating. Furthermore, the two mutations were demonstrated with 132 HAPE patients and 127 non-HAPE individuals. Unfortunately, all the haplogroups were not associated with the HAPE haplogroups.

Conclusions

We provided evidence of differences in mtDNA polymorphism frequencies between HAPE and non-HAPE Han Chinese. Genotypes of mtDNA 16172C and 16519C were correlated with HAPE susceptibility, indicating the role of the mitochondrial genome in the pathogenesis of HAPE.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12890-021-01791-1.

High-altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high-altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high-altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein-coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high-altitude vertebrates than in low-altitude vertebrates. The seven rapidly evolving genes were shared by the high-altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high-altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high-altitude vertebrates was higher than in low-altitude vertebrates as their evolution requires more energy in a high-altitude environment. Our study demonstrates the high-altitude environment (low atmospheric O2 levels) drives vertebrate evolution in mtDNA PCGs.

Mitochondrial genetic variation in human bioenergetics, adaptation, and adult disease

OBJECTIVES

Mitochondria are critical for the survival of eukaryotic organisms due to their ability to produce cellular energy, which drives virtually all aspects of host biology. However, the effects of mitochondrial DNA (mtDNA) variation in relation to disease etiology and adaptation within contemporary global human populations remains incompletely understood.

METHODS

To develop a more holistic understanding of the role of mtDNA diversity in human adaptation, health, and disease, we investigated mitochondrial biology and bioenergetics. More specifically, we synthesized details from studies of mitochondrial function and variation in the context of haplogroup background, climatic adaptation, and oxidative disease.

RESULTS

The majority of studies show that mtDNA variation arose during modern human dispersal around the world. Some of these variants appear to have been positively selected for their adaptiveness in colder climates, with these sequence changes having implications for tissue-specific function and thermogenic capacity. In addition, many variants modulating energy production are also associated with damaging metabolic byproducts and mitochondrial dysfunction, which, in turn, are implicated in the onset and severity of several different adult mitochondrial diseases. Thus, mtDNA variation that governs bioenergetics, metabolism, and thermoregulation may potentially have adverse consequences for human health, depending on the genetic background and context in which it occurs.

CONCLUSIONS

Our review suggests that the mitochondrial research field would benefit from independently replicating mtDNA haplogroup-phenotype associations across global populations, incorporating potentially confounding environmental, demographic, and disease covariates into studies of mtDNA variation, and extending association-based studies to include analyses of complete mitogenomes and assays of mitochondrial function.

The Complete Mitogenome of Elymus sibiricus and Insights Into Its Evolutionary Pattern Based on Simple Repeat Sequences of Seed Plant Mitogenomes

The most intriguing characteristics of plant mitochondrial genomes (mitogenomes) include their high variation in both sequence and structure, the extensive horizontal gene transfer (HGT), and the important role they play in hypoxic adaptation. However, the investigation of the mechanisms of hypoxic adaptation and HGT in plant mitochondria remains challenging due to the limited number of sequenced mitogenomes and non-coding nature of the transferred DNA. In this study, the mitogenome of Elymus sibiricus (Gramineae, Triticeae), a perennial grass species native to the Qinghai-Tibet plateau (QTP), was de novo assembled and compared with the mitogenomes of eight Gramineae species. The unique haplotype composition and higher TE content compared to three other Triticeae species may be attributed to the long-term high-altitude plateau adaptability of E. sibiricus. We aimed to discover the connection between mitogenome simple sequence repeats (SSRs) (mt-SSRs) and HGT. Therefore, we predicted and annotated the mt-SSRs of E. sibiricus along with the sequencing of 87 seed plants. The clustering result based on all of the predicted compound mitogenome SSRs (mt-c-SSRs) revealed an expected synteny within systematic taxa and also inter-taxa. The mt-c-SSRs were annotated to 11 genes, among which "(ATA)3agtcaagtcaag (AAT)3" occurred in the nad5 gene of 8 species. The above-mentioned results further confirmed the HGT of mitogenomes sequences even among distant species from the aspect of mt-c-SSRs. Two genes, nad4 and nad7, possessed a vast number of SSRs in their intron regions across the seed plant mitogenomes. Furthermore, five pairs of SSRs developed from the mitogenome of E. sibiricus could be considered as potential markers to distinguish between the species E. sibiricus and its related sympatric species E. nutans.

Mitochondria Are Fundamental for the Emergence of Metazoans: On Metabolism, Genomic Regulation, and the Birth of Complex Organisms

Out of many intracellular bacteria, only the mitochondria and chloroplasts abandoned their independence billions of years ago and became endosymbionts within the host eukaryotic cell. Consequently, one cannot grow eukaryotic cells without their mitochondria, and the mitochondria cannot divide outside of the cell, thus reflecting interdependence. Here, we argue that such interdependence underlies the fundamental role of mitochondrial activities in the emergence of metazoans. Several lines of evidence support our hypothesis: (a) Differentiation and embryogenesis rely on mitochondrial function; (b) mitochondrial metabolites are primary precursors for epigenetic modifications (such as methyl and acetyl), which are critical for chromatin remodeling and gene expression, particularly during differentiation and embryogenesis; and (c) mitonuclear coregulation adapted to accommodate both housekeeping and tissue-dependent metabolic needs. We discuss the evolution of the unique mitochondrial genetic system, mitochondrial metabolites, mitonuclear coregulation, and their critical roles in the emergence of metazoans and in human disorders.

Analysis of very important pharmacogene variants in the Tibetan population from China

Personalized medicine, the treatment best suited for an individual, is a hot field of clinical research in the world. Many recent studies have shown that genetic variations have a great influence on the treatment. This study aimed to identify the distribution differences of very important pharmacogene (VIP) variants between the Tibetan population and the other 26 populations from the 1000 Genomes project. Based on the PharmGKB database, we successfully genotyped 50 VIP variants located in 27 genes in the Tibetan population. We also compared the genotype frequencies of VIP variants between Tibetan population and the other 26 populations. Without adjustment, the Chi-square test showed that the only significant variant between Tibetans and every other group was rs1801159 in dihydropyrimidine dehydrogenase (DPYD), followed by rs1800566 in NAD(P)H quinone dehydrogenase 1 (NQO1) and rs1051296 in solute carrier family 19 member 1 (SLC19A1). After Bonferroni's multiple adjustments, the genotype frequencies distribution of DPYD rs1801159 was found to be different in Tibetans compared to the other 26 groups, apart from ACB and ASW. Moreover, genetic structure/F-statistics (Fst) analysis and the phylogenetic tree illustrated that Tibetans had a closer affinity with CDX, CHB, CHS, JPT and KHV. Our data will complement pharmacogenomics information of the Tibetan population and provide theoretical support for the realization of individualized medical treatment for Tibetans in the future.

Mitogenomic Perspectives on the Adaptation to Extreme Alkaline Environment of Amur ide (Leuciscus waleckii)

Amur ide (Leuciscus waleckii, Family Cyprinidae) is widely distributed in Northeast Asia. L. waleckii usually inhabits freshwater environments but can also survive in the Lake Dali Nur, one of the most extreme aquatic environments on the earth, with an alkalinity up to 50 mmol/L (pH 9.6). To investigate mechanisms of mitogenomic evolution underlying adaptation to extreme environments, we determined 30 complete mitogenomes that included Lake Dali Nur (alkaline environment, AL) population and Amur basin (freshwater environment, FW) population. Through phylogenetic and divergence time analysis, we found that AL and FW populations forming distinct two groups which were consistent with geographic divergence (the formation of Lake Dali Nur). In addition, we found that almost of the windows exhibited higher nucleotide diversity in FW population (avg 0.0046) than AL population (avg 0.0012). This result indicated that severe environment selection had remarkably reduced the genetic diversity of mitogenome in AL population and suggested that severe environment selection had remarkably reduced the genetic diversity of mitogenome in the AL population. Compared with the FW population (ω = 0.064), the AL population (ω = 0.092) had a larger mean ω (dN/dS ratios) value for the 13 concatenated mitochondrial protein-coding genes, indicating that the high alkaline tolerated group had accumulated more nonsynonymous mutations. These nonsynonymous mutations had resulted in slightly beneficial amino acid changes that allowed adaption to the severe conditions. This study provides an additional view to decipher the adaptive mitogenome evolution of L. waleckii of the high alkaline environment.

A battle for transmission: the cooperative and selfish animal mitochondrial genomes

The mitochondrial genome is an evolutionarily persistent and cooperative component of metazoan cells that contributes to energy production and many other cellular processes. Despite sharing the same host as the nuclear genome, the multi-copy mitochondrial DNA (mtDNA) follows very different rules of replication and transmission, which translate into differences in the patterns of selection. On one hand, mtDNA is dependent on the host for its transmission, so selections would favour genomes that boost organismal fitness. On the other hand, genetic heterogeneity within an individual allows different mitochondrial genomes to compete for transmission. This intra-organismal competition could select for the best replicator, which does not necessarily give the fittest organisms, resulting in mito-nuclear conflict. In this review, we discuss the recent advances in our understanding of the mechanisms and opposing forces governing mtDNA transmission and selection in bilaterians, and what the implications of these are for mtDNA evolution and mitochondrial replacement therapy.

Mitochondria and the Origin of Species: Bridging Genetic and Ecological Perspectives on Speciation Processes

Mitochondria have been known to be involved in speciation through the generation of Dobzhansky-Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky-Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.

Mitogenome evolution in ladybirds: Potential association with dietary adaptation

Dietary shifts can alter the relative availability of different nutrients and are therefore associated with metabolic adaptation in animals. The Coccinellidae (ladybirds) exhibits three major types of feeding habits and provides a useful model to study the effects of dietary changes on the evolution of mitogenomes, which encode proteins directly involved in energy metabolism. Here, mitogenomes of three coccinellid species were newly sequenced. These data were combined with other ten previously sequenced coccinellid mitogenomes to explore the relationship between mitogenome evolution and diets. Our results indicate that mitogenomic data can be effectively used to resolve phylogenetic relationships of Coccinellidae. Strong codon usage bias in coccinellid mitogenomes was predominantly determined by nucleotide composition. The 13 mitochondrial protein-coding genes (PCGs) globally evolved under negative constraints, with some PCGs showing a stronger purifying selection. Six PCGs (nad3, nad4L, and nad5 from Complex I; cox1 and cox3 from Complex IV; and atp6 from Complex V) displayed signs of positive selection. Of these, adaptive changes in cox3 were potentially associated with metabolic differences resulting from dietary shifts in Coccinellidae. Our results provide insights into the adaptive evolution of coccinellid mitogenomes in response to both dietary shifts and other life history traits.

The complete mitochondrial genome of Calyptogena marissinica (Heterodonta: Veneroida: Vesicomyidae): Insight into the deep-sea adaptive evolution of vesicomyids

The deep-sea chemosynthetic environment is one of the most extreme environments on the Earth, with low oxygen, high hydrostatic pressure and high levels of toxic substances. Species of the family Vesicomyidae are among the dominant chemosymbiotic bivalves found in this harsh habitat. Mitochondria play a vital role in oxygen usage and energy metabolism; thus, they may be under selection during the adaptive evolution of deep-sea vesicomyids. In this study, the mitochondrial genome (mitogenome) of the vesicomyid bivalve Calyptogena marissinica was sequenced with Illumina sequencing. The mitogenome of C. marissinica is 17,374 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL) and 22 transfer RNA genes. All of these genes are encoded on the heavy strand. Some special elements, such as tandem repeat sequences, “G(A)_nT” motifs and AT-rich sequences, were observed in the control region of the C. marissinica mitogenome, which is involved in the regulation of replication and transcription of the mitogenome and may be helpful in adjusting the mitochondrial energy metabolism of organisms to adapt to the deep-sea chemosynthetic environment. The gene arrangement of protein-coding genes was identical to that of other sequenced vesicomyids. Phylogenetic analyses clustered C. marissinica with previously reported vesicomyid bivalves with high support values. Positive selection analysis revealed evidence of adaptive change in the mitogenome of Vesicomyidae. Ten potentially important adaptive residues were identified, which were located in cox1, cox3, cob, nad2, nad4 and nad5. Overall, this study sheds light on the mitogenomic adaptation of vesicomyid bivalves that inhabit the deep-sea chemosynthetic environment.

Whole Mitochondrial DNA Sequencing Analysis in 47 Han Populations in Southwest China

Background

Mitochondrial DNA (mtDNA) sequencing has been used in many areas, including forensic genetics. Due to the rapid development of sequencing technology, whole mtDNA sequencing is now possible and may be used in epidemiological and forensic studies. This study aimed to use whole mtDNA sequencing to investigate 47 Chongqing Han populations in southwest China and the diversity in the mtGenome reference data.

Material/Methods

The mtDNA of 47 Chongqing Han populations was generated using the Ion Torrent Personal Genome Machine (PGM) system. The extent of the effects of the mtDNA on the subpopulations was investigated and compared with six other populations from published studies. Pairwise fixation index (FST), a measure of population differentiation due to genetic structure, were calculated. Analysis of molecular variance (AMOVA) was performed, and 1257 hypervariable region data sets were added to the principal component analysis (PCA).

Results

The whole mtDNA sequencing data of 47 southwest Chinese Han populations were successfully recovered. Expanding the sequencing rage increased the discrimination power of mtDNA from three-times to 25-times based on different populations. The subpopulation effects showed 20 times the differences in match probability when compared with south China regions.

Conclusions

Whole mtDNA sequencing distinguished between individuals from 47 Chongqing Han populations in southwest China and has potential applications that include high-quality forensic identification.

Mitochondrial DNA sequencing reveals association of variants and haplogroup M33a2'3 with High altitude pulmonary edema susceptibility in Indian male lowlanders

High Altitude Pulmonary Edema (HAPE) is a threatening disorder caused due to acute exposure to high altitude above 3000 m. Apart from multiple factors involved, the genetic factors also play an important function in the pathogenesis of HAPE. This study aims to evaluate the role of mtDNA polymorphism and their association with haplogroup in understanding the etiology of HAPE. In this study, all the HAPE susceptible and acclimatized control subjects could be classified into nine haplogroups pertaining mostly to Macrohaplogroup M and U. The frequency of haplogroup M was significantly higher in HAPE susceptibles whereas the haplogroup M33a2'3 was found only in HAPE susceptibles. The variant G4491A and A4944G of MT-ND2, A14002G of MT-ND5, and C8562T of MT-ATP8, were definition site of haplogroup M33a2'3. The frequency of A10398G of MT-ND3, A8701G of MT-ATP6 and C14766T of MT-CYB genes were significantly higher in HAPE susceptibles. mtDNA copy number also plays a significant synergistic role in HAPE susceptibility. Our findings suggests that variants in MT-ND2 and MT-ND5 were predicted to confer decreased protein stability in HAPE susceptibles and in particular, highly conserved variants G4491A, A4944G and A14002G associated with haplogroup M33a2'3 may be the primary cause of susceptibility to HAPE in Indian male lowlanders.

Convergent Evolution of Mitochondrial Genes in Deep-Sea Fishes

Deep seas have extremely harsh conditions including high hydrostatic pressure, total darkness, cold, and little food and oxygen. The adaptations of fishes to deep-sea environment apparently have occurred independently many times. The genetic basis of adaptation for obtaining their energy remains unknown. Mitochondria play a central role in aerobic respiration. Analyses of the available 2,161 complete mitochondrial genomes of 1,042 fishes, including 115 deep-sea species, detect signals of positive selection in mitochondrial genes in nine branches of deep-sea fishes. Aerobic metabolism yields much more energy per unit of source material than anaerobic metabolism. The adaptive evolution of the mtDNA may reflect that aerobic metabolism plays a more important role than anaerobic metabolism in deep-sea fishes, whose energy sources (food) are extremely limited. This strategy maximizes the usage of energy sources. Eleven mitochondrial genes have convergent/parallel amino acid changes between branches of deep-sea fishes. Thus, these amino acid sites may be functionally important in the acquisition of energy, and reflect convergent evolution during their independent invasion of the harsh deep-sea ecological niche.

Positive Selection Drove the Adaptation of Mitochondrial Genes to the Demands of Flight and High-Altitude Environments in Grasshoppers

The molecular evolution of mitochondrial genes responds to changes in energy requirements and to high altitude adaptation in animals, but this has not been fully explored in invertebrates. The evolution of atmospheric oxygen content from high to low necessarily affects the energy requirements of insect movement. We examined 13 mitochondrial protein-coding genes (PCGs) of grasshoppers to test whether the adaptive evolution of genes involved in energy metabolism occurs in changes in atmospheric oxygen content and high altitude adaptation. Our molecular evolutionary analysis of the 13 PCGs in 15 species of flying grasshoppers and 13 related flightless grasshoppers indicated that, similar to previous studies, flightless grasshoppers have experienced relaxed selection. We found evidence of significant positive selection in the genes ATP8, COX3, ND2, ND4, ND4L, ND5, and ND6 in flying lineages. This results suggested that episodic positive selection allowed the mitochondrial genes of flying grasshoppers to adapt to increased energy demands during the continuous reduction of atmospheric oxygen content. Our analysis of five grasshopper endemic to the Tibetan Plateau and 13 non-Tibetan grasshoppers indicated that, due to positive selection, more non-synonymous nucleotide substitutions accumulated in Tibetan grasshoppers than in non-Tibetan grasshoppers. We also found evidence for significant positive selection in the genes ATP6, ND2, ND3, ND4, and ND5 in Tibetan lineages. Our results thus strongly suggest that, in grasshoppers, positive selection drives mitochondrial genes to better adapt both to the energy requirements of flight and to the high altitude of the Tibetan Plateau.

Genetic polymorphisms of the drug-metabolizing enzyme CYP2J2 in a Tibetan population

As an important metabolic enzyme, it is necessary to investigate the genetic polymorphisms of CYP2J2 among healthy Tibetan individuals. Genetic polymorphisms of CYP2J2 could affect enzyme activity and lead to differences among individual responses to drugs.We sequenced the whole gene of CYP2J2 in 100 unrelated, healthy Tibetan volunteers from the Tibet Autonomous Region and screened for genetic variants in the promoters, introns, exons, and the 3'-UTR regions.We detected 4 novel genetic polymorphisms of the CYP2J2 gene. The allelic frequencies of CYP2D6*1 and *7 were 0.955 and 0.045, respectively. CYP2D6*1/*7 decreased the activity of CYP2J2 and was expressed in 9% of the sample population.Our results provided basic data about CYP2J2 polymorphisms in a Tibetan population, suggested that the enzymatic activities of CYP2J2 might be different within the ethnic group, and offered a theoretical basis for individualized medical treatment and drug genomics studies.

Analysis of mitochondrial DNA sequence and copy number variation across five high-altitude species and their low-altitude relatives

High-altitude inhospitable environments impose a formidable life challenge for the local animals. Training and exposure to high-altitude environments produce both distinct physiological and phenotypic characteristics. The mitochondrion, an organelle crucial for the energy production, plays an important role in hypoxia adaptation. In this study, we investigated the mitochondrial DNA (mtDNA) polymorphism and copy number variation between the population pairs from distinct altitudes across the multi-species. Higher mitochondrial DNA control region's genetic diversity is conspicuous in high-altitude animals versus low-altitude relatives. We also found an accordant decrease of mtDNA copy number in most of the tissues from high-altitude animals. Compared to mammals, chickens have significantly distinct mitogenomic characteristics, and more significant changes in the skeletal muscle mtDNA copy number between high- and low-altitude individuals. Our study catches a snapshot of the biological similarities and differences in the mitochondrial high-altitude acclimation across the species.

Mitochondrial genomes uncover the maternal history of the Pamir populations

The Pamirs, among the world's highest mountains in Central Asia, are one of homelands with the most extreme high altitude for several ethnic groups. The settlement history of modern humans on the Pamirs remains still opaque. Herein, we have sequenced the mitochondrial DNA (mtDNA) genomes of 382 individuals belonging to eight populations from the Pamirs and the surrounding lowlands in Central Asia. We construct the Central Asian (including both highlanders and lowlanders) mtDNA haplogroup tree at the highest resolution. All the matrilineal components are assigned into the defined mtDNA haplogroups in East and West Eurasians. No basal lineages that directly emanate from the Eurasian founder macrohaplogroups M, N, and R are found. Our data support the origin of Central Asian being the result of East-West Eurasian admixture. The coalescence ages for more than 93% mtDNA lineages in Central Asians are dated after the last glacial maximum (LGM). The post-LGM and/or later dispersals/admixtures play dominant roles in shaping the maternal gene pool of Central Asians. More importantly, our analyses reveal the mtDNA heterogeneity in the Pamir highlanders, not only between the Turkic Kyrgyz and the Indo-European Tajik groups, but also among three highland Tajiks. No evidence supports positive selection or relaxation of selective constraints in the mtDNAs of highlanders as compared to that of lowlanders. Our results suggest a complex history for the peopling of Pamirs by multiple waves of migrations from various genetic resources during different time scales.

Measuring high-altitude adaptation

High altitudes (>8,000 ft or 2,500 m) provide an experiment of nature for measuring adaptation and the physiological processes involved. Studies conducted over the past ~25 years in Andeans, Tibetans, and, less often, Ethiopians show varied but distinct O₂ transport traits from those of acclimatized newcomers, providing indirect evidence for genetic adaptation to high altitude. Short-term (acclimatization, developmental) and long-term (genetic) responses to high altitude exhibit a temporal gradient such that, although all influence O₂ content, the latter also improve O₂ delivery and metabolism. Much has been learned concerning the underlying physiological processes, but additional studies are needed on the regulation of blood flow and O₂ utilization. Direct evidence of genetic adaptation comes from single-nucleotide polymorphism (SNP)-based genome scans and whole genome sequencing studies that have identified gene regions acted upon by natural selection. Efforts have begun to understand the connections between the two with Andean studies on the genetic factors raising uterine blood flow, fetal growth, and susceptibility to Chronic Mountain Sickness and Tibetan studies on genes serving to lower hemoglobin and pulmonary arterial pressure. Critical for future studies will be the selection of phenotypes with demonstrable effects on reproductive success, the calculation of actual fitness costs, and greater inclusion of women among the subjects being studied. The well-characterized nature of the O₂ transport system, the presence of multiple long-resident populations, and relevance for understanding hypoxic disorders in all persons underscore the importance of understanding how evolutionary adaptation to high altitude has occurred.NEW & NOTEWORTHY Variation in O₂ transport characteristics among Andean, Tibetan, and, when available, Ethiopian high-altitude residents supports the existence of genetic adaptations that improve the distribution of blood flow to vital organs and the efficiency of O₂ utilization. Genome scans and whole genome sequencing studies implicate a broad range of gene regions. Future studies are needed using phenotypes of clear relevance for reproductive success for determining the mechanisms by which naturally selected genes are acting.

Detecting signatures of selection within the Tibetan sheep mitochondrial genome

Tibetan sheep, a Chinese indigenous breed, are mainly distributed in plateau and mountain-valley areas at a terrestrial elevation between 2260 and 4100 m. The herd is genetically distinct from the other domestic sheep and undergoes acclimatization to adapt to the hypoxic environment. To date, whether the mitochondrial DNA modification of Tibetan sheep shares the same feature as the other domestic breed remains unknown. In this study, we compared the whole mitogenome sequences from 32 Tibetan sheep, 22 domestic sheep and 24 commercial sheep to identify the selection signatures of hypoxic-tolerant in Tibetan sheep. Nucleotide diversity analysis using the sliding window method showed that the highest level of nucleotide diversity was observed in the control region with a peak value of π = 0.05215, while the lowest π value was detected in the tRNAs region. qPCR results showed that the relative mtDNA copy number in Tibetan sheep was significantly lower than that in Suffolk sheep. None of the mutations in 12S rRNA were fixed in Tibetan sheep, which indicated that there has been less artificial selection in this herd than the other domestic and commercial breeds. Although one site (1277G) might undergo the purifying selection, it was not identified as the breed-specific allele in Tibetan sheep. We proposed that nature selection was the main drive during the domestication of Tibetan sheep and single mutation (or locus) could not reveal the signature of selection as for the high diversity in the mitogenome of Tibetan sheep.

Mitochondrial haplogroup M9a1a1c1b is associated with hypoxic adaptation in the Tibetans

While hypoxic environment at high altitude remains a major challenge for travelers from low-altitude areas, Tibetans have adapted to the high-altitude environment. Mitochondria are the energy conversion and supplement centers in eukaryotic cells. In recent years, studies have found that the diversity of the mitochondrial genome may have a role in the adaptation to hypoxia in Tibetans. In this study, mitochondrial haplogroup classification and variant genotyping were performed in Tibetan and Han Chinese populations living at different altitudes. The frequencies of mitochondrial haplogroups B and M7 in the high-altitude population were significantly lower compared with those in the low-altitude population (P=0.003 and 0.029, respectively), whereas the frequencies of haplogroups G and M9a1a1c1b in the high-altitude group were significantly higher compared with those in the low-altitude group (P=0.01 and 0.002, respectively). The frequencies of T3394C and G7697A, which are the definition sites of haplogroup M9a1a1c1b, were significantly higher in the high-altitude group compared with that in the low-altitude group (P=0.012 and 0.02, respectively). Our results suggest that mitochondrial haplogroups B and M7 are associated with inadaptability to hypoxic environments, whereas haplogroups G and M9a1a1c1b may be associated with hypoxic adaptation. In particular, the T3394C and G7697A variants on haplogroup M9a1a1c1b may be the primary cause of adaptation to hypoxia.

Mitogenomic perspectives on the origin of Tibetan loaches and their adaptation to high altitude

Tibetan loaches are the largest group of Tibetan fishes and are well adapted to the Tibetan Plateau. To investigate the origin of Tibetan loaches and their adaptations to the Tibetan Plateau, we determined 32 complete mitochondrial genomes that included 29 Tibetan loach species, two Barbatula species and Schistura longus. By combining these newly determined sequences with other previously published mitochondrial genomes, we assembled a large mitogenomic data set (11,433 bp) of 96 species in the superfamily Cobitoidea, to investigate the phylogenetic status of the genus Triplophysa. The resulting phylogeny strongly supported that the genus Triplophysa forms a monophyletic group within Nemacheilidae. Our molecular dating time suggests that the lineage leading to the Tibetan loaches and other loaches diverged approximately 23.5 Ma, which falls within the period of recent major uplifts of the Tibetan Plateau in the Early Miocene. Selection analyses revealed that the mitochondrial protein-coding genes of Tibetan loaches have larger ratios of nonsynonymous to synonymous substitutions than do those of non-Tibetan loaches, indicating that Tibetan loaches accumulated more nonsynonymous mutations than non-Tibetan loaches and exhibited rapid evolution. Two positively selected sites were identified in the ATP8 and ND1 genes.

Analysis of mitochondrial DNA in Tibetan gastric cancer patients at high altitude

The highest risk areas of gastric cancer are currently Japan, Korea and China; Qinghai, a high-altitude area, has one of the highest gastric cancer rates in China. The incidence of gastric cancer is higher in the Tibetan ethnic group compared to that in the Han ethnic group in Qinghai. This study was conducted to determine the clinical characteristics of mitochondrial DNA (mtDNA) mutations and copy numbers among Tibetans with gastric cancer residing at high altitudes and investigate the association between adaptations to hypoxic conditions and oncogenesis. A total of 23 Tibetan gastric cancer patients and 40 matched controls were recruited in this study. Leukocyte mtDNA genes and copy numbers were analyzed. The haplogroups were classified based on mitochondrial gene sequences. A total of 56.5% of the study participants had used alcohol at some point in their lives and 73.9% were positive for Helicobacter pylori (H. pylori). Eight mutations in 8 mitochondrial genes were identified in 43.4% of the Tibetan cancer patient group. There were no significant differences in leukocyte mtDNA copy number levels based on smoking status, alchohol consumption, obesity or H. pylori infection between the control and cancer groups. Statistical differences were also not found between gastric cancer patients with and those without mtDNA mutations. The majority of Tibetan patients with gastric cancer belonged to the mitochondrial haplogroup M9. In conclusion, Tibetans with gastric cancer residing at high altitudes exhibited a wide spectrum of mtDNA mutations. However, leukocyte mtDNA copy numbers in stage II gastric cancer were not statistically different compared to those in healthy Tibetans.

Mito-nuclear co-evolution: the positive and negative sides of functional ancient mutations

Most cell functions are carried out by interacting factors, thus underlying the functional importance of genetic interactions between genes, termed epistasis. Epistasis could be under strong selective pressures especially in conditions where the mutation rate of one of the interacting partners notably differs from the other. Accordingly, the order of magnitude higher mitochondrial DNA (mtDNA) mutation rate as compared to the nuclear DNA (nDNA) of all tested animals, should influence systems involving mitochondrial-nuclear (mito-nuclear) interactions. Such is the case of the energy producing oxidative phosphorylation (OXPHOS) and mitochondrial translational machineries which are comprised of factors encoded by both the mtDNA and the nDNA. Additionally, the mitochondrial RNA transcription and mtDNA replication systems are operated by nDNA-encoded proteins that bind mtDNA regulatory elements. As these systems are central to cell life there is strong selection toward mito-nuclear co-evolution to maintain their function. However, it is unclear whether (A) mito-nuclear co-evolution befalls only to retain mitochondrial functions during evolution or, also, (B) serves as an adaptive tool to adjust for the evolving energetic demands as species' complexity increases. As the first step to answer these questions we discuss evidence of both negative and adaptive (positive) selection acting on the mtDNA and nDNA-encoded genes and the effect of both types of selection on mito-nuclear interacting factors. Emphasis is given to the crucial role of recurrent ancient (nodal) mutations in such selective events. We apply this point-of-view to the three available types of mito-nuclear co-evolution: protein-protein (within the OXPHOS system), protein-RNA (mainly within the mitochondrial ribosome), and protein-DNA (at the mitochondrial replication and transcription machineries).

Detecting mitochondrial signatures of selection in wild Tibetan pigs and domesticated pigs

Selection in genomic regions is prevalent in mammals; however, the effects of selection on the mitogenome are not clearly understood. We determined the complete mitochondrial DNA (mtDNA) sequences from six wild Tibetan pigs from the Tibetan plateau and four domestic pig breeds from the lowland of neighboring southwest China. Nucleotide diversity analysis using the sliding window method showed that the nucleotide diversity of wild Tibetan pigs in most regions of the mitogenome was higher than that of domestic pigs. The 12 s ribosomal RNA showed relatively lower nucleotide diversity in Tibetan pigs, suggesting purifying selection of these genes during high-altitude adaptation. More non-synonymous nucleotide substitutions in the ATP6 were found in wild Tibetan pigs, indicating adaptive selection in Tibetan pigs. The results suggested distinct impacts of natural selection and artificial selection upon the mitogenome, especially the mitochondrial signatures of adaptive evolution in wild Tibetan pigs under natural selection.

Developmental functional adaptation to high altitude: review

Various approaches have been used to understand the origins of the functional traits that characterize the Andean high-altitude native. Based on the conceptual framework of developmental functional adaptation which postulates that environmental influences during the period of growth and development have long lasting effects that may be expressed during adulthood, we initiated a series of studies addressed at determining the pattern of physical growth and the contribution of growth and development to the attainment of full functional adaptation to high-altitude of low and high altitude natives living under rural and urban conditions. Current research indicate that: (a) the pattern of growth at high altitude due to limited nutritional resources, physical growth in body size is delayed but growth in lung volumes is accelerated because of hypoxic stress); (b) low-altitude male and female urban natives can attain a full functional adaptation to high altitude by exposure to high-altitude hypoxia during the period of growth and development; (c) both experimental studies on animals and comparative human studies indicate that exposure to high altitude during the period of growth and development results in the attainment of a large residual lung volume; (d) this developmentally acquired enlarged residual lung volume and its associated increase in alveolar area when combined with the increased tissue capillarization and moderate increase in red blood cells and hemoglobin concentration contributes to the successful functional adaptation of the Andean high-altitude native to hypoxia; and (e) any specific genetic traits that are related to the successful functional adaptation of Andean high-altitude natives have yet to be identified.

Functional recurrent mutations in the human mitochondrial phylogeny: dual roles in evolution and disease

Mutations frequently reoccur in the human mitochondrial DNA (mtDNA). However, it is unclear whether recurrent mtDNA nodal mutations (RNMs), that is, recurrent mutations in stems of unrelated phylogenetic nodes, are functional and hence selectively constrained. To answer this question, we performed comprehensive parsimony and maximum likelihood analyses of 9,868 publicly available whole human mtDNAs revealing 1,606 single nodal mutations (SNMs) and 679 RNMs. We then evaluated the potential functionality of synonymous, nonsynonymous and RNA SNMs and RNMs. For synonymous mutations, we have implemented the Codon Adaptation Index. For nonsynonymous mutations, we assessed evolutionary conservation, and employed previously described pathogenicity score assessment tools. For RNA genes’ mutations, we designed a bioinformatic tool which compiled evolutionary conservation and potential effect on RNA structure. While comparing the functionality scores of nonsynonymous and RNA SNMs and RNMs with those of disease-causing mtDNA mutations, we found significant difference (P < 0.001). However, 24 RNMs and 67 SNMs had comparable values with disease-causing mutations reflecting their potential function thus being the best candidates to participate in adaptive events of unrelated lineages. Strikingly, some functional RNMs occurred in unrelated mtDNA lineages that independently altered susceptibility to the same diseases, thus suggesting common functionality. To our knowledge, this is the most comprehensive analysis of selective signatures in the mtDNA not only within proteins but also within RNA genes. For the first time, we discover virtually all positively selected RNMs in our phylogeny while emphasizing their dual role in past evolutionary events and in disease today.

Human adaptation to the hypoxia of high altitude: the Tibetan paradigm from the pregenomic to the postgenomic era

The Tibetan Plateau is one of the highest regions on Earth. Tibetan highlanders are adapted to life and reproduction in a hypoxic environment and possess a suite of distinctive physiological traits. Recent studies have identified genomic loci that have undergone natural selection in Tibetans. Two of these loci, EGLN1 and EPAS1, encode major components of the hypoxia-inducible factor transcriptional system, which has a central role in oxygen sensing and coordinating an organism's response to hypoxia, as evidenced by studies in humans and mice. An association between genetic variants within these genes and hemoglobin concentration in Tibetans at high altitude was demonstrated in some of the studies (8, 80, 96). Nevertheless, the functional variants within these genes and the underlying mechanisms of action are still not known. Furthermore, there are a number of other possible phenotypic traits, besides hemoglobin concentration, upon which natural selection may have acted. Integration of studies at the genomic level with functional molecular studies and studies in systems physiology has the potential to provide further understanding of human evolution in response to high-altitude hypoxia. The Tibetan paradigm provides further insight on the role of the hypoxia-inducible factor system in humans in relation to oxygen homeostasis.

Mitochondrial DNA response to high altitude: a new perspective on high-altitude adaptation

Mitochondria are the energy metabolism centers of the cell. More than 95% of cellular energy is produced by mitochondrial oxidative phosphorylation. Hypoxia affects a wide range of energy generation and consumption processes in animals. The most important mechanisms limiting ATP consumption increase the efficiency of ATP production and accommodate the reduced production of ATP by the body. All of these mechanisms relate to changes in mitochondrial function. Mitochondrial function can be affected by variations in mitochondrial DNA, including polymorphisms, content changes, and deletions. These variations play an important role in acclimatization or adaptation to hypoxia. In this paper, the association between mitochondrial genome sequences and high-altitude adaptation is reviewed.