Interlaboratory evaluation of automated, multiplexed peptide immunoaffinity enrichment coupled to multiple reaction monitoring mass spectrometry for quantifying proteins in plasma

The inability to quantify large numbers of proteins in tissues and biofluids with high precision, sensitivity, and throughput is a major bottleneck in biomarker studies. We previously demonstrated that coupling immunoaffinity enrichment using anti-peptide antibodies (SISCAPA) to multiple reaction monitoring mass spectrometry (MRM-MS) produces Immunoprecipitation MRM-MS (immuno-MRM-MS) assays that can be multiplexed to quantify proteins in plasma with high sensitivity, specificity, and precision. Here we report the first systematic evaluation of the interlaboratory performance of multiplexed (8-plex) immuno-MRM-MS in three independent labs. A staged study was carried out in which the effect of each processing and analysis step on assay coefficient of variance, limit of detection, limit of quantification, and recovery was evaluated. Limits of detection were at or below 1 ng/ml for the assayed proteins in 30 μl of plasma. Assay reproducibility was acceptable for verification studies, with median intra- and interlaboratory coefficients of variance above the limit of quantification of 11% and <14%, respectively, for the entire immuno-MRM-MS assay process, including enzymatic digestion of plasma. Trypsin digestion and its requisite sample handling contributed the most to assay variability and reduced the recovery of target peptides from digested proteins. Using a stable isotope-labeled protein as an internal standard instead of stable isotope-labeled peptides to account for losses in the digestion process nearly doubled assay accuracy for this while improving assay precision 5%. Our results demonstrate that multiplexed immuno-MRM-MS can be made reproducible across independent laboratories and has the potential to be adopted widely for assaying proteins in matrices as complex as plasma.     

ARF6-Dependent Regulation of P2Y Receptor Traffic and Function in Human Platelets

Adenosine diphosphate (ADP) is a critical regulator of platelet activation, mediating its actions through two G protein-coupled receptors, the P2Y(1) and P2Y(12) purinoceptors. Recently, we demonstrated that P2Y(1) and P2Y(12) purinoceptor activities are rapidly and reversibly modulated in human platelets, revealing that the underlying mechanism requires receptor internalization and subsequent trafficking as an essential part of this process. In this study we investigated the role of the small GTP-binding protein ADP ribosylation factor 6 (ARF6) in the internalization and function of P2Y(1) and P2Y(12) purinoceptors in human platelets. ARF6 has been implicated in the internalization of a number of GPCRs, although its precise molecular mechanism in this process remains unclear. In this study we show that activation of either P2Y(1) or P2Y(12) purinoceptors can stimulate ARF6 activity. Further blockade of ARF6 function either in cell lines or human platelets blocks P2Y purinoceptor internalization. This blockade of receptor internalization attenuates receptor resensitization. Furthermore, we demonstrate that Nm23-H1, a nucleoside diphosphate (NDP) kinase regulated by ARF6 which facilitates dynamin-dependent fission of coated vesicles during endocytosis, is also required for P2Y purinoceptor internalization. These data describe a novel function of ARF6 in the internalization of P2Y purinoceptors and demonstrate the integral importance of this small GTPase upon platelet ADP receptor function.     

Metagenomics of the Svalbard reindeer rumen microbiome reveals abundance of polysaccharide utilization loci

Lignocellulosic biomass remains a largely untapped source of renewable energy predominantly due to its recalcitrance and an incomplete understanding of how this is overcome in nature. We present here a compositional and comparative analysis of metagenomic data pertaining to a natural biomass-converting ecosystem adapted to austere arctic nutritional conditions, namely the rumen microbiome of Svalbard reindeer (Rangifer tarandus platyrhynchus). Community analysis showed that deeply-branched cellulolytic lineages affiliated to the Bacteroidetes and Firmicutes are dominant, whilst sequence binning methods facilitated the assemblage of metagenomic sequence for a dominant and novel Bacteroidales clade (SRM-1). Analysis of unassembled metagenomic sequence as well as metabolic reconstruction of SRM-1 revealed the presence of multiple polysaccharide utilization loci-like systems (PULs) as well as members of more than 20 glycoside hydrolase and other carbohydrate-active enzyme families targeting various polysaccharides including cellulose, xylan and pectin. Functional screening of cloned metagenome fragments revealed high cellulolytic activity and an abundance of PULs that are rich in endoglucanases (GH5) but devoid of other common enzymes thought to be involved in cellulose degradation. Combining these results with known and partly re-evaluated metagenomic data strongly indicates that much like the human distal gut, the digestive system of herbivores harbours high numbers of deeply branched and as-yet uncultured members of the Bacteroidetes that depend on PUL-like systems for plant biomass degradation.     

Short-lived, transitory cell-cell interactions foster migration-dependent aggregation

During embryonic development, motile cells aggregate into cohesive groups, which give rise to tissues and organs. The role of cell migration in regulating aggregation is unclear. The current paradigm for aggregation is based on an equilibrium model of differential cell adhesivity to neighboring cells versus the underlying substratum. In many biological contexts, however, dynamics is critical. Here, we provide evidence that multicellular aggregation dynamics involves both local adhesive interactions and transport by cell migration. Using time-lapse video microscopy, we quantified the duration of cell-cell contacts among migrating cells that collided and adhered to another cell. This lifetime of cell-cell interactions exhibited a monotonic decreasing dependence on substratum adhesivity. Parallel quantitative measurements of cell migration speed revealed that across the tested range of adhesive substrata, the mean time needed for cells to migrate and encounter another cell was greater than the mean adhesion lifetime, suggesting that aggregation dynamics may depend on cell motility instead of the local differential adhesivity of cells. Consistent with this hypothesis, aggregate size exhibited a biphasic dependence on substratum adhesivity, matching the trend we observed for cell migration speed. Our findings suggest a new role for cell motility, alongside differential adhesion, in regulating developmental aggregation events and motivate new design principles for tuning aggregation dynamics in tissue engineering applications.     

Role of DDAH-1 in lipid peroxidation product-mediated inhibition of endothelial NO generation

Altered nitric oxide (NO) biosynthesis is thought to play a role in the initiation and progression of atherosclerosis and may contribute to increased risk seen in other cardiovascular diseases. It is hypothesized that altered NO bioavailability may result from an increase in endogenous NO synthase (NOS) inhibitors, asymmetric dimethly araginine (ADMA), and N(G)-monomethyl arginine, which are normally metabolized by dimethyarginine dimethylamine hydrolase (DDAH). Lipid hydroperoxides and their degradation products are generated during inflammation and oxidative stress and have been implicated in the pathogenesis of cardiovascular disorders. Here, we show that the lipid hydroperoxide degradation product 4-hydroxy-2-nonenal (4-HNE) causes a dose-dependent decrease in NO generation from bovine aortic endothelial cells, accompanied by a decrease in DDAH enzyme activity. The inhibitory effects of 4-HNE (50 microM) on endothelial NO production were partially reversed with L-Arg supplementation (1 mM). Overexpression of human DDAH-1 along with antioxidant supplementation completely restored endothelial NO production following exposure to 4-HNE (50 microM). These results demonstrate a critical role for the endogenous methylarginines in the pathogenesis of endothelial dysfunction. Because lipid hydroperoxides and their degradation products are known to be involved in atherosclerosis, modulation of DDAH and methylarginines may serve as a novel therapeutic target in the treatment of cardiovascular disorders associated with oxidative stress.     

Inferring the contribution of sexual reproduction,migration and off‐season survival to the temporal maintenance of microbial populations: a case study on the wheat fungal pathogen Puccinia striiformis f.sp. tritici

Understanding the mode of temporal maintenance of plant pathogens is an important domain of microbial ecology research. Due to the inconspicuous nature of microbes, their temporal maintenance cannot be studied directly through tracking individuals and their progeny. Here, we suggest a series of population genetic analyses on molecular marker variation in temporally spaced samples to infer about the relative contribution of sexual reproduction, off‐season survival and migration to the temporal maintenance of pathogen populations. We used the proposed approach to investigate the temporal maintenance of wheat yellow rust pathogen, Puccinia striiformis f.sp. tritici (PST), in the Himalayan region of Pakistan. Multilocus microsatellite genotyping of PST isolates revealed high genotypic diversity and recombinant population structure across all locations, confirming the existence of sexual reproduction in this region. The genotypes were assigned to four genetic groups, revealing a clear differentiation between zones with and without Berberis spp., the alternate host of PST, with an additional subdivision within the Berberis zone. The lack of any differentiation between samples across two sampling years, and the very infrequent resampling of multilocus genotypes over years at a given location was consistent with limited over‐year clonal survival, and a limited genetic drift. The off‐season oversummering population in the Berberis zone, likely to be maintained locally, served as a source of migrants contributing to the temporal maintenance in the non‐Berberis zone. Our study hence demonstrated the contribution of both sexual recombination and off‐season oversummering survival to the temporal maintenance of the pathogen. These new insights into the population biology of PST highlight the general usefulness of the analytical approach proposed.     

Predicting the likely response of data‐poor ecosystems to climate change using space‐for‐time substitution across domains

Predicting ecological response to climate change is often limited by a lack of relevant local data from which directly applicable mechanistic models can be developed. This limits predictions to qualitative assessments or simplistic rules of thumb in data‐poor regions, making management of the relevant systems difficult. We demonstrate a method for developing quantitative predictions of ecological response in data‐poor ecosystems based on a space‐for‐time substitution, using distant, well‐studied systems across an inherent climatic gradient to predict ecological response. Changes in biophysical data across the spatial gradient are used to generate quantitative hypotheses of temporal ecological responses that are then tested in a target region. Transferability of predictions among distant locations, the novel outcome of this method, is demonstrated via simple quantitative relationships that identify direct and indirect impacts of climate change on physical, chemical and ecological variables using commonly available data sources. Based on a limited subset of data, these relationships were demonstrably plausible in similar yet distant (>2000 km) ecosystems. Quantitative forecasts of ecological change based on climate‐ecosystem relationships from distant regions provides a basis for research planning and informed management decisions, especially in the many ecosystems for which there are few data. This application of gradient studies across domains – to investigate ecological response to climate change – allows for the quantification of effects on potentially numerous, interacting and complex ecosystem components and how they may vary, especially over long time periods (e.g. decades). These quantitative and integrated long‐term predictions will be of significant value to natural resource practitioners attempting to manage data‐poor ecosystems to prevent or limit the loss of ecological value. The method is likely to be applicable to many ecosystem types, providing a robust scientific basis for estimating likely impacts of future climate change in ecosystems where no such method currently exists.     

ROVER variant caller: read-pair overlap considerate variant-calling software applied to PCR-based massively parallel sequencing datasets

Background

We recently described Hi-Plex, a highly multiplexed PCR-based target-enrichment system for massively parallel sequencing (MPS), which allows the uniform definition of library size so that subsequent paired-end sequencing can achieve complete overlap of read pairs. Variant calling from Hi-Plex-derived datasets can thus rely on the identification of variants appearing in both reads of read-pairs, permitting stringent filtering of sequencing chemistry-induced errors. These principles underly ROVER software (derived from Read Overlap PCR-MPS variant caller), which we have recently used to report the screening for genetic mutations in the breast cancer predisposition gene PALB2. Here, we describe the algorithms underlying ROVER and its usage.

Results

ROVER enables users to quickly and accurately identify genetic variants from PCR-targeted, overlapping paired-end MPS datasets. The open-source availability of the software and threshold tailorability enables broad access for a range of PCR-MPS users.

Methods

ROVER is implemented in Python and runs on all popular POSIX-like operating systems (Linux, OS X). The software accepts a tab-delimited text file listing the coordinates of the target-specific primers used for targeted enrichment based on a specified genome-build. It also accepts aligned sequence files resulting from mapping to the same genome-build. ROVER identifies the amplicon a given read-pair represents and removes the primer sequences by using the mapping co-ordinates and primer co-ordinates. It considers overlapping read-pairs with respect to primer-intervening sequence. Only when a variant is observed in both reads of a read-pair does the signal contribute to a tally of read-pairs containing or not containing the variant. A user-defined threshold informs the minimum number of, and proportion of, read-pairs a variant must be observed in for a ‘call’ to be made. ROVER also reports the depth of coverage across amplicons to facilitate the identification of any regions that may require further screening.

Conclusions

ROVER can facilitate rapid and accurate genetic variant calling for a broad range of PCR-MPS users.     

Sub-alpine amphibian distributions related to species palatability to non-native salmonids in the Klamath mountains of northern California

The goal of this study was to examine how introduced trout influence the distributions and abundances of a sub‐alpine amphibian assemblage whose members display a variety of different life‐history and defence strategies. Our study was conducted in the sub‐alpine lentic habitats of three wilderness areas that form the core of the Klamath‐Siskiyou Bioregion of northern California, a biodiversity ‘hotspot’ that supports the highest diversity of sub‐alpine, lentic‐breeding amphibians in the western USA. These wilderness areas contain no native fishes, but all have been populated with non‐native trout for recreational fishing. Five of the eight amphibian species that occur in this region were sufficiently common to use in our study; these included one that breeds in both temporary and permanent waters and is palatable to fish (Pacific treefrog, Pseudacris regilla), two that breed primarily in permanent waters and are unpalatable to fish (western toad, Bufo boreas, and rough‐skinned newt, Taricha granulosa), and two that breed primarily in permanent waters and are palatable to fish (Cascades frog, Rana cascadae, and long‐toed salamander, Ambystoma macrodactylum). Based on life histories and predator defence strategies (i.e. palatable or not), we predicted that the three palatable species would likely be negatively correlated with introduced trout, but with P. regilla less impacted because of its use of both temporary and permanent waters. We predicted that B. boreas and T. granulosa would not be significantly correlated with introduced trout due to the lack of any predator/prey interactions between them. We surveyed 728 pond, lake, or wet meadow sites during the summers of 1999–2002, using timed gill‐net sets to measure trout occurrence and relative density, and visual encounter surveys to determine amphibian presence and abundance. We used semiparametric logistic regression models to quantify the effect of trout presence/absence and density on the probability of finding amphibian species in a water body while accounting for variation within and among the various lentic habitats sampled. The distributions of P. regilla, A. macrodactylum and R. cascadae were strongly negatively correlated with trout presence across all three wilderness areas. Ambystoma macrodactylum was 44 times more likely to be found in lakes without fish than in lakes with fish. Rana cascadae and P. regilla were 3.7 and 3.0 times more likely, respectively, to be found in fishless than fish‐containing waters. In contrast, the two unpalatable species were either uncorrelated (T. granulosa) or positively correlated (B. boreas) with fish presence. We found that the relative density of fish (catch per unit effort) was negatively correlated with the combined abundances of the three palatable amphibians, and also with both the length and the condition of the fish themselves. Our results are consistent with a compelling body of evidence that introduced fishes greatly alter the aquatic community structure of mountain lakes, ponds, and wet meadows.     

Isolation by distance and gene flow in the Eurasian badger (Meles meles) at both a local and broad scale

Eurasian badgers, Meles meles, have been shown to possess limited genetic population structure within Europe; however, field studies have detected high levels of philopatry, which are expected to increase population structure. Population structure will be a consequence of both contemporary dispersal and historical processes, each of which is expected to be evident at a different scale. Therefore, to gain a greater understanding of gene flow in the badger, we examined microsatellite diversity both among and within badger populations, focusing on populations from the British Isles and western Europe. We found that while populations differed in their allelic diversity, the British Isles displayed a similar degree of diversity to the rest of western Europe. The lower genetic diversity occurring in Ireland, Norway and Scotland was more likely to have resulted from founder effects rather than contemporary population density. While there was significant population structure (F ST = 0.19), divergence among populations was generally well explained by geographic distance (P < 0.0001) across the entire range studied of more than 3000 km. Transient effects from the Pleistocene appear to have been replaced by a strong pattern of genetic isolation by distance across western Europe, suggestive of colonization from a single refugium. Analysis of individuals within British populations through Mantel tests and spatial autocorrelation demonstrated that there was significant local population structure across 3-30 km, confirming that dispersal is indeed restricted. The isolation by distance observed among badger populations across western Europe is likely to be a consequence of this restricted local dispersal.