Evolving Sensitivity Balances Boolean Networks

We investigate the sensitivity of Boolean Networks (BNs) to mutations. We are interested in Boolean Networks as a model of Gene Regulatory Networks (GRNs). We adopt Ribeiro and Kauffman’s Ergodic Set and use it to study the long term dynamics of a BN. We define the sensitivity of a BN to be the mean change in its Ergodic Set structure under all possible loss of interaction mutations. Insilico experiments were used to selectively evolve BNs for sensitivity to losing interactions. We find that maximum sensitivity was often achievable and resulted in the BNs becoming topologically balanced, i.e. they evolve towards network structures in which they have a similar number of inhibitory and excitatory interactions. In terms of the dynamics, the dominant sensitivity strategy that evolved was to build BNs with Ergodic Sets dominated by a single long limit cycle which is easily destabilised by mutations. We discuss the relevance of our findings in the context of Stem Cell Differentiation and propose a relationship between pluripotent stem cells and our evolved sensitive networks

Cell cycle-independent phospho-regulation of Fkh2 during hyphal growth regulates Candida albicans pathogenesis.

The opportunistic human fungal pathogen, Candida albicans, undergoes morphological and transcriptional adaptation in the switch from commensalism to pathogenicity. Although previous gene-knockout studies have identified many factors involved in this transformation, it remains unclear how these factors are regulated to coordinate the switch. Investigating morphogenetic control by post-translational phosphorylation has generated important regulatory insights into this process, especially focusing on coordinated control by the cyclin-dependent kinase Cdc28. Here we have identified the Fkh2 transcription factor as a regulatory target of both Cdc28 and the cell wall biosynthesis kinase Cbk1, in a role distinct from its conserved function in cell cycle progression. In stationary phase yeast cells 2D gel electrophoresis shows that there is a diverse pool of Fkh2 phospho-isoforms. For a short window on hyphal induction, far before START in the cell cycle, the phosphorylation profile is transformed before reverting to the yeast profile. This transformation does not occur when stationary phase cells are reinoculated into fresh medium supporting yeast growth. Mass spectrometry and mutational analyses identified residues phosphorylated by Cdc28 and Cbk1. Substitution of these residues with non-phosphorylatable alanine altered the yeast phosphorylation profile and abrogated the characteristic transformation to the hyphal profile. Transcript profiling of the phosphorylation site mutant revealed that the hyphal phosphorylation profile is required for the expression of genes involved in pathogenesis, host interaction and biofilm formation. We confirmed that these changes in gene expression resulted in corresponding defects in pathogenic processes. Furthermore, we identified that Fkh2 interacts with the chromatin modifier Pob3 in a phosphorylation-dependent manner, thereby providing a possible mechanism by which the phosphorylation of Fkh2 regulates its specificity. Thus, we have discovered a novel cell cycle-independent phospho-regulatory event that subverts a key component of the cell cycle machinery to a role in the switch from commensalism to pathogenicity

A large genome-wide association study of age-related macular degeneration highlights contributions of rare and common variants.

This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/ng.3448Advanced age-related macular degeneration (AMD) is the leading cause of blindness in the elderly, with limited therapeutic options. Here we report on a study of >12 million variants, including 163,714 directly genotyped, mostly rare, protein-altering variants. Analyzing 16,144 patients and 17,832 controls, we identify 52 independently associated common and rare variants (P < 5 × 10(-8)) distributed across 34 loci. Although wet and dry AMD subtypes exhibit predominantly shared genetics, we identify the first genetic association signal specific to wet AMD, near MMP9 (difference P value = 4.1 × 10(-10)). Very rare coding variants (frequency <0.1%) in CFH, CFI and TIMP3 suggest causal roles for these genes, as does a splice variant in SLC16A8. Our results support the hypothesis that rare coding variants can pinpoint causal genes within known genetic loci and illustrate that applying the approach systematically to detect new loci requires extremely large sample sizes.We thank all participants of all the studies included for enabling this research by their participation in these studies. Computer resources for this project have been provided by the high-performance computing centers of the University of Michigan and the University of Regensburg. Group-specific acknowledgments can be found in the Supplementary Note. The Center for Inherited Diseases Research (CIDR) Program contract number is HHSN268201200008I. This and the main consortium work were predominantly funded by 1X01HG006934-01 to G.R.A. and R01 EY022310 to J.L.H

Sloan Digital Sky Survey IV: Mapping the Milky Way, Nearby Galaxies, and the Distant Universe

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratios in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially resolved spectroscopy for thousands of nearby galaxies (median $z\sim 0.03$). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between $z\sim 0.6$ and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGNs and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5 m Sloan Foundation Telescope at the Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5 m du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in 2016 July

The Fifteenth Data Release of the Sloan Digital Sky Surveys: First Release of MaNGA-derived Quantities, Data Visualization Tools, and Stellar Library

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (2014 July–2017 July). This is the third data release for SDSS-IV, and the 15th from SDSS (Data Release Fifteen; DR15). New data come from MaNGA—we release 4824 data cubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g., stellar and gas kinematics, emission-line and other maps) from the MaNGA Data Analysis Pipeline, and a new data visualization and access tool we call "Marvin." The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper, we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials, and examples of data use. Although SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020–2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data

Mortality from gastrointestinal congenital anomalies at 264 hospitals in 74 low-income, middle-income, and high-income countries: a multicentre, international, prospective cohort study

Background: Congenital anomalies are the fifth leading cause of mortality in children younger than 5 years globally. Many gastrointestinal congenital anomalies are fatal without timely access to neonatal surgical care, but few studies have been done on these conditions in low-income and middle-income countries (LMICs). We compared outcomes of the seven most common gastrointestinal congenital anomalies in low-income, middle-income, and high-income countries globally, and identified factors associated with mortality. // Methods: We did a multicentre, international prospective cohort study of patients younger than 16 years, presenting to hospital for the first time with oesophageal atresia, congenital diaphragmatic hernia, intestinal atresia, gastroschisis, exomphalos, anorectal malformation, and Hirschsprung's disease. Recruitment was of consecutive patients for a minimum of 1 month between October, 2018, and April, 2019. We collected data on patient demographics, clinical status, interventions, and outcomes using the REDCap platform. Patients were followed up for 30 days after primary intervention, or 30 days after admission if they did not receive an intervention. The primary outcome was all-cause, in-hospital mortality for all conditions combined and each condition individually, stratified by country income status. We did a complete case analysis. // Findings: We included 3849 patients with 3975 study conditions (560 with oesophageal atresia, 448 with congenital diaphragmatic hernia, 681 with intestinal atresia, 453 with gastroschisis, 325 with exomphalos, 991 with anorectal malformation, and 517 with Hirschsprung's disease) from 264 hospitals (89 in high-income countries, 166 in middle-income countries, and nine in low-income countries) in 74 countries. Of the 3849 patients, 2231 (58·0%) were male. Median gestational age at birth was 38 weeks (IQR 36–39) and median bodyweight at presentation was 2·8 kg (2·3–3·3). Mortality among all patients was 37 (39·8%) of 93 in low-income countries, 583 (20·4%) of 2860 in middle-income countries, and 50 (5·6%) of 896 in high-income countries (p<0·0001 between all country income groups). Gastroschisis had the greatest difference in mortality between country income strata (nine [90·0%] of ten in low-income countries, 97 [31·9%] of 304 in middle-income countries, and two [1·4%] of 139 in high-income countries; p≤0·0001 between all country income groups). Factors significantly associated with higher mortality for all patients combined included country income status (low-income vs high-income countries, risk ratio 2·78 [95% CI 1·88–4·11], p<0·0001; middle-income vs high-income countries, 2·11 [1·59–2·79], p<0·0001), sepsis at presentation (1·20 [1·04–1·40], p=0·016), higher American Society of Anesthesiologists (ASA) score at primary intervention (ASA 4–5 vs ASA 1–2, 1·82 [1·40–2·35], p<0·0001; ASA 3 vs ASA 1–2, 1·58, [1·30–1·92], p<0·0001]), surgical safety checklist not used (1·39 [1·02–1·90], p=0·035), and ventilation or parenteral nutrition unavailable when needed (ventilation 1·96, [1·41–2·71], p=0·0001; parenteral nutrition 1·35, [1·05–1·74], p=0·018). Administration of parenteral nutrition (0·61, [0·47–0·79], p=0·0002) and use of a peripherally inserted central catheter (0·65 [0·50–0·86], p=0·0024) or percutaneous central line (0·69 [0·48–1·00], p=0·049) were associated with lower mortality. // Interpretation: Unacceptable differences in mortality exist for gastrointestinal congenital anomalies between low-income, middle-income, and high-income countries. Improving access to quality neonatal surgical care in LMICs will be vital to achieve Sustainable Development Goal 3.2 of ending preventable deaths in neonates and children younger than 5 years by 2030

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects

Functionality and sensing in Boolean networks

The main theme of this thesis is investigating how a cell’s biological function relates to the topology of its Gene Regulatory Network (GRN). In this context, the limits a biological function places on evolution are examined and also whether genetic networks can evolve the capacity to sense internal mutations. GRNs are modelled using Threshold Boolean Networks (TBNs), abstracting away details so that certain computational approaches become viable. For instance in Chapter 3, all possible TBNs that attain a specified functional path (of the form {v(t)}T t=0) through the expression state space are exhaustively found from a possible 3N2 TBNs, where N is the number of genes (nodes) in the network. This allows for the detailed examination of the complete neutral evolutionary space of a given functional path. It is demonstrated that the major quantities of interest, such as the connectivity of this neutral space under point mutations, the mutational and noise robustness of the TBNs in this space and even the number of networks all depend strongly on the duration T of the paths. The neutral space is found to disintegrate rapidly into disconnected components as T is increased. The effect of more exotic functional forms is also investigated. Chapter 4 focuses on evolving networks which are sensitive to deletion mutations. It is found that increased sensitivity is readily evolvable in TBNs, with the networks evolving to be more topologically balanced (they possess a similar number of excitatory and inhibitory interactions). Networks are only found to achieve maximal sensitivity through attaining long limit cycles. The study of sensitivity is extended to static populations of TBNs in Chapter 5 and the question is asked about whether a population of cells can develop the capacity to sense the presence of a mutant among them. The multicellular framework is also used to investigate the effect of intercellular connectivity on the dynamics. It is found that the greater the intercellular connectivity the more uniform the expression patterns are between cells. Chapter 6 applies the general Ergodic Set (ES) [Ribeiro and Kauffman, 2007] concept to stem cell differentiation and cancer. An alternative hypothesis to that in [Serra et al., 2010] is proposed about how to model stem cell differentiation using ESs. Coupled with results from Chapter 4, I suggest that under this new hypothesis, pluripotent stem cells will correspond to more sensitive TBNs and that differentiated stem cells will correspond to more robust TBNs

Bacterial drug tolerance under clinical conditions is governed by anaerobic adaptation but not anaerobic respiration

Non-inherited antibiotic resistance is a phenomenon whereby a subpopulation of genetically identical bacteria displays phenotypic tolerance to antibiotics. We show that, compared to E. coli, the clinically relevant genus Burkholderia displays much higher levels of cells that tolerate ceftazidime. Measuring the dynamics of the formation of drug-tolerant cells, under conditions that mimic in vivo infections, we show that in Burkholderia oxygen levels affect the formation of these cells. The drug-tolerant cells are characterised by an anaerobic metabolic signature and can be eliminated by oxygenation of the system, or by the addition of nitrate. The transcriptome profile suggest that these cells are not dormant persister cells, and are likely to be drug-tolerant as a consequence of the up-regulation of anaerobic nitrate respiration, efflux pumps, beta lactamases and stress response proteins. These findings have important implications for the treatment of chronic bacterial infections and the methodologies and conditions which are used to study drug-tolerant and persister cells in vitro