Molecular taxonomy and phylogenetic affinities of two groundwater amphipods, Crangonyx islandicus and Crymostygius thingvallensis, endemic to Iceland

The amphipod superfamily Crangonyctoidea is distributed exclusively in freshwater habitats worldwide and is characteristic of subterranean habitats. Two members of the family, Crangonyx islandicus and Crymostygius thingvallensis, are endemic to Iceland and were recently discovered in groundwater underneath lava fields. Crangonyx islandicus belongs to a well-known genus with representatives both in North America and in Eurasia. Crymostygius thingvallensis defines a new family, Crymostygidae. Considering the incongruences observed recently between molecular and morphological taxonomy within subterranean species, we aim to assess the taxonomical status of the two species using molecular data. Additionally, the study contributes to the phylogenetic relationships among several crangonyctoidean species and specifically among species from four genera of the family Crangonyctidae. Given the available data we consider how the two Icelandic species could have colonized Iceland, by comparing geographical origin of the species with the phylogeny. Regions of two nuclear (18S and 28S rRNA) and two mitochondrial genes (16S rRNA and COI) for 20 different species of three families of the Crangonyctoidea were sequenced. Four different methods were used to align the RNA gene sequences and phylogenetic trees were constructed using bayesian and maximum likelihood analysis. The Crangonyctidae monophyly is supported. Crangonyx islandicus appeared more closely related to species from the Nearctic region. Crymostygius thingvallensis is clearly divergent from the other species of Crangonyctoidea. Crangonyx and Synurella genera are clearly polyphyletic and showed a geographical association, being split into a Nearctic and a Palearctic group. This research confirms that the studied species of Crangonyctidae share a common ancestor, which was probably widespread in the Northern hemisphere well before the break up of Laurasia. The Icelandic species are of particular interest since Iceland emerged after the separation of Eurasia and North America, is geographically isolated and has repeatedly been covered by glaciers during the Ice Age. The close relation between Crangonyx islandicus and North American species supports the hypothesis of the Trans-Atlantic land bridge between Greenland and Iceland which might have persisted until 6 million years ago. The status of the family Crymostygidae is supported, whereas Crangonyx islandicus might represent a new genus. As commonly observed in subterranean animals, molecular and morphological taxonomy led to different conclusions, probably due to convergent evolution of morphological traits. Our molecular analysis suggests that the family Crangonyctidae needs taxonomic revisions.     

A Computational Module Assembled from Different Protease Family Motifs Identifies PI PLC from Bacillus cereus as a Putative Prolyl Peptidase with a Serine Protease Scaffold

Proteolytic enzymes have evolved several mechanisms to cleave peptide bonds. These distinct types have been systematically categorized in the MEROPS database. While a BLAST search on these proteases identifies homologous proteins, sequence alignment methods often fail to identify relationships arising from convergent evolution, exon shuffling, and modular reuse of catalytic units. We have previously established a computational method to detect functions in proteins based on the spatial and electrostatic properties of the catalytic residues (CLASP). CLASP identified a promiscuous serine protease scaffold in alkaline phosphatases (AP) and a scaffold recognizing a β-lactam (imipenem) in a cold-active Vibrio AP. Subsequently, we defined a methodology to quantify promiscuous activities in a wide range of proteins. Here, we assemble a module which encapsulates the multifarious motifs used by protease families listed in the MEROPS database. Since APs and proteases are an integral component of outer membrane vesicles (OMV), we sought to query other OMV proteins, like phospholipase C (PLC), using this search module. Our analysis indicated that phosphoinositide-specific PLC from Bacillus cereus is a serine protease. This was validated by protease assays, mass spectrometry and by inhibition of the native phospholipase activity of PI-PLC by the well-known serine protease inhibitor AEBSF (IC50 = 0.018 mM). Edman degradation analysis linked the specificity of the protease activity to a proline in the amino terminal, suggesting that the PI-PLC is a prolyl peptidase. Thus, we propose a computational method of extending protein families based on the spatial and electrostatic congruence of active site residues.     

Mixed Model with Correction for Case-Control Ascertainment Increases Association Power

We introduce a liability-threshold mixed linear model (LTMLM) association statistic for case-control studies and show that it has a well-controlled false-positive rate and more power than existing mixed-model methods for diseases with low prevalence. Existing mixed-model methods suffer a loss in power under case-control ascertainment, but no solution has been proposed. Here, we solve this problem by using a χ² score statistic computed from posterior mean liabilities (PMLs) under the liability-threshold model. Each individual’s PML is conditional not only on that individual’s case-control status but also on every individual’s case-control status and the genetic relationship matrix (GRM) obtained from the data. The PMLs are estimated with a multivariate Gibbs sampler; the liability-scale phenotypic covariance matrix is based on the GRM, and a heritability parameter is estimated via Haseman-Elston regression on case-control phenotypes and then transformed to the liability scale. In simulations of unrelated individuals, the LTMLM statistic was correctly calibrated and achieved higher power than existing mixed-model methods for diseases with low prevalence, and the magnitude of the improvement depended on sample size and severity of case-control ascertainment. In a Wellcome Trust Case Control Consortium 2 multiple sclerosis dataset with >10,000 samples, LTMLM was correctly calibrated and attained a 4.3% improvement (p = 0.005) in χ² statistics over existing mixed-model methods at 75 known associated SNPs, consistent with simulations. Larger increases in power are expected at larger sample sizes. In conclusion, case-control studies of diseases with low prevalence can achieve power higher than that in existing mixed-model methods.     

Dynamics fingerprint and inherent asymmetric flexibility of a cold-adapted homodimeric enzyme. A case study of the Vibrio alkaline phosphatase

Background

Protein dynamics influence protein function and stability and modulate conformational changes. Such motions depend on the underlying networks of intramolecular interactions and communicating residues within the protein structure. Here, we provide the first characterization of the dynamic fingerprint of the dimeric alkaline phosphatase (AP) from the cold-adapted Vibrio strain G15-21 (VAP), which is among the APs with the highest known k_cat at low temperatures.

Methods

Multiple all-atom explicit solvent molecular dynamics simulations were employed in conjunction with different metrics to analyze the dynamical patterns and the paths of intra- and intermolecular communication.

Results

Interactions and coupled motions at the interface between the two VAP subunits have been characterized, along with the networks of intramolecular interactions. It turns out a low number of intermolecular interactions and coupled motions, which result differently distributed in the two monomers. The paths of long-range communication mediated from the catalytic residues to distal sites were also characterized, pointing out a different information flow in the two subunits.

Conclusions

A pattern of asymmetric flexibility is evident in the two identical subunits of the VAP dimer that is intimately linked to a different distribution of intra- and intermolecular interactions. The asymmetry was also evident in pairs of cross-correlated residues during the dynamics.

General significance

The results here discussed provide a structural rationale to the half-of-site mechanism previously proposed for VAP and other APs, as well as a general framework to characterize asymmetric dynamics in homomeric enzymes.     

Effect of Genetic Variation in a Drosophila Model of Diabetes-Associated Misfolded Human Proinsulin

The identification and validation of gene–gene interactions is a major challenge in human studies. Here, we explore an approach for studying epistasis in humans using a Drosophila melanogaster model of neonatal diabetes mellitus. Expression of the mutant preproinsulin (hINS^C96Y) in the eye imaginal disc mimics the human disease: it activates conserved stress-response pathways and leads to cell death (reduction in eye area). Dominant-acting variants in wild-derived inbred lines from the Drosophila Genetics Reference Panel produce a continuous, highly heritable distribution of eye-degeneration phenotypes in a hINS^C96Y background. A genome-wide association study (GWAS) in 154 sequenced lines identified a sharp peak on chromosome 3L, which mapped to a 400-bp linkage block within an intron of the gene sulfateless (sfl). RNAi knockdown of sfl enhanced the eye-degeneration phenotype in a mutant-hINS-dependent manner. RNAi against two additional genes in the heparan sulfate (HS) biosynthetic pathway (ttv and botv), in which sfl acts, also modified the eye phenotype in a hINS^C96Y-dependent manner, strongly suggesting a novel link between HS-modified proteins and cellular responses to misfolded proteins. Finally, we evaluated allele-specific expression difference between the two major sfl-intronic haplotypes in heterozygtes. The results showed significant heterogeneity in marker-associated gene expression, thereby leaving the causal mutation(s) and its mechanism unidentified. In conclusion, the ability to create a model of human genetic disease, map a QTL by GWAS to a specific gene, and validate its contribution to disease with available genetic resources and the potential to experimentally link the variant to a molecular mechanism demonstrate the many advantages Drosophila holds in determining the genetic underpinnings of human disease.     

Characterisation and expression of the paired box protein 7 (Pax7) gene in polymorphic Arctic charr (Salvelinus alpinus)

Arctic charr (Salvelinus alpinus L.) from Lake Thingvallavatn, Iceland occur as four distinct morphs: large benthivorous (LB), dwarf benthivorous (DB), piscivorous (PI) and planktonivorous (PL). The morphs differ with respect to body size, head morphology, growth rate, and life history. The aim of this study was to investigate the paired box protein 7 (Pax7) gene as a candidate for such polymorphisms due to its importance in cranio-facial, skeletal muscle, and central nervous system development. No variation in coding and intronic sequences was found between morphs. We identified 10 alternate Pax7 isoforms with insertions/deletions: a four-residue (GNRT) deletion, a GEASS insertion truncated by the first serine residue (GEAS), and a thirteen-residue insertion (GQYA/TGPEYVYCGT). The latter insertion with a threonine (T) contains a putative casein kinase II (CK-2) phosphorylation site. Pax7 spatial expression patterns were identical in embryos of DB-, LB-, and PL-morphs, and were similar to those described for zebrafish Pax7c, but a difference in temporal expression for segmentation was observed between DB and LB morphs. At the end of segmentation, novel expression was observed in the mandibular region as two bilateral domains. The potential role of multiple alternative splicing of the Pax7 gene for the generation of different Arctic charr morphs is briefly discussed.     

Microscopic rate-constants for substrate binding and acylation in cold-adaptation of trypsin I from Atlantic cod

Temperature imposes limits on where life can thrive and this is evident in the evolution of the basic structural properties of proteins. Cold-adaptation of enzymes is one example, where the catalytic rate constant (k(cat)) is increased compared with hot-acclimated homologous under identical assay conditions. Trypsin I from Atlantic cod (Gadus morhua) has catalytic efficiency (k(cat)/K(m)) for amide hydrolysis that is 17-fold larger than observed for bovine trypsin. Here, the individual rate-constants for association of substrate (k(1)), dissociation of substrate (k(-1)), and acylation of the enzyme (k(2)) have been determined using benzoyl-Arg-p-nitroanilide or benzyloxycarbonyl-Gly-Pro-Arg-p-nitroanilide as substrates. Rather unexpectedly, by far the largest difference (37-fold increase) was observed in k(1), the rate constant for binding of substrate. The cold-adaptation of the dissociation and catalytic steps were not as prominent (increased by 3.7-fold). The length of substrate did have an effect by increasing the reaction rate by 70-fold, and again, the step most affected was the initial binding-step.     

Islands of tractability for parsimony haplotyping

We study the parsimony approach to haplotype inference, which calls for finding a set of haplotypes of minimum cardinality that explains an input set of genotypes. We prove that the problem is APX-hard even in very restricted cases. On the positive side, we identify islands of tractability for the problem, by focusing on instances with specific structure of haplotype sharing among the input genotypes. We exploit the structure of those instance to give polynomial and constant-approximation algorithms to the problem. We also show that the general parsimony haplotyping problem is fixed parameter tractable.     

Natural hybrids in Atlantic eels (Anguilla anguilla, A. rostrata): evidence for successful reproduction and fluctuating abundance in space and time

The outcome of natural hybridization is highly variable and depends on the nonexclusive effects of both pre- and post-mating reproductive barriers. The objective of this study was to address three specific questions regarding the dynamics of hybridization between the American and European eels (Anguilla rostrata and Anguilla anguilla). Using 373 AFLP loci, 1127 eels were genotyped, representing different life stages from both continents, as well as multiple Icelandic locations. We first evaluated the extent of hybridization and tested for the occurrence of hybrids beyond the first generation. Second, we tested whether hybrids were randomly distributed across continents and among Icelandic sampling sites. Third, we tested for a difference in the proportion of hybrids between glass eel and yellow eel stages in Iceland. Our results provided evidence for (i) an overall hybrid proportion of 15.5% in Iceland, with values ranging from 6.7% to 100% depending on life stages and locations; (ii) the existence of hybrids beyond the first generation; (iii) a nonrandom geographic distribution of hybrids in the North Atlantic; and (iv) a higher proportion of first and later generation hybrids in yellow eels compared to glass eels, as well as a significant latitudinal gradient in the proportion of hybrids in Icelandic freshwater. We propose that the combined effect of both differential survival of hybrids and variation in hybridization rate through time best explain these patterns. We discuss the possibility that climate change, which is impacting many environmental features in the North Atlantic, may have a determinant effect on the outcome of natural hybridization in Atlantic eels.