Estimating divergence times in large phylogenetic trees

A new method, PATHd8, for estimating ultrametric trees from trees with edge (branch) lengths proportional to the number of substitutions is proposed. The method allows for an arbitrary number of reference nodes for time calibration, each defined either as absolute age, minimum age, or maximum age, and the tree need not be fully resolved. The method is based on estimating node ages by mean path lengths from the node to the leaves but correcting for deviations from a molecular clock suggested by reference nodes. As opposed to most existing methods allowing substitution rate variation, the new method smoothes substitution rates locally, rather than simultaneously over the whole tree, thus allowing for analysis of very large trees. The performance of PATHd8 is compared with other frequently used methods for estimating divergence times. In analyses of three separate data sets, PATHd8 gives similar divergence times to other methods, the largest difference being between crown group ages, where unconstrained nodes get younger ages when analyzed with PATHd8. Overall, chronograms obtained from other methods appear smoother, whereas PATHd8 preserves more of the heterogeneity seen in the original edge lengths. Divergence times are most evenly spread over the chronograms obtained from the Bayesian implementation and the clock-based Langley-Fitch method, and these two methods produce very similar ages for most nodes. Evaluations of PATHd8 using simulated data suggest that PATHd8 is slightly less precise compared with penalized likelihood, but it gives more sensible answers for extreme data sets. A clear advantage with PATHd8 is that it is more or less instantaneous even with trees having several thousand leaves, whereas other programs often run into problems when analyzing trees with hundreds of leaves. PATHd8 is implemented in freely available software.     

Repeated inhalations of diesel exhaust particles and oxidatively damaged DNA in young oxoguanine DNA glycosylase (OGG1) deficient mice

DNA repair may prevent increased levels of oxidatively damaged DNA from prolonged oxidative stress induced by, e.g. exposure to diesel exhaust particles (DEP). We studied oxidative damage to DNA in broncho-alveolar lavage cells, lungs, and liver after 4 × 1.5 h inhalations of DEP (20 mg/m³) in Ogg1^- / -  and wild type (WT) mice with similar extent of inflammation. DEP exposure increased lung levels of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) in Ogg1^- / -  mice, whereas no effect on 8-oxodG or oxidized purines in terms of formamidopyrimidine DNA glycosylase (FPG) sites was observed in WT mice. In both unexposed and exposed Ogg1^- / -  mice the level of FPG sites in the lungs was 3-fold higher than in WT mice. The high basal level of FPG sites in Ogg1^- / -  mice probably saturated the assay and prevented detection of DEP-generated damage. In conclusion, Ogg1^- / -  mice have elevated pulmonary levels of FPG sites and accumulate genomic 8-oxodG after repeated inhalations of DEP.     

A conserved modified wobble nucleoside (mcm5s2U) in lysyl-tRNA is required for viability in yeast

Transfer RNAs specific for Gln, Lys, and Glu from all organisms (except Mycoplasma) and organelles have a 2-thiouridine derivative (xm(5)s(2)U) as wobble nucleoside. These tRNAs read the A- and G-ending codons in the split codon boxes His/Gln, Asn/Lys, and Asp/Glu. In eukaryotic cytoplasmic tRNAs the conserved constituent (xm(5)-) in position 5 of uridine is 5-methoxycarbonylmethyl (mcm(5)). A protein (Tuc1p) from yeast resembling the bacterial protein TtcA, which is required for the synthesis of 2-thiocytidine in position 32 of the tRNA, was shown instead to be required for the synthesis of 2-thiouridine in the wobble position (position 34). Apparently, an ancient member of the TtcA family has evolved to thiolate U34 in tRNAs of organisms from the domains Eukarya and Archaea. Deletion of the TUC1 gene together with a deletion of the ELP3 gene, which results in the lack of the mcm(5) side chain, removes all modifications from the wobble uridine derivatives of the cytoplasmic tRNAs specific for Gln, Lys, and Glu, and is lethal to the cell. Since excess of the unmodified form of these three tRNAs rescued the double mutant elp3 tuc1, the primary function of mcm(5)s(2)U34 seems to be to improve the efficiency to read the cognate codons rather than to prevent mis-sense errors. Surprisingly, overexpression of the mcm(5)s(2)U-lacking tRNA(Lys) alone was sufficient to restore viability of the double mutant.     

Structure probing of tmRNA in distinct stages of trans-translation

Ribosomes stalled on problematic mRNAs in bacterial cells can be rescued by transfer-messenger RNA (tmRNA), its helper protein (small protein B, SmpB), and elongation factor Tu (EF-Tu) through a mechanism called trans-translation. In this work we used lead(II) footprinting to probe the interactions of tmRNA with SmpB and other components of the translation machinery at different steps of the trans-translation cycle. Ribosomes with a short nascent peptide stalled on a truncated mRNA were reacted with Ala-tmRNA*EF-Tu*GTP, SmpB, and other translation components to initiate and execute trans-translation. Free tmRNA was probed with lead(II) acetate with and without SmpB, and ribosome bound tmRNA was probed in one of four different trans-translation states stabilized by antibiotic addition or selective exclusion of translation components. For comparison, we also analyzed lead(II) cleavage patterns of tmRNA in vivo in a wild-type as well as in an SmpB-deficient Escherichia coli strain. We observed some specific cleavages/protections in tmRNA for the individual steps of trans-translation, but the overall tmRNA conformation appeared to be similar in the stages analyzed. Our findings suggest that, in vivo, a dominant fraction of tmRNA is in complex with SmpB and that, in vitro, SmpB remains tmRNA bound at the initial steps of trans-translation.     

Predicting survival from microarray data--a comparative study

MOTIVATION: Survival prediction from gene expression data and other high-dimensional genomic data has been subject to much research during the last years. These kinds of data are associated with the methodological problem of having many more gene expression values than individuals. In addition, the responses are censored survival times. Most of the proposed methods handle this by using Cox's proportional hazards model and obtain parameter estimates by some dimension reduction or parameter shrinkage estimation technique. Using three well-known microarray gene expression data sets, we compare the prediction performance of seven such methods: univariate selection, forward stepwise selection, principal components regression (PCR), supervised principal components regression, partial least squares regression (PLS), ridge regression and the lasso. RESULTS: Statistical learning from subsets should be repeated several times in order to get a fair comparison between methods. Methods using coefficient shrinkage or linear combinations of the gene expression values have much better performance than the simple variable selection methods. For our data sets, ridge regression has the overall best performance. AVAILABILITY: Matlab and R code for the prediction methods are available at http://www.med.uio.no/imb/stat/bmms/software/microsurv/.     

Antifungal compounds from cultures of dairy propionibacteria type strains

Antifungal compounds from cultures of five type strains of dairy propionibacteria, as well as from the cultivation medium, were studied. Cell-free supernatants and medium were fractionated by C(18) solid phase extraction. The aqueous 95% acetonitrile fractions were analyzed by GC-MS or subjected to reversed-phase HPLC, to identify, quantify or isolate antifungal substances. The resulting HPLC fractions were screened for antifungal activity against the mold Aspergillus fumigatus and the yeast Rhodotorula mucilaginosa. Active fractions were further separated by HPLC and the structures of the compounds were determined by spectroscopic and chromatographic methods. All five strains produced 3-phenyllactic acid, at concentrations ranging from 1.0 microg mL(-1) (Propionibacterium freudenreichii ssp. shermanii) to 15.1 microg mL(-1) (Propionibacterium thoenii), and at L/D -ratios ranging from 2 : 3 (Propionibacterium acidipropionici) to 9 : 1 (Propionibacterium freudenreichii). A number of active compounds found in cultures of propionibacteria were also present in noninoculated growth medium: two antifungal diketopiperazines, cyclo(L-Phe-L-Pro) and cyclo(L-Ile-L-Pro), and seven antifungal linear peptides. Three of the linear peptides corresponded to sequences found in the medium component casein, suggesting their origin from this component, whereas the diketopiperazines were suggested to be formed from medium peptides by heat treatment.     

Fouling-resistant surfaces of tropical sea stars

Qualitative evidence suggests sea stars are free of fouling organisms; however the presence of fouling-resistant surfaces of sea stars has not previously been documented. Field surveys were conducted in northern Queensland, Australia, during the wet and dry seasons and several tropical sea star species were examined for surface-associated micro- and macro-organisms. Mean bacterial abundances on seven sea star species were approximately 10(4) to 10(5) cells cm(-2) during both seasons. There were no consistent trends in bacterial abundances with season, species and aboral positions on sea star arms. No common generalist fouling organisms, such as algae, barnacles, serpulid polychaetes, bryozoans and ascidians, were found on any specimens of 12 sea star species. However, low numbers of parasitic and commensal macro-organisms were found on six sea star species. The gastropods Parvioris fulvescens, Asterolamia hians, Thyca (Granulithyca) nardoafrianti and Thyca crystallina were found exclusively on the sea stars Archaster typicus, Astropecten indicus, Nardoa pauciforis and Linckia laevigata, respectively. The shrimp Periclimenes soror was only found on Acanthaster planci, and the polychaete Ophiodromus sp. on A. typicus. The copepods Stellicola illgi and Paramolgus sp. were only found on L. laevigata and Echinaster luzonicus, respectively. As no common generalist fouling organisms were discovered, sea stars offer an excellent model to investigate the mechanisms driving fouling-resistant surfaces and the selective settlement of specialist invertebrates.     

EphA4-dependent axon guidance is mediated by the RacGAP alpha2-chimaerin

Neuronal network formation in the developing nervous system is dependent on the accurate navigation of nerve cell axons and dendrites, which is controlled by attractive and repulsive guidance cues. Ephrins and their cognate Eph receptors mediate many repulsive axonal guidance decisions by intercellular interactions resulting in growth cone collapse and axon retraction of the Eph-presenting neuron. We show that the Rac-specific GTPase-activating protein alpha2-chimaerin binds activated EphA4 and mediates EphA4-triggered axonal growth cone collapse. alpha-Chimaerin mutant mice display a phenotype similar to that of EphA4 mutant mice, including aberrant midline axon guidance and defective spinal cord central pattern generator activity. Our results reveal an alpha-chimaerin-dependent signaling pathway downstream of EphA4, which is essential for axon guidance decisions and neuronal circuit formation in vivo.     

A dynamic mass balance model for phosphorus fluxes and concentrations in coastal areas

This paper presents a general, process-based mass balance model (CoastMab) for total phosphorus (TP) in defined coastal areas (at the ecosystem scale). The model is based on ordinary differential equations and calculates inflow, outflow and internal fluxes on a monthly basis. It consists of four compartments: surface water, deep water, erosion/transportation areas for fine sediments and accumulation areas for fine sediments. The separation between surface water and deep water is not done based on water temperature, but on sedimentological criteria instead (from the theoretical wave base). There are algorithms for all major internal TP fluxes (sedimentation, resuspension, diffusion, mixing and burial). Validations were performed using data from 21 different Baltic coastal areas. The results show that the model predicts monthly TP in water and chlorophyll a very well (generally within the uncertainty bands of the empirical data). The model has also been put through sensitivity tests, which show that the most important factor regulating the predictions of the model is generally the TP concentration in the sea beyond the coast. The model is simple to apply, since all driving variables may be accessed from maps or monitoring programs. The driving variables include coastal area, section area (between the defined coastal area and the adjacent sea), mean and maximum depths, latitude (used to predict water temperatures, stratification and mixing), salinity and TP concentration in the sea. Many of the model structures are general and could be used for areas other than those included in this study, e.g., for open coasts, estuaries or tidal coasts, as well as for other substances than phosphorus.     

Androgen dependent regulation of protein kinase A subunits in prostate cancer cells

Neuroendocrine (NE) cells may play a role in prostate cancer progression. Both androgen deprivation and cAMP are well known inducers of NE differentiation (NED) in the prostate. Gene-expression profiling of LNCaP cells, incubated in androgen stripped medium, showed that the Cbeta isoform of PKA is up-regulated during NE differentiation. Furthermore, by using semi-quantitative RT-PCR and immunoblotting analysis, we observed that the Cbeta splice variants are differentially regulated during this process. Whereas the Cbeta2 splice variant is down-regulated in growth arrested LNCaP cells, the Cbeta1, Cbeta3 and Cbeta4 variants, as well as the RIIbeta subunit of PKA, are induced in NE-like LNCaP cells. The opposite effect of Cbeta expression could be mimicked by androgen stimulation, implying the Cbeta gene of PKA as a putative new target gene for the androgen receptor in prostate cancer. Moreover, to investigate expression of PKA subunits during prostate cancer progression, we did immunoblotting of several prostatic cell lines and normal and tumor tissue from prostate cancer patients. Interestingly, multiple Cbeta subunits were also observed in human prostate specimens, and the Cbeta2 variant was up-regulated in tumor cells. In conclusion, it seems that the Cbeta isoforms play different roles in proliferation and differentiation and could therefore be potential markers for prostate cancer progression.