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. ^

Piecing together evolution of the vertebrate endocrine system

One of the great challenges in biology is to understand how particular complex morphological and physiological characters originated in specific evolutionary lineages. In this article, we address the origin of the vertebrate hypothalamic-pituitary-peripheral gland (H-P-PG) endocrine system, a complex network of specialized tissues, ligands and receptors. Analysis of metazoan nucleotide and protein sequences reveals a patchwork pattern of H-P-PG gene conservation between vertebrates and closely related invertebrates (ascidians). This is consistent with a model of how the vertebrate H-P-PG endocrine system could have emerged in relatively few steps by gene family expansion and by regulatory and structural modifications to genes that are present in a chordate ancestor. Some of these changes might have resulted in new connections between metabolic or signaling pathways, such as the bridging of 'synthesis islands' to form an efficient system for steroid hormone synthesis.

Phylogeny, molecular and fossil dating, and biogeographic history of Annonaceae and Myristicaceae (Magnoliales)

Annonaceae and Myristicaceae, the two largest families of Magnoliales, are pantropical groups of uncertain geographic history. The most recent morphological and molecular phylogenetic analyses identify the Asian-American genus Anaxagorea as sister to all other Annonaceae and the ambavioids, consisting of small genera endemic to South America, Africa, Madagascar, and Asia, as a second branch. However, most genera form a large clade in which the basal lines are African, and South American and Asian taxa are more deeply nested. Although it has been suggested that Anaxagorea was an ancient Laurasian line, present data indicate that this genus is basically South American. These considerations may mean that the family as a whole began its radiation in Africa and South America in the Late Cretaceous, when the South Atlantic was narrower, and several lines dispersed from Africa-Madagascar into Laurasia as the Tethys closed in the Tertiary. This scenario is consistent with the occurrence of annonaceous seeds in the latest Cretaceous of Nigeria and the Eocene of England and with molecular dating of the family. Based on distribution of putatively primitive taxa in Madagascar and derived taxa in Asia, it has been suggested that Myristicaceae had a similar history. Phylogenetic analyses of Myristicaceae, using morphology and several plastid regions, confirm that the ancestral area was Africa-Madagascar and that Asian taxa are derived. However, Myristicaceae as a whole show strikingly lower molecular divergence than Annonaceae, indicating either a much younger age or a marked slowdown in molecular evolution. The fact that the oldest diagnostic fossils of Myristicaceae are Miocene seeds might be taken as evidence that Myristicaceae are much younger than Annonaceae, but this is implausible in requiring transoceanic dispersal of their large, animal-dispersed seeds.

Limited effect of anthropogenic habitat fragmentation on molecular diversity in a rain forest skink, Gnypetoscincus queenslandiae

To examine the effects of recent habitat fragmentation, we assayed genetic diversity in a rain forest endemic lizard, the prickly forest skink (Gnypetoscincus queenslandiae), from seven forest fragments and five sites in continuous forest on the Atherton tableland of northeastern Queensland, Australia. The rain forest in this region was fragmented by logging and clearing for dairy farms in the early 1900s and most forest fragments studied have been isolated for 50-80 years or nine to 12 skink generations. We genotyped 411 individuals at nine microsatellite DNA loci and found fewer alleles per locus in prickly forest skinks from small rain forest fragments and a lower ratio of allele number to allele size range in forest fragments than in continuous forest, indicative of a decrease in effective population size. In contrast, and as expected for populations with small neighbourhood sizes, neither heterozygosity nor variance in allele size differed between fragments and sites in continuous forests. Considering measures of among population differentiation, there was no increase in F-ST among fragments and a significant isolation by distance pattern was identified across all 12 sites. However, the relationship between genetic (F-ST) and geographical distance was significantly stronger for continuous forest sites than for fragments, consistent with disruption of gene flow among the latter. The observed changes in genetic diversity within and among populations are small, but in the direction predicted by the theory of genetic erosion in recently fragmented populations. The results also illustrate the inherent difficulty in detecting genetic consequences of recent habitat fragmentation, even in genetically variable species, and especially when effective population size and dispersal rates are low.

Speciation of Antechinus stuartii and A. subtropicus (Marsupialia : Dasyuridae) in eastern Australia: molecular and morphological evidence

This paper evaluates the systematic status of the Antechinus populations of northern New South Wales and southern Queensland, using a combined morphological and molecular (allozymes and mitochondrial DNA) approach. Analysis of the d-loop section of the mitochondrial DNA control region revealed two highly supported clades within A. stuartii sensu lato that were sympatric in the Border Ranges of northern New South Wales. However, genetic distances between these clades were small ( approximately 3%), indicating that time of divergence was probably relatively recent. Allozyme electrophoresis also showed very small differences between clades/ species. Analyses of cranial and dental characters showed that the members of each of these clades differed morphologically and that the clades corresponded to A. stuartii and the recently described A. subtropicus. The combined results support the species status of A. stuartii and A. subtropicus, and suggest that speciation was likely a result of a recent vicariant event.

Molecular tools for studying puberty in the grey mullet, Mugil cephalus: cloning of cDNAs that regulate reproductive function

The cDNAs coding for the brain GnRHs (AY373449-51), pituitary GH, SL and PRL, and liver IGFs (AY427954-5) were isolated. Partial cDNA sequences of the brain (Cyp19b) and gonadal (Cyp19a) aromatases have also been obtained. These tools would be utilized to study the endocrine regulation of puberty in the grey mullet.