How duplicated transcription regulators can diversify to govern the expression of nonoverlapping sets of genes.

The duplication of transcription regulators can elicit major regulatory network rearrangements over evolutionary timescales. However, few examples of duplications resulting in gene network expansions are understood in molecular detail. Here we show that four Candida albicans transcription regulators that arose by successive duplications have differentiated from one another by acquiring different intrinsic DNA-binding specificities, different preferences for half-site spacing, and different associations with cofactors. The combination of these three mechanisms resulted in each of the four regulators controlling a distinct set of target genes, which likely contributed to the adaption of this fungus to its human host. Our results illustrate how successive duplications and diversification of an ancestral transcription regulator can underlie major changes in an organism's regulatory circuitry.

A framework for integrating directed and undirected annotations to build explanatory models of cis-eQTL data.

A longstanding goal of regulatory genetics is to understand how variants in genome sequences lead to changes in gene expression. Here we present a method named Bayesian Annotation Guided eQTL Analysis (BAGEA), a variational Bayes framework to model cis-eQTLs using directed and undirected genomic annotations. We used BAGEA to integrate directed genomic annotations with eQTL summary statistics from tissues of various origins. This analysis revealed epigenetic marks that are relevant for gene expression in different tissues and cell types. We estimated the predictive power of the models that were fitted based on directed genomic annotations. This analysis showed that, depending on the underlying eQTL data used, the directed genomic annotations could predict up to 1.5% of the variance observed in the expression of genes with top nominal eQTL association p-values < 10-7. For genes with estimated effect sizes in the top 25% quantile, up to 5% of the expression variance could be predicted. Based on our results, we recommend the use of BAGEA for the analysis of cis-eQTL data to reveal annotations relevant to expression biology.

Evolution of increased complexity in a molecular machine.

Many cellular processes are carried out by molecular 'machines'-assemblies of multiple differentiated proteins that physically interact to execute biological functions. Despite much speculation, strong evidence of the mechanisms by which these assemblies evolved is lacking. Here we use ancestral gene resurrection and manipulative genetic experiments to determine how the complexity of an essential molecular machine--the hexameric transmembrane ring of the eukaryotic V-ATPase proton pump--increased hundreds of millions of years ago. We show that the ring of Fungi, which is composed of three paralogous proteins, evolved from a more ancient two-paralogue complex because of a gene duplication that was followed by loss in each daughter copy of specific interfaces by which it interacts with other ring proteins. These losses were complementary, so both copies became obligate components with restricted spatial roles in the complex. Reintroducing a single historical mutation from each paralogue lineage into the resurrected ancestral proteins is sufficient to recapitulate their asymmetric degeneration and trigger the requirement for the more elaborate three-component ring. Our experiments show that increased complexity in an essential molecular machine evolved because of simple, high-probability evolutionary processes, without the apparent evolution of novel functions. They point to a plausible mechanism for the evolution of complexity in other multi-paralogue protein complexes.

PhyloBot: A Web Portal for Automated Phylogenetics, Ancestral Sequence Reconstruction, and Exploration of Mutational Trajectories.

The method of phylogenetic ancestral sequence reconstruction is a powerful approach for studying evolutionary relationships among protein sequence, structure, and function. In particular, this approach allows investigators to (1) reconstruct and "resurrect" (that is, synthesize in vivo or in vitro) extinct proteins to study how they differ from modern proteins, (2) identify key amino acid changes that, over evolutionary timescales, have altered the function of the protein, and (3) order historical events in the evolution of protein function. Widespread use of this approach has been slow among molecular biologists, in part because the methods require significant computational expertise. Here we present PhyloBot, a web-based software tool that makes ancestral sequence reconstruction easy. Designed for non-experts, it integrates all the necessary software into a single user interface. Additionally, PhyloBot provides interactive tools to explore evolutionary trajectories between ancestors, enabling the rapid generation of hypotheses that can be tested using genetic or biochemical approaches. Early versions of this software were used in previous studies to discover genetic mechanisms underlying the functions of diverse protein families, including V-ATPase ion pumps, DNA-binding transcription regulators, and serine/threonine protein kinases. PhyloBot runs in a web browser, and is available at the following URL: http://www.phylobot.com. The software is implemented in Python using the Django web framework, and runs on elastic cloud computing resources from Amazon Web Services. Users can create and submit jobs on our free server (at the URL listed above), or use our open-source code to launch their own PhyloBot server.

Identification and characterization of a previously undescribed family of sequence-specific DNA-binding domains.

Sequence-specific DNA-binding proteins are among the most important classes of gene regulatory proteins, controlling changes in transcription that underlie many aspects of biology. In this work, we identify a transcriptional regulator from the human fungal pathogen Candida albicans that binds DNA specifically but has no detectable homology with any previously described DNA- or RNA-binding protein. This protein, named White-Opaque Regulator 3 (Wor3), regulates white-opaque switching, the ability of C. albicans to switch between two heritable cell types. We demonstrate that ectopic overexpression of WOR3 results in mass conversion of white cells to opaque cells and that deletion of WOR3 affects the stability of opaque cells at physiological temperatures. Genome-wide chromatin immunoprecipitation of Wor3 and gene expression profiling of a wor3 deletion mutant strain indicate that Wor3 is highly integrated into the previously described circuit regulating white-opaque switching and that it controls a subset of the opaque transcriptional program. We show by biochemical, genetic, and microfluidic experiments that Wor3 binds directly to DNA in a sequence-specific manner, and we identify the set of cis-regulatory sequences recognized by Wor3. Bioinformatic analyses indicate that the Wor3 family arose more recently in evolutionary time than most previously described DNA-binding domains; it is restricted to a small number of fungi that include the major fungal pathogens of humans. These observations show that new families of sequence-specific DNA-binding proteins may be restricted to small clades and suggest that current annotations--which rely on deep conservation--underestimate the fraction of genes coding for transcriptional regulators.

Patch and matrix striatonigral neurons differentially regulate locomotion.

Striatonigral neurons, known to promote locomotion, reside in both the patch and matrix compartments of the dorsal striatum. However, their compartment-specific contributions to locomotion remain largely unexplored. Using molecular identifier Kremen1 and Calb1 , we showed in mouse models that patch and matrix striatonigral neurons exert opposite influences on locomotion. Matrix striatonigral neurons reduced their activity before the cessation of self-paced locomotion, while patch striatonigral neuronal activity increased, suggesting an inhibitory function. Indeed, optogenetic activation of patch striatonigral neurons suppressed ongoing locomotion with reduced striatal dopamine release, contrasting with the locomotion-promoting effect of matrix striatonigral neurons, which showed an initial increase in dopamine release. Furthermore, genetic deletion of the GABA-B receptor in Aldehyde dehydrogenase 1A1-positive (ALDH1A1 + ) nigrostriatal dopaminergic neurons completely abolished the locomotion-suppressing effect of patch striatonigral neurons. Our findings unravel a compartment-specific mechanism governing locomotion in the dorsal striatum, where patch striatonigral neurons suppress locomotion by inhibiting ALDH1A1 + nigrostriatal dopaminergic neurons.

Patch and matrix striatonigral neurons differentially regulate locomotion.

The striatonigral neurons are known to promote locomotion1,2. These neurons reside in both the patch (also known as striosome) and matrix compartments of the dorsal striatum3-5. However, the specific contribution of patch and matrix striatonigral neurons to locomotion remain largely unexplored. Using molecular identifier Kringle-Containing Protein Marking the Eye and the Nose (Kremen1) and Calbidin (Calb1)6, we showed in mouse models that patch and matrix striatonigral neurons exert opposite influence on locomotion. While a reduction in neuronal activity in matrix striatonigral neurons precedes the cessation of locomotion, fiber photometry recording during self-paced movement revealed an unexpected increase of patch striatonigral neuron activity, indicating an inhibitory function. Indeed, optogenetic activation of patch striatonigral neurons suppressed locomotion, contrasting with the locomotion-promoting effect of matrix striatonigral neurons. Consistently, patch striatonigral neuron activation markedly inhibited dopamine release, whereas matrix striatonigral neuron activation initially promoted dopamine release. Moreover, the genetic deletion of inhibitory GABA-B receptor Gabbr1 in Aldehyde dehydrogenase 1A1-positive (ALDH1A1+) nigrostriatal dopaminergic neurons (DANs) completely abolished the locomotion-suppressing effect caused by activating patch striatonigral neurons. Together, our findings unravel a compartment-specific mechanism governing locomotion in the dorsal striatum, where patch striatonigral neurons suppress locomotion by inhibiting the activity of ALDH1A1+ nigrostriatal DANs.

Unravelling cell type-specific responses to Parkinson's Disease at single cell resolution.

Parkinson's Disease (PD) is the second most common neurodegenerative disorder. The pathological hallmark of PD is loss of dopaminergic neurons and the presence of aggregated α-synuclein, primarily in the substantia nigra pars compacta (SNpc) of the midbrain. However, the molecular mechanisms that underlie the pathology in different cell types is not currently understood. Here, we present a single nucleus transcriptome analysis of human post-mortem SNpc obtained from 15 sporadic Parkinson's Disease (PD) cases and 14 Controls. Our dataset comprises ∼84K nuclei, representing all major cell types of the brain, allowing us to obtain a transcriptome-level characterization of these cell types. Importantly, we identify multiple subpopulations for each cell type and describe specific gene sets that provide insights into the differing roles of these subpopulations. Our findings reveal a significant decrease in neuronal cells in PD samples, accompanied by an increase in glial cells and T cells. Subpopulation analyses demonstrate a significant depletion of tyrosine hydroxylase (TH) enriched astrocyte, microglia and oligodendrocyte populations in PD samples, as well as TH enriched neurons, which are also depleted. Moreover, marker gene analysis of the depleted subpopulations identified 28 overlapping genes, including those associated with dopamine metabolism (e.g., ALDH1A1, SLC6A3 & SLC18A2). Overall, our study provides a valuable resource for understanding the molecular mechanisms involved in dopaminergic neuron degeneration and glial responses in PD, highlighting the existence of novel subpopulations and cell type-specific gene sets.

Dilemmas in parenteral glucose delivery and approach to glucose monitoring and interpretation in the neonate.

Glucose control continues to be challenging for intensivists, in particular in high-risk neonates. Many factors play a role in glucose regulation including intrinsic and extrinsic factors. Optimal targets for euglycemia are debatable with uncertain short and long-term effects. Glucose measurement technology has continued to advance over the past decade; unfortunately, the availability of these advanced devices outside of research continues to be problematic. Treatment approaches should be individualized depending on etiology, symptoms, and neonatal conditions. Glucose infusions should be titrated based upon variations in organ glucose uptake, co-morbidities and postnatal development. In this article we summarize the most common dilemmas encountered in the NICU: ranges for euglycemia, physiological differences, approach for glucose measurements, monitoring and best strategies to control parenteral glucose delivery.

Systematic comparative analysis of strand-specific RNA-seq library preparation methods for low input samples.

Despite the recent precipitous decline in the cost of genome sequencing, library preparation for RNA-seq is still laborious and expensive for applications such as high throughput screening. Limited availability of RNA generated by some experimental workflows poses an additional challenge and increases the cost of RNA library preparation. In a search for low cost, automation-compatible RNA library preparation kits that maintain strand specificity and are amenable to low input RNA quantities, we systematically tested two recent commercial technologies-Swift RNA and Swift Rapid RNA, presently offered by Integrated DNA Technologies (IDT) -alongside the Illumina TruSeq stranded mRNA, the de facto standard workflow for bulk transcriptomics. We used the Universal Human Reference RNA (UHRR) (composed of equal quantities of total RNA from 10 human cancer cell lines) to benchmark gene expression in these kits, at input quantities ranging between 10 to 500 ng. We found normalized read counts between all treatment groups to be in high agreement. Compared to the Illumina TruSeq stranded mRNA kit, both Swift RNA library kits offer shorter workflow times enabled by their patented Adaptase technology. We also found the Swift RNA kit to produce the fewest number of differentially expressed genes and pathways directly attributable to input mRNA amount.

Robustness of ancestral sequence reconstruction to phylogenetic uncertainty.

Ancestral sequence reconstruction (ASR) is widely used to formulate and test hypotheses about the sequences, functions, and structures of ancient genes. Ancestral sequences are usually inferred from an alignment of extant sequences using a maximum likelihood (ML) phylogenetic algorithm, which calculates the most likely ancestral sequence assuming a probabilistic model of sequence evolution and a specific phylogeny--typically the tree with the ML. The true phylogeny is seldom known with certainty, however. ML methods ignore this uncertainty, whereas Bayesian methods incorporate it by integrating the likelihood of each ancestral state over a distribution of possible trees. It is not known whether Bayesian approaches to phylogenetic uncertainty improve the accuracy of inferred ancestral sequences. Here, we use simulation-based experiments under both simplified and empirically derived conditions to compare the accuracy of ASR carried out using ML and Bayesian approaches. We show that incorporating phylogenetic uncertainty by integrating over topologies very rarely changes the inferred ancestral state and does not improve the accuracy of the reconstructed ancestral sequence. Ancestral state reconstructions are robust to uncertainty about the underlying tree because the conditions that produce phylogenetic uncertainty also make the ancestral state identical across plausible trees; conversely, the conditions under which different phylogenies yield different inferred ancestral states produce little or no ambiguity about the true phylogeny. Our results suggest that ML can produce accurate ASRs, even in the face of phylogenetic uncertainty. Using Bayesian integration to incorporate this uncertainty is neither necessary nor beneficial.

Reconstructed evolutionary history of the yeast septins Cdc11 and Shs1.

Septins are GTP-binding proteins conserved across metazoans. They can polymerize into extended filaments and, hence, are considered a component of the cytoskeleton. The number of individual septins varies across the tree of life-yeast (Saccharomyces cerevisiae) has seven distinct subunits, a nematode (Caenorhabditis elegans) has two, and humans have 13. However, the overall geometric unit (an apolar hetero-octameric protomer and filaments assembled there from) has been conserved. To understand septin evolutionary variation, we focused on a related pair of yeast subunits (Cdc11 and Shs1) that appear to have arisen from gene duplication within the fungal clade. Either Cdc11 or Shs1 occupies the terminal position within a hetero-octamer, yet Cdc11 is essential for septin function and cell viability, whereas Shs1 is not. To discern the molecular basis of this divergence, we utilized ancestral gene reconstruction to predict, synthesize, and experimentally examine the most recent common ancestor ("Anc.11-S") of Cdc11 and Shs1. Anc.11-S was able to occupy the terminal position within an octamer, just like the modern subunits. Although Anc.11-S supplied many of the known functions of Cdc11, it was unable to replace the distinct function(s) of Shs1. To further evaluate the history of Shs1, additional intermediates along a proposed trajectory from Anc.11-S to yeast Shs1 were generated and tested. We demonstrate that multiple events contributed to the current properties of Shs1: (1) loss of Shs1-Shs1 self-association early after duplication, (2) co-evolution of heterotypic Cdc11-Shs1 interaction between neighboring hetero-octamers, and (3) eventual repurposing and acquisition of novel function(s) for its C-terminal extension domain. Thus, a pair of duplicated proteins, despite constraints imposed by assembly into a highly conserved multi-subunit structure, could evolve new functionality via a complex evolutionary pathway.

Characterization of recombinant fructose-1,6-bisphosphatase gene mutations: evidence of inhibition/activation of FBPase protein by gene mutation.

Specific residues of the highly regulated fructose-1,6-bisphosphatase (FBPase) enzyme serve as important contributors to the catalytic activity of the enzyme. Previous clinical studies exploring the genetic basis of hypoglycemia revealed two significant mutations in the coding region of the FBPase gene in patients with hypoglycemia, linking the AMP-binding site to the active site of the enzyme. In the present study, a full kinetic analysis of similar mutants was performed. Kinetic results of mutants Y164A and M177A revealed an approximate two to three-fold decrease in inhibitory constants (Ki's) for natural inhibitors AMP and fructose-2,6-bisphosphate (F2,6-BP) compared with the Wild-type enzyme (WT). A separate mutation (M248D) was performed in the active site of the enzyme to investigate whether the enzyme could be activated. This mutant displayed an approximate seven-fold increase in Ki for F2,6-BP. Interfacial mutants L56A and L73A exhibited an increase in Ki for F2,6-BP by approximately five-fold. Mutations in the AMP-binding site (K112A and Y113A) demonstrated an eight to nine-fold decrease in AMP inhibition. Additionally, mutant M248D displayed a four-fold decrease in its apparent Michelis constant (Km), and a six-fold increase in catalytic efficiency (CE). The importance-and medical relevance-of specific residues for FBPase structural/functional relationships in both the catalytic site and AMP-binding site is discussed.

Haploinsufficiency leads to neurodegeneration in C9ORF72 ALS/FTD human induced motor neurons.

An intronic GGGGCC repeat expansion in C9ORF72 is the most common cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), but the pathogenic mechanism of this repeat remains unclear. Using human induced motor neurons (iMNs), we found that repeat-expanded C9ORF72 was haploinsufficient in ALS. We found that C9ORF72 interacted with endosomes and was required for normal vesicle trafficking and lysosomal biogenesis in motor neurons. Repeat expansion reduced C9ORF72 expression, triggering neurodegeneration through two mechanisms: accumulation of glutamate receptors, leading to excitotoxicity, and impaired clearance of neurotoxic dipeptide repeat proteins derived from the repeat expansion. Thus, cooperativity between gain- and loss-of-function mechanisms led to neurodegeneration. Restoring C9ORF72 levels or augmenting its function with constitutively active RAB5 or chemical modulators of RAB5 effectors rescued patient neuron survival and ameliorated neurodegenerative processes in both gain- and loss-of-function C9ORF72 mouse models. Thus, modulating vesicle trafficking was able to rescue neurodegeneration caused by the C9ORF72 repeat expansion. Coupled with rare mutations in ALS2, FIG4, CHMP2B, OPTN and SQSTM1, our results reveal mechanistic convergence on vesicle trafficking in ALS and FTD.

Evolution of an ancient protein function involved in organized multicellularity in animals.

To form and maintain organized tissues, multicellular organisms orient their mitotic spindles relative to neighboring cells. A molecular complex scaffolded by the GK protein-interaction domain (GKPID) mediates spindle orientation in diverse animal taxa by linking microtubule motor proteins to a marker protein on the cell cortex localized by external cues. Here we illuminate how this complex evolved and commandeered control of spindle orientation from a more ancient mechanism. The complex was assembled through a series of molecular exploitation events, one of which - the evolution of GKPID's capacity to bind the cortical marker protein - can be recapitulated by reintroducing a single historical substitution into the reconstructed ancestral GKPID. This change revealed and repurposed an ancient molecular surface that previously had a radically different function. We show how the physical simplicity of this binding interface enabled the evolution of a new protein function now essential to the biological complexity of many animals.

Building Data-Driven Pathways From Routinely Collected Hospital Data: A Case Study on Prostate Cancer.

BACKGROUND: Routinely collected data in hospitals is complex, typically heterogeneous, and scattered across multiple Hospital Information Systems (HIS). This big data, created as a byproduct of health care activities, has the potential to provide a better understanding of diseases, unearth hidden patterns, and improve services and cost. The extent and uses of such data rely on its quality, which is not consistently checked, nor fully understood. Nevertheless, using routine data for the construction of data-driven clinical pathways, describing processes and trends, is a key topic receiving increasing attention in the literature. Traditional algorithms do not cope well with unstructured processes or data, and do not produce clinically meaningful visualizations. Supporting systems that provide additional information, context, and quality assurance inspection are needed. OBJECTIVE: The objective of the study is to explore how routine hospital data can be used to develop data-driven pathways that describe the journeys that patients take through care, and their potential uses in biomedical research; it proposes a framework for the construction, quality assessment, and visualization of patient pathways for clinical studies and decision support using a case study on prostate cancer. METHODS: Data pertaining to prostate cancer patients were extracted from a large UK hospital from eight different HIS, validated, and complemented with information from the local cancer registry. Data-driven pathways were built for each of the 1904 patients and an expert knowledge base, containing rules on the prostate cancer biomarker, was used to assess the completeness and utility of the pathways for a specific clinical study. Software components were built to provide meaningful visualizations for the constructed pathways. RESULTS: The proposed framework and pathway formalism enable the summarization, visualization, and querying of complex patient-centric clinical information, as well as the computation of quality indicators and dimensions. A novel graphical representation of the pathways allows the synthesis of such information. CONCLUSIONS: Clinical pathways built from routinely collected hospital data can unearth information about patients and diseases that may otherwise be unavailable or overlooked in hospitals. Data-driven clinical pathways allow for heterogeneous data (ie, semistructured and unstructured data) to be collated over a unified data model and for data quality dimensions to be assessed. This work has enabled further research on prostate cancer and its biomarkers, and on the development and application of methods to mine, compare, analyze, and visualize pathways constructed from routine data. This is an important development for the reuse of big data in hospitals.

Mapping of neurokinin-like immunoreactivity in the human brainstem.

BACKGROUND: Using an indirect immunoperoxidase technique, we have studied the distribution of immunoreactive fibers and cell bodies containing neurokinin in the adult human brainstem with no prior history of neurological or psychiatric disease. RESULTS: Clusters of immunoreactive cell bodies and high densities of neurokinin-immunoreactive fibers were located in the periaqueductal gray, the dorsal motor nucleus of the vagus and in the reticular formation of the medulla, pons and mesencephalon. Moreover, immunoreactive cell bodies were found in the inferior colliculus, the raphe obscurus, the nucleus prepositus hypoglossi, and in the midline of the anterior medulla oblongata. In general, immunoreactive fibers containing neurokinin were observed throughout the whole brainstem. In addition to the nuclei mentioned above, the highest densities of such immunoreactive fibers were located in the spinal trigeminal nucleus, the lateral reticular nucleus, the nucleus of the solitary tract, the superior colliculus, the substantia nigra, the nucleus ambiguus, the gracile nucleus, the cuneate nucleus, the motor hypoglossal nucleus, the medial and superior vestibular nuclei, the nucleus prepositus hypoglossi and the interpeduncular nucleus. CONCLUSION: The widespread distribution of immunoreactive structures containing neurokinin in the human brainstem indicates that neurokinin might be involved in several physiological mechanisms, acting as a neurotransmitter and/or neuromodulator.

General practice training environment and its impact on preparedness.

INTRODUCTION: The notion of preparedness for practice is poorly defined in medical education literature. It is unclear what preparedness means and how the training environment impacts on preparedness for practice. OBJECTIVES: This paper aims to explore the meaning that GP trainees and newly qualified GPs attach to the notion of preparedness, and to examine the ways in which they perceive their training environment to impact on preparedness. METHODS: We used a qualitative interpretive approach and conducted 27 in-depth semi-structured interviews with 15 newly qualified GPs and 12 GP trainees at the end of their training. RESULTS: Two central categories describing preparedness emerged; 'confidence' and 'adaptability'. Inclusive training practices, characterised by non-hierarchical relationships between the doctors, particularly vis-à-vis trainees, were reported to be more 'progressive' and were better at preparing trainees. The way the training practice can impact on preparedness can be explained drawing on Lave and Wenger's theory of 'situated learning'. The role of the trainer was also pivotal in preparing trainees. Supervision tailored to trainees' needs, and guided decision making enhanced confidence of trainees in their ability to work independently in the future. CONCLUSIONS: We suggest that for GP trainees to be better prepared it is not enough to extend GP training; rather it is important that GP trainees' time is spent in inclusive training environments.

Ancestral resurrection reveals evolutionary mechanisms of kinase plasticity.

Protein kinases have evolved diverse specificities to enable cellular information processing. To gain insight into the mechanisms underlying kinase diversification, we studied the CMGC protein kinases using ancestral reconstruction. Within this group, the cyclin dependent kinases (CDKs) and mitogen activated protein kinases (MAPKs) require proline at the +1 position of their substrates, while Ime2 prefers arginine. The resurrected common ancestor of CDKs, MAPKs, and Ime2 could phosphorylate substrates with +1 proline or arginine, with preference for proline. This specificity changed to a strong preference for +1 arginine in the lineage leading to Ime2 via an intermediate with equal specificity for proline and arginine. Mutant analysis revealed that a variable residue within the kinase catalytic cleft, DFGx, modulates +1 specificity. Expansion of Ime2 kinase specificity by mutation of this residue did not cause dominant deleterious effects in vivo. Tolerance of cells to new specificities likely enabled the evolutionary divergence of kinases.