Differential transcriptional responses to Ebola and Marburg virus infection in bat and human cells

The unprecedented outbreak of Ebola in West Africa resulted in over 28,000 cases and 11,000 deaths, underlining the need for a better understanding of the biology of this highly pathogenic virus to develop specific counter strategies. Two filoviruses, the Ebola and Marburg viruses, result in a severe and often fatal infection in humans. However, bats are natural hosts and survive filovirus infections without obvious symptoms. The molecular basis of this striking difference in the response to filovirus infections is not well understood. We report a systematic overview of differentially expressed genes, activity motifs and pathways in human and bat cells infected with the Ebola and Marburg viruses, and we demonstrate that the replication of filoviruses is more rapid in human cells than in bat cells. We also found that the most strongly regulated genes upon filovirus infection are chemokine ligands and transcription factors. We observed a strong induction of the JAK/STAT pathway, of several genes encoding inhibitors of MAP kinases (DUSP genes) and of PPP1R15A, which is involved in ER stress-induced cell death. We used comparative transcriptomics to provide a data resource that can be used to identify cellular responses that might allow bats to survive filovirus infections.Additional co-authors: Andreas J. Gruber, Franziska Hufsky, Henrike Indrischek, Sabina Kanton, Jörg Linde, Nelly Mostajo, Roman Ochsenreiter, Konstantin Riege, Lorena Rivarola-Duarte, Abdullah H. Sahyoun, Sita J. Saunders, Stefan E. Seemann, Andrea Tanzer, Bertram Vogel, Michael T. Wolfinger, Rolf Backofen, Jan Gorodkin, Ivo Grosse, Ivo Hofacker, Steve Hoffmann, Christoph Kaleta, Peter F. Stadler, Stephan Becker, and Manja Marz

Hydrogen-Bond Driven Loop-Closure Kinetics in Unfolded Polypeptide Chains

Characterization of the length dependence of end-to-end loop-closure kinetics in unfolded polypeptide chains provides an understanding of early steps in protein folding. Here, loop-closure in poly-glycine-serine peptides is investigated by combining single-molecule fluorescence spectroscopy with molecular dynamics simulation. For chains containing more than 10 peptide bonds loop-closing rate constants on the 20–100 nanosecond time range exhibit a power-law length dependence. However, this scaling breaks down for shorter peptides, which exhibit slower kinetics arising from a perturbation induced by the dye reporter system used in the experimental setup. The loop-closure kinetics in the longer peptides is found to be determined by the formation of intra-peptide hydrogen bonds and transient β-sheet structure, that accelerate the search for contacts among residues distant in sequence relative to the case of a polypeptide chain in which hydrogen bonds cannot form. Hydrogen-bond-driven polypeptide-chain collapse in unfolded peptides under physiological conditions found here is not only consistent with hierarchical models of protein folding, that highlights the importance of secondary structure formation early in the folding process, but is also shown to speed up the search for productive folding events

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis

Revealing a brain network endophenotype in families with idiopathic generalised epilepsy

Idiopathic generalised epilepsy (IGE) has a genetic basis. The mechanism of seizure expression is not fully known, but is assumed to involve large-scale brain networks. We hypothesised that abnormal brain network properties would be detected using EEG in patients with IGE, and would be manifest as a familial endophenotype in their unaffected first-degree relatives. We studied 117 participants: 35 patients with IGE, 42 unaffected first-degree relatives, and 40 normal controls, using scalp EEG. Graph theory was used to describe brain network topology in five frequency bands for each subject. Frequency bands were chosen based on a published Spectral Factor Analysis study which demonstrated these bands to be optimally robust and independent. Groups were compared, using Bonferroni correction to account for nonindependent measures and multiple groups. Degree distribution variance was greater in patients and relatives than controls in the 6-9 Hz band (p = 0.0005, p = 0.0009 respectively). Mean degree was greater in patients than healthy controls in the 6-9 Hz band (p = 0.0064). Clustering coefficient was higher in patients and relatives than controls in the 6-9 Hz band (p = 0.0025, p = 0.0013). Characteristic path length did not differ between groups. No differences were found between patients and unaffected relatives. These findings suggest brain network topology differs between patients with IGE and normal controls, and that some of these network measures show similar deviations in patients and in unaffected relatives who do not have epilepsy. This suggests brain network topology may be an inherited endophenotype of IGE, present in unaffected relatives who do not have epilepsy, as well as in affected patients. We propose that abnormal brain network topology may be an endophenotype of IGE, though not in itself sufficient to cause epilepsy

Bayesian inference of accurate population sizes and FRET efficiencies from single diffusing biomolecules.

It is of significant biophysical interest to obtain accurate intramolecular distance information and population sizes from single-molecule Förster resonance energy transfer (smFRET) data obtained from biomolecules in solution. Experimental methods of increasing cost and complexity are being developed to improve the accuracy and precision of data collection. However, the analysis of smFRET data sets currently relies on simplistic, and often arbitrary methods, for the selection and denoising of fluorescent bursts. Although these methods are satisfactory for the analysis of simple, low-noise systems with intermediate FRET efficiencies, they display systematic inaccuracies when applied to more complex systems. We have developed an inference method for the analysis of smFRET data from solution studies based on rigorous model-based Bayesian techniques. We implement a Monte Carlo Markov chain (MCMC) based algorithm that simultaneously estimates population sizes and intramolecular distance information directly from a raw smFRET data set, with no intermediate event selection and denoising steps. Here, we present both our parametric model of the smFRET process and the algorithm developed for data analysis. We test the algorithm using a combination of simulated data sets and data from dual-labeled DNA molecules. We demonstrate that our model-based method systematically outperforms threshold-based techniques in accurately inferring both population sizes and intramolecular distances.This is the final published version. It's also available from ACS in Analytical Chemistry: http://pubs.acs.org/doi/pdf/10.1021/ac501188r

Cooperation of local motions in the Hsp90 molecular chaperone ATPase mechanism

The Hsp90 chaperone is a central node of protein homeostasis activating a large number of diverse client proteins. Hsp90 functions as a molecular clamp that closes and opens in response to the binding and hydrolysis of ATP. Crystallographic studies define distinct conformational states of the mechanistic core implying structural changes that have not yet been observed in solution. Here, we engineered one-nanometer fluorescence probes based on photo-induced electron transfer into yeast Hsp90 to observe these motions. We found that the ATPase activity of the chaperone was reflected in the kinetics of specific structural rearrangements at remote positions that acted cooperatively. Nanosecond single-molecule fluorescence fluctuation analysis uncovered that critical structural elements that undergo rearrangement are mobile on a sub-millisecond time scale. We identified a two-step mechanism for lid closure over the nucleotide-binding pocket. The activating co-chaperone Aha1 mobilizes the lid of apo Hsp90, suggesting an early role in the catalytic cycle

Genome-Wide Copy Number Variation in Epilepsy: Novel Susceptibility Loci in Idiopathic Generalized and Focal Epilepsies

Epilepsy is one of the most common neurological disorders in humans with a prevalence of 1% and a lifetime incidence of 3%. Several genes have been identified in rare autosomal dominant and severe sporadic forms of epilepsy, but the genetic cause is unknown in the vast majority of cases. Copy number variants (CNVs) are known to play an important role in the genetic etiology of many neurodevelopmental disorders, including intellectual disability (ID), autism, and schizophrenia. Genome-wide studies of copy number variation in epilepsy have not been performed. We have applied whole-genome oligonucleotide array comparative genomic hybridization to a cohort of 517 individuals with various idiopathic, non-lesional epilepsies. We detected one or more rare genic CNVs in 8.9% of affected individuals that are not present in 2,493 controls; five individuals had two rare CNVs. We identified CNVs in genes previously implicated in other neurodevelopmental disorders, including two deletions in AUTS2 and one deletion in CNTNAP2. Therefore, our findings indicate that rare CNVs are likely to contribute to a broad range of generalized and focal epilepsies. In addition, we find that 2.9% of patients carry deletions at 15q11.2, 15q13.3, or 16p13.11, genomic hotspots previously associated with ID, autism, or schizophrenia. In summary, our findings suggest common etiological factors for seemingly diverse diseases such as ID, autism, schizophrenia, and epilepsy

Structural Heterogeneity and Quantitative FRET Efficiency Distributions of Polyprolines through a Hybrid Atomistic Simulation and Monte Carlo Approach

Förster Resonance Energy Transfer (FRET) experiments probe molecular distances via distance dependent energy transfer from an excited donor dye to an acceptor dye. Single molecule experiments not only probe average distances, but also distance distributions or even fluctuations, and thus provide a powerful tool to study biomolecular structure and dynamics. However, the measured energy transfer efficiency depends not only on the distance between the dyes, but also on their mutual orientation, which is typically inaccessible to experiments. Thus, assumptions on the orientation distributions and averages are usually made, limiting the accuracy of the distance distributions extracted from FRET experiments. Here, we demonstrate that by combining single molecule FRET experiments with the mutual dye orientation statistics obtained from Molecular Dynamics (MD) simulations, improved estimates of distances and distributions are obtained. From the simulated time-dependent mutual orientations, FRET efficiencies are calculated and the full statistics of individual photon absorption, energy transfer, and photon emission events is obtained from subsequent Monte Carlo (MC) simulations of the FRET kinetics. All recorded emission events are collected to bursts from which efficiency distributions are calculated in close resemblance to the actual FRET experiment, taking shot noise fully into account. Using polyproline chains with attached Alexa 488 and Alexa 594 dyes as a test system, we demonstrate the feasibility of this approach by direct comparison to experimental data. We identified cis-isomers and different static local environments as sources of the experimentally observed heterogeneity. Reconstructions of distance distributions from experimental data at different levels of theory demonstrate how the respective underlying assumptions and approximations affect the obtained accuracy. Our results show that dye fluctuations obtained from MD simulations, combined with MC single photon kinetics, provide a versatile tool to improve the accuracy of distance distributions that can be extracted from measured single molecule FRET efficiencies

Genomic and transcriptomic changes complement each other in the pathogenesis of sporadic Burkitt lymphoma

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing