Reconstructing an ultrametric galled phylogenetic network from a distance matrix

Given a distance matrix M that specifies the pairwise evolutionary distances between n species, the phylogenetic tree reconstruction problem asks for an edge-weighted phylogenetic tree that satisfies M, if one exists. We study some extensions of this problem to rooted phylogenetic networks. Our main result is an O(n(2) log n)-time algorithm for determining whether there is an ultrametric galled network that satisfies M, and if so, constructing one. In fact, if such an ultrametric galled network exists, our algorithm is guaranteed to construct one containing the minimum possible number of nodes with more than one parent (hybrid nodes). We also prove that finding a largest possible submatrix M' of M such that there exists an ultrametric galled network that satisfies M' is NP-hard. Furthermore, we show that given an incomplete distance matrix (i.e. where some matrix entries are missing), it is also NP-hard to determine whether there exists an ultrametric galled network which satisfies it.     

Dissection of a metastatic gene expression signature into distinct components

Background

Metastasis, the process whereby cancer cells spread, is in part caused by an incompletely understood interplay between cancer cells and the surrounding stroma. Gene expression studies typically analyze samples containing tumor cells and stroma. Samples with less than 50% tumor cells are generally excluded, thereby reducing the number of patients that can benefit from clinically relevant signatures.

Results

For a head-neck squamous cell carcinoma (HNSCC) primary tumor expression signature that predicts the presence of lymph node metastasis, we first show that reduced proportions of tumor cells results in decreased predictive accuracy. To determine the influence of stroma on the predictive signature and to investigate the interaction between tumor cells and the surrounding microenvironment, we used laser capture microdissection to divide the metastatic signature into six distinct components based on tumor versus stroma expression and on association with the metastatic phenotype. A strikingly skewed distribution of metastasis associated genes is revealed.

Conclusion

Dissection of predictive signatures into different components has implications for design of expression signatures and for our understanding of the metastatic process. Compared to primary tumors that have not formed metastases, primary HNSCC tumors that have metastasized are characterized by predominant down-regulation of tumor cell specific genes and exclusive up-regulation of stromal cell specific genes. The skewed distribution agrees with poor signature performance on samples that contain less than 50% tumor cells. Methods for reducing tumor composition bias that lead to greater predictive accuracy and an increase in the types of samples that can be included are presented.     

Inhibition of Cardiomyocytes Differentiation of Mouse Embryonic Stem Cells by CD38/cADPR/Ca2+ Signaling Pathway

Cyclic adenosine diphosphoribose (cADPR) is an endogenous Ca²⁺ mobilizing messenger that is formed by ADP-ribosyl cyclases from nicotinamide adenine dinucleotide (NAD). The main ADP-ribosyl cyclase in mammals is CD38, a multi-functional enzyme and a type II membrane protein. Here we explored the role of CD38-cADPR-Ca²⁺ in the cardiomyogenesis of mouse embryonic stem (ES) cells. We found that the mouse ES cells are responsive to cADPR and possess the key components of the cADPR signaling pathway. In vitro cardiomyocyte (CM) differentiation of mouse ES cells was initiated by embryoid body (EB) formation. Interestingly, beating cells appeared earlier and were more abundant in CD38 knockdown EBs than in control EBs. Real-time RT-PCR and Western blot analyses further showed that the expression of several cardiac markers, including GATA4, MEF2C, NKX2.5, and α-MLC, were increased markedly in CD38 knockdown EBs than those in control EBs. Similarly, FACS analysis showed that more cardiac Troponin T-positive CMs existed in CD38 knockdown or 8-Br-cADPR, a cADPR antagonist, treated EBs compared with that in control EBs. On the other hand, overexpression of CD38 in mouse ES cells significantly inhibited CM differentiation. Moreover, CD38 knockdown ES cell-derived CMs possess the functional properties characteristic of normal ES cell-derived CMs. Last, we showed that the CD38-cADPR pathway negatively modulated the FGF4-Erks1/2 cascade during CM differentiation of ES cells, and transiently inhibition of Erk1/2 blocked the enhanced effects of CD38 knockdown on the differentiation of CM from ES cells. Taken together, our data indicate that the CD38-cADPR-Ca²⁺ signaling pathway antagonizes the CM differentiation of mouse ES cells.     

Molecular docking and enzyme kinetic studies of dihydrotanshinone on metabolism of a model CYP2D6 probe substrate in human liver microsomes

The effects of Danshen and its active components (tanshinone I, tanshinone IIA, dihydrotanshinone and cryptotanshinone) on CYP2D6 activity was investigated by measuring the metabolism of a model CYP2D6 probe substrate, dextromethorphan to dextrorphan in human pooled liver microsomes. The ethanolic extract of crude Danshen (6.25-100 μg/ml) decreased dextromethorphan O-demethylation in vitro (IC(50)=23.3 μg/ml) and the water extract of crude Danshen (0.0625-1 mg/ml) showed no inhibition. A commercially available Danshen pill (31.25-500 μg/ml) also decreased CYP2D6 activity (IC(50)=265.8 μg/ml). Among the tanshinones, only dihydrotanshinone significantly inhibited CYP2D6 activity (IC(50)=35.4 μM), compared to quinidine, a specific CYP2D6 inhibitor (IC(50)=0.9 μM). Crytotanshinone, tanshinone I and tanshinone IIA produced weak inhibition, with IC(20) of 40.8 μM, 16.5 μM and 61.4 μM, respectively. Water soluble components such as salvianolic acid B and danshensu did not affect CYP2D6-mediated metabolism. Enzyme kinetics studies showed that inhibition of CYP2D6 activity by the ethanolic extract of crude Danshen and dihydrotanshinone was concentration-dependent, with K(i) values of 4.23 μg/ml and 2.53 μM, respectively, compared to quinidine, K(i)=0.41 μM. Molecular docking study confirmed that dihydrotanshinone and tanshinone I interacted with the Phe120 amino acid residue in the active cavity of CYP2D6 through Pi-Pi interaction, but did not interact with Glu216 and Asp301, the key residues for substrate binding. The logarithm of free binding energy of dihydrotanshinone (-7.6 kcal/mol) to Phe120 was comparable to quinidine (-7.0 kcal/mol) but greater than tanshinone I (-5.4 kcal/mol), indicating dihydrotanshinone has similar affinity to quinidine in binding to the catalytic site on CYP2D6.     

CFTR mediates bicarbonate-dependent activation of miR-125b in preimplantation embryo development

Although HCO₃⁻ is known to be required for early embryo development, its exact role remains elusive. Here we report that HCO₃⁻ acts as an environmental cue in regulating miR-125b expression through CFTR-mediated influx during preimplantation embryo development. The results show that the effect of HCO₃⁻ on preimplantation embryo development can be suppressed by interfering the function of a HCO₃⁻-conducting channel, CFTR, by a specific inhibitor or gene knockout. Removal of extracellular HCO₃⁻ or inhibition of CFTR reduces miR-125b expression in 2 cell-stage mouse embryos. Knockdown of miR-125b mimics the effect of HCO₃⁻ removal and CFTR inhibition, while injection of miR-125b precursor reverses it. Downregulation of miR-125b upregulates p53 cascade in both human and mouse embryos. The activation of miR-125b is shown to be mediated by sAC/PKA-dependent nuclear shuttling of NF-κB. These results have revealed a critical role of CFTR in signal transduction linking the environmental HCO₃⁻ to activation of miR-125b during preimplantation embryo development and indicated the importance of ion channels in regulation of miRNAs.     

H5-type influenza virus hemagglutinin is functionally recognized by the natural killer-activating receptor NKp44

Antiviral immune defenses involve natural killer (NK) cells. We previously showed that the NK-activating receptor NKp44 is involved in the functional recognition of H1-type influenza virus strains by NK cells. In the present study, we investigated the interaction of NKp44 and the hemagglutinin of a primary influenza virus H5N1 isolate. Here we show that recombinant NKp44 recognizes H5-expressing cells and specifically interacts with soluble H5 hemagglutinin. H5-pseudotyped lentiviral particles bind to NK cells expressing NKp44. Following interaction with target cells expressing H5, pseudotyped lentiviral particles, or membrane-associated H5, NK cells show NKp44-mediated induced activity. These findings indicate that NKp44-H5 interactions induce functional NK activation.     

Dynamics of a Cytokine Storm

Six volunteers experienced severe inflammatory response during the Phase I clinical trial of a monoclonal antibody that was designed to stimulate a regulatory T cell response. Soon after the trial began, each volunteer experienced a “cytokine storm”, a dramatic increase in cytokine concentrations. The monoclonal antibody, TGN1412, raised serum concentrations of both pro- and anti-inflammatory cytokines το very hiγh values during the first day, while lymphocyte and monocyte concentrations plummeted. Because the subjects were healthy and had no prior indications of immune deficiency, this event provided an unusual opportunity to study the dynamic interactions of cytokines and other measured parameters. Here, the response histories of nine cytokines have been modeled by a set of linear ordinary differential equations. A general search procedure identifies parameters of the model, whose response fits the data well during the five-day measurement period. The eighteenth-order model reveals plausible cause-and-effect relationships among the cytokines, showing how each cytokine induces or inhibits other cytokines. It suggests that perturbations in IL2, IL8, and IL10 have the most significant inductive effect, while IFN-γ and IL12 have the greatest inhibiting effect on other cytokine concentrations. Although TNF-α is a major pro-inflammatory factor, IFN-γ and three other cytokines have faster initial and median response to TGN1412 infusion. Principal-component analysis of the data reveals three clusters of similar cytokine responses: [TNF-α, IL1, IL10], [IFN-γ, IL2, IL4, IL8, and IL12], and [IL6]. IL1, IL6, IL10, and TNF-α have the highest degree of variability in response to uncertain initial conditions, exogenous effects, and parameter estimates. This study illuminates details of a cytokine storm event, and it demonstrates the value of linear modeling for interpreting complex, coupled biological system dynamics from empirical data.     

Genome Sequence of Pseudomonas mendocina DLHK,Isolated from a Biotrickling Reactor

Pseudomonas mendocina DLHK is an aerobic bacterium isolated from a biotrickling reactor which can remove nitric oxide, a common air pollutant from combustion exhaust gas. Here, we present the draft genome of Pseudomonas mendocina DLHK.