Gulliver, a long terminal repeat retrotransposon from the genome of the oriental blood fluke Schistosoma japonicum

We characterized the consensus sequence and structure of a long terminal repeat (LTR) retrotransposon from the genome of the human blood fluke, Schistosoma japonicum, and have earned this element, Gulliver. The full length, consensus Gulliver LTR retrotransposon was 4788 bp, and it was flanked at its 5'- and 3'-ends by LTRs of 259 bp. Each LTR included RNA polymerase II promoter sequences, a CAAT signal and a TATA box, Gulliver exhibited features characteristic of a functional LTR retrotransposon including two read through (termination) ORFs encoding retroviral gag and pol proteins of 312 and 1071 amino acid residues, respectively. The gag ORF encoded motifs conserved in nucleic acid binding proteins, while the pol ORF encoded conserved domains of aspartic protease, reverse transcriptase (RT), RNaseH and integrase, in that order, a pol pattern conserved in the gypsy lineage of LTR retrotransposons. Whereas the sequence and structure of Gulliver was similar to that of gypsy, phylogenetic analysis revealed that Gulliver did not group particularly closely with the gypsy family. Rather, its closest relatives were a LTR retrotransposon from Caenorhabditis elegans, mag from Bombyx mori and, to a lesser extent, easel from the salmon Oncorhynchus keta. Dot blot hybridizations indicated that Gulliver was present at between 100 and several thousand copies in the S. japonicum genome, and Southern hybridization analysis suggested its probable presence in the genome of Schistosoma mansoni. Transcripts encoding the RT domain of Gulliver were detected by RT-PCR in larval and adult stages of S. japonicum, indicating that (at least) the RT domain of Gulliver is transcribed. This is the first report of the sequence and structure of an LTR retrotransposon from any schistosome or indeed from any species belonging to the phylum Platyhelminthes. (C) 2001 Elsevier Science B.V. All rights reserved.

Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda)

The complete nucleotide sequence of the mitochondrial (mt) DNA molecule of the liverfluke, Fasciola hepatica (phylum Platyhelminthes, class Trematoda, family Fasciolidae), was determined, It comprises 14462 bp, contains 12 protein-encoding, 2 ribosomal and 22 transfer RNA genes, and is the second complete flatworm (and the first trematode) mitochondrial sequence to be described in detail. All of the genes are transcribed from the same strand. Of the genes typically found in mitochondrial genomes of eumetazoans, only atp8 is absent. The nad4L and nad4 genes overlap by 40 nt. Most intergenic sequences are very short. Two larger non-coding regions are present. The longer one (817 nt) is located between trnG and cox3 and consists of 8 identical tandem repeats of 85 nt, rich in G and C, followed by 1 imperfect repeat. The shorter non-coding region (187 nt) exhibits no special features and is separated from the longer region by trnG. The gene arrangement resembles that of some other trematodes including the eastern Asian Schistosoma species (and cyclophyllidean cestode species) but it is strikingly different from that of the African schistosomes, represented by Schistosoma mansoni. The genetic code is as inferred previously for flatworms. Transfer RNA genes range in length from 58 to 70 nt, their products producing characteristic 'clover leaf' structures, except for tRNA(S-VON) and tRNA(S-AGN) lacking the DHU arm.

The draft genome of Ciona intestinalis: Insights into chordate and vertebrate origins

The first chordates appear in the fossil record at the time of the Cambrian explosion, nearly 550 million years ago. The modern ascidian tadpole represents a plausible approximation to these ancestral chordates. To illuminate the origins of chordate and vertebrates, we generated a draft of the protein-coding portion of the genome of the most studied ascidian, Ciona intestinalis. The Ciona genome contains similar to16,000 protein-coding genes, similar to the number in other invertebrates, but only half that found in vertebrates. Vertebrate gene families are typically found in simplified form in Ciona, suggesting that ascidians contain the basic ancestral complement of genes involved in cell signaling and development. The ascidian genome has also acquired a number of lineage-specific innovations, including a group of genes engaged in cellulose metabolism that are related to those in bacteria and fungi.

The human zoo: Endogenous retroviruses in the human genome

The main focus of the human genome sequencing project has been gene discovery, but a great additional benefit is that it offers the chance to examine the large proportion of the genome that does not contain human genes. The nature of this ‘noncoding’ DNA is poorly understood, both as an evolutionary question (how did it get there?) and in the functional sense (what is it doing now?). Much of the noncoding DNA is derived from retroviruses that have inserted their DNA into the genome. The availability of complete genomic sequences will revolutionize studies of the number and location of endogenous retroviruses, their role in genome evolution, and their contribution to human disease.

Inhibition of tubulin assembly and covalent binding to microtubular protein by valproic acid glucuronide in vitro

Acyl glucuronides are reactive metabolites of carboxylate drugs, able to undergo a number of reactions in vitro and in vivo, including isomerization via intramolecular rearrangement and covalent adduct formation with proteins. The intrinsic reactivity of a particular acyl glucuronide depends upon the chemical makeup of the drug moiety. The least reactive acyl glucuronide yet reported is valproic acid acyl glucuronide (VPA-G), which is the major metabolite of the antiepileptic agent valproic acid (VPA). In this study, we showed that both VPA-G and its rearrangement isomers (iso-VPA-G) interacted with bovine brain microtubular protein (MTP, comprised of 85% tubulin and 15% microtubule associated proteins {MAPs}). MTP was incubated with VPA, VPA-G and iso-VPA-G for 2 h at room temperature and pH 7.5 at various concentrations up to 4 mM. VPA-G and iso-VPA-G caused dose-dependent inhibition of assembly of MTP into microtubules, with 50% inhibition (IC50) values of 1.0 and 0.2 mM respectively, suggesting that iso-VPA-G has five times more inhibitory potential than VPA-G. VPA itself did not inhibit microtubule formation except at very high concentrations (greater than or equal to2 mM). Dialysis to remove unbound VPA-G and iso-VPA-G (prior to the assembly assay) diminished inhibition while not removing it. Comparison of covalent binding of VPA-G and iso-VPA-G (using [C-14]-labelled species) showed that adduct formation was much greater for iso-vTA-G. When [C-14]-iso-VPA-G was reacted with MTP in the presence of sodium cyanide (to stabilize glycation adducts), subsequent separation into tubulin and MAPs fractions by ion exchange chromatography revealed that 78 and 22% of the covalent binding occurred with the MAPs and tubulin fractions respectively. These experiments support the notion of both covalent and reversible binding playing parts in the inhibition of microtubule formation from MTP (though the acyl glucuronide of VPA is less important than its rearrangement isomers in this regard), and that both tubulin and (perhaps more importantly) MAPs form adducts with acyl glucuronides. (C) 2002 Elsevier Science Inc. All rights reserved.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

Cadherin-directed actin assembly: E-cadherin physically associates with the Arp2/3 complex to direct actin assembly in nascent adhesive contacts

Cadherin cell adhesion molecules are major determinants of tissue patterning which function in cooperation with the actin cytoskeleton [1-4]. In the context of stable adhesion [1], cadherin/catenin complexes are often envisaged to passively scaffold onto cortical actin filaments. However, cadherins also form dynamic adhesive contacts during wound healing and morphogenesis [2]. Here actin polymerization has been proposed to drive cell surfaces together [5], although F-actin reorganization also occurs as cell contacts mature [6]. The interaction between cadherins and actin is therefore likely to depend on the functional state of adhesion. We sought to analyze the relationship between cadherin homophilic binding and cytoskeletal activity during early cadherin adhesive contacts. Dissecting the specific effect of cadherin ligation alone on actin regulation is difficult in native cell-cell contacts, due to the range of juxtacrine signals that can arise when two cell surfaces adhere [7]. We therefore activated homophilic ligation using a specific functional recombinant protein. We report the first evidence that E-cadherin associates with the Arp2/3 complex actin nucleator and demonstrate that cadherin binding can exert an active, instructive influence on cells to mark sites for actin assembly at the cell surface.

pido, a non-long terminal repeat retrotransposon of the chicken repeat 1 family from the genome of the Oriental blood fluke, Schistosoma japonicum

A newly described non-long terminal repeat (non-LTR) retrotransposon element was isolated from the genome of the Oriental schistosome, Schistosoma japonicum. At least 1000 partial copies of the element, which was named pido, were dispersed throughout the genome of S. japonicum. As is usual with non-LTR retrotransposons, it is expected that many pido elements will be 5'-truncated. A consensus sequence of 3564 bp of the truncated pido element was assembled from several genomic fragments that contained pido-hybridizing sequences. The sequence encoded part of the first open reading frame (ORF), the entire second ORF and, at its 3'-terminus, a tandemly repetitive, A-rich (TA(6)TA(5)TA(8)) tail, The ORF1 of pido encoded a nucleic acid binding protein and ORF2 encoded a retroviral-like polyprotein that included apurinic/apyrimidinic endonuclease (EN) and reverse transcriptase (RT) domains, in that order. Based on its sequence and structure, and phylogenetic analyses of both the RT and EN domains, pido belongs to the chicken repeat 1 (CR1)-like lineage of elements known from the chicken, turtle, puffer fish, mosquitoes and other taxa. pido shared equal similarity with CRI from chicken, an uncharacterized retrotransposon from Caenorhabditis elegans and SR1 (a non-LTR retrotransposon) from the related blood fluke Schistosoma mansoni; the level of similarity between pido and SR1 indicated that these two schistosome retrotransposons were related but not orthologous. The findings indicate that schistosomes have been colonized by at least two discrete CRI-like elements. Whereas pido did not appear to have a tight target site specificity, at least one copy of pido has inserted into the 3'-untranslated region of a protein-encoding gene (GeriBank AW736757) of as yet unknown identity. mRNA encoding the RT of pido was detected by reverse transcription-polymerase chain reaction in the egg, miracidium. and adult developmental stages of S. japonicum, indicating that the RT domain was transcribed and suggesting that pido was replicating actively and mobile within the S. japonicum genome. (C) 2002 Elsevier Science B.V. All rights reserved.

Molecular and functional analyses of Kunjin virus infectious cDNA clones demonstrate the essential roles for NS2A in virus assembly and for a nonconservative residue in NS3 in RNA replication

A number of full-length cDNA clones of Kunjin virus (KUN) were previously prepared; it was shown that two of them, pAKUN and FLSDX, differed in specific infectivities of corresponding in vitro transcribed RNAs by similar to100,000-fold (A. A. Khromykh et al., J. Virol. 72:7270-7279, 1998). In this study, we analyzed a possible genetic determinant(s) of the observed differences in infectivity initially by sequencing the entire cDNAs of both clones and comparing them with the published sequence of the parental KUN strain MRM61C. We found six common amino acid residues in both cDNA clones that were different from those in the published MRM61C sequence but were similar to those in the published sequences of other flaviviruses from the same subgroup. pAKUN clone had four additional codon changes, i.e., Ile59 to Asn and Arg175 to Lys in NS2A and Tyr518 to His and Ser557 to Pro in NS3. Three of these substitutions except the previously shown marker mutation, Arg175 to Lys in NS2A, reverted to the wild-type sequence in the virus eventually recovered from pAKUN RNA-transfected BHK cells, demonstrating the functional importance of these residues in viral replication and/or viral assembly. Exchange of corresponding DNA fragments between pAKUN and FLSDX clones and site-directed mutagenesis revealed that the Tyr518-to-His mutation in NS3 was responsible for an similar to5-fold decrease in specific infectivity of transcribed RNA, while the Ile59-to-Asn mutation in NS2A completely blocked virus production. Correction of the Asn59 in pAKUN NS2A to the wild-type lie residue resulted in complete restoration of RNA infectivity. Replication of KUN replicon RNA with an Ile59-to-Asn substitution in NS2A and with a Ser557-to-Pro substitution in NS3 was not affected, while the Tyr518-to-His substitution in NS3 led to severe inhibition of RNA replication. The impaired function of the mutated NS2A in production of infectious virus was complemented in trans by the helper wild-type NS2A produced from the KUN replicon RNA. However, replicon RNA with mutated NS2A could not be packaged in trans by the KUN structural proteins. The data demonstrated essential roles for the KUN nonstructural protein NS2A in virus assembly and for NS3 in RNA replication and identified specific single-amino-acid residues involved in these functions.

Significance in replication of the terminal nucleotides of the Flavivirus genome

Point mutations that resulted in a substitution of the conserved 3'-penultimate cytidine in genomic RNA or the RNA negative strand of the self-amplifying replicon of the Flavivirus Kunjin virus completely blocked in vivo replication. Similarly, substitutions of the conserved 3'-terminal uridine in the RNA negative or positive strand completely blocked replication or caused much-reduced replication, respectively. The same preference for cytidine in the 3'-terminal dinucleotide was noted in reports of the in vitro activity of the RNA-dependent RNA polymerase (RdRp) for the other genera of Flaviviridae that also employ a double-stranded RNA (dsRNA) template to initiate asymmetric semiconservative RNA positive-strand synthesis. The Kunjin virus replicon results were interpreted in the context of a proposed model for initiation of RNA synthesis based on the solved crystal structure of the RdRp of phi6 bacteriophage, which also replicates efficiently using a dsRNA template with conserved 3'-penultimate cytidines and a 3'-terminal pyrimidine. A previously untested substitution of the conserved pentanucleotide at the top of the 3'-terminal stem-loop of all Flavivirus species also blocked detectable in vivo replication of the Kunjin virus replicon RNA.

The value of idiosyncratic markers and changes to conserved tRNA sequences from the mitochondrial genome of hard ticks (Acari: Ixodida: Ixodidae) for phylogenetic inference

Idiosyncratic markers are features of genes and genomes that are so unusual that it is unlikely that they evolved more than once in a lineage of organisms. Here we explore further the potential of idiosyncratic markers and changes to typically conserved tRNA sequences for phylogenetic inference. Hard ticks were chosen as the model group because their phylogeny has been studied extensively. Fifty-eight candidate markers from hard ticks ( family Ixodidae) and 22 markers from the subfamily Rhipicephalinae sensu lato were mapped onto phylogenies of these groups. Two of the most interesting markers, features of the secondary structure of two different tRNAs, gave strong support to the hypothesis that species of the Prostriata ( Ixodes spp.) are monophyletic. Previous analyses of genes and morphology did not strongly support this relationship, instead suggesting that the Prostriata is paraphyletic with respect to the Metastriata ( the rest of the hard ticks). Parallel or convergent evolution was not found in the arrangements of mitochondrial genes in ticks nor were there any reversals to the ancestral arthropod character state. Many of the markers identified were phylogenetically informative, whereas others should be informative with study of additional taxa. Idiosyncratic markers and changes to typically conserved nucleotides in tRNAs that are phylogenetically informative were common in this data set, and thus these types of markers might be found in other organisms.

The highly rearranged mitochondrial genome of the plague thrips, Thrips imaginis (Insecta : thysanoptera): Convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes

To help understand the mechanisms of gene rearrangement in the mitochondrial (mt) genomes of hemipteroid insects, we sequenced the mt genome of the plague thrips, Thrips imaginis (Thysanoptera). This genome is circular, 15,407 by long, and has many unusual features, including (1) rRNA genes inverted and distant from one another, (2) an extra gene for tRNA-Ser, (3) a tRNA-Val lacking a D-arm, (4) two pseudo-tRNA genes, (5) duplicate control regions, and (6) translocations and/or inversions of 24 of the 37 genes. The mechanism of rRNA gene transcription in T. imaginis may be different from that of other arthropods since the two rRNA genes have inverted and are distant from one another. Further, the rRNA genes are not adjacent or even close to either of the two control regions. Tandem duplication and deletion is a plausible model for the evolution of duplicate control regions and for the gene translocations, but intramitochondrial recombination may account for the gene inversions in T. imaginis. All the 18 genes between control regions #1 and #2 have translocated and/or inverted, whereas only six of the 20 genes outside this region have translocated and/or inverted. Moreover, the extra tRNA gene and the two pseudo-tRNA genes are either in this region or immediately adjacent to one of the control regions. These observations suggest that tandem duplication and deletion may be facilitated by the duplicate control regions and may have occurred a number of times in the lineage leading to T. imaginis. T. imaginis shares two novel gene boundaries with a lepidopsocid species from another order of hemipteroid insects, the Psocoptera. The evidence available suggests that these shared gene boundaries evolved by convergence and thus are not informative for the interordinal phylogeny of hemipteroid insects. We discuss the potential of hemipteroid insects as a model system for studies of the evolution of animal rut genomes and outline some fundamental questions that may be addressed with this system.

Prospects for whole genome linkage disequilibrium mapping in domestic dog breeds

Linkage disequilibrium (LD) mapping is commonly used as a fine mapping tool in human genome mapping and has been used with some success for initial disease gene isolation in certain isolated inbred human populations. An understanding of the population history of domestic dog breeds suggests that LID mapping could be routinely utilized in this species for initial genome-wide scans. Such an approach offers significant advantages over traditional linkage analysis. Here, we demonstrate, using canine copper toxicosis in the Bedlington terrier as the model, that LID mapping could be reasonably expected to be a useful strategy in low-resolution, genome-wide scans in pure-bred dogs. Significant LID was demonstrated over distances up to 33.3 cM. It is very unlikely, for a number of reasons discussed, that this result could be extrapolated to the rest of the genome. It is, however, consistent with the expectation given the population structure of canine breeds and, in this breed at least, with the hypothesis that it may be possible to utilize LID in a genome-wide scan. In this study, LD mapping confirmed the location of the copper toxicosis in Bedlington terrier gene (CT-BT) and was able to do so in a population that was refractory to traditional linkage analysis.