Equifinality and preservation potential of complex eskers

Eskers are useful for reconstructing meltwater drainage systems of glaciers and ice sheets. However, our process understanding of eskers suffers from a disconnect between sporadic detailed morpho‐sedimentary investigations of abundant large‐scale ancient esker systems, and a small number of modern analogues where esker formation has been observed. This paper presents the results of detailed field and high‐resolution remote sensing studies into two esker systems that have recently emerged at Hørbyebreen, Svalbard, and one at Breiðamerkurjökull, Iceland. Despite the different glaciological settings (polythermal valley glacier vs. active temperate piedmont lobe), in all cases a distinctive planform morphology has developed, where ridges are orientated in two dominant directions corresponding to the direction of ice flow and the shape of the ice margin. These two orientations in combination form a cross‐cutting and locally rectilinear pattern. One set of ridges at Hørbyebreen is a hybrid of eskers and geometric ridges formed during a surge and/or jökulhlaup event. The other sets of ridges are eskers formed time‐transgressively at a retreating ice margin. The similar morphology of esker complexes formed in different ways on both glacier forelands implies equifinality, meaning that care should be taken when interpreting Quaternary esker patterns. The eskers at Hørbyebreen contain substantial ice‐cores with a high ice:sediment ratio, suggesting that they would be unlikely to survive after ice melt. The Breiðamerkurjökull eskers emerged from terrain characterized by buried ice that has melted out. Our observations lead us to conclude that eskers may reflect a wide range of processes at dynamic ice margins, including significant paraglacial adjustments. This work, as well as previous studies, confirms that constraints on esker morphology include: topographic setting (e.g. confined valley or broad plain); sediment and meltwater availability (including surges and jökulhlaups); position of formation (supraglacial, englacial or subglacial); and ice‐marginal dynamics such as channel abandonment, the formation of outwash heads or the burial and/or exhumation of dead ice

Physico-chemical foundations underpinning microarray and next-generation sequencing experiments

Hybridization of nucleic acids on solid surfaces is a key process involved in high-throughput technologies such as microarrays and, in some cases, next-generation sequencing (NGS). A physical understanding of the hybridization process helps to determine the accuracy of these technologies. The goal of a widespread research program is to develop reliable transformations between the raw signals reported by the technologies and individual molecular concentrations from an ensemble of nucleic acids. This research has inputs from many areas, from bioinformatics and biostatistics, to theoretical and experimental biochemistry and biophysics, to computer simulations. A group of leading researchers met in Ploen Germany in 2011 to discuss present knowledge and limitations of our physico-chemical understanding of high-throughput nucleic acid technologies. This meeting inspired us to write this summary, which provides an overview of the state-of-the-art approaches based on physico-chemical foundation to modeling of the nucleic acids hybridization process on solid surfaces. In addition, practical application of current knowledge is emphasized

Sequencing by Hybridization of Long Targets

Sequencing by Hybridization (SBH) reconstructs an n-long target DNA sequence from its biochemically determined l-long subsequences. In the standard approach, the length of a uniformly random sequence that can be unambiguously reconstructed is limited to due to repetitive subsequences causing reconstruction degeneracies. We present a modified sequencing method that overcomes this limitation without the need for different types of biochemical assays and is robust to error

HMMSplicer: A Tool for Efficient and Sensitive Discovery of Known and Novel Splice Junctions in RNA-Seq Data

Background: High-throughput sequencing of an organism’s transcriptome, or RNA-Seq, is a valuable and versatile new strategy for capturing snapshots of gene expression. However, transcriptome sequencing creates a new class of alignment problem: mapping short reads that span exon-exon junctions back to the reference genome, especially in the case where a splice junction is previously unknown. Methodology/Principal Findings: Here we introduce HMMSplicer, an accurate and efficient algorithm for discovering canonical and non-canonical splice junctions in short read datasets. HMMSplicer identifies more splice junctions than currently available algorithms when tested on publicly available A. thaliana, P. falciparum, and H. sapiens datasets without a reduction in specificity. Conclusions/Significance: HMMSplicer was found to perform especially well in compact genomes and on genes with low expression levels, alternative splice isoforms, or non-canonical splice junctions. Because HHMSplicer does not rely on prebuilt gene models, the products of inexact splicing are also detected. For H. sapiens, we find 3.6 % of 39 splice sites and 1.4% of 59 splice sites are inexact, typically differing by 3 bases in either direction. In addition, HMMSplicer provides a score for every predicted junction allowing the user to set a threshold to tune false positive rates depending on the needs of the experiment. HMMSplicer is implemented in Python. Code and documentation are freely available a

Anatomy of terminal moraine segments and implied lake stability on Ngozumpa Glacier, Nepal, from electrical resistivity tomography (ERT)

This research was supported financially by the European Commission FP7-MC-IEF (PIEF-GA-2012-330805), the University Centre in Svalbard (UNIS), National Geographic Society GRANT #W135-10.Moraine-dammed lakes at debris-covered glaciers are becoming increasingly common and pose significant outburst flood hazards if the dam is breached. While moraine subsurface structure and internal processes are likely to influence dam stability, only few sites have so far been investigated. We conducted electrical resistivity tomography (ERT) surveys at two sites on the terminal moraine complex of the Ngozumpa Glacier, Nepal, to aid assessment of future terminus stability. The resistivity signature of glacier ice at the site (100-15 kΩ m) is more consistent with values measured from cold glacier ice and while this may be feasible, uncertainties in the data inversion introduce ambiguity to this thermal interpretation. However, the ERT data does provide a significant improvement to our knowledge of the subsurface characteristics at these sites, clearly showing the presence (or absence) of glacier ice. Our interpretation is that of a highly complex latero-terminal moraine, resulting from interaction between previous glacier advance, recession and outburst flooding. If the base-level Spillway Lake continues to expand to a fully formed moraine-dammed glacial lake, the degradation of the ice core could have implications for glacial lake outburst risk.Publisher PDFPeer reviewe

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038−2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pKa values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pKa values in closer agreement to those obtained experimentally for the more slowly converging residues

Accurate Estimates of Microarray Target Concentration from a Simple Sequence-Independent Langmuir Model

Background: Microarray technology is a commonly used tool for assessing global gene expression. Many models for estimation of target concentration based on observed microarray signal have been proposed, but, in general, these models have been complex and platform-dependent. Principal Findings: We introduce a universal Langmuir model for estimation of absolute target concentration from microarray experiments. We find that this sequence-independent model, characterized by only three free parameters, yields excellent predictions for four microarray platforms, including Affymetrix, Agilent, Illumina and a custom-printed microarray. The model also accurately predicts concentration for the MAQC data sets. This approach significantly reduces the computational complexity of quantitative target concentration estimates. Conclusions: Using a simple form of the Langmuir isotherm model, with a minimum of parameters and assumptions, and without explicit modeling of individual probe properties, we were able to recover absolute transcript concentrations with high R 2 on four different array platforms. The results obtained here suggest that with a ‘‘spiked-in’ ’ concentration serie

Distinct Cytoplasmic and Nuclear Functions of the Stress Induced Protein DDIT3/CHOP/GADD153

DDIT3, also known as GADD153 or CHOP, encodes a basic leucine zipper transcription factor of the dimer forming C/EBP family. DDIT3 is known as a key regulator of cellular stress response, but its target genes and functions are not well characterized. Here, we applied a genome wide microarray based expression analysis to identify DDIT3 target genes and functions. By analyzing cells carrying tamoxifen inducible DDIT3 expression constructs we show distinct gene expression profiles for cells with cytoplasmic and nuclear localized DDIT3. Of 175 target genes identified only 3 were regulated by DDIT3 in both cellular localizations. More than two thirds of the genes were downregulated, supporting a role for DDIT3 as a dominant negative factor that could act by either cytoplasmic or nuclear sequestration of dimer forming transcription factor partners. Functional annotation of target genes showed cell migration, proliferation and apoptosis/survival as the most affected categories. Cytoplasmic DDIT3 affected more migration associated genes, while nuclear DDIT3 regulated more cell cycle controlling genes. Cell culture experiments confirmed that cytoplasmic DDIT3 inhibited migration, while nuclear DDIT3 caused a G1 cell cycle arrest. Promoters of target genes showed no common sequence motifs, reflecting that DDIT3 forms heterodimers with several alternative transcription factors that bind to different motifs. We conclude that expression of cytoplasmic DDIT3 regulated 94 genes. Nuclear translocation of DDIT3 regulated 81 additional genes linked to functions already affected by cytoplasmic DDIT3. Characterization of DDIT3 regulated functions helps understanding its role in stress response and involvement in cancer and degenerative disorders

Switch from antagonist to agonist of the androgen receptor blocker bicalutamide is associated with prostate tumour progression in a new model system

Advanced prostate cancer is treated by androgen ablation and/or androgen receptor (AR) antagonists. In order to investigate the mechanisms relevant to the development of therapy-resistant tumours, we established a new tumour model which closely resembles the situation in patients who receive androgen ablation therapy. Androgen-sensitive LNCaP cells were kept in androgen-depleted medium for 87 passages. The new LNCaP cell subline established in this manner, LNCaP-abl, displayed a hypersensitive biphasic proliferative response to androgen until passage 75. Maximal proliferation of LNCaP-abl cells was achieved at 0.001 nM of the synthetic androgen methyltrienolone (R1881), whereas 0.01 nM of this compound induced the same effect in parental cells. At later passages (> 75), androgen exerted an inhibitory effect on growth of LNCaP-abl cells. The non-steroidal anti-androgen bicalutamide stimulated proliferation of LNCaP-abl cells. AR protein expression in LNCaP-abl cells increased approximately fourfold. The basal AR transcriptional activity was 30-fold higher in LNCaP-abl than in LNCaP cells. R1881 stimulated reporter gene activity in LNCaP-abl cells even at 0.01 nM, whereas 0.1 nM of R1881 was needed for induction of the same level of reporter gene activity in LNCaP cells. Bicalutamide that acts as a pure antagonist in parental LNCaP cells showed agonistic effects on AR transactivation activity in LNCaP-abl cells and was not able to block the effects of androgen in these cells. The non-steroidal AR blocker hydroxyflutamide exerted stimulatory effects on AR activity in both LNCaP and LNCaP-abl cells; however, the induction of reporter gene activity by hydroxyflutamide was 2.4- to 4-fold higher in the LNCaP-abl subline. The changes in AR activity were associated neither with a new alteration in AR cDNA sequence nor with amplification of the AR gene. Growth of LNCaP-abl xenografts in nude mice was stimulated by bicalutamide and repressed by testosterone. In conclusion, our results show for the first time that the non-steroidal anti-androgen bicalutamide acquires agonistic properties during long-term androgen ablation. These findings may have repercussions on the natural course of prostate cancer with androgen deprivation and on strategies of therapeutic intervention. © 1999 Cancer Research Campaig

Design, Validation and Annotation of Transcriptome-Wide Oligonucleotide Probes for the Oligochaete Annelid Eisenia fetida

High density oligonucleotide probe arrays have increasingly become an important tool in genomics studies. In organisms with incomplete genome sequence, one strategy for oligo probe design is to reduce the number of unique probes that target every non-redundant transcript through bioinformatic analysis and experimental testing. Here we adopted this strategy in making oligo probes for the earthworm Eisenia fetida, a species for which we have sequenced transcriptome-scale expressed sequence tags (ESTs). Our objectives were to identify unique transcripts as targets, to select an optimal and non-redundant oligo probe for each of these target ESTs, and to annotate the selected target sequences. We developed a streamlined and easy-to-follow approach to the design, validation and annotation of species-specific array probes. Four 244K-formatted oligo arrays were designed using eArray and were hybridized to a pooled E. fetida cRNA sample. We identified 63,541 probes with unsaturated signal intensities consistently above the background level. Target transcripts of these probes were annotated using several sequence alignment algorithms. Significant hits were obtained for 37,439 (59%) probed targets. We validated and made publicly available 63.5K oligo probes so the earthworm research community can use them to pursue ecological, toxicological, and other functional genomics questions. Our approach is efficient, cost-effective and robust because it (1) does not require a major genomics core facility; (2) allows new probes to be easily added and old probes modified or eliminated when new sequence information becomes available, (3) is not bioinformatics-intensive upfront but does provide opportunities for more in-depth annotation of biological functions for target genes; and (4) if desired, EST orthologs to the UniGene clusters of a reference genome can be identified and selected in order to improve the target gene specificity of designed probes. This approach is particularly applicable to organisms with a wealth of EST sequences but unfinished genome