Invertebrate Gonadotropin-Releasing Hormone Receptor Signaling and Its Relevant Biological Actions.

Gonadotropin-releasing hormones (GnRHs) play pivotal roles in reproduction via the hypothalamus-pituitary-gonad axis (HPG axis) in vertebrates. GnRHs and their receptors (GnRHRs) are also conserved in invertebrates lacking the HPG axis, indicating that invertebrate GnRHs do not serve as "gonadotropin-releasing factors" but, rather, function as neuropeptides that directly regulate target tissues. All vertebrate and urochordate GnRHs comprise 10 amino acids, whereas amphioxus, echinoderm, and protostome GnRH-like peptides are 11- or 12-residue peptides. Intracellular calcium mobilization is the major second messenger for GnRH signaling in cephalochordates, echinoderms, and protostomes, while urochordate GnRHRs also stimulate cAMP production pathways. Moreover, the ligand-specific modulation of signal transduction via heterodimerization between GnRHR paralogs indicates species-specific evolution in Ciona intestinalis. The characterization of authentic or putative invertebrate GnRHRs in various tissues and their in vitro and in vivo activities indicate that invertebrate GnRHs are responsible for the regulation of both reproductive and nonreproductive functions. In this review, we examine our current understanding of and perspectives on the primary sequences, tissue distribution of mRNA expression, signal transduction, and biological functions of invertebrate GnRHs and their receptors.

PathwAX: a web server for network crosstalk based pathway annotation.

Pathway annotation of gene lists is often used to functionally analyse biomolecular data such as gene expression in order to establish which processes are activated in a given experiment. Databases such as KEGG or GO represent collections of how genes are known to be organized in pathways, and the challenge is to compare a given gene list with the known pathways such that all true relations are identified. Most tools apply statistical measures to the gene overlap between the gene list and pathway. It is however problematic to avoid false negatives and false positives when only using the gene overlap. The pathwAX web server (http://pathwAX.sbc.su.se/) applies a different approach which is based on network crosstalk. It uses the comprehensive network FunCoup to analyse network crosstalk between a query gene list and KEGG pathways. PathwAX runs the BinoX algorithm, which employs Monte-Carlo sampling of randomized networks and estimates a binomial distribution, for estimating the statistical significance of the crosstalk. This results in substantially higher accuracy than gene overlap methods. The system was optimized for speed and allows interactive web usage. We illustrate the usage and output of pathwAX.

An expanded Notch-Delta model exhibiting long-range patterning and incorporating MicroRNA regulation.

Notch-Delta signaling is a fundamental cell-cell communication mechanism that governs the differentiation of many cell types. Most existing mathematical models of Notch-Delta signaling are based on a feedback loop between Notch and Delta leading to lateral inhibition of neighboring cells. These models result in a checkerboard spatial pattern whereby adjacent cells express opposing levels of Notch and Delta, leading to alternate cell fates. However, a growing body of biological evidence suggests that Notch-Delta signaling produces other patterns that are not checkerboard, and therefore a new model is needed. Here, we present an expanded Notch-Delta model that builds upon previous models, adding a local Notch activity gradient, which affects long-range patterning, and the activity of a regulatory microRNA. This model is motivated by our experiments in the ascidian Ciona intestinalis showing that the peripheral sensory neurons, whose specification is in part regulated by the coordinate activity of Notch-Delta signaling and the microRNA miR-124, exhibit a sparse spatial pattern whereby consecutive neurons may be spaced over a dozen cells apart. We perform rigorous stability and bifurcation analyses, and demonstrate that our model is able to accurately explain and reproduce the neuronal pattern in Ciona. Using Monte Carlo simulations of our model along with miR-124 transgene over-expression assays, we demonstrate that the activity of miR-124 can be incorporated into the Notch decay rate parameter of our model. Finally, we motivate the general applicability of our model to Notch-Delta signaling in other animals by providing evidence that microRNAs regulate Notch-Delta signaling in analogous cell types in other organisms, and by discussing evidence in other organisms of sparse spatial patterns in tissues where Notch-Delta signaling is active.

Identification of testis-specific ubiquitin-conjugating enzyme in the ascidian Ciona intestinalis.

The ubiquitin-proteasome system is known to play a key role in fertilization in ascidians, sea urchins, and mammals. To obtain insights into the ubiquitin-conjugating enzymes (Ube2) involved in reproductive systems, we systematically explored Ube2 enzymes expressed in the testis of the ascidian Ciona intestinalis. Here, we report cDNA cloning and characterization of a novel type of Ube2r (Ci0100152677) that is capable of making a thiolester bond with ubiquitin. Northern analysis, whole-mount in situ hybridization and immunocytochemistry indicate that this enzyme is exclusively expressed in the testis, mainly in the germ cells during the late stage of spermatogenesis, and is localized in the sperm head and tail, suggesting possible participation in fertilization or spermatogenesis/spermiogenesis.

MGEScan-non-LTR: computational identification and classification of autonomous non-LTR retrotransposons in eukaryotic genomes.

Computational methods for genome-wide identification of mobile genetic elements (MGEs) have become increasingly necessary for both genome annotation and evolutionary studies. Non-long terminal repeat (non-LTR) retrotransposons are a class of MGEs that have been found in most eukaryotic genomes, sometimes in extremely high numbers. In this article, we present a computational tool, MGEScan-non-LTR, for the identification of non-LTR retrotransposons in genomic sequences, following a computational approach inspired by a generalized hidden Markov model (GHMM). Three different states represent two different protein domains and inter-domain linker regions encoded in the non-LTR retrotransposons, and their scores are evaluated by using profile hidden Markov models (for protein domains) and Gaussian Bayes classifiers (for linker regions), respectively. In order to classify the non-LTR retrotransposons into one of the 12 previously characterized clades using the same model, we defined separate states for different clades. MGEScan-non-LTR was tested on the genome sequences of four eukaryotic organisms, Drosophila melanogaster, Daphnia pulex, Ciona intestinalis and Strongylocentrotus purpuratus. For the D. melanogaster genome, MGEScan-non-LTR found all known 'full-length' elements and simultaneously classified them into the clades CR1, I, Jockey, LOA and R1. Notably, for the D. pulex genome, in which no non-LTR retrotransposon has been annotated, MGEScan-non-LTR found a significantly larger number of elements than did RepeatMasker, using the current version of the RepBase Update library. We also identified novel elements in the other two genomes, which have only been partially studied for non-LTR retrotransposons.

Integrative assessment of coastal pollution in a Ría coastal system (Galicia, NW Spain): correspondence between sediment chemistry and toxicity.

Elutriate embryo-larval bioassays with sea-urchins (Paracentrotus lividus) and ascidians (Ciona intestinalis) were conducted concurrently with trace metal analyses as part of an integrative assessment of sediment pollution at Ría de Pontevedra (Galicia, NW Spain). High metal contents in sediments were found in localised areas from the inner part of the estuary indicating a clear anthropogenic influence. In particular, very high Cu, Zn and Pb levels were found at sites P2 and P3, which were also the most toxic to the embryo-larval bioassays. Sediment quality guidelines were used to help in the ecological interpretation of sediment chemistry data and to identify pollutants of concern. Cu and Zn in P3 were consistently above the effects range median (ERM) values, which seem to be good predictors of toxicity to sea-urchin and ascidian embryos. A toxic unit approach, based on published EC(50) values and metal levels in elutriates, was used to assess the harmful ecological effects associated to sediment chemistry. Toxicity detected in P3 may be explained on the basis of the toxic unit model; however, the high toxicity detected at P2 may be attributable not only to the metals quantified in the analyses but also to unmeasured organic pollutants. Multidimensional scaling applied independently to the toxicology and chemistry data resulted in a good agreement between both types of configurations. Moreover, the Mantel test revealed a significant correlation (r(M)=0.481; p=0.019) between metal concentrations and toxicity data profiles, supporting the correspondence between configurations.

Markov chain-based promoter structure modeling for tissue-specific expression pattern prediction.

Transcriptional regulation is the first level of regulation of gene expression and is therefore a major topic in computational biology. Genes with similar expression patterns can be assumed to be co-regulated at the transcriptional level by promoter sequences with a similar structure. Current approaches for modeling shared regulatory features tend to focus mainly on clustering of cis-regulatory sites. Here we introduce a Markov chain-based promoter structure model that uses both shared motifs and shared features from an input set of promoter sequences to predict candidate genes with similar expression. The model uses positional preference, order, and orientation of motifs. The trained model is used to score a genomic set of promoter sequences: high-scoring promoters are assumed to have a structure similar to the input sequences and are thus expected to drive similar expression patterns. We applied our model on two datasets in Caenorhabditis elegans and in Ciona intestinalis. Both computational and experimental verifications indicate that this model is capable of predicting candidate promoters driving similar expression patterns as the input-regulatory sequences. This model can be useful for finding promising candidate genes for wet-lab experiments and for increasing our understanding of transcriptional regulation.

Accuracy assessment of diploid consensus sequences.

If the origins of fragments are known in genome sequencing projects, it is straightforward to reconstruct diploid consensus sequences. In reality, however, this is not true. Although there are proposed methods to reconstruct haplotypes from genome sequencing projects, an accuracy assessment is required to evaluate the confidence of the estimated diploid consensus sequences. In this paper, we define the confidence score of diploid consensus sequences. It requires the calculation of the likelihood of an assembly. To calculate the likelihood, we propose a linear time algorithm with respect to the number of polymorphic sites. The likelihood calculation and confidence score are used for further improvements of haplotype estimation in two directions. One direction is that low-scored phases are disconnected. The other direction is that, instead of using nominal frequency 1/2, the haplotype frequency is estimated to reflect the actual contribution of each haplotype. Our method was evaluated on the simulated data whose polymorphism rate (1.2 percent) was based on Ciona intestinalis. As a result, the high accuracy of our algorithm was indicated: The true positive rate of the haplotype estimation was greater than 97 percent.

Isolation and phylogeny of novel cytochrome P450 genes from tunicates (Ciona spp.): a CYP3 line in early deuterostomes?

Cytochromes P450 (CYPs) form a gene superfamily involved in the biotransformation of numerous endogenous and exogenous natural and synthetic compounds. In humans, CYP3A4 is regarded as one of the most important CYPs due to its abundance in liver and its capacity to metabolize more than 50% of all clinically used drugs. It has been suggested that all CYP3s arose from a common ancestral gene lineage that diverged between 800 and 1100 million years ago, before the deuterostome-protostome split. While CYP3s are well known in mammals and have been described in lower vertebrates, they have not been reported in non-vertebrate deuterostomes. Members of the genus Ciona belong to the tunicates, whose lineage is thought to be the most basal among the chordates, and from which the vertebrate line diverged. Here we describe the cloning, exon-intron structure, phylogeny, and estimated expression of four novel genes from Ciona intestinalis. We also describe the gene structure and phylogeny of homologous genes in Ciona savignyi. Comparing these genes with other members of the CYP clan 3, show that the Ciona sequences bear remarkable similarity to vertebrate CYP3A genes, and may be an early deuterostome CYP3 line.

Toxicity assessment of the antifouling compound zinc pyrithione using early developmental stages of the ascidian Ciona intestinalis.

This study investigated the toxicity of zinc pyrithione (Zpt) on the early stages of development of the ascidian Ciona intestinalis. Larval morphological abnormalities were studied after the exposure of C. intestinalis embryos at different stages of development. The median effective concentrations (EC50) ranged from 226-590 nM. The larval settlement stage was the most sensitive to Zpt. Toxic effects of Zpt on larval settlement were detected at 9 nM (EC10). The inhibition of C. intestinalis embryonic development was also used to study the loss of toxicity in Zpt solutions exposed to direct sunlight and laboratory UV light. The results showed that the toxicity of Zpt solutions decreased but did not disappear after 4 h exposure to direct sunlight (EC50 = 484 nM) or UV light (EC50 = 453 nM), compared to control Zpt solutions prepared in dark conditions. On the basis of the present data, predicted no effect concentrations of Zpt to C. intestinalis larvae are lower than predicted environmental concentrations of Zpt in certain polluted areas and therefore, may pose a risk to C. intestinalis populations.

Real-time quantitative PCR for assessment of abundance of Pseudoalteromonas species in marine samples.

A real-time quantitative PCR (RTQ-PCR) method for measuring the abundance of Pseudoalteromonas species in marine samples is presented. PCR primers targeting a Pseudoalteromonas-specific region of the 16S rRNA gene were tested at three different levels using database searches (in silico), a selection of pure cultures (in vitro), and a combined denaturing gradient gel electrophoresis and cloning approach on environmental DNA (in situ). The RTQ-PCR method allowed for the detection of SYBR Green fluorescence from double-stranded DNA over a linear range spanning six orders of magnitude. The detection limit was determined as 1.4 fg of target DNA (1,000 gene copies) measured in the presence of 20 ng of nontarget DNA from salmon testes. In this study, we discuss the importance of robust post-PCR analyses to overcome pitfalls in RTQ-PCR when samples from different complex marine habitats are analyzed and compared on a nonroutine basis. Representatives of the genus Pseudoalteromonas were detected in samples from all investigated habitats, suggesting a widespread distribution of this genus across many marine habitats (e.g., seawater, rocks, macroalgae, and marine animals). Three sample types were analyzed by RTQ-PCR to determine the relative abundance of Pseudoalteromonas ribosomal DNA (rDNA) compared to the total abundance of eubacterial rDNA. The rDNA fractions of Pseudoalteromonas compared to all Eubacteria were 1.55% on the green alga Ulva lactuca, 0.10% on the tunicate Ciona intestinalis, and 0.06% on the green alga Ulvaria fusca.

Assessment of coastal marine pollution in Galicia (NW Iberian Peninsula); metal concentrations in seawater, sediments and mussels (Mytilus galloprovincialis) versus embryo-larval bioassays using Paracentrotus lividus and Ciona intestinalis.

Sediments from three Galician Rias were tested for toxicity using sea-urchin and ascidian sediment elutriate embryo-larval bioassays. Trace metal contents in seawater, sediments and mussels were also determined and subjected to multidimensional scaling methods which grouped stations according to chemical contamination. High metal contents were found in seawater, sediments and mussels from the Ria of Pontevedra, and moderate levels were detected in the Ria of Vigo and Ria of Arousa. The results revealed that samples assessed as toxic, according to the sea-urchin and ascidian embryo-larval bioassays, were among the most polluted by trace metals. A good agreement was reported between ordination plots resulting from applying multidimensional scaling to the chemical data, and the results of the biological endpoints tested.

Integrative assessment of marine pollution in Galician estuaries using sediment chemistry, mussel bioaccumulation, and embryo-larval toxicity bioassays.

An integrative assessment of environmental quality was carried out in selected sites along the Galician coast (NW Iberian Peninsula) combining analytical chemistry of seawater and sediments, bioaccumulation in the marine mussel, and embryo-larval sediment toxicity bioassays, in order to link biological and chemical criteria for the assessment of coastal pollution. Maximum values of Hg and Cu in seawater, sediment and mussels, were found in the inner part of Ria of Pontevedra, while maximum levels of organics (polychlorinated biphenyls, hexachlorobenzene and aldrin) were found in mussels from A Coruña. Outstanding values of Cu, Pb and Zn have been found in seawater and sediment from a single site, P3, which also was the most toxic in the embryo-larval bioassays performed with four different phyla of marine organisms: mollusks, echinoderms, arthropods and chordates. Sediment quality effects range-median values provided a valuable reference to predict biological effects from sediment chemistry data, while effects range-low values were too conservative. Sediment toxicity could also be predicted by using a toxic-unit model based on published EC50 values for trace metals and mobilization factors independently obtained from measurements of metal contents in sediments and their elutriates. When chemical and toxicological data are independently used to arrange sampling sites by using non-metric multidimensional scaling, a remarkable degree of concordance between both types of configurations could be observed.

Toxicity of Hg, Cu, Cd, and Cr on early developmental stages of Ciona intestinalis (Chordata, Ascidiacea) with potential application in marine water quality assessment.

The toxicity of mercury, copper, cadmium and chromium on sperm viability, fertilisation, embryogenesis and larval attachment of Ciona intestinalis was examined. Fertilisation rate (FR) showed a small decrease even at the highest metal concentration tested. The median effective concentrations (EC50) reducing rates of embryogenesis and larval attachment by 50% were 54 microg Hg/l (0.27 microM), 46 microg Cu/l (0.72 microM), 838 microg Cd/l (7.46 microM), 10,318 microg Cr/l (198 microM), and 35 microg Hg/l (0.18 microM), 34 microg Cu/l (0.54 microM) and 11,755 microg Cr/l (226 microM), respectively. Therefore, Hg is three times more toxic than Cu (on a molar basis), ca. 30 times more toxic than Cd and ca. 1000 times more toxic than Cr to early stages of C. intestinalis. Rates of larval attachment and embryogenesis were the most sensitive endpoints, although the latter is more advisable for routine assessment of seawater quality because of its greater simplicity. In addition to bivalves and sea-urchins, ascidian embryos can provide biological criteria for seawater quality standards taking into account the sensitivity of a chordate and contributing to the detection of harmful chemicals with no marked effect on the species currently in use in seawater quality bioassays.