Adaptive evolution of cytochrome c oxidase subunit VIII in anthropoid primates

Astrocytes Drive Divergent Metabolic Gene Expression in Humans and Chimpanzees

The human brain utilizes ∼20% of all of the body's metabolic resources, while chimpanzee brains use <10%. Although previous work shows significant differences in metabolic gene expression between the brains of primates, we have yet to fully resolve the contribution of distinct brain cell types. To investigate cell type-specific interspecies differences in brain gene expression, we conducted RNA-seq on neural progenitor cells, neurons, and astrocytes generated from induced pluripotent stem cells from humans and chimpanzees. Interspecies differential expression analyses revealed that twice as many genes exhibit differential expression in astrocytes (12.2% of all genes expressed) than neurons (5.8%). Pathway enrichment analyses determined that astrocytes, rather than neurons, diverged in expression of glucose and lactate transmembrane transport, as well as pyruvate processing and oxidative phosphorylation. These findings suggest that astrocytes may have contributed significantly to the evolution of greater brain glucose metabolism with proximity to humans.

Disentangling Positive Selection from Relaxed Selection in Animal Mitochondrial Genomes

AbstractDisentangling different types of selection is a common goal in molecular evolution. Elevated dN/dS ratios (the ratio of nonsynonymous to synonymous substitution rates) in focal lineages are often interpreted as signs of positive selection. Paradoxically, relaxed purifying selection can also result in elevated dN/dS ratios, but tests to distinguish these two causes are seldomly implemented. Here, we reevaluated seven case studies describing elevated dN/dS ratios in animal mitochondrial DNA (mtDNA) and their accompanying hypotheses regarding selection. They included flightless lineages versus flighted lineages in birds, bats, and insects and physiological adaptations in snakes, two groups of electric fishes, and primates. We found that elevated dN/dS ratios were often not caused by the predicted mechanism, and we sometimes found strong support for the opposite mechanism. We discuss reasons why energetic hypotheses may be confounded by other selective forces acting on mtDNA and caution against overinterpreting singular molecular signals, including elevated dN/dS ratios.

The mtDNA mutation spectrum in the PolG mutator mouse reveals germline and somatic selection

BACKGROUND

Mitochondrial DNA (mtDNA) codes for products necessary for electron transport and mitochondrial gene translation. mtDNA mutations can lead to human disease and influence organismal fitness. The PolG mutator mouse lacks mtDNA proofreading function and rapidly accumulates mtDNA mutations, making it a model for examining the causes and consequences of mitochondrial mutations. Premature aging in PolG mice and their physiology have been examined in depth, but the location, frequency, and diversity of their mtDNA mutations remain understudied. Identifying the locations and spectra of mtDNA mutations in PolG mice can shed light on how selection shapes mtDNA, both within and across organisms.

RESULTS

Here, we characterized somatic and germline mtDNA mutations in brain and liver tissue of PolG mice to quantify mutation count (number of unique mutations) and frequency (mutation prevalence). Overall, mtDNA mutation count and frequency were the lowest in the D-loop, where an mtDNA origin of replication is located, but otherwise uniform across the mitochondrial genome. Somatic mtDNA mutations have a higher mutation count than germline mutations. However, germline mutations maintain a higher frequency and were also more likely to be silent. Cytosine to thymine mutations characteristic of replication errors were the plurality of basepair changes, and missense C to T mutations primarily resulted in increased protein hydrophobicity. Unlike wild type mice, PolG mice do not appear to show strand asymmetry in mtDNA mutations. Indel mutations had a lower count and frequency than point mutations and tended to be short, frameshift deletions.

CONCLUSIONS

Our results provide strong evidence that purifying selection plays a major role in the mtDNA of PolG mice. Missense mutations were less likely to be passed down in the germline, and they were less likely to spread to high frequencies. The D-loop appears to have resistance to mutations, either through selection or as a by-product of replication processes. Missense mutations that decrease hydrophobicity also tend to be selected against, reflecting the membrane-bound nature of mtDNA-encoded proteins. The abundance of mutations from polymerase errors compared with reactive oxygen species (ROS) damage supports previous studies suggesting ROS plays a minimal role in exacerbating the PolG phenotype, but our findings on strand asymmetry provide discussion for the role of polymerase errors in wild type organisms. Our results provide further insight on how selection shapes mtDNA mutations and on the aging mechanisms in PolG mice.

Molecular and Functional Effects of Loss of Cytochrome c Oxidase Subunit 8A

In this work we studied molecular and functional effects of the loss of the smallest nuclear encoded subunit of cytochrome c oxidase COX8A in fibroblasts from a patient with a homozygous splice site mutation and in CRISPR/Cas9 genome-edited HEK293T cells. In both cellular model systems, between 20 to 30% of the residual enzymatic activity of cytochrome c oxidase (COX) was detectable. In immunoblots of BN-PAGE separated mitochondria from both cellular models almost no monomers and dimers of the fully assembled COX could be visualized. Interestingly, supercomplexes of COX formed with complex III and also with complexes I and III retained considerable immunoreactivity, while nearly no immunoreactivity attributable to subassemblies was found. That indicates that COX lacking subunit 8A is stabilized in supercomplexes, while monomers and dimers are rapidly degraded. With transcriptome analysis by 3'-RNA sequencing we failed to detect in our cellular models of COX8A deficiency transcriptional changes of genes involved in the mitochondrial unfolded protein response (mtUPR) and the integrated stress response (ISR). Thus, our data strongly suggest that the smallest subunit of cytochrome c oxidase COX8A is required for maintenance of the structural stability of COX monomers and dimers.

Not all mitochondrial DNAs are made equal and the nucleus knows it

The oxidative phosphorylation (OXPHOS) system is the only structure in animal cells with components encoded by two genomes, maternally transmitted mitochondrial DNA (mtDNA), and biparentally transmitted nuclear DNA (nDNA). MtDNA‐encoded genes have to physically assemble with their counterparts encoded in the nucleus to build together the functional respiratory complexes. Therefore, structural and functional matching requirements between the protein subunits of these molecular complexes are rigorous. The crosstalk between nDNA and mtDNA needs to overcome some challenges, as the nuclear‐encoded factors have to be imported into the mitochondria in a correct quantity and match the high number of organelles and genomes per mitochondria that encode and synthesize their own components locally. The cell is able to sense the mito‐nuclear match through changes in the activity of the OXPHOS system, modulation of the mitochondrial biogenesis, or reactive oxygen species production. This implies that a complex signaling cascade should optimize OXPHOS performance to the cellular‐specific requirements, which will depend on cell type, environmental conditions, and life stage. Therefore, the mitochondria would function as a cellular metabolic information hub integrating critical information that would feedback the nucleus for it to respond accordingly. Here, we review the current understanding of the complex interaction between mtDNA and nDNA.

New Insights on the Evolution of the Sweet Taste Receptor of Primates Adapted to Harsh Environments

Simple Summary

A sense of taste is vital to an animal’s fitness. It enables animals to select and ingest beneficial foods and avoid harmful substances in their daily lives. There have been relatively few studies regarding the evolution of the taste receptor gene for sweetness, particularly in regard to endemic Bornean primates. We constructed TAS1R2 gene phylogenies for 20 anthropoid primate species using four different methods as well as established the phylogenic time divergence. The phylogenetic analysis successfully separated the primates into their taxonomic groups and as per their dietary preferences. Of note, the estimated time of divergence for the primate speciation pattern in this study was more recent than the previously published estimates. This difference may have been due to environmental changes, such as food scarcity and climate change, during the late Miocene epoch, which likely forced primates to adapt their dietary preferences. These findings establish a foundation for further investigations.

Abstract

Taste perception is an essential function that provides valuable dietary and sensory information, which is crucial for the survival of animals. Studies into the evolution of the sweet taste receptor gene (TAS1R2) are scarce, especially for Bornean endemic primates such as Nasalis larvatus (proboscis monkey), Pongo pygmaeus (Bornean orangutan), and Hylobates muelleri (Muller’s Bornean gibbon). Primates are the perfect taxa to study as they are diverse dietary feeders, comprising specialist folivores, frugivores, gummivores, herbivores, and omnivores. We constructed phylogenetic trees of the TAS1R2 gene for 20 species of anthropoid primates using four different methods (neighbor-joining, maximum parsimony, maximum-likelihood, and Bayesian) and also established the time divergence of the phylogeny. The phylogeny successfully separated the primates into their taxonomic groups as well as by their dietary preferences. Of note, the reviewed time of divergence estimation for the primate speciation pattern in this study was more recent than the previously published estimates. It is believed that this difference may be due to environmental changes, such as food scarcity and climate change, during the late Miocene epoch, which forced primates to change their dietary preferences. These findings provide a starting point for further investigation.

Role of cytochrome c oxidase nuclear-encoded subunits in health and disease

Cytochrome c oxidase (COX), the terminal enzyme of mitochondrial electron transport chain, couples electron transport to oxygen with generation of proton gradient indispensable for the production of vast majority of ATP molecules in mammalian cells. The review summarizes current knowledge of COX structure and function of nuclear-encoded COX subunits, which may modulate enzyme activity according to various conditions. Moreover, some nuclear-encoded subunits posess tissue-specific and development-specific isoforms, possibly enabling fine-tuning of COX function in individual tissues. The importance of nuclear-encoded subunits is emphasized by recently discovered pathogenic mutations in patients with severe mitopathies. In addition, proteins substoichiometrically associated with COX were found to contribute to COX activity regulation and stabilization of the respiratory supercomplexes. Based on the summarized data, a model of three levels of quaternary COX structure is postulated. Individual structural levels correspond to subunits of the i) catalytic center, ii) nuclear-encoded stoichiometric subunits and iii) associated proteins, which may constitute several forms of COX with varying composition and differentially regulated function.

The atypical subunit composition of respiratory complexes I and IV is associated with original extra structural domains in Euglena gracilis

In mitochondrial oxidative phosphorylation, electron transfer from NADH or succinate to oxygen by a series of large protein complexes in the inner mitochondrial membrane (complexes I-IV) is coupled to the generation of an electrochemical proton gradient, the energy of which is utilized by complex V to generate ATP. In Euglena gracilis, a non-parasitic secondary green alga related to trypanosomes, these respiratory complexes totalize more than 40 Euglenozoa-specific subunits along with about 50 classical subunits described in other eukaryotes. In the present study the Euglena proton-pumping complexes I, III, and IV were purified from isolated mitochondria by a two-steps liquid chromatography approach. Their atypical subunit composition was further resolved and confirmed using a three-steps PAGE analysis coupled to mass spectrometry identification of peptides. The purified complexes were also observed by electron microscopy followed by single-particle analysis. Even if the overall structures of the three oxidases are similar to the structure of canonical enzymes (e.g. from mammals), additional atypical domains were observed in complexes I and IV: an extra domain located at the tip of the peripheral arm of complex I and a "helmet-like" domain on the top of the cytochrome c binding region in complex IV.

Evidence of a Conserved Molecular Response to Selection for Increased Brain Size in Primates

The adaptive significance of human brain evolution has been frequently studied through comparisons with other primates. However, the evolution of increased brain size is not restricted to the human lineage but is a general characteristic of primate evolution. Whether or not these independent episodes of increased brain size share a common genetic basis is unclear. We sequenced and de novo assembled the transcriptome from the neocortical tissue of the most highly encephalized nonhuman primate, the tufted capuchin monkey (Cebus apella). Using this novel data set, we conducted a genome-wide analysis of orthologous brain-expressed protein coding genes to identify evidence of conserved gene–phenotype associations and species-specific adaptations during three independent episodes of brain size increase. We identify a greater number of genes associated with either total brain mass or relative brain size across these six species than show species-specific accelerated rates of evolution in individual large-brained lineages. We test the robustness of these associations in an expanded data set of 13 species, through permutation tests and by analyzing how genome-wide patterns of substitution co-vary with brain size. Many of the genes targeted by selection during brain expansion have glutamatergic functions or roles in cell cycle dynamics. We also identify accelerated evolution in a number of individual capuchin genes whose human orthologs are associated with human neuropsychiatric disorders. These findings demonstrate the value of phenotypically informed genome analyses, and suggest at least some aspects of human brain evolution have occurred through conserved gene–phenotype associations. Understanding these commonalities is essential for distinguishing human-specific selection events from general trends in brain evolution.

mtDNA polymorphism and metabolic inhibition affect sperm performance in conplastic mice

Whereas a broad link exists between nucleotide substitutions in the mitochondrial genome (mtDNA) and a range of metabolic pathologies, exploration of the effect of specific mtDNA genotypes is on-going. Mitochondrial DNA mutations are of particular relevance for reproductive traits, since they are expected to have profound effects on male specific processes as a result of the strict maternal inheritance of mtDNA. Sperm motility is crucially dependent on ATP in most systems studied. However, the importance of mitochondrial function in the production of the ATP necessary for sperm function remains uncertain. In this study, we test the effect of mtDNA polymorphisms upon mouse sperm performance and bioenergetics by using five conplastic inbred strains that share the same nuclear background while differing in their mitochondrial genomes. We found that, while genetic polymorphisms across distinct mtDNA haplotypes are associated with modification in sperm progressive velocity, this effect is not related to ATP production. Furthermore, there is no association between the number of mtDNA polymorphisms and either (a) the magnitude of sperm performance decrease, or (b) performance response to specific inhibition of the main sperm metabolic pathways. The observed variability between strains may be explained in terms of additive effects of single nucleotide substitutions on mtDNA coding sequences, which have been stabilized through genetic drift in the different laboratory strains. Alternatively, the decreased sperm performance might have arisen from the disruption of the nuclear DNA/mtDNA interactions that have coevolved during the radiation of Mus musculus subspecies.

Tissue- and Condition-Specific Isoforms of Mammalian Cytochrome c Oxidase Subunits: From Function to Human Disease

Cytochrome c oxidase (COX) is the terminal enzyme of the electron transport chain and catalyzes the transfer of electrons from cytochrome c to oxygen. COX consists of 14 subunits, three and eleven encoded, respectively, by the mitochondrial and nuclear DNA. Tissue- and condition-specific isoforms have only been reported for COX but not for the other oxidative phosphorylation complexes, suggesting a fundamental requirement to fine-tune and regulate the essentially irreversible reaction catalyzed by COX. This article briefly discusses the assembly of COX in mammals and then reviews the functions of the six nuclear-encoded COX subunits that are expressed as isoforms in specialized tissues including those of the liver, heart and skeletal muscle, lung, and testes: COX IV-1, COX IV-2, NDUFA4, NDUFA4L2, COX VIaL, COX VIaH, COX VIb-1, COX VIb-2, COX VIIaH, COX VIIaL, COX VIIaR, COX VIIIH/L, and COX VIII-3. We propose a model in which the isoforms mediate the interconnected regulation of COX by (1) adjusting basal enzyme activity to mitochondrial capacity of a given tissue; (2) allosteric regulation to adjust energy production to need; (3) altering proton pumping efficiency under certain conditions, contributing to thermogenesis; (4) providing a platform for tissue-specific signaling; (5) stabilizing the COX dimer; and (6) modulating supercomplex formation.

Mitogenomic analysis in European brown hare (Lepus europaeus) proposes genetic and functional differentiation between the distinct lineages

European brown hare is a small game species spreading across Europe to Asia Minor, with important economic traits. Population genetics studies using mitochondrial DNA markers have revealed the existence of two major phylogeographic lineages, the European and the Anatolian. European lineage is further divided in the European type halpogroup and south-eastern European type halpogroup, while Anatolian consists only by the Anatolian/Middle Eastern type halpogroup. All three haplogroups show a discrete geographical distribution, with an overlapping zone forming in North-East Greece and Bulgaria, forming a contact zone. Despite the existence of a contact zone, European haplotype was never detected in Anatolia and vice versa, proposing the presence of genetic barriers responsible for this phenomenon. In this study, we analyzed the whole mitochondrial genomes of specimens originating from both lineages, aiming to detect the genetic and functional differentiation of the oxidative phosphorylation complexes that are encoded by mtDNA that could lead gradually to the reproductive isolation of the lineages.

SNP Identification through Transcriptome Analysis of the European Brown Hare (Lepus europaeus): Cellular Energetics and Mother's Curse

The European brown hare (Lepus europaeus, Pallas 1778) is an important small game species in Europe. Due to its size and position in the food chain, as well as its life history, phenotypic variation and the relatively recent speciation events, brown hare plays an important role in the structure of various ecosystems and has emerged as an important species for population management and evolutionary studies. In order to identify informative SNPs for such studies, heart and liver tissues of three samples from the European lineage and a three-sample pool from the Anatolian lineage were subjected to RNA-Sequencing analysis. This effort resulted in 9496 well-assembled protein-coding sequences with close homology to human. After applying very stringent filtering criteria, 66185 polymorphic sites were identified in 7665 genes/cds and 2050 of those polymorphic sites are potentially capable of distinguishing the European from the Anatolian lineage. From these distinguishing mutations we focused on those in genes that are involved in cellular energy production, namely the glycolysis, Krebs cycle and the OXPHOS machinery. A selected set of SNPs was also validated by Sanger sequencing. By simulating the three European individuals as one pool, no substantial informative-SNP identification was lost, making it a cost-efficient approach. To our knowledge this is the first attempt to correlate the differentiation in both nuclear and mitochondrial genome between the two different lineages of L. europaeus with the observed spatial partitioning of the lineages of the species, proposing a possible mechanism that is maintaining the reproductive isolation of the lineages.

Loss of the smallest subunit of cytochrome c oxidase, COX8A, causes Leigh-like syndrome and epilepsy

Isolated cytochrome c oxidase (complex IV) deficiency is one of the most frequent respiratory chain defects in humans and is usually caused by mutations in proteins required for assembly of the complex. Mutations in nuclear-encoded structural subunits are very rare. In a patient with Leigh-like syndrome presenting with leukodystrophy and severe epilepsy, we identified a homozygous splice site mutation in COX8A, which codes for the ubiquitously expressed isoform of subunit VIII, the smallest nuclear-encoded subunit of complex IV. The mutation, affecting the last nucleotide of intron 1, leads to aberrant splicing, a frame-shift in the highly conserved exon 2, and decreased amount of the COX8A transcript. The loss of the wild-type COX8A protein severely impairs the stability of the entire cytochrome c oxidase enzyme complex and manifests in isolated complex IV deficiency in skeletal muscle and fibroblasts, similar to the frequent c.845_846delCT mutation in the assembly factor SURF1 gene. Stability and activity of complex IV could be rescued in the patient's fibroblasts by lentiviral expression of wild-type COX8A. Our findings demonstrate that COX8A is indispensable for function of human complex IV and its mutation causes human disease.

Conservative and compensatory evolution in oxidative phosphorylation complexes of angiosperms with highly divergent rates of mitochondrial genome evolution

Interactions between nuclear and mitochondrial gene products are critical for eukaryotic cell function. Nuclear genes encoding mitochondrial-targeted proteins (N-mt genes) experience elevated rates of evolution, which has often been interpreted as evidence of nuclear compensation in response to elevated mitochondrial mutation rates. However, N-mt genes may be under relaxed functional constraints, which could also explain observed increases in their evolutionary rate. To disentangle these hypotheses, we examined patterns of sequence and structural evolution in nuclear- and mitochondrial-encoded oxidative phosphorylation proteins from species in the angiosperm genus Silene with vastly different mitochondrial mutation rates. We found correlated increases in N-mt gene evolution in species with fast-evolving mitochondrial DNA. Structural modeling revealed an overrepresentation of N-mt substitutions at positions that directly contact mutated residues in mitochondrial-encoded proteins, despite overall patterns of conservative structural evolution. These findings support the hypothesis that selection for compensatory changes in response to mitochondrial mutations contributes to the elevated rate of evolution in N-mt genes. We discuss these results in light of theories implicating mitochondrial mutation rates and mitonuclear coevolution as drivers of speciation and suggest comparative and experimental approaches that could take advantage of heterogeneity in rates of mtDNA evolution across eukaryotes to evaluate such theories.

The Roles of Compensatory Evolution and Constraint in Aminoacyl tRNA Synthetase Evolution

Mitochondrial protein translation requires interactions between transfer RNAs encoded by the mitochondrial genome (mt-tRNAs) and mitochondrial aminoacyl tRNA synthetase proteins (mt-aaRS) encoded by the nuclear genome. It has been argued that animal mt-tRNAs have higher deleterious substitution rates relative to their nuclear-encoded counterparts, the cytoplasmic tRNAs (cyt-tRNAs). This dynamic predicts elevated rates of compensatory evolution of mt-aaRS that interact with mt-tRNAs, relative to aaRS that interact with cyt-tRNAs (cyt-aaRS). We find that mt-aaRS do evolve at significantly higher rates (exemplified by higher d_N and d_N/d_S) relative to cyt-aaRS, across mammals, birds, and Drosophila. While this pattern supports a model of compensatory evolution, the level at which a gene is expressed is a more general predictor of protein evolutionary rate. We find that gene expression level explains 10–56% of the variance in aaRS d_N/d_S, and that cyt-aaRS are more highly expressed in addition to having lower d_N/d_S values relative to mt-aaRS, consistent with more highly expressed genes being more evolutionarily constrained. Furthermore, we find no evidence of positive selection acting on either class of aaRS protein, as would be expected under a model of compensatory evolution. Nevertheless, the signature of faster mt-aaRS evolution persists in mammalian, but not bird or Drosophila, lineages after controlling for gene expression, suggesting some additional effect of compensatory evolution for mammalian mt-aaRS. We conclude that gene expression is the strongest factor governing differential amino acid substitution rates in proteins interacting with mitochondrial versus cytoplasmic factors, with important differences in mt-aaRS molecular evolution among taxonomic groups.

The subunit composition and function of mammalian cytochrome c oxidase

Cytochrome c oxidase (COX) from mammals and birds is composed of 13 subunits. The three catalytic subunits I-III are encoded by mitochondrial DNA, the ten nuclear-coded subunits (IV, Va, Vb, VIa, VIb, VIc, VIIa, VIIb, VIIc, VIII) by nuclear DNA. The nuclear-coded subunits are essentially involved in the regulation of oxygen consumption and proton translocation by COX, since their removal or modification changes the activity and their mutation causes mitochondrial diseases. Respiration, the basis for ATP synthesis in mitochondria, is differently regulated in organs and species by expression of tissue-, developmental-, and species-specific isoforms for COX subunits IV, VIa, VIb, VIIa, VIIb, and VIII, but the holoenzyme in mammals is always composed of 13 subunits. Various proteins and enzymes were shown, e.g., by co-immunoprecipitation, to bind to specific COX subunits and modify its activity, but these interactions are reversible, in contrast to the tightly bound 13 subunits. In addition, the formation of supercomplexes with other oxidative phosphorylation complexes has been shown to be largely variable. The regulatory complexity of COX is increased by protein phosphorylation. Up to now 18 phosphorylation sites have been identified under in vivo conditions in mammals. However, only for a few phosphorylation sites and four nuclear-coded subunits could a specific function be identified. Research on the signaling pathways leading to specific COX phosphorylations remains a great challenge for understanding the regulation of respiration and ATP synthesis in mammalian organisms. This article reviews the function of the individual COX subunits and their isoforms, as well as proteins and small molecules interacting and regulating the enzyme.

Mitochondrial-nuclear interactions: compensatory evolution or variable functional constraint among vertebrate oxidative phosphorylation genes?

Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate dS and nonsynonymous substitution rate dN) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d(N) of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d(N) between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

The first Chameleon transcriptome: comparative genomic analysis of the OXPHOS system reveals loss of COX8 in Iguanian lizards

Recently, we found dramatic mitochondrial DNA divergence of Israeli Chamaeleo chamaeleon populations into two geographically distinct groups. We aimed to examine whether the same pattern of divergence could be found in nuclear genes. However, no genomic resource is available for any chameleon species. Here we present the first chameleon transcriptome, obtained using deep sequencing (SOLiD). Our analysis identified 164,000 sequence contigs of which 19,000 yielded unique BlastX hits. To test the efficacy of our sequencing effort, we examined whether the chameleon and other available reptilian transcriptomes harbored complete sets of genes comprising known biochemical pathways, focusing on the nDNA-encoded oxidative phosphorylation (OXPHOS) genes as a model. As a reference for the screen, we used the human 86 (including isoforms) known structural nDNA-encoded OXPHOS subunits. Analysis of 34 publicly available vertebrate transcriptomes revealed orthologs for most human OXPHOS genes. However, OXPHOS subunit COX8 (Cytochrome C oxidase subunit 8), including all its known isoforms, was consistently absent in transcriptomes of iguanian lizards, implying loss of this subunit during the radiation of this suborder. The lack of COX8 in the suborder Iguania is intriguing, since it is important for cellular respiration and ATP production. Our sequencing effort added a new resource for comparative genomic studies, and shed new light on the evolutionary dynamics of the OXPHOS system.

Parallel evolution of IDH2 gene in cetaceans, primates and bats

Cetaceans and primates both have large brains that require large amounts of aerobic energy metabolism. In bats, the cost of flight makes locomotion energetically demanding. These mammalian groups may represent three independent evolutionary origins of an energy-demanding lifestyle in mammals. IDH2 encodes an enzyme in the tricarboxylic acid cycle in the mitochondrion, which plays a key role in aerobic energy metabolism. In this study, we cloned and sequenced this gene in two cetaceans, and 19 bat species, and compared the data with available primate sequences to test its evolution. We found significant signals of parallel evolution in this gene among these three groups. Parallel evolution of this gene may reflect their parallel evolution towards a higher demand for energy.

An Incompatibility between a mitochondrial tRNA and its nuclear-encoded tRNA synthetase compromises development and fitness in Drosophila

Mitochondrial transcription, translation, and respiration require interactions between genes encoded in two distinct genomes, generating the potential for mutations in nuclear and mitochondrial genomes to interact epistatically and cause incompatibilities that decrease fitness. Mitochondrial-nuclear epistasis for fitness has been documented within and between populations and species of diverse taxa, but rarely has the genetic or mechanistic basis of these mitochondrial–nuclear interactions been elucidated, limiting our understanding of which genes harbor variants causing mitochondrial–nuclear disruption and of the pathways and processes that are impacted by mitochondrial–nuclear coevolution. Here we identify an amino acid polymorphism in the Drosophila melanogaster nuclear-encoded mitochondrial tyrosyl–tRNA synthetase that interacts epistatically with a polymorphism in the D. simulans mitochondrial-encoded tRNA^Tyr to significantly delay development, compromise bristle formation, and decrease fecundity. The incompatible genotype specifically decreases the activities of oxidative phosphorylation complexes I, III, and IV that contain mitochondrial-encoded subunits. Combined with the identity of the interacting alleles, this pattern indicates that mitochondrial protein translation is affected by this interaction. Our findings suggest that interactions between mitochondrial tRNAs and their nuclear-encoded tRNA synthetases may be targets of compensatory molecular evolution. Human mitochondrial diseases are often genetically complex and variable in penetrance, and the mitochondrial–nuclear interaction we document provides a plausible mechanism to explain this complexity.

The ancient symbiosis between two prokaryotes that gave rise to the eukaryotic cell has required genomic cooperation for at least a billion years. Eukaryotic cells respire through the coordinated expression of their nuclear and mitochondrial genomes, both of which encode the proteins and RNAs required for mitochondrial transcription, translation, and aerobic respiration. Genetic interactions between these genomes are hypothesized to influence the effects of mitochondrial mutations on disease and drive mitochondrial–nuclear coevolution. Here we characterize the molecular cause and the cellular and organismal consequences of a mitochondrial–nuclear interaction in Drosophila between naturally occurring mutations in a mitochondrial tRNA and a nuclear-encoded tRNA synthetase. These mutations have little effect on their own; but, when combined, they severely compromise development and reproduction. tRNA synthetases attach the appropriate amino acid onto their cognate tRNA, and this reaction is required for efficient and accurate protein synthesis. We show that disruption of this interaction compromises mitochondrial function, providing hypotheses for the variable penetrance of diseases associated with mitochondrial tRNAs and for which pathways and processes are likely to be affected by mitochondrial–nuclear interactions.

Evolution of the couple cytochrome c and cytochrome c oxidase in primates

Mitochondrial energy metabolism has been affected by a broad set of ancient and recent evolutionary events. The oldest example is the endosymbiosis theory that led to mitochondria and a recently proposed example is adaptation to cold climate by anatomically modern human lineages. Mitochondrial energy metabolism has also been associated with an important area in anthropology and evolutionary biology, brain enlargement in human evolution. Indeed, several studies have pointed to the need for a major metabolic rearrangement to supply a sufficient amount of energy for brain development in primates.The genes encoding for the coupled cytochrome c (Cyt c) and cytochrome c oxidase (COX, complex IV, EC 1.9.3.1) seem to have an exceptional pattern of evolution in the anthropoid lineage. It has been proposed that this evolution was linked to the rearrangement of energy metabolism needed for brain enlargement. This hypothesis is reinforced by the fact that the COX enzyme was proposed to have a large role in control of the respiratory chain and thereby global energy production.After summarizing major events that occurred during the evolution of COX and cytochrome c on the primate lineage, we review the different evolutionary forces that could have influenced primate COX evolution and discuss the probable causes and consequences of this evolution. Finally, we discuss and review the co-occurring primate phenotypic evolution.

Is the COI barcoding gene involved in speciation through intergenomic conflict?

We here test the proposition that changes in the barcoding region of COI are commonly involved in speciation through intergenomic conflict. We demonstrate that this is unlikely given that even with incomplete taxon sampling, 78-90% of closely-related animal species have identical COI amino acid sequences. In addition, in those cases where amino acid substitutions between closely related species are observed, the inter- and intra-specific substitution patterns are very similar and/or lack consistent differences in the number, position and type of amino acid change. Overall, we conclude that there is little evidence for a widespread involvement of the barcoding gene in speciation.

Silencing, positive selection and parallel evolution: busy history of primate cytochromes C

Cytochrome c (cyt c) participates in two crucial cellular processes, energy production and apoptosis, and unsurprisingly is a highly conserved protein. However, previous studies have reported for the primate lineage (i) loss of the paralogous testis isoform, (ii) an acceleration and then a deceleration of the amino acid replacement rate of the cyt c somatic isoform, and (iii) atypical biochemical behavior of human cyt c. To gain insight into the cause of these major evolutionary events, we have retraced the history of cyt c loci among primates. For testis cyt c, all primate sequences examined carry the same nonsense mutation, which suggests that silencing occurred before the primates diversified. For somatic cyt c, maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses yielded the same tree topology. The evolutionary analyses show that a fast accumulation of non-synonymous mutations (suggesting positive selection) occurred specifically on the anthropoid lineage root and then continued in parallel on the early catarrhini and platyrrhini stems. Analysis of evolutionary changes using the 3D structure suggests they are focused on the respiratory chain rather than on apoptosis or other cyt c functions. In agreement with previous biochemical studies, our results suggest that silencing of the cyt c testis isoform could be linked with the decrease of primate reproduction rate. Finally, the evolution of cyt c in the two sister anthropoid groups leads us to propose that somatic cyt c evolution may be related both to COX evolution and to the convergent brain and body mass enlargement in these two anthropoid clades.

Relaxed clocks and inferences of heterogeneous patterns of nucleotide substitution and divergence time estimates across whales and dolphins (Mammalia: Cetacea)

Various nucleotide substitution models have been developed to accommodate among lineage rate heterogeneity, thereby relaxing the assumptions of the strict molecular clock. Recently developed "uncorrelated relaxed clock" and "random local clock" (RLC) models allow decoupling of nucleotide substitution rates between descendant lineages and are thus predicted to perform better in the presence of lineage-specific rate heterogeneity. However, it is uncertain how these models perform in the presence of punctuated shifts in substitution rate, especially between closely related clades. Using cetaceans (whales and dolphins) as a case study, we test the performance of these two substitution models in estimating both molecular rates and divergence times in the presence of substantial lineage-specific rate heterogeneity. Our RLC analyses of whole mitochondrial genome alignments find evidence for up to ten clade-specific nucleotide substitution rate shifts in cetaceans. We provide evidence that in the uncorrelated relaxed clock framework, a punctuated shift in the rate of molecular evolution within a subclade results in posterior rate estimates that are either misled or intermediate between the disparate rate classes present in baleen and toothed whales. Using simulations, we demonstrate abrupt changes in rate isolated to one or a few lineages in the phylogeny can mislead rate and age estimation, even when the node of interest is calibrated. We further demonstrate how increasing prior age uncertainty can bias rate and age estimates, even while the 95% highest posterior density around age estimates decreases; in other words, increased precision for an inaccurate estimate. We interpret the use of external calibrations in divergence time studies in light of these results, suggesting that rate shifts at deep time scales may mislead inferences of absolute molecular rates and ages.

Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.

Stability of mitochondrial membrane proteins in terrestrial vertebrates predicts aerobic capacity and longevity

The cellular energy produced by mitochondria is a fundamental currency of life. However, the extent to which mitochondrial (mt) performance (power and endurance) is adapted to habitats and life strategies of vertebrates is not well understood. A global analysis of mt genomes revealed that hydrophobicity (HYD) of mt membrane proteins (MMPs) is much lower in terrestrial vertebrates than in fishes and shows a strong negative correlation with serine/threonine composition (STC). Here, we present evidence that this systematic feature of MMPs was crucial for the evolution of large terrestrial vertebrates with high aerobic capacity. An Arrhenius-type equation gave positive correlations between STC and maximum life span (MLS) in terrestrial vertebrates (with a few exceptions relating to the lifestyle of small animals with a high resting metabolic rate [RMR]) and negative correlations in secondary marine vertebrates, such as cetaceans and alligators (which returned from land to water, utilizing buoyancy with increased body size). In particular, marked STC increases in primates (especially hominoids) among placentals were associated with very high MLS values. We connected these STC increases in MMPs with greater stability of respiratory complexes by estimating the degradation of the Arrhenius plot given by accelerating mtRMR up to mt maximum metabolic rate. Both mtRMR and HYD in terrestrial vertebrates decreased with increasing body mass. Decreases in mtRMR raise MMP stability when high mobility is not required, whereas decreased HYD may weaken this stability under the hydrophobic environment of lipid bilayer. High maximal metabolic rates (5–10 RMR), which we postulate require high MMP mobility, presumably render MMPs more unstable. A marked rise in STC may therefore be essential to stabilize MMPs, perhaps as dynamic supercomplexes, via hydrogen bonds associated with serine/threonine motifs.

Cytochrome c oxidase: evolution of control via nuclear subunit addition

According to theory, present eukaryotic cells originated from a beneficial association between two free-living cells. Due to this endosymbiotic event the pre-eukaryotic cell gained access to oxidative phosphorylation (OXPHOS), which produces more than 15 times as much ATP as glycolysis. Because cellular ATP needs fluctuate and OXPHOS both requires and produces entities that can be toxic for eukaryotic cells such as ROS or NADH, we propose that the success of endosymbiosis has largely depended on the regulation of endosymbiont OXPHOS. Several studies have presented cytochrome c oxidase as a key regulator of OXPHOS; for example, COX is the only complex of mammalian OXPHOS with known tissue-specific isoforms of nuclear encoded subunits. We here discuss current knowledge about the origin of nuclear encoded subunits and the appearance of different isozymes promoted by tissue and cellular environments such as hypoxia. We also review evidence for recent selective pressure acting on COX among vertebrates, particularly in primate lineages, and discuss the unique pattern of co-evolution between the nuclear and mitochondrial genomes. Finally, even though the addition of nuclear encoded subunits was a major event in eukaryotic COX evolution, this does not lead to emergence of a more efficient COX, as might be expected from an anthropocentric point of view, for the "higher" organism possessing large brains and muscles. The main function of these subunits appears to be "only" to control the activity of the mitochondrial subunits. We propose that this control function is an as yet under appreciated key point of evolution. Moreover, the importance of regulating energy supply may have caused the addition of subunits encoded by the nucleus in a process comparable to a "domestication scenario" such that the host tends to control more and more tightly the ancestral activity of COX performed by the mtDNA encoded subunits.

Adaptive evolution of energy metabolism genes and the origin of flight in bats

Bat flight poses intriguing questions about how flight independently developed in mammals. Flight is among the most energy-consuming activities. Thus, we deduced that changes in energy metabolism must be a primary factor in the origin of flight in bats. The respiratory chain of the mitochondrial produces 95% of the adenosine triphosphate (ATP) needed for locomotion. Because the respiratory chain has a dual genetic foundation, with genes encoded by both the mitochondrial and nuclear genomes, we examined both genomes to gain insights into the evolution of flight within mammals. Evidence for positive selection was detected in 23.08% of the mitochondrial-encoded and 4.90% of nuclear-encoded oxidative phosphorylation (OXPHOS) genes, but in only 2.25% of the nuclear-encoded nonrespiratory genes that function in mitochondria or 1.005% of other nuclear genes in bats. To address the caveat that the two available bat genomes are of only draft quality, we resequenced 77 OXPHOS genes from four species of bats. The analysis of the resequenced gene data are in agreement with our conclusion that a significantly higher proportion of genes involved in energy metabolism, compared with background genes, show evidence of adaptive evolution specific on the common ancestral bat lineage. Both mitochondrial and nuclear-encoded OXPHOS genes display evidence of adaptive evolution along the common ancestral branch of bats, supporting our hypothesis that genes involved in energy metabolism were targets of natural selection and allowed adaptation to the huge change in energy demand that were required during the origin of flight.

Both noncoding and protein-coding RNAs contribute to gene expression evolution in the primate brain

Despite striking differences in cognition and behavior between humans and our closest primate relatives, several studies have found little evidence for adaptive change in protein-coding regions of genes expressed primarily in the brain. Instead, changes in gene expression may underlie many cognitive and behavioral differences. Here, we used digital gene expression: tag profiling (here called Tag-Seq, also called DGE:tag profiling) to assess changes in global transcript abundance in the frontal cortex of the brains of 3 humans, 3 chimpanzees, and 3 rhesus macaques. A substantial fraction of transcripts we identified as differentially transcribed among species were not assayed in previous studies based on microarrays. Differentially expressed tags within coding regions are enriched for gene functions involved in synaptic transmission, transport, oxidative phosphorylation, and lipid metabolism. Importantly, because Tag-Seq technology provides strand-specific information about all polyadenlyated transcripts, we were able to assay expression in noncoding intragenic regions, including both sense and antisense noncoding transcripts (relative to nearby genes). We find that many noncoding transcripts are conserved in both location and expression level between species, suggesting a possible functional role. Lastly, we examined the overlap between differential gene expression and signatures of positive selection within putative promoter regions, a sign that these differences represent adaptations during human evolution. Comparative approaches may provide important insights into genes responsible for differences in cognitive functions between humans and nonhuman primates, as well as highlighting new candidate genes for studies investigating neurological disorders.

Adaptive threonine increase in transmembrane regions of mitochondrial proteins in higher primates

BACKGROUND

The mitochondrial (mt) gene tree of placental mammals reveals a very strong acceleration of the amino acid (AA) replacement rate and a change in AA compositional bias in the lineage leading to the higher primates (simians), in contrast to the nuclear gene tree. Whether this acceleration and compositional bias were caused by adaptive evolution at the AA level or directional mutation pressure at the DNA level has been vigorously debated.

METHODOLOGY/PRINCIPAL FINDINGS

Our phylogenetic analysis indicates that the rate acceleration in the simian lineage is accompanied by a marked increase in threonine (Thr) residues in the transmembrane helix regions of mt DNA-encoded proteins. This Thr increase involved the replacement of hydrophobic AAs in the membrane interior. Even after accounting for lack of independence due to phylogeny, a regression analysis reveals a statistical significant positive correlation between Thr composition and longevity in primates.

CONCLUSION/SIGNIFICANCE

Because crucial roles of Thr and Ser in membrane proteins have been proposed to be the formation of hydrogen bonds enhancing helix-helix interactions, the Thr increase detected in the higher primates might be adaptive by serving to reinforce stability of mt proteins in the inner membrane. The correlation between Thr composition in the membrane interior and the longevity of animals is striking, especially because some mt functions are thought to be involved in aging.

Evidence implicating the candidate schizophrenia/bipolar disorder susceptibility gene G72 in mitochondrial function

G72 is a strong candidate susceptibility gene for schizophrenia and bipolar disorder, whose function remains enigmatic. Here we show that one splicing isoform of the gene (LG72) encodes for a mitochondrial protein. We also provide convergent lines of evidence that increase of endogenous or exogenous G72 levels promotes robust mitochondrial fragmentation in mammalian cell lines and primary neurons, which proceeds in a manner that does not depend on induction of apoptosis or alteration in mitochondrial transmembrane potential. Finally, we show that increase in G72 levels in immature primary neurons is accompanied by a marked increase in dendritic arborization. By contrast, we failed to confirm the originally proposed functional interaction between G72 and D-amino acid oxidase (DAO) in two tested cell lines. Our results suggest an alternative role for G72 in modulating mitochondrial function.

Comparative genomics search for losses of long-established genes on the human lineage

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.

Molecular evolution of the cytochrome c oxidase subunit 5A gene in primates

BACKGROUND

Many electron transport chain (ETC) genes show accelerated rates of nonsynonymous nucleotide substitutions in anthropoid primate lineages, yet in non-anthropoid lineages the ETC proteins are typically highly conserved. Here, we test the hypothesis that COX5A, the ETC gene that encodes cytochrome c oxidase subunit 5A, shows a pattern of anthropoid-specific adaptive evolution, and investigate the distribution of this protein in catarrhine brains.

RESULTS

In a dataset comprising 29 vertebrate taxa, including representatives from all major groups of primates, there is nearly 100% conservation of the COX5A amino acid sequence among extant, non-anthropoid placental mammals. The most recent common ancestor of these species lived about 100 million years (MY) ago. In contrast, anthropoid primates show markedly elevated rates of nonsynonymous evolution. In particular, branch site tests identify five positively selected codons in anthropoids, and ancestral reconstructions infer that substitutions in these codons occurred predominantly on stem lineages (anthropoid, ape and New World monkey) and on the human terminal branch. Examination of catarrhine brain samples by immunohistochemistry characterizes for the first time COX5A protein distribution in the primate neocortex, and suggests that the protein is most abundant in the mitochondria of large-size projection neurons. Real time quantitative PCR supports previous microarray results showing COX5A is expressed in cerebral cortical tissue at a higher level in human than in chimpanzee or gorilla.

CONCLUSION

Taken together, these results suggest that both protein structural and gene regulatory changes contributed to COX5A evolution during humankind's ancestry. Furthermore, these findings are consistent with the hypothesis that adaptations in ETC genes contributed to the emergence of the energetically expensive anthropoid neocortex.

Encephalization, emergent properties, and psychiatry: a minicolumnar perspective

The focus of the authors' attention is the consequence of brain growth understood in terms of the development of networks of cortical cell minicolumns, the elemental information-processing units of the brain. The authors view cortical growth, encephalization, and the emergence of higher cognitive functions in humans as the consequence of an increase in the number of minicolumns and their connections. Encephalization has proceeded via weak linkages of canonical circuits, which facilitate the emergence of novel cortical functions. In addition to reframing the evolution of mind, this perspective provides a conceptual framework for a better understanding of the origin and maladaptive nature of certain psychiatric conditions.

Rate of evolution in brain-expressed genes in humans and other primates

Brain-expressed genes are known to evolve slowly in mammals. Nevertheless, since brains of higher primates have evolved rapidly, one might expect acceleration in DNA sequence evolution in their brain-expressed genes. In this study, we carried out full-length cDNA sequencing on the brain transcriptome of an Old World monkey (OWM) and then conducted three-way comparisons among (i) mouse, OWM, and human, and (ii) OWM, chimpanzee, and human. Although brain-expressed genes indeed appear to evolve more rapidly in species with more advanced brains (apes > OWM > mouse), a similar lineage effect is observable for most other genes. The broad inclusion of genes in the reference set to represent the genomic average is therefore critical to this type of analysis. Calibrated against the genomic average, the rate of evolution among brain-expressed genes is probably lower (or at most equal) in humans than in chimpanzee and OWM. Interestingly, the trend of slow evolution in coding sequence is no less pronounced among brain-specific genes, vis-à-vis brain-expressed genes in general. The human brain may thus differ from those of our close relatives in two opposite directions: (i) faster evolution in gene expression, and (ii) a likely slowdown in the evolution of protein sequences. Possible explanations and hypotheses are discussed.

When calibrated against the genomic average, the rate of evolution among brain-expressed genes in humans is probably lower than or equal to that of other closely related primates.

Whether comparing morphology or cognitive ability, it is clear that the human brain has evolved rapidly relative to that of other primates. But the extent to which genes expressed in the brain also reflect this overall pattern is unclear. To address this question, it's necessary to measure any variations in the DNA sequences of these genes between human and chimpanzee. And, to do this as accurately as possible, it's also important to require an appropriate reference group to act as a benchmark against which the differences can be measured. We therefore compared publicly available genomic sequences of chimps and humans with complementary DNA sequences of several thousand genes expressed in the brain of another closely related primate—the macaque, an Old World monkey—as well as the more distantly related mouse. Our analyses of the rates of protein evolution in these species suggest that genes expressed in the human brain have in fact slowed down in their evolution since the split between human and chimpanzee, contrary to some previously published reports. We suggest that advanced brains are driven primarily by the increasing complexity in the network of gene interactions. As a result, brain-expressed genes are constrained in their sequence evolution, although their expression levels may change rapidly.

Why do we still have a maternally inherited mitochondrial DNA? Insights from evolutionary medicine

The human cell is a symbiosis of two life forms, the nucleus-cytosol and the mitochondrion. The nucleus-cytosol emphasizes structure and its genes are Mendelian, whereas the mitochondrion specializes in energy and its mitochondrial DNA (mtDNA) genes are maternal. Mitochondria oxidize calories via oxidative phosphorylation (OXPHOS) to generate a mitochondrial inner membrane proton gradient (DeltaP). DeltaP then acts as a source of potential energy to produce ATP, generate heat, regulate reactive oxygen species (ROS), and control apoptosis, etc. Interspecific comparisons of mtDNAs have revealed that the mtDNA retains a core set of electron and proton carrier genes for the proton-translocating OXPHOS complexes I, III, IV, and V. Human mtDNA analysis has revealed these genes frequently contain region-specific adaptive polymorphisms. Therefore, the mtDNA with its energy controlling genes may have been retained to permit rapid adaptation to new environments.

The jewels of our genome: the search for the genomic changes underlying the evolutionarily unique capacities of the human brain

The recent publication of the initial sequence and analysis of the chimp genome allows us, for the first time, to compare our genome with that of our closest living evolutionary relative. With more primate genome sequences being pursued, and with other genome-wide, cross-species comparative techniques emerging, we are entering an era in which we will be able to carry out genomic comparisons of unprecedented scope and detail. These studies should yield a bounty of new insights about the genes and genomic features that are unique to our species as well as those that are unique to other primate lineages, and may begin to causally link some of these to lineage-specific phenotypic characteristics. The most intriguing potential of these new approaches will be in the area of evolutionary neurogenomics and in the possibility that the key human lineage–specific (HLS) genomic changes that underlie the evolution of the human brain will be identified. Such new knowledge should provide fresh insights into neuronal development and higher cognitive function and dysfunction, and may possibly uncover biological mechanisms for information storage, analysis, and retrieval never previously seen.

Adaptive selection of mitochondrial complex I subunits during primate radiation

Mammalian oxidative phosphorylation (OXPHOS) complexes I, III, IV and V are assembled from both mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) encoded subunits, with complex I encompassing 39 nDNA and seven mtDNA subunits. Yet the sequence variation of the mtDNA genes is more than ten fold greater than that of the nDNA encoded genes of the OXPHOS complexes and the mtDNA proteins have been found to be influenced by positive (adaptive) selection. To maintain a functional complex I, nDNA and mtDNA subunits must interact, implying that certain nDNA complex I genes may also have been influenced by positive selection. To determine if positive selection has influenced nDNA complex I genes, we analyzed the DNA sequences of all of the nDNA and mtDNA encoded complex I subunits from orangutan, gorilla, chimpanzee, human and all available vertebrate sequences. This revealed that three nDNA complex I genes (NDUFC2, NDUFA1, and NDUFA4) had significantly increased amino acid substitution rates by both PAML and Z-test, suggesting that they have been subjected to adaptive selection during primate radiation. Since all three of these subunits reside in the membrane domain of complex I along with the mtDNA subunits, we compared amino acid changes in these three nDNA genes with those of the mtDNA genes across species. Changes in the nDNA NDUFC2 cysteine 39 were found to correlate with those in the mtDNA ND5 cysteine 330. Therefore, adaptive selection has influenced some nDNA complex I genes and nDNA and mtDNA complex I genes may have co-evolved.

Snakes as agents of evolutionary change in primate brains

Current hypotheses that use visually guided reaching and grasping to explain orbital convergence, visual specialization, and brain expansion in primates are open to question now that neurological evidence reveals no correlation between orbital convergence and the visual pathway in the brain that is associated with reaching and grasping. An alternative hypothesis proposed here posits that snakes were ultimately responsible for these defining primate characteristics. Snakes have a long, shared evolutionary existence with crown-group placental mammals and were likely to have been their first predators. Mammals are conservative in the structures of the brain that are involved in vigilance, fear, and learning and memory associated with fearful stimuli, e.g., predators. Some of these areas have expanded in primates and are more strongly connected to visual systems. However, primates vary in the extent of brain expansion. This variation is coincident with variation in evolutionary co-existence with the more recently evolved venomous snakes. Malagasy prosimians have never co-existed with venomous snakes, New World monkeys (platyrrhines) have had interrupted co-existence with venomous snakes, and Old World monkeys and apes (catarrhines) have had continuous co-existence with venomous snakes. The koniocellular visual pathway, arising from the retina and connecting to the lateral geniculate nucleus, the superior colliculus, and the pulvinar, has expanded along with the parvocellular pathway, a visual pathway that is involved with color and object recognition. I suggest that expansion of these pathways co-occurred, with the koniocellular pathway being crucially involved (among other tasks) in pre-attentional visual detection of fearful stimuli, including snakes, and the parvocellular pathway being involved (among other tasks) in protecting the brain from increasingly greater metabolic demands to evolve the neural capacity to detect such stimuli quickly. A diet that included fruits or nectar (though not to the exclusion of arthropods), which provided sugars as a neuroprotectant, may have been a required preadaptation for the expansion of such metabolically active brains. Taxonomic differences in evolutionary exposure to venomous snakes are associated with similar taxonomic differences in rates of evolution in cytochrome oxidase genes and in the metabolic activity of cytochrome oxidase proteins in at least some visual areas in the brains of primates. Raptors that specialize in eating snakes have larger eyes and greater binocularity than more generalized raptors, and provide non-mammalian models for snakes as a selective pressure on primate visual systems. These models, along with evidence from paleobiogeography, neuroscience, ecology, behavior, and immunology, suggest that the evolutionary arms race begun by constrictors early in mammalian evolution continued with venomous snakes. Whereas other mammals responded by evolving physiological resistance to snake venoms, anthropoids responded by enhancing their ability to detect snakes visually before the strike.

Rapid Adaptive Evolution of the Tumor Suppressor Gene Pten in an Insect Lineage

The Pten gene was initially identified in humans as a tumor suppressor. It has since been shown to play important roles in the control of cell size, cell motility, apoptosis, and organ size, and it has also been implicated in aging. Pten is highly conserved among organisms as diverse as nematodes, insects, and vertebrates. In contrast, a phylogenetic analysis by maximum likelihood of a 133-amino acid region showed an average nonsynonymous-to-synonymous rate ratio of 10.4 for Pten in the lineage leading to parasitoid wasps of the Nasonia genus, indicating very strong positive selection. A previous study identified Pten as a potential QTL candidate gene for differences in male wing size in Nasonia. Most of the amino acid replacements that occurred in the Nasonia lineage cluster in a small region of the protein surface, suggesting that they might be involved in an interaction between Pten and another protein. The phenotypic changes due to Pten are not yet known, although it is not associated with known differences in male wing size. Introgression of Pten from one species to another does affect longevity, but a causal relationship is not established.

Asymmetrical Evolution of Cytochrome bd Subunits

Functionally linked genes generally evolve at similar rates and the knowledge of this particular feature of genomic evolution has been used as the basis for the phylogenetic profiling method. We illustrate here an exception to this rule in the evolution of the cytochrome bd complex. This is a two-component oxidase complex, with the subunits I and II known to be widely present in bacteria. The subunits within the cytochrome bd complex are under the same evolutionary pressure and most likely behave in the same evolutionary manner. However, the sequence similarity of genes encoding subunit II varies considerably across species. Genes encoding subunit II evolve 1.2 times faster on most of the branches of their phylogeny than subunit I genes. Furthermore, the genes encoding subunit II in Oceanobacillus iheyensis, Bacillus halodurans, and Staphylococcus species do not have detectable homologues within E. coli due to their large divergence. Together, the two subunits of cytochrome bd reveal an interesting example of an asymmetric pattern of evolutionary change.

Rapid nonsynonymous evolution of the iron-sulfur protein in anthropoid primates

Cytochrome c (CYC) and 9 of the 13 subunits of cytochrome c oxidase (complex IV; COX) were previously shown to have accelerated rates of nonsynonymous substitution in anthropoid primates. Cytochrome b, the mtDNA encoded subunit of ubiquinol-cytochrome c reductase (complex III), also showed an accelerated nonsynonymous substitution rate in anthropoid primates but rate information about the nuclear encoded subunits of complex III has been lacking. We now report that phylogenetic and relative rates analysis of a nuclear encoded catalytically active subunit of complex III, the iron-sulfur protein (ISP), shows an accelerated rate of amino acid replacement similar to cytochrome b. Because both ISP and subunit 9, whose function is not directly related to electron transport, are produced by cleavage into two subunits of the initial translation product of a single gene, it is probable that these two subunits of complex III have essentially identical underlying rates of mutation. Nevertheless, we find that the catalytically active ISP has an accelerated rate of amino acid replacement in anthropoid primates whereas the catalytically inactive subunit 9 does not.

Adaptive evolution of the human fatty acid synthase gene: support for the cancer selection and fat utilization hypotheses?

Cancer may act as the etiological agent for natural selection in some genes. This selective pressure would act to reduce the success of neoplastic lineages over normal cell lineages in individuals of reproductive age. In addition, human's relatively larger brain and longer lifespan may have also acted as a selective force requiring new genotypes. One of the most important proteins in both processes is the fatty acid synthase (FAS) gene involved in fatty acid biosynthesis. A variety of other proteins, including PTEN, MAPK1, SREBP1, SREBP2 and PI are also involved in the regulation of fatty acid biosynthesis. We have specifically analysed variability in selective pressure across all these genes in human, mouse and other vertebrates. We have found that the FAS gene alone has signatures indicative of adaptive evolution. We did not find any signatures of adaptive evolution in any of the other proteins. In the FAS gene, we have detected an excess of non-synonymous over synonymous substitutions in approximately 6% of sites in the human lineage. Contrastingly, the substitution process at these sites in other available vertebrates and mammals indicates strong purifying selection. This is likely to reflect a functional shift in human FAS and correlates well with previously observed changes in FAS biochemical activities. We speculate that the role played by FAS either in cancer development or in human brain development has created this selective pressure, although we cannot rule out the various other functions of FAS.