Andrew B. Feldman

Next-Generation Sequencing of Human Mitochondrial Reference Genomes Uncovers High Heteroplasmy Frequency

We describe methods for rapid sequencing of the entire human mitochondrial genome (mtgenome), which involve long-range PCR for specific amplification of the mtgenome, pyrosequencing, quantitative mapping of sequence reads to identify sequence variants and heteroplasmy, as well as de novo sequence assembly. These methods have been used to study 40 publicly available HapMap samples of European (CEU) and African (YRI) ancestry to demonstrate a sequencing error rate <5.63×10−4, nucleotide diversity of 1.6×10−3 for CEU and 3.7×10−3 for YRI, patterns of sequence variation consistent with earlier studies, but a higher rate of heteroplasmy varying between 10% and 50%. These results demonstrate that next-generation sequencing technologies allow interrogation of the mitochondrial genome in greater depth than previously possible which may be of value in biology and medicine.

A single equivalent moving dipole model: an efficient approach for localizing sites of origin of ventricular electrical activation.

AbstractWe propose a new method for guiding catheter ablation procedures to abolish sites of origin of arrhythmias. This method models both cardiac electrical activity and current pulses delivered from the tip of the ablation catheter with a single equivalent moving dipole (SEMD). The SEMD parameters are obtained from analysis of body surface potentials. In this paper we examine the feasibility of this method by evaluating the performance of an inverse algorithm we developed to localize the SEMD from the surface potentials. In computer simulations realistic levels of measurement noise led to uncertainties in SEMD location ∼0.005 cm. Dipole orientation randomization contributed to increased uncertainty (0.04 cm) in SEMD location only when boundary effects were included. In ventricular pacing swine studies, we found that the SEMD model accurately accounted for electrocardiographic wave forms and that measurement noise led to an uncertainty of approximately 0.04 cm in the SEMD at 15 ms after the pacing spike. We have also found that the algorithm we developed to identify the SEMD parameters yielded positions for two spatially separated pacing sites that maintained their direction and were very close to their physical separation. These results suggest that the SEMD method may potentially be used to guide radio-frequency ablation procedures.

Genomic Signatures of Strain Selection and Enhancement in Bacillus atrophaeus var. globigii ,a Historical Biowarfare Simulant

Background Despite the decades-long use of Bacillus atrophaeus var. globigii (BG) as a simulant for biological warfare (BW) agents, knowledge of its genome composition is limited. Furthermore, the ability to differentiate signatures of deliberate adaptation and selection from natural variation is lacking for most bacterial agents. We characterized a lineage of BGwith a long history of use as a simulant for BW operations, focusing on classical bacteriological markers, metabolic profiling and whole-genome shotgun sequencing (WGS). Results Archival strains and two “present day” type strains were compared to simulant strains on different laboratory media. Several of the samples produced multiple colony morphotypes that differed from that of an archival isolate. To trace the microevolutionary history of these isolates, we obtained WGS data for several archival and present-day strains and morphotypes. Bacillus-wide phylogenetic analysis identified B. subtilis as the nearest neighbor to B. atrophaeus. The genome of B. atrophaeus is, on average, 86% identical to B. subtilis on the nucleotide level. WGS of variants revealed that several strains were mixed but highly related populations and uncovered a progressive accumulation of mutations among the “military” isolates. Metabolic profiling and microscopic examination of bacterial cultures revealed enhanced growth of “military” isolates on lactate-containing media, and showed that the “military” strains exhibited a hypersporulating phenotype. Conclusions Our analysis revealed the genomic and phenotypic signatures of strain adaptation and deliberate selection for traits that were desirable in a simulant organism. Together, these results demonstrate the power of whole-genome and modern systems-level approaches to characterize microbial lineages to develop and validate forensic markers for strain discrimination and reveal signatures of deliberate adaptation.

Establishing Drug Resistance in Microorganisms by Mass Spectrometry

AbstractA rapid method to determine drug resistance in bacteria based on mass spectrometry is presented. In it, a mass spectrum of an intact microorganism grown in drug-containing stable isotope-labeled media is compared with a mass spectrum of the intact microorganism grown in non-labeled media without the drug present. Drug resistance is determined by predicting characteristic mass shifts of one or more microorganism biomarkers using bioinformatics algorithms. Observing such characteristic mass shifts indicates that the microorganism is viable even in the presence of the drug, thus incorporating the isotopic label into characteristic biomarker molecules. The performance of the method is illustrated on the example of intact E. coli, grown in control (unlabeled) and 13C-labeled media, and analyzed by MALDI TOF MS. Algorithms for data analysis are presented as well. Figureᅟ

Flexibility in surface-exposed loops in a virus capsid mediates escape from antibody neutralization.

ABSTRACT New human norovirus strains emerge every 2 to 3 years, partly due to mutations in the viral capsid that allow escape from antibody neutralization and herd immunity. To understand how noroviruses evolve antibody resistance, we investigated the structural basis for the escape of murine norovirus (MNV) from antibody neutralization. To identify specific residues in the MNV-1 protruding (P) domain of the capsid that play a role in escape from the neutralizing monoclonal antibody (MAb) A6.2, 22 recombinant MNVs were generated with amino acid substitutions in the A′B′ and E′F′ loops. Six mutations in the E′F′ loop (V378F, A382K, A382P, A382R, D385G, and L386F) mediated escape from MAb A6.2 neutralization. To elucidate underlying structural mechanisms for these results, the atomic structure of the A6.2 Fab was determined and fitted into the previously generated pseudoatomic model of the A6.2 Fab/MNV-1 virion complex. Previously, two distinct conformations, A and B, of the atomic structures of the MNV-1 P domain were identified due to flexibility in the two P domain loops. A superior stereochemical fit of the A6.2 Fab to the A conformation of the MNV P domain was observed. Structural analysis of our observed escape mutants indicates changes toward the less-preferred B conformation of the P domain. The shift in the structural equilibrium of the P domain toward the conformation with poor structural complementarity to the antibody strongly supports a unique mechanism for antibody escape that occurs via antigen flexibility instead of direct antibody-antigen binding. IMPORTANCE Human noroviruses cause the majority of all nonbacterial gastroenteritis worldwide. New epidemic strains arise in part by mutations in the viral capsid leading to escape from antibody neutralization. Herein, we identify a series of point mutations in a norovirus capsid that mediate escape from antibody neutralization and determine the structure of a neutralizing antibody. Fitting of the antibody structure into the virion/antibody complex identifies two conformations of the antibody binding domain of the viral capsid: one with a superior fit and the other with an inferior fit to the antibody. These data suggest a unique mode of antibody neutralization. In contrast to other viruses that largely escape antibody neutralization through direct disruption of the antibody-virus interface, we identify mutations that acted indirectly by limiting the conformation of the antibody binding loop in the viral capsid and drive the antibody binding domain into the conformation unable to be bound by the antibody.

Statistical accuracy of a moving equivalent dipole method to identify sites of origin of cardiac electrical activation

While radio frequency (RF) catheter ablation (RCA) procedures for treating ventricular arrhythmias have evolved significantly over the past several years, the use of RCA has been limited to treating slow ventricular tachycardias (VTs). In this paper, we present preliminary results from computer and animal studies to evaluate the accuracy of an algorithm that uses the single equivalent moving dipole (SEMD) model in an infinite homogeneous volume conductor to guide the RF catheter to the site of origin of the arrhythmia. Our method involves measuring body surface electrocardiographic (ECG) signals generated by arrhythmic activity and by bipolar current pulses emanating from a catheter tip, and representing each of them by a SEMD model source at each instant of the cardiac cycle, thus enabling rapid repositioning of the catheter tip requiring only a few cycles of the arrhythmia. We found that the SEMD model accurately reproduced body surface ECG signals with a correlation coefficients >0.95. We used a variety of methods to estimate the uncertainty of the SEMD parameters due to measurement noise and found that at the time when the arrhythmia is mostly localized during the cardiac cycle, the estimates of the uncertainty of the spatial SEMD parameters (from ECG signals) are between 1 and 3 mm. We used pacing data from spatially separated epicardial sites in a swine model as surrogates for focal ventricular arrhythmic sources and found that the spatial SEMD estimates of the two pacing sites agreed with both their physical separation and orientation with respect to each other. In conclusion, our algorithm to estimate the SEMD parameters from body surface ECG can potentially be a useful method for rapidly positioning the catheter tip to the arrhythmic focus during an RCA procedure.

Rapid, targeted and culture-free viral infectivity assay in drop-based microfluidics

A key viral property is infectivity, and its accurate measurement is crucial for the understanding of viral evolution, disease and treatment. Currently viral infectivity is measured using plaque assays, which involve prolonged culturing of host cells, and whose measurement is unable to differentiate between specific strains and is prone to low number fluctuation. We developed a rapid, targeted and culture-free infectivity assay using high-throughput drop-based microfluidics. Single infectious viruses are incubated in a large number of picoliter drops with host cells for one viral replication cycle followed by in-drop gene-specific amplification to detect infection events. Using murine noroviruses (MNV) as a model system, we measure their infectivity and determine the efficacy of a neutralizing antibody for different variants of MNV. Our results are comparable to traditional plaque-based assays and plaque reduction neutralization tests. However, the fast, low-cost, highly accurate genomic-based assay promises to be a superior method for drug screening and isolation of resistant viral strains. Moreover our technique can be adapted to measuring the infectivity of other pathogens, such as bacteria and fungi.

Isolation and Analysis of Rare Norovirus Recombinants from Coinfected Mice Using Drop-Based Microfluidics

ABSTRACT Human noroviruses (HuNoVs) are positive-sense RNA viruses that can cause severe, highly infectious gastroenteritis. HuNoV outbreaks are frequently associated with recombination between circulating strains. Strain genotyping and phylogenetic analyses show that noroviruses often recombine in a highly conserved region near the junction of the viral polyprotein (open reading frame 1 [ORF1]) and capsid (ORF2) genes and occasionally within the RNA-dependent RNA polymerase (RdRP) gene. Although genotyping methods are useful for tracking changes in circulating viral populations, they report only the dominant recombinant strains and do not elucidate the frequency or range of recombination events. Furthermore, the relatively low frequency of recombination in RNA viruses has limited studies to cell culture or in vitro systems, which do not reflect the complexities and selective pressures present in an infected organism. Using two murine norovirus (MNV) strains to model coinfection, we developed a microfluidic platform to amplify, detect, and recover individual recombinants following in vitro and in vivo coinfection. One-step reverse transcriptase PCR (RT-PCR) was performed in picoliter drops with primers that identified the wild-type and recombinant progenies and scanned for recombination breakpoints at ∼1-kb intervals. We detected recombination between MNV strains at multiple loci spanning the viral protease, RdRP, and capsid ORFs and isolated individual recombinant RNA genomes that were present at a frequency of 1/300,000 or higher. This study is the first to examine norovirus recombination following coinfection of an animal and suggests that the exchange of RNA among viral genomes in an infected host occurs in multiple locations and is an important driver of genetic diversity. IMPORTANCE RNA viruses increase diversity and escape host immune barriers by genomic recombination. Studies using a number of viral systems indicate that recombination occurs via template switching by the virus-encoded RNA-dependent RNA polymerase (RdRP). However, factors that govern the frequency and positions of recombination in an infected organism remain largely unknown. This work leverages advances in the applied physics of drop-based microfluidics to isolate and sequence rare recombinants arising from the coinfection of mice with two distinct strains of murine norovirus. This study is the first to detect and analyze norovirus recombination in an animal model.

Applicability of the single equivalent point dipole model to represent a spatially distributed bio-electrical source

Although the single equivalent point dipole model has been used to represent well-localised bio-electrical sources, in realistic situations the source is distributed. Consequently, position estimates of point dipoles determined by inverse algorithms suffer from systematic error due to the non-exact applicability of the inverse model. In realistic situations, this systematic error cannot be avoided, a limitation that is independent of the complexity of the torso model used. This study quantitatively investigates the intrinsic limitations in the assignment of a location to the equivalent dipole due to distributed electrical source. To simulate arrhythmic activity in the heart, a model of a wave of depolarisation spreading from a focal source over the surface of a spherical shell is used. The activity is represented by a sequence of concentric belt sources (obtained by slicing the shell with a sequence of parallel plane pairs), with constant dipole moment per unit length (circumferentially) directed parallel to the propagation direction. The distributed source is represented by N dipoles at equal arc lengths along the belt. The sum of the dipole potentials is calculated at predefined electrode locations. The inverse problem involves finding a single equivalent point dipole that best reproduces the electrode potentials due to the distributed source. The inverse problem is implemented by minimising the χ2 per degree of freedom. It is found that the trajectory traced by the equivalent dipole is sensitive to the location of the spherical shell relative to the fixed electrodes. It is shown that this trajectory does not coincide with the sequence of geometrical centres of the consecutive belt sources. For distributed sources within a bounded spherical medium, displaced from the spheres centre by 40% of the spheres radius, it is found that the error in the equivalent dipole location varies from 3 to 20% for sources with size between 5 and 50% of the spheres radius. Finally, a method is devised to obtain the size of the distributed source during the cardiac cycle.

Artifact-Free Quantification and Sequencing of Rare Recombinant Viruses by Using Drop-Based Microfluidics.

Recombination is an important driver in the evolution of viruses and thus is key to understanding viral epidemics and improving strategies to prevent future outbreaks. Characterization of rare recombinant subpopulations remains technically challenging because of artifacts such as artificial recombinants, known as chimeras, and amplification bias. To overcome this, we have developed a high‐throughput microfluidic technique with a second verification step in order to amplify and sequence single recombinant viruses with high fidelity in picoliter drops. We obtained the first artifact‐free estimate of in vitro recombination rate between murine norovirus strains MNV‐1 and WU20 co‐infecting a cell (Prec=3.3×10−4±2×10−5) for a 1205 nt region. Our approach represents a time‐ and cost‐effective improvement over current methods, and can be adapted for genomic studies requiring artifact‐ and bias‐free selective amplification, such as microbial pathogens, or rare cancer cells.