Reducing the worst case running times of a family of RNA and CFG problems, using Valiant's approach

Background

RNA secondary structure prediction is a mainstream bioinformatic domain, and is key to computational analysis of functional RNA. In more than 30 years, much research has been devoted to defining different variants of RNA structure prediction problems, and to developing techniques for improving prediction quality. Nevertheless, most of the algorithms in this field follow a similar dynamic programming approach as that presented by Nussinov and Jacobson in the late 70's, which typically yields cubic worst case running time algorithms. Recently, some algorithmic approaches were applied to improve the complexity of these algorithms, motivated by new discoveries in the RNA domain and by the need to efficiently analyze the increasing amount of accumulated genome-wide data.

Results

We study Valiant's classical algorithm for Context Free Grammar recognition in sub-cubic time, and extract features that are common to problems on which Valiant's approach can be applied. Based on this, we describe several problem templates, and formulate generic algorithms that use Valiant's technique and can be applied to all problems which abide by these templates, including many problems within the world of RNA Secondary Structures and Context Free Grammars.

Conclusions

The algorithms presented in this paper improve the theoretical asymptotic worst case running time bounds for a large family of important problems. It is also possible that the suggested techniques could be applied to yield a practical speedup for these problems. For some of the problems (such as computing the RNA partition function and base-pair binding probabilities), the presented techniques are the only ones which are currently known for reducing the asymptotic running time bounds of the standard algorithms.     

Monitoring influenza activity in the United States: a comparison of traditional surveillance systems with Google Flu Trends

Background

Google Flu Trends was developed to estimate US influenza-like illness (ILI) rates from internet searches; however ILI does not necessarily correlate with actual influenza virus infections.

Methods and Findings

Influenza activity data from 2003–04 through 2007–08 were obtained from three US surveillance systems: Google Flu Trends, CDC Outpatient ILI Surveillance Network (CDC ILI Surveillance), and US Influenza Virologic Surveillance System (CDC Virus Surveillance). Pearson''s correlation coefficients with 95% confidence intervals (95% CI) were calculated to compare surveillance data. An analysis was performed to investigate outlier observations and determine the extent to which they affected the correlations between surveillance data. Pearson''s correlation coefficient describing Google Flu Trends and CDC Virus Surveillance over the study period was 0.72 (95% CI: 0.64, 0.79). The correlation between CDC ILI Surveillance and CDC Virus Surveillance over the same period was 0.85 (95% CI: 0.81, 0.89). Most of the outlier observations in both comparisons were from the 2003–04 influenza season. Exclusion of the outlier observations did not substantially improve the correlation between Google Flu Trends and CDC Virus Surveillance (0.82; 95% CI: 0.76, 0.87) or CDC ILI Surveillance and CDC Virus Surveillance (0.86; 95%CI: 0.82, 0.90).

Conclusions

This analysis demonstrates that while Google Flu Trends is highly correlated with rates of ILI, it has a lower correlation with surveillance for laboratory-confirmed influenza. Most of the outlier observations occurred during the 2003–04 influenza season that was characterized by early and intense influenza activity, which potentially altered health care seeking behavior, physician testing practices, and internet search behavior.     

Human telomeres maintain their overhang length at senescence

Normal human cells in culture enter replicative senescence after a finite number of population doublings. The exact molecular mechanisms triggering the growth arrest are poorly understood. A recent report on the disappearance of the G-rich 3' telomeric overhang in senescent cells led to the hypothesis that loss of the 3' G-rich overhang is the molecular signal that triggers senescence. Here, we describe a quantitative assay to measure the length of the G-rich 3' telomeric overhangs from cultured cells. Using both this assay and the conventional nondenaturing hybridization assay for measuring G-rich overhangs, we show that normal human fibroblasts can maintain their overhangs at senescence. Furthermore, cells do not lose their overhangs when they bypass senescence after the inactivation of p53 and Rb. We thus conclude that a global reduction in overhang length is not the molecular signal that triggers replicative senescence.     

Neural circuitry that governs Drosophila male courtship behavior

Male-specific fruitless (fru) products (Fru(M)) are both necessary and sufficient to "hardwire" the potential for male courtship behavior into the Drosophila nervous system. Fru(M) is expressed in approximately 2% of neurons in the male nervous system, but not in the female. We have targeted the insertion of GAL4 into the fru locus, allowing us to visualize and manipulate the Fru(M)-expressing neurons in the male as well as their counterparts in the female. We present evidence that these neurons are directly and specifically involved in male courtship behavior and that at least some of them are interconnected in a circuit. This circuit includes olfactory neurons required for the behavioral response to sex pheromones. Anatomical differences in this circuit that might account for the dramatic differences in male and female sexual behavior are not apparent.     

Telomere-end processing the terminal nucleotides of human chromosomes

Mammalian telomeres end in single-stranded, G-rich 3' overhangs resulting from both the "end-replication problem" (the inability of DNA polymerase to replicate the very end of the telomeres) and postreplication processing. Telomeric G-rich overhangs are precisely defined in ciliates; the length and the terminal nucleotides are fixed. Human telomeres have very long overhangs that are heterogeneous in size (35-600 nt), indicating that their processing must differ in some respects from model organisms. We developed telomere-end ligation protocols that allowed us to identify the terminal nucleotides of both the C-rich and the G-rich telomere strands. Up to approximately 80% of the C-rich strands terminate in CCAATC-5', suggesting that after replication a nuclease with high specificity or constrained action acts on the C strand. In contrast, the G-terminal nucleotide was less precise than Tetrahymena and Euplotes but still had a bias that changed as a function of telomerase expression.     

Lesion bypass by human DNA polymerase mu reveals a template-dependent, sequence-independent nucleotidyl transferase activity

DNA polymerase mu (pol mu), which is related to terminal deoxynucleotidyl transferase and DNA polymerase beta, is thought to be involved in non-homologous end joining and V(D)J recombination. Pol mu is induced by ionizing radiation and exhibits low fidelity. Analysis of translesion replication by purified human pol mu revealed that it bypasses a synthetic abasic site with high efficiency, using primarily a misalignment mechanism. It can also replicate across two tandem abasic sites, using the same mechanism. Pol mu extends primers whose 3'-terminal nucleotides are located opposite the abasic site. Most remarkably, this extension occurs via a mode of nucleotidyl transferase activity, which does not depend on the sequence of the template. This is not due to simple terminal nucleotidyl transferase activity, because pol mu is unable to add dNTPs to an oligo(dT)29 primer or to a blunt end duplex oligonucleotide under standard conditions. Thus, pol mu is a dual mode DNA-synthesizing enzyme, which can act as either a classical DNA polymerase or as a non-canonical, template-dependent, but sequence-independent nucleotidyl transferase. To our knowledge, this is the first report on a DNA-synthesizing enzyme with such properties. These activities may be required for its function in non-homologous end joining in the processing of DNA ends prior to ligation.     

The matrilineal ancestry of Ashkenazi Jewry: portrait of a recent founder event

Both the extent and location of the maternal ancestral deme from which the Ashkenazi Jewry arose remain obscure. Here, using complete sequences of the maternally inherited mitochondrial DNA (mtDNA), we show that close to one-half of Ashkenazi Jews, estimated at 8,000,000 people, can be traced back to only 4 women carrying distinct mtDNAs that are virtually absent in other populations, with the important exception of low frequencies among non-Ashkenazi Jews. We conclude that four founding mtDNAs, likely of Near Eastern ancestry, underwent major expansion(s) in Europe within the past millennium.     

Near-infrared fluorescent deoxyglucose analogue for tumor optical imaging in cell culture and living mice

2-deoxy-2-[18F]fluoro-D-glucose ([18F]FDG) has extensively been used for clinical diagnosis, staging, and therapy monitoring of cancer and other diseases. Nonradioactive glucose analogues enabling the screening of the glucose metabolic rate of tumors are of particular interest for anticancer drug development. A nonradioactive fluorescent deoxyglucose analogue may have many applications for both imaging of tumors and monitoring therapeutic efficacy of drugs in living animals and may eventually translate to clinical applications. We found that a fluorescent 2-deoxyglucose analogue, 2-[N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino]-2-deoxy-D-glucose (2-NBDG), can be delivered in several tumor cells via the glucose transporters (GLUTs). We therefore conjugated D-glucosamine with a near-infrared (NIR) fluorphor Cy5.5 and tested the feasibility of the Cy5.5-D-glucosamine (Cy5.5-2DG) conjugate for NIR fluorescence imaging of tumors in a preclinical xenograft animal model. Cy5.5-2DG was prepared by conjugating Cy5.5 monofunctional N-hydroxysuccinimide ester (Cy5.5-NHS) and D-glucosamine followed by high-performance liquid chromatography purification. The accumulation of Cy5.5-2DG and Cy5.5-NHS in different tumor cell lines at 37 and 4 degrees C were imaged using a fluorescence microscope. Tumor targeting and retention of Cy5.5-2DG and Cy5.5-NHS in a subcutaneous U87MG glioma and A375M melanoma tumor model were evaluated and quantified by a Xenogen IVIS 200 optical cooled charged-coupled device system. Fluorescence microscopy imaging shows that Cy5.5-2DG and Cy5.5-NHS are taken up and trapped by a variety of tumor cell lines at 37 degrees C incubation, while they exhibit marginal uptake at 4 degrees C. The tumor cell uptake of Cy5.5-2DG cannot be blocked by the 50 mM D-glucose, suggesting that Cy5.5-2DG may not be delivered in tumor cells by GLUTs. U87MG and A375M tumor localization was clearly visualized in living mice with both NIR fluorescent probes. Tumor/muscle contrast was clearly visible as early as 30 min postinjection (pi), and the highest U87MG tumor/muscle ratios of 2.81 +/- 0.10 and 3.34 +/- 0.23 were achieved 24 h pi for Cy5.5-2DG and Cy5.5-NHS, respectively. While as a comparison, the micropositron emission tomography imaging study shows that [18F]FDG preferentially localizes to the U87MG tumor, with resulting tumor/muscle ratios ranging from 3.89 to 4.08 after 30 min to 2 h postadministration of the probe. In conclusion, the NIR fluorescent glucose analogues, Cy5.5-2DG and Cy5.5-NHS, both demonstrate tumor-targeting abilities in cell culture and living mice. More studies are warranted to further explore their application for optical tumor imaging. To develop NIR glucose analogues with the ability to target GLUTs/hexokinase, it is highly important to select NIR dyes with a reasonable molecular size.