Mitochondrial DNA associations with East Asian metabolic syndrome.

Mitochondrial dysfunction has repeatedly been reported associated with type 2 diabetes mellitus (T2DM) and metabolic syndrome (MS), as have mitochondrial DNA (mtDNA) tRNA and duplication mutations and mtDNA haplogroup lineages. We identified 19 Taiwanese T2DM and MS pedigrees from Taiwan, with putative matrilineal transmission, one of which harbored the pathogenic mtDNA tRNALeu(UUR) nucleotide (nt) 3243A>G mutation on the N9a3 haplogroup background. We then recruited three independent Taiwanese cohorts, two from Taipei (N = 498, mean age 52 and N = 1002, mean age 44) and one from a non-urban environment (N = 501, mean age 57). All three cohorts were assessed for an array of metabolic parameters, their mtDNA haplogroups determined, and the haplogroups correlated with T2DM/MS phenotypes. Logistic regression analysis revealed that mtDNA haplogroups D5, F4, and N9a conferred T2DM protection, while haplogroups F4 and N9a were risk factors for hypertension (HTN), and F4 was a risk factor for obesity (OB). Additionally, the 5263C>T (ND2 A165V) variant commonly associated with F4 was associated with hypertension (HTN). Cybrids were prepared with macro-haplogroup N (defined by variants m.ND3 10398A (114T) and m.ATP6 8701A (59T)) haplogroups B4 and F1 mtDNAs and from macro-haplogroup M (variants m.ND3 10398G (114A) and m.ATP6 8701G (59A)) haplogroup M9 mtDNAs. Additionally, haplogroup B4 and F1 cybrids were prepared with and without the mtDNA variant in ND1 3394T>C (Y30H) reported to be associated with T2DM. Assay of mitochondria complex I in these cybrids revealed that macro-haplogroup N cybrids had lower activity than M cybrids, that haplogroup F cybrids had lower activity than B4 cybrids, and that the ND1 3394T>C (Y30H) variant reduced complex I on both the B4 and F1 background but with very different cumulative effects. These data support the hypothesis that functional mtDNA variants may contribute to the risk of developing T2DM and MS.

Concurrent action of purifying selection and gene conversion results in extreme conservation of the major stress-inducible Hsp70 genes in mammals.

Several evolutionary mechanisms alter the fate of mutations and genes within populations based on their exhibited functional effects. To understand the underlying mechanisms involved in the evolution of the cellular stress response, a very conserved mechanism in the course of organismal evolution, we studied the patterns of natural genetic variation and functional consequences of polymorphisms of two stress-inducible Hsp70 genes. These genes, HSPA1A and HSPA1B, are major orchestrators of the cellular stress response and are associated with several human diseases. Our phylogenetic analyses revealed that the duplication of HSPA1A and HSPA1B originated in a lineage proceeding to placental mammals, and henceforth they remained in conserved synteny. Additionally, analyses of synonymous and non-synonymous changes suggest that purifying selection shaped the HSPA1 gene diversification, while gene conversion resulted in high sequence conservation within species. In the human HSPA1-cluster, the vast majority of mutations are synonymous and specific genic regions are devoid of mutations. Furthermore, functional characterization of several human polymorphisms revealed subtle differences in HSPA1A stability and intracellular localization. Collectively, the observable patterns of HSPA1A-1B variation describe an evolutionary pattern, in which purifying selection and gene conversion act simultaneously and conserve a major orchestrator of the cellular stress response.

Association Between Mitochondrial DNA Haplogroup Variation and Autism Spectrum Disorders.

Importance: Autism spectrum disorders (ASD) are characterized by impairments in social interaction, communication, and repetitive or restrictive behavior. Although multiple physiologic and biochemical studies have reported defects in mitochondrial oxidative phosphorylation in patients with ASD, the role of mitochondrial DNA (mtDNA) variation has remained relatively unexplored. Objective: To assess what impact mitochondrial lineages encompassing ancient mtDNA functional polymorphisms, termed haplogroups, have on ASD risk. Design, Setting, and Participants: In this cohort study, individuals with autism and their families were studied using the Autism Genetic Resource Exchange cohort genome-wide association studies data previously generated at the Children's Hospital of Philadelphia. From October 2010 to January 2017, we analyzed the data and used the mtDNA single-nucleotide polymorphisms interrogated by the Illumina HumanHap 550 chip to determine the mtDNA haplogroups of the individuals. Taking into account the familial structure of the Autism Genetic Resource Exchange data, we then determined whether the mtDNA haplogroups correlate with ASD risk. Main Outcomes and Measures: Odds ratios of mitochondrial haplogroup as predictors of ASD risk. Results: Of 1624 patients with autism included in this study, 1299 were boys (80%) and 325 were girls (20%). Families in the Autism Genetic Resource Exchange collection (933 families, encompassing 4041 individuals: 1624 patients with ASD and 2417 healthy parents and siblings) had been previously recruited in the United States with no restrictions on age, sex, race/ethnicity, or socioeconomic status. Relative to the most common European haplogroup HHV, European haplogroups I, J, K, O-X, T, and U were associated with increased risk of ASD, as were Asian and Native American haplogroups A and M, with odds ratios ranging from 1.55 (95% CI, 1.16-2.06) to 2.18 (95% CI, 1.59-3) (adjusted P < .04). Hence, mtDNA haplogroup variation is an important risk factor for ASD. Conclusions and Relevance: Because haplogroups I, J, K, O-X, T, and U encompass 55% of the European population, mtDNA lineages must make a significant contribution to overall ASD risk.

Functional diversification and specialization of cytosolic 70-kDa heat shock proteins.

A fundamental question in molecular evolution is how protein functional differentiation alters the ability of cells and organisms to cope with stress and survive. To answer this question we used two paralogous Hsp70s from mouse and explored whether these highly similar cytosolic molecular chaperones, which apart their temporal expression have been considered functionally interchangeable, are differentiated with respect to their lipid-binding function. We demonstrate that the two proteins bind to diverse lipids with different affinities and therefore are functionally specialized. The observed lipid-binding patterns may be related with the ability of both Hsp70s to induce cell death by binding to a particular plasma-membrane lipid, and the potential of only one of them to promote cell survival by binding to a specific lysosomal-membrane lipid. These observations reveal that two seemingly identical proteins differentially modulate cellular adaptation and survival by having acquired specialized functions via sequence divergence. Therefore, this study provides an evolutionary paradigm, where promiscuity, specificity, sub- and neo-functionalization orchestrate one of the most conserved systems in nature, the cellular stress-response.

Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease.

The unorthodox genetics of the mtDNA is providing new perspectives on the etiology of the common "complex" diseases. The maternally inherited mtDNA codes for essential energy genes, is present in thousands of copies per cell, and has a very high mutation rate. New mtDNA mutations arise among thousands of other mtDNAs. The mechanisms by which these "heteroplasmic" mtDNA mutations come to predominate in the female germline and somatic tissues is poorly understood, but essential for understanding the clinical variability of a range of diseases. Maternal inheritance and heteroplasmy also pose major challengers for the diagnosis and prevention of mtDNA disease.

mtDNA Variation and Analysis Using Mitomap and Mitomaster.

The Mitomap database of human mitochondrial DNA (mtDNA) information has been an important compilation of mtDNA variation for researchers, clinicians and genetic counselors for the past twenty-five years. The Mitomap protocol shows how users may look up human mitochondrial gene loci, search for public mitochondrial sequences, and browse or search for reported general population nucleotide variants as well as those reported in clinical disease. Within Mitomap is the powerful sequence analysis tool for human mitochondrial DNA, Mitomaster. The Mitomaster protocol gives step-by-step instructions showing how to submit sequences to identify nucleotide variants relative to the rCRS, to determine the haplogroup, and to view species conservation. User-supplied sequences, GenBank identifiers and single nucleotide variants may be analyzed.

Heteroplasmy of mouse mtDNA is genetically unstable and results in altered behavior and cognition.

Maternal inheritance of mtDNA is the rule in most animals, but the reasons for this pattern remain unclear. To investigate the consequence of overriding uniparental inheritance, we generated mice containing an admixture (heteroplasmy) of NZB and 129S6 mtDNAs in the presence of a congenic C57BL/6J nuclear background. Analysis of the segregation of the two mtDNAs across subsequent maternal generations revealed that proportion of NZB mtDNA was preferentially reduced. Ultimately, this segregation process produced NZB-129 heteroplasmic mice and their NZB or 129 mtDNA homoplasmic counterparts. Phenotypic comparison of these three mtDNA lines demonstrated that the NZB-129 heteroplasmic mice, but neither homoplasmic counterpart, had reduced activity, food intake, respiratory exchange ratio; accentuated stress response; and cognitive impairment. Therefore, admixture of two normal but different mouse mtDNAs can be genetically unstable and can produce adverse physiological effects, factors that may explain the advantage of uniparental inheritance of mtDNA.

Mitochondrial DNA variant associated with Leber hereditary optic neuropathy and high-altitude Tibetans.

The distinction between mild pathogenic mtDNA mutations and population polymorphisms can be ambiguous because both are homoplasmic, alter conserved functions, and correlate with disease. One possible explanation for this ambiguity is that the same variant may have different consequences in different contexts. The NADH dehydrogenase subunit 1 (ND1) nucleotide 3394 T > C (Y30H) variant is such a case. This variant has been associated with Leber hereditary optic neuropathy and it reduces complex I activity and cellular respiration between 7% and 28% on the Asian B4c and F1 haplogroup backgrounds. However, complex I activity between B4c and F1 mtDNAs, which harbor the common 3394T allele, can also differ by 30%. In Asia, the 3394C variant is most commonly associated with the M9 haplogroup, which is rare at low elevations but increases in frequency with elevation to an average of 25% of the Tibetan mtDNAs (odds ratio = 23.7). In high-altitude Tibetan and Indian populations, the 3394C variant occurs on five different macrohaplogroup M haplogroup backgrounds and is enriched on the M9 background in Tibet and the C4a4 background on the Indian Deccan Plateau (odds ratio = 21.9). When present on the M9 background, the 3394C variant is associated with a complex I activity that is equal to or higher than that of the 3394T variant on the B4c and F1 backgrounds. Hence, the 3394C variant can either be deleterious or beneficial depending on its haplogroup and environmental context. Thus, this mtDNA variant fulfills the criteria for a common variant that predisposes to a "complex" disease.

Phylogenetic Profiles Reveal Structural and Functional Determinants of Lipid-binding.

One of the major challenges in the genomic era is annotating structure/function to the vast quantities of sequence information now available. Indeed, most of the protein sequence database lacks comprehensive annotation, even when experimental evidence exists. Further, within structurally resolved and functionally annotated protein domains, additional functionalities contained in these domains are not apparent. To add further complication, small changes in the amino-acid sequence can lead to profound changes in both structure and function, underscoring the need for rapid and reliable methods to analyze these types of data. Phylogenetic profiles provide a quantitative method that can relate the structural and functional properties of proteins, as well as their evolutionary relationships. Using all of the structurally resolved Src-Homology-2 (SH2) domains, we demonstrate that knowledge-bases can be used to create single-amino acid phylogenetic profiles which reliably annotate lipid-binding. Indeed, these measures isolate the known phosphotyrosine and hydrophobic pockets as integral to lipid-binding function. In addition, we determined that the SH2 domain of Tec family kinases bind to lipids with varying affinity and specificity. Simulating mutations in Bruton's tyrosine kinase (BTK) that cause X-Linked Agammaglobulinemia (XLA) predict that these mutations alter lipid-binding, which we confirm experimentally. In light of these results, we propose that XLA-causing mutations in the SH3-SH2 domain of BTK alter lipid-binding, which could play a causative role in the XLA-phenotype. Overall, our study suggests that the number of lipid-binding proteins is drastically underestimated and, with further development, phylogenetic profiles can provide a method for rapidly increasing the functional annotation of protein sequences.

Origins and evolution of the formin multigene family that is involved in the formation of actin filaments.

In eukaryotes, the assembly and elongation of unbranched actin filaments is controlled by formins, which are long, multidomain proteins. These proteins are important for dynamic cellular processes such as determination of cell shape, cell division, and cellular interaction. Yet, no comprehensive study has been done about the origins and evolution of this gene family. We therefore performed extensive phylogenetic and motif analyses of the formin genes by examining 597 prokaryotic and 53 eukaryotic genomes. Additionally, we used three-dimensional protein structure data in an effort to uncover distantly related sequences. Our results suggest that the formin homology 2 (FH2) domain, which promotes the formation of actin filaments, is a eukaryotic innovation and apparently originated only once in eukaryotic evolution. Despite the high degree of FH2 domain sequence divergence, the FH2 domains of most eukaryotic formins are predicted to assume the same fold and thus have similar functions. The formin genes have experienced multiple taxon-specific duplications and followed the birth-and-death model of evolution. Additionally, the formin genes experienced taxon-specific genomic rearrangements that led to the acquisition of unrelated protein domains. The evolutionary diversification of formin genes apparently increased the number of formin's interacting molecules and consequently contributed to the development of a complex and precise actin assembly mechanism. The diversity of formin types is probably related to the range of actin-based cellular processes that different cells or organisms require. Our results indicate the importance of gene duplication and domain acquisition in the evolution of the eukaryotic cell and offer insights into how a complex system, such as the cytoskeleton, evolved.

Ancient origin of the new developmental superfamily DANGER.

Developmental proteins play a pivotal role in the origin of animal complexity and diversity. We report here the identification of a highly divergent developmental protein superfamily (DANGER), which originated before the emergence of animals (approximately 850 million years ago) and experienced major expansion-contraction events during metazoan evolution. Sequence analysis demonstrates that DANGER proteins diverged via multiple mechanisms, including amino acid substitution, intron gain and/or loss, and recombination. Divergence for DANGER proteins is substantially greater than for the prototypic member of the superfamily (Mab-21 family) and other developmental protein families (e.g., WNT proteins). DANGER proteins are widely expressed and display species-dependent tissue expression patterns, with many members having roles in development. DANGER1A, which regulates the inositol trisphosphate receptor, promotes the differentiation and outgrowth of neuronal processes. Regulation of development may be a universal function of DANGER family members. This family provides a model system to investigate how rapid protein divergence contributes to morphological complexity.

Origin and evolution of the chicken leukocyte receptor complex.

In mammals, the cell surface receptors encoded by the leukocyte receptor complex (LRC) regulate the activity of T lymphocytes and B lymphocytes, as well as that of natural killer cells, and thus provide protection against pathogens and parasites. The chicken genome encodes many Ig-like receptors that are homologous to the LRC receptors. The chicken Ig-like receptor (CHIR) genes are members of a large monophyletic gene family and are organized into genomic clusters, which are in conserved synteny with the mammalian LRC. One-third of CHIR genes encode polypeptide molecules that contain both activating and inhibitory motifs. These genes are present in different phylogenetic groups, suggesting that the primordial CHIR gene could have encoded both types of motifs in a single molecule. In contrast to the mammalian LRC genes, the CHIR genes with similar function (inhibition or activation) are evolutionarily closely related. We propose that, in addition to recombination, single nucleotide substitutions played an important role in the generation of receptors with different functions. Structural models and amino acid analyses of the CHIR proteins reveal the presence of different types of Ig-like domains in the same phylogenetic groups, as well as sharing of conserved residues and conserved changes of residues between different CHIR groups and between CHIRs and LRCs. Our data support the notion that the CHIR gene clusters are regions homologous to the mammalian LRC gene cluster and favor a model of evolution by repeated processes of birth and death (expansion-contraction) of the Ig-like receptor genes.

Retrospective study on poliomyelitis cases in northern Greece.

Here we report the results of a retrospective study on the epidemiological characteristics and genetic relationships of the virus isolates responsible for the last poliomyelitis cases in Greece. The last wild poliomyelitis case in Greece was detected in 1996, and the last vaccine-related strain was isolated in 1998. The whole of Europe, including Greece, is now considered to be polio-free.