Characterization of mucus-associated proteins from abalone (Haliotis) - candidates for chemical signaling

Living in groups is a widespread phenomenon in the animal kingdom. For free-spawning aquatic animals, such as the abalone (Haliotis), being in the close proximity to potential mating partners enhances reproductive success. In this study, we investigated whether chemical cues could be present in abalone mucus that enable species-specific aggregation. A comparative MS analysis of mucus obtained from trailing or fixed stationary Haliotis asinina, and from seawater surrounding aggregations, indicated that water-soluble biomolecules are present and that these can stimulate sensory activity in conspecifics. Purified extracts of trail mucus contain at least three small proteins [termed H.?asinina mucus-associated proteins (Has-MAPs)-1-3], which readily diffuse into the surrounding seawater and evoke a robust cephalic tentacle response in conspecifics. Mature Has-MAP-1 is approximately 9.9 kDa in size, and has a glycine-rich N-terminal region. Has-MAP-2 is approximately 6.2 kDa in size, and has similarities to schistosomin, a protein that is known to play a role in mollusc reproduction. The mature Has-MAP-3 is approximately 12.5 kDa in size, and could only be identified within trail mucus of animals outside of the reproductive season. All three Has-MAP genes are expressed at high levels within secretory cells of the juvenile abalone posterior pedal gland, consistent with a role in scent marking. We infer from these results that abalone mucus-associated proteins are candidate chemical cues that could provide informational cues to conspecifics living in close proximity and, given their apparent stability and hydrophilicity, animals further afield.     

A polynomial time algorithm for calculating the probability of a ranked gene tree given a species tree

ABSTRACT: BACKGROUND: The ancestries of genes form gene trees which do not necessarily have the same topology as the species tree due to incomplete lineage sorting. Available algorithms determining the probability of a gene tree given a species tree require exponential computational runtime. RESULTS: In this paper, we provide a polynomial time algorithm to calculate the probability of a ranked gene tree topology for a given species tree, where a ranked tree topology is a tree topology with the internal vertices being ordered. The probability of a gene tree topology can thus be calculated in polynomial time if the number of orderings of the internal vertices is a polynomial number. However, the complexity of calculating the probability of a gene tree topology with an exponential number of rankings for a given species tree remains unknown. CONCLUSIONS: Polynomial algorithms for calculating ranked gene tree probabilities may become useful in developing methodology to infer species trees based on a collection of gene trees, leading to a more accurate reconstruction of ancestral species relationships.     

Multiple merger gene genealogies in two species: Monophyly, paraphyly, and polyphyly for two examples of Lambda coalescents

Probabilities of monophyly, paraphyly, and polyphyly of two-species gene genealogies are computed for modest sample sizes and compared for two different Λ coalescent processes. Coalescent processes belonging to the Λ coalescent family admit asynchronous multiple mergers of active ancestral lineages. Assigning a timescale to the time of divergence becomes a central issue when different populations have different coalescent processes running on different timescales. Clade probabilities in single populations are also computed, which can be useful for testing for taxonomic distinctiveness of an observed set of monophyletic lineages. The coalescence rates of multiple merger coalescent processes are functions of coalescent parameters. The effect of coalescent parameters on the probabilities studied depends on the coalescent process, and if the population is ancestral or derived. The probability of reciprocal monophyly tends to be somewhat lower, when associated with a Λ coalescent, under the null hypothesis that two groups come from the same population. However, even for fairly recent divergence times, the probability of monophyly tends to be higher as a function of the number of generations for coalescent processes that admit multiple mergers, and is sensitive to the parameter of one of the example processes.     

The role of MAPK signaling in patterning and establishing axial symmetry in the gastropod Haliotis asinina

Gastropods are members of the Spiralia, a diverse group of invertebrates that share a common early developmental program, which includes spiral cleavage and a larval trochophore stage. The spiral cleavage program results in the division of the embryo into four quadrants. Specification of the dorsal (D) quadrant is intimately linked with body plan organization and in equally cleaving gastropods occurs when one of the vegetal macromeres makes contact with overlying micromeres and receives an inductive signal that activates a MAPK signaling cascade. Following the induction of the 3D macromere, the embryo begins to gastrulate and assumes a bilateral cleavage pattern. Here we inhibit MAPK activation in 3D with U0126 and examine its effect on the formation and patterning of the trochophore, using a suite of territory-specific markers. The head (pretrochal) region appears to maintain quadri-radial symmetry in U0126-treated embryos, supporting a role for MAPK signaling in 3D in establishing dorsoventral polarity in this region. Posterior (posttrochal) structures - larval musculature, shell and foot - fail to develop in MAPK inhibited trochophores. Inhibition of 3D specification by an alternative method - monensin treatment - yields similar abnormal trochophores. However, genes that are normally expressed in the ectodermal structures (shell and foot) are detected in U0126- and monensin-perturbed larvae in patterns that suggest that this region has latent dorsoventral polarity that is manifested even in the absence of D quadrant specification.     

Conservation and diversity in the immunity regions of wild phages with the immunity specificity of phage lambda

The gene regulatory circuitry of phage lambda is among the best-understood circuits. Much of the circuitry centres around the immunity region, which includes genes for two repressors, CI and Cro, and their cis-acting sites. Related phages, termed lambdoid phages, have different immunity regions, but similar regulatory circuitry and genome organization to that of lambda, and show a mosaic organization, arising by recombination between lambdoid phages. We sequenced the immunity regions of several wild phages with the immunity specificity of lambda, both to determine whether natural variation exists in regulation, and to analyse conservation and variability in a region rich in well-studied regulatory elements. CI, Cro and their cis-acting sites are almost identical to those in lambda, implying that regulatory mechanisms controlled by the immunity region are conserved. A segment adjacent to one of the operator regions is also conserved, and may be a novel regulatory element. In most isolates, different alleles of two regulatory proteins (N and CII) flank the immunity region; possibly the lysis-lysogeny decision is more variable among isolates. Extensive mosaicism was observed for several elements flanking the immunity region. Very short sequence elements or microhomologies were also identified. Our findings suggest mechanisms by which fine-scale mosaicism arises.     

Mitochondrial diversity of early-branching metazoa is revealed by the complete mt genome of a haplosclerid demosponge

The first mitochondrial (mt) genomes of demosponges have recently been sequenced and appear to be markedly different from published eumetazoan mt genomes. Here we show that the mt genome of the haplosclerid demosponge Amphimedon queenslandica has features that it shares with both demosponges and eumetazoans. Although the A. queenslandica mt genome has typical demosponge features, including size, long noncoding regions, and bacterialike rRNA genes, it lacks atp9, which is found in the other demosponges sequenced to date. We found strong evidence of a recent transposon-mediated transfer of atp9 to the nuclear genome. In addition, A. queenslandica bears an incomplete tRNA set, unusual amino acid deletion patterns, and a putative control region. Furthermore, the arrangement of mt rRNA genes differs from that of other demosponges. These genes evolve at significantly higher rates than observed in other demosponges, similar to previously observed nuclear rRNA gene rates in other haplosclerid demosponges.     

Identifying and quantifying metabolites by scoring peaks of GC-MS data

Background

Metabolomics is one of most recent omics technologies. It has been applied on fields such as food science, nutrition, drug discovery and systems biology. For this, gas chromatography-mass spectrometry (GC-MS) has been largely applied and many computational tools have been developed to support the analysis of metabolomics data. Among them, AMDIS is perhaps the most used tool for identifying and quantifying metabolites. However, AMDIS generates a high number of false-positives and does not have an interface amenable for high-throughput data analysis. Although additional computational tools have been developed for processing AMDIS results and to perform normalisations and statistical analysis of metabolomics data, there is not yet a single free software or package able to reliably identify and quantify metabolites analysed by GC-MS.

Results

Here we introduce a new algorithm, PScore, able to score peaks according to their likelihood of representing metabolites defined in a mass spectral library. We implemented PScore in a R package called MetaBox and evaluated the applicability and potential of MetaBox by comparing its performance against AMDIS results when analysing volatile organic compounds (VOC) from standard mixtures of metabolites and from female and male mice faecal samples. MetaBox reported lower percentages of false positives and false negatives, and was able to report a higher number of potential biomarkers associated to the metabolism of female and male mice.

Conclusions

Identification and quantification of metabolites is among the most critical and time-consuming steps in GC-MS metabolome analysis. Here we present an algorithm implemented in a R package, which allows users to construct flexible pipelines and analyse metabolomics data in a high-throughput manner.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0374-2) contains supplementary material, which is available to authorized users.