Bayesian and Maximum Likelihood Phylogenetic Analyses of Protein Sequence Data Under Relative Branch-Length Differences and Model Violation

Background: Bayesian phylogenetic inference holds promise as an alternative to maximum likelihood, particularly for large molecular-sequence data sets. We have investigated the performance of Bayesian inference with empirical and simulated protein-sequence data under conditions of relative branch-length differences and model violation. Results: With empirical protein-sequence data, Bayesian posterior probabilities provide more-generous estimates of subtree reliability than does the nonparametric bootstrap combined with maximum likelihood inference, reaching 100% posterior probability at bootstrap proportions around 80%. With simulated 7-taxon protein-sequence datasets, Bayesian posterior probabilities are somewhat more generous than bootstrap proportions, but do not saturate. Compared with likelihood, Bayesian phylogenetic inference can be as or more robust to relative branch-length differences for datasets of this size, particularly when among-sites rate variation is modeled using a gamma distribution. When the (known) correct model was used to infer trees, Bayesian inference recovered the (known) correct tree in 100% of instances in which one or two branches were up to 20-fold longer than the others. At ratios more extreme than 20-fold, topological accuracy of reconstruction degraded only slowly when only one branch was of relatively greater length, but more rapidly when there were two such branches. Under an incorrect model of sequence change, inaccurate trees were sometimes observed at less extreme branch-length ratios, and (particularly for trees with single long branches) such trees tended to be more inaccurate. The effect of model violation on accuracy of reconstruction for trees with two long branches was more variable, but gamma-corrected Bayesian inference nonetheless yielded more-accurate trees than did either maximum likelihood or uncorrected Bayesian inference across the range of conditions we examined. Assuming an exponential Bayesian prior on branch lengths did not improve, and under certain extreme conditions significantly diminished, performance. The two topology-comparison metrics we employed, edit distance and Robinson-Foulds symmetric distance, yielded different but highly complementary measures of performance. Conclusions: Our results demonstrate that Bayesian inference can be relatively robust against biologically reasonable levels of relative branch-length differences and model violation, and thus may provide a promising alternative to maximum likelihood for inference of phylogenetic trees from protein-sequence data

Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation

BACKGROUND: Bayesian phylogenetic inference holds promise as an alternative to maximum likelihood, particularly for large molecular-sequence data sets. We have investigated the performance of Bayesian inference with empirical and simulated protein-sequence data under conditions of relative branch-length differences and model violation. RESULTS: With empirical protein-sequence data, Bayesian posterior probabilities provide more-generous estimates of subtree reliability than does the nonparametric bootstrap combined with maximum likelihood inference, reaching 100% posterior probability at bootstrap proportions around 80%. With simulated 7-taxon protein-sequence datasets, Bayesian posterior probabilities are somewhat more generous than bootstrap proportions, but do not saturate. Compared with likelihood, Bayesian phylogenetic inference can be as or more robust to relative branch-length differences for datasets of this size, particularly when among-sites rate variation is modeled using a gamma distribution. When the (known) correct model was used to infer trees, Bayesian inference recovered the (known) correct tree in 100% of instances in which one or two branches were up to 20-fold longer than the others. At ratios more extreme than 20-fold, topological accuracy of reconstruction degraded only slowly when only one branch was of relatively greater length, but more rapidly when there were two such branches. Under an incorrect model of sequence change, inaccurate trees were sometimes observed at less extreme branch-length ratios, and (particularly for trees with single long branches) such trees tended to be more inaccurate. The effect of model violation on accuracy of reconstruction for trees with two long branches was more variable, but gamma-corrected Bayesian inference nonetheless yielded more-accurate trees than did either maximum likelihood or uncorrected Bayesian inference across the range of conditions we examined. Assuming an exponential Bayesian prior on branch lengths did not improve, and under certain extreme conditions significantly diminished, performance. The two topology-comparison metrics we employed, edit distance and Robinson-Foulds symmetric distance, yielded different but highly complementary measures of performance. CONCLUSIONS: Our results demonstrate that Bayesian inference can be relatively robust against biologically reasonable levels of relative branch-length differences and model violation, and thus may provide a promising alternative to maximum likelihood for inference of phylogenetic trees from protein-sequence data

Anaplastic Lymphoma Kinase Spares Organ Growth during Nutrient Restriction in Drosophila

SummaryDeveloping animals survive periods of starvation by protecting the growth of critical organs at the expense of other tissues. Here, we use Drosophila to explore the as yet unknown mechanisms regulating this privileged tissue growth. As in mammals, we observe in Drosophila that the CNS is more highly spared than other tissues during nutrient restriction (NR). We demonstrate that anaplastic lymphoma kinase (Alk) efficiently protects neural progenitor (neuroblast) growth against reductions in amino acids and insulin-like peptides during NR via two mechanisms. First, Alk suppresses the growth requirement for amino acid sensing via Slimfast/Rheb/TOR complex 1. And second, Alk, rather than insulin-like receptor, primarily activates PI3-kinase. Alk maintains PI3-kinase signaling during NR as its ligand, Jelly belly (Jeb), is constitutively expressed from a glial cell niche surrounding neuroblasts. Together, these findings identify a brain-sparing mechanism that shares some regulatory features with the starvation-resistant growth programs of mammalian tumors.PaperCli

Community Impacts of Prosopis Juliflora Invasion: Biogeographic and Congeneric Comparisons

We coordinated biogeographical comparisons of the impacts of an exotic invasive tree in its native and non-native ranges with a congeneric comparison in the non-native range. Prosopis juliflora is taxonomically complicated and with P. pallida forms the P. juliflora complex. Thus we sampled P. juliflora in its native Venezuela, and also located two field sites in Peru, the native range of Prosopis pallida. Canopies of Prosopis juliflora, a native of the New World but an invader in many other regions, had facilitative effects on the diversity of other species in its native Venezuela, and P. pallida had both negative and positive effects depending on the year, (overall neutral effects) in its native Peru. However, in India and Hawaii, USA, where P. juliflora is an aggressive invader, canopy effects were consistently and strongly negative on species richness. Prosopis cineraria, a native to India, had much weaker effects on species richness in India than P. juliflora. We carried out multiple congeneric comparisons between P. juliflora and P. cineraria, and found that soil from the rhizosphere of P. juliflora had higher extractable phosphorus, soluble salts and total phenolics than P. cineraria rhizosphere soils. Experimentally applied P. juliflora litter caused far greater mortality of native Indian species than litter from P. cineraria. Prosopis juliflora leaf leachate had neutral to negative effects on root growth of three common crop species of north-west India whereas P. cineraria leaf leachate had positive effects. Prosopis juliflora leaf leachate also had higher concentrations of total phenolics and L-tryptophan than P. cineraria, suggesting a potential allelopathic mechanism for the congeneric differences. Our results also suggest the possibility of regional evolutionary trajectories among competitors and that recent mixing of species from different trajectories has the potential to disrupt evolved interactions among native species

Conformational changes and CO2-induced channel gating in connexin26

Connexins form large-pore channels that function either as dodecameric gap junctions or hexameric hemichannels to allow the regulated movement of small molecules and ions across cell membranes. Opening or closing of the channels is controlled by a variety of stimuli, and dysregulation leads to multiple diseases. An increase in the partial pressure of carbon dioxide (PCO2) has been shown to cause connexin26 (Cx26) gap junctions to close. Here, we use cryoelectron microscopy (cryo-EM) to determine the structure of human Cx26 gap junctions under increasing levels of PCO2. We show a correlation between the level of PCO2 and the size of the aperture of the pore, governed by the N-terminal helices that line the pore. This indicates that CO2 alone is sufficient to cause conformational changes in the protein. Analysis of the conformational states shows that movements at the N terminus are linked to both subunit rotation and flexing of the transmembrane helices

Cryo-EM structural analysis of FADD:Caspase-8 complexes defines the catalytic dimer architecture for co-ordinated control of cell fate.

Regulated cell death is essential in development and cellular homeostasis. Multi-protein platforms, including the Death-Inducing Signaling Complex (DISC), co-ordinate cell fate via a core FADD:Caspase-8 complex and its regulatory partners, such as the cell death inhibitor c-FLIP. Here, using electron microscopy, we visualize full-length procaspase-8 in complex with FADD. Our structural analysis now reveals how the FADD-nucleated tandem death effector domain (tDED) helical filament is required to orientate the procaspase-8 catalytic domains, enabling their activation via anti-parallel dimerization. Strikingly, recruitment of c-FLIPS into this complex inhibits Caspase-8 activity by altering tDED triple helix architecture, resulting in steric hindrance of the canonical tDED Type I binding site. This prevents both Caspase-8 catalytic domain assembly and tDED helical filament elongation. Our findings reveal how the plasticity, composition and architecture of the core FADD:Caspase-8 complex critically defines life/death decisions not only via the DISC, but across multiple key signaling platforms including TNF complex II, the ripoptosome, and RIPK1/RIPK3 necrosome

Cryo-EM structural analysis of FADD:Caspase-8 complexes defines the catalytic dimer architecture for co-ordinated control of cell fate.

Regulated cell death is essential in development and cellular homeostasis. Multi-protein platforms, including the Death-Inducing Signaling Complex (DISC), co-ordinate cell fate via a core FADD:Caspase-8 complex and its regulatory partners, such as the cell death inhibitor c-FLIP. Here, using electron microscopy, we visualize full-length procaspase-8 in complex with FADD. Our structural analysis now reveals how the FADD-nucleated tandem death effector domain (tDED) helical filament is required to orientate the procaspase-8 catalytic domains, enabling their activation via anti-parallel dimerization. Strikingly, recruitment of c-FLIPS into this complex inhibits Caspase-8 activity by altering tDED triple helix architecture, resulting in steric hindrance of the canonical tDED Type I binding site. This prevents both Caspase-8 catalytic domain assembly and tDED helical filament elongation. Our findings reveal how the plasticity, composition and architecture of the core FADD:Caspase-8 complex critically defines life/death decisions not only via the DISC, but across multiple key signaling platforms including TNF complex II, the ripoptosome, and RIPK1/RIPK3 necrosome

Open and closed conformations of a sub-80 kDa Chagas vaccine candidate defined by a cryo-EM led integrative approach

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, remains a significant global public health concern. It affects an estimated eight million individuals worldwide, with the majority remaining undiagnosed. Despite its profound health impact in both endemic and non-endemic areas, no vaccine is available, and the existing therapies are outdated, producing severe side effects.The 80kDa prolyl oligopeptidase of Trypanosoma cruzi (TcPOP) has been recently identified as a leading candidate for Chagas vaccine development. However, its three-dimensional structure has remained elusive for almost two decades since its discovery. We report the first three-dimensional structure of TcPOP in open and closed conformation, at a resolution of 3.0 and 2.5 Angstroms respectively, determined using single-particle cryo-electron microscopy. Multiple conformations were observed and were further characterized, using plasmonic optical tweezers.To assess the immunogenic potential of TcPOP, we immunized mice and evaluated both polyclonal and monoclonal responses against the TcPOP antigen and its homologues. The results revealed unexpected cross-reactivity across prolyl POPs from other closely related parasites, but intriguingly, not towards the human homologue.Altogether, our findings provide critical structural insights necessary to understand the immunogenicity of TcPOP for future Chagas vaccine development and diagnostic applications.Additionally, our integrative approach indicated that stage-tilted acquisition can yield biologically relevant information for challenging sub-80kDa proteins and could adequately resolve the cryoEM structures. Consequently, this comprehensive strategy can significantly enhance the success rate in determining the structures of proteins that present challenges in characterization

Cryo-EM structure of human Pol κ bound to DNA and mono-ubiquitylated PCNA

Y-family DNA polymerase κ (Pol κ) can replicate damaged DNA templates to rescue stalled replication forks. Access of Pol κ to DNA damage sites is facilitated by its interaction with the processivity clamp PCNA and is regulated by PCNA mono-ubiquitylation. Here, we present cryo-EM reconstructions of human Pol κ bound to DNA, an incoming nucleotide, and wild type or mono-ubiquitylated PCNA (Ub-PCNA). In both reconstructions, the internal PIP-box adjacent to the Pol κ Polymerase-Associated Domain (PAD) docks the catalytic core to one PCNA protomer in an angled orientation, bending the DNA exiting the Pol κ active site through PCNA, while Pol κ C-terminal domain containing two Ubiquitin Binding Zinc Fingers (UBZs) is invisible, in agreement with disorder predictions. The ubiquitin moieties are partly flexible and extend radially away from PCNA, with the ubiquitin at the Pol κ-bound protomer appearing more rigid. Activity assays suggest that, when the internal PIP-box interaction is lost, Pol κ is retained on DNA by a secondary interaction between the UBZs and the ubiquitins flexibly conjugated to PCNA. Our data provide a structural basis for the recruitment of a Y-family TLS polymerase to sites of DNA damage.This research was supported by King Abdullah University of Science and Technology through core funding (to S.M.H.) and the Competitive Research Award Grant CRG8 URF/1/4036‐01‐01 (to S.M.H. and A.D.B.), and by the Wellcome Trust (to A.D.B.). R.C. acknowledges funding from the MINECO (CTQ2016-78636-P) and to AGAUR, (2017 SGR 324). The MD project has been carried out using CSUC resources. We acknowledge The Midlands Regional Cryo-EM Facility at the Leicester Institute of Structural and Chemical Biology (LISCB), major funding from MRC (MC_PC_17136). We thank Christos Savva (LISCB, University of Leicester) for his help in cryo-EM data collection and advice on data processing.Peer reviewe

Automatically Harnessing Sparse Acceleration

Sparse linear algebra is central to many scientific programs, yet compilers fail to optimize it well. High-performance libraries are available, but adoption costs are significant. Moreover, libraries tie programs into vendor-specific software and hardware ecosystems, creating non-portable code. In this paper, we develop a new approach based on our specification Language for implementers of Linear Algebra Computations (LiLAC). Rather than requiring the application developer to (re)write every program for a given library, the burden is shifted to a one-off description by the library implementer. The LiLAC-enabled compiler uses this to insert appropriate library routines without source code changes. LiLAC provides automatic data marshaling, maintaining state between calls and minimizing data transfers. Appropriate places for library insertion are detected in compiler intermediate representation, independent of source languages. We evaluated on large-scale scientific applications written in FORTRAN; standard C/C++ and FORTRAN benchmarks; and C++ graph analytics kernels. Across heterogeneous platforms, applications and data sets we show speedups of 1.1$\times$ to over 10$\times$ without user intervention.Comment: Accepted to CC 202