Activity of ancestral restriction factors against ancient retroviruses

Analysis of TRIM5α and APOBEC3G genes suggests that these two restriction factors underwent strong positive selection throughout primate evolution. This pressure was possibly imposed by ancient exogenous retroviruses, of which endogenous retroviruses are remnants. Our study aims to assess in vitro the activity of these factors against ancient retroviruses by reconstructing their ancestral gag sequences, as well as the ancestral TRIM5α and APOBEC3G for primates. Based on evolutionary genomics approach, we reconstructed ancestors of the two largest families of human endogenous retroviruses (HERV), namely HERV-K and HERV-H, as well as primate ancestral TRIM5α and APOBEC3G variants. The oldest TRIM5α sequence was the catarhinne TRIM5α, common ancestor of Old World monkeys and hominoids, dated from 25 million years ago (mya). From the oldest, to the youngest, ancestral TRIM5α variants showed less restriction of HIV-1 in vitro [1]. Likewise three ancestral APOBEC3Gs sequences common to hominoids (18 mya), Old World monkeys, and catarhinnes (25 mya) were reconstructed. All ancestral APOBEC3G variants inhibited efficiently HIV-1Δvif in vitro, compared to modern APOBEC3Gs. The ability of Vif proteins (HIV-1, HIV-2, SIVmac and SIVagm) to counteract their activity tallied with the residue 128 on ancestral APOBEC3Gs. Moreover we are attempting to reconstruct older ancestral sequences of both restriction factors by using prosimian orthologue sequences. An infectious onemillion- years-old HERV-KCON previously reconstituted was shown to be resistant to modern TRIM5α and APOBEC3G [2]. Our ancestral TRIM5α and APOBEC3G variants were inactive against HERV-KCON. Besides we reconstructed chimeric HERV-K bearing ancestral capsids (up to 7 mya) that resulted in infectious viruses resistant to modern and ancestral TRIM5α. Likewise HERV-K viruses bearing ancestral nucleocapsids will be tested for ancestral and modern APOBEC3G restriction. In silico reconstruction and structural modeling of ancestral HERV-H capsids resulted in structures homologous to that of the gammaretrovirus MLV. Thus we are attempting to construct chimeric MLV virus bearing HERV-H ancestral capsids. These chimeric ancestral HERVs will be tested for infectivity and restriction by ancestral TRIM5α. Similarly chimeric MLV viruses bearing ancestral HERV-H nucleocapsids will be reconstructed and tested for APOBEC3G restriction.

Evolution of C(4) phosphoenolpyruvate carboxykinase in grasses, from genotype to phenotype.

C(4) photosynthesis is an adaptation over the classical C(3) pathway that has evolved multiple times independently. These convergences are accompanied by strong variations among the independent C(4) lineages. The decarboxylating enzyme used to release CO(2) around Rubisco particularly differs between C(4) species, a criterion used to distinguish three distinct biochemical C(4) subtypes. The phosphoenolpyruvate carboxykinase (PCK) serves as a primary decarboxylase in a minority of C(4) species. This enzyme is also present in C(3) plants, where it is responsible for nonphotosynthetic functions. The genetic changes responsible for the evolution of C(4)-specific PCK are still unidentified. Using phylogenetic analyses on PCK sequences isolated from C(3) and C(4) grasses, this study aimed at resolving the evolutionary history of C(4)-specific PCK enzymes. Four independent evolutions of C(4)-PCK were shown to be driven by positive selection, and nine C(4)-adaptive sites underwent parallel genetic changes in different C(4) lineages. C(4)-adaptive residues were also observed in C(4) species from the nicotinamide adenine dinucleotide phosphate-malic enzyme (NADP-ME) subtype and particularly in all taxa where a PCK shuttle was previously suggested to complement the NADP-ME pathway. Acquisitions of C(4)-specific PCKs were mapped on a species tree, which revealed that the PCK subtype probably appeared at the base of the Chloridoideae subfamily and was then recurrently lost and secondarily reacquired at least three times. Linking the genotype to subtype phenotype shed new lights on the evolutionary transitions between the different C(4) subtypes.

Pathogenesis-related proteins and the jasmonate signaling pathway

Résumé Etant une importante source d'énergie, les plantes sont constamment attaquées par des pathogènes. Ne pouvant se mouvoir, elles ont développé des systèmes de défense sophistiqués afin de lutter contre ces prédateurs. Parmi ces systèmes, les voies de signalisation mettant en jeu des éliciteurs endog8nes tels que les jasmonates permettent d'induire la production de protéines de défense telles que les protéines dites "liées à la pathogénèse". Les gènes codant pour ces protéines appartiennent à des familles multigéniques. Le premier but de cette thèse est d'évaluer le nombre de ces gènes dans le génome d'Arabidopsis thaliana et d'estimer la part de ce système de défense, dépendant de la voie de signalisation des jasmonates. Nous avons défini un cluster de seulement 1S gènes sur 266, "liés à la pathogénèse", exclusivement régulés par les jasmonates. De multiples membres des familles des lectines de type jacaline et des inhibiteurs de trypsines semblent dépendre du jasmonate. Présente dans tous les systèmes immunitaires des eucaryotes, la famille des défensines est une famille très intéressante. Chez Arabidopsis thaliana, 317 protéines similaires aux défensines ont été définies, cependant seulement 15 défensines (PDF) sont bien annotées. Ces 15 défensines sont séparées en deux groupes dont un semble avoir évolué plus récemment. Le second but de cette thèse est d'étudier ce groupe de défensines à l'aide de la bioinformatique et des techniques de biologie moléculaire (gêne rapporteur, PCR en temps réel). Nous avons montré que ce groupe contenait une défensine acide intéressante, PDF1.5, qui semblait avoir subi une sélection positive. Cette protéine n'avait encore jamais été étudiée. Contrairement à ce que nous pensions, nous avons établi que cette protéine pouvait avoir une activité biologique liée à la défense. Ce travail de thèse a permis de préciser le nombre de gènes "liées à la pathogénèse" induits par la voie des jasmonates et d'apporter des éléments de réponse sur la question de la redondance des gènes de défense. En conclusion, même si de nombreuses familles de gènes intervenant dans la défense sont bien définies chez Arabidopsis, il reste encore de nombreuses études à faire sur chacun de ces membres. Abstract Being an important source of energy, plants are constantly attacked by herbivores and pathogens. As sessile organisms, they have developed sophisticated defense responses to cope with attack. Among these responses, signalling pathways, using endogenous elicitors including jasmonates (JA), allow the plant to induce the production of defense proteins such as pathogenesis-related (PR) proteins. The genes encoding these proteins belong to multigenic families. The first goal of this thesis was to evaluate the number of PR genes in the genome of Arabidopsis thaliana and estimate how much of this plant defense system was dependent on the jasmonate signaling pathway in leaves. Surprisingly a cluster of only 1S genes out of 2ó6 PR genes was exclusively regulated by JA. Multiple members of the jacalin lectin and trypsin inhibitor gene families were shown to be regulated by JA. Present in all eukaryotic immune systems, defensins are an attractive PR family to study. In Arabidopsis thaliana, 317 defensin-related proteins have been found but just 1S defensins (i.e. PDF family) are well annotated. These defensins are split into 2 groups. One of these groups may have appeared and diversified recently. The second goal of this thesis was to study this defensin gene group combining bioinformatic, reporter gene and quantitative PCR techniques. We have shown that this group contains an interesting acidic defensin, PDF1.S, which seems to have undergone positive selection. No information was known on this protein. We have established that this protein may have a biological activity in plant defense. This thesis allowed us to define the number of PR genes induced by the jasmonate pathway and gave initial leads to explain the redundancy of the PR genes in the genome of Arabidopsis. In conclusion, even if many defense gene families are already defined in the Arabidopsis genome, much work remains to be done on individual members.

Association of hepatitis C virus NS5B variants with resistance to new antiviral drugs among untreated patients

Mutations located in the 109-amino acid fragment of NS5B are typically associated with resistance to interferon (IFN) and ribavirin (RIB) and to new antiviral drugs. The prevalence of these mutations was examined in 69 drug-naïve individuals with hepatitis C virus (HCV) infections in Rio de Janeiro, Brazil. Mutations related to non-response to IFN/RIB were observed in all subtypes studied (1a, 1b, 2b, 3a and 4). The most common mutation was Q309R, present in all subtypes, except subtype 2b with frequency above 20%. D244N was detected only in subtype 3a and A333E was detected only in subtype 2b. We did not detect the S282T, S326G or T329I mutations in any of the samples analysed. Of note, the C316N mutation, previously related to a new non-nucleoside compound (HCV796 and AG-021541), was observed in only eight of 33 (24%) samples from subtype 1b. Site 316 was under positive selection in this HCV variant. Our data highlight the presence of previously described resistance mutations in HCV genotypes from drug-naïve patients.

Emergence of young human genes after a burst of retroposition in primates.

The origin of new genes through gene duplication is fundamental to the evolution of lineage- or species-specific phenotypic traits. In this report, we estimate the number of functional retrogenes on the lineage leading to humans generated by the high rate of retroposition (retroduplication) in primates. Extensive comparative sequencing and expression studies coupled with evolutionary analyses and simulations suggest that a significant proportion of recent retrocopies represent bona fide human genes. We estimate that at least one new retrogene per million years emerged on the human lineage during the past approximately 63 million years of primate evolution. Detailed analysis of a subset of the data shows that the majority of retrogenes are specifically expressed in testis, whereas their parental genes show broad expression patterns. Consistently, most retrogenes evolved functional roles in spermatogenesis. Proteins encoded by X chromosome-derived retrogenes were strongly preserved by purifying selection following the duplication event, supporting the view that they may act as functional autosomal substitutes during X-inactivation of late spermatogenesis genes. Also, some retrogenes acquired a new or more adapted function driven by positive selection. We conclude that retroduplication significantly contributed to the formation of recent human genes and that most new retrogenes were progressively recruited during primate evolution by natural and/or sexual selection to enhance male germline function.

Survival-independent function of NF-kappaB/Rel during late stages of thymocyte differentiation.

Transcription factors of the NF-kappaB/Rel family are important mediators of extracellular signals. Their implication in positive selection of thymocytes is suggested by a defective thymic development in transgenic mice that over-express IkappaB in thymocytes. These mice exhibit an accumulation of an unusually prominent population of TCRhigh/CD4/CD8 double positive cells in the thymus and a dramatic reduction of CD4+ and CD8+ cells in the periphery. The present study addresses the role of NF-kappaB in survival and differentiation processes of maturing thymocytes using IkappaB/bcl-2 and IkappaB/HY double-transgenic mice. Neither the introduction of the anti-apoptosis gene bcl-2 nor the positively selecting background in female HY transgenic mice resulted in a rescue of the maturational defects observed in the thymus of IkappaB transgenic mice. Thus, rather than promoting survival the main role of NF-kappaB/Rel proteins during positive selection of thymocytes appears to be the mediation of differentiation signals.

Evidence for an episodic model of protein sequence evolution.

The evolution of protein function appears to involve alternating periods of conservative evolution and of relatively rapid change. Evidence for such episodic evolution, consistent with some theoretical expectations, comes from the application of increasingly sophisticated models of evolution to large sequence datasets. We present here some of the recent methods to detect functional shifts, using amino acid or codon models. Both provide evidence for punctual shifts in patterns of amino acid conservation, including the fixation of key changes by positive selection. Although a link to gene duplication, a presumed source of functional changes, has been difficult to establish, this episodic model appears to apply to a wide variety of proteins and organisms.

Oseltamivir-resistant influenza A(H1N1)pdm2009 strains found in Brazil are endowed with permissive mutations, which compensate the loss of fitness imposed by antiviral resistance

The 2009 pandemic influenza A virus outbreak led to the systematic use of the neuraminidase (NA) inhibitor oseltamivir (OST). Consequently, OST-resistant strains, carrying the mutation H275Y, emerged in the years after the pandemics, with a prevalence of 1-2%. Currently, OST-resistant strains have been found in community settings, in untreated individuals. To spread in community settings, H275Y mutants must contain additional mutations, collectively called permissive mutations. We display the permissive mutations in NA of OST-resistant A(H1N1)pdm09 virus found in Brazilian community settings. The NAs from 2013 are phylogenetically distinct from those of 2012, indicating a tendency of positive selection of NAs with better fitness. Some previously predicted permissive mutations, such as V241I and N369K, found in different countries, were also detected in Brazil. Importantly, the change D344N, also predicted to compensate loss of fitness imposed by H275Y mutation, was found in Brazil, but not in other countries in 2013. Our results reinforce the notion that OST-resistant A(H1N1)pdm09 strains with compensatory mutations may arise in an independent fashion, with samples being identified in different states of Brazil and in different countries. Systematic circulation of these viral strains may jeopardise the use of the first line of anti-influenza drugs in the future.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Vitellogenin Underwent Subfunctionalization to Acquire Caste and Behavioral Specific Expression in the Harvester Ant Pogonomyrmex barbatus.

The reproductive ground plan hypothesis (RGPH) proposes that the physiological pathways regulating reproduction were co-opted to regulate worker division of labor. Support for this hypothesis in honeybees is provided by studies demonstrating that the reproductive potential of workers, assessed by the levels of vitellogenin (Vg), is linked to task performance. Interestingly, contrary to honeybees that have a single Vg ortholog and potentially fertile nurses, the genome of the harvester ant Pogonomyrmex barbatus harbors two Vg genes (Pb_Vg1 and Pb_Vg2) and nurses produce infertile trophic eggs. P. barbatus, thus, provides a unique model to investigate whether Vg duplication in ants was followed by subfunctionalization to acquire reproductive and non-reproductive functions and whether Vg reproductive function was co-opted to regulate behavior in sterile workers. To investigate these questions, we compared the expression patterns of P. barbatus Vg genes and analyzed the phylogenetic relationships and molecular evolution of Vg genes in ants. qRT-PCRs revealed that Pb_Vg1 is more highly expressed in queens compared to workers and in nurses compared to foragers. By contrast, the level of expression of Pb_Vg2 was higher in foragers than in nurses and queens. Phylogenetic analyses show that a first duplication of the ancestral Vg gene occurred after the divergence between the poneroid and formicoid clades and subsequent duplications occurred in the lineages leading to Solenopsis invicta, Linepithema humile and Acromyrmex echinatior. The initial duplication resulted in two Vg gene subfamilies preferentially expressed in queens and nurses (subfamily A) or in foraging workers (subfamily B). Finally, molecular evolution analyses show that the subfamily A experienced positive selection, while the subfamily B showed overall relaxation of purifying selection. Our results suggest that in P. barbatus the Vg gene underwent subfunctionalization after duplication to acquire caste- and behavior- specific expression associated with reproductive and non-reproductive functions, supporting the validity of the RGPH in ants.

Senescence in cell oxidative status in two bird species with contrasting life expectancy.

Oxidative stress occurs when the production of reactive oxygen species (ROS) by an organism exceeds its capacity to mitigate the damaging effects of the ROS. Consequently, oxidative stress hypotheses of ageing argue that a decline in fecundity and an increase in the likelihood of death with advancing age reported at the organism level are driven by gradual disruption of the oxidative balance at the cellular level. Here, we measured erythrocyte resistance to oxidative stress in the same individuals over several years in two free-living bird species with contrasting life expectancy, the great tit (known maximum life expectancy is 15.4 years) and the Alpine swift (26 years). In both species, we found evidence for senescence in cell resistance to oxidative stress, with patterns of senescence becoming apparent as subjects get older. In the Alpine swift, there was also evidence for positive selection on cell resistance to oxidative stress, the more resistant subjects being longer lived. The present findings of inter-individual selection and intra-individual deterioration in cell oxidative status at old age in free-living animals support a role for oxidative stress in the ageing of wild animals.

Birth and adaptive evolution of a hominoid gene that supports high neurotransmitter flux.

The enzyme glutamate dehydrogenase (GDH) is important for recycling the chief excitatory neurotransmitter, glutamate, during neurotransmission. Human GDH exists in housekeeping and brain-specific isotypes encoded by the genes GLUD1 and GLUD2, respectively. Here we show that GLUD2 originated by retroposition from GLUD1 in the hominoid ancestor less than 23 million years ago. The amino acid changes responsible for the unique brain-specific properties of the enzyme derived from GLUD2 occurred during a period of positive selection after the duplication event.

Evolutionary and Functional Analyses of the Interaction between the Myeloid Restriction Factor SAMHD1 and the Lentiviral Vpx Protein.

SAMHD1 has recently been identified as an HIV-1 restriction factor operating in myeloid cells. As a countermeasure, the Vpx accessory protein from HIV-2 and certain lineages of SIV have evolved to antagonize SAMHD1 by inducing its ubiquitin-proteasome-dependent degradation. Here, we show that SAMHD1 experienced strong positive selection episodes during primate evolution that occurred in the Catarrhini ancestral branch prior to the separation between hominoids (gibbons and great apes) and Old World monkeys. The identification of SAMHD1 residues under positive selection led to mapping the Vpx-interaction domain of SAMHD1 to its C-terminal region. Importantly, we found that while SAMHD1 restriction activity toward HIV-1 is evolutionarily maintained, antagonism of SAMHD1 by Vpx is species-specific. The distinct evolutionary signature of SAMHD1 sheds light on the development of its antiviral specificity.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

A generalized mechanistic codon model.

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.