Efficient linkage mapping using exome capture and extreme QTL in schistosome parasites

Background

Identification of parasite genes that underlie traits such as drug resistance and host specificity is challenging using classical linkage mapping approaches. Extreme QTL (X-QTL) methods, originally developed by rodent malaria and yeast researchers, promise to increase the power and simplify logistics of linkage mapping in experimental crosses of schistosomes (or other helminth parasites), because many 1000s of progeny can be analysed, phenotyping is not required, and progeny pools rather than individuals are genotyped. We explored the utility of this method for mapping a drug resistance gene in the human parasitic fluke Schistosoma mansoni.

Results

We staged a genetic cross between oxamniquine sensitive and resistant parasites, then between two F1 progeny, to generate multiple F2 progeny. One group of F2s infecting hamsters was treated with oxamniquine, while a second group was left untreated. We used exome capture to reduce the size of the genome (from 363 Mb to 15 Mb) and exomes from pooled F2 progeny (treated males, untreated males, treated females, untreated females) and the two parent parasites were sequenced to high read depth (mean = 95-366×) and allele frequencies at 14,489 variants compared. We observed dramatic enrichment of alleles from the resistant parent in a small region of chromosome 6 in drug-treated male and female pools (combined analysis: = 11.07, p = 8.74 × 10^-29). This region contains Smp_089320 a gene encoding a sulfotransferase recently implicated in oxamniquine resistance using classical linkage mapping methods.

Conclusions

These results (a) demonstrate the utility of exome capture for generating reduced representation libraries in Schistosoma mansoni, and (b) provide proof-of-principle that X-QTL methods can be successfully applied to an important human helminth. The combination of these methods will simplify linkage analysis of biomedically or biologically important traits in this parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-617) contains supplementary material, which is available to authorized users.     

Molecular phylogeny of the major arthropod groups indicates polyphyly of crustaceans and a new hypothesis for the origin of hexapods

A phylogeny of the arthropods was inferred from analyses of amino acid sequences derived from the nuclear genes encoding elongation factor-1 alpha and the largest subunit of RNA polymerase II using maximum- parsimony, neighbor-joining, and maximum-likelihood methods. Analyses of elongation factor-1 alpha from 17 arthropods and 4 outgroup taxa recovered many arthropod clades supported by previous morphological studies, including Diplopoda, Myriapoda, Insecta, Hexapoda, Branchiopoda (Crustacea), Araneae, Tetrapulmonata, Arachnida, Chelicerata, and Malacostraca (Crustacea). However, counter to previous studies, elongation factor-1 alpha placed Malacostraca as sister group to the other arthropods. Branchiopod crustaceans were found to be more closely related to hexapods and myriapods than to malacostracan crustaceans. Sequences for RNA polymerase II were obtained from 11 arthropod taxa and were analyzed separately and in combination with elongation factor-1 alpha. Results from these analyses were concordant with those derived from elongation factor-1 alpha alone and provided support for a Hexapoda/Branchiopoda clade, thus arguing against the monophyly of the traditionally defined Atelocerata (Hexapoda + Myriapoda).      

Statistical methods of DNA sequence analysis: detection of intragenic recombination or gene conversion

Simple but exact statistical tests for detecting a cluster of associated nucleotide changes in DNA are presented. The tests are based on the linear distribution of a set of s sites among a total of n sites, where the s sites may be the variable sites, sites of insertion/deletion, or categorized in some other way. These tests are especially useful for detecting gene conversion and intragenic recombination in a sample of DNA sequences. In this case, the sites of interest are those that correspond to particular ways of splitting the sequences into two groups (e.g., sequences A and D vs. sequences B, C, and E-J). Each such split is termed a phylogenetic partition. Application of these methods to a well-documented case of gene conversion in human gamma-globin genes shows that sites corresponding to two of the three observed partitions are significantly clustered, whereas application to hominoid mitochondrial DNA sequences--among which no recombination is expected to occur--shows no evidence of such clustering. This indicates that clustering of partition-specific sites is largely due to intragenic recombination or gene conversion. Alternative hypotheses explaining the observed clustering of sites, such as biased selection or mutation, are discussed.      

Acute effects of prolonged intermittent low-intensity isometric warm-up schemes on jump,sprint, and agility performance in collegiate soccer players

The aim of the present study was to compare the effects of different warm-up interventions on jump, sprint and agility performance in collegiate soccer players. Twenty-one healthy male college soccer players (age: 20.14 ± 1.65 years; body height: 179.9 ± 8.34 cm; body mass: 74.4 ± 13.0 kg; % body fat: 9.45 ± 4.8) participated in the study. Subjects underwent four different randomized warm-up protocols separated by at least 48 hours. The warm-up schemes were: 1. no conditioning contraction protocol (NCC); 2. dynamic stretching (DS); 3. prolonged intermittent low-intensity isometric exercise (ST); and, 4. ST with an additional external load equal to 30% of body weight (ST + 30% BW). All interventions were preceded by a general warm-up. Results from one-way repeated measures ANOVA demonstrated a significant difference in countermovement jump (CMJ) at F(3,60) = 10.2, ηρ² = 0.337, p < 0.01. Post hoc analysis revealed a significant difference in CMJ performance in DS when compared to NCC and ST + 30% BW. No significant difference in CMJ was observed between DS and ST. CMJ scores in NCC, ST, and ST + 30% BW were non-significant. There was a significant difference in speed; F(3, 60) = 6.61, ηρ² = 0.248, p < 0.01. Post hoc analysis revealed significantly better time in DS than NCC and ST. However, no difference in speed was observed between DS and ST + 30% BW. Similarly, speed was similar in NCC, ST and ST + 30% BW. A significant difference in agility performance was also observed; F(3, 60) = 24.1, ηρ²= 0.546, p < 0.01. Post hoc analysis revealed significantly greater performance gains in DS than NCC. No significant difference in agility was observed in DS, ST and ST + 30% BW. In conclusion, a prolonged intermittent low-intensity isometric protocol using bodyweight only showed similar benefits with dynamic stretching in countermovement jump performance. When the same isometric condition with additional load equal to 30% of bodyweight was applied, effects in speed and agility were similar to dynamic stretching.     

Rapid screening for phenotype-genotype associations by linear transformations of genomic evaluations

Background

Currently, association studies are analysed using statistical mixed models, with marker effects estimated by a linear transformation of genomic breeding values. The variances of marker effects are needed when performing the tests of association. However, approaches used to estimate the parameters rely on a prior variance or on a constant estimate of the additive variance. Alternatively, we propose a standardized test of association using the variance of each marker effect, which generally differ among each other. Random breeding values from a mixed model including fixed effects and a genomic covariance matrix are linearly transformed to estimate the marker effects.

Results

The standardized test was neither conservative nor liberal with respect to type I error rate (false-positives), compared to a similar test using Predictor Error Variance, a method that was too conservative. Furthermore, genomic predictions are solved efficiently by the procedure, and the p-values are virtually identical to those calculated from tests for one marker effect at a time. Moreover, the standardized test reduces computing time and memory requirements.The following steps are used to locate genome segments displaying strong association. The marker with the highest − log(p-value) in each chromosome is selected, and the segment is expanded one Mb upstream and one Mb downstream of the marker. A genomic matrix is calculated using the information from those markers only, which is used as the variance-covariance of the segment effects in a model that also includes fixed effects and random genomic breeding values. The likelihood ratio is then calculated to test for the effect in every chromosome against a reduced model with fixed effects and genomic breeding values. In a case study with pigs, a significant segment from chromosome 6 explained 11% of total genetic variance.

Conclusions

The standardized test of marker effects using their own variance helps in detecting specific genomic regions involved in the additive variance, and in reducing false positives. Moreover, genome scanning of candidate segments can be used in meta-analyses of genome-wide association studies, as it enables the detection of specific genome regions that affect an economically relevant trait when using multiple populations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-246) contains supplementary material, which is available to authorized users.     

Thrombin promotes actin stress fiber formation in RPE through Rho/ROCK‐mediated MLC phosphorylation

The retinal pigment epithelium (RPE) forms the outer blood–retina barrier (BRB). Most retinal diseases involve BRB breakdown, whereupon thrombin contained in serum directly contacts the RPE. Thrombin is known to promote actin stress fiber formation, an important determinant in eye diseases involving the epithelial–mesenchymal transition (EMT) and migration of RPE cells, such as proliferative vitreoretinopathy. We analyzed thrombin effect on signaling pathways leading to myosin light chain (MLC) phosphorylation and actin stress fiber formation in primary cultures of rat RPE cells, in order to support a role for thrombin in RPE transdifferentiation. MLC phosphorylation was measured by Western blot; actin cytoskeleton was visualized using immunofluorescent phalloidin, and Rho GTPase activation was assessed by ELISA. We showed that thrombin/PAR‐1 induces the time‐ and dose‐dependent phosphorylation of MLC through the activation of Rho/ROCK and myosin light chain kinase (MLCK). ROCK increased phospho‐MLC by phosphorylating MLC and by inhibiting MLC phosphatase. Thrombin effect was abolished by the ROCK inhibitor Y‐27632, whereas MLCK inhibitor ML‐7 and PLC‐β inhibitor U73122 attenuated MLC phosphorylation by ≈50%, suggesting the activation of MLCK by PLC‐β‐mediated calcium increase. Additionally, thrombin‐induced MLC phosphorylation was blocked by the inhibitory PKCζ pseudosubstrate, wortmannin, and LY294002, indicating IP₃/PKCζ involvement in the control of MLC phosphorylation. Moreover, we demonstrated that thrombin effect on MLC induces actin stress fiber formation, since this effect was prevented by inhibiting the pathways leading to MLC phosphorylation. We conclude that thrombin stimulation of MLC phosphorylation and actin stress fiber formation may be involved in thrombin‐induced RPE cell transformation subsequent to BRB dysfunction. J. Cell. Physiol. 226: 414–423, 2011. © 2010 Wiley‐Liss, Inc.     

Conceptual-level workflow modeling of scientific experiments using NMR as a case study

Background  

Scientific workflows improve the process of scientific experiments by making computations explicit, underscoring data flow, and emphasizing the participation of humans in the process when intuition and human reasoning are required. Workflows for experiments also highlight transitions among experimental phases, allowing intermediate results to be verified and supporting the proper handling of semantic mismatches and different file formats among the various tools used in the scientific process. Thus, scientific workflows are important for the modeling and subsequent capture of bioinformatics-related data. While much research has been conducted on the implementation of scientific workflows, the initial process of actually designing and generating the workflow at the conceptual level has received little consideration.