Genetic evidence and historical theories of the Asian and African origins of the present Malagasy population.

The origin of the Malagasy population has been a subject of speculation since the 16th century. Contributions of African, Asian, Indian, Melanesian, Arabic and Persian populations have been suggested based on physical and cultural anthropology, oral tradition, linguistics and later also by archaeology. In the mid-20th century, increased knowledge of heredity rules and technical progress enabled the identification of African and Asian populations as main contributors. Recent access to the genomic landscape of Madagascar demonstrated pronounced regional variability in the relative contributions of these two ancestries, yet with significant presence of both African and Asian components throughout Madagascar. This article reviews the extent to which genetic results have settled historical questions concerning the origin of the Malagasy population. After an overview of the early literature, the genetic results of the 20th and 21th centuries are discussed and then complemented by the latest results in genome-wide analyses. While there is still much uncertainty regarding when, how and the circumstances under which the ancestors of the modern Malagasy population arrived on the island, we propose a scenario based on historical texts and genomic results.

Genomic landscape of human diversity across Madagascar.

Although situated ∼400 km from the east coast of Africa, Madagascar exhibits cultural, linguistic, and genetic traits from both Southeast Asia and Eastern Africa. The settlement history remains contentious; we therefore used a grid-based approach to sample at high resolution the genomic diversity (including maternal lineages, paternal lineages, and genome-wide data) across 257 villages and 2,704 Malagasy individuals. We find a common Bantu and Austronesian descent for all Malagasy individuals with a limited paternal contribution from Europe and the Middle East. Admixture and demographic growth happened recently, suggesting a rapid settlement of Madagascar during the last millennium. However, the distribution of African and Asian ancestry across the island reveals that the admixture was sex biased and happened heterogeneously across Madagascar, suggesting independent colonization of Madagascar from Africa and Asia rather than settlement by an already admixed population. In addition, there are geographic influences on the present genomic diversity, independent of the admixture, showing that a few centuries is sufficient to produce detectable genetic structure in human populations.

Modeling the respiratory chain complexes with biothermokinetic equations - the case of complex I.

The mitochondrial respiratory chain plays a crucial role in energy metabolism and its dysfunction is implicated in a wide range of human diseases. In order to understand the global expression of local mutations in the rate of oxygen consumption or in the production of adenosine triphosphate (ATP) it is useful to have a mathematical model in which the changes in a given respiratory complex are properly modeled. Our aim in this paper is to provide thermodynamics respecting and structurally simple equations to represent the kinetics of each isolated complexes which can, assembled in a dynamical system, also simulate the behavior of the respiratory chain, as a whole, under a large set of different physiological and pathological conditions. On the example of the reduced nicotinamide adenine dinucleotide (NADH)-ubiquinol-oxidoreductase (complex I) we analyze the suitability of different types of rate equations. Based on our kinetic experiments we show that very simple rate laws, as those often used in many respiratory chain models, fail to describe the kinetic behavior when applied to a wide concentration range. This led us to adapt rate equations containing the essential parameters of enzyme kinetic, maximal velocities and Henri-Michaelis-Menten like-constants (KM and KI) to satisfactorily simulate these data.

From in silico to in spectro kinetics of respiratory complex I.

An enzyme's activity is the consequence of its structure. The stochastic approach we developed to study the functioning of the respiratory complexes is based upon their 3D structure and their physical and chemical properties. Consequently it should predict their kinetic properties. In this paper we compare the predictions of our stochastic model derived for the complex I with a number of experiments performed with a large range of complex I substrates and products. A good fit was found between the experiments and the prediction of our stochastic approach. We show that, due to the spatial separation of the two half redox reactions (NADH/NAD and Q/QH(2)), the kinetics cannot necessarily obey a simple mechanism (ordered or ping-pong for instance). A plateau in the kinetics is observed at high substrates concentrations, well evidenced in the double reciprocal plots, which is explained by the limiting rate of quinone reduction as compared with the oxidation of NADH at the other end of complex I. Moreover, we show that the set of the seven redox reactions in between the two half redox reactions (NADH/NAD and Q/QH(2)) acts as an electron buffer. An inhibition of complex I activity by quinone is observed at high concentration of this molecule, which cannot be explained by a simple stochastic model based on the known structure. We hypothesize that the distance between the catalytic site close to N2 (iron/sulfur redox center that transfers electrons to quinone) and the membrane forces the quinone/quinol to take several positions in between these sites. We represent these possible positions by an extra site necessarily occupied by the quinone/quinol molecules on their way to the redox site. With this hypothesis, we are able to fit the kinetic experiments over a large range of substrates and products concentrations. The slow rate constants derived for the transition between the two sites could be an indication of a conformational change of the enzyme during the quinone/quinol movement. This article is part of a Special Issue entitled: 17th European Bioenergetics Conference (EBEC 2012).

Mitochondrial energetic metabolism-some general principles.

In nonphotosynthetic organisms, mitochondria are the power plant of the cell, emphasizing their great potentiality for adenosine triphosphate (ATP) synthesis from the redox span between nutrients and oxygen. Also of great importance is their role in the maintenance of the cell redox balance. Even though crystallographic structures of respiratory complexes, ATP synthase, and ATP/adenosine diphosphate (ADP) carrier are now quite well known, the coupling between ATP synthesis and cell redox state remains a controversial issue. In this review, we will present some of the processes that allow a modular coupling between ATP synthesis and redox state. Furthermore, we will present some theoretical approaches of this highly integrated system.

Canadian Occupational Performance Measure: Benefits and Limitations Highlighted Using the Delphi Method and Principal Component Analysis.

Introduction: The objective of this study was to establish a baseline of current use in practice of the Canadian Occupational Performance Measure (COPM) by consulting 33 expert French occupational therapists, who trained in this method between 2012 and 2017 and use of the COPM with their clients. The areas of health intervention are pediatrics, psychiatry, neurology, and geriatrics. An email invitation to participate in the research was therefore sent to 113 occupational therapists. We received 33 responses. Methods: A novel mixed method study combined a Delphi method with a lexical analysis of experts' responses and principal component analysis (PCA). Results: In the last Delphi round, the consensus of the expert group was attained on 31 benefits and 1 limitation, confirming the generally positive influence of the COPM in French health services. Discussion. The COPM was clearly identified as a tool that supports occupational therapy clinical reasoning, facilitates team decision-making for care pathways, and enables people with disabilities and health conditions to make decisions for their care. Conclusion: The Delphi study revealed that the COPM appears to be well adapted to French culture and should be widely incorporated into preregistration training.

The loss of biodiversity in Madagascar is contemporaneous with major demographic events.

Only 400 km off the coast of East Africa, the island of Madagascar is one of the last large land masses to have been colonized by humans. While many questions surround the human occupation of Madagascar, recent studies raise the question of human impact on endemic biodiversity and landscape transformation. Previous genetic and linguistic analyses have shown that the Malagasy population has emerged from an admixture that happened during the last millennium, between Bantu-speaking African populations and Austronesian-speaking Asian populations. By studying the sharing of chromosome segments between individuals (IBD determination), local ancestry information, and simulated genetic data, we inferred that the Malagasy ancestral Asian population was isolated for more than 1,000 years with an effective size of just a few hundred individuals. This isolation ended around 1,000 years before present (BP) by admixture with a small African population. Around the admixture time, there was a rapid demographic expansion due to intrinsic population growth of the newly admixed population, which coincides with extensive changes in Madagascar's landscape and the extinction of all endemic large-bodied vertebrates. Therefore, our approach can provide new insights into past human demography and associated impacts on ecosystems.

Measurement of kinetic parameters of human platelet DNA polymerase gamma.

Synthesis of mitochondrial DNA is performed by DNA polymerase gamma. Mutations in POLG, the gene encoding the catalytic subunit of DNA polymerase gamma, are a major cause of neurological disease. A large proportion of patients carry rare nucleotide substitutions leading to single amino acid changes. Confirming that these replacements are pathogenic can be problematic without biochemical evidence. Here, we provide a hands-on protocol for an in vitro kinetic assay of DNA polymerase gamma which allows assessment of the K(m) and V(max) for the incoming nucleotide of the polymerization reaction. To avoid measurement of contaminating nuclear DNA polymerases, platelet extracts are used since platelets do not contain a nucleus. Moreover, platelets have the advantage of being obtainable relatively non-invasively. Polymerization activity is determined by measurement of the incorporation of radioactive thymidine 5'-triphosphate (dTTP) on the homopolymeric RNA substrate poly(rA).oligo(dT)(12-18). To further minimize nuclear DNA polymerase activity, aphidicolin, an inhibitor of most nuclear DNA polymerases, is included in the reaction. In addition, reactions are carried out in the absence and presence of the competitive inhibitor of DNA polymerase gamma, 2',3'-dideoxythymidine 5'-triphosphate (ddTTP), to allow calculation of the ddTTP-sensitive incorporation. With this method, platelets from healthy control subjects extracted with 3% Triton X-100 showed a K(m) for dTTP of 1.42 microM and a V(max) of 0.83 pmol min(-1)mg(-1).

African Gene Flow Reduces Beta-Ionone Anosmia/Hyposmia Prevalence in Admixed Malagasy Populations.

While recent advances in genetics make it possible to follow the genetic exchanges between populations and their phenotypic consequences, the impact of the genetic exchanges on the sensory perception of populations has yet to be explored. From this perspective, the present study investigated the consequences of African gene flow on odor perception in a Malagasy population with a predominantly East Asian genetic background. To this end, we combined psychophysical tests with genotype data of 235 individuals who were asked to smell the odorant molecule beta-ionone (ßI). Results showed that in this population the ancestry of the OR5A1 gene significantly influences the ability to detect ßI. At the individual level, African ancestry significantly protects against specific anosmia/hyposmia due to the higher frequency of the functional gene (OR ratios = 14, CI: 1.8-110, p-value = 0.012). At the population level, African introgression decreased the prevalence of specific anosmia/hyposmia to this odorous compound. Taken together, these findings validate the conjecture that in addition to cultural exchanges, genetic transfer may also influence the sensory perception of the population in contact.

Virtual mitochondrion: towards an integrated model of oxidative phosphorylation complexes and beyond.

The modelling of OXPHOS (oxidative phosphorylation) in order to integrate all kinetic and thermodynamic aspects of chemiosmotic theory has a long history. We briefly review this history and show how new ways of modelling are required to integrate a local model of the individual respiratory complexes into a global model of OXPHOS and, beyond that, into a reliable overall model of central metabolism.

Smell and taste changes are early indicators of the COVID-19 pandemic and political decision effectiveness.

In response to the COVID-19 pandemic, many governments have taken drastic measures to avoid an overflow of intensive care units. Accurate metrics of disease spread are critical for the reopening strategies. Here, we show that self-reports of smell/taste changes are more closely associated with hospital overload and are earlier markers of the spread of infection of SARS-CoV-2 than current governmental indicators. We also report a decrease in self-reports of new onset smell/taste changes as early as 5 days after lockdown enforcement. Cross-country comparisons demonstrate that countries that adopted the most stringent lockdown measures had faster declines in new reports of smell/taste changes following lockdown than a country that adopted less stringent lockdown measures. We propose that an increase in the incidence of sudden smell and taste change in the general population may be used as an indicator of COVID-19 spread in the population.

Mitochondrial energetic metabolism: a simplified model of TCA cycle with ATP production.

Mitochondria play a central role in cellular energetic metabolism. The essential parts of this metabolism are the tricarboxylic acid (TCA) cycle, the respiratory chain and the adenosine triphosphate (ATP) synthesis machinery. Here a simplified model of these three metabolic components with a limited set of differential equations is presented. The existence of a steady state is demonstrated and results of numerical simulations are presented. The relevance of a simple model to represent actual in vivo behavior is discussed.

Strong selection during the last millennium for African ancestry in the admixed population of Madagascar.

While admixed populations offer a unique opportunity to detect selection, the admixture in most of the studied populations occurred too recently to produce conclusive signals. By contrast, Malagasy populations originate from admixture between Asian and African populations that occurred ~27 generations ago, providing power to detect selection. We analyze local ancestry across the genomes of 700 Malagasy and identify a strong signal of recent positive selection, with an estimated selection coefficient >0.2. The selection is for African ancestry and affects 25% of chromosome 1, including the Duffy blood group gene. The null allele at this gene provides resistance to Plasmodium vivax malaria, and previous studies have suggested positive selection for this allele in the Malagasy population. This selection event also influences numerous other genes implicated in immunity, cardiovascular diseases, and asthma and decreases the Asian ancestry genome-wide by 10%, illustrating the role played by selection in recent human history.

Comprehensive mathematical model of oxidative phosphorylation valid for physiological and pathological conditions.

We developed a mathematical model of oxidative phosphorylation (OXPHOS) that allows for a precise description of mitochondrial function with respect to the respiratory flux and the ATP production. The model reproduced flux-force relationships under various experimental conditions (state 3 and 4, uncoupling, and shortage of respiratory substrate) as well as time courses, exhibiting correct P/O ratios. The model was able to reproduce experimental threshold curves for perturbations of the respiratory chain complexes, the F1 F0 -ATP synthase, the ADP/ATP carrier, the phosphate/OH carrier, and the proton leak. Thus, the model is well suited to study complex interactions within the OXPHOS system, especially with respect to physiological adaptations or pathological modifications, influencing substrate and product affinities or maximal catalytic rates. Moreover, it could be a useful tool to study the role of OXPHOS and its capacity to compensate or enhance physiopathologies of the mitochondrial and cellular energy metabolism.

Modeling of mitochondria bioenergetics using a composable chemiosmotic energy transduction rate law: theory and experimental validation.

Mitochondrial bioenergetic processes are central to the production of cellular energy, and a decrease in the expression or activity of enzyme complexes responsible for these processes can result in energetic deficit that correlates with many metabolic diseases and aging. Unfortunately, existing computational models of mitochondrial bioenergetics either lack relevant kinetic descriptions of the enzyme complexes, or incorporate mechanisms too specific to a particular mitochondrial system and are thus incapable of capturing the heterogeneity associated with these complexes across different systems and system states. Here we introduce a new composable rate equation, the chemiosmotic rate law, that expresses the flux of a prototypical energy transduction complex as a function of: the saturation kinetics of the electron donor and acceptor substrates; the redox transfer potential between the complex and the substrates; and the steady-state thermodynamic force-to-flux relationship of the overall electro-chemical reaction. Modeling of bioenergetics with this rate law has several advantages: (1) it minimizes the use of arbitrary free parameters while featuring biochemically relevant parameters that can be obtained through progress curves of common enzyme kinetics protocols; (2) it is modular and can adapt to various enzyme complex arrangements for both in vivo and in vitro systems via transformation of its rate and equilibrium constants; (3) it provides a clear association between the sensitivity of the parameters of the individual complexes and the sensitivity of the system's steady-state. To validate our approach, we conduct in vitro measurements of ETC complex I, III, and IV activities using rat heart homogenates, and construct an estimation procedure for the parameter values directly from these measurements. In addition, we show the theoretical connections of our approach to the existing models, and compare the predictive accuracy of the rate law with our experimentally fitted parameters to those of existing models. Finally, we present a complete perturbation study of these parameters to reveal how they can significantly and differentially influence global flux and operational thresholds, suggesting that this modeling approach could help enable the comparative analysis of mitochondria from different systems and pathological states. The procedures and results are available in Mathematica notebooks at http://www.igb.uci.edu/tools/sb/mitochondria-modeling.html.

Mutation rate switch inside Eurasian mitochondrial haplogroups: impact of selection and consequences for dating settlement in Europe.

R-lineage mitochondrial DNA represents over 90% of the European population and is significantly present all around the planet (North Africa, Asia, Oceania, and America). This lineage played a major role in migration "out of Africa" and colonization in Europe. In order to determine an accurate dating of the R lineage and its sublineages, we analyzed 1173 individuals and complete mtDNA sequences from Mitomap. This analysis revealed a new coalescence age for R at 54.500 years, as well as several limitations of standard dating methods, likely to lead to false interpretations. These findings highlight the association of a striking under-accumulation of synonymous mutations, an over-accumulation of non-synonymous mutations, and the phenotypic effect on haplogroup J. Consequently, haplogroup J is apparently not a Neolithic group but an older haplogroup (Paleolithic) that was subjected to an underestimated selective force. These findings also indicated an under-accumulation of synonymous and non-synonymous mutations localized on coding and non-coding (HVS1) sequences for haplogroup R0, which contains the major haplogroups H and V. These new dates are likely to impact the present colonization model for Europe and confirm the late glacial resettlement scenario.

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.