VDJServer: A Cloud-Based Analysis Portal and Data Commons for Immune Repertoire Sequences and Rearrangements.

Background: Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation. Results: VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity. Conclusion: VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.

VDJPipe: a pipelined tool for pre-processing immune repertoire sequencing data.

BACKGROUND: Pre-processing of high-throughput sequencing data for immune repertoire profiling is essential to insure high quality input for downstream analysis. VDJPipe is a flexible, high-performance tool that can perform multiple pre-processing tasks with just a single pass over the data files. RESULTS: Processing tasks provided by VDJPipe include base composition statistics calculation, read quality statistics calculation, quality filtering, homopolymer filtering, length and nucleotide filtering, paired-read merging, barcode demultiplexing, 5' and 3' PCR primer matching, and duplicate reads collapsing. VDJPipe utilizes a pipeline approach whereby multiple processing steps are performed in a sequential workflow, with the output of each step passed as input to the next step automatically. The workflow is flexible enough to handle the complex barcoding schemes used in many immunosequencing experiments. Because VDJPipe is designed for computational efficiency, we evaluated this by comparing execution times with those of pRESTO, a widely-used pre-processing tool for immune repertoire sequencing data. We found that VDJPipe requires <10% of the run time required by pRESTO. CONCLUSIONS: VDJPipe is a high-performance tool that is optimized for pre-processing large immune repertoire sequencing data sets.

VDJML: a file format with tools for capturing the results of inferring immune receptor rearrangements.

BACKGROUND: The genes that produce antibodies and the immune receptors expressed on lymphocytes are not germline encoded; rather, they are somatically generated in each developing lymphocyte by a process called V(D)J recombination, which assembles specific, independent gene segments into mature composite genes. The full set of composite genes in an individual at a single point in time is referred to as the immune repertoire. V(D)J recombination is the distinguishing feature of adaptive immunity and enables effective immune responses against an essentially infinite array of antigens. Characterization of immune repertoires is critical in both basic research and clinical contexts. Recent technological advances in repertoire profiling via high-throughput sequencing have resulted in an explosion of research activity in the field. This has been accompanied by a proliferation of software tools for analysis of repertoire sequencing data. Despite the widespread use of immune repertoire profiling and analysis software, there is currently no standardized format for output files from V(D)J analysis. Researchers utilize software such as IgBLAST and IMGT/High V-QUEST to perform V(D)J analysis and infer the structure of germline rearrangements. However, each of these software tools produces results in a different file format, and can annotate the same result using different labels. These differences make it challenging for users to perform additional downstream analyses. RESULTS: To help address this problem, we propose a standardized file format for representing V(D)J analysis results. The proposed format, VDJML, provides a common standardized format for different V(D)J analysis applications to facilitate downstream processing of the results in an application-agnostic manner. The VDJML file format specification is accompanied by a support library, written in C++ and Python, for reading and writing the VDJML file format. CONCLUSIONS: The VDJML suite will allow users to streamline their V(D)J analysis and facilitate the sharing of scientific knowledge within the community. The VDJML suite and documentation are available from https://vdjserver.org/vdjml/ . We welcome participation from the community in developing the file format standard, as well as code contributions.

owlcpp: a C++ library for working with OWL ontologies.

BACKGROUND: The increasing use of ontologies highlights the need for a library for working with ontologies that is efficient, accessible from various programming languages, and compatible with common computational platforms. RESULTS: We developed owlcpp, a library for storing and searching RDF triples, parsing RDF/XML documents, converting triples into OWL axioms, and reasoning. The library is written in ISO-compliant C++ to facilitate efficiency, portability, and accessibility from other programming languages. Internally, owlcpp uses the Raptor RDF Syntax library for parsing RDF/XML and the FaCT++ library for reasoning. The current version of owlcpp is supported under Linux, OSX, and Windows platforms and provides an API for Python. CONCLUSIONS: The results of our evaluation show that, compared to other commonly used libraries, owlcpp is significantly more efficient in terms of memory usage and searching RDF triple stores. owlcpp performs strict parsing and detects errors ignored by other libraries, thus reducing the possibility of incorrect semantic interpretation of ontologies. owlcpp is available at http://owl-cpp.sf.net/ under the Boost Software License, Version 1.0.

MSPrecise: A molecular diagnostic test for multiple sclerosis using next generation sequencing.

BACKGROUND: We have previously demonstrated that cerebrospinal fluid-derived B cells from early relapsing-remitting multiple sclerosis (RRMS) patients that express a VH4 gene accumulate specific replacement mutations. These mutations can be quantified as a score that identifies such patients as having or likely to convert to RRMS. Furthermore, we showed that next generation sequencing is an efficient method for obtaining the sequencing information required by this mutation scoring tool, originally developed using the less clinically viable single-cell Sanger sequencing. OBJECTIVE: To determine the accuracy of MSPrecise, the diagnostic test that identifies the presence of the RRMS-enriched mutation pattern from patient cerebrospinal fluid B cells. METHODS: Cerebrospinal fluid cell pellets were obtained from RRMS and other neurological disease (OND) patient cohorts. VH4 gene segments were amplified, sequenced by next generation sequencing and analyzed for mutation score. RESULTS: The diagnostic test showed a sensitivity of 75% on the RRMS cohort and a specificity of 88% on the OND cohort. The accuracy of the test in identifying RRMS patients or patients that will develop RRMS is 84%. CONCLUSION: MSPrecise exhibits good performance in identifying patients with RRMS irrespective of time with RRMS.

The antibody genetics of multiple sclerosis: comparing next-generation sequencing to sanger sequencing.

We previously identified a distinct mutation pattern in the antibody genes of B cells isolated from cerebrospinal fluid (CSF) that can identify patients who have relapsing-remitting multiple sclerosis (RRMS) and patients with clinically isolated syndromes who will convert to RRMS. This antibody gene signature (AGS) was developed using Sanger sequencing of single B cells. While potentially helpful to patients, Sanger sequencing is not an assay that can be practically deployed in clinical settings. In order to provide AGS evaluations to patients as part of their diagnostic workup, we developed protocols to generate AGS scores using next-generation DNA sequencing (NGS) on CSF-derived cell pellets without the need to isolate single cells. This approach has the potential to increase the coverage of the B-cell population being analyzed, reduce the time needed to generate AGS scores, and may improve the overall performance of the AGS approach as a diagnostic test in the future. However, no investigations have focused on whether NGS-based repertoires will properly reflect antibody gene frequencies and somatic hypermutation patterns defined by Sanger sequencing. To address this issue, we isolated paired CSF samples from eight patients who either had MS or were at risk to develop MS. Here, we present data that antibody gene frequencies and somatic hypermutation patterns are similar in Sanger and NGS-based antibody repertoires from these paired CSF samples. In addition, AGS scores derived from the NGS database correctly identified the patients who initially had or subsequently converted to RRMS, with precision similar to that of the Sanger sequencing approach. Further investigation of the utility of the AGS in predicting conversion to MS using NGS-derived antibody repertoires in a larger cohort of patients is warranted.

Haplotype Association Mapping Identifies a Candidate Gene Region in Mice Infected With Staphylococcus aureus.

Exposure to Staphylococcus aureus has a variety of outcomes, from asymptomatic colonization to fatal infection. Strong evidence suggests that host genetics play an important role in susceptibility, but the specific host genetic factors involved are not known. The availability of genome-wide single nucleotide polymorphism (SNP) data for inbred Mus musculus strains means that haplotype association mapping can be used to identify candidate susceptibility genes. We applied haplotype association mapping to Perlegen SNP data and kidney bacterial counts from Staphylococcus aureus-infected mice from 13 inbred strains and detected an associated block on chromosome 7. Strong experimental evidence supports the result: a separate study demonstrated the presence of a susceptibility locus on chromosome 7 using consomic mice. The associated block contains no genes, but lies within the gene cluster of the 26-member extended kallikrein gene family, whose members have well-recognized roles in the generation of antimicrobial peptides and the regulation of inflammation. Efficient mixed-model association (EMMA) testing of all SNPs with two alleles and located within the gene cluster boundaries finds two significant associations: one of the three polymorphisms defining the associated block and one in the gene closest to the block, Klk1b11. In addition, we find that 7 of the 26 kallikrein genes are differentially expressed between susceptible and resistant mice, including the Klk1b11 gene. These genes represent a promising set of candidate genes influencing susceptibility to Staphylococcus aureus.

The protease domain increases the translocation stepping efficiency of the hepatitis C virus NS3-4A helicase.

Hepatitis C virus (HCV) NS3 protein has two enzymatic activities of helicase and protease that are essential for viral replication. The helicase separates the strands of DNA and RNA duplexes using the energy from ATP hydrolysis. To understand how ATP hydrolysis is coupled to helicase movement, we measured the single turnover helicase translocation-dissociation kinetics and the pre-steady-state P(i) release kinetics on single-stranded RNA and DNA substrates of different lengths. The parameters of stepping were determined from global fitting of the two types of kinetic measurements into a computational model that describes translocation as a sequence of coupled hydrolysis-stepping reactions. Our results show that the HCV helicase moves with a faster rate on single stranded RNA than on DNA. The HCV helicase steps on the RNA or DNA one nucleotide at a time, and due to imperfect coupling, not every ATP hydrolysis event produces a successful step. Comparison of the helicase domain (NS3h) with the protease-helicase (NS3-4A) shows that the most significant contribution of the protease domain is to improve the translocation stepping efficiency of the helicase. Whereas for NS3h, only 20% of the hydrolysis events result in translocation, the coupling for NS3-4A is near-perfect 93%. The presence of the protease domain also significantly reduces the stepping rate, but it doubles the processivity. These effects of the protease domain on the helicase can be explained by an improved allosteric cross-talk between the ATP- and nucleic acid-binding sites achieved by the overall stabilization of the helicase domain structure.

Experimental and computational analysis of DNA unwinding and polymerization kinetics.

DNA unwinding and polymerization are complex processes involving many intermediate species in the reactions. Our understanding of these processes is limited because the rates of the reactions or the existence of intermediate species is not apparent without specially designed experimental techniques and data analysis procedures. In this chapter we describe how pre-steady state and single-turnover measurements analyzed by model-based methods can be used for estimating the elementary rate constants. Using the hexameric helicase and the DNA polymerase from bacteriophage T7 as model systems, we provide stepwise procedures for measuring the kinetics of the reactions they catalyze based on radioactivity and fluorescence. We also describe analysis of the experimental measurements using publicly available models and software gfit ( http://gfit.sf.net ).

Model-based global analysis of heterogeneous experimental data using gfit.

Regression analysis is indispensible for quantitative understanding of biological systems and for developing accurate computational models. By applying regression analysis, one can validate models and quantify components of the system, including ones that cannot be observed directly. Global (simultaneous) analysis of all experimental data available for the system produces the most informative results. To quantify components of a complex system, the dataset needs to contain experiments of different types performed under a broad range of conditions. However, heterogeneity of such datasets complicates implementation of the global analysis. Computational models continuously evolve to include new knowledge and to account for novel experimental data, creating the demand for flexible and efficient analysis procedures. To address these problems, we have developed gfit software to globally analyze many types of experiments, to validate computational models, and to extract maximum information from the available experimental data.

Mechanism of ATP-driven PCNA clamp loading by S. cerevisiae RFC.

Circular clamps tether polymerases to DNA, serving as essential processivity factors in genome replication, and function in other critical cellular processes as well. Clamp loaders catalyze clamp assembly onto DNA, and the question of how these proteins construct a topological link between a clamp and DNA, especially the mechanism by which ATP is utilized for the task, remains open. Here we describe pre-steady-state analysis of ATP hydrolysis, proliferating cell nuclear antigen (PCNA) clamp opening, and DNA binding by Saccharomyces cerevisiae replication factor C (RFC), and present the first kinetic model of a eukaryotic clamp-loading reaction validated by global data analysis. ATP binding to multiple RFC subunits initiates a slow conformational change in the clamp loader, enabling it to bind and open PCNA and to bind DNA as well. PCNA opening locks RFC into an active state, and the resulting RFC.ATP.PCNA((open)) intermediate is ready for the entry of DNA into the clamp. DNA binding commits RFC to ATP hydrolysis, which is followed by PCNA closure and PCNA.DNA release. This model enables quantitative understanding of the multistep mechanism of a eukaryotic clamp loader and furthermore facilitates comparative analysis of loaders from diverse organisms.

Mechanism of cadmium-mediated inhibition of Msh2-Msh6 function in DNA mismatch repair.

The observation that Cadmium (Cd(2+)) inhibits Msh2-Msh6, which is responsible for identifying base pair mismatches and other discrepancies in DNA, has led to the proposal that selective targeting of this protein and consequent suppression of DNA repair or apoptosis promote the carcinogenic effects of the heavy metal toxin. It has been suggested that Cd(2+) binding to specific sites on Msh2-Msh6 blocks its DNA binding and ATPase activities. To investigate the mechanism of inhibition, we measured Cd(2+) binding to Msh2-Msh6, directly and by monitoring changes in protein structure and enzymatic activity. Global fitting of the data to a multiligand binding model revealed that binding of about 100 Cd(2+) ions per Msh2-Msh6 results in its inactivation. This finding indicates that the inhibitory effect of Cd(2+) occurs via a nonspecific mechanism. Cd(2+) and Msh2-Msh6 interactions involve cysteine sulfhydryl groups, and the high Cd(2+):Msh2-Msh6 ratio implicates other ligands such as histidine, aspartate, glutamate, and the peptide backbone as well. Our study also shows that cadmium inactivates several unrelated enzymes similarly, consistent with a nonspecific mechanism of inhibition. Targeting of a variety of proteins, including Msh2-Msh6, in this generic manner would explain the marked broad-spectrum impact of Cd(2+) on biological processes. We propose that the presence of multiple nonspecific Cd(2+) binding sites on proteins and their propensity to change conformation on interaction with Cd(2+) are critical determinants of the susceptibility of corresponding biological systems to cadmium toxicity.

Multiplexed RNA trafficking in oligodendrocytes and neurons.

In oligodendrocytes and neurons genetic information is transmitted from the nucleus to dendrites in the form of RNA granules. Here we describe how transport of multiple different RNA molecules in individual granules is analogous to the process of multiplexing in telecommunications. In both cases multiple messages are combined into a composite signal for transmission on a single carrier. Multiplexing provides a mechanism to coordinate local expression of ensembles of genes in myelin in oligodendrocytes and at synapses in neurons.

Multiplexed dendritic targeting of alpha calcium calmodulin-dependent protein kinase II, neurogranin, and activity-regulated cytoskeleton-associated protein RNAs by the A2 pathway.

In neurons, many different RNAs are targeted to dendrites where local expression of the encoded proteins mediates synaptic plasticity during learning and memory. It is not known whether each RNA follows a separate trafficking pathway or whether multiple RNAs are targeted to dendrites by the same pathway. Here, we show that RNAs encoding alpha calcium calmodulin-dependent protein kinase II, neurogranin, and activity-regulated cytoskeleton-associated protein are coassembled into the same RNA granules and targeted to dendrites by the same cis/trans-determinants (heterogeneous nuclear ribonucleoprotein [hnRNP] A2 response element and hnRNP A2) that mediate dendritic targeting of myelin basic protein RNA by the A2 pathway in oligodendrocytes. Multiplexed dendritic targeting of different RNAs by the same pathway represents a new organizing principle for coordinating gene expression at the synapse.

Significant proportions of nuclear transport proteins with reduced intracellular mobilities resolved by fluorescence correlation spectroscopy.

Nuclear transport requires freely diffusing nuclear transport proteins to facilitate movement of cargo molecules through the nuclear pore. We analyzed dynamic properties of importin alpha, importin beta, Ran and NTF2 in nucleus, cytoplasm and at the nuclear pore of neuroblastoma cells using fluorescence correlation spectroscopy. Mobile components were quantified by global fitting of autocorrelation data from multiple cells. Immobile components were quantified by analysis of photobleaching kinetics. Wild-type Ran was compared to various mutant Ran proteins to identify components representing GTP or GDP forms of Ran. Untreated cells were compared to cells treated with nocodazole or latrunculin to identify components associated with cytoskeletal elements. The results indicate that freely diffusing importin alpha, importin beta, Ran and NTF2 are in dynamic equilibrium with larger pools associated with immobile binding partners such as microtubules in the cytoplasm. These findings suggest that formation of freely diffusing nuclear transport intermediates is in competition with binding to immobile partners. Variation in concentrations of freely diffusing nuclear transport intermediates among cells indicates that the nuclear transport system is sufficiently robust to function over a wide range of conditions.

DNA synthesis provides the driving force to accelerate DNA unwinding by a helicase.

Helicases are molecular motors that use the energy of nucleoside 5'-triphosphate (NTP) hydrolysis to translocate along a nucleic acid strand and catalyse reactions such as DNA unwinding. The ring-shaped helicase of bacteriophage T7 translocates along single-stranded (ss)DNA at a speed of 130 bases per second; however, T7 helicase slows down nearly tenfold when unwinding the strands of duplex DNA. Here, we report that T7 DNA polymerase, which is unable to catalyse strand displacement DNA synthesis by itself, can increase the unwinding rate to 114 base pairs per second, bringing the helicase up to similar speeds compared to its translocation along ssDNA. The helicase rate of stimulation depends upon the DNA synthesis rate and does not rely on specific interactions between T7 DNA polymerase and the carboxy-terminal residues of T7 helicase. Efficient duplex DNA synthesis is achieved only by the combined action of the helicase and polymerase. The strand displacement DNA synthesis by the DNA polymerase depends on the unwinding activity of the helicase, which provides ssDNA template. The rapid trapping of the ssDNA bases by the DNA synthesis activity of the polymerase in turn drives the helicase to move forward through duplex DNA at speeds similar to those observed along ssDNA.

A Brownian motor mechanism of translocation and strand separation by hepatitis C virus helicase.

Helicases translocate along their nucleic acid substrates using the energy of ATP hydrolysis and by changing conformations of their nucleic acid-binding sites. Our goal is to characterize the conformational changes of hepatitis C virus (HCV) helicase at different stages of ATPase cycle and to determine how they lead to translocation. We have reported that ATP binding reduces HCV helicase affinity for nucleic acid. Now we identify the stage of the ATPase cycle responsible for translocation and unwinding. We show that a rapid directional movement occurs upon helicase binding to DNA in the absence of ATP, resulting in opening of several base pairs. We propose that HCV helicase translocates as a Brownian motor with a simple two-stroke cycle. The directional movement step is fueled by single-stranded DNA binding energy while ATP binding allows for a brief period of random movement that prepares the helicase for the next cycle.

The DNA-unwinding mechanism of the ring helicase of bacteriophage T7.

Helicases are motor proteins that use the chemical energy of NTP hydrolysis to drive mechanical processes such as translocation and nucleic acid strand separation. Bacteriophage T7 helicase functions as a hexameric ring to drive the replication complex by separating the DNA strands during genome replication. Our studies show that T7 helicase unwinds DNA with a low processivity, and the results indicate that the low processivity is due to ring opening and helicase dissociating from the DNA during unwinding. We have measured the single-turnover kinetics of DNA unwinding and globally fit the data to a modified stepping model to obtain the unwinding parameters. The comparison of the unwinding properties of T7 helicase with its translocation properties on single-stranded (ss)DNA has provided insights into the mechanism of strand separation that is likely to be general for ring helicases. T7 helicase unwinds DNA with a rate of 15 bp/s, which is 9-fold slower than the translocation speed along ssDNA. T7 helicase is therefore primarily an ssDNA translocase that does not directly destabilize duplex DNA. We propose that T7 helicase achieves DNA unwinding by its ability to bind ssDNA because it translocates unidirectionally, excluding the complementary strand from its central channel. The results also imply that T7 helicase by itself is not an efficient helicase and most likely becomes proficient at unwinding when it is engaged in a replication complex.

The functional interaction of the hepatitis C virus helicase molecules is responsible for unwinding processivity.

Although helicases participate in virtually every cellular process involving nucleic acids, the details of their mechanism including the role of interaction between the subunits remains unclear. Here we study the unwinding kinetics of the helicase from hepatitis C virus using DNA substrates with a range of tail and duplex lengths. The binding of the helicase to the substrates was characterized by electron microscopy and fluorimetric titrations. Depending on the length of the ssDNA tail, one or more helicase molecules can be loaded on the DNA. Unwinding was measured under single-turnover conditions, and the results show that a monomer is active on short duplexes yet multiple molecules are needed to unwind long duplexes. Thus, increasing the ssDNA tail length increases the unwinding efficiency. The unwinding kinetics was modeled as a stepwise process performed by single or multiple helicase molecules. The model programmed in MATLAB was used for global fitting of the kinetics, yielding values for the rate of unwinding, processivity, cooperativity, step size, and occlusion site. The results indicate that a single hepatitis C virus helicase molecule unwinds DNA with a low processivity. The multiple helicase molecules present on the DNA substrate show functional cooperativity and unwind with greater efficiency, although they bind and release the substrate non-cooperatively, and the ATPase cycle of the helicase molecules is not coordinated. The functional interaction model explains the efficient unwinding by multiple helicases and is generally applicable.

Fluorescence correlation spectroscopy and quantitative cell biology.

Fluorescence correlation spectroscopy (FCS) analyzes fluctuations in fluorescence within a small observation volume. Autocorrelation analysis of FCS fluctuation data can be used to measure concentrations, diffusion properties, and kinetic constants for individual fluorescent molecules. Photon count histogram analysis of fluorescence fluctuation data can be used to study oligomerization of individual fluorescent molecules. If the FCS observation volume is positioned inside a living cell, these parameters can be measured in vivo. FCS can provide the requisite quantitative data for analysis of molecular interaction networks underlying complex cell biological processes.