Bayesian modelling of compositional heterogeneity in molecular phylogenetics

In molecular phylogenetics, standard models of sequence evolution generally assume that sequence composition remains constant over evolutionary time. However, this assumption is violated in many datasets which show substantial heterogeneity in sequence composition across taxa. We propose a model which allows compositional heterogeneity across branches, and formulate the model in a Bayesian framework. Specifically, the root and each branch of the tree is associated with its own composition vector whilst a global matrix of exchangeability parameters applies everywhere on the tree. We encourage borrowing of strength between branches by developing two possible priors for the composition vectors: one in which information can be exchanged equally amongst all branches of the tree and another in which more information is exchanged between neighbouring branches than between distant branches. We also propose a Markov chain Monte Carlo (MCMC) algorithm for posterior inference which uses data augmentation of substitutional histories to yield a simple complete data likelihood function that factorises over branches and allows Gibbs updates for most parameters. Standard phylogenetic models are not informative about the root position. Therefore a significant advantage of the proposed model is that it allows inference about rooted trees. The position of the root is fundamental to the biological interpretation of trees, both for polarising trait evolution and for establishing the order of divergence among lineages. Furthermore, unlike some other related models from the literature, inference in the model we propose can be carried out through a simple MCMC scheme which does not require problematic dimension-changing moves. We investigate the performance of the model and priors in analyses of two alignments for which there is strong biological opinion about the tree topology and root position.

Plastid establishment did not require a chlamydial partner

Primary plastids descend from the cyanobacterial endosymbiont of an ancient eukaryotic host, but the initial selective drivers that stabilized the association between these two cells are still unclear. One hypothesis that has achieved recent prominence suggests that the first role of the cyanobiont was in energy provision for a host cell whose reserves were being depleted by an intracellular chlamydial pathogen. A pivotal claim is that it was chlamydial proteins themselves that converted otherwise unusable cyanobacterial metabolites into host energy stores. We test this hypothesis by investigating the origins of the key enzymes using sophisticated phylogenetics. Here we show a mosaic origin for the relevant pathway combining genes with host, cyanobacterial or bacterial ancestry, but we detect no strong case for Chlamydiae to host transfer under the best-fitting models. Our conclusion is that there is no compelling evidence from gene trees that Chlamydiae played any role in establishing the primary plastid endosymbiosis.

Primary plastids descend from an endosymbiosis involving cyanobacteria, an ancient eukaryotic host and, possibly, a chlamydial pathogen. Here, Domman and colleagues use sophisticated phylogenetic methods to show that Chlamydiae did not play a role in establishing the primary plastid endosymbiosis.

Archaeal "dark matter" and the origin of eukaryotes

Current hypotheses about the history of cellular life are mainly based on analyses of cultivated organisms, but these represent only a small fraction of extant biodiversity. The sequencing of new environmental lineages therefore provides an opportunity to test, revise, or reject existing ideas about the tree of life and the origin of eukaryotes. According to the textbook three domains hypothesis, the eukaryotes emerge as the sister group to a monophyletic Archaea. However, recent analyses incorporating better phylogenetic models and an improved sampling of the archaeal domain have generally supported the competing eocyte hypothesis, in which core genes of eukaryotic cells originated from within the Archaea, with important implications for eukaryogenesis. Given this trend, it was surprising that a recent analysis incorporating new genomes from uncultivated Archaea recovered a strongly supported three domains tree. Here, we show that this result was due in part to the use of a poorly fitting phylogenetic model and also to the inclusion by an automated pipeline of genes of putative bacterial origin rather than nucleocytosolic versions for some of the eukaryotes analyzed. When these issues were resolved, analyses including the new archaeal lineages placed core eukaryotic genes within the Archaea. These results are consistent with a number of recent studies in which improved archaeal sampling and better phylogenetic models agree in supporting the eocyte tree over the three domains hypothesis.

Evolution: rooting the eukaryotic tree of life

The root of the eukaryotic tree is a major unresolved question in evolutionary biology. A recent study marshals mitochondrial genes to place that root between the enigmatic Excavates and all other eukaryotes, providing an interesting new perspective on early eukaryotic evolution.

Reduction and expansion in microsporidian genome evolution: new insights from comparative genomics

Microsporidia are an abundant group of obligate intracellular parasites of other eukaryotes, including immunocompromised humans, but the molecular basis of their intracellular lifestyle and pathobiology are poorly understood. New genomes from a taxonomically broad range of microsporidians, complemented by published expression data, provide an opportunity for comparative analyses to identify conserved and lineage-specific patterns of microsporidian genome evolution that have underpinned this success. In this study, we infer that a dramatic bottleneck in the last common microsporidian ancestor (LCMA) left a small conserved core of genes that was subsequently embellished by gene family expansion driven by gene acquisition in different lineages. Novel expressed protein families represent a substantial fraction of sequenced microsporidian genomes and are significantly enriched for signals consistent with secretion or membrane location. Further evidence of selection is inferred from the gain and reciprocal loss of functional domains between paralogous genes, for example, affecting transport proteins. Gene expansions among transporter families preferentially affect those that are located on the plasma membrane of model organisms, consistent with recruitment to plug conserved gaps in microsporidian biosynthesis and metabolism. Core microsporidian genes shared with other eukaryotes are enriched in orthologs that, in yeast, are highly expressed, highly connected, and often essential, consistent with strong negative selection against further reduction of the conserved gene set since the LCMA. Our study reveals that microsporidian genome evolution is a highly dynamic process that has balanced constraint, reductive evolution, and genome expansion during adaptation to an extraordinarily successful obligate intracellular lifestyle.

An archaeal origin of eukaryotes supports only two primary domains of life

The discovery of the Archaea and the proposal of the three-domains 'universal' tree, based on ribosomal RNA and core genes mainly involved in protein translation, catalysed new ideas for cellular evolution and eukaryotic origins. However, accumulating evidence suggests that the three-domains tree may be incorrect: evolutionary trees made using newer methods place eukaryotic core genes within the Archaea, supporting hypotheses in which an archaeon participated in eukaryotic origins by founding the host lineage for the mitochondrial endosymbiont. These results provide support for only two primary domains of life--Archaea and Bacteria--because eukaryotes arose through partnership between them.

Subfunctionalization of cyprinid hypoxia-inducible factors for roles in development and oxygen sensing

Among vertebrates, teleost fishes have evolved the most impressive adaptations to variable oxygen tensions in water (Shoubridge and Hochachka 1980; Nilsson and Randall 2010). Under conditions of oxygen deprivation (hypoxia), major changes in gene expression are mediated by hypoxia-inducible factors (HIF alpha). Here we show that hif alpha genes were duplicated in the teleost specific whole-genome duplication. Although one of each paralogous gene pair was lost in most teleosts, both copies were retained in cyprinids. Computational analyses suggest that these duplicates have become subfunctionalized with complementary changes in coding and regulatory sequences within each paralogous gene pair. We tested our predictions with comparisons of hif alpha transcription in zebrafish, a cyprinid, and sturgeon, an outgroup that diverged from teleosts before the duplication event. Our experiments revealed distinct transcriptional profiles in the cyprinid duplicates: while one of each paralogous pair maintained the ancestral developmental response, the other was more sensitive to changes in oxygen tension. These results demonstrate the subfunctionalization of cyprinid hif alpha paralogs for specialized roles in development and the hypoxic stress response.

A congruent phylogenomic signal places eukaryotes within the Archaea

Determining the relationships among the major groups of cellular life is important for understanding the evolution of biological diversity, but is difficult given the enormous time spans involved. In the textbook ‘three domains’ tree based on informational genes, eukaryotes and Archaea share a common ancestor to the exclusion of Bacteria. However, some phylogenetic analyses of the same data have placed eukaryotes within the Archaea, as the nearest relatives of different archaeal lineages. We compared the support for these competing hypotheses using sophisticated phylogenetic methods and an improved sampling of archaeal biodiversity. We also employed both new and existing tests of phylogenetic congruence to explore the level of uncertainty and conflict in the data. Our analyses suggested that much of the observed incongruence is weakly supported or associated with poorly fitting evolutionary models. All of our phylogenetic analyses, whether on small subunit and large subunit ribosomal RNA or concatenated protein-coding genes, recovered a monophyletic group containing eukaryotes and the TACK archaeal superphylum comprising the Thaumarchaeota, Aigarchaeota, Crenarchaeota and Korarchaeota. Hence, while our results provide no support for the iconic three-domain tree of life, they are consistent with an extended eocyte hypothesis whereby vital components of the eukaryotic nuclear lineage originated from within the archaeal radiation.

Proteome-wide analysis of functional divergence in bacteria: exploring a host of ecological adaptations

Functional divergence is the process by which new genes and functions originate through the modification of existing ones. Both genetic and environmental factors influence the evolution of new functions, including gene duplication or changes in the ecological requirements of an organism. Novel functions emerge at the expense of ancestral ones and are generally accompanied by changes in the selective forces at constrained protein regions. We present software capable of analyzing whole proteomes, identifying putative amino acid replacements leading to functional change in each protein and performing statistical tests on all tabulated data. We apply this method to 750 complete bacterial proteomes to identify high-level patterns of functional divergence and link these patterns to ecological adaptations. Proteome-wide analyses of functional divergence in bacteria with different ecologies reveal a separation between proteins involved in information processing (Ribosome biogenesis etc.) and those which are dependent on the environment (energy metabolism, defense etc.). We show that the evolution of pathogenic and symbiotic bacteria is constrained by their association with the host, and also identify unusual events of functional divergence even in well-studied bacteria such as Escherichia coli. We present a description of the roles of phylogeny and ecology in functional divergence at the level of entire proteomes in bacteria.

Aldosterone does not modify gene expression in human endothelial cells

The toxic effects of aldosterone on the vasculature, and in particular on the endothelial layer, have been proposed as having an important role in the cardiovascular pathology observed in mineralocorticoid-excess states. In order to characterize the genomic molecular mechanisms driving the aldosterone-induced endothelial dysfunction, we performed an expression microarray on transcripts obtained from both human umbilical vein endothelial cells and human coronary artery endothelial cells stimulated with 10 - 7 M aldosterone for 18 h. The results were then subjected to qRT-PCR confirmation, also including a group of genes known to be involved in the control of the endothelial function or previously described as regulated by aldosterone. The state of activation of the mineralocorticoid receptor was investigated by means of a luciferase-reporter assay using a plasmid encoding a mineralocorticoid and glucocorticoid-sensitive promoter. Aldosterone did not determine any significant change in gene expression in either cell type both in the microarray and in the qRT-PCR analysis. The luciferase-reporter assay showed no activation of the mineralocorticoid receptor following aldosterone stimulation. The status of nonfunctionality of the mineralocorticoid receptor expressed in cultured human umbilical and coronary artery endothelial cells does not allow aldosterone to modify gene expression and provides evidence against either a beneficial or harmful genomic effect of aldosterone on healthy endothelial cells.

Outcome of ranibizumab treatment in neovascular age related macula degeneration in eyes with baseline visual acuity better than 6/12

BACKGROUND

The beneficial effect of intravitreal ranibizumab in the treatment of neovascular age-related macula degeneration (nAMD) is well known. Outcome data for eyes presenting with visual acuity better than 6/12 is limited.

AIMS

To assess the effect of baseline vision on outcome in ranibizumab-treated nAMD eyes, including a subgroup with baseline vision ≥6/12 (<0.30  logmar).

DESIGN

Prospective, consecutive and interventional case series.

METHODS

A consecutive cohort of patients treated with intravitreal ranibizumab for nAMD with 52-week follow-up were studied. Patients who had received previous treatment for nAMD were excluded. Eyes were stratified according to baseline logmar visual acuity into four groups: <0.30 (>6/12), 0.30-0.59 (6/12-6/24), 0.60-0.99 (6/24-6/60) and 1.00-1.20 (6/60-6/96). Intravitreal ranibizumab (0.5 mg in 0.05 ml) was administered in three loading monthly doses followed by PRN dosing according to optical coherence tomography (OCT) findings.

RESULTS

A total of 615 eyes were studied including 88 eyes with baseline vision <0.30. The mean change in logmar letters at 52 weeks was +5.5 (entire study group), -0.5 (<0.30 subgroup), +2.2 (0.30-0.59 subgroup), +6.5 (0.60-0.99 subgroup) and +15.3 (1.00-1.20 subgroup). In the <0.30 subgroup, 60 of 88 eyes (68%) had best-corrected visual acuity (BCVA) equal to or better than baseline and 82 of 88 eyes (93%) lost <15 letters at 52 weeks. Within this subgroup 56 of 67 eyes (84%) maintained UK driving standard BCVA visual acuity over the study period.

CONCLUSIONS

This study provides evidence that intravitreal ranibizumab treatment stabilises good vision in nAMD presenting with vision better than 6/12 over 52 weeks follow-up.

Informational gene phylogenies do not support a fourth domain of life for nucleocytoplasmic large DNA viruses

Mimivirus is a nucleocytoplasmic large DNA virus (NCLDV) with a genome size (1.2 Mb) and coding capacity ( 1000 genes) comparable to that of some cellular organisms. Unlike other viruses, Mimivirus and its NCLDV relatives encode homologs of broadly conserved informational genes found in Bacteria, Archaea, and Eukaryotes, raising the possibility that they could be placed on the tree of life. A recent phylogenetic analysis of these genes showed the NCLDVs emerging as a monophyletic group branching between Eukaryotes and Archaea. These trees were interpreted as evidence for an independent "fourth domain" of life that may have contributed DNA processing genes to the ancestral eukaryote. However, the analysis of ancient evolutionary events is challenging, and tree reconstruction is susceptible to bias resulting from non-phylogenetic signals in the data. These include compositional heterogeneity and homoplasy, which can lead to the spurious grouping of compositionally-similar or fast-evolving sequences. Here, we show that these informational gene alignments contain both significant compositional heterogeneity and homoplasy, which were not adequately modelled in the original analysis. When we use more realistic evolutionary models that better fit the data, the resulting trees are unable to reject a simple null hypothesis in which these informational genes, like many other NCLDV genes, were acquired by horizontal transfer from eukaryotic hosts. Our results suggest that a fourth domain is not required to explain the available sequence data.

The effect of chaperonin buffering on protein evolution

Molecular chaperones are highly conserved and ubiquitous proteins that help other proteins in the cell to fold. Pioneering work by Rutherford and Lindquist suggested that the chaperone Hsp90 could buffer (i.e., suppress) phenotypic variation in its client proteins and that alternate periods of buffering and expression of these variants might be important in adaptive evolution. More recently, Tokuriki and Tawfik presented an explicit mechanism for chaperone-dependent evolution, in which the Escherichia coli chaperonin GroEL facilitated the folding of clients that had accumulated structurally destabilizing but neofunctionalizing mutations in the protein core. But how important an evolutionary force is chaperonin-mediated buffering in nature? Here, we address this question by modeling the per-residue evolutionary rate of the crystallized E. coli proteome, evaluating the relative contributions of chaperonin buffering, functional importance, and structural features such as residue contact density. Previous findings suggest an interaction between codon bias and GroEL in limiting the effects of misfolding errors. Our results suggest that the buffering of deleterious mutations by GroEL increases the evolutionary rate of client proteins. We then examine the evolutionary fate of GroEL clients in the Mycoplasmas, a group of bacteria containing the only known organisms that lack chaperonins. We show that GroEL was lost once in the common ancestor of a monophyletic subgroup of Mycoplasmas, and we evaluate the effect of this loss on the subsequent evolution of client proteins, providing evidence that client homologs in 11 Mycoplasma species have lost their obligate dependency on GroEL for folding. Our analyses indicate that individual molecules such as chaperonins can have significant effects on proteome evolution through their modulation of protein folding.

Effect of length of stay in intensive care unit on hospital and long-term mortality of critically ill adult patients

BACKGROUND

Critical illness leading to prolonged length of stay (LOS) in an intensive care unit (ICU) is associated with significant mortality and resource utilization. This study assessed the independent effect of ICU LOS on in-hospital and long-term mortality after hospital discharge.

METHODS

Clinical and mortality data of 22 298 patients, aged 16 yr and older, admitted to ICU between 1987 and 2002 were included in this linked-data cohort study. Cox's regression with restricted cubic spline function was used to model the effect of LOS on in-hospital and long-term mortality after adjusting for age, gender, acute physiology score (APS), maximum number of organ failures, era of admission, elective admission, Charlson's co-morbidity index, and diagnosis. The variability each predictor explained was calculated by the percentage of the chi(2) statistic contribution to the total chi(2) statistic.

RESULTS

Most hospital deaths occurred within the first few days of ICU admission. Increasing LOS in ICU was not associated with an increased risk of in-hospital mortality after adjusting for other covariates, but was associated with an increased risk of long-term mortality after hospital discharge. The variability on the long-term mortality effect associated with ICU LOS (2.3%) appeared to reach a plateau after the first 10 days in ICU and was not as important as age (35.8%), co-morbidities (18.6%), diagnosis (10.9%), and APS (3.6%).

CONCLUSIONS

LOS in ICU was not an independent risk factor for in-hospital mortality, but it had a small effect on long-term mortality after hospital discharge after adjustment for other risk factors.

Two chaperonin systems in bacterial genomes with distinct ecological roles

Bacterial chaperonins are essential to cell viability and have a role in endosymbiosis, which leads to increased biological complexity. However, the extent to which chaperonins promote ecological innovation is unknown. We screened 622 bacterial genomes for genes encoding chaperonins, and found archaeal-like chaperonins in bacteria that inhabit archaeal ecological niches. We found that chaperonins encoded in pathogenic bacteria are the most functionally divergent. We identified the molecular basis of the dramatic structural changes in mitochondrial GROEL, a highly derived chaperonin gene. Our analysis suggests that chaperonins are important capacitors of evolutionary and ecological change.

Effect of an episode of critical illness on subsequent hospitalisation: a linked data study

Healthcare utilisation can affect quality of life and is important in assessing the cost-effectiveness of medical interventions. A clinical database was linked to two Australian state administrative databases to assess the difference in incidence of healthcare utilisation of 19,921 patients who survived their first episode of critical illness. The number of hospital admissions and days of hospitalisation per patient-year was respectively 150% and 220% greater after than before an episode of critical illness (assessed over the same time period). This was the case regardless of age or type of surgery (i.e. cardiac vs non-cardiac). After adjusting for the ageing effect of the cohort as a whole, there was still an unexplained two to four-fold increase in hospital admissions per patient-year after an episode of critical illness. We conclude that an episode of critical illness is a robust predictor of subsequent healthcare utilisation.

Genome-wide functional divergence after the symbiosis of proteobacteria with insects unraveled through a novel computational approach

Symbiosis has been among the most important evolutionary steps to generate biological complexity. The establishment of symbiosis required an intimate metabolic link between biological systems with different complexity levels. The strict endo-cellular symbiotic bacteria of insects are beautiful examples of the metabolic coupling between organisms belonging to different kingdoms, a eukaryote and a prokaryote. The host (eukaryote) provides the endosymbiont (prokaryote) with a stable cellular environment while the endosymbiont supplements the host's diet with essential metabolites. For such communication to take place, endosymbionts' genomes have suffered dramatic modifications and reconfigurations of proteins' functions. Two of the main modifications, loss of genes redundant for endosymbiotic bacteria or the host and bacterial genome streamlining, have been extensively studied. However, no studies have accounted for possible functional shifts in the endosymbiotic proteomes. Here, we develop a simple method to screen genomes for evidence of functional divergence between two species clusters, and we apply it to identify functional shifts in the endosymbiotic proteomes. Despite the strong effects of genetic drift in the endosymbiotic systems, we unexpectedly identified genes to be under stronger selective constraints in endosymbionts of aphids and ants than in their free-living bacterial relatives. These genes are directly involved in supplementing the host's diet with essential metabolites. A test of functional divergence supports a strong relationship between the endosymbiosis and the functional shifts of proteins involved in the metabolic communication with the insect host. The correlation between functional divergence in the endosymbiotic bacterium and the ecological requirements of the host uncovers their intimate biochemical and metabolic communication and provides insights on the role of symbiosis in generating species diversity.

Biological complexity has emerged on earth by the combination of living forms. This combination, called symbiosis, had to overcome the problems caused by the uncoupled metabolisms of the organisms involved. One way to do so was through the loss of genes that were no longer needed for the endosymbiont in the protected cellular environment provided by the host. Another step necessary to adjust both metabolisms was through the change in the function of bacterial proteins to perform new roles in the symbiotic system. In this article, we test such events in symbiotic systems involving an insect and a bacterium by developing a new and simple method to identify proteome-wide functional shifts. Our results show that most of the functional changes occurred at genes involved in metabolic communication with the host and are correlated with the host's ecological traits.

Conceptual basis and methodology of the SOPHIA study

To clarify the role of gene polymorphisms on the effect of losartan and losartan plus hydrochlorothiazide on blood pressure (primary end point) and on cardiac, vascular and metabolic phenotypes (secondary end point) after 4, 8, 12, 16 and 48 weeks treatment, an Italian collaborative study - The Study of the Pharmacogenomics in Italian hypertensive patients treated with the Angiotensin receptor blocker losartan (SOPHIA) - on never-treated essential hypertensives (n = 800) was planned. After an 8 week run-in, losartan 50 mg once daily will be given and doubled to 100 mg at week +4 if blood pressure is more than 140/90 mmHg. Hydroclorothiazide 25 mg once daily at week +8 and amlodipine 5 mg at week +16 will be added if blood pressure is more than 140/90 mmHg. Cardiac mass (echocardiography), carotid intima-media thickness, 24 h ambulatory blood pressure, homeostatic model assessment (HOMA) index, microalbuminuria, plasma renin activity and aldosterone, endogenous lithium clearance, brain natriuretic peptide and losartan metabolites will be evaluated. Genes of the renin-angiotensin-aldosterone system, salt sensitivity, the beta-adrenergic system and losartan metabolism will be studied (Illumina custom arrays). A whole-genome scan will also be performed in half of the study cohort (1M array, Illumina 500 GX beadstation).