ALBERT

Description: A novel mechanism of collisionless heating in large planar arrays of small inductive coils operated at radio frequencies is presented. In contrast to the well-known case of non-local heating related to the transversal conductivity, when the electrons move perpendicular to the planar coil, we investigate the problem of electrons moving in a plane parallel to the coils. Two types of periodic structures are studied. Resonance velocities where heating is efficient are calculated analytically by solving the Vlasov equation. Certain scaling parameters are identified. The concept is further investigated by a single particle simulation based on the ergodic principle and combined with a Monte Carlo code allowing for collisions with Argon atoms. Resonances, energy exchange, and distribution functions are obtained. The analytical results are confirmed by the numerical simulation. Pressure and electric field dependences are studied. Stochastic heating is found to be most efficient when the electron mean free path exceeds the size of a single coil cell. Then the mean energy increases approximately exponentially with the electric field amplitude.

Description: A pulsed slow-positron beam generated by an electron linear accelerator was directly used for positron annihilation lifetime spectroscopy without any positron storage devices. A waveform digitizer was introduced to simultaneously capture multiple gamma-ray signals originating from positron annihilation events during a single accelerator pulse. The positron pulse was chopped and bunched with the chopper signals also sent to the waveform digitizer. Time differences between the annihilation gamma-ray and chopper peaks were calculated and accumulated as lifetime spectra in a computer. The developed technique indicated that positron annihilation lifetime spectroscopy can be performed in a 20 μ s time window at a pulse repetition rate synchronous with the linear accelerator. Lifetime spectra of a Kapton sheet and a thermally grown SiO 2 layer on Si were successfully measured. Synchronization of positron lifetime measurements with pulsed ion irradiation was demonstrated by this technique.

Description: Surface nuclear magnetic resonance (surface-NMR) is a promising technique for exploring shallow subsurface aquifer structures. Surface-NMR can be applied in environments that are characterized as a 1-D layered Earth. The technique utilizes a single loop and is referred to as magnetic resonance sounding. The technique referred to as magnetic resonance tomography (MRT) allows complex 2-D aquifer structures to be explored. Currently, MRT requires multiple loops and a roll along measurement scheme, which causes long survey time. We propose a loop layout using an elongated transmitter and an in-loop receiver arrays (ETRA) to conduct a 2-D survey with just one measurement. We present a comprehensive comparison between the new layout and the common approaches based on sensitivity and resolution analyses and show synthetic and field data. The results show that ETRA generates subsurface images at sufficient resolution with significantly lower survey times than other loop layouts.

Description: Measurements of ground deformation can be used to identify and interpret geophysical processes occurring at volcanoes. Most studies rely on a single geodetic technique, or fit a geophysical model to the results of multiple geodetic techniques. Here we present a methodology that combines GPS, Total Station measurements and InSAR into a single reference frame to produce an integrated 3-D geodetic velocity surface without any prior geophysical assumptions. The methodology consists of five steps: design of the network, acquisition and processing of the data, spatial integration of the measurements, time series computation and finally the integration of spatial and temporal measurements. The most significant improvements of this method are (1) the reduction of the required field time, (2) the unambiguous detection of outliers, (3) an increased measurement accuracy and (4) the construction of a 3-D geodetic velocity field. We apply this methodology to ongoing motion on Arenal's western flank. Integration of multiple measurement techniques at Arenal volcano revealed a deformation field that is more complex than that described by individual geodetic techniques, yet remains consistent with previous studies. This approach can be applied to volcano monitoring worldwide and has the potential to be extended to incorporate other geodetic techniques and to study transient deformation.

Description: Structured metal surfaces could support spoof surface plasmon polaritons (SPPs), the dispersion of which is determined by the cutoff condition of guided modes in the nanostructures. We show that we can achieve split spoof SPPs by breaking the degeneracy of guided helical modes in concentric nanostructures via the classic analogue of the Zeeman effect. This split effect is shown to be observable from the spectra of enhanced electromagnetic transmission. Spin-sensitive enhanced electromagnetic transmission and the associated characteristics of field are investigated. Transmission branches versus parallel wavevector can be satisfactorily fitted by using the dispersion of spoof SPPs.

Description: Light-induced degradation (LID) is a deleterious effect in crystalline silicon, which is considered to originate from recombination-active boron-oxygen complexes and/or copper-related defects. Although LID in both cases appears as a fast initial decay followed by a second slower degradation, we show that the time constant of copper-related degradation increases with increasing boron concentration in contrast to boron-oxygen LID. Temperature-dependent analysis reveals that the defect formation is limited by copper diffusion. Finally, interface defect density measurements confirm that copper-related LID is dominated by recombination in the wafer bulk.

Description: In this paper, we report the existence of intervalence charge transfer (IVCT) luminescence in Yb-doped fluorite-type crystals associated with Yb 2+ –Yb 3+ mixed valence pairs. By means of embedded cluster, wave function theory ab initio calculations, we show that the widely studied, very broad band, anomalous emission of Yb 2+ -doped CaF 2 and SrF 2 , usually associated with impurity-trapped excitons, is, rather, an IVCT luminescence associated with Yb 2+ –Yb 3+ mixed valence pairs. The IVCT luminescence is very efficiently excited by a two-photon upconversion mechanism where each photon provokes the same strong 4 f 14 –1 A 1g → 4 f 13 ( 2 F 7/2 )5 de g –1 T 1u absorption in the Yb 2+ part of the pair: the first one, from the pair ground state; the second one, from an excited state of the pair whose Yb 3+ moiety is in the higher 4 f 13 ( 2 F 5/2 ) multiplet. The Yb 2+ –Yb 3+ → Yb 3+ –Yb 2+ IVCT emission consists of an Yb 2+ 5 de g → Yb 3+ 4 f 7/2 charge transfer accompanied by a 4 f 7/2 → 4 f 5/2 deexcitation within the Yb 2+ 4 f 13 subshell: [ 2 F 5/2 5 de g , 2 F 7/2 ] → [ 2 F 7/2 ,4 f 14 ]. The IVCT vertical transition leaves the oxidized and reduced moieties of the pair after electron transfer very far from their equilibrium structures; this explains the unexpectedly large band width of the emission band and its low peak energy, because the large reorganization energies are subtracted from the normal emission. The IVCT energy diagrams resulting from the quantum mechanical calculations explain the different luminescent properties of Yb-doped CaF 2 , SrF 2 , BaF 2 , and SrCl 2 : the presence of IVCT luminescence in Yb-doped CaF 2 and SrF 2 ; its coexistence with regular 5 d -4 f emission in SrF 2 ; its absence in BaF 2 and SrCl 2 ; the quenching of all emissions in BaF 2 ; and the presence of additional 5 d –4 f emissions in SrCl 2 which are absent in SrF 2 . They also allow to interpret and reproduce recent experiments on transient photoluminescence enhancement in Yb 2+ -doped CaF 2 and SrF 2 , the appearance of Yb 2+ 4 f –5 d absorption bands in the excitation spectra of the IR Yb 3+ emission in partly reduced CaF 2 :Yb 3+ samples, and to identify the broadband observed in the excitation spectrum of the so far called anomalous emission of SrF 2 :Yb 2+ as an IVCT absorption, which corresponds to an Yb 2+ 4 f 5/2 → Yb 3+ 4 f 7/2 electron transfer.

Description: A classical limit of quantum dynamics can be defined by compensation of the quantum potential in the time-dependent Schrödinger equation. The quantum potential is a non-local quantity, defined in the trajectory-based form of the Schrödinger equation, due to Madelung, de Broglie, and Bohm, which formally generates the quantum-mechanical features in dynamics. Selective inclusion of the quantum potential for the degrees of freedom deemed “quantum,” defines a hybrid quantum/classical dynamics, appropriate for molecular systems comprised of light and heavy nuclei. The wavefunction is associated with all of the nuclei, and the Ehrenfest, or mean-field, averaging of the force acting on the classical degrees of freedom, typical of the mixed quantum/classical methods, is avoided. The hybrid approach is used to examine evolution of light/heavy systems in the harmonic and double-well potentials, using conventional grid-based and approximate quantum-trajectory time propagation. The approximate quantum force is defined on spatial domains, which removes unphysical coupling of the wavefunction fragments corresponding to distinct classical channels or configurations. The quantum potential, associated with the quantum particle, generates forces acting on both quantum and classical particles to describe the backreaction.

Description: In the present work, we redefine and generalize the action principle for dissipative systems proposed by Riewe by fixing the mathematical inconsistencies present in the original approach. In order to formulate a quadratic Lagrangian for non-conservative systems, the Lagrangian functions proposed depend on mixed integer order and fractional order derivatives. As examples, we formulate a quadratic Lagrangian for a particle under a frictional force proportional to the velocity and to the classical problem of an accelerated point charge.

Description: We consider a very natural generalization of quantum theory by letting the dimension of the Bloch ball be not necessarily three. We analyze bipartite state spaces where each of the components has a d -dimensional Euclidean ball as state space. In addition to this, we impose two very natural assumptions: the continuity and reversibility of dynamics and the possibility of characterizing bipartite states by local measurements. We classify all these bipartite state spaces and prove that, except for the quantum two-qubit state space, none of them contains entangled states. Equivalently, in any of these non-quantum theories, interacting dynamics is impossible. This result reveals that “existence of entanglement” is the requirement with minimal logical content which singles out quantum theory from our family of theories.

Description: A covariant algorithm for deriving the conserved quantities for natural Hamiltonian systems is combined with the non-relativistic framework of Eisenhart, and of Duval, in which the classical trajectories arise as geodesics in a higher dimensional space-time, realized by Brinkmann manifolds. Conserved quantities which are polynomial in the momenta can be built using time-dependent conformal Killing tensors with flux. The latter are associated with terms proportional to the Hamiltonian in the lower dimensional theory and with spectrum generating algebras for higher dimensional quantities of order 1 and 2 in the momenta. Illustrations of the general theory include the Runge-Lenz vector for planetary motion with a time-dependent gravitational constant G ( t ), motion in a time-dependent electromagnetic field of a certain form, quantum dots, the Hénon-Heiles and Holt systems, respectively, providing us with Killing tensors of rank that ranges from one to six.

Description: We demonstrate both theoretically and experimentally that the microwave-shielding effectiveness of a double-layer metallic mesh with a submillimeter period can be improved by increasing the separation between the two mesh layers (without affecting transmittance). This double-layer mesh consists of two layers of square aluminum mesh separated by a quartz-glass substrate. By increasing the substrate's optical thickness from zero to λ / 4   of the shielding band's upper frequency, the shielding of the double-layer mesh improves considerably, owing to the increased reflectivity of the double-layer mesh with increasing separation in the low-frequency band. A Ku-band shielding effectiveness of over 32 dB is observed for the double-layer mesh with a normalized visible transmittance greater than 91%. It is found that the electromagnetic shielding effectiveness is enhanced by over 7 dB (80.0% energy attenuation) across the Ku-band, compared with that of a single-layer mesh, while the optical transmittances are almost identical for both tested structures. Such an enhancement permits the design of high-transparency optical elements with stronger microwave shielding that can be achieved using single-layer metallic mesh.

Description: This letter performed polarized microscopic laser Raman scattering spectroscopy on the curved edges of transferred epitaxial graphene on SiO 2 /Si. The intensity ratio between the parallel and perpendicular polarized D band is evolved, providing a spectroscopy-based technique to probe the atomic-scale edge structures in graphene. A detailed analysis procedure for non-ideal disordered curved edges of graphene is developed combining the atomic-scale zigzag and armchair edge structures along with some point defects. These results could provide valuable information of the realistic edges of graphene at the atomic-scale that can strongly influence the performance of graphene-based nanodevices.

Description: We measure and quantify properties of galactic outflows and diffuse gas at z ≥ 1 in cosmological hydrodynamical simulations. Our novel subresolution model, Multi-Phase Particle Integrator (MUPPI), implements supernova feedback using fully local gas properties, where the wind velocity and mass loading are not given as input. We find the following trends at z  = 2 by analysing central galaxies having a stellar mass higher than 10 9 M . The outflow velocity and mass outflow rate ( $\dot{M}_{\rm out}$ ) exhibit positive correlations with galaxy mass and with the star formation rate (SFR). However, most of the relations present a large scatter. The outflow mass loading factor () is between 0.2 and 10. The comparison effective model generates a constant outflow velocity, and a negative correlation of with halo mass. The number fraction of galaxies where outflow is detected decreases at lower redshifts, but remains more than 80 per cent over z  = 1–5. The outflow velocity correlation with SFR becomes flatter at z  = 1, and displays a negative correlation with halo mass in massive galaxies. Our study demonstrates that both the MUPPI and effective models produce significant outflows at ~1/10 of the virial radius; at the same time shows that the properties of outflows generated can be different from the input speed and mass loading in the effective model. Our MUPPI model, using local properties of gas in the subresolution recipe, is able to develop galactic outflows whose properties correlate with global galaxy properties , and consistent with observations.

Description: In this paper, we explore the effects of electrostatic parametric amplification on a high quality factor (Q 〉 100 000) encapsulated disk resonator gyroscope (DRG), fabricated in 〈100〉 silicon. The DRG was operated in the n = 2 degenerate wineglass mode at 235 kHz, and electrostatically tuned so that the frequency split between the two degenerate modes was less than 100 mHz. A parametric pump at twice the resonant frequency is applied to the sense axis of the DRG, resulting in a maximum scale factor of 156.6  μ V/(°/s), an 8.8× improvement over the non-amplified performance. When operated with a parametric gain of 5.4, a minimum angle random walk of 0.034°/√h and bias instability of 1.15°/h are achieved, representing an improvement by a factor of 4.3× and 1.5×, respectively.

Description: We propose a control element for a Josephson spin valve. It is a complex Josephson device containing ferromagnetic (F) layer in the weak-link area consisting of two regions, representing 0 and π Josephson junctions, respectively. The valve's state is defined by mutual orientations of the F-layer magnetization vector and boundary line between 0 and π sections of the device. We consider possible implementation of the control element by introduction of a thin normal metal layer in a part of the device area. By means of theoretical simulations, we study properties of the valve's structure as well as its operation, revealing such advantages as simplicity of control, high characteristic frequency, and good legibility of the basic states.

Description: Thin films of the lead-free ferroelectric Na 0.5 Bi 0.5 TiO 3 grown on thin-film Pt electrodes supported by SrTiO 3 substrates have a complex microstructure consisting of crystalline grains with three distinct major crystallographic orientations. The piezoelectric response measured in spatially separated sub-micron grains using time-resolved synchrotron x-ray microdiffraction is highly inhomogeneous even among grains sharing the same major orientation. The piezoelectric coefficient d 33 varies by nearly a factor of two in a series of areas sharing the 〈001〉 orientation. The piezoelectric inhomogeneity is linked to the peculiar microstructure of the film, arising from local variations in the stress imposed by surrounding grains with different crystallographic orientations and differing directions of the ferroelectric remnant polarization. A systematic nonlinearity of the piezoelectric strain is observed in applied electric fields with small magnitudes in all regions, consistent with the coexistence of domains of differing polarization direction at zero applied electric field.

Description: The advent in high-throughput-sequencing (HTS) technologies has revolutionized conventional biodiversity research by enabling parallel capture of DNA sequences possessing species-level diagnosis. However, polymerase chain reaction (PCR)-based implementation is biased by the efficiency of primer binding across lineages of organisms. A PCR-free HTS approach will alleviate this artefact and significantly improve upon the multi-locus method utilizing full mitogenomes. Here we developed a novel multiplex sequencing and assembly pipeline allowing for simultaneous acquisition of full mitogenomes from pooled animals without DNA enrichment or amplification. By concatenating assemblies from three de novo assemblers, we obtained high-quality mitogenomes for all 49 pooled taxa, with 36 species 〉15 kb and the remaining 〉10 kb, including 20 complete mitogenomes and nearly all protein coding genes (99.6%). The assembly quality was carefully validated with Sanger sequences, reference genomes and conservativeness of protein coding genes across taxa. The new method was effective even for closely related taxa, e.g. three Drosophila spp., demonstrating its broad utility for biodiversity research and mito-phylogenomics. Finally, the in silico simulation showed that by recruiting multiple mito-loci, taxon detection was improved at a fixed sequencing depth. Combined, these results demonstrate the plausibility of a multi-locus mito-metagenomics approach as the next phase of the current single-locus metabarcoding method.

Description: The efficacy and the mutation spectrum of genome editing methods can vary substantially depending on the targeted sequence. A simple, quick assay to accurately characterize and quantify the induced mutations is therefore needed. Here we present TIDE, a method for this purpose that requires only a pair of PCR reactions and two standard capillary sequencing runs. The sequence traces are then analyzed by a specially developed decomposition algorithm that identifies the major induced mutations in the projected editing site and accurately determines their frequency in a cell population. This method is cost-effective and quick, and it provides much more detailed information than current enzyme-based assays. An interactive web tool for automated decomposition of the sequence traces is available. TIDE greatly facilitates the testing and rational design of genome editing strategies.

Description: NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients ( r 2 ) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.

Description: Bicyclic oxazaphospholidine monomers were used to prepare a series of phosphorothioate (PS)-modified gapmer antisense oligonucleotides (ASOs) with control of the chirality of each of the PS linkages within the 10-base gap. The stereoselectivity was determined to be 98% for each coupling. The objective of this work was to study how PS chirality influences biophysical and biological properties of the ASO including binding affinity ( T m ), nuclease stability, activity in vitro and in vivo , RNase H activation and cleavage patterns (both human and E. coli ) in a gapmer context. Compounds that had nine or more Sp-linkages in the gap were found to be poorly active in vitro , while compounds with uniform Rp-gaps exhibited activity very similar to that of the stereo-random parent ASOs. Conversely, when tested in vivo , the full Rp-gap compound was found to be quickly metabolized resulting in low activity. A total of 31 ASOs were prepared with control of the PS chirally of each linkage within the gap in an attempt to identify favorable Rp/Sp positions. We conclude that a mix of Rp and Sp is required to achieve a balance between good activity and nuclease stability.

Description: In an isogenic cell population, phenotypic heterogeneity among individual cells is common and critical for survival of the population under different environment conditions. DNA modification is an important epigenetic factor that can regulate phenotypic heterogeneity. The single molecule real-time (SMRT) sequencing technology provides a unique platform for detecting a wide range of DNA modifications, including N6-methyladenine (6-mA), N4-methylcytosine (4-mC) and 5-methylcytosine (5-mC). Here we present qDNAmod, a novel bioinformatic tool for genome-wide quantitative profiling of intercellular heterogeneity of DNA modification from SMRT sequencing data. It is capable of estimating proportion of isogenic haploid cells, in which the same loci of the genome are differentially modified. We tested the reliability of qDNAmod with the SMRT sequencing data of Streptococcus pneumoniae strain ST556. qDNAmod detected extensive intercellular heterogeneity of DNA methylation (6-mA) in a clonal population of ST556. Subsequent biochemical analyses revealed that the recognition sequences of two type I restriction–modification (R-M) systems are responsible for the intercellular heterogeneity of DNA methylation initially identified by qDNAmod. qDNAmod thus represents a valuable tool for studying intercellular phenotypic heterogeneity from genome-wide DNA modification.

Description: Combinatorial transcription factor (TF) binding is essential for cell-type-specific gene regulation. However, much remains to be learned about the mechanisms of TF interactions, including to what extent constrained spacing and orientation of interacting TFs are critical for regulatory element activity. To examine the relative prevalence of the ‘enhanceosome’ versus the ‘TF collective’ model of combinatorial TF binding, a comprehensive analysis of TF binding site sequences in large scale datasets is necessary. We developed a motif-pair discovery pipeline to identify motif co-occurrences with preferential distance(s) between motifs in TF-bound regions. Utilizing a compendium of 289 mouse haematopoietic TF ChIP-seq datasets, we demonstrate that haematopoietic-related motif-pairs commonly occur with highly conserved constrained spacing and orientation between motifs. Furthermore, motif clustering revealed specific associations for both heterotypic and homotypic motif-pairs with particular haematopoietic cell types. We also showed that disrupting the spacing between motif-pairs significantly affects transcriptional activity in a well-known motif-pair—E-box and GATA, and in two previously unknown motif-pairs with constrained spacing—Ets and Homeobox as well as Ets and E-box. In this study, we provide evidence for widespread sequence-specific TF pair interaction with DNA that conforms to the ‘enhanceosome’ model, and furthermore identify associations between specific haematopoietic cell-types and motif-pairs.

Description: In silico tools have been developed to predict variants that may have an impact on pre-mRNA splicing. The major limitation of the application of these tools to basic research and clinical practice is the difficulty in interpreting the output. Most tools only predict potential splice sites given a DNA sequence without measuring splicing signal changes caused by a variant. Another limitation is the lack of large-scale evaluation studies of these tools. We compared eight in silico tools on 2959 single nucleotide variants within splicing consensus regions (scSNVs) using receiver operating characteristic analysis. The Position Weight Matrix model and MaxEntScan outperformed other methods. Two ensemble learning methods, adaptive boosting and random forests, were used to construct models that take advantage of individual methods. Both models further improved prediction, with outputs of directly interpretable prediction scores. We applied our ensemble scores to scSNVs from the Catalogue of Somatic Mutations in Cancer database. Analysis showed that predicted splice-altering scSNVs are enriched in recurrent scSNVs and known cancer genes. We pre-computed our ensemble scores for all potential scSNVs across the human genome, providing a whole genome level resource for identifying splice-altering scSNVs discovered from large-scale sequencing studies.

Description: In mammals, RNA interference is primarily a post-transcriptional mechanism. Evidence has accumulated for additional role in transcriptional gene silencing (TGS) but the question for a good paradigm for small interfering antigene RNA (agRNA)-induced chromatin modification remains unanswered. Here, we show that SETDB1, a histone H3-lysine 9 (H3K9)-specific methyltransferase, cooperates with Argonaute-2 (AGO2) and plays an essential role in agRNA-induced TGS. The androgen receptor ( AR ) gene was transcriptionally silenced by agRNA targeted to its promoter, and we show that this repression was mitigated by knockdown of SETDB1 or AGO2 . Chromatin immunoprecipitation demonstrated that agRNA-driven AGO2 was first targeted to the AR promoter, followed by SETDB1. SIN3A and HDAC1/2, the components of the SIN3-HDAC complex, immunoprecipitated with SETDB1, and localized at the agRNA-targeted promoter. Agreeing with the presence of SETDB1, trimethyl-H3K9 was enriched in the AR promoter. Both EZH2 and trimethyl-H3K27 were also present in the targeted locus; accordingly, EZH2 immunoprecipitated with SETDB1. DNA methylation level was not significantly changed, suggesting the absence of de novo methylating activity in agRNA-induced AR promoter. Our results demonstrate that SETDB1, together with AGO2, plays an essential role in TGS through recruiting chromatin remodeler and/or other modifiers, consequently creating a repressive chromatin milieu at the targeted promoter.

Description: LEF/TCFs direct the final step in Wnt/β-catenin signalling by recruiting β-catenin to genes for activation of transcription. Ancient, non-vertebrate TCFs contain two DNA binding domains, a High Mobility Group box for recognition of the Wnt Response Element (WRE; 5'-CTTTGWWS-3') and the C-clamp domain for recognition of the GC-rich Helper motif (5'-RCCGCC-3'). Two vertebrate TCFs (TCF-1/TCF7 and TCF-4/TCF7L2) use the C-clamp as an alternatively spliced domain to regulate cell-cycle progression, but how the C-clamp influences TCF binding and activity genome-wide is not known. Here, we used a doxycycline inducible system with ChIP-seq to assess how the C-clamp influences human TCF1 binding genome-wide. Metabolic pulse-labeling of nascent RNA with 4'Thiouridine was used with RNA-seq to connect binding to the Wnt transcriptome. We find that the C-clamp enables targeting to a greater number of gene loci for stronger occupancy and transcription regulation. The C-clamp uses Helper sites concurrently with WREs for gene targeting, but it also targets TCF1 to sites that do not have readily identifiable canonical WREs. The coupled ChIP-seq/4'Thiouridine-seq analysis identified new Wnt target genes, including additional regulators of cell proliferation. Thus, C-clamp containing isoforms of TCFs are potent transcriptional regulators with an expanded transcriptome directed by C-clamp-Helper site interactions.

Description: Chromatin modifiers and histone modifications are components of a chromatin-signaling network involved in transcription and its regulation. The interactions between chromatin modifiers and histone modifications are often unknown, are based on the analysis of few genes or are studied in vitro . Here, we apply computational methods to recover interactions between chromatin modifiers and histone modifications from genome-wide ChIP-Seq data. These interactions provide a high-confidence backbone of the chromatin-signaling network. Many recovered interactions have literature support; others provide hypotheses about yet unknown interactions. We experimentally verified two of these predicted interactions, leading to a link between H4K20me1 and members of the Polycomb Repressive Complexes 1 and 2. Our results suggest that our computationally derived interactions are likely to lead to novel biological insights required to establish the connectivity of the chromatin-signaling network involved in transcription and its regulation.

Description: As the biomedical impact of small RNAs grows, so does the need to understand competing structural alternatives for regions of functional interest. Suboptimal structure analysis provides significantly more RNA base pairing information than a single minimum free energy prediction. Yet computational enhancements like Boltzmann sampling have not been fully adopted by experimentalists since identifying meaningful patterns in this data can be challenging. Profiling is a novel approach to mining RNA suboptimal structure data which makes the power of ensemble-based analysis accessible in a stable and reliable way. Balancing abstraction and specificity, profiling identifies significant combinations of base pairs which dominate low-energy RNA secondary structures. By design, critical similarities and differences are highlighted, yielding crucial information for molecular biologists. The code is freely available via http://gtfold.sourceforge.net/profiling.html .

Description: The stress-sensitive restriction-modification (RM) system CglI from Corynebacterium glutamicum and the homologous NgoAVII RM system from Neisseria gonorrhoeae FA1090 are composed of three genes: a DNA methyltransferase (M.CglI and M.NgoAVII), a putative restriction endonuclease (R.CglI and R.NgoAVII, or R-proteins) and a predicted DEAD-family helicase/ATPase (N.CglI and N.NgoAVII or N-proteins). Here we report a biochemical characterization of the R- and N-proteins. Size-exclusion chromatography and SAXS experiments reveal that the isolated R.CglI, R.NgoAVII and N.CglI proteins form homodimers, while N.NgoAVII is a monomer in solution. Moreover, the R.CglI and N.CglI proteins assemble in a complex with R 2 N 2 stoichiometry. Next, we show that N-proteins have ATPase activity that is dependent on double-stranded DNA and is stimulated by the R-proteins. Functional ATPase activity and extensive ATP hydrolysis (~170 ATP/s/monomer) are required for site-specific DNA cleavage by R-proteins. We show that ATP-dependent DNA cleavage by R-proteins occurs at fixed positions (6–7 nucleotides) downstream of the asymmetric recognition sequence 5'-GCCGC-3'. Despite similarities to both Type I and II restriction endonucleases, the CglI and NgoAVII enzymes may employ a unique catalytic mechanism for DNA cleavage.

Description: Emerging evidence points to roles for tRNA modifications and tRNA abundance in cellular stress responses. While isolated instances of stress-induced tRNA degradation have been reported, we sought to assess the effects of stress on tRNA levels at a systems level. To this end, we developed a next-generation sequencing method that exploits the paucity of ribonucleoside modifications at the 3'-end of tRNAs to quantify changes in all cellular tRNA molecules. Application of this tRNA-seq method to Saccharomyces cerevisiae identified all 76 expressed unique tRNA species out of 295 coded in the yeast genome, including all isoacceptor variants, with highly precise relative (fold-change) quantification of tRNAs. In studies of stress-induced changes in tRNA levels, we found that oxidation (H 2 O 2 ) and alkylation (methylmethane sulfonate, MMS) stresses induced nearly identical patterns of up- and down-regulation for 58 tRNAs. However, 18 tRNAs showed opposing changes for the stresses, which parallels our observation of signature reprogramming of tRNA modifications caused by H 2 O 2 and MMS. Further, stress-induced degradation was limited to only a small proportion of a few tRNA species. With tRNA-seq applicable to any organism, these results suggest that translational control of stress response involves a contribution from tRNA abundance.

Description: During meiosis programmed DNA double-strand breaks (DSBs) are repaired by homologous recombination using the sister chromatid or the homologous chromosome (homolog) as a template. This repair results in crossover (CO) and non-crossover (NCO) recombinants. Only CO formation between homologs provides the physical linkages guiding correct chromosome segregation, which are essential to produce healthy gametes. The factors that determine the CO/NCO decision are still poorly understood. Using Schizosaccharomyces pombe as a model we show that the Rad51/Dmc1-paralog complexes Rad55-Rad57 and Rdl1-Rlp1-Sws1 together with Swi5-Sfr1 play a major role in antagonizing both the FANCM-family DNA helicase/translocase Fml1 and the RecQ-type DNA helicase Rqh1 to limit hybrid DNA formation and promote Mus81-Eme1-dependent COs. A common attribute of these protein complexes is an ability to stabilize the Rad51/Dmc1 nucleoprotein filament, and we propose that it is this property that imposes constraints on which enzymes gain access to the recombination intermediate, thereby controlling the manner in which it is processed and resolved.

Description: To study target sequence specificity, selectivity, and reaction kinetics of Streptococcus pyogenes Cas9 activity, we challenged libraries of random variant targets with purified Cas9::guide RNA complexes in vitro . Cleavage kinetics were nonlinear, with a burst of initial activity followed by slower sustained cleavage. Consistent with other recent analyses of Cas9 sequence specificity, we observe considerable (albeit incomplete) impairment of cleavage for targets mutated in the PAM sequence or in ‘seed’ sequences matching the proximal 8 bp of the guide. A second target region requiring close homology was located at the other end of the guide::target duplex (positions 13–18 relative to the PAM). Sequences flanking the guide+PAM region had measurable (albeit modest) effects on cleavage. In addition, the first-base Guanine constraint commonly imposed by gRNA expression systems has little effect on overall cleavage efficiency. Taken together, these studies provide an in vitro understanding of the complexities of Cas9–gRNA interaction and cleavage beyond the general paradigm of site determination based on the ‘seed’ sequence and PAM.

Description: Transcripts possessing a 5'-triphosphate are a hallmark of viral transcription and can trigger the host antiviral response. 5'-triphosphates are also found on common host transcripts transcribed by RNA polymerase III (RNAP III), yet how these transcripts remain non-immunostimulatory is incompletely understood. Most microRNAs (miRNAs) are 5'-monophosphorylated as a result of sequential endonucleolytic processing by Drosha and Dicer from longer RNA polymerase II (RNAP II)-transcribed primary transcripts. In contrast, bovine leukemia virus (BLV) expresses subgenomic RNAP III transcripts that give rise to miRNAs independent of Drosha processing. Here, we demonstrate that each BLV pre-miRNA is directly transcribed by RNAP III from individual, compact RNAP III type II genes. Thus, similar to manmade RNAP III-generated short hairpin RNAs (shRNAs), the BLV pre-miRNAs are initially 5'-triphosphorylated. Nonetheless, the derivative 5p miRNAs and shRNA-generated 5p small RNAs (sRNAs) possess a 5'-monophosphate. Our enzymatic characterization and small RNA sequencing data demonstrate that BLV 5p miRNAs are co-terminal with 5'-triphosphorylated miRNA precursors (pre-miRNAs). Thus, these results identify a 5'-tri-phosphatase activity that is involved in the biogenesis of BLV miRNAs and shRNA-generated sRNAs. This work advances our understanding of retroviral miRNA and shRNA biogenesis and may have implications regarding the immunostimulatory capacity of RNAP III transcripts.

Description: Within the field of synthetic biology, a rational design of genetic parts should include a causal understanding of their input-output responses—the so-called transfer function—and how to tune them. However, a commonly adopted strategy is to fit data to Hill-shaped curves without considering the underlying molecular mechanisms. Here we provide a novel mathematical formalization that allows prediction of the global behavior of a synthetic device by considering the actual information from the involved biological parts. This is achieved by adopting an enzymology-like framework, where transfer functions are described in terms of their input affinity constant and maximal response. As a proof of concept, we characterize a set of Lux homoserine-lactone-inducible genetic devices with different levels of Lux receptor and signal molecule. Our model fits the experimental results and predicts the impact of the receptor's ribosome-binding site strength, as a tunable parameter that affects gene expression. The evolutionary implications are outlined.

Description: The restriction endonuclease (REase) NgoAVII is composed of two proteins, R.NgoAVII and N.NgoAVII, and shares features of both Type II restriction enzymes and Type I/III ATP-dependent restriction enzymes (see accompanying paper Zaremba et al. , 2014). Here we present crystal structures of the R.NgoAVII apo-protein and the R.NgoAVII C-terminal domain bound to a specific DNA. R.NgoAVII is composed of two domains: an N-terminal nucleolytic PLD domain; and a C-terminal B3-like DNA-binding domain identified previously in BfiI and EcoRII REases, and in plant transcription factors. Structural comparison of the B3-like domains of R.NgoAVII, EcoRII, BfiI and the plant transcription factors revealed a conserved DNA-binding surface comprised of N- and C-arms that together grip the DNA. The C-arms of R.NgoAVII, EcoRII, BfiI and plant B3 domains are similar in size, but the R.NgoAVII N-arm which makes the majority of the contacts to the target site is much longer. The overall structures of R.NgoAVII and BfiI are similar; however, whilst BfiI has stand-alone catalytic activity, R.NgoAVII requires an auxiliary cognate N.NgoAVII protein and ATP hydrolysis in order to cleave DNA at the target site. The structures we present will help formulate future experiments to explore the molecular mechanisms of intersubunit crosstalk that control DNA cleavage by R.NgoAVII and related endonucleases.

Description: Distinct translational initiation mechanisms between prokaryotes and eukaryotes limit the exploitation of prokaryotic riboswitch repertoire for regulatory RNA circuit construction in mammalian application. Here, we explored programmed ribosomal frameshifting (PRF) as the regulatory gene expression platform for engineered ligand-responsive RNA devices in higher eukaryotes. Regulation was enabled by designed ligand-dependent conformational rearrangements of the two cis-acting RNA motifs of opposite activity in -1 PRF. Particularly, RNA elements responsive to trans-acting ligands can be tailored to modify co-translational RNA refolding dynamics of a hairpin upstream of frameshifting site to achieve reversible and adjustable -1 PRF attenuating activity. Combined with a ligand-responsive stimulator, synthetic RNA devices for synergetic translational-elongation control of gene expression can be constructed. Due to the similarity between co-transcriptional RNA hairpin folding and co-translational RNA hairpin refolding, the RNA-responsive ligand repertoire provided in prokaryotic systems thus becomes accessible to gene-regulatory circuit construction for synthetic biology application in mammalian cells.

Description: Cutaneous photosensitization is a common side effect of drug treatment and can be associated with an increased skin cancer risk. The immunosuppressant azathioprine, the fluoroquinolone antibiotics and vemurafenib—a BRAF inhibitor used to treat metastatic melanoma—are all recognized clinical photosensitizers. We have compared the effects of UVA radiation on cultured human cells treated with 6-thioguanine (6-TG, a DNA-embedded azathioprine surrogate), the fluoroquinolones ciprofloxacin and ofloxacin and vemurafenib. Despite widely different structures and modes of action, each of these drugs potentiated UVA cytotoxicity. UVA photoactivation of 6-TG, ciprofloxacin and ofloxacin was associated with the generation of singlet oxygen that caused extensive protein oxidation. In particular, these treatments were associated with damage to DNA repair proteins that reduced the efficiency of nucleotide excision repair. Although vemurafenib was also highly phototoxic to cultured cells, its effects were less dependent on singlet oxygen. Highly toxic combinations of vemurafenib and UVA caused little protein carbonylation but were nevertheless inhibitory to nucleotide excision repair. Thus, for three different classes of drugs, photosensitization by at least two distinct mechanisms is associated with reduced protection against potentially mutagenic and carcinogenic DNA damage.

Description: The first level of genome packaging in eukaryotic cells involves the formation of dense nucleosome arrays, with DNA coverage near 90% in yeasts. How cells achieve such high coverage within a short time, e.g. after DNA replication, remains poorly understood. It is known that random sequential adsorption of impenetrable particles on a line reaches high density extremely slowly, due to a jamming phenomenon. The nucleosome-shifting action of remodeling enzymes has been proposed as a mechanism to resolve such jams. Here, we suggest two biophysical mechanisms which assist rapid filling of DNA with nucleosomes, and we quantitatively characterize these mechanisms within mathematical models. First, we show that the ‘softness’ of nucleosomes, due to nucleosome breathing and stepwise nucleosome assembly, significantly alters the filling behavior, speeding up the process relative to ‘hard’ particles with fixed, mutually exclusive DNA footprints. Second, we explore model scenarios in which the progression of the replication fork could eliminate nucleosome jamming, either by rapid filling in its wake or via memory of the parental nucleosome positions. Taken together, our results suggest that biophysical effects promote rapid nucleosome filling, making the reassembly of densely packed nucleosomes after DNA replication a simpler task for cells than was previously thought.

Description: Abasic (AP) lesions are the most frequent type of damages occurring in cellular DNA. Here we describe the conformational effects of AP sites substituted for 2'-deoxyadenosine in the first ( ap7 ), second ( ap13 ) or third ( ap19 ) loop of the quadruplex formed in K + by the human telomere DNA 5'-d[AG 3 (TTAG 3 ) 3 ]. CD spectra and electrophoresis reveal that the presence of AP sites does not hinder the formation of intramolecular quadruplexes. NMR spectra show that the structural heterogeneity is substantially reduced in ap7 and ap19 as compared to that in the wild-type. These two ( ap7 and ap19 ) sequences are shown to adopt the hybrid-1 and hybrid-2 quadruplex topology, respectively, with AP site located in a propeller-like loop. All three studied sequences transform easily into parallel quadruplex in dehydrating ethanol solution. Thus, the AP site in any loop region facilitates the formation of the propeller loop. Substitution of all adenines by AP sites stabilizes the parallel quadruplex even in the absence of ethanol. Whereas guanines are the major determinants of quadruplex stability, the presence or absence of loop adenines substantially influences quadruplex folding. The naturally occurring adenine-lacking sites in the human telomere DNA can change the quadruplex topology in vivo with potentially vital biological consequences.

Description: The DNA sequence preferences of nearly all sequence specific DNA binding proteins are influenced by the identities of bases that are not directly contacted by protein. Discrimination between non-contacted base sequences is commonly based on the differential abilities of DNA sequences to allow narrowing of the DNA minor groove. However, the factors that govern the propensity of minor groove narrowing are not completely understood. Here we show that the differential abilities of various DNA sequences to support formation of a highly ordered and stable minor groove solvation network are a key determinant of non-contacted base recognition by a sequence-specific binding protein. In addition, disrupting the solvent network in the non-contacted region of the binding site alters the protein's ability to recognize contacted base sequences at positions 5–6 bases away. This observation suggests that DNA solvent interactions link contacted and non-contacted base recognition by the protein.

Description: We give a bordism-theoretic characterization of those closed almost contact $(2q{+ }1)$ -manifolds (with $q\geq 2$ ) that admit a Stein fillable contact structure. Our method is to apply Eliashberg's $h$ -principle for Stein manifolds in the setting of Kreck's modified surgery. As an application, we show that any simply connected almost contact 7-manifold with torsion-free second homotopy group is Stein fillable. We also discuss the Stein fillability of exotic spheres and examine subcritical Stein fillability.

Description: Casson-type invariants emerging from Donaldson theory over certain negative-definite four-manifolds were recently suggested by Teleman. These are defined by an algebraic count of points in a zero-dimensional moduli space of flat instantons. Motivated by the cobordism programme of proving Witten's conjecture, we use a moduli space of ${\rm PU}(2)$ Seiberg–Witten monopoles to exhibit an oriented one-dimensional cobordism of the instanton moduli space to the empty space. The Casson-type invariant must therefore vanish.

Description: The dynamical and stationary behaviors of a fourth-order equation in the unit ball with clamped boundary conditions and a singular reaction term are investigated. The equation arises in the modeling of microelectromechanical systems and includes a positive voltage parameter $\lambda$ . It is shown that there is a threshold value $\lambda _* 〉 0$ of the voltage parameter such that no radially symmetric stationary solution exists for $\lambda 〉 \lambda _* $ , while at least two such solutions exist for $\lambda \in (0,\lambda _* )$ . Local and global well-posedness results are obtained for the corresponding hyperbolic and parabolic evolution problems as well as the occurrence of finite time singularities when $\lambda 〉 \lambda _* $ .

Description: In a recent article, Serov et al. [J. Appl. Phys. 116 , 034507 (2014)] claim: “This study represents the first time that the high-field behavior in graphene on a substrate was investigated taking into account intrinsic graphene properties,” ignoring the most recent anisotropic distribution function [V. K. Arora et al ., J. Appl. Phys. 112 , 114330 (2012)] also published in J. Appl. Phys., targeting the same experimental data [V. E. Dorgan et al. , Appl. Phys. Lett. 97 , 082112 (2010)]. The claim of Serov et al. of being first is refuted and many shortcomings of the hydrodynamic model for a highly quantum and degenerate graphene nanolayer are pointed out.

Description: Fe thin films were deposited on sodium chloride (NaCl) substrates using magnetron sputtering to investigate means of texture control in free standing metal films. The Fe thin films were studied using transmission electron microscopy equipped with automated crystallographic orientation microscopy. Using this technique, the microstructure of each film was characterized in order to elucidate the effects of altering deposition parameters. The natural tendency for Fe films grown on (100) NaCl is to form a randomly oriented nanocrystalline microstructure. By careful selection of substrate and deposition conditions, it is possible to drive the texture of the film toward a single (100) orientation while retaining the nanocrystalline microstructure.

Description: The current paper reports successful syntheses of Ni–TiN composite coatings by pulse electrodeposition. The effect of pulse frequency on the microstructures, nanomechanical, and wear properties of the coatings was investigated using transmission electron microscopy, X–ray diffraction, nanoindenter, scanning electron microscopy, and wear test instrument. The results showed that the Ni–TiN composite coating prepared at the pulse frequency of 100 Hz showed the presence of a less number of TiN particles and some degrees of aggregation in micro-regions. By contrast, in the Ni–TiN coating deposited at the pulse frequency of 500 Hz, the TiN particles were large in number and dispersed homogeneously, thereby, offering the coating a uniform and fine structure. The average grain diameters of Ni and TiN in the coating prepared at 100 Hz were 154.7 and 44.8 nm, respectively, whereas those for the coating prepared at 500 Hz were 67.3 and 25.9 nm, respectively. The maximum TiN content in the Ni-TiN coating deposited at 800 Hz was approximately 10.5 wt. %. The maximum microhardness and the Young's modulus values for the Ni–TiN composite coatings deposited at 800 Hz were 35.7 GPa and 167.4 GPa, respectively. Furthermore, the Ni–TiN composite coating prepared at 100 Hz had more severe damages, whereas the morphologies of worn surface of the coatings deposited at 500 Hz and 800 Hz were smooth and only a few small pits appeared on the surface.

Description: Excited states and the ground state of the diatomic molecule RbSr were calculated by post Hartree-Fock molecular orbital theory up to 22 000 cm −1 . We applied a multireference configuration interaction calculation based on multiconfigurational self-consistent field wave functions. Both methods made use of effective core potentials and core polarization potentials. Potential energy curves, transition dipole moments, and permanent electric dipole moments were determined for RbSr and could be compared with other recent calculations. We found a good agreement with experimental spectra, which have been obtained recently by helium nanodroplet isolation spectroscopy. For the lowest two asymptotes (Rb (5 s 2 S) + Sr (5 s 4 d 3 P°) and Rb (5 p 2 P°) + Sr (5 s 2 1 S)), which exhibit a significant spin-orbit coupling, we included relativistic effects by two approaches, one applying the Breit-Pauli Hamiltonian to the multireference configuration interaction wave functions, the other combining a spin-orbit Hamiltonian and multireference configuration interaction potential energy curves. Using the results for the relativistic potential energy curves that correspond to the Rb (5 s 2 S) + Sr (5 s 4 d 3 P°) asymptote, we have simulated dispersed fluorescence spectra as they were recently measured in our lab. The comparison with experimental data allows to benchmark both methods and demonstrate that spin-orbit coupling has to be included for the lowest states of RbSr.

Description: Ozonolysis of alkenes, a principle non-photolytic source of atmospheric OH radicals, proceeds through unimolecular decay of energized carbonyl oxide intermediates, known as Criegee intermediates. In this work, cold dimethyl-substituted Criegee intermediates are vibrationally activated in the CH stretch overtone region to drive the 1,4 hydrogen transfer reaction that leads to OH radical products. IR excitation of (CH 3 ) 2 COO reveals the vibrational states with sufficient oscillator strength, coupling to the reaction coordinate, and energy to surmount the effective barrier (≤ 16.0 kcal mol −1 ) to reaction. Insight on the dissociation dynamics is gleaned from homogeneous broadening of the spectral features, indicative of rapid intramolecular vibrational energy redistribution and/or reaction, as well as the quantum state distribution of the OH X 2 Π (v = 0) products. The experimental results are compared with complementary electronic structure calculations, which provide the IR absorption spectrum and geometric changes along the intrinsic reaction coordinate. Additional theoretical analysis reveals the vibrational modes and couplings that permit (CH 3 ) 2 COO to access to the transition state region for reaction. The experimental and theoretical results are compared with an analogous recent study of the IR activation of syn -CH 3 CHOO and its unimolecular decay to OH products [F. Liu, J. M. Beames, A. S. Petit, A. B. McCoy, and M. I. Lester, Science345, 1596 (2014)].

Description: In the present work, the structural and dynamic properties of liquid and undercooled boron are investigated by means of ab initio molecular dynamics simulation. Our results show that both liquid and undercooled states present a well pronounced short-range order (SRO) mainly due to the formation of inverted umbrella structural units. Moreover, we observe the development of a medium-range order (MRO) in the undercooling regime related to the increase of inverted umbrella structural units and of their interconnection as the temperature decreases. We also evidence that this MRO leads to a partial crystallization in the β-rhombohedral crystal below T = 1900 K. Finally, we discuss the role played by the SRO and MRO in the nearly Arrhenius evolution of the diffusion and the non-Arrhenius temperature dependence of the shear viscosity, in agreement with the experiment.

Description: The properties of dipolar cubic lattices are studied and the paradox of how to obtain a long range polarization in such lattices is resolved by choosing a proper shape of the total system. It has been shown that large but finite number of aligned dipoles prefer to exist as needle shaped macroscopic particles [M. Yoon and D. Tománek, J. Phys.: Condens. Matter22, 455105 (2010)]. The total energy for a particle in such a system has one short range contribution depending on the packing (the chosen lattice) and one long range term depending on the dipole density of the system. We show that the latter term corresponds exactly to the polarization term from a dielectric medium embedding a sphere of the considered system. There is no need to include a dielectric medium in this modeling and the “dielectric stabilization” is generated solely by the dipoles of the system.

Description: A practical hyperdynamics method is proposed to accelerate systems with highly endothermic and exothermic reactions such as hydrocarbon pyrolysis and oxidation reactions. In this method, referred to as the “adaptive hyperdynamics (AHD) method,” the bias potential parameters are adaptively updated according to the change in potential energy. The approach is intensively examined for JP-10 (exo-tetrahydrodicyclopentadiene) pyrolysis simulations using the ReaxFF reactive force field. Valid boost parameter ranges are clarified as a result. It is shown that AHD can be used to model pyrolysis at temperatures as low as 1000 K while achieving a boost factor of around 10 5 .

Description: E -crotonic acid was isolated in cryogenic solid N 2 and xenon matrices, and subjected to Laser ultraviolet (UV) and near-infrared (NIR) irradiations. In the deposited matrices, the two low-energy cis C–O E -cc and E -ct conformers, which are the only forms significantly populated in the gas phase, were observed. UV irradiation (λ= 250 nm) of the compound in N 2 matrix allows for experimental detection, not just of the two low-energy cis C–O isomers of Z -crotonic acid previously observed in the experiments carried out in argon matrix [ Z -cc and Z -ct; R. Fausto, A. Kulbida, and O. Schrems, J. Chem. Soc., Faraday Trans.91, 3755–3770 (1995)] but also of the never observed before high-energy forms of both E - and Z -crotonic acids bearing the carboxylic acid group in the trans arrangement ( E -tc and Z -tc conformers). In turn, NIR irradiation experiments in the N 2 matrix allow to produce the high-energy E -tc trans C–O conformer in a selective way, from the initially deposited E -cc form. The vibrational signatures of all the 6 rotameric structures of the crotonic acids experimentally observed, including those of the new trans C–O forms, were determined and the individual spectra fully assigned, also with support of theoretically obtained data. On the other hand, as found before for the compound isolated in argon matrix, the experiments performed in xenon matrix failed to experimental detection of the trans C–O forms. This demonstrates that in noble gas matrices these forms are not stable long enough to allow for their observation by steady state spectroscopy techniques. In these matrices, the trans C–O forms convert spontaneously into their cis C–O counterparts, by tunnelling. Some mechanistic details of the studied processes were extracted and discussed.

Description: We describe a novel approach for the calculation of local electric dipole moments for periodic systems. Since the position operator is ill-defined in periodic systems, maximally localized Wannier functions based on the Berry-phase approach are usually employed for the evaluation of local contributions to the total electric dipole moment of the system. We propose an alternative approach: within a subsystem-density functional theory based embedding scheme, subset electric dipole moments are derived without any additional localization procedure, both for hybrid and non-hybrid exchange–correlation functionals. This opens the way to a computationally efficient evaluation of local electric dipole moments in (molecular) periodic systems as well as their rigorous splitting into atomic electric dipole moments. As examples, Infrared spectra of liquid ethylene carbonate and dimethyl carbonate are presented, which are commonly employed as solvents in Lithium ion batteries.

Description: Long-duration gamma-ray bursts (GRBs) are thought to come from the core collapse of Wolf–Rayet stars. Whereas their stellar masses M * have a rather narrow distribution, the population of GRBs is very diverse, with gamma-ray luminosities L spanning several orders of magnitude. This suggests the existence of a ‘hidden’ stellar variable whose burst-to-burst variation leads to a spread in L . Whatever this hidden variable is, its variation should not noticeably affect the shape of GRB light curves, which display a constant luminosity (in a time-average sense) followed by a sharp drop at the end of the burst seen with Swift /XRT. We argue that such a hidden variable is progenitor star's large-scale magnetic flux. Shortly after the core collapse, most of stellar magnetic flux accumulates near the black hole (BH) and remains there. The flux extracts BH rotational energy and powers jets of roughly a constant luminosity, L j . However, once BH mass accretion rate $\dot{M}$ falls below ~ L j / c 2 , the flux becomes dynamically important and diffuses outwards, with the jet luminosity set by the rapidly declining mass accretion rate, $L_{\rm j}\sim \dot{M}c^2$ . This provides a potential explanation for the sharp end of GRBs and the universal shape of their light curves. During the GRB, gas infall translates spatial variation of stellar magnetic flux into temporal variation of L j . We make use of the deviations from constancy in L j to perform stellar magnetic flux ‘tomography’. Using this method, we infer the presence of magnetized tori in the outer layers of progenitor stars for GRB 920513 and GRB 940210.

Description: We make use of a large set of fast simulations of an intensity mapping experiment with characteristics similar to those expected of the Square Kilometre Array in order to study the viability and limits of blind foreground subtraction techniques. In particular, we consider three different approaches: polynomial fitting, principal component analysis (PCA) and independent component analysis (ICA). We review the motivations and algorithms for the three methods, and show that they can all be described, using the same mathematical framework, as different approaches to the blind source separation problem. We study the efficiency of foreground subtraction both in the angular and radial (frequency) directions, as well as the dependence of this efficiency on different instrumental and modelling parameters. For well-behaved foregrounds and instrumental effects, we find that foreground subtraction can be successful to a reasonable level on most scales of interest. We also quantify the effect that the cleaning has on the recovered signal and power spectra. Interestingly, we find that the three methods yield quantitatively similar results, with PCA and ICA being almost equivalent.

Description: Observed galaxies with high stellar masses or in dense environments have low specific star formation rates, i.e. they are quenched. Based on cosmological hydrodynamic simulations that include a prescription where quenching occurs in regions dominated by hot (〉10 5.4  K) gas, we argue that this hot gas quenching in haloes 〉10 12 M drives both mass quenching (i.e. central quenching) and environment quenching (i.e. satellite quenching). These simulations reproduce a broad range of locally observed trends among quenching, halo mass, stellar mass, environment, and distance to halo centre. Mass quenching is independent of environment because ~10 12 –10 13 M ‘mass quenching haloes’ inhabit a large range of environments. On the other hand, environment quenching is independent of stellar mass because galaxies of all stellar masses may live in dense environments as satellites of groups and clusters. As in observations, the quenched fraction of satellites increases with halo mass and decreases with distance to the centre of the group or cluster. We investigate pre-processing in group haloes, ejected former satellites, and hot gas that extends beyond the virial radius. The agreement of our model with key observational trends suggests that hot gas in massive haloes plays a leading role in quenching low-redshift galaxies.

Description: We used Arecibo Observatory and the Green Bank Telescope to observe OH in 12 early-type galaxies with known reservoirs of dense gas. We present three new detections of OH in absorption in the 1667 MHz line. One objective of our survey was to find evidence of molecular outflows, but our sensitivity and the strength of the OH absorption were insufficient to detect outflows. The detected sources have infrared luminosities and dust temperatures among the lowest of any galaxy detected in OH absorption. The ratio L HCN / L CO , a measure of the dense gas fraction in galaxies, is a powerful selector of OH megamasers for galaxies with high infrared luminosity. In early-type galaxies, which have much lower infrared luminosities, L HCN / L CO is also a promising tool for discovering OH, but in absorption rather than in maser emission. In addition to dense molecular gas, a radio continuum source and a favourable line of sight to the observer are likely key factors in detecting OH absorbers.

Description: The initial mass function (IMF) of the Arches cluster, which was formed a few million years ago in the harsh environment of the Galactic Centre (GC), has long been a target of interest to those who study the GC and the theory of star formation. The distinct star-forming conditions in the GC might have caused the cluster to have a shallower slope or an elevated lower mass cutoff in its IMF. But its mass function (MF) has been revealed only down to 1–2 M (the lower limit of resolved stars), and the low-end MF of the Arches is still unknown. To estimate the unresolved part of the Arches MF, we have devised a novel photometric method that involves the histogram of pixel intensities in the observed image, which contains information on the unresolved, faint stars. By comparing the pixel intensity histograms (PIHs) of numerous artificial images constructed from model IMFs with the observed PIH, we find that the best-fitting model IMF for the Arches cluster has a cutoff mass less than or similar to 0.1 M and a shape very close to that of the Kroupa MF. Our findings imply that the IMF of the Arches cluster is similar to those found in the Galactic disc.

Description: The collision-less transfer of momentum and energy from explosive debris plasma to magnetized background plasma is a salient feature of various astrophysical and space environments. While much theoretical and computational work has investigated collision-less coupling mechanisms and relevant parameters, an experimental validation of the results demands the measurement of the complex, collective electric fields associated with debris-background plasma interaction. Emission spectroscopy offers a non-interfering diagnostic of electric fields via the Stark effect. A unique experiment at the University of California, Los Angeles, that combines the Large Plasma Device (LAPD) and the Phoenix laser facility has investigated the marginally super-Alfvénic, quasi-perpendicular expansion of a laser-produced carbon (C) debris plasma through a preformed, magnetized helium (He) background plasma via emission spectroscopy. Spectral profiles of the He II 468.6 nm line measured at the maximum extent of the diamagnetic cavity are observed to intensify, broaden, and develop equally spaced modulations in response to the explosive C debris, indicative of an energetic electron population and strong oscillatory electric fields. The profiles are analyzed via time-dependent Stark effect models corresponding to single-mode and multi-mode monochromatic (single frequency) electric fields, yielding temporally resolved magnitudes and frequencies. The proximity of the measured frequencies to the expected electron plasma frequency suggests the development of the electron beam-plasma instability, and a simple saturation model demonstrates that the measured magnitudes are feasible provided that a sufficiently fast electron population is generated during C debris–He background interaction. Potential sources of the fast electrons, which likely correspond to collision-less coupling mechanisms, are briefly considered.

Description: It is shown that co-linear injection of electrons or positrons into the wakefield of the self-modulating particle beam is possible and ensures high energy gain. The witness beam must co-propagate with the tail part of the driver, since the plasma wave phase velocity there can exceed the light velocity, which is necessary for efficient acceleration. If the witness beam is many wakefield periods long, then the trapped charge is limited by beam loading effects. The initial trapping is better for positrons, but at the acceleration stage a considerable fraction of positrons is lost from the wave. For efficient trapping of electrons, the plasma boundary must be sharp, with the density transition region shorter than several centimeters. Positrons are not susceptible to the initial plasma density gradient.

Description: The foundational theory for dusty plasmas is the dust charging theory that provides the dust potential and charge arising from the dust interaction with a plasma. The most widely used dust charging theory for negatively charged dust particles is the so-called orbital motion limited (OML) theory, which predicts the dust potential and heat collection accurately for a variety of applications, but was previously found to be incapable of evaluating the dust charge and plasma response in any situation. Here, we report a revised OML formulation that is able to predict the plasma response and hence the dust charge. Numerical solutions of the new OML model show that the widely used Whipple approximation of dust charge-potential relationship agrees with OML theory in the limit of small dust radius compared with plasma Debye length, but incurs large (order-unity) deviation from the OML prediction when the dust size becomes comparable with or larger than plasma Debye length. This latter case is expected for the important application of dust particles in a tokamak plasma.

Description: We report on the controlled removal of an amorphous Se capping layer from Bi 2 Te 3 and Bi 2 Se 3 topological insulators. We show that the Se coalesces into micron-sized islands before desorbing from the surface at a temperature of ∼150 °C. In situ Auger Electron Spectroscopy reveals that Se replaces a significant fraction of the Te near the top surface of the Bi 2 Te 3 . Rutherford Backscattering Spectrometry and Transmission Electron Microscopy show that after heating, Se has been incorporated in the Bi 2 Te 3 lattice down to ∼7 nm from its top surface while remaining iso-structural.

Description: A system of generating and receiving orbital angular momentum (OAM) radio beams, which are collectively formed by two circular array antennas (CAAs) and effectively optimized by two intensity controlled masks, is proposed and experimentally investigated. The scheme is effective in blocking of the unwanted OAM modes and enhancing the power of received radio signals, which results in the capacity gain of system and extended transmission distance of the OAM radio beams. The operation principle of the intensity controlled masks, which can be regarded as both collimator and filter, is feasible and simple to realize. Numerical simulations of intensity and phase distributions at each key cross-sectional plane of the radio beams demonstrate the collimated results. The experimental results match well with the theoretical analysis and the receive distance of the OAM radio beam at radio frequency (RF) 20 GHz is extended up to 200 times of the wavelength of the RF signals, the measured distance is 5 times of the original measured distance. The presented proof-of-concept experiment demonstrates the feasibility of the system.

Description: We demonstrate chemical doping of a topological insulator Bi 2 Se 3 using ion implantation. Ion beam-induced structural damage was characterized using grazing incidence X-ray diffraction and transmission electron microscopy. Ion damage was reversed using a simple thermal annealing step. Carrier-type conversion was achieved using ion implantation followed by an activation anneal in Bi 2 Se 3 thin films. These two sets of experiments establish the feasibility of ion implantation for chemical modification of Bi 2 Se 3 , a prototypical topological insulator. Ion implantation can, in principle, be used for any topological insulator. The direct implantation of dopants should allow better control over carrier concentrations for the purposes of achieving low bulk conductivity. Ion implantation also enables the fabrication of inhomogeneously doped structures, which in turn should make possible new types of device designs.

Description: Scanning near-field photoluminescence (PL) spectroscopy was applied to study spatial variations of emission spectra of Al x Ga 1−x N epilayers with 0.6 ≤ x ≤ 0.7 . PL spectra were found to be spatially uniform with peak wavelength standard deviations of only ∼2 meV and ratios between peak intensity standard deviations and average peak intensity values of 0.06. The observed absence of correlation between the PL peak wavelength and intensity shows that spatial distribution of nonradiative recombination centers is not related to band potential fluctuations. Our results demonstrate that the homogeneous broadening and the random cation distribution primarily determine PL linewidths for layers grown under optimized conditions.

Description: Using the ideal magnetohydrodynamic model, we calculate the temporal evolution of initial ripples on the boundaries of a planar plasma slab that is subjected to the magneto-Rayleigh-Taylor instability. The plasma slab consists of three regions. We assume that in each region the plasma density is constant with an arbitrary value and the magnetic field is also constant with an arbitrary magnitude and an arbitrary direction parallel to the interfaces. Thus, the instability may be driven by a combination of magnetic pressure and kinetic pressure. The general dispersion relation is derived, together with the feedthrough factor between the two interfaces. The temporal evolution is constructed from the superposition of the eigenmodes. Previously established results are recovered in the various limits. Numerical examples are given on the temporal evolution of ripples on the interfaces of the finite plasma slab.

Description: We demonstrate direct evidence that the strain variation induced by local lattice distortion exists in the surface layers of SnO 2 nanowires by coupled scanning transmission electron microscopy and digital image correlation techniques. First-principles calculations suggest that surface reduction and subsurface oxygen vacancies account for such vigorous wavelike strain. Our study revealed that the localized change of surface atomistic configuration was responsible for the observed reduction of elastic modulus and hardness of SnO 2 nanowires, as well as the superior sensing properties of SnO 2 nanowire network.

Description: Complementary logic based on tunnel field-effect transistors (TFETs) would drastically reduce power consumption thanks to the TFET's potential to obtain a sub-60 mV/dec subthreshold swing (SS). However, p-type TFETs typically do not meet the performance of n-TFETs for direct bandgap III-V configurations. The p-TFET SS stays well above 60 mV/dec, due to the low density of states in the conduction band. We therefore propose a source configuration in which a highly doped region is maintained only near the tunnel junction. In the remaining part of the source, the hot carriers in the exponential tail of the Fermi-Dirac distribution are blocked by reducing the doping degeneracy, either with a source section with a lower doping concentration or with a heterostructure. We apply this concept to n-p-i-p configurations consisting of In 0.53 Ga 0.47 As and an InP-InAs heterostructure. 15-band quantum mechanical simulations predict that the configurations with our source design can obtain sub-60 mV/dec SS, with an on-current comparable to the conventional source design.

Description: We design and experimentally realize an ultra-broad band metamaterial-based acoustic absorption material. Unlike traditional acoustic absorbers, the designed device features a simple configuration unrestricted by the material type and does not require extra sound-absorbing materials, suggesting the potential to have simultaneously structural-stiffness and environmental-friendliness. Analytical analyses are provided to explain such distinct characteristics, which are revealed to stem from the localization and dissipation of waves with different frequencies at particular spatial positions. This is also demonstrated both numerically and experimentally. Our results may offer possible designs for various applications such as noise reduction and making underwater anechoic materials.

Description: The absolute photoluminescence (PL) quantum yield (QY) of multilayers of Silicon nanocrystals (SiNCs) separated by SiO 2 barriers were thoroughly studied as function of the barrier thickness, excitation wavelength, and temperature. By mastering the plasma-enhanced chemical vapor deposition growth, we produce a series of samples with the same size-distribution of SiNCs but variable interlayer barrier distance. These samples enable us to clearly demonstrate that the increase of barrier thickness from ∼1 to larger than 2 nm induces doubling of the PL QY value, which corresponds to the change of number of close neighbors in the hcp structure. The temperature dependence of PL QY suggests that the PL QY changes are due to a thermally activated transport of excitation into non-radiative centers in dark NCs or in the matrix. We estimate that dark NCs represent about 68% of the ensemble of NCs. The PL QY excitation spectra show no significant changes upon changing the barrier thickness and no clear carrier multiplication effects. The dominant effect is the gradual decrease of the PL QY with increasing excitation photon energy.

Description: High-quality BN-Graphene-BN nanoribbon capacitors with double side-gates of graphene have been experimentally realized. The double side-gates can effectively modulate the electronic properties of graphene nanoribbon capacitors. By applying anti-symmetric side-gate voltages, we observed significant upward shifting and flattening of the V-shaped capacitance curve near the charge neutrality point. Symmetric side-gate voltages, however, only resulted in tilted upward shifting along the opposite direction of applied gate voltages. These modulation effects followed the behavior of graphene nanoribbons predicted theoretically for metallic side-gate modulation. The negative quantum capacitance phenomenon predicted by numerical simulations for graphene nanoribbons modulated by graphene side-gates was not observed, possibly due to the weakened interactions between the graphene nanoribbon and side-gate electrodes caused by the Ga + beam etching process.

Description: GaAs/Al-GaAs core-shell nanowires fabricated by molecular beam epitaxy contain quantum confining structures susceptible of producing narrow photoluminescence (PL) and single photons. The nanoscale chemical mapping of these structures is analyzed in 3D by atom probe tomography (APT). The study allows us to confirm that Al atoms tend to segregate within the AlGaAs shells towards the vertices of the hexagons defining the nanowire cross section. We also find strong alloy fluctuations remaining AlGaAs shell, leading occasionally to the formation of quantum dots (QDs). The PL emission energies predicted in the framework of a 3D effective mass model for a QD analyzed by APT and the PL spectra measured on other nanowires from the same growth batch are consistent within the experimental uncertainties.

Description: Rectilinear deposition of elongated DNA molecules was achieved by the forced dewetting of a DNA solution droplet over a nanograting. Uncoiling of double stranded DNA is made by the conjunction of both DNA terminal anchoring on a functionalized substrate and capillary force acting throughout the forced dewetting of a DNA solution droplet. The deposition over a nanograting allows the molecule to be uncoiled on the edges of the grooves and to maintain a rectilinear conformation. This DNA deposition technique uses transparent nanograting obtained by laser interference lithography and has been developed for the specific need in observation dsDNA molecules in extended conformation.

Description: Perfectly lattice-matched magnetic tunnel junctions (MTJs) consisting of a Heusler alloy B 2-Co 2 FeAl (CFA) electrode and a cation-disorder spinel (Mg-Al-O) barrier were fabricated by sputtering and plasma oxidation. We achieved a large tunnel magnetoresistance (TMR) ratio of 228% at room temperature (RT) (398% at 5 K) for the epitaxial CFA/MgAl-O/CoFe(001) MTJ, in which the effect of lattice defects on TMR ratios is excluded. With inserting a ultrathin (≤1.5 nm) CoFe layer between the CFA and Mg-Al-O, the TMR ratio further increased up to 280% at RT (453% at 5 K), which reflected the importance of controlling barrier-electrode interface states other than the lattice matching.

Description: A major goal of nanotechnology is to develop the capability to arrange matter at will by placing individual atoms at desired locations in a predetermined configuration to build a nanostructure with specific properties or function. The scanning tunneling microscope has demonstrated the ability to arrange the basic building blocks of matter, single atoms, in two-dimensional configurations. An array of various nanostructures has been assembled, which display the quantum mechanics of quantum confined geometries. The level of human interaction needed to physically locate the atom and bring it to the desired location limits this atom assembly technology. Here we report the use of autonomous atom assembly via path planning technology; this allows atomically perfect nanostructures to be assembled without the need for human intervention, resulting in precise constructions in shorter times. We demonstrate autonomous assembly by assembling various quantum confinement geometries using atoms and molecules and describe the benefits of this approach.

Description: We developed a silicon avalanche photodiode (Si-APD) linear-array detector for use in nuclear resonant scattering experiments using synchrotron X-rays. The Si-APD linear array consists of 64 pixels (pixel size: 100 × 200 μm 2 ) with a pixel pitch of 150 μm and depletion depth of 10 μm. An ultrafast frontend circuit allows the X-ray detector to obtain a high output rate of 〉10 7 cps per pixel. High-performance integrated circuits achieve multichannel scaling over 1024 continuous time bins with a 1 ns resolution for each pixel without dead time. The multichannel scaling method enabled us to record a time spectrum of the 14.4 keV nuclear radiation at each pixel with a time resolution of 1.4 ns (FWHM). This method was successfully applied to nuclear forward scattering and nuclear small-angle scattering on 57 Fe.

Description: We describe a nanosecond time-resolved fluorescence spectrometer that acquires fluorescence decay waveforms from each well of a 384-well microplate in 3 min with signal-to-noise exceeding 400 using direct waveform recording. The instrument combines high-energy pulsed laser sources (5–10 kHz repetition rate) with a photomultiplier and high-speed digitizer (1 GHz) to record a fluorescence decay waveform after each pulse. Waveforms acquired from rhodamine or 5-((2-aminoethyl)amino) naphthalene-1-sulfonic acid dyes in a 384-well plate gave lifetime measurements 5- to 25-fold more precise than the simultaneous intensity measurements. Lifetimes as short as 0.04 ns were acquired by interleaving with an effective sample rate of 5 GHz. Lifetime measurements resolved mixtures of single-exponential dyes with better than 1% accuracy. The fluorescence lifetime plate reader enables multiple-well fluorescence lifetime measurements with an acquisition time of 0.5 s per well, suitable for high-throughput fluorescence lifetime screening applications.

Description: Diagnostic for investigating and distinguishing different laser ion acceleration mechanisms has been developed and successfully tested. An ion separation wide angle spectrometer can simultaneously investigate three important aspects of the laser plasma interaction: (1) acquire angularly resolved energy spectra for two ion species, (2) obtain ion energy spectra for multiple species, separated according to their charge to mass ratio, along selected axes, and (3) collect laser radiation reflected from and transmitted through the target and propagating in the same direction as the ion beam. Thus, the presented diagnostic constitutes a highly adaptable tool for accurately studying novel acceleration mechanisms in terms of their angular energy distribution, conversion efficiency, and plasma density evolution.

Description: Motivation: Mapping of high-throughput sequencing data and other bulk sequence comparison applications have motivated a search for high-efficiency sequence alignment algorithms. The bit-parallel approach represents individual cells in an alignment scoring matrix as bits in computer words and emulates the calculation of scores by a series of logic operations composed of AND, OR, XOR, complement, shift and addition. Bit-parallelism has been successfully applied to the longest common subsequence (LCS) and edit-distance problems, producing fast algorithms in practice. Results: We have developed BitPAl, a bit-parallel algorithm for general, integer-scoring global alignment. Integer-scoring schemes assign integer weights for match, mismatch and insertion/deletion. The BitPAl method uses structural properties in the relationship between adjacent scores in the scoring matrix to construct classes of efficient algorithms, each designed for a particular set of weights. In timed tests, we show that BitPAl runs 7–25 times faster than a standard iterative algorithm. Availability and implementation: Source code is freely available for download at http://lobstah.bu.edu/BitPAl/BitPAl.html . BitPAl is implemented in C and runs on all major operating systems. Contact : jloving@bu.edu or yhernand@bu.edu or gbenson@bu.edu Supplementary information : Supplementary data are available at Bioinformatics online.

Description: : Next-generation sequencing (NGS) has a large potential in HIV diagnostics, and genotypic prediction models have been developed and successfully tested in the recent years. However, albeit being highly accurate, these computational models lack computational efficiency to reach their full potential. In this study, we demonstrate the use of graphics processing units (GPUs) in combination with a computational prediction model for HIV tropism. Our new model named gCUP, parallelized and optimized for GPU, is highly accurate and can classify 〉175 000 sequences per second on an NVIDIA GeForce GTX 460. The computational efficiency of our new model is the next step to enable NGS technologies to reach clinical significance in HIV diagnostics. Moreover, our approach is not limited to HIV tropism prediction, but can also be easily adapted to other settings, e.g. drug resistance prediction. Availability and implementation: The source code can be downloaded at http://www.heiderlab.de Contact: d.heider@wz-straubing.de

Description: : We present a new method to incrementally construct the FM-index for both short and long sequence reads, up to the size of a genome. It is the first algorithm that can build the index while implicitly sorting the sequences in the reverse (complement) lexicographical order without a separate sorting step. The implementation is among the fastest for indexing short reads and the only one that practically works for reads of averaged kilobases in length. Availability and implementation: https://github.com/lh3/ropebwt2 Contact: hengli@broadinstitute.org

Description: : AliView is an alignment viewer and editor designed to meet the requirements of next-generation sequencing era phylogenetic datasets. AliView handles alignments of unlimited size in the formats most commonly used, i.e. FASTA, Phylip, Nexus, Clustal and MSF. The intuitive graphical interface makes it easy to inspect, sort, delete, merge and realign sequences as part of the manual filtering process of large datasets. AliView also works as an easy-to-use alignment editor for small as well as large datasets. Availability and implementation: AliView is released as open-source software under the GNU General Public License, version 3.0 (GPLv3), and is available at GitHub ( www.github.com/AliView ). The program is cross-platform and extensively tested on Linux, Mac OS X and Windows systems. Downloads and help are available at http://ormbunkar.se/aliview Contact: anders.larsson@ebc.uu.se Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: The ability to accurately read the order of nucleotides in DNA and RNA is fundamental for modern biology. Errors in next-generation sequencing can lead to many artifacts, from erroneous genome assemblies to mistaken inferences about RNA editing. Uneven coverage in datasets also contributes to false corrections. Result: We introduce Trowel, a massively parallelized and highly efficient error correction module for Illumina read data. Trowel both corrects erroneous base calls and boosts base qualities based on the k -mer spectrum. With high-quality k -mers and relevant base information, Trowel achieves high accuracy for different short read sequencing applications.The latency in the data path has been significantly reduced because of efficient data access and data structures. In performance evaluations, Trowel was highly competitive with other tools regardless of coverage, genome size read length and fragment size. Availability and implementation: Trowel is written in C++ and is provided under the General Public License v3.0 (GPLv3). It is available at http://trowel-ec.sourceforge.net . Contact: euncheon.lim@tue.mpg.de or weigel@tue.mpg.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : The application of protein–protein docking in large-scale interactome analysis is a major challenge in structural bioinformatics and requires huge computing resources. In this work, we present MEGADOCK 4.0, an FFT-based docking software that makes extensive use of recent heterogeneous supercomputers and shows powerful, scalable performance of 〉97% strong scaling. Availability and Implementation: MEGADOCK 4.0 is written in C++ with OpenMPI and NVIDIA CUDA 5.0 (or later) and is freely available to all academic and non-profit users at: http://www.bi.cs.titech.ac.jp/megadock . Contact: akiyama@cs.titech.ac.jp Supplementary information: Supplementary data are available at Bioinformatics online

Description: Motivation: The identification of active transcriptional regulatory elements is crucial to understand regulatory networks driving cellular processes such as cell development and the onset of diseases. It has recently been shown that chromatin structure information, such as DNase I hypersensitivity (DHS) or histone modifications, significantly improves cell-specific predictions of transcription factor binding sites. However, no method has so far successfully combined both DHS and histone modification data to perform active binding site prediction. Results: We propose here a method based on hidden Markov models to integrate DHS and histone modifications occupancy for the detection of open chromatin regions and active binding sites. We have created a framework that includes treatment of genomic signals, model training and genome-wide application. In a comparative analysis, our method obtained a good trade-off between sensitivity versus specificity and superior area under the curve statistics than competing methods. Moreover, our technique does not require further training or sequence information to generate binding location predictions. Therefore, the method can be easily applied on new cell types and allow flexible downstream analysis such as de novo motif finding. Availability and implementation: Our framework is available as part of the Regulatory Genomics Toolbox. The software information and all benchmarking data are available at http://costalab.org/wp/dh-hmm . Contact: ivan.costa@rwth-aachen.de or eduardo.gusmao@rwth-aachen.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: A proper target or marker is essential in any diagnosis (e.g. an infection or cancer). An ideal diagnostic target should be both conserved in and unique to the pathogen. Currently, these targets can only be identified manually, which is time-consuming and usually error-prone. Because of the increasingly frequent occurrences of emerging epidemics and multidrug-resistant ‘superbugs’, a rapid diagnostic target identification process is needed. Results: A new method that can identify uniquely conserved regions (UCRs) as candidate diagnostic targets for a selected group of organisms solely from their genomic sequences has been developed and successfully tested. Using a sequence-indexing algorithm to identify UCRs and a k -mer integer-mapping model for computational efficiency, this method has successfully identified UCRs within the bacteria domain for 15 test groups, including pathogenic, probiotic, commensal and extremophilic bacterial species or strains. Based on the identified UCRs, new diagnostic primer sets were designed, and their specificity and efficiency were tested by polymerase chain reaction amplifications from both pure isolates and samples containing mixed cultures. Availability and implementation: The UCRs identified for the 15 bacterial species are now freely available at http://ucr.synblex.com . The source code of the programs used in this study is accessible at http://ucr.synblex.com/bacterialIdSourceCode.d.zip Contact: yazhousun@synblex.com Supplementary Information: Supplementary data are available at Bioinformatics online.

Description: Motivation: A popular method for classification of protein domain movements apportions them into two main types: those with a ‘hinge’ mechanism and those with a ‘shear’ mechanism. The intuitive assignment of domain movements to these classes has limited the number of domain movements that can be classified in this way. Furthermore, whether intended or not, the term ‘shear’ is often interpreted to mean a relative translation of the domains. Results: Numbers of occurrences of four different types of residue contact changes between domains were optimally combined by logistic regression using the training set of domain movements intuitively classified as hinge and shear to produce a predictor for hinge and shear. This predictor was applied to give a 10-fold increase in the number of examples over the number previously available with a high degree of precision. It is shown that overall a relative translation of domains is rare, and that there is no difference between hinge and shear mechanisms in this respect. However, the shear set contains significantly more examples of domains having a relative twisting movement than the hinge set. The angle of rotation is also shown to be a good discriminator between the two mechanisms. Availability and implementation: Results are free to browse at http://www.cmp.uea.ac.uk/dyndom/interface/ . Contact: sjh@cmp.uea.ac.uk . Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Recent studies on human disease have revealed that aberrant interaction between proteins probably underlies a substantial number of human genetic diseases. This suggests a need to investigate disease inheritance mode using interaction, and based on which to refresh our conceptual understanding of a series of properties regarding inheritance mode of human disease. Results: We observed a strong correlation between the number of protein interactions and the likelihood of a gene causing any dominant diseases or multiple dominant diseases, whereas no correlation was observed between protein interaction and the likelihood of a gene causing recessive diseases. We found that dominant diseases are more likely to be associated with disruption of important interactions. These suggest inheritance mode should be understood using protein interaction. We therefore reviewed the previous studies and refined an interaction model of inheritance mode, and then confirmed that this model is largely reasonable using new evidences. With these findings, we found that the inheritance mode of human genetic diseases can be predicted using protein interaction. By integrating the systems biology perspectives with the classical disease genetics paradigm, our study provides some new insights into genotype–phenotype correlations. Contact: haodapeng@ems.hrbmu.edu.cn or biofomeng@hotmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : Recently, several high profile studies collected cell viability data from panels of cancer cell lines treated with many drugs applied at different concentrations. Such drug sensitivity data for cancer cell lines provide suggestive treatments for different types and subtypes of cancer. Visualization of these datasets can reveal patterns that may not be obvious by examining the data without such efforts. Here we introduce Drug/Cell-line Browser (DCB), an online interactive HTML5 data visualization tool for interacting with three of the recently published datasets of cancer cell lines/drug-viability studies. DCB uses clustering and canvas visualization of the drugs and the cell lines, as well as a bar graph that summarizes drug effectiveness for the tissue of origin or the cancer subtypes for single or multiple drugs. DCB can help in understanding drug response patterns and prioritizing drug/cancer cell line interactions by tissue of origin or cancer subtype. Availability and implementation: DCB is an open source Web-based tool that is freely available at: http://www.maayanlab.net/LINCS/DCB Contact: avi.maayan@mssm.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: As an important subtype of structural variations, chromosomal translocation is associated with various diseases, especially cancers, by disrupting gene structures and functions. Traditional methods for identifying translocations are time consuming and have limited resolutions. Recently, a few studies have employed next-generation sequencing (NGS) technology for characterizing chromosomal translocations on human genome, obtaining high-throughput results with high resolutions. However, these studies are mainly focused on mechanism-specific or site-specific translocation mapping. In this study, we conducted a comprehensive genome-wide analysis on the characterization of human chromosomal material exchange with regard to the chromosome translocations. Using NGS data of 1,481 subjects from the 1000 Genomes Project, we identified 15,349,092 translocated DNA fragment pairs, ranging from 65 to 1,886 bp and with an average size of approximately 102 bp. On average, each individual genome carried about 10,364 pairs, covering approximately 0.069% of the genome. We identified 16 translocation hot regions, among which two regions did not contain repetitive fragments. Results of our study overlapped with a majority of previous results, containing approximately 79% of approximately 2,340 translocations characterized in three available translocation databases. In addition, our study identified five novel potential recurrent chromosomal material exchange regions with greater than 20% detection rates. Our results will be helpful for an accurate characterization of translocations in human genomes, and contribute as a resource for future studies of the roles of translocations in human disease etiology and mechanisms.

Description: Fractured rocks are known to exhibit seismic anisotropy and shear wave splitting (SWS). SWS is commonly used for fractured rock characterization and has been shown to be sensitive to fluid type. The presence of partial liquid/gas saturation is also known to affect the elastic properties of rocks. The combined effect of both fractures and partial liquid/gas saturation is still unknown. Using synthetic, silica-cemented sandstones with aligned penny-shaped voids, we conducted laboratory ultrasonic experiments to investigate the effect fractures aligned at an oblique angle to wave propagation would have on SWS under partial liquid/gas saturation conditions. The result for the fractured rock shows a saturation dependence which can be explained by combining a fractured rock model and a partial saturation model. At high to full water saturation values, SWS decreases as a result of the fluid bulk modulus effect on the quasi-shear wave. This bulk modulus effect is frequency dependent as a result of wave-induced fluid flow mechanisms, which would in turn lead to frequency dependent SWS. This result suggests the possible use of SWS for discriminating between full liquid saturation and partial liquid/gas saturation.

Description: The effect of network density and geometric distribution on kinematic non-linear source inversion is investigated by inverting synthetic ground motions from a buried strike-slip fault ( M w 6.5), that have been generated by dynamic spontaneous rupture modelling. For the inversion, we use a physics-based regularized Yoffe function as slip velocity function. We test three different cases of station network geometry: (i) single station, varying azimuth and epicentral distance; (ii) multistation circular configurations, that is stations at similar distances from the fault, and regularly spaced around the fault; (iii) irregular multistation configurations using different numbers of stations. Our results show: (1) single station tests suggest that it may be possible to obtain a relatively good source model even using a single station. The best source model using a single station is obtained with stations at which amplitude ratios between three components are not large. We infer that both azimuthal angle and source-to-station distance play an important role in the design of optimal seismic network for source inversion. (2) Multistation tests show that the quality of the inverted source systematically correlates neither with the number of stations, nor with waveform misfit. (3) Waveform misfit has a direct correlation with the number of stations, resulting in overfitting the observed data without any systematic improvement of the source. It suggests that the best source model is not necessarily derived from the model with minimum waveform misfit. (4) A seismic network with a small number of well-spaced stations around the fault may be sufficient to obtain acceptable source inversion.

Description: Fermat's interferometric principle is used to compute interior transmission traveltimes pq from exterior transmission traveltimes sp and sq . Here, the exterior traveltimes are computed for sources s on a boundary B that encloses a volume V of interior points p and q . Once the exterior traveltimes are computed, no further ray tracing is needed to calculate the interior times pq . Therefore this interferometric approach can be more efficient than explicitly computing interior traveltimes pq by ray tracing. Moreover, the memory requirement of the traveltimes is reduced by one dimension, because the boundary B is of one fewer dimension than the volume V . An application of this approach is demonstrated with interbed multiple (IM) elimination. Here, the IMs in the observed data are predicted from the migration image and are subsequently removed by adaptive subtraction. This prediction is enabled by the knowledge of interior transmission traveltimes pq computed according to Fermat's interferometric principle. We denote this principle as the ‘traveltime holographic principle’, by analogy with the holographic principle in cosmology where information in a volume is encoded on the region's boundary.

Description: Earthquake hypocentre locations are crucial in many domains of application (academic and industrial) as seismic event location maps are commonly used to delineate faults or fractures. The interpretation of these maps depends on location accuracy and on the reliability of the associated uncertainties. The largest contribution to location and uncertainty errors is due to the fact that the velocity model errors are usually not correctly taken into account. We propose a new Bayesian formulation that integrates properly the knowledge on the velocity model into the formulation of the probabilistic earthquake location. In this work, the velocity model uncertainties are first estimated with a Bayesian tomography of active shot data. We implement a sampling Monte Carlo type algorithm to generate velocity models distributed according to the posterior distribution. In a second step, we propagate the velocity model uncertainties to the seismic event location in a probabilistic framework. This enables to obtain more reliable hypocentre locations as well as their associated uncertainties accounting for picking and velocity model uncertainties. We illustrate the tomography results and the gain in accuracy of earthquake location for two synthetic examples and one real data case study in the context of induced microseismicity.

Description: Clinical and neuropathological similarities between dementia with Lewy bodies (DLB), Parkinson's and Alzheimer's diseases (PD and AD, respectively) suggest that these disorders may share etiology. To test this hypothesis, we have performed an association study of 54 genomic regions, previously implicated in PD or AD, in a large cohort of DLB cases and controls. The cohort comprised 788 DLB cases and 2624 controls. To minimize the issue of potential misdiagnosis, we have also performed the analysis including only neuropathologically proven DLB cases (667 cases). The results show that the APOE is a strong genetic risk factor for DLB, confirming previous findings, and that the SNCA and SCARB2 loci are also associated after a study-wise Bonferroni correction, although these have a different association profile than the associations reported for the same loci in PD. We have previously shown that the p.N370S variant in GBA is associated with DLB, which, together with the findings at the SCARB2 locus, suggests a role for lysosomal dysfunction in this disease. These results indicate that DLB has a unique genetic risk profile when compared with the two most common neurodegenerative diseases and that the lysosome may play an important role in the etiology of this disorder. We make all these data available.

Description: Mitochondrial DNA mutations at MT-ATP6 gene are relatively common in individuals suffering from striatal necrosis syndromes. These patients usually do not show apparent histochemical and/or biochemical signs of oxidative phosphorylation dysfunction. Because of this, MT-ATP6 is not typically analyzed in many other mitochondrial disorders that have not been previously associated to mutations in this gene. To correct this bias, we have performed a screening of the MT-ATP6 gene in a large collection of patients suspected of suffering different mitochondrial DNA (mtDNA) disorders. In three cases, biochemical, molecular-genetics and other analyses in patient tissues and cybrids were also carried out. We found three new pathologic mutations. Two of them in patients showing phenotypes that have not been commonly associated to mutations in the MT-ATP6 gene. These results remark the importance of sequencing the MT-ATP6 gene in patients with striatal necrosis syndromes, but also within other mitochondrial pathologies. This gene should be sequenced at least in all those patients suspected of suffering an mtDNA disorder disclosing normal results for histochemical and biochemical analyses of respiratory chain.

Description: Immunoglobulin-like domain containing receptor 1 ( ILDR1 ) is a poorly characterized gene that was first identified in lymphoma cells. Recently, ILDR1 has been found to be responsible for autosomal recessive hearing impairment DFNB42. Patients with ILDR1 mutations cause bilateral non-progressive moderate-to-profound sensorineural hearing impairment. However, the etiology and mechanism of ILDR1 -related hearing loss remains to be elucidated. In order to uncover the pathology of DFNB42 deafness, we used the morpholino injection technique to establish an ildr1b -morphant zebrafish model. Ildr1b -morphant zebrafish displayed defective hearing and imbalanced swimming, and developmental delays were seen in the semicircular canals of the inner ear. The gene expression profile and real-time PCR revealed down-regulation of atp1b2b (encoding Na + /K + transporting, beta 2b polypeptide) in ildr1b -morphant zebrafish. We found that injection of atp1b2b mRNA into ildr1b -knockdown zebrafish could rescue the phenotype of developmental delay of the semicircular canals. Moreover, ildr1b -morphant zebrafish had reduced numbers of lateral line neuromasts due to the disruption of lateral line primordium migration. In situ hybridization showed the involvement of attenuated FGF signaling and the chemokine receptor 4b ( cxcr4b ) and chemokine receptor 7b ( cxcr7b ) in posterior lateral line primordium of ildr1b -morphant zebrafish. We concluded that Ildr1b is crucial for the development of the inner ear and the lateral line system. This study provides the first evidence for the mechanism of Ildr1b on hearing in vivo and sheds light on the pathology of DFNB42.