The four hexamerin genes in the honey bee: structure, molecular evolution and function deduced from expression patterns in queens, workers and drones

Background: Hexamerins are hemocyanin-derived proteins that have lost the ability to bind copper ions and transport oxygen; instead, they became storage proteins. The current study aimed to broaden our knowledge on the hexamerin genes found in the honey bee genome by exploring their structural characteristics, expression profiles, evolution, and functions in the life cycle of workers, drones and queens. Results: The hexamerin genes of the honey bee (hex 70a, hex 70b, hex 70c and hex 110) diverge considerably in structure, so that the overall amino acid identity shared among their deduced protein subunits varies from 30 to 42%. Bioinformatics search for motifs in the respective upstream control regions (UCRs) revealed six overrepresented motifs including a potential binding site for Ultraspiracle (Usp), a target of juvenile hormone (JH). The expression of these genes was induced by topical application of JH on worker larvae. The four genes are highly transcribed by the larval fat body, although with significant differences in transcript levels, but only hex 110 and hex 70a are re-induced in the adult fat body in a caste-and sex-specific fashion, workers showing the highest expression. Transcripts for hex 110, hex 70a and hex70b were detected in developing ovaries and testes, and hex 110 was highly transcribed in the ovaries of egg-laying queens. A phylogenetic analysis revealed that HEX 110 is located at the most basal position among the holometabola hexamerins, and like HEX 70a and HEX 70c, it shares potential orthology relationship with hexamerins from other hymenopteran species. Conclusions: Striking differences were found in the structure and developmental expression of the four hexamerin genes in the honey bee. The presence of a potential binding site for Usp in the respective 5' UCRs, and the results of experiments on JH level manipulation in vivo support the hypothesis of regulation by JH. Transcript levels and patterns in the fat body and gonads suggest that, in addition to their primary role in supplying amino acids for metamorphosis, hexamerins serve as storage proteins for gonad development, egg production, and to support foraging activity. A phylogenetic analysis including the four deduced hexamerins and related proteins revealed a complex pattern of evolution, with independent radiation in insect orders.

Molecular evolution and phylogeny of elapid snake venom three-finger toxins

Animal venom components are of considerable interest to researchers across a wide variety of disciplines, including molecular biology, biochemistry, medicine, and evolutionary genetics. The three-finger family of snake venom peptides is a particularly interesting and biochemically complex group of venom peptides, because they are encoded by a large multigene family and display a diverse array of functional activities. In addition, understanding how this complex and highly varied multigene family evolved is an interesting question to researchers investigating the biochemical diversity of these peptides and their impact on human health. Therefore, the purpose of our study was to investigate the long-term evolutionary patterns exhibited by these snake venom toxins to understand the mechanisms by which they diversified into a large, biochemically diverse, multigene family. Our results show a much greater diversity of family members than was previously known, including a number of subfamilies that did not fall within any previously identified groups with characterized activities. In addition, we found that the long-term evolutionary processes that gave rise to the diversity of three-finger toxins are consistent with the birth-and-death model of multigene family evolution. It is anticipated that this three-finger toxin toolkit will prove to be useful in providing a clearer picture of the diversity of investigational ligands or potential therapeutics available within this important family.

High rates of molecular evolution in hantaviruses

Hantaviruses are rodent-borne Bunyaviruses that infect the Arvicolinae, Murinae, and Sigmodontinae subfamilies of Muridae. The rate of molecular evolution in the hantaviruses has been previously estimated at approximately 10(-7) nucleotide substitutions per site, per year (substitutions/site/year), based on the assumption of codivergence and hence shared divergence times with their rodent hosts. If substantiated, this would make the hantaviruses among the slowest evolving of all RNA viruses. However, as hantaviruses replicate with an RNA-dependent RNA polymerase, with error rates in the region of one mutation per genome replication, this low rate of nucleotide substitution is anomalous. Here, we use a Bayesian coalescent approach to estimate the rate of nucleotide substitution from serially sampled gene sequence data for hantaviruses known to infect each of the 3 rodent subfamilies: Araraquara virus ( Sigmodontinae), Dobrava virus ( Murinae), Puumala virus ( Arvicolinae), and Tula virus ( Arvicolinae). Our results reveal that hantaviruses exhibit shortterm substitution rates of 10(-2) to 10(-4) substitutions/site/year and so are within the range exhibited by other RNA viruses. The disparity between this substitution rate and that estimated assuming rodent-hantavirus codivergence suggests that the codivergence hypothesis may need to be reevaluated.

Molecular clocks in reptiles: Life history influences rate of molecular evolution

Life history has been implicated as a determinant of variation in rate of molecular evolution amongst vertebrate species because of a negative correlation between bode size and substitution rate for many Molecular data sets. Both the generality and the cause of the negative bode size trend have been debated, and the validity of key studies has been questioned (particularly concerning the failure to account for phylogenetic bias). In this study, a comparative method has been used to test for an association between a range of life-history variables-such as body size age at maturity, and clutch size-and DNA substitution rate for three genes (NADH4, cytochrome b, and c-mos). A negative relationship between body size and rate of molecular evolution was found for phylogenetically independent pairs of reptile species spanning turtles. lizards. snakes, crocodile, and tuatara. Although this Study was limited by the number of comparisons for which both sequence and lite-history data were available, the results, suggest that a negative bode size trend in rate of molecular evloution may be a general feature of reptile molecular evolution. consistent with similar studies of mammals and birds. This observation has important implications for uncovering the mechanisms of molecular evolution and warns against assuming that related lineages will share the same substitution rate (a local molecular clock) in order to date evolutionary divergences from DNA sequences.

Testing the link between the latitudinal gradient in species richness and rates of molecular evolution

Numerous hypotheses have been proposed to explain latitudinal gradients in species richness, but all are subject to ongoing debate. Here we examine Rohde's (1978, 1992) hypothesis, which proposes that climatic conditions at low latitudes lead to elevated rates of speciation. This hypothesis predicts that rates of molecular evolution should increase towards lower latitudes, but this prediction has never been tested. We discuss potential links between rates of molecular evolution and latitudinal diversity gradients, and present the first test of latitudinal variation in rates of molecular evolution. Using 45 phylogenetically independent, latitudinally separated pairs of bird species and higher taxa, we compare rates of evolution of two mitochondrial genes and DNA-DNA hybridization distances. We find no support for an effect of latitude on rate of molecular evolution. This result casts doubt on the generality of a key component of Rohde's hypothesis linking climate and speciation.

The Role of the Humoral Immune Response in the Molecular Evolution of the Envelope C2, V3 and C3 Regions in Chronically HIV-2 Infected Patients

BACKGROUND: This study was designed to investigate, for the first time, the short-term molecular evolution of the HIV-2 C2, V3 and C3 envelope regions and its association with the immune response. Clonal sequences of the env C2V3C3 region were obtained from a cohort of eighteen HIV-2 chronically infected patients followed prospectively during 2-4 years. Genetic diversity, divergence, positive selection and glycosylation in the C2V3C3 region were analysed as a function of the number of CD4+ T cells and the anti-C2V3C3 IgG and IgA antibody reactivity RESULTS: The mean intra-host nucleotide diversity was 2.1% (SD, 1.1%), increasing along the course of infection in most patients. Diversity at the amino acid level was significantly lower for the V3 region and higher for the C2 region. The average divergence rate was 0.014 substitutions/site/year, which is similar to that reported in chronic HIV-1 infection. The number and position of positively selected sites was highly variable, except for codons 267 and 270 in C2 that were under strong and persistent positive selection in most patients. N-glycosylation sites located in C2 and V3 were conserved in all patients along the course of infection. Intra-host variation of C2V3C3-specific IgG response over time was inversely associated with the variation in nucleotide and amino acid diversity of the C2V3C3 region. Variation of the C2V3C3-specific IgA response was inversely associated with variation in the number of N-glycosylation sites. CONCLUSION: The evolutionary dynamics of HIV-2 envelope during chronic aviremic infection is similar to HIV-1 implying that the virus should be actively replicating in cellular compartments. Convergent evolution of N-glycosylation in C2 and V3, and the limited diversification of V3, indicates that there are important functional constraints to the potential diversity of the HIV-2 envelope. C2V3C3-specific IgG antibodies are effective at reducing viral population size limiting the number of virus escape mutants. The C3 region seems to be a target for IgA antibodies and increasing N-linked glycosylation may prevent HIV-2 envelope recognition by these antibodies. Our results provide new insights into the biology of HIV-2 and its relation with the human host and may have important implications for vaccine design.

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.

Molecular evolution of the Drosophila olfactory system

The olfactory system is an attractive model to study the genetic mechanisms underlying evolution of the nervous system. This sensory system mediates the detection and behavioural responses to an enormous diversity of volatile chemicals in the environment and displays rapid evolution, as species acquire, modify and discard olfactory receptors and circuits to adapt to new olfactory stimuli. Drosophilids provide an attractive model to study these processes. The availability of 12 sequenced genomes of Drosophila species occupying diverse ecological niches provides a rich resource for genomic analyses. Moreover, one of these species, Drosophila melanogaster, is amenable to a powerful combination of genetic and electrophysiological analyses. D. melanogaster has two distinct families of olfactory receptors to detect odours, the well-characterised Odorant Receptors (ORs) and the recently identified lonotropic Receptors (IRs). In my thesis, I have provided new insights into the genetic mechanisms underlying olfactory system evolution through three distinct, but interrelated projects. First, I performed a comparative genomic analysis of the IR repertoire in 12 sequenced Drosophila species, which has revealed that the olfactory IRs are highly conserved across species. By contrast, a large fraction of IRs that are not expressed in the olfactory system - and which may be gustatory receptors - are much more variable in sequence and gene copy number. Second, to identify ligands for IR expressing olfactory sensory neurons, I have performed an electrophysiological screen in D. melanogaster using a panel of over 160 odours. I found that the IRs respond to a number of amines, aldehydes and acids, contrasting with the chemical specificity of the OR repertoire, which is mainly tuned to esters, alcohols and ketones. Finally, the identification of ligands for IRs in this species allowed me to investigate in detail the molecular and functional evolution of a tandem array of IRs, IR75a/IR75b/IR75c, in D. sechellia. This species is endemic to the Seychelles archipelago and highly specialised to breed on the fruits of Morinda citrifolia, which is repulsive and toxic for other Drosophila species. These studies led me to discover that receptor loss, changes in receptor specificity and changes in receptor expression have likely played an important role during the evolution of these IRs in D. sechellia. These changes may explain, in part, the unique chemical ecology of this species. - Le système olfactif est un excellent modèle pour étudier les mécanismes génétiques impliqués dans l'étude de l'évolution du système nerveux. Ce système sensoriel permet la détection de nombreux composés volatils présents dans l'environnement et est à la base des réponses comportementales. Il est propre à chaque espèce et évolue rapidement en modifiant ou en éliminant des récepteurs et leurs circuits olfactifs correspondants pour s'adapter à de nouvelles odeurs. Pour étudier le système olfactif et son évolution, nous avons décidé d'utiliser la drosophile comme modèle. Le séquençage complet de 12 souches de drosophiles habitant différentes niches écologiques permet une analyse génomique conséquente. De plus, l'une de ces espèces Drosophila melanogaster permet la combinaison d'analyses génétiques et électrophysiologiques. En effet, D. melanogaster possède 2 familles distinctes de récepteurs olfactifs qui permettent la détection d'odeurs: les récepteurs olfactifs (ORs) étant les mieux caractérisés et les récepteurs ionotropiques (IRs), plus récemment identifiés. Au cours de ma thèse, j'ai apporté des nouvelles connaissances qui m'ont permis de mieux comprendre les mécanismes génétiques à la base de l'évolution du système olfactif au travers de trois projets différents, mais interdépendants. Premièrement, j'ai réalisé une analyse génomique comparative de l'ensemble des IRs dans les 12 souches de drosophiles séquencées jusqu'à présent. Ceci a montré que les récepteurs olfactifs IRs sont hautement conservés parmi l'ensemble de ces espèces. Au contraire, une grande partie des IRs qui ne sont pas exprimés dans le système olfactif, et qui semblent être des récepteurs gustatifs, sont beaucoup plus variables dans leur séquence et dans le nombre de copie de gènes. Deuxièmement, pour identifier les ligands des récepteurs IRs exprimés par les neurones sensoriels olfactifs, j'ai réalisé une étude électrophysiologique chez D. melanogaster e η testant l'effet de plus de 160 composés chimiques sur les IRs. J'ai trouvé que les IRs répondent à un nombre d'amines, d'aldéhydes et d'acides, contrairement aux récepteurs olfactifs ORs qui eux répondent principalement aux esthers, alcools et cétones. Finalement, l'identification de ligands pour les IRs dans ces espèces m'a permis d'étudier en détail l'évolution fonctionnelle et moléculaire des IR75a/IR75b/IR75c dans D. sechellia. Cette espèce est endémique de l'archipel des Seychelles et se nourrit spécifiquement du fruit Morinda citrifolia qui est répulsif et toxique pour d'autres souches de drosophiles. Ces études m'ont poussé à découvrir que, la perte de IR75a, le changement dans la spécificité de IR75b ainsi que le changement dans l'expression de IR75c ont probablement joué un rôle important dans l'évolution des IRs chez D. sechellia. Ces changements peuvent expliquer, en partie, l'écologie chimique propre à cette espèce. Résumé français large public Le système olfactif permet aux animaux de détecter des milliers de molécules odorantes, les aidant ainsi à trouver de la nourriture, à distinguer si elle est fraîche ou avariée, à trouver des partenaires sexuels, ainsi qu'à éviter les prédateurs. Selon l'environnement et le mode de vie des espèces, le système olfactif doit détecter des odeurs très diverses ; en effet, un moustique qui recherche du sang humain pour se nourrir doit détecter des odeurs bien différentes d'une abeille qui recherche des fleurs. Dans ma thèse, j'ai essayé de comprendre comment les systèmes olfactifs d'une espèce évoluent pour s'adapter aux exigences induites par son environnement. Un très bon modèle pour étudier cela est la drosophile dont les différentes espèces se nichent dans des habitats très divers. Pour ce faire, j'ai étudié les récepteurs olfactifs de différentes espèces de la drosophile. Ces récepteurs sont des protéines qui se lient à des odeurs spécifiques. Lorsqu'ils se lient, ils activent un neurone qui envoie un signal électrique au cerveau. Ce signal est ensuite traité par ce dernier qui indique à la mouche si l'odeur est attractive ou répulsive. J'ai identifié les récepteurs olfactifs de plusieurs espèces de drosophile et étudié s'il y avait des différences entre elles. La plupart des récepteurs sont similaires entre les espèces, cependant dans l'une d'entre elles, certains récepteurs sont différents. Ce fait est particulièrement intéressant car cette espèce de drosophile se nourrit de fruits que les autres espèces n'apprécient pas. Comme nous ne savons pas quels récepteurs se lient à quelles odeurs, j'ai testé un grand nombre de composants odorants. Ceci m'a permis de constater que, effectivement, certains changements produits dans ces récepteurs expliquent pourquoi cette espèce aime particulièrement ces fruits. En outre, mes résultats contribuent à mieux comprendre les changements génétiques qui sont impliqués dans l'évolution du système olfactif.

Genome-wide analysis of adaptive molecular evolution in the carnivorous plant Utricularia gibba

The genome of the bladderwort Utricularia gibba provides an unparalleled opportunity to uncover the adaptive landscape of an aquatic carnivorous plant with unique phenotypic features such as absence of roots, development of water-filled suction bladders, and a highly ramified branching pattern. Despite its tiny size, the U. gibba genome accommodates approximately as many genes as other plant genomes. To examine the relationship between the compactness of its genome and gene turnover, we compared the U. gibba genome with that of four other eudicot species, defining a total of 17,324 gene families (orthogroups). These families were further classified as either 1) lineage-specific expanded/contracted or 2) stable in size. The U. gibba-expanded families are generically related to three main phenotypic features: 1) trap physiology, 2) key plant morphogenetic/developmental pathways, and 3) response to environmental stimuli, including adaptations to life in aquatic environments. Further scans for signatures of protein functional specialization permitted identification of seven candidate genes with amino acid changes putatively fixed by positive Darwinian selection in the U. gibba lineage. The Arabidopsis orthologs of these genes (AXR, UMAMIT41, IGS, TAR2, SOL1, DEG9, and DEG10) are involved in diverse plant biological functions potentially relevant for U. gibba phenotypic diversification, including 1) auxin metabolism and signal transduction, 2) flowering induction and floral meristem transition, 3) root development, and 4) peptidases. Taken together, our results suggest numerous candidate genes and gene families as interesting targets for further experimental confirmation of their functional and adaptive roles in the U. gibba's unique lifestyle and highly specialized body plan.

Molecular Evolution of Metazoan Hypoxia-inducing Factors

Most metazoans rely on aerobic energy production, which is dependent on adequate oxygen supply. In the case of reduced oxygen supply (hypoxia), the most profound changes in gene expression are mediated by transcription factors named hypoxia-inducible factors (HIF alpha). These proteins are post-translationally regulated by prolyl-4-hydroxylase (PHD) enzymes that are direct “sensors” of cellular oxygen levels. This thesis examines the molecular evolution of metazoan HIF systems. In early metazoans the HIF system emerged from pre-existing PHD oxygen sensors and early bHLH-PAS transcription factors. In invertebrates our analysis revealed an unexpected diversity of PHD genes and HIF alpha sequence characteristics. An early branching vertebrate, the epaulette shark (Hemiscyllium ocellatum) was chosen for sequencing and hypoxia preconditioning studies of HIF alpha and PHD genes. As no quantitative PCR reference genes were available, this thesis includes the first study of reference genes in cartilaginous fish species. Applying multiple statistical analysis we also discoveredthat commonly used reference gene software may perform poorly with some data sets. Novel reference genes allowed accurate measurements of the mRNAlevels of the studied target genes. Cartilaginous fishes have three genomic duplicates of both HIF alpha and PHD genes like mammals and teleost fishes. Combining functional divergence and selection analyses it was possible to describe how sequence changes in both HIF alpha and PHD duplicates may have contributed to the differential oxygen sensitivityof HIF alphas. Additionally, novel teleost HIF-1 alpha sequences were produced and used to reveal the molecular evolution of HIF-1 alpha in this lineage rich with hypoxia tolerant species.

Bayesian coalescent analysis reveals a high rate of molecular evolution in GB virus C

GB virus C/hepatitis G (GBV-C) is an RNA virus of the family Flaviviridae. Despite replicating with an RNA-dependent RNA polymerase, some previous estimates of rates of evolutionary change in GBV-C suggest that it fixes mutations at the anomalously low rate of similar to 100(-7) nucleotide substitution per site, per year. However, these estimates were largely based on the assumption that GBV-C and its close relative GBV-A (New World monkey GB viruses) codiverged with their primate hosts over millions of years. Herein, we estimated the substitution rate of GBV-C using the largest set of dated GBV-C isolates compiled to date and a Bayesian coalescent approach that utilizes the year of sampling and so is independent of the assumption of codivergence. This revealed a rate of evolutionary change approximately four orders of magnitude higher than that estimated previously, in the range of 10(-2) to 10(-3) sub/site/year, and hence in line with those previously determined for RNA viruses in general and the Flaviviridae in particular. In addition, we tested the assumption of host-virus codivergence in GBV-A by performing a reconciliation analysis of host and virus phylogenies. Strikingly, we found no statistical evidence for host-virus codivergence in GBV-A, indicating that substitution rates in the GB viruses should not be estimated from host divergence times.

Hepatitis C virus molecular evolution: Transmission, disease progression and antiviral therapy

Hepatitis C virus (HCV) infection represents an important public health problem worldwide. Reduction of HCV morbidity and mortality is a current challenge owned to several viral and host factors. Virus molecular evolution plays an important role in HCV transmission, disease progression and therapy outcome. The high degree of genetic heterogeneity characteristic of HCV is a key element for the rapid adaptation of the intrahost viral population to different selection pressures (e.g., host immune responses and antiviral therapy). HCV molecular evolution is shaped by different mechanisms including a high mutation rate, genetic bottlenecks, genetic drift, recombination, temporal variations and compartmentalization. These evolutionary processes constantly rearrange the composition of the HCV intrahost population in a staging manner. Remarkable advances in the understanding of the molecular mechanism controlling HCV replication have facilitated the development of a plethora of direct-acting antiviral agents against HCV. As a result, superior sustained viral responses have been attained. The rapidly evolving field of anti-HCV therapy is expected to broad its landscape even further with newer, more potent antivirals, bringing us one step closer to the interferon-free era. (C) 2014 Baishideng Publishing Group Inc. All rights reserved.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.