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.

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.

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.

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 F₁ F₀ -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 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).

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.

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.

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).

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.