Molecular cloning of natriuretic peptides from the heart of reptiles : loss of ANP in diapsid reptiles and birds

Atrial natriuretic peptide (ANP) and B-type NP (BNP) are hormones involved in homeostatic control of body fluid and cardiovascular regulation. Both ANP and BNP have been cloned from the heart of mammals, amphibians, and teleost fishes, while an additional cardiac peptide, ventricular NP, has been found in selected species of teleost fish. However, in chicken, BNP is the primary cardiac peptide identified thus far. In contrast, the types of NP/s present in the reptilian heart are unknown, representing a considerable gap in our understanding of NP evolution. In the present study, we cloned and sequenced a BNP cDNA from the atria of representative species of reptile, including crocodile, lizard, snake, and tortoise. In addition, we cloned BNP from the pigeon atria. The reptilian and pigeon BNP cDNAs had ATTTA repeats in the 3′ untranslated region, as observed in all vertebrate BNP mRNAs. A high sequence homology was evident when comparing reptile and pigeon preproBNP with the previously identified chicken preproBNP. In particular, the predicted mature BNP-29 was identical between crocodile, tortoise, and chicken, with pigeon having a single amino acid substitution; lizard and snake BNP had seven and nine substitutions, respectively. Furthermore, an ANP cDNA could only be cloned from the tortoise atria. Since ANP was not isolated from the heart of any non-chelonian reptile and appears to be absent in birds, we propose that the ANP gene has been lost after branching of the turtles in the amniote line. This data provides new avenues for research on NP function in reptiles.

Molecular biology of natriuretic peptides in reptiles and birds

The heart produces natriuretic peptides that are critical regulators of blood pressure and renal function. This study examined the molecular evolution of natriuretic peptides in vertebrates and discovered novel forms of the peptides in birds. The research outcomes advanced our knowledge of the importance of these peptides in animal physiology.

Evolution of lactation : ancient origin and extreme adaptations of the lactation system

Lactation, an important characteristic of mammalian reproduction, has evolved by exploiting a diversity of strategies across mammals. Comparative genomics and transcriptomics experiments have now allowed a more in-depth analysis of the molecular evolution of lactation. Milk cell and mammary gland genomic studies have started to reveal conserved milk proteins and other components of the lactation system of monotreme, marsupial, and eutherian lineages. These analyses confirm the ancient origin of the lactation system and provide useful insight into the function of specific milk proteins in the control of lactation. These studies also illuminate the role of milk in the regulation of growth and development of the young beyond simple nutritive aspects.

Molecular characterization and phylogenetic analysis of JussuMP-I: A RGD-P-III class hemorrhagic metalloprotease from Bothrops jararacussu snake venom

Snake venom metalloproteases (SVMPs) embody zinc-dependent multidomain enzymes responsible for a relevant pathophysiology in envenomation. including local and systemic hemorrhage. The molecular features responsible for hemorrhagic potency of SVMPs have been associated with their multidomains structures which can target these proteins them to several receptors of different tissues and cellular types. BjussuMP-I. a SVMP isolated from the Bothrops jararacussu venom, has been characterized as a P-III hemorrhagic metalloprotease. The complete cDNA sequence of BjussuMP-I with 1641bp encodes open reading frames of 547 amino acid residues, which conserve the common domains of P-III high molecular weight hemorrhagic metalloproteases: (i) pre-pro-peptide, (ii) metalloprotease, (iii) disintegrin-like and (iv) rich cysteine domain. BjussuMP-I induced lyses in fibrin clots and inhibited collagen- and ADP-induced platelet aggregation. We are reporting, for the first time, the primary structure of an RGD-P-III class snake venom metalloprotease. A phylogenetic analysis of the BjussuMP-1 metalloprotease/catalytic domain was performed to get new insights into the molecular evolution of the metalloproteases. A theoretical molecular model of this domain was built through folding recognition (threading) techniques and refined by molecular dynamics simulation. Then, the final BjussuMP-I catalytic domain model was compared to other SVMPs and Reprolysin family proteins in order to identify eventual structural differences, which could help to understand the biochemical activities of these enzymes. The presence of large hydrophobic areas and some conserved surface charge-positive residues were identified as important features of the SVMPs and other metalloproteases. (C) 2006 Elsevier B.V. All rights reserved.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

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

Molecular phylogenetic relationships and phenotypic diversity in miniaturized toadlets, genus Brachycephalus (Amphibia: Anura: Brachycephalidae)

Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18. mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate. © 2011 Elsevier Inc.

Chromosome mapping of repetitive sequences in Anostomidae species: Implications for genomic and sex chromosome evolution

Background: Members of the Anostomidae family provide an interesting model system for the study of the influence of repetitive elements on genome composition, mainly because they possess numerous heterochromatic segments and a peculiar system of female heterogamety that is restricted to a few species of the Leporinus genus. The aim of this study was to isolate and identify important new repetitive DNA elements in Anostomidae through restriction enzyme digestion, followed by cloning, characterisation and chromosome mapping of this fragment. To identify repetitive elements in other Leporinus species and expand on studies of repetitive elements in Anostomidae, hybridisation experiments were also performed using previously described probes of LeSpeI repetitive elements. Results: The 628-base pair (bp) LeSpeII fragment was hybridised to metaphase cells of L. elongatus individuals as well as those of L. macrocephalus, L. obtusidens, L. striatus, L. lacustris, L. friderici, Schizodon borellii and S. isognathus. In L. elongatus, both male and female cells contained small clusters of LeSpeII repetitive elements dispersed on all of the chromosomes, with enrichment near most of the terminal portions of the chromosomes. In the female sex chromosomes of L. elongatus (Z2,Z2/W1W 2), however, this repeated element was absent. In the remaining species, a dispersed pattern of hybridisation was observed on all chromosomes irrespective of whether or not they were sex chromosomes. The repetitive element LeSpeI produced positive hybridisations signals only in L. elongatus, L. macrocephalus and L. obtusidens, i.e., species with differentiated sex chromosomes. In the remaining species, the LeSpeI element did not produce hybridisation signals. Conclusions: Results are discussed in terms of the effects of repetitive sequences on the differentiation of the Anostomidae genome, especially with respect to sex chromosome evolution. LeSpeII showed hybridisation patterns typical of Long Interspersed Elements (LINEs). The differential distribution of this element may be linked to sex chromosome differentiation in L. elongatus species. The relationship between sex chromosome specificity and the LeSpeI element is confirmed in the species L. elongatus, L. macrocephalus and L. obtusidens. © 2012 da Silva et al.; licensee BioMed Central Ltd.

Circulation of Different Lineages of Dengue Virus 2, Genotype American/Asian in Brazil: Dynamics and Molecular and Phylogenetic Characterization

The American/Asian genotype of Dengue virus type 2 (DENV-2) was introduced into the Americas in the 80′s. Although there is no data showing when this genotype was first introduced into Brazil, it was first detected in Brazil in 1990. After which the virus spread throughout the country and major epidemics occurred in 1998, 2007/08 and 2010. In this study we sequenced 12 DENV-2 genomes obtained from serum samples of patients with dengue fever residing in São José do Rio Preto, São Paulo (SJRP/SP), Brazil, in 2008. The whole open reading frame or envelope sequences were used to perform phylogenetic, phylogeographic and evolutionary analyses. Isolates from SJRP/SP were grouped within one lineage (BR3) close to isolates from Rio de Janeiro, Brazil. Isolates from SJRP were probably introduced there at least in 2007, prior to its detection in the 2008 outbreak. DENV-2 circulation in Brazil is characterized by the introduction, displacement and circulation of three well-defined lineages in different times, most probably from the Caribbean. Thirty-seven unique amino acid substitutions were observed among the lineages, including seven amino acid differences in domains I to III of the envelope protein. Moreover, we dated here, for the first time, the introduction of American/Asian genotype into Brazil (lineage BR1) to 1988/89, followed by the introduction of lineages BR2 (1998-2000) and BR3 (2003-05). Our results show a delay between the introduction and detection of DENV-2 lineages in Brazil, reinforcing the importance and need for surveillance programs to detect and trace the evolution of these viruses. Additionally, Brazilian DENV-2 differed in genetic diversity, date of introduction and geographic origin and distribution in Brazil, and these are important factors for the evolution, dynamics and control of dengue. © 2013 Drumond et al.

On Hepatitis C Virus Evolution: The Interaction between Virus and Host towards Treatment Outcome

Background:Hepatitis C is a disease spread throughout the world. Hepatitis C virus (HCV), the etiological agent of this disease, is a single-stranded positive RNA virus. Its genome encodes a single precursor protein that yields ten proteins after processing. NS5A, one of the non-structural viral proteins, is most associated with interferon-based therapy response, the approved treatment for hepatitis C in Brazil. HCV has a high mutation rate and therefore high variability, which may be important for evading the immune system and response to therapy. The aim of this study was to analyze the evolution of NS5A quasispecies before, during, and after treatment in patients infected with HCV genotype 3a who presented different therapy responses.Methods:Viral RNA was extracted, cDNA was synthesized, the NS5A region was amplified and cloned, and 15 clones from each time-point were sequenced. The sequences were analyzed for evolutionary history, genetic diversity and selection.Results:This analysis shows that the viral population that persists after treatment for most non-responder patients is present in before-treatment samples, suggesting it is adapted to evade treatment. In contrast, the population found in before treatment samples from most end-of-treatment responder patients either are selected out or appears in low frequency after relapse, therefore changing the population structure. The exceptions illustrate the uniqueness of the evolutionary process, and therefore the treatment resistance process, in each patient.Conclusion:Although evolutionary behavior throughout treatment showed that each patient presented different population dynamics unrelated to therapy outcome, it seems that the viral population from non-responders that resists the treatment already had strains that could evade therapy before it started. © 2013 Bittar et al.

PRP8 intein in cryptic species of Histoplasma capsulatum: Evolution and phylogeny

The PRP8 intein is the most widespread intein among the Kingdom Fungi. This genetic element occurs within the prp8 gene, and is transcribed and translated simultaneously with the gene. After translation, the intein excises itself from the Prp8 protein by an autocatalytic splicing reaction, subsequently joining the N and C terminals of the host protein, which retains its functional conformation. Besides the splicing domain, some PRP8 inteins also have a homing endonuclease (HE) domain which, if functional, makes the intein a mobile element capable of becoming fixed in a population. This work aimed to study (1) The occurrence of this intein in Histoplasma capsulatum isolates (n=. 99) belonging to different cryptic species collected in diverse geographical locations, and (2) The functionality of the endonuclease domains of H. capsulatum PRP8 inteins and their phylogenetic relationship among the cryptic species. Our results suggest that the PRP8 intein is fixed in H. capsulatum populations and that an admixture or a probable ancestral polymorphism of the PRP8 intein sequences is responsible for the apparent paraphyletic pattern of the LAmA clade which, in the intein phylogeny, also encompasses sequences from LAmB isolates. The PRP8 intein sequences clearly separate the different cryptic species, and may serve as an additional molecular typing tool, as previously proposed for other fungi genus, such as Cryptococcus and Paracoccidioides. © 2013 Elsevier B.V.

Evolution of the Metazoan Mitochondrial Replicase

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

Epidemiologia molecular do vírus da Hepatite C: análise comparativa de diferentes regiões subgenômicas aplicadas a estudos de associação genética

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A molecular phylogeny of bivalve mollusks: ancient radiations and divergences as revealed by mitochondrial genes

The main scope of my PhD is the reconstruction of the large-scale bivalve phylogeny on the basis of four mitochondrial genes, with samples taken from all major groups of the class. To my knowledge, it is the first attempt of such a breadth in Bivalvia. I decided to focus on both ribosomal and protein coding DNA sequences (two ribosomal encoding genes -12s and 16s -, and two protein coding ones - cytochrome c oxidase I and cytochrome b), since either bibliography and my preliminary results confirmed the importance of combined gene signals in improving evolutionary pathways of the group. Moreover, I wanted to propose a methodological pipeline that proved to be useful to obtain robust results in bivalves phylogeny. Actually, best-performing taxon sampling and alignment strategies were tested, and several data partitioning and molecular evolution models were analyzed, thus demonstrating the importance of molding and implementing non-trivial evolutionary models. In the line of a more rigorous approach to data analysis, I also proposed a new method to assess taxon sampling, by developing Clarke and Warwick statistics: taxon sampling is a major concern in phylogenetic studies, and incomplete, biased, or improper taxon assemblies can lead to misleading results in reconstructing evolutionary trees. Theoretical methods are already available to optimize taxon choice in phylogenetic analyses, but most involve some knowledge about genetic relationships of the group of interest, or even a well-established phylogeny itself; these data are not always available in general phylogenetic applications. The method I proposed measures the "phylogenetic representativeness" of a given sample or set of samples and it is based entirely on the pre-existing available taxonomy of the ingroup, which is commonly known to investigators. Moreover, it also accounts for instability and discordance in taxonomies. A Python-based script suite, called PhyRe, has been developed to implement all analyses.

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

Evolution of cis-regulatory elements in a repetitive sequence adjacent to a sea urchin aboral ectoderm -specific gene

The creation, preservation, and degeneration of cis-regulatory elements controlling developmental gene expression are fundamental genome-level evolutionary processes about which little is known. In this study, critical differences in cis-regulatory elements controlling the expression of the sea urchin aboral ectoderm-specific spec genes were identified and explored. In genomes of species within the Strongylocentrotidae family, multiple copies of a repetitive sequence element termed RSR were present, but RSRs were not detected in genomes of species outside Strongylocentrotidae. RSRs are invariably associated with spec genes, and in Strongylocentrotus purpuratus, the spec2a RSR functioned as a transcriptional enhancer displaying greater activity than RSRs from the spec1 or spec2c paralogs. Single base-pair differences at two cis-regulatory elements within the spec2a RSR greatly increased the binding affinities of four transcription factors: SpCCAAT-binding factor at one element and SpOtx, SpGoosecoid, and SpGATA-E at another. The cis-regulatory elements to which SpCCAAT-binding factor, SpOtx, SpGoosecoid, and SpGATA-E bound were recent evolutionary acquisitions that could act either to activate or repress transcription, depending on the cell type. These elements were found in the spec2a RSR ortholog in Strongylocentrotus pallidus but not in the RSR orthologs of Strongylocentrotus droebachiensis or Hemicentrotus pulcherrimus. These results indicate that spec genes exhibit a dynamic pattern of cis-regulatory element evolution while stabilizing selection preserves their aboral ectoderm expression domain. ^