Influence of base stacking and hydrogen bonding on the fluorescence of 2-aminopurine and pyrrolocytosine in nucleic acids

Fluorescent nucleobase analogues are used extensively to probe the structure and dynamics of nucleic acids. The fluorescence of the adenine analogue 2-aminopurine and the cytosine analogue pyrrolocytosine is significantly quenched when the bases are located in regions of double-stranded nucleic acids. To allow more detailed structural information to be obtained from fluorescence studies using these bases, we have studied the excited-state properties of the bases at the CIS and TDB3LYP level in hydrogen-bonded and base-stacked complexes. The results reveal that the first excited state (the fluorescent state) of a hydrogen-bonded complex containing 2-aminopurine and thymine is just the first excited state of 2-aminopurine alone. However, the same cannot be said for structures in which 2-aminopurine is base stacked with other nucleobases. Stacking causes the molecular orbitals involved in the fluorescence transition to spread over more than one base. The predicted rate for the fluorescence transition is reduced, thus reducing the fluorescence quantum yield. The decrease in radiative rate varies with the stacking arrangement (e.g., A- or B-form DNA) and with the identity of the nucleobase with which 2-aminopurine is stacked. Stacking 2-aminopurine between two guanine moieties is shown to significantly decrease the energy gap between the first and second excited states. We do not find reliable evidence for a low-energy charge-transfer state in any of the systems that were studied. In the case of pyrrolocytosine, base stacking was found to reduce the oscillator strength for the fluorescence transition, but very little spreading of molecular orbitals across more than one base was observed.     

Quantitative trait loci for locomotor behavior in Drosophila melanogaster

Locomotion is an integral component of most animal behaviors and many human diseases and disorders are associated with locomotor deficits, but little is known about the genetic basis of natural variation in locomotor behavior. Locomotion is a complex trait, with variation attributable to the joint segregation of multiple interacting quantitative trait loci (QTL), with effects that are sensitive to the environment. We assessed variation in a component of locomotor behavior (locomotor reactivity) in a population of 98 recombinant inbred lines of Drosophila melanogaster and mapped four QTL affecting locomotor reactivity by linkage to polymorphic roo transposable element insertion sites. We used complementation tests of deficiencies to fine map these QTL to 12 chromosomal regions and complementation tests of mutations to identify 13 positional candidate genes affecting locomotor reactivity, including Dopa decarboxylase (Ddc), which catalyzes the final step in the synthesis of serotonin and dopamine. Linkage disequilibrium mapping in a population of 164 second chromosome substitution lines derived from a single natural population showed that polymorphisms at Ddc were associated with naturally occurring genetic variation in locomotor behavior. These data implicate variation in the synthesis of bioamines as a factor contributing to natural variation in locomotor reactivity.     

Multiple sequence alignment accuracy and phylogenetic inference

Phylogenies are often thought to be more dependent upon the specifics of the sequence alignment rather than on the method of reconstruction. Simulation of sequences containing insertion and deletion events was performed in order to determine the role that alignment accuracy plays during phylogenetic inference. Data sets were simulated for pectinate, balanced, and random tree shapes under different conditions (ultrametric equal branch length, ultrametric random branch length, nonultrametric random branch length). Comparisons between hypothesized alignments and true alignments enabled determination of two measures of alignment accuracy, that of the total data set and that of individual branches. In general, our results indicate that as alignment error increases, topological accuracy decreases. This trend was much more pronounced for data sets derived from more pectinate topologies. In contrast, for balanced, ultrametric, equal branch length tree shapes, alignment inaccuracy had little average effect on tree reconstruction. These conclusions are based on average trends of many analyses under different conditions, and any one specific analysis, independent of the alignment accuracy, may recover very accurate or inaccurate topologies. Maximum likelihood and Bayesian, in general, outperformed neighbor joining and maximum parsimony in terms of tree reconstruction accuracy. Results also indicated that as the length of the branch and of the neighboring branches increase, alignment accuracy decreases, and the length of the neighboring branches is the major factor in topological accuracy. Thus, multiple-sequence alignment can be an important factor in downstream effects on topological reconstruction.     

Development of a screening assay for surrogate markers of CHK1 inhibitor-induced cell cycle release

Chk1 is a key regulator of the S and G2/M checkpoints and is activated following DNA damage by agents such as the topoisomerase I inhibitor camptothecin (CPT). It has been proposed that Chk1 inhibitors used in combination with such a DNA damaging agent to treat tumors would potentiate cytotoxicity and increase the therapeutic index, particularly in tumors lacking functional p53. The aim of this study was to determine whether gene expression analysis could be used to inform lead optimization of a novel series of Chk1 inhibitors. The candidate small-molecule Chk1 inhibitors were used in combination with CPT to identify potential markers of functional Chk1 inhibition, as well as resulting cell cycle progression, using cDNA-based microarrays. Differential expression of several of these putative marker genes was further validated by RT-PCR for use as a medium-throughput assay. In the presence of DNA damage, Chk1 inhibitors altered CPT-dependent effects on the expression of cell cycle and DNA repair genes in a manner consistent with a Chk1-specific mechanism of action. Furthermore, differential expression of selected marker genes, cyclin E2, EGR1, and DDIT3, was dose dependent for Chk1 inhibition. RT-PCR results for these genes following treatment with a panel of Chk1 inhibitors showed a strong correlation between marker gene response and the ability of each compound to abrogate cell cycle arrest in situ following CPT-induced DNA damage. These results demonstrate the utility of global expression analysis to identify surrogate markers, providing an alternative method for rapid compound characterization to support advancement decisions in early drug discovery.     

The effects of chronic immune stimulation on muscle growth in rainbow trout

Successful production of aquaculture species depends on efficient growth with low susceptibility to disease. Therefore, selection programs have focused on rapid growth combined with disease resistance. However, chronic immune stimulation diminishes muscle growth (a syndrome referred to as cachexia), and decreases growth efficiency in production animals, including rainbow trout. In mammals, recent results show that increased levels of pro-inflammatory cytokines, such as those seen during an immune assault, specifically target myosin and MyoD and inhibit muscle growth. This suggests that increased disease resistance in fish, a desired trait for production, may actually decrease the growth of muscle, the main aquacultural commodity. To test this possibility, a rainbow trout model of cachexia was developed and characterized. A six-week study was conducted in which rainbow trout were chronically immune stimulated by repeated injections of LPS. Growth indices were monitored, and whole body and muscle proximate analyses, real-time PCR, and Western blotting were conducted to examine the resulting cachectic phenotype. Muscle ratio was decreased in fish chronically immunostimulated, however expression levels of MyoD2 and myosin were not decreased compared to fish that were not immunostimulated, indicating that while muscle accretion was altered, the mechanism by which it occurred was somewhat different than that characterized in mammals. Microarray analysis was used to compare gene expression in fish that had been chronically immunostimulated versus those that had not to identify possible alternative mechanisms of cachexia in fish.     

MDC1 maintains genomic stability by participating in the amplification of ATM-dependent DNA damage signals

MDC1 functions in checkpoint activation and DNA repair following DNA damage. To address the physiological role of MDC1, we disrupted the MDC1 gene in mice. MDC1-/- mice recapitulated many phenotypes of H2AX-/- mice, including growth retardation, male infertility, immune defects, chromosome instability, DNA repair defects, and radiation sensitivity. At the molecular level, H2AX, MDC1, and ATM form a positive feedback loop, with MDC1 directly mediating the interaction between H2AX and ATM. MDC1 binds phosphorylated H2AX through its BRCT domain and ATM through its FHA domain. Through these interactions, MDC1 accumulates activated ATM flanking the sites of DNA damage, facilitating further ATM-dependent phosphorylation of H2AX and the amplification of DNA damage signals. In the absence of MDC1, many downstream ATM signaling events are defective. These results suggest that MDC1, as a signal amplifier of the ATM pathway, is vital in controlling proper DNA damage response and maintaining genomic stability.     

Microbial community dynamics in a humic lake: differential persistence of common freshwater phylotypes

In an effort to better understand the factors contributing to patterns in freshwater bacterioplankton community composition and diversity, we coupled automated ribosomal intergenic spacer analysis (ARISA) to analysis of 16S ribosomal RNA (rRNA) gene sequences to follow the persistence patterns of 46 individual phylotypes over 3 years in Crystal Bog Lake. Additionally, we sought to identify linkages between the observed phylotype variations and known chemical and biological drivers. Sequencing of 16S rRNA genes obtained from the water column indicated the presence of phylotypes associated with the Actinobacteria, Bacteroidetes, Firmicutes, Proteobacteria, TM7 and Verrucomicrobia phyla, as well as phylotypes with unknown affiliation. Employment of the 16S rRNA gene/ARISA method revealed that specific phylotypes varied independently of the entire bacterial community dynamics. Actinobacteria, which were present on greater than 95% of sampling dates, did not share the large temporal variability of the other identified phyla. Examination of phylotype relative abundance patterns (inferred using ARISA fragment relative fluorescence) revealed a strong correlation between the dominant phytoplankton succession and the relative abundance patterns of the majority of individual phylotypes. Further analysis revealed covariation among unique phylotypes, which formed several distinct bacterial assemblages correlated with particular phytoplankton communities. These data indicate the existence of unique persistence patterns for different common freshwater phylotypes, which may be linked to the presence of dominant phytoplankton species.     

The statistics of identifying differentially expressed genes in Expresso and TM4: a comparison

Background  

Analysis of DNA microarray data takes as input spot intensity measurements from scanner software and returns differential expression of genes between two conditions, together with a statistical significance assessment. This process typically consists of two steps: data normalization and identification of differentially expressed genes through statistical analysis. The Expresso microarray experiment management system implements these steps with a two-stage, log-linear ANOVA mixed model technique, tailored to individual experimental designs. The complement of tools in TM4, on the other hand, is based on a number of preset design choices that limit its flexibility. In the TM4 microarray analysis suite, normalization, filter, and analysis methods form an analysis pipeline. TM4 computes integrated intensity values (IIV) from the average intensities and spot pixel counts returned by the scanner software as input to its normalization steps. By contrast, Expresso can use either IIV data or median intensity values (MIV). Here, we compare Expresso and TM4 analysis of two experiments and assess the results against qRT-PCR data.     

Quality control for terms and definitions in ontologies and taxonomies

Background  

Ontologies and taxonomies are among the most important computational resources for molecular biology and bioinformatics. A series of recent papers has shown that the Gene Ontology (GO), the most prominent taxonomic resource in these fields, is marked by flaws of certain characteristic types, which flow from a failure to address basic ontological principles. As yet, no methods have been proposed which would allow ontology curators to pinpoint flawed terms or definitions in ontologies in a systematic way.