Conservation of hepatitis C virus nonstructural protein 3 amino acid sequence in viral isolates during liver transplantation

As a consequence of selective pressure exerted by the immune response during hepatitis C virus (HCV) infection, a high rate of nucleotide mutations in the viral genome is observed which leads to the emergence of viral escape mutants. The aim of this study was to evaluate the evolution of the amino acid (aa) sequence of the HCV nonstructural protein 3 (NS3) in viral isolates after liver transplantation. Six patients with HCV-induced liver disease undergoing liver transplantation (LT) were followed up for sequence analysis. Hepatitis C recurrence was observed in all patients after LT. The rate of synonymous (dS) nucleotide substitutions was much higher than that of nonsynonymous (dN) ones in the NS3 encoding region. The high values of the dS/dN ratios suggest no sustained adaptive evolution selection pressure and, therefore, absence of specific NS3 viral populations. Clinical genotype assignments were supported by phylogenetic analysis. Serial samples from each patient showed lower mean nucleotide genetic distance when compared with samples of the same HCV genotype and subtype. The NS3 samples studied had an N-terminal aa sequence with several differences as compared with reference ones, mainly in genotype 1b-infected patients. After LT, as compared with the sequences before, a few reverted aa substitutions and several established aa substitutions were observed at the N-terminal of NS3. Sites described to be involved in important functions of NS3, notably those of the catalytic triad and zinc binding, remained unaltered in terms of aa sequence. Rare or frequent aa substitutions occurred indiscriminately in different positions. Several cytotoxic T lymphocyte epitopes described for HCV were present in our 1b samples. Nevertheless, the deduced secondary structure of the NS3 protease showed a few alterations in samples from genotype 3a patients, but none were seen in 1b cases. Our data, obtained from patients under important selective pressure during LT, show that the NS3 protease remains well conserved, mainly in HCV 3a patients. It reinforces its potential use as an antigenic candidate for further studies aiming at the development of a protective immune response.

Adaptation from standing genetic variation and from mutation

Dissertation presented to obtain the Ph.D degree in Evolutionary Biology

Evolutionary dynamics of Chlamydia trachomatis genome and identification of molecular patterns of hypothetical protein coding genes

The obligate intracellular bacterium Chlamydia trachomatis is a human pathogen of major public health significance. Strains can be classified into 15 main serovars (A to L3) that preferentially cause ocular infections (A-C), genital infections (D-K) or lymphogranuloma venereum (LGV) (L1-L3), but the molecular basis behind their distinct tropism, ecological success and pathogenicity is not welldefined. Most chlamydial research demands culture in eukaryotic cell lines, but it is not known if stains become laboratory adapted. By essentially using genomics and transcriptomics, we aimed to investigate the evolutionary patterns underlying the adaptation of C. trachomatis to the different human tissues, given emphasis to the identification of molecular patterns of genes encoding hypothetical proteins, and to understand the adaptive process behind the C. trachomatis in vivo to in vitro transition. Our results highlight a positive selection-driven evolution of C. trachomatis towards nichespecific adaptation, essentially targeting host-interacting proteins, namely effectors and inclusion membrane proteins, where some of them also displayed niche-specific expression patterns. We also identified potential "ocular-specific" pseudogenes, and pointed out the major gene targets of adaptive mutations associated with LGV infections. We further observed that the in vivo-derived genetic makeup of C. trachomatis is not significantly compromised by its long-term laboratory propagation. In opposition, its introduction in vitro has the potential to affect the phenotype, likely yielding virulence attenuation. In fact, we observed a "genital-specific" rampant inactivation of the virulence gene CT135, which may impact the interpretation of data derived from studies requiring culture. Globally, the findings presented in this Ph.D. thesis contribute for the understanding of C.trachomatis adaptive evolution and provides new insights into the biological role of C. trachomatishypothetical proteins. They also launch research questions for future functional studies aiming toclarify the determinants of tissue tropism, virulence or pathogenic dissimilarities among C. trachomatisstrains.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection.

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

Influence of learning on range expansion and adaptation to novel habitats.

Learning has been postulated to 'drive' evolution, but its influence on adaptive evolution in heterogeneous environments has not been formally examined. We used a spatially explicit individual-based model to study the effect of learning on the expansion and adaptation of a species to a novel habitat. Fitness was mediated by a behavioural trait (resource preference), which in turn was determined by both the genotype and learning. Our findings indicate that learning substantially increases the range of parameters under which the species expands and adapts to the novel habitat, particularly if the two habitats are separated by a sharp ecotone (rather than a gradient). However, for a broad range of parameters, learning reduces the degree of genetically-based local adaptation following the expansion and facilitates maintenance of genetic variation within local populations. Thus, in heterogeneous environments learning may facilitate evolutionary range expansions and maintenance of the potential of local populations to respond to subsequent environmental changes.

Evolutionary trajectories of primate genes involved in HIV pathogenesis.

The current availability of five complete genomes of different primate species allows the analysis of genetic divergence over the last 40 million years of evolution. We hypothesized that the interspecies differences observed in susceptibility to HIV-1 would be influenced by the long-range selective pressures on host genes associated with HIV-1 pathogenesis. We established a list of human genes (n = 140) proposed to be involved in HIV-1 biology and pathogenesis and a control set of 100 random genes. We retrieved the orthologous genes from the genome of humans and of four nonhuman primates (Pan troglodytes, Pongo pygmaeus abeli, Macaca mulatta, and Callithrix jacchus) and analyzed the nucleotide substitution patterns of this data set using codon-based maximum likelihood procedures. In addition, we evaluated whether the candidate genes have been targets of recent positive selection in humans by analyzing HapMap Phase 2 single-nucleotide polymorphisms genotyped in a region centered on each candidate gene. A total of 1,064 sequences were used for the analyses. Similar median K(A)/K(S) values were estimated for the set of genes involved in HIV-1 pathogenesis and for control genes, 0.19 and 0.15, respectively. However, genes of the innate immunity had median values of 0.37 (P value = 0.0001, compared with control genes), and genes of intrinsic cellular defense had K(A)/K(S) values around or greater than 1.0 (P value = 0.0002). Detailed assessment allowed the identification of residues under positive selection in 13 proteins: AKT1, APOBEC3G, APOBEC3H, CD4, DEFB1, GML, IL4, IL8RA, L-SIGN/CLEC4M, PTPRC/CD45, Tetherin/BST2, TLR7, and TRIM5alpha. A number of those residues are relevant for HIV-1 biology. The set of 140 genes involved in HIV-1 pathogenesis did not show a significant enrichment in signals of recent positive selection in humans (intraspecies selection). However, we identified within or near these genes 24 polymorphisms showing strong signatures of recent positive selection. Interestingly, the DEFB1 gene presented signatures of both interspecies positive selection in primates and intraspecies recent positive selection in humans. The systematic assessment of long-acting selective pressures on primate genomes is a useful tool to extend our understanding of genetic variation influencing contemporary susceptibility to HIV-1.

Odorant Receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura Lineages

Background: In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection. Methodology/Principal Findings: We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ~37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster ¿OR65c and OR67a ¿and D. pseudoobscura ¿OR65b1 and OR67c. Conclusions: The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.

Natural genetic variation in Arabidopsis: tools, traits and prospects for evolutionary ecology.

BACKGROUND: The model plant Arabidopsis thaliana (Arabidopsis) shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, the potential of Arabidopsis for molecular genetic analysis of this natural variation has increased dramatically in recent years. SCOPE: Advanced genomics has accelerated molecular phylogenetic analysis and gene identification by quantitative trait loci (QTL) mapping and/or association mapping in Arabidopsis. In particular, QTL mapping utilizing natural accessions is now becoming a major strategy of gene isolation, offering an alternative to artificial mutant lines. Furthermore, the genomic information is used by researchers to uncover the signature of natural selection acting on the genes that contribute to phenotypic variation. The evolutionary significance of such genes has been evaluated in traits such as disease resistance and flowering time. However, although molecular hallmarks of selection have been found for the genes in question, a corresponding ecological scenario of adaptive evolution has been difficult to prove. Ecological strategies, including reciprocal transplant experiments and competition experiments, and utilizing near-isogenic lines of alleles of interest will be a powerful tool to measure the relative fitness of phenotypic and/or allelic variants. CONCLUSIONS: As the plant model organism, Arabidopsis provides a wealth of molecular background information for evolutionary genetics. Because genetic diversity between and within Arabidopsis populations is much higher than anticipated, combining this background information with ecological approaches might well establish Arabidopsis as a model organism for plant evolutionary ecology.

The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.

A Novel Isoform of the Hepatic Antimicrobial Peptide, Hepcidin (Hepc-CB1), from a Deep-Sea Fish, the Spinyjaw Greeneye Chlorophthalmus bicornis (Norman, 1939): Molecular Characterisation and Phylogeny

Hepcidin is cysteine-rich short peptide of innate immune system of fishes, equipped to perform prevention and proliferation of invading pathogens like bacteria and viruses by limiting iron availability and activating intracellular cascades. Hepcidins are diverse in teleost fishes, due to the varied aquatic environments including exposure to pathogens, oxygenation and iron concentration. In the present study, we report a 87-amino acid (aa) preprohepcidin (Hepc-CB1) with a signal peptide of 24 aa, a prodomain of 39 aa and a bioactive mature peptide of 24 aa from the gill mRNA transcripts of the deep-sea fish spinyjaw greeneye, Chlorophthalmus bicornis. Molecular characterisation and phylogenetic analysis categorised the peptide to HAMP2-like group with a mature peptide of 2.53 kDa; a net positive charge (?3) and capacity to form b-hairpin-like structure configured by 8 conserved cysteines. The present work provides new insight into the mass gene duplication events and adaptive evolution of hepcidin isoforms with respect to environmental influences and positive Darwinian selection. This work reports a novel hepcidin isoform under the group HAMP2 from a nonacanthopterygian deep-sea fish, C. bicornis

Determination of mitochondrial genetic diversity in mammals

Mitochondrial DNA (mtDNA) is one of the most Popular population genetic markers. Its relevance as an indicator Of Population size and history has recently been questioned by several large-scale studies in animals reporting evidence for recurrent adaptive evolution, at least in invertebrates. Here we focus on mammals, a more restricted taxonomic group for which the issue of mtDNA near neutrality is crucial. By analyzing the distribution of mtDNA diversity across species and relating 4 to allozyme diversity, life-history traits, and taxonomy, we show that (i) mtDNA in mammals (toes not reject the nearly neutral model; (ii) mtDNA diversity, however, is unrelated to any of the 14 life-history and ecological variables that we analyzed, including body mass, geographic range, and The World Conservation Union (IUCN) categorization; (iii) mtDNA diversity is highly variable between mammalian orders and families; (iv) this taxonomic effect is most likely explained by variations of mutation rate between lineages. These results are indicative of a strong stochasticity of effective population size in mammalian species. They Suggest that, even in the absence of selection, mtDNA genetic diversity is essentially unpredictable, knowing species biology, and probably uncorrelated to species abundance.

EVOLUTIONARY BIOLOGY IN BIODIVERSITY SCIENCE, CONSERVATION, AND POLICY: A CALL TO ACTION

Evolutionary biologists have long endeavored to document how many species exist on Earth, to understand the processes by which biodiversity waxes and wanes, to document and interpret spatial patterns of biodiversity, and to infer evolutionary relationships. Despite the great potential of this knowledge to improve biodiversity science, conservation, and policy, evolutionary biologists have generally devoted limited attention to these broader implications. Likewise, many workers in biodiversity science have underappreciated the fundamental relevance of evolutionary biology. The aim of this article is to summarize and illustrate some ways in which evolutionary biology is directly relevant We do so in the context of four broad areas: (1) discovering and documenting biodiversity, (2) understanding the causes of diversification, (3) evaluating evolutionary responses to human disturbances, and (4) implications for ecological communities, ecosystems, and humans We also introduce bioGENESIS, a new project within DIVERSITAS launched to explore the potential practical contributions of evolutionary biology In addition to fostering the integration of evolutionary thinking into biodiversity science, bioGENESIS provides practical recommendations to policy makers for incorporating evolutionary perspectives into biodiversity agendas and conservation. We solicit your involvement in developing innovative ways of using evolutionary biology to better comprehend and stem the loss of biodiversity.

Role of the triad N46, S106 and T107 and the surface charges in the determination of the acidic pH optimum of digestive lysozymes from Musca domestica

Structures of digestive lysozymes 1 and 2 from housefly (MdL1 and MdL2) show that S106-T107 delimit a polar pocket around E32 (catalytic acid/base) and N46 contributes to the positioning of 050 (catalytic nucleophile), whereas those residues are replaced by V109-A110 and D48 in the non-digestive lysozyme from hen egg-white (HEWL). Further analyses revealed that MdL1 and MdL2 surfaces are less positively charged than HEWL surface. To verify the relevance of these differences to the acidic pH optimum of digestive lysozymes it was determined that pKas of the catalytic residues of the triple mutant MdL2 (N46D-S106V-T107A) are similar to HEWL pKas and higher than those for MdL2. In agreement, triple mutant MdL2 and HEWL exhibits the same pH optimum upon methylumbelliferylchitotrioside. In addition to that, the introduction of six basic residues on MdL1 surface increased by 1 unit the pH optimum for the activity upon bacterial walls. Thus, the acidic pH optimum for MdL2 and MdL1 activities upon methylumbelliferylchitotrioside is determined by the presence of N46, S106 and T107 in the environment of their catalytic residues, which favors pKas reduction. Conversely, acidic pH optimum upon bacterial walls is determined by a low concentration of positive charges on the MdL2 and MdL1 surfaces. (C) 2010 Elsevier Inc. All rights reserved.

Thermal relationships and exercise physiology in anuran amphibians: Integration and evolutionary implications

Thermal and water balance are coupled in anurans, and species with particularly permeable skin avoid overheating more effectively than minimizing variance of body temperature. In turn, temperature affects muscle performance in several ways, so documenting the mean and variance of body temperature of active frogs can help explain variation in behavioral performance. The two types of activities studied in most detail, jumping and calling, differ markedly in duration and intensity, and there are distinct differences in the metabolic profile and fiber type of the supporting muscles. Characteristics of jumping and calling also vary significantly among species, and these differences have a number of implications that we discuss in some detail throughout this paper. One question that emerges from this topic is whether anuran species exhibit activity temperatures that match the temperature range over which they perform best. Although this seems the case, thermal preferences are variable and may not necessarily reflect typical activity temperatures. The performance versus temperature curves and the thermal limits for anuran activity reflect the thermal ecology of species more than their systematic position. Anuran thermal physiology, therefore, seems to be phenotypically plastic and susceptible to adaptive evolution. Although generalizations regarding the mechanistic basis of such adjustments are not yet possible, recent attempts have been made to reveal the mechanistic basis of acclimation and acclimatization. (C) 2007 Elsevier B.V. All rights reserved.