Curious bivalves: Systematic utility and unusual properties of anomalodesmatan mitochondrial genomes.

Mitogenomic trees for Bivalvia have proved problematic in the past, but several highly divergent lineages were missing from these analyses and increased representation of these groups may yet improve resolution. Here, we add seven new sequences from the Anomalodesmata and one unidentified semelid species (Bryopa lata, Euciroa cf. queenslandica, Laternula elliptica, Laternula truncata, Lyonsia norwegica, Myadora brevis, Tropidomya abbreviata, "Abra" sp.). We show that relationships in a mitogenomic tree for the Class are improved by the addition of seven anomalodesmatans from this highly divergent clade, but are still not completely consistent with relationships recovered in studies of nuclear genes. We suggest that some anomalous relationships (for instance the non-monophyly of Bivalvia) may be partially explained by compositional heterogeneity in the mitogenome and suggest that the addition of more taxa may help resolve both this effect and possible instances of long branch attraction. We also identify several curious features about anomalodesmatan mitogenomes. For example, many protein-coding gene boundaries are poorly defined in marine bivalves, but particularly so in anomalodesmatans, primarily due to non-conserved boundary sequences. The use of transcriptomic and genomic data together enabled better definition of gene boundaries, the identification of possible pseudogenes and suggests that most genes are translated monocistronically, which contrasts with many other studies. We also identified a possible case of gene duplication of ND5 in Myadora brevis (Myochamidae). Mitogenome size in the Anomalodesmata ranges from very small compact molecules, with the smallest for Laternula elliptica (Laternulidae) only 14,622bp, to Bryopa lata (Clavagellidae) which is at least 31,969bp long and may be >40,000bp. Finally, sampled species show a high degree of sequence divergence and variable gene order, although intraspecific variation in Laternula elliptica is very low.

Mitochondrial Genomes of Anopheles (Kerteszia) (Diptera: Culicidae) From the Atlantic Forest, Brazil.

Mitochondrial genome sequences are widely used as molecular markers for phylogenetic studies of mosquito species complexes, such as the Anopheles albitarsis complex. Except for a few studies that employed a limited number of nuclear or mitochondrial loci to address the genetic structure and species status of Anopheles cruzii, Anopheles bellator, and Anopheles homunculus, little is known about genetic markers that can be employed in studies focusing on Kerteszia species. The complete mitochondrial genomes of seven specimens of An. bellator, An. cruzii, An. homunculus, and Anopheles laneanus were sequenced using long-range polymerase chain reaction and Illumina sequencing. The mitochondrial genomes varied from 15,446 to 15,738 bp in length and contained 37 genes (13 protein-encoding genes, 2 rRNA genes [12S rRNA and 16S rRNA] and 22 tRNA genes), and the AT-rich control region, as all do other Anopheles mitochondrial genomes sequenced to date. Specimens from four populations of An. cruzii showed differences in codon composition.

The complete mitochondrial genome of a turbinid vetigastropod from MiSeq Illumina sequencing of genomic DNA and steps towards a resolved gastropod phylogeny.

A need to increase sampling of mitochondrial genomes for Vetigastropoda has been identified as an important step towards resolving relationships within the Gastropoda. We used shotgun sequencing of genomic DNA, using an Illumina MiSeq, to obtain the first mitochondrial genome for the vetigastropod family Turbinidae, doubling the number of genomes for the species-rich superfamily Trochoidea. This method avoids the necessity of finding suitable primers for long PCRs or primer-walking amplicons, resulting in a timely and cost-effective method for obtaining whole mitochondrial genomes from ethanol-preserved tissue samples. Bayesian analysis of amino acid variation for all available gastropod genomes including the new turbinid mtgenome produced a well resolved tree with high nodal support for most nodes. Major clades within Gastropoda were recovered with strong support, with the exception of Littorinimorpha, which was polyphyletic. We confirm here that mitogenomics is a useful tool for molluscan phylogenetics, especially when using powerful new models of amino acid evolution, but recognise that increased taxon sampling is still required to resolve existing differences between nuclear and mitochondrial gene trees.

Detection of a new mumps virus genotype during parotitis epidemic of 2006-2007 in the state of São Paulo, Brazil.

Nucleotide sequence analyses of the SH gene of 18 mumps virus isolates collected in the 2006-2007 parotitis epidemic in the state of São Paulo identified a new genotype, designated genotype M. This new designation fulfills all the parameters required to define a new mumps virus genotype. The parameters were established by an expert panel in collaboration with the World Health Organization (WHO) in 2005. This information will enhance the mumps virus surveillance program both at the national and global levels.

Structure and content of the Entamoeba histolytica genome.

The intestinal parasite Entamoeba histolytica is one of the first protists for which a draft genome sequence has been published. Although the genome is still incomplete, it is unlikely that many genes are missing from the list of those already identified. In this chapter we summarise the features of the genome as they are currently understood and provide previously unpublished analyses of many of the genes.

High-throughput structural biology in drug discovery: protein kinases.

Structural biology is an invaluable tool in modern drug discovery, providing key insights into the interactions of small-molecule drugs with their protein targets. As in many aspects of the drug discovery process, significant synergies can be realized in structural biology by the contemporaneous pursuit of many target proteins from a single structural and functional class. We will review some of those synergies here using the example of the protein kinases--an important class of drug targets that has recently been the subject of intensive study. We conclude by discussing some of the technical advances in X-ray crystallography that have enabled implementation of high-throughput structural biology as applied to drug lead optimization.

Comparison of EST libraries from seven beetle species: towards a framework for phylogenomics of the Coleoptera.

Relatively little is known about Coleoptera genes and genomes and how these compare in different taxa. We describe here the construction, DNA sequencing and sequence comparisons of cDNA libraries from seven beetle species. A total of 6717 bacterial colonies were screened for cDNA insert containing plasmids and 2784 size selected clones were 5'- and 3'-end sequenced to produce 1620 assembled sequences. Similarity comparisons with existing protein sequence databases revealed that 65.1% had matches (E < 10(-4)) in other organisms, with greater numbers of matches in Drosophila melanogaster than Caenorhabditis elegans and Saccharomyces cerevisiae databases. tBlastX comparisons also revealed numerous similarity hits (E < 10(-20)) in intra- and interlibrary comparisons. These results show the potential of small cDNA libraries for discovery and comparative analysis of genes useful for phylogenomic and functional studies.

Crystal structure of the macrocycle-forming thioesterase domain of the erythromycin polyketide synthase: versatility from a unique substrate channel.

As the first structural elucidation of a modular polyketide synthase (PKS) domain, the crystal structure of the macrocycle-forming thioesterase (TE) domain from the 6-deoxyerythronolide B synthase (DEBS) was solved by a combination of multiple isomorphous replacement and multiwavelength anomalous dispersion and refined to an R factor of 24.1% to 2.8-A resolution. Its overall tertiary architecture belongs to the alpha/beta-hydrolase family, with two unusual features unprecedented in this family: a hydrophobic leucine-rich dimer interface and a substrate channel that passes through the entire protein. The active site triad, comprised of Asp-169, His-259, and Ser-142, is located in the middle of the substrate channel, suggesting the passage of the substrate through the protein. Modeling indicates that the active site can accommodate and orient the 6-deoxyerythronolide B precursor uniquely, while at the same time shielding the active site from external water and catalyzing cyclization by macrolactone formation. The geometry and organization of functional groups explain the observed substrate specificity of this TE and offer strategies for engineering macrocycle biosynthesis. Docking of a homology model of the upstream acyl carrier protein (ACP6) against the TE suggests that the 2-fold axis of the TE dimer may also be the axis of symmetry that determines the arrangement of domains in the entire DEBS. Sequence conservation suggests that all TEs from modular polyketide synthases have a similar fold, dimer 2-fold axis, and substrate channel geometry.

Iron hydrogenases and the evolution of anaerobic eukaryotes.

Hydrogenases, oxygen-sensitive enzymes that can make hydrogen gas, are key to the function of hydrogen-producing organelles (hydrogenosomes), which occur in anaerobic protozoa scattered throughout the eukaryotic tree. Hydrogenases also play a central role in the hydrogen and syntrophic hypotheses for eukaryogenesis. Here, we show that sequences related to iron-only hydrogenases ([Fe] hydrogenases) are more widely distributed among eukaryotes than reports of hydrogen production have suggested. Genes encoding small proteins which contain conserved structural features unique to [Fe] hydrogenases were identified on all well-surveyed aerobic eukaryote genomes. Longer sequences encoding [Fe] hydrogenases also occur in the anaerobic eukaryotes Entamoeba histolytica and Spironucleus barkhanus, both of which lack hydrogenosomes. We also identified a new [Fe] hydrogenase sequence from Trichomonas vaginalis, bringing the total of [Fe] hydrogenases reported for this organism to three, all of which may function within its hydrogenosomes. Phylogenetic analysis and hypothesis testing using likelihood ratio tests and parametric bootstrapping suggest that the [Fe] hydrogenases in anaerobic eukaryotes are not monophyletic. Iron-only hydrogenases from Entamoeba, Spironucleus, and Trichomonas are plausibly monophyletic, consistent with the hypothesis that a gene for [Fe] hydrogenase was already present on the genome of the common, perhaps also anaerobic, ancestor of these phylogenetically distinct eukaryotes. Trees where the [Fe] hydrogenase from the hydrogenosomal ciliate Nyctotherus was constrained to be monophyletic with the other eukaryote sequences were rejected using a likelihood ratio test of monophyly. In most analyses, the Nyctotherus sequence formed a sister group with a [Fe] hydrogenase on the genome of the eubacterium Desulfovibrio vulgaris. Thus, it is possible that Nyctotherus obtained its hydrogenosomal [Fe] hydrogenase from a different source from Trichomonas for its hydrogenosomes. We find no support for the hypothesis that components of the Nyctotherus [Fe] hydrogenase fusion protein derive from the mitochondrial respiratory chain.

The structural basis for tRNA recognition and pseudouridine formation by pseudouridine synthase I.

Pseudouridine synthases catalyze the isomerization of specific uridines to pseudouridine in a variety of RNAs, yet the basis for recognition of the RNA sites or how they catalyze this reaction is unknown. The crystal structure of pseudouridine synthase I from Escherichia coli, which, for example, modifies positions 38, 39 and/or 40 in tRNA, reveals a dimeric protein that contains two positively charged, RNA-binding clefts along the surface of the protein. Each cleft contains a highly conserved aspartic acid located at its center. The structural domains have a topological similarity to those of other RNA-binding proteins, though the mode of interaction with tRNA appears to be unique. The structure suggests that a dimeric enzyme is required for binding transfer RNA and subsequent pseudouridine formation.

Compositional bias may affect both DNA-based and protein-based phylogenetic reconstructions.

It is now well-established that compositional bias in DNA sequences can adversely affect phylogenetic analysis based on those sequences. Phylogenetic analyses based on protein sequences are generally considered to be more reliable than those derived from the corresponding DNA sequences because it is believed that the use of encoded protein sequences circumvents the problems caused by nucleotide compositional biases in the DNA sequences. There exists, however, a correlation between AT/GC bias at the nucleotide level and content of AT- and GC-rich codons and their corresponding amino acids. Consequently, protein sequences can also be affected secondarily by nucleotide compositional bias. Here, we report that DNA bias not only may affect phylogenetic analysis based on DNA sequences, but also drives a protein bias which may affect analyses based on protein sequences. We present a striking example where common phylogenetic tools fail to recover the correct tree from complete animal mitochondrial protein-coding sequences. The data set is very extensive, containing several thousand sites per sequence, and the incorrect phylogenetic trees are statistically very well supported. Additionally, neither the use of the LogDet/paralinear transform nor removal of positions in the protein alignment with AT- or GC-rich codons allowed recovery of the correct tree. Two taxa with a large compositional bias continually group together in these analyses, despite a lack of close biological relatedness. We conclude that even protein-based phylogenetic trees may be misleading, and we advise caution in phylogenetic reconstruction using protein sequences, especially those that are compositionally biased.

Nucleotide composition bias affects amino acid content in proteins coded by animal mitochondria.

We show that in animal mitochondria homologous genes that differ in guanine plus cytosine (G + C) content code for proteins differing in amino acid content in a manner that relates to the G + C content of the codons. DNA sequences were analyzed using square plots, a new method that combines graphical visualization and statistical analysis of compositional differences in both DNA and protein. Square plots divide codons into four groups based on first and second position A + T (adenine plus thymine) and G + C content and indicate differences in amino acid content when comparing sequences that differ in G + C content. When sequences are compared using these plots, the amino acid content is shown to correlate with the nucleotide bias of the genes. This amino acid effect is shown in all protein-coding genes in the mitochondrial genome, including cox I, cox II, and cyt b, mitochondrial genes which are commonly used for phylogenetic studies. Furthermore, nucleotide content differences are shown to affect the content of all amino acids with A + T- and G + C-rich codons. We speculate that phylogenetic analysis of genes so affected may tend erroneously to indicate relatedness (or lack thereof) based only on amino acid content.

Unbiased estimation of symmetrical directional mutation pressure from protein-coding DNA.

The most generally applicable procedure for obtaining estimates of the symmetrical, or strandnonspecific, directional mutation pressure (microD) on protein-coding DNA sequences is to determine the G+C content at synonymous codon sites (Psyn), and to divide Psyn by twice the arithmetic mean of the G+C content at synonymous codon sites of a large number of randomly generated, synonymously coding DNA sequences (Psyn). Unfortunately, the original procedure yields biased estimates of Psyn and microD and is computationally expensive. We here present a fast procedure for estimating unbiased microD values. The procedure employs direct calculation of Psyn (approximately Psyn) and two normalization procedures, one for Psyn < or = Psyn and another for Psyn > or = Psyn. The normalization removes a bias sometimes caused by codons specifying arginine, asparagine, isoleucine, and leucine. Consequently, comparison of protein-coding genes that are translated using different genetic codes is facilitated.

Mechanism-based inhibition of thymidylate synthase by 5-(trifluoromethyl)-2'-deoxyuridine 5'-monophosphate.

Thymidylate synthase (TS) from Lactobacillus casei is inhibited by 5-(trifluoromethyl)-2'-deoxyuridine 5'-monophosphate (CF3dUMP). CF3dUMP binds to the active site of TS in the absence of 5,10-methylenetetrahydrofolate, and attack of the catalytic nucleophile cysteine 198 at C6 of the pyrimidine leads to activation of the trifluoromethyl group and release of fluoride ion. Subsequently, the activated heterocycle reacts with a nucleophile of the enzyme to form a moderately stable covalent complex. Proteolytic digestion of TS treated with [2'-3H]CF3dUMP, followed by sequencing of the labeled peptides, revealed that tyrosine 146 and cysteine 198 are covalently bound to the inhibitor in the enzyme-inhibitor complex. The presence of dithiothreitol (DTT) or beta-mercaptoethanol resulted in the breakdown of the covalent complex, and products from the breakdown of the complex were isolated and characterized. The three-dimensional structure of the enzyme-inhibitor complex was determined by X-ray crystallography, clearly demonstrating covalent attachment of the nucleotide to tyrosine 146. A chemical reaction mechanism for the inhibition of TS by CF3dUMP is presented that is consistent with the kinetic, biochemical, and structural results.

Refined structures of substrate-bound and phosphate-bound thymidylate synthase from Lactobacillus casei.

Crystal structures of two crystal forms of the complex of Lactobacillus casei (TS) with its substrate dUMP have been solved and refined at 2.55 A resolution. The two crystal forms differ by approximately 5% in the c-axis length. The TS-dUMP complexes are symmetric dimers with dUMP bound equivalently in both active sites. dUMP is non-covalently bound in the same conformation as in ternary complexes of TS with dUMP and cofactor or cofactor analogs. The same hydrogen bonds are made between TS and substrate in the binary and ternary complexes. We have also determined the 2.36 A crystal structure of phosphate-bound L. casei TS. This structure has been refined to an R-factor of 19.3% with highly constrained geometry. Refinement has revealed the locations of all residues in the protein, including the disordered residues 90 to 119, which are part of an insert found only in the L. casei and Staphylococcus aureus transposon Tn4003 TS sequences. The 2.9 A multiple isomorphous replacement (MIR) structure of L. casei TS in a complex with its substrate dUMP has been refined to a crystallographic R-factor of 15.5%. Reducing agents were withheld from crystallization solutions during MIR structure determination to allow heavy-metal labeling of the cysteine residues. Therefore, the active-site cysteine residue in this structure is oxidized and the dUMP is found at half-occupancy in the active site. No significant conformational difference was found between the phosphate-bound and dUMP-bound structures. The TS-dUMP structures were better ordered than the phosphate-bound TS or the oxidized TS-dUMP, particularly Arg23, which is clearly hydrogen-bonded to the phosphate group of dUMP. A large and a small P6(1)22 crystal form are observed for both phosphate-bound and dUMP-bound L. casei TS. The small cell forms of the phosphate-bound and dUMP-bound enzyme are isomorphous, whereas the cell constants of the larger cell form change slightly when dUMP is bound (c = 240 A versus c = 243 A). For both liganded and unliganded enzyme, conversion from the small to the large crystal form sometimes occurs spontaneously, and the crystal packing changes at a single interface. Conversion may be the result of a small change in pH in the mother liquor surrounding the crystal. A single intermolecular contact between symmetry-related Asp287 residues is disrupted on going from the small to the large crystal form.

Construction of a dihydrofolate reductase-deficient mutant of Escherichia coli by gene replacement.

The dihydrofolate reductase (fol) gene in Escherichia coli has been deleted and replaced by a selectable marker. Verification of the delta fol::kan strain has been accomplished using genetic and biochemical criteria, including Southern analysis of the chromosomal DNA. The delta fol::kan mutation is stable in E. coli K549 [thyA polA12 (Ts)] and can be successfully transduced to other E. coli strains providing they have mutations in their thymidylate synthetase (thyA) genes. A preliminary investigation of the relationship between fol and thyA gene expression suggests that a Fol- cell (i.e., a dihydrofolate reductase deficiency phenotype) is not viable unless thymidylate synthetase activity is concurrently eliminated. This observation indicates that either the nonproductive accumulation of dihydrofolate or the depletion of tetrahydrofolate cofactor pools is lethal in a Fol- ThyA+ strain. Strains containing the thyA delta fol::kan lesions require the presence of Fol end products for growth, and these lesions typically increase the doubling time of the strain by a factor of 2.5 in rich medium.

Immunological relationships among vertebrate lysine-rich histones.

1. The relationship between sequence homology and immunological cross-reaction was investigated by enzyme-linked immunosorbent assay and immunoblotting using polyclonal antisera against lysine-rich histones (LRH) of known sequence, chicken H1 and H5, trout Hl and Xenopus H1s. 2. The order of immunological relatedness was consistent with known homologies among these LRH and goose H5, but quantitative correlations reflected varied locations of antigenic determinants. 3. When immunoblotting was extended to LRH from eight more vertebrates, it was evident that avian H5, mammalian H1o and anuran H1s form a sub-class, to which turtle H1s may belong, that urodele erythrocytes contain no H5-like histone and that fish "H5" is more like H1 than the H5/H1s/H1o subclass.

Predictability of sequence homologies among lysine-rich histones by immunological distance.

The relationship between immunological distance (I.D.) measured by microcomplement fixation and amino acid sequence difference for lysine-rich histones was tested using antisera to lysine-rich histones of known sequence, chicken H1 and H5, goose H5, and trout H1 as well as to trout H5. The best relationship between I.D. (y) and percent sequence difference (x) for lysine-rich histones, y = 2x, applies as well to other histones of known sequence but it differs from y = 5x, reported for other proteins and often used for histones. Although deviations indicate that I.D. is a poor predictor of primary sequence differences among histones, it suggests that trout H5 is more closely related to H1 than to chicken H5.