ECOdrug: a database connecting drugs and conservation of their targets across species.

Pharmaceuticals are designed to interact with specific molecular targets in humans and these targets generally have orthologs in other species. This provides opportunities for the drug discovery community to use alternative model species for drug development. It also means, however, there is potential for mode of action related effects in non-target wildlife species as many pharmaceuticals reach the environment through patient use and manufacturing wastes. Acquiring insight in drug target ortholog predictions across species and taxonomic groups has proven difficult because of the lack of an optimal strategy and because necessary information is spread across multiple and diverse sources and platforms. We introduce a new research platform tool, ECOdrug, that reliably connects drugs to their protein targets across divergent species. It harmonizes ortholog predictions from multiple sources via a simple user interface underpinning critical applications for a wide range of studies in pharmacology, ecotoxicology and comparative evolutionary biology. ECOdrug can be used to identify species with drug targets and identify drugs that interact with those targets. As such, it can be applied to support intelligent targeted drug safety testing by ensuring appropriate and relevant species are selected in ecological risk assessments. ECOdrug is freely accessible and available at: http://www.ecodrug.org.

The novel isoxazoline ectoparasiticide lotilaner (Credelio™): a non-competitive antagonist specific to invertebrates γ-aminobutyric acid-gated chloride channels (GABACls).

BACKGROUND: The isoxazolines are a novel class of parasiticides that are potent inhibitors of γ-aminobutyric acid (GABA)-gated chloride channels (GABACls) and, to a lesser extent, of inhibitory glutamate-gated chloride channels (GluCls). Lotilaner (Credelio™), a novel representative of this chemical class, is currently evaluated for its excellent ectoparasiticide properties. METHODS: In this study, we investigated the molecular mode of action and pharmacology of lotilaner. We report the successful gene identification, cDNA cloning and functional expression in Xenopus oocytes of Drosohpila melanogaster (wild type and dieldrin/fipronil-resistant forms), Lepeophtheirus salmonis (an ectoparasite copepod crustacean of salmon), Rhipicephalus microplus and Canis lupus familiaris GABACls. Automated Xenopus oocyte two-electrode voltage clamp electrophysiology was used to assess GABACls functionality and to compare ion channel inhibition by lotilaner with that of established insecticides addressing GABACls as targets. RESULTS: In these assays, we demonstrated that lotilaner is a potent non-competitive antagonist of insects (fly) GABACls. No cross-resistance with dieldrin or fipronil resistance mutations was detected, suggesting that lotilaner might bind to a site at least partly different from the one bound by known GABACl blockers. Using co-application experiments, we observed that lotilaner antagonism differs significantly from the classical open channel blocker fipronil. We finally confirmed for the first time that isoxazoline compounds are not only powerful antagonists of GABACls of acari (ticks) but also of crustaceans (sea lice), while no activity on a dog GABAA receptor was observed up to a concentration of 10 µM. CONCLUSIONS: Together, these results demonstrate that lotilaner is a non-competitive antagonist specific to invertebrate's γ-aminobutyric acid-gated chloride channels (GABACls). They contribute to our understanding of the mode of action of this new ectoparasiticide compound.

Genetic assessment of meiobenthic community composition and spatial distribution in coastal sediments along northern Gulf of Mexico.

Meiobenthic (meiofauna and micro-eukaryotes) organisms are important contributors to ecosystem functioning in aquatic environments through their roles in nutrient transport, sediment stability, and food web interactions. Despite their ecological importance, information pertaining to variation of these communities at various spatial and temporal scales is not widely known. Many studies in the Gulf of Mexico (GOM) have focused either on deep sea or continental shelf areas, while little attention has been paid to bays and coastal regions. Herein, we take a holistic approach by using high-throughput sequencing approaches to examine spatial variation in meiobenthic communities within Alabama bays and the coastal northern GOM region. Sediment samples were collected along three transects (Mississippi Sound: MS, FOCAL: FT, and Orange Beach: OB) from September 2010 to April 2012 and community composition was determined by metabarcoding the V9 hypervariable region of the nuclear18S rRNA gene. Results showed that Stramenopiles (diatoms), annelids, arthropods (copepods), and nematodes were the dominate groups within samples, while there was presence of other phyla throughout the dataset. Location played a larger role than time sampled in community composition. However, samples were collected over a short temporal scale. Samples clustered in reference to transect, with the most eastern transect (OB) having a distinct community composition in comparison to the other two transects (MS and FT). Communities also differed in reference to region (Bay versus Shelf). Bulk density and percent inorganic carbon were the only measured environmental factors that were correlated with community composition.

Ernest Everett Just: Egg and embryo as excitable systems.

Ernest Everett Just (1883-1941) was an African American embryologist of international standing whose research interests lay in the area of fertilization and early development in marine invertebrates. Perhaps best known for his discovery of the dynamical and structural blocks to polyspermy that sweep over the egg upon fertilization, E. E. Just also was the first to associate cell surface changes with stages of embryonic development. He was deeply familiar with the natural history of the animals whose eggs he studied, and his knowledge of natural settings led him to emphasize the importance of using laboratory conditions that closely match those in nature. Based on more than 30 years of work, he came to believe that it was the cell surface that played the most critical role in development, heredity, and evolution. He promoted a holistic view of cells and organisms in opposition to the gene-centric view that was becoming more prevalent with the rise of genetics, but rejected the vitalism espoused by some biologists of his era, calling instead for "a physics and chemistry in a new dimension …superimposed upon the now known physics and chemistry" to account for biological phenomena. Just's incisive critique of genetic reductionism finds echoes in contemporary multiscale, systems approaches in biology. His speculations on the relationship between developmental and evolutionary mechanisms resonate with today's evolutionary developmental biology. After a brief biographical sketch, this paper outlines and discusses some of Just's scientific contributions, and shows how his ideas remain relevant today.

Dicyema Pax6 and Zic: tool-kit genes in a highly simplified bilaterian.

BACKGROUND: Dicyemid mesozoans (Phylum Dicyemida) are simple (8-40-cell) cephalopod endoparasites. They have neither body cavities nor differentiated organs, such as nervous and gastrointestinal systems. Whether dicyemids are intermediate between Protozoa and Metazoa (as represented by their "Mesozoa" classification) or degenerate species of more complex metazoans is controversial. Recent molecular phylogenetic studies suggested that they are simplified bilaterians belonging to the Lophotrochozoa. We cloned two genes developmentally critical in bilaterian animals (Pax6 and Zic), together with housekeeping genes (actin, fructose-bisphosphate aldolase, and ATP synthase beta subunit) from a dicyemid to reveal whether their molecular phylogeny supported the "simplification" hypothesis, and to clarify evolutionary changes in dicyemid gene structure and expression profiles. RESULTS: Genomic/cDNA sequence analysis showed that 1) the Pax6 molecular phylogeny and Zic intron positions supported the idea of dicyemids as reduced bilaterians; 2) the aa sequences deduced from the five genes were highly divergent; and 3) Dicyema genes contained very short introns of uniform length. In situ hybridization analyses revealed that Zic genes were expressed in hermaphroditic gonads, and Pax6 was expressed weakly throughout the developmental stages of the 2 types of embryo and in the hermaphroditic gonads. CONCLUSION: The accelerated evolutionary rates and very short and uniform intron may represent a part of Dicyema genomic features. The presence and expression of the two tool-kit genes (Pax6 and Zic) in Dicyema suggests that they can be very versatile genes even required for the highly reduced bilaterian like Dicyema. Dicyemids may be useful models of evolutionary body plan simplification.

Ancient biomolecules from deep ice cores reveal a forested southern Greenland.

It is difficult to obtain fossil data from the 10% of Earth's terrestrial surface that is covered by thick glaciers and ice sheets, and hence, knowledge of the paleoenvironments of these regions has remained limited. We show that DNA and amino acids from buried organisms can be recovered from the basal sections of deep ice cores, enabling reconstructions of past flora and fauna. We show that high-altitude southern Greenland, currently lying below more than 2 kilometers of ice, was inhabited by a diverse array of conifer trees and insects within the past million years. The results provide direct evidence in support of a forested southern Greenland and suggest that many deep ice cores may contain genetic records of paleoenvironments in their basal sections.

Retroposition and evolution of the DNA-binding motifs of YY1, YY2 and REX1.

YY1 is a DNA-binding transcription factor found in both vertebrates and invertebrates. Database searches identified 62 YY1 related sequences from all the available genome sequences ranging from flying insects to human. These sequences are characterized by high levels of sequence conservation, ranging from 66% to 100% similarity, in the zinc finger DNA-binding domain of the predicted proteins. Phylogenetic analyses uncovered duplication events of YY1 in several different lineages, including flies, fish and mammals. Retroposition is responsible for generating one duplicate in flies, PHOL from PHO, and two duplicates in placental mammals, YY2 and Reduced Expression 1 (REX1) from YY1. DNA-binding motif studies have demonstrated that YY2 still binds to the same consensus sequence as YY1 but with much lower affinity. In contrast, REX1 binds to DNA motifs divergent from YY1, but the binding motifs of REX1 and YY1 share some similarity at their core regions (5'-CCAT-3'). This suggests that the two duplicates, YY2 and REX1, although generated through similar retroposition events have undergone different selection schemes to adapt to new roles in placental mammals. Overall, the conservation of YY2 and REX1 in all placental mammals predicts that each duplicate has co-evolved with some unique features of eutherian mammals.

Insights into the molecular evolution of the PDZ/LIM family and identification of a novel conserved protein motif.

The PDZ and LIM domain-containing protein family is encoded by a diverse group of genes whose phylogeny has currently not been analyzed. In mammals, ten genes are found that encode both a PDZ- and one or several LIM-domains. These genes are: ALP, RIL, Elfin (CLP36), Mystique, Enigma (LMP-1), Enigma homologue (ENH), ZASP (Cypher, Oracle), LMO7 and the two LIM domain kinases (LIMK1 and LIMK2). As conventional alignment and phylogenetic procedures of full-length sequences fell short of elucidating the evolutionary history of these genes, we started to analyze the PDZ and LIM domain sequences themselves. Using information from most sequenced eukaryotic lineages, our phylogenetic analysis is based on full-length cDNA-, EST-derived- and genomic- PDZ and LIM domain sequences of over 25 species, ranging from yeast to humans. Plant and protozoan homologs were not found. Our phylogenetic analysis identifies a number of domain duplication and rearrangement events, and shows a single convergent event during evolution of the PDZ/LIM family. Further, we describe the separation of the ALP and Enigma subfamilies in lower vertebrates and identify a novel consensus motif, which we call 'ALP-like motif' (AM). This motif is highly-conserved between ALP subfamily proteins of diverse organisms. We used here a combinatorial approach to define the relation of the PDZ and LIM domain encoding genes and to reconstruct their phylogeny. This analysis allowed us to classify the PDZ/LIM family and to suggest a meaningful model for the molecular evolution of the diverse gene architectures found in this multi-domain family.

Hidden diversity and host specificity in cycliophorans: a phylogeographic analysis along the North Atlantic and Mediterranean Sea.

In order to elucidate the evolutionary history and the population structure of the members of the phylum Cycliophora, which live commensally on three species of lobsters, we studied sequence variation in the mitochondrial gene cyctochrome c oxidase subunit I. Overall 242 sequences from 16 locations on both coasts of the North Atlantic, including the North Sea and the Mediterranean, were analysed, revealing 28 haplotypes, with a maximum sequence divergence of 16.6%. Total genetic diversity was high (h = 0.8322, pi = 0.0898), as it was for the commensals on Homarus americanus (17 haplotypes, h = 0.7506, pi = 0.0504). However, it was low for commensals on Nephrops norvegicus (6 haplotypes, h = 0.3899, pi = 0.0035), and intermediate for cycliophorans on Homarus gammarus (5 haplotypes, h = 0.3020, pi = 0.0140). Although two of the host lobsters co-inhabit the coastal waters of Europe, a strong genetic structure (78.45% of the observed genetic variation) was detected among populations on all host species, indicating the existence of a reproductively isolated species on each lobster. In addition, genetic structure over long distances exists among populations on each host species. Such patterns can be explained by the limited dispersal ability of the cycliophoran chordoid larva. Demographic and phylogenetic analyses suggest old and possibly cryptic populations present on H. americanus and H. gammarus, while the latter may have experienced recent bottlenecks, perhaps during Pleistocene glaciations. Populations on N. norvegicus appear to be of recent origin.

DIRS-1 and the other tyrosine recombinase retrotransposons.

DIRS-1 is a retroelement from the slime mold Dictyostelium discoideum. Until recently only two related retrotransposons had been described: PAT from the nematode Panagrellus redivivus and Prt1 from the zygomycete fungus Phycomyces blakesleeanus. Analyses of the reverse transcriptase sequences encoded by these three elements suggested that they were closely related to each other and more distantly related to the Ty3/gypsy Long Terminal Repeat (LTR) retroelements. They have several unusual structural features that distinguish them from typical LTR elements. For instance, they each encode a tyrosine recombinase (YR), but not a DDE-type integrase or an aspartic protease. Although the DIRS-1-related elements are bordered by terminal repeats these differ from typical LTRs in a number of ways. In DIRS-1, for example, the terminal repeats are inverted (complementary), non-identical in sequence, and the outer edges of the terminal sequences are repeated (adjacent to each other) in the internal region. PAT has so-called "split" direct repeats in which the unrelated terminal sequences appear as direct repeats adjacent to each other in the internal region. The only repetition displayed by Prt1 is the presence of short inverted terminal repeats, but the sequenced copy of this element is believed to be a truncated version of an element with a structure resembling DIRS-1. The unusual structure of the terminal repeats of the DIRS1-like elements appears to be related to their replication via free circular intermediates. Site-specific recombination is believed to integrate the circle without creating duplications of the target sites. In recognition of these important distinctions it is proposed that the retrotransposons that encode tyrosine recombinases be called the tyrosine recombinase (or YR) retrotransposons. Recently a large number of additional YR retrotransposons have been described, including elements from fungi (zygomycetes and basidiomycetes), plants (green algae) and a wide range of animals including nematodes, insects, sea urchins, fish and amphibia, while remnants of elements related to DIRS-1 occur in the human genome. The complete set of YR retrotransposons can be divided into two major groups, the DIRS elements and the Ngaro elements, the two groups forming distinct clades on phylogenetic trees based on alignments of RT/RH and recombinase sequences, and also having some structural distinctions. A third group of transposable elements, which we call Cryptons, also carry tyrosine recombinases. These elements do not encode a reverse transcriptase and so are believed to be DNA transposons not retrotransposons. They have been detected in several pathogenic fungi, including the basidiomycete Cryptococcus neoformans, and the ascomycetes Coccidioides posadasii and Histoplasma capsulatum. Sequence comparisons suggest that the Crypton YRs are related to those of the YR retrotransposons. We suggest that the YR retrotransposons arose from the combination of a Crypton-like YR DNA transposon and the RT/RH encoding sequence of a retrotransposon.

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.

Tolerance to biodegraded crude oil in marine invertebrate embryos and larvae is associated with expression of a multixenobiotic resistance transporter.

The toxicity of water-soluble fractions of biodegraded crude oil (BWSF) to embryos and larvae of two marine invertebrates, the white sea urchin (Lytechinus anamesus) and the fat innkeeper (Urechis caupo), was studied. Santa Barbara Channel crude oil was artificially weathered and subjected to biodegradation using a mixed microbe culture obtained from natural oil seep sites. The degradation culture inoculated with seep sediment microbes accumulated 43.7 microg/l water-soluble hydrocarbons. In contrast water-soluble fractions from the non-degraded cultures (NWSF) only accumulated 3.05 microg/l. BWSF proved deleterious to Lytechinus embryo development at low concentrations (EC50 = 0.33 mg/l) but was essentially non-toxic to Urechis embryos/larvae up to 3.0 mg/l. An established mechanism for handling of a wide array of xenobiotics in Urechis embryos is the multixenobiotoic resistance transporter multixenobiotic response (MXR, also known as multidrug resistance, MDR). This mechanism is primarily mediated by ATP-dependent, efflux pumps that extrude a wide array of xenobiotic compounds. In this study, we show that Lytechinus larvae do not appear to express MXR efflux protein nor MXR mediated dye efflux capacity. In contrast, BWSF acts as a competitive inhibitor of MXR transport-mediated dye efflux in Urechis larvae. These results suggest that MXR may be an important mechanism for extrusion of the by-products of crude oil degradation by microbes, and that the level of its expression may determine the susceptibility of organisms to degraded oil hydrocarbons.

Pacemaker oscillations in heart and brain: a key role for hyperpolarization-activated cation channels.

Rhythmic activity of single cells or multicellular networks is a common feature of all organisms. The oscillatory activity is characterized by time intervals of several seconds up to many hours. Cellular rhythms govern the beating of the heart, the swimming behavior of sperm, cycles of sleep and wakefulness, breathing, and the release of hormones. Many neurons in the brain and cardiac cells are characterized by endogenous rhythmic activity, which relies on a complex interplay between several distinct ion channels. In particular, one type of ion channel plays a prominent role in the control of rhythmic electrical activity since it determines the frequency of the oscillations. The activity of the channels is thus setting the "pace" of the oscillations; therefore, these channels are often referred to as "pacemaker" channels. Despite their obvious important physiological function, it was not until recently that genes encoding pacemaker channels have been identified. Because both hyperpolarization and cyclic nucleotides are key elements that control their activity, pacemaker channels have now been designated hyperpolarization-activated and cyclic nucleotide-gated (HCN) channels. The molecular identification of the channels and the upcoming studies on their properties in heterologous systems will certainly enhance our understanding of "pacemaking" in physiological systems. This review gives a brief insight into the physiological importance of these channels and sums up what we have learned since the first cloning of genes succeeded (for recent reviews, see also Clapham 1998; Luthi and McCormick 1998a; Biel et al. 1999; Ludwig, Zong, Hofmann, et al. 1999; Santoro and Tibbs 1999).

A novel A-isoform-like inositol 1,4,5-trisphosphate 3-kinase from chicken erythrocytes exhibits alternative splicing and conservation of intron positions between vertebrates and invertebrates.

Based on the partial peptide sequence of inositol 1,4, 5-trisphosphate 3-kinase purified with 135 000-fold enrichment from chicken erythrocytes, cDNA-fragments were cloned by RT-PCR using degenerate oligonucleotides. Subsequent hybridization screening of an embryonic chicken cDNA library and 5'-RACE yielded a cDNA-contig of 2418 bp, encoding a 452 amino acid protein. The amino acid sequence shows the highest degree of homology with A-isoforms of inositol 1,4,5-trisphosphate 3-kinase (65% identities), whereas homology towards B and C isoforms was lower (57% and 52% amino acid identities respectively). These findings reveal a new tissue-specific pattern of A-isoform expression, a form which so far has only been found in brain and testes. Two overlapping lambda-genomic clones for chicken inositol 1,4,5-trisphosphate 3-kinase, isolated by hybridization screening, covered 18 499 bp of genomic sequence. This contig included four exons: three of them were present in all cDNA clones, whereas one was only represented in a single cDNA clone. In addition, the sequence of the latter differed from the other cDNAs by an in-frame deletion of 72 bp within the coding region for the catalytic domain of the enzyme. This divergent cDNA suggests the existence of alternative splice products, at least in embryonic tissue.A comparison of the position of introns, with the respective introns known from the corresponding gene from Caenorhabditis elegans, revealed a high degree of conservation of intron positions between vertebrates and invertebrates. Functional data for the enzyme suggests that the conserved exons represent defined functional protein modules.

Nuclear receptors, nuclear-receptor factors, and nuclear-receptor-like orphans form a large paralog cluster in Homo sapiens.

We studied a human protein paralog cluster formed by 38 nonredundant sequences taken from the Swiss-Prot database and its supplement, TrEMBL. These sequences include nuclear receptors, nuclear-receptor factors and nuclear-receptor-like orphans. Working separately with both the central cysteine-rich DNA-binding domain and the carboxy-terminal ligand-binding domain, we performed multialignment analyses that included drawings of paralog trees. Our results show that the cluster is highly multibranched, with considerable differences in the amino acid sequence in the ligand-binding domain (LBD), and 17 proximal subbranches which are identifiable and fully coincident when independent trees from both domains are compared. We identified the six recently proposed subfamilies as groups of neighboring clusters in the LBD paralog tree. We found similarities of 80%-100% for the N-terminal transactivation domain among mammalian ortholog receptors, as well as some paralog resemblances within diverse subbranches. Our studies suggest that during the evolutionary process, the three domains were assembled in a modular fashion with a nonshuffled modular fusion of the LBD. We used the EMBL server PredictProtein to make secondary-structure predictions for all 38 LBD subsequences. Amino acid residues in the multialigned homologous domains--taking the beginning of helix H3 of the human retinoic acid receptor-gamma as the initial point of reference--were substituted with H or E, which identify residues predicted to be helical or extended, respectively. The result was a secondary structure multialignment with the surprising feature that the prediction follows a canonical pattern of alignable alpha-helices with some short extended elements in between, despite the fact that a number of subsequences resemble each other by less than 25% in terms of the similarity index. We also identified the presence of a binary patterning in all of the predicted helices that were conserved throughout the 38-sequence sample. Our results fit well with a recently proposed evolutionary model that combines protein secondary structure and amino acid replacement. We propose a new hypothesis for molecular evolution, in which chaperones--acting as an endogenous cellular device for selection--play a crucial role in preserving protein secondary structure.

Organization of non-vertebrate globin genes.

The organization of non-vertebrate globin genes exhibits substantially more variability than the three-exon, two-intron structure of the vertebrate globin genes. (1) The structures of genes of the single-domain globin chains of the annelid Lumbricus and the mollusc Anadara, and the globin gene coding for the two-domain chains of the clam Barbatia, are similar to the vertebrate plan. (2) Genes for single-domain chains exist in bacteria and protozoa. Although the globin gene is highly expressed in the bacterium Vitreoscilla, the putative globin gene hmp in E. coli, which codes for a chimeric protein whose N-terminal moiety of 139 residues contains 67 residues identical to the Vitreoscilla globin, may be either unexpressed or expressed at very low levels, despite the presence of normal regulatory sequences. The DNA sequence of the globin gene of the protozoan Paramecium, determined recently by Yamauchi and collaborators, appears to consist of two exons separated by a short intron. (3) Among the lower eukaryotes, the yeasts Saccharomyces and Candida have chimeric proteins consisting of N-terminal globin and C-terminal flavoprotein moieties of about the same size. The structure of the gene for the chimeric protein of Saccharomyces exhibits no introns. According to Riggs, the presence of chimeric proteins in E. coli and other prokaryotes, such as Alcaligenes and Rhizobium, as well as in yeasts, suggests a previously unrecognized evolutionary pathway for hemoglobin, namely that of a multipurpose heme-binding domain attached to a variety of unrelated proteins with diverse functions. (4) The published globin gene sequences of the insect larva Chironomus have an intron-less structure and are present as clusters of multiple copies; the expression of the globin genes is tissue and developmental stage-specific. Furthermore, the expression of many of these genes has not yet been demonstrated despite the presence of apparently normal regulatory sequences in the two flanking regions. Unexpectedly, Bergtrom and collaborators have recently shown that at least three Ctt globin II beta genes contain putative introns. (5) Pohajdak and collaborators have found a seven-exon and six-intron structure for the globin gene of the nematode Pseudoterranova which codes for a two-domain globin chain. Although the second and fourth introns of the N-terminal domain correspond to the two introns found in vertebrate globin genes, the position of the third intron is close to that of the central intron in plant hemoglobins.(ABSTRACT TRUNCATED AT 400 WORDS)