A hemifused complex is the hub in a network of pathways to membrane fusion.

Membrane fusion is a critical step for many essential processes, from neurotransmission to fertilization. For over 40 years, protein-free fusion driven by calcium or other cationic species has provided a simplified model of biological fusion, but the mechanisms remain poorly understood. Cation-mediated membrane fusion and permeation are essential in their own right to drug delivery strategies based on cell-penetrating peptides or cation-bearing lipid nanoparticles. Experimental studies suggest calcium drives anionic membranes to a hemifused intermediate that constitutes a hub in a network of pathways, but the pathway selection mechanism is unknown. Here we develop a mathematical model that identifies the network hub as a highly dynamic hemifusion complex. Multivalent cations drive expansion of this high-tension hemifusion interface between interacting vesicles during a brief transient. The fate of this interface determines the outcome, either fusion, dead-end hemifusion, or vesicle lysis. The model reproduces the unexplained finding that calcium-driven fusion of vesicles with planar membranes typically stalls at hemifusion, and we show the equilibrated hemifused state is a novel lens-shaped complex. Thus, membrane fusion kinetics follow a stochastic trajectory within a network of pathways, with outcome weightings set by a hemifused complex intermediate.

Model of influenza virus acidification.

Internal acidification of the influenza virus, mediated by the M2 proton channel, is a key step in its life cycle. The interior M1 protein shell dissolves at pH~5.5 to 6.0, allowing the release of vRNA to the cytoplasm upon fusion of the viral envelope with the endosomal membrane. Previous models have described the mechanisms and rate constants of M2-mediated transport but did not describe the kinetics of pH changes inside the virus or consider exterior pH changes due to endosome maturation. Therefore, we developed a mathematical model of M2-mediated virion acidification. We find that ~32,000 protons are required to acidify a typically-sized virion. Predicted acidification kinetics were consistent with published in vitro experiments following internal acidification. Finally, we applied the model to the in vivo situation. For all rates of endosomal maturation considered, internal acidification lagged ~1 min behind endosomal acidification to pH 6. For slow endosomal maturation requiring several minutes or more, internal and endosomal pH decay together in pseudo-equilibrium to the late endosomal pH~5.0. For fast endosomal maturation (≲2 min), a lag of tens of seconds continued toward the late endosomal pH. Recent experiments suggest in vivo maturation is in this "fast" regime where lag is considerable. We predict that internal pH reaches the threshold for M1 shell solvation just before the external pH triggers membrane fusion mediated by the influenza protein hemagglutinin, critical because outward proton diffusion through a single small fusion pore is faster than the collective M2-mediated transport inward.

Entropic forces drive self-organization and membrane fusion by SNARE proteins.

SNARE proteins are the core of the cell's fusion machinery and mediate virtually all known intracellular membrane fusion reactions on which exocytosis and trafficking depend. Fusion is catalyzed when vesicle-associated v-SNAREs form trans-SNARE complexes ("SNAREpins") with target membrane-associated t-SNAREs, a zippering-like process releasing ∼65 kT per SNAREpin. Fusion requires several SNAREpins, but how they cooperate is unknown and reports of the number required vary widely. To capture the collective behavior on the long timescales of fusion, we developed a highly coarse-grained model that retains key biophysical SNARE properties such as the zippering energy landscape and the surface charge distribution. In simulations the ∼65-kT zippering energy was almost entirely dissipated, with fully assembled SNARE motifs but uncomplexed linker domains. The SNAREpins self-organized into a circular cluster at the fusion site, driven by entropic forces that originate in steric-electrostatic interactions among SNAREpins and membranes. Cooperative entropic forces expanded the cluster and pulled the membranes together at the center point with high force. We find that there is no critical number of SNAREs required for fusion, but instead the fusion rate increases rapidly with the number of SNAREpins due to increasing entropic forces. We hypothesize that this principle finds physiological use to boost fusion rates to meet the demanding timescales of neurotransmission, exploiting the large number of v-SNAREs available in synaptic vesicles. Once in an unfettered cluster, we estimate ≥15 SNAREpins are required for fusion within the ∼1-ms timescale of neurotransmitter release.

Metagenomic analysis of medicinal Cannabis samples; pathogenic bacteria, toxigenic fungi, and beneficial microbes grow in culture-based yeast and mold tests.

Background: The presence of bacteria and fungi in medicinal or recreational Cannabis poses a potential threat to consumers if those microbes include pathogenic or toxigenic species. This study evaluated two widely used culture-based platforms for total yeast and mold (TYM) testing marketed by 3M Corporation and Biomérieux, in comparison with a quantitative PCR (qPCR) approach marketed by Medicinal Genomics Corporation. Methods: A set of 15 medicinal Cannabis samples were analyzed using 3M and Biomérieux culture-based platforms and by qPCR to quantify microbial DNA. All samples were then subjected to next-generation sequencing and metagenomics analysis to enumerate the bacteria and fungi present before and after growth on culture-based media. Results: Several pathogenic or toxigenic bacterial and fungal species were identified in proportions of >5% of classified reads on the samples, including Acinetobacter baumannii, Escherichia coli, Pseudomonas aeruginosa, Ralstonia pickettii, Salmonella enterica, Stenotrophomonas maltophilia, Aspergillus ostianus, Aspergillus sydowii, Penicillium citrinum and Penicillium steckii. Samples subjected to culture showed substantial shifts in the number and diversity of species present, including the failure of Aspergillus species to grow well on either platform. Substantial growth of Clostridium botulinum and other bacteria were frequently observed on one or both of the culture-based TYM platforms. The presence of plant growth promoting (beneficial) fungal species further influenced the differential growth of species in the microbiome of each sample. Conclusions: These findings have important implications for the Cannabis and food safety testing industries.

Cholesterol Increases the Openness of SNARE-Mediated Flickering Fusion Pores.

Flickering of fusion pores during exocytotic release of hormones and neurotransmitters is well documented, but without assays that use biochemically defined components and measure single-pore dynamics, the mechanisms remain poorly understood. We used total internal reflection fluorescence microscopy to quantify fusion-pore dynamics in vitro and to separate the roles of soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins and lipid bilayer properties. When small unilamellar vesicles bearing neuronal v-SNAREs fused with planar bilayers reconstituted with cognate t-SNARES, lipid and soluble cargo transfer rates were severely reduced, suggesting that pores flickered. From the lipid release times we computed pore openness, the fraction of time the pore is open, which increased dramatically with cholesterol. For most lipid compositions tested, SNARE-mediated and nonspecifically nucleated pores had similar openness, suggesting that pore flickering was controlled by lipid bilayer properties. However, with physiological cholesterol levels, SNAREs substantially increased the fraction of fully open pores and fusion was so accelerated that there was insufficient time to recruit t-SNAREs to the fusion site, consistent with t-SNAREs being preclustered by cholesterol into functional docking and fusion platforms. Our results suggest that cholesterol opens pores directly by reducing the fusion-pore bending energy, and indirectly by concentrating several SNAREs into individual fusion events.

Expanded genetic codes in next generation sequencing enable decontamination and mitochondrial enrichment.

We have developed a PCR method, coined Déjà vu PCR, that utilizes six nucleotides in PCR with two methyl specific restriction enzymes that respectively digest these additional nucleotides. Use of this enzyme-and-nucleotide combination enables what we term a "DNA diode", where DNA can advance in a laboratory in only one direction and cannot feedback into upstream assays. Here we describe aspects of this method that enable consecutive amplification with the introduction of a 5th and 6th base while simultaneously providing methylation dependent mitochondrial DNA enrichment. These additional nucleotides enable a novel DNA decontamination technique that generates ephemeral and easy to decontaminate DNA.

Evolution of the hemifused intermediate on the pathway to membrane fusion.

The pathway to membrane fusion in synthetic and biological systems is thought to pass through hemifusion, in which the outer leaflets are fused while the inner leaflets engage in a hemifusion diaphragm (HD). Fusion has been proposed to be completed by lysis of the expanded HD that matures from a localized stalklike initial connection. However, the process that establishes the expanded HD is poorly understood. Here we mathematically modeled hemifusion of synthetic vesicles, where hemifusion and fusion are most commonly driven by calcium and membrane tension. The model shows that evolution of the hemifused state is driven by these agents and resisted by interleaflet frictional and tensile stresses. Predicted HD growth rates depend on tension and salt concentration, and agree quantitatively with experimental measurements. For typical conditions, we predict that HDs expand at ~30 µm(2)/s, reaching a final equilibrium area ~7% of the vesicle area. Key model outputs are the evolving HD tension and area during the growth transient, properties that may determine whether HD lysis occurs. Applying the model to numerous published experimental studies that reported fusion, our results are consistent with a final fusion step in which the HD ruptures due to super-lysis HD membrane tensions.

Accurate variant detection across non-amplified and whole genome amplified DNA using targeted next generation sequencing.

BACKGROUND: Many hypothesis-driven genetic studies require the ability to comprehensively and efficiently target specific regions of the genome to detect sequence variations. Often, sample availability is limited requiring the use of whole genome amplification (WGA). We evaluated a high-throughput microdroplet-based PCR approach in combination with next generation sequencing (NGS) to target 384 discrete exons from 373 genes involved in cancer. In our evaluation, we compared the performance of six non-amplified gDNA samples from two HapMap family trios. Three of these samples were also preamplified by WGA and evaluated. We tested sample pooling or multiplexing strategies at different stages of the tested targeted NGS (T-NGS) workflow. RESULTS: The results demonstrated comparable sequence performance between non-amplified and preamplified samples and between different indexing strategies [sequence specificity of 66.0% ± 3.4%, uniformity (coverage at 0.2× of the mean) of 85.6% ± 0.6%]. The average genotype concordance maintained across all the samples was 99.5% ± 0.4%, regardless of sample type or pooling strategy. We did not detect any errors in the Mendelian patterns of inheritance of genotypes between the parents and offspring within each trio. We also demonstrated the ability to detect minor allele frequencies within the pooled samples that conform to predicted models. CONCLUSION: Our described PCR-based sample multiplex approach and the ability to use WGA material for NGS may enable researchers to perform deep resequencing studies and explore variants at very low frequencies and cost.

The hemifused state on the pathway to membrane fusion.

Fusion of compartments enclosed by membrane bilayers enables secretion and other vital cellular processes and is widely studied in model synthetic membrane systems. Experiments suggest the fusion pathway passes through a hemifused intermediate where only outer monolayers are fused. Here we show membrane tension and divalent cations drive vesicles to hemifused equilibrium with expanded hemifusion diaphragms (HDs) where inner monolayers engage. Predicted HD sizes agree with recent measurements of Nikolaus et al. [Biophys. J. 98, 1192 (2010).]. The fusion pathway is completed by HD lysis provided HD tension is sufficiently high.

Pipeline for large-scale microdroplet bisulfite PCR-based sequencing allows the tracking of hepitype evolution in tumors.

Cytosine methylation provides an epigenetic level of cellular plasticity that is important for development, differentiation and cancerogenesis. We adopted microdroplet PCR to bisulfite treated target DNA in combination with second generation sequencing to simultaneously assess DNA sequence and methylation. We show measurement of methylation status in a wide range of target sequences (total 34 kb) with an average coverage of 95% (median 100%) and good correlation to the opposite strand (rho = 0.96) and to pyrosequencing (rho = 0.87). Data from lymphoma and colorectal cancer samples for SNRPN (imprinted gene), FGF6 (demethylated in the cancer samples) and HS3ST2 (methylated in the cancer samples) serve as a proof of principle showing the integration of SNP data and phased DNA-methylation information into "hepitypes" and thus the analysis of DNA methylation phylogeny in the somatic evolution of cancer.

Application of microdroplet PCR for large-scale targeted bisulfite sequencing.

Cytosine methylation of DNA CpG dinucleotides in gene promoters is an epigenetic modification that regulates gene transcription. While many methods exist to interrogate methylation states, few current methods offer large-scale, targeted, single CpG resolution. We report an approach combining bisulfite treatment followed by microdroplet PCR with next-generation sequencing to assay the methylation state of 50 genes in the regions 1 kb upstream of and downstream from their transcription start sites. This method yielded 96% coverage of the targeted CpGs and demonstrated high correlation between CpG island (CGI) DNA methylation and transcriptional regulation. The method was scaled to interrogate the methylation status of 77,674 CpGs in the promoter regions of 2100 genes in primary CD4 T cells. The 2100 gene library yielded 97% coverage of all targeted CpGs and 99% of the target amplicons.

Management of open apex in a central incisor using mineral trioxide aggregate.

UNLABELLED: Trauma to an immature tooth can result in pulpal devitalization and arrested apexogenesis, resulting in open apices; this is most common in anterior teeth of the permanent dentition. The attainment of an apical seal in such cases is challenging. The management of an open apex in a central incisor in a 16-year-old male patient is described in this report. The satisfactory result achieved suggests that mineral trioxide aggregate (MTA) can be used successfully as a filling material in teeth with open apices. CLINICAL RELEVANCE: Preservation of anterior teeth following trauma is essential, especially in young patients. A good option is the use of mineral trioxide aggregate in the treatment of open apices.

Microdroplet-based PCR enrichment for large-scale targeted sequencing.

Targeted enrichment of specific loci of the human genome is a promising approach to enable sequencing-based studies of genetic variation in large populations. Here we describe an enrichment approach based on microdroplet PCR, which enables 1.5 million amplifications in parallel. We sequenced six samples enriched by microdroplet or traditional singleplex PCR using primers targeting 435 exons of 47 genes. Both methods generated similarly high-quality data: 84% of the uniquely mapping reads fell within the targeted sequences; coverage was uniform across approximately 90% of targeted bases; sequence variants were called with >99% accuracy; and reproducibility between samples was high (r(2) = 0.9). We scaled the microdroplet PCR to 3,976 amplicons totaling 1.49 Mb of sequence, sequenced the resulting sample with both Illumina GAII and Roche 454, and obtained data with equally high specificity and sensitivity. Our results demonstrate that microdroplet technology is well suited for processing DNA for massively parallel enrichment of specific subsets of the human genome for targeted sequencing.

Model of SNARE-mediated membrane adhesion kinetics.

SNARE proteins are conserved components of the core fusion machinery driving diverse membrane adhesion and fusion processes in the cell. In many cases micron-sized membranes adhere over large areas before fusion. Reconstituted in vitro assays have helped isolate SNARE mechanisms in small membrane adhesion-fusion and are emerging as powerful tools to study large membrane systems by use of giant unilamellar vesicles (GUVs). Here we model SNARE-mediated adhesion kinetics in SNARE-reconstituted GUV-GUV or GUV-supported bilayer experiments. Adhesion involves many SNAREs whose complexation pulls apposing membranes into contact. The contact region is a tightly bound rapidly expanding patch whose growth velocity v(patch) increases with SNARE density Gamma(snare). We find three patch expansion regimes: slow, intermediate, fast. Typical experiments belong to the fast regime where v(patch) ~ (Gamma(snare)(2/3) depends on SNARE diffusivities and complexation binding constant. The model predicts growth velocities ~10 - 300 microm/s. The patch may provide a close contact region where SNAREs can trigger fusion. Extending the model to a simple description of fusion, a broad distribution of fusion times is predicted. Increasing SNARE density accelerates fusion by boosting the patch growth velocity, thereby providing more complexes to participate in fusion. This quantifies the notion of SNAREs as dual adhesion-fusion agents.

Short-term calorie restriction in male mice feminizes gene expression and alters key regulators of conserved aging regulatory pathways.

BACKGROUND: Calorie restriction (CR) is the only intervention known to extend lifespan in a wide range of organisms, including mammals. However, the mechanisms by which it regulates mammalian aging remain largely unknown, and the involvement of the TOR and sirtuin pathways (which regulate aging in simpler organisms) remain controversial. Additionally, females of most mammals appear to live longer than males within species; and, although it remains unclear whether this holds true for mice, the relationship between sex-biased and CR-induced gene expression remains largely unexplored. METHODOLOGY/PRINCIPAL FINDINGS: We generated microarray gene expression data from livers of male mice fed high calorie or CR diets, and we find that CR significantly changes the expression of over 3,000 genes, many between 10- and 50-fold. We compare our data to the GenAge database of known aging-related genes and to prior microarray expression data of genes expressed differently between male and female mice. CR generally feminizes gene expression and many of the most significantly changed individual genes are involved in aging, hormone signaling, and p53-associated regulation of the cell cycle and apoptosis. Among the genes showing the largest and most statistically significant CR-induced expression differences are Ddit4, a key regulator of the TOR pathway, and Nnmt, a regulator of lifespan linked to the sirtuin pathway. Using western analysis we confirmed post-translational inhibition of the TOR pathway. CONCLUSIONS: Our data show that CR induces widespread gene expression changes and acts through highly evolutionarily conserved pathways, from microorganisms to mammals, and that its life-extension effects might arise partly from a shift toward a gene expression profile more typical of females.

Rapid whole-genome mutational profiling using next-generation sequencing technologies.

Forward genetic mutational studies, adaptive evolution, and phenotypic screening are powerful tools for creating new variant organisms with desirable traits. However, mutations generated in the process cannot be easily identified with traditional genetic tools. We show that new high-throughput, massively parallel sequencing technologies can completely and accurately characterize a mutant genome relative to a previously sequenced parental (reference) strain. We studied a mutant strain of Pichia stipitis, a yeast capable of converting xylose to ethanol. This unusually efficient mutant strain was developed through repeated rounds of chemical mutagenesis, strain selection, transformation, and genetic manipulation over a period of seven years. We resequenced this strain on three different sequencing platforms. Surprisingly, we found fewer than a dozen mutations in open reading frames. All three sequencing technologies were able to identify each single nucleotide mutation given at least 10-15-fold nominal sequence coverage. Our results show that detecting mutations in evolved and engineered organisms is rapid and cost-effective at the whole-genome level using new sequencing technologies. Identification of specific mutations in strains with altered phenotypes will add insight into specific gene functions and guide further metabolic engineering efforts.

Systematic identification of mammalian regulatory motifs' target genes and functions.

We developed an algorithm, Lever, that systematically maps metazoan DNA regulatory motifs or motif combinations to sets of genes. Lever assesses whether the motifs are enriched in cis-regulatory modules (CRMs), predicted by our PhylCRM algorithm, in the noncoding sequences surrounding the genes. Lever analysis allows unbiased inference of functional annotations to regulatory motifs and candidate CRMs. We used human myogenic differentiation as a model system to statistically assess greater than 25,000 pairings of gene sets and motifs or motif combinations. We assigned functional annotations to candidate regulatory motifs predicted previously and identified gene sets that are likely to be co-regulated via shared regulatory motifs. Lever allows moving beyond the identification of putative regulatory motifs in mammalian genomes, toward understanding their biological roles. This approach is general and can be applied readily to any cell type, gene expression pattern or organism of interest.

High-throughput quantitative polymerase chain reaction in picoliter droplets.

Limiting dilution PCR has become an increasingly useful technique for the detection and quantification of rare species in a population, but the limit of detection and accuracy of quantification are largely determined by the number of reactions that can be analyzed. Increased throughput may be achieved by reducing the reaction volume and increasing processivity. We have designed a high-throughput microfluidic chip that encapsulates PCR reagents in millions of picoliter droplets in a continuous oil flow. The oil stream conducts the droplets through alternating denaturation and annealing zones, resulting in rapid (55-s cycles) and efficient PCR amplification. Inclusion of fluorescent probes in the PCR reaction mix permits the amplification process to be monitored within individual droplets at specific locations within the microfluidic chip. We show that amplification of a 245-bp adenovirus product can be detected and quantified in 35 min at starting template concentrations as low as 1 template molecule/167 droplets (0.003 pg/microL). The frequencies of positive reactions over a range of template concentrations agree closely with the frequencies predicted by Poisson statistics, demonstrating both the accuracy and sensitivity of this platform for limiting dilution and digital PCR applications.

Role of N-linked oligosaccharide flexibility in mannose phosphorylation of lysosomal enzyme cathepsin L.

Mannose phosphorylation of N-linked oligosaccharides by UDP-GlcNAc:lysosomal enzyme N-acetylglucosamine-1-phosphotransferase is a key step in the targeting of lysosomal enzymes in mammalian cells and tissues. The selectivity of this process is determined by lysine-based phosphorylation signals shared by lysosomal enzymes of diverse structure and function. By introducing new glycosylation sites at several locations on the surface of mouse procathepsin L and modeling oligosaccharide conformations for sites that are phosphorylated, it was shown that the inherent flexibility of N-linked oligosaccharides can account for the specificity of the transferase for oligosaccharides at different locations on the protein. By using this approach, the physical relationship between the lysine-based signal and the site of phosphorylation of mannose residues was determined. The analysis also revealed the existence of additional independent lysine-based phosphorylation signals on procathepsin L, which account for the low level of phosphorylation observed when the primary Lys-54/Lys-99 signal is ablated. Mutagenesis of residues that surround Lys-54 and Lys-99 and demonstration of mannose phosphorylation of a glycosylated derivative of green fluorescent protein provide strong evidence that the cathepsin L phosphorylation signal is a simple structure composed of as few as two well placed lysine residues.