Prospecção de marcadores microssatélites e análise da estrutura populacional em Atta laevigata (Formicidae: Attini)

Pós-graduação em Ciências Biológicas (Biologia Celular e Molecular) - IBRC

Vertical transmission of hepatitis C virus: A tale of multiple outcomes

Globally, hepatitis C virus (HCV) infection affects approximately 130 million people and 3 million new infections occur annually. HCV is also recognized as an important cause of chronic liver disease in children. The absence of proofreading properties of the HCV RNA polymerase leads to a highly error prone replication process, allowing HCV to escape host immune response. The adaptive nature of HCV evolution dictates the outcome of the disease in many ways. Here, we investigated the molecular evolution of HCV in three unrelated children who acquired chronic HCV infection as a result of mother-to-child transmission, two of whom were also coinfected with HIV-1. The persistence of discrete HCV variants and their population structure were assessed using median joining network and Bayesian approaches. While patterns of viral evolution clearly differed between subjects, immune system dysfunction related to HIV coinfection or persistent HCV seronegativity stand as potential mechanisms to explain the lack of molecular evolution observed in these three cases. In contrast, treatment of HCV infection with PegIFN, which did not lead to sustained virologic responses in all 3 cases, was not associated with commensurate variations in the complexity of the variant spectrum. Finally, the differences in the degree of divergence suggest that the mode of transmission of the virus was not the main factor driving viral evolution. (C) 2013 Elsevier B. V. All rights reserved.

Análise de genes envolvidos na neurotransmissão, formação e manutenção sináptica em indivíduos com transtornos do espectro do Autismo

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Horizontal gene transfer confers fermentative metabolism in the respiratory-deficient plant trypanosomatid Phytomonas serpens

Among trypanosomatids, the genus Phytomonas is the only one specifically adapted to infect plants. These hosts provide a particular habitat with a plentiful supply of carbohydrates. Phytomonas sp. lacks a cytochrome-mediated respiratory chain and Krebs cycle, and ATP production relies predominantly on glycolysis. We have characterised the complete gene encoding a putative pyruvate/indolepyruvate decarboxylase (PDC/IPDC) (548 amino acids) of P. serpens, that displays high amino acid sequence similarity with phytobacteria and Leishmania enzymes. No orthologous PDC/IPDC genes were found in Trypanosoma cruzi or T. brucei. Conservation of the PDC/IPDC gene sequence was verified in 14 Phytomonas isolates. A phylogenetic analysis shows that Phytomonas protein is robustly monophyletic with Leishmania spp. and C. fasciculata enzymes. In the trees this clade appears as a sister group of indolepyruvate decarboxylases of gamma-proteobacteria. This supports the proposition that a horizontal gene transfer event from a donor phytobacteria to a recipient ancestral trypanosome has occurred prior to the separation between Phytomonas. Leishmania and Crithidia. We have measured the PDC activity in P. serpens cell extracts. The enzyme has a Km value for pyruvate of 1.4 mM. The acquisition of a PDC, a key enzyme in alcoholic fermentation, explains earlier observations that ethanol is one of the major end-products of glucose catabolism under aerobic and anaerobic conditions. This represents an alternative and necessary route to reoxidise part of the NADH produced in the highly demanding glycolytic pathway and highlights the importance of this type of event in metabolic adaptation. (C) 2012 Elsevier B.V. All rights reserved.

Avian coronavirus Spike Glycoprotein Ectodomain Shows a Low Codon Adaptation to Gallus gallus with Virus-Exclusive Codons in Strategic Amino Acids Positions

This is a study on the Avian coronavirus IBV and chicken host-relationship from the codon usage point of view based on fifty-nine non-redundant IBV S1 sequences (nt 1-507) from strains detected worldwide and chicken tissue-specific protein genes sequences from IBV-replicating sites. The effective number of codons (ENC) values ranged from 36 to 47.8, indicating a high-to-moderate codon usage bias. The highest IBV codon adaptation index (CAI) value was 0.7, indicating a distant virus versus host synonymous codons usage. The ENC x GC3 % curve indicates that both mutational pressure and natural selection are the driving forces on codon usage pattern in S1. The low CAI values agree with a low S protein expression and considering that S protein is a determinant for attachment and neutralization, this could be a further mechanism besides mRNA transcription attenuation for a low expression of this protein leading to an immune camouflage.

Das Globingen-Repertoire von Amphibien und Teleostiern

Im Rahmen der vorliegenden Dissertation wurde die molekulare Evolution von Globinen in Amphibien und Teleostiern untersucht und Analysen zur Genexpression ausgewählter Globine durchgeführt. Die bisher besonders für die neueren Mitglieder der Superfamilie der Globine – Neuroglobin und Cytoglobin – schwerpunktmäßig in Mammaliern erbrachten Daten sollten durch die Analyse in Amphibien und Teleostiern auf ihre generelle Gültigkeit für Vertebraten überprüft werden. Die Analysen zur Genexpression wurden sowohl in silico, basierend auf genomischen wie EST-Daten, als auch experimentell durch qualitative und quantitative RT-PCR-Nachweise durchgeführt. Die mRNA-Lokalisation wurde durch in situ-Hybridisierungen an Gewebeschnitten beziehungsweise durch Whole mount in situ-Hybridisierung an ganzen Embryonen detektiert. In einem ersten Teil der Arbeit wurde das Globin-Repertoire von Xenopus tropicalis umfassend analysiert. Die Expressionsanalyse der gefundenen Globine umfasste nicht nur adulte Tiere, sonder erstmals auch detailliert die Entwicklungsstadien eines Vertebraten. Dabei wurde festgestellt, dass die vorwiegend neuronale Expression des streng konservierten Neuroglobins ein generelles Charakteristikum aller Tetrapoden ist und bereits in der frühembryonalen Entwicklung auftritt. Auch für das als Einzelkopie im Amphibiengenom vertretene Cytoglobin konnte eine strenge Sequenzkonservierung gezeigt werden. Das Expressionsmuster des Amphibien-Cytoglobins stimmte mit dem aus Mammaliern bekannten überein und zeigte konservierte Charakteristika dieses Globins bei Tetrapoden auf. Die Analyse des Xenopus-Genoms ergab zudem, dass Krallenfrösche nicht über Myoglobin verfügen. Genomische Vergleiche syntäner Genregionen ließen auf Rearrangements in diesem Genombereich im Verlauf der Evolution schließen, in deren Folge das Myoglobingen in den Krallenfröschen deletiert wurde. Die Hämoglobine wurden in Xenopus tropicalis erstmals in einem Amphibium umfassend analysiert. Die Gene zeigten demnach eine geclusterte Anordnung: der tropische Krallenfrosch verfügte über je ein funktionelles α- bzw. β-adultes und sieben bzw. vier α- bzw. β-larvale Hämoglobine, die während der Entwicklung bzw. in adulten Tieren charakteristisch exprimiert wurden. Die Analyse der Hämoglobine hinsichtlich ihrer Lage in einem Cluster, ihrer phylogenetischen Relation zueinander und nicht zuletzt ihres Expressionsmusters ließen Rückschlüsse auf ihre Evolution zu. Zusätzlich zu diesen bereits bekannten Globinen konnte im Rahmen dieser Dissertation das Globingen-Repertoirs von Xenopus um zwei weitere, bisher unbekannte Globine erweitert werden. Diese wurden entsprechend ihrer bisher unbekannten Funktion als GlobinX und GlobinY bezeichnet. Während GlobinY bisher ausschließlich in Amphibien nachgewiesen werden konnte, wurde GlobinX zudem in Teleostiern detektiert und repräsentiert damit ein auf Anamnia beschränktes Globin. Die rekombinante Proteinexpression von Neuroglobin, Cytoglobin, GlobinX und GlobinY des tropischen Krallenfrosches zeigte ein hexakoordiniertes Bindungsschema dieser Globine in ihrem Deoxy-Zustand. In einem zweiten Teil dieser Dissertation wurden Neuroglobin und Cytoglobin in Teleostiern untersucht und die Analyse für diese zwei Gene somit über die Tetrapoden hinaus auf den gesamten Stammbaum der Vertebraten ausgedehnt. Dabei wurde deutlich, dass die vorwiegend neuronale Expression des seit 420 Millionen Jahren streng konservierten Neuroglobins ein generelles Merkmal dieses Globins in allen Vertebraten ist. Der in Amphibien und Teleostiern erbrachte und mit Ergebnissen in Mammaliern übereinstimmende Nachweis von Neuroglobin in neuronalen Geweben mit einem hohen Stoffwechsel lässt derzeit eine Funktion dieses Globins im Sauerstoffmetabolimus als wahrscheinlich erscheinen. Ob Neuroglobin dabei als kurzzeitiger Sauerstoffspeicher, O2-Transoprter oder aber in der Detoxifikation reaktiver Sauerstoff- bzw. Stickstoffspezies agiert, bleibt zu untersuchen. Für Cytoglobin konnte eine offenbar alle Teleostier betreffende Genduplikation nachgewiesen werden. Phylogenetische Analysen zeigen die Monophylie der Vertebraten-Cytoglobine. Der Vergleich der paralogen Cytoglobine der Teleostier mit dem syntänen Genombereich des humanen Cytoglobins zeigte die wahrscheinliche Entstehung der Fisch-Cytoglobine durch eine Genomduplikation in einem Vorfahren aller Teleostier vor etwa 300-450 Millionen Jahren. Die paralogen Cytoglobine zeigten in Danio rerio und Tetraodon nigroviridis differierende, charakteristische Expressionsmuster, die mit der Theorie der Subfunktionalisierung von Genen in Folge eines Duplikationsereignisses kompatibel sind. Die Analyse zeigte, dass Cygb-1 prädominant in Gehirn und Herz exprimiert wurde, Cygb-2 hingegen bevorzugt in Gehirn und Auge. Dies bestätigte indirekt die Hypothese, nach der das Cytoglobin der Mammalier zwei unterschiedliche Funktionen in differenten Geweben wahrnimmt. Die rekombinante Expression von Cygb-1 des Zebrabärblings zeigte zudem, das auch dieses Globin in seiner Deoxy-Form über ein hexakoordiniertes Bindungsschema verfügt.

Caratterizzazione genomica di virus influenzali suini in Italia mediante next generation sequencing

I sottotipi H1N1, H1N2 e H3N2 di influenza A virus sono largamente diffusi nella popolazione suina di tutto il mondo. Nel presente lavoro è stato sviluppato un protocollo di sequenziamento di c.d. nuova generazione, su piattaforma Ion Torrent PGM, idoneo per l’analisi di tutti i virus influenzali suini (SIV). Per valutare l’evoluzione molecolare dei SIV italiani, sono stati sequenziati ed analizzati mediante analisi genomica e filogenetica un totale di sessantadue ceppi di SIV appartenenti ai sottotipi H1N1, H1N2 e H3N2, isolati in Italia dal 1998 al 2014. Sono stati evidenziati in sei campioni due fenomeni di riassortimento: tutti i SIV H1N2 esaminati presentavano una neuraminidasi di derivazione umana, diversa da quella dei SIV H1N2 circolanti in Europa, inoltre l’emoagglutinina (HA) di due isolati H1N2 era originata dal riassortimento con un SIV H1N1 avian-like. L’analisi molecolare dell’HA ha permesso di rivelare un’inserzione di due amminoacidi in quattro SIV H1N1 pandemici e una delezione di due aminoacidi in quattro SIV H1N2, entrambe a livello del sito di legame con il recettore cellulare. E’ stata inoltre evidenziata un’elevata omologia di un SIV H1N1 con ceppi europei isolati negli anni ’80, suggerendo la possibile origine vaccinale di questo virus. E’ stato possibile, in aggiunta, applicare il nuovo protocollo sviluppato per sequenziare un virus influenzale aviare altamente patogeno trasmesso all’uomo, direttamente da campione biologico. La diversità genetica nei SIV esaminati in questo studio conferma l’importanza di un continuo monitoraggio della costellazione genomica dei virus influenzali nella popolazione suina.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.

Patterns of positive selection in seven ant genomes

The evolution of ants is marked by remarkable adaptations that allowed the development of very complex social systems. To identify how ant-specific adaptations are associated with patterns of molecular evolution, we searched for signs of positive selection on amino-acid changes in proteins. We identified 24 functional categories of genes which were enriched for positively selected genes in the ant lineage. We also reanalyzed genome-wide datasets in bees and flies with the same methodology, to check whether positive selection was specific to ants or also present in other insects. Notably, genes implicated in immunity were enriched for positively selected genes in the three lineages, ruling out the hypothesis that the evolution of hygienic behaviors in social insects caused a major relaxation of selective pressure on immune genes. Our scan also indicated that genes implicated in neurogenesis and olfaction started to undergo increased positive selection before the evolution of sociality in Hymenoptera. Finally, the comparison between these three lineages allowed us to pinpoint molecular evolution patterns that were specific to the ant lineage. In particular, there was ant-specific recurrent positive selection on genes with mitochondrial functions, suggesting that mitochondrial activity was improved during the evolution of this lineage. This might have been an important step toward the evolution of extreme lifespan that is a hallmark of ants.

Primary structure and evolutionary relationship between the adult alpha-globin genes and their 5'-flanking regions of Xenopus laevis and Xenopus tropicalis

To investigate the evolution of globin genes in the genus Xenopus, we have determined the primary structure of the related adult alpha I- and alpha II-globin genes of X. laevis and of the adult alpha-globin gene of X. tropicalis, including their 5'-flanking regions. All three genes are comprised of three exons and two introns at homologous positions. The exons are highly conserved and code for 141 amino acids. By contrast, the corresponding introns vary in length and show considerable divergence. Comparison of 900 bp of the 5'-flanking region revealed that the X. tropicalis gene contains a conserved proximal 310-bp promoter sequence, comprised of the canonical TATA and CCAAT motifs at homologous positions, and five conserved elements in the same order and at similar positions as previously shown for the corresponding genes of X. laevis. We therefore conclude that these conserved upstream elements may represent regulatory sequences for cell-specific regulation of the adult Xenopus globin genes.

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

Detecting the signature of natural selection with microsatellites

Natural selection is one of the major factors in the evolution of all organisms. Detecting the signature of natural selection has been a central theme in evolutionary genetics. With the availability of microsatellite data, it is of interest to study how natural selection can be detected with microsatellites. ^ The overall aim of this research is to detect signatures of natural selection with data on genetic variation at microsatellite loci. The null hypothesis to be tested is the neutral mutation theory of molecular evolution, which states that different alleles at a locus have equivalent effects on fitness. Currently used tests of this hypothesis based on data on genetic polymorphism in natural populations presume that mutations at the loci follow the infinite allele/site models (IAM, ISM), in the sense that at each site at most only one mutation event is recorded, and each mutation leads to an allele not seen before in the population. Microsatellite loci, which are abundant in the genome, do not obey these mutation models, since the new alleles at such loci can be created either by contraction or expansion of tandem repeat sizes of core motifs. Since the current genome map is mainly composed of microsatellite loci and this class of loci is still most commonly studied in the context of human genome diversity, this research explores how the current test procedures for testing the neutral mutation hypothesis should be modified to take into account a generalized model of forward-backward stepwise mutations. In addition, recent literature also suggested that past demographic history of populations, presence of population substructure, and varying rates of mutations across loci all have confounding effects for detecting signatures of natural selection. ^ The effects of the stepwise mutation model and other confounding factors on detecting signature of natural selection are the main results of the research. ^

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.