Genomic organization, nucleotide sequence analysis of the core histone genes cluster in Chlamys farreri and molecular evolution assessment of the H2A and H2B

This work represents the nucleotide sequence of the core histone gene cluster from scallop Chlamys farreri. The tandemly repeated unit of 5671 bp containing a copy of the four core histone genes H4, H2B, H2A and H3 was amplified and identified by the techniques of homology cloning and genomic DNA walking. All the histone genes in the cluster had the structures in their 3' flanking region which related to the evolution of histone gene expression patterns throughout the cell cycle, including two different termination signals, the hairpin structure and at least one AATAAA polyadenylation signal. In their 5' region, the transcription initiation sites with a conserved sequence of 5'-PyATTCPu-3' known as the CAP site were present in all genes except to H2B, generally 37-45 bp upstream of the start code. Canonical TATA and CAAT boxes were identified only in certain histone genes. In the case of the promoters of H2B and H2A genes, there was a 5'-GATCC-3' element, which had been found to be essential to start transcription at the appropriate site. After this element, in the promoter of H2B, there was another sequence, 5'-GGATCGAAACGTTC-3', which was similar to the consensus sequence of 5'-GGAATAAACGTATTC-3' corresponding to the H2B-specific promoter element. The presence of enhancer sequences (5'-TGATATATG-3') was identified from the H4 and H3 genes, matching perfectly with the consensus sequence defined for histone genes. There were several slightly more complex repetitive DNA in the intergene regions. The presence of the series of conserved sequences and reiterated sequences was consistent with the view that mollusc histone gene cluster arose by duplicating of an ancestral precursor histone gene, the birth-and-death evolution model with strong purifying selection enabled the histone cluster less variation and more conserved function. Meanwhile, the H2A and the H2B were demonstrated to be potential good marks for phylogenetic analysis. All the results will be contributed to the characterization of repeating histone gene families in molluscs.

Molecular evolution of proadrenomedullin N-terminal 20 peptide (PAMP):evidence for gene co-option

Posttranslational processing of proadrenomedullin generates two biologically active peptides, adrenomedullin (AM) and proadrenomedullin N-terminal 20 peptide (PAMP). Sequence comparison of homologous proadrenomedullin genes in vertebrate evolution shows a high degree of stability in the reading frame for AM, whereas PAMP sequence changes rapidly. Here we investigate the functional significance of PAMP phylogenetic variation studying two of PAMP's better characterized physiological activities, angiogenic potential and antimicrobial capability, with synthetic peptides carrying the predicted sequence for human, mouse, chicken, and fish PAMP. All tested peptides induced angiogenesis when compared with untreated controls, but chicken and fish PAMP, which lack terminal amidation, were apparently less angiogenic than their human and mouse homologs. Confirming the role of amidation in angiogenesis, Gly-extended and free acid variants of human PAMP produced responses similar to the natural nonamidated peptides. In contrast, antimicrobial activity was restricted to human PAMP, indicating that this function may have been acquired at a late time during the evolution of PAMP. Interestingly, free acid human PAMP retained antimicrobial activity whereas the Gly-extended form did not. This fact may reflect the need for maintaining a tightly defined structural conformation in the pore-forming mechanism proposed for these antimicrobial agents. The evolution of PAMP provides an example of an angiogenic peptide that developed antimicrobial capabilities without losing its original function.

Biodiversity and molecular evolution of malacostraca

No actual cenário de perda acelerada de biodiversidade, o nosso conhecimento dos ecossistemas marinhos, apesar da sua extensão e complexidade, continua muito inferior ao dos ecossistemas terrestres. A classe Malacostraca (Arthropoda, Crustacea), um grupo dos mais representativos nos ecossistemas marinhos, apresenta um elevado nível de diversidade morfológica e ecológica, mas difícil sua identificação ao nível de espécie requer frequentemente a ajuda de especialistas em taxonomia. A utilização recente do “barcoding” (código de barras do ADN), revelou ser um método rápido e eficaz para a identificação de espécies em diversos grupos de metazoários, incluindo os Malacostraca. No âmbito desta tese foi construída uma base de dados de código de barras de ADN envolvendo 132 espécies de Malacostraca vários locais de amostragem no Atlântico Nordeste e Mediterrâneo. As sequências de ADN mitocondrial provenientes de 601 espécimes formaram, em 95% dos casos, grupos congruentes com as identificações baseadas em características morfológicas. No entanto, foi detectado polimorfismo em seis casos e a divergência intra-específica foi elevada em exemplares pertencentes a duas espécies morfológicas, sugerindo, neste caso, a ocorrência de especiação críptica. Este estudo confirma a utilidade do código de barras de ADN para a identificação de Malacostraca marinhos. Apesar do sucesso obtido, este método apresenta alguns problemas, como por exemplo a possível amplificação de pseudogenes. A ocorrência de pseudogenes e as possíveisabordagens para a detecção e resolução deste tipo de problemas são discutidas com base em casos de estudo: análises dos códigos de barras ADN na espécie Goneplax rhomboides (Crustacea, Decapoda). A análise dos códigos de barras ADN revelou ainda grupos prioritários de decápodes para estudos taxonómicos e sistemáticos, nomeadamente os decápodes dos géneros Plesionika e Pagurus. Neste âmbito são discutidas as relações filogenéticas entre espécies seleccionadas dos géneros Plesionika e Pagurus. Este trabalho aponta para várias questões no âmbito da biodiversidade e evolução molecular da classe Malacostraca que carecem de um maior esclarecimento, podendo ser considerado como a base para estudo futuros. Análises filogenéticas adicionais integrando dados morfológicos e moleculares de um maior número de espécies e de famílias deverão certamente conduzir a uma melhor avaliação da biodiversidade e da evolução dentro da classe.

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.

Design, realization and characterization of automated millifluidic bioreactors for investigating the molecular evolution of lytic or lysogenic vector phages infecting engineered host e. Coli

In questa tesi viene presentato un bioreattore in grado di mantenere nel tempo condizioni biologiche tali che consentano di massimizzare i cicli di evoluzione molecolare di vettori di clonazione fagici: litico (T7) o lisogeno (M13). Verranno quindi introdtti concetti legati alla Teoria della Quasispecie e alla relazione tra errori di autoreplicazione e pressioni selettive naturali o artificiali su popolazioni di virus: il modello naturale del sistema evolutivo. Tuttavia, mantenere delle popolazioni di virus significa formire loro un substrato dove replicare. Per fare ciò, altri gruppi di ricerca hanno giá sviluppato complessi e costosi prototipi di macchinari per la crescita continua di popolazioni batteriche: i compartimenti dei sistemi evolutivi. Il bioreattore, oggetto di questo lavoro, fa parte del progetto europeo Evoprog: general purpose programmable machine evolution on a chip (Jaramillo’s Lab, University of Warwick) che, utilizzando tecnologie fagiche e regolazioni sintetiche esistenti, sará in grado di produrre funzionalità biocomputazionali di due ordini di grandezza più veloci rispetto alle tecniche convenzionali, riducendo allo stesso tempo i costi complessivi. Il primo prototipo consiste in uno o piú fermentatori, dove viene fatta crescere la cultura batterica in condizioni ottimizzate di coltivazione continua, e in un cellstat, un volume separato, dove avviene solo la replicazione dei virus. Entrambi i volumi sono di pochi millilitri e appropriatamente interconnessi per consentire una sorta di screening continuo delle biomolecole prodotte all’uscita. Nella parte finale verranno presentati i risultati degli esperimenti preliminari, a dimostrazione dell’affidabilità del prototipo costruito e dei protocolli seguiti per la sterilizzazione e l’assemblaggio del bioreattore. Gli esperimenti effettuati dimostrano il successo di due coltivazioni virali continue e una ricombinazione in vivo di batteriofagi litici o lisogeni ingegnerizzati. La tesi si conclude valutando i futuri sviluppi e i limiti del sistema, tenendo in considerazione, in particolare, alcune applicazioni rivolte agli studi di una terapia batteriofagica.

Adaptive Molecular Evolution in the Opsin Genes of Rapidly Speciating Cichlid Species

Cichlid fish inhabit a diverse range of environments that vary in the spectral content of light available for vision. These differences should result in adaptive selective pressure on the genes involved in visual sensitivity, the opsin genes. This study examines the evidence for differential adaptive molecular evolution in East African cichlid opsin genes due to gross differences in environmental light conditions. First, we characterize the selective regime experienced by cichlid opsin genes using a likelihood ratio test format, comparing likelihood models with different constraints on the relative rates of amino acid substitution, across sites. Second, we compare turbid and clear lineages to determine if there is evidence of differences in relative rates of substitution. Third, we present evidence of functional diversification and its relationship to the photic environment among cichlid opsin genes. We report statistical evidence of positive selection in all cichlid opsin genes, except short wavelength–sensitive 1 and short wavelength–sensitive 2b. In all genes predicted to be under positive selection, except short wavelength–sensitive 2a, we find differences in selective pressure between turbid and clear lineages. Potential spectral tuning sites are variable among all cichlid opsin genes; however, patterns of substitution consistent with photic environment–driven evolution of opsin genes are observed only for short wavelength–sensitive 1 opsin genes. This study identifies a number of promising candidate-tuning sites for future study by site-directed mutagenesis. This work also begins to demonstrate the molecular evolutionary dynamics of cichlid visual sensitivity and its relationship to the photic environment.

Molecular evolution of primate immunodeficiency viruses and hepatitis delta virus

Primate immunodeficiency viruses, or lentiviruses (HIV-1, HIV-2, and SIV), and hepatitis delta virus (HDV) are RNA viruses characterized by rapid evolution. Infection by primate immunodeficiency viruses usually results in the development of acquired immunodeficiency syndrome (AIDS) in humans and AIDS-like illnesses in Asian macaques. Similarly, hepatitis delta virus infection causes hepatitis and liver cancer in humans. These viruses are heterogeneous within an infected patient and among individuals. Substitution rates in the virus genomes are high and vary in different lineages and among sites. Methods of phylogenetic analysis were applied to study the evolution of primate lentiviruses and the hepatitis delta virus. The following results have been obtained: (1) The substitution rate varies among sites of primate lentivirus genes according to the two parameter gamma distribution, with the shape parameter $\alpha$ being close to 1. (2) Primate immunodeficiency viruses fall into species-specific lineages. Therefore, viral transmissions across primate species are not as frequent as suggested by previous authors. (3) Primate lentiviruses have acquired or lost their pathogenicity several times in the course of evolution. (4) Evidence was provided for multiple infections of a North American patient by distinct HIV-1 strains of the B subtype. (5) Computer simulations indicate that the probability of committing an error in testing HIV transmission depends on the number of virus sequences and their length, the divergence times among sequences, and the model of nucleotide substitution. (6) For future investigations of HIV-1 transmissions, using longer virus sequences and avoiding the use of distant outgroups is recommended. (7) Hepatitis delta virus strains are usually related according to the geographic region of isolation. (8) Evolution of HDV is characterized by the rate of synonymous substitution being lower than the nonsynonymous substitution rate and the rate of evolution of the noncoding region. (9) There is a strong preference for G and C nucleotides at the third codon positions of the HDV coding region. ^

Population and molecular evolution studies on the Melanocortin 1 receptor locus implicate its role in human pigmentation variation

Human pigmentation is a complex trait with the observed variation caused by the varied production of eumelanin (brown/black melanins) and phaeomelanin (red/yellow melanins) by the melanocytes. The melanocortin 1 receptor (MC1R), a G protein-coupled receptor expressed in the melanocytes, is a regulator eu- and phaeomelanin synthesis, and MC1R mutations causing skin and coat color changes are known in many mammals. To understand the role of MC1R in human pigmentation variation, I have sequenced the MC1R gene in 121 individuals sampled from world populations. In addition, I have sequenced the MC1R gene in common and pygmy chimpanzees, gorilla, orangutan, and baboon to study the evolution of MC1R and to infer the ancestral human MC1R sequence. The ancestral MC1R sequence is observed in all 25 African individuals studied, but at lower frequencies in the other populations examined, especially in East and Southeast Asians. The Arg163Gln variant is absent in the Africans studied, almost absent in Europeans, and at a low frequency in Indians, but is at an exceptionally high frequency (70%) in East and Southeast Asians. To further evaluate the role of MC1R variants in human pigmentation variation, I have combined these molecular evolution and population studies with functional assays on MC1R variants and primate MC1Rs. ^

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Molecular evolution of integrins: Genes encoding integrin β subunits from a coral and a sponge

The integrin family of cell surface receptors is strongly conserved in higher animals, but the evolutionary history of integrins is obscure. We have identified and sequenced cDNAs encoding integrin β subunits from a coral (phylum Cnidaria) and a sponge (Porifera), indicating that these proteins existed in the earliest stages of metazoan evolution. The coral βCn1 and, especially, the sponge βPo1 sequences are the most divergent of the “β1-class” integrins and share a number of features not found in any other vertebrate or invertebrate integrins. Perhaps the greatest difference from other β subunits is found in the third and fourth repeats of the cysteine-rich stalk, where the generally conserved spacings between cysteines are highly variable, but not similar, in βCn1 and βPo1. Alternatively spliced cDNAs, containing a stop codon about midway through the full-length translated sequence, were isolated from the sponge library. These cDNAs appear to define a boundary between functional domains, as they would encode a protein that includes the globular ligand-binding head but would be missing the stalk, transmembrane, and cytoplasmic domains. These and other sequence comparisons with vertebrate integrins are discussed with respect to models of integrin structure and function.

Tracking molecular evolution of photosynthesis by characterization of a major photosynthesis gene cluster from Heliobacillus mobilis

A DNA sequence has been obtained for a 35.6-kb genomic segment from Heliobacillus mobilis that contains a major cluster of photosynthesis genes. A total of 30 ORFs were identified, 20 of which encode enzymes for bacteriochlorophyll and carotenoid biosynthesis, reaction-center (RC) apoprotein, and cytochromes for cyclic electron transport. Donor side electron-transfer components to the RC include a putative RC-associated cytochrome c553 and a unique four-large-subunit cytochrome bc complex consisting of Rieske Fe-S protein (encoded by petC), cytochrome b6 (petB), subunit IV (petD), and a diheme cytochrome c (petX). Phylogenetic analysis of various photosynthesis gene products indicates a consistent grouping of oxygenic lineages that are distinct and descendent from anoxygenic lineages. In addition, H. mobilis was placed as the closest relative to cyanobacteria, which form a monophyletic origin to chloroplast-based photosynthetic lineages. The consensus of the photosynthesis gene trees also indicates that purple bacteria are the earliest emerging photosynthetic lineage. Our analysis also indicates that an ancient gene-duplication event giving rise to the paralogous bchI and bchD genes predates the divergence of all photosynthetic groups. In addition, our analysis of gene duplication of the photosystem I and photosystem II core polypeptides supports a “heterologous fusion model” for the origin and evolution of oxygenic photosynthesis.