Structural and functional comparison of magnesium transporters throughout evolution.

Magnesium (Mg2+) is the most prevalent divalent intracellular cation. As co-factor in many enzymatic reactions, Mg2+ is essential for protein synthesis, energy production, and DNA stability. Disturbances in intracellular Mg2+ concentrations, therefore, unequivocally result in delayed cell growth and metabolic defects. To maintain physiological Mg2+ levels, all organisms rely on balanced Mg2+ influx and efflux via Mg2+ channels and transporters. This review compares the structure and the function of prokaryotic Mg2+ transporters and their eukaryotic counterparts. In prokaryotes, cellular Mg2+ homeostasis is orchestrated via the CorA, MgtA/B, MgtE, and CorB/C Mg2+ transporters. For CorA, MgtE, and CorB/C, the motifs that form the selectivity pore are conserved during evolution. These findings suggest that CNNM proteins, the vertebrate orthologues of CorB/C, also have Mg2+ transport capacity. Whereas CorA and CorB/C proteins share the gross quaternary structure and functional properties with their respective orthologues, the MgtE channel only shares the selectivity pore with SLC41 Na+/Mg2+ transporters. In eukaryotes, TRPM6 and TRPM7 Mg2+ channels provide an additional Mg2+ transport mechanism, consisting of a fusion of channel with a kinase. The unique features these TRP channels allow the integration of hormonal, cellular, and transcriptional regulatory pathways that determine their Mg2+ transport capacity. Our review demonstrates that understanding the structure and function of prokaryotic magnesiotropic proteins aids in our basic understanding of Mg2+ transport.

Does mitochondrial DNA evolution in metazoa drive the origin of new mitochondrial proteins?

Most eukaryotic cells contain mitochondria with a genome that evolved from their α-proteobacterial ancestor. In the course of eukaryotic evolution, the mitochondrial genome underwent a dramatic reduction in size, caused by the loss and translocation of genes. This required adjustments in mitochondrial gene expression mechanisms and resulted in a complex collaborative system of mitochondrially encoded transfer RNAs and ribosomal RNAs with nuclear encoded proteins to express the mitochondrial encoded oxidative phosphorylation (OXPHOS) proteins. In this review, we examine mitochondrial gene expression from an evolutionary point of view: to what extent can we correlate changes in the mitochondrial genome in the evolutionary lineage leading to human with the origin of new nuclear encoded proteins. We dated the evolutionary origin of mitochondrial proteins that interact with mitochondrial DNA or its RNA and/or protein products in a systematic manner and compared them with documented changes in the mitochondrial DNA. We find anecdotal but accumulating evidence that metazoan RNA-interacting proteins arose in conjunction with changes of the mitochondrial DNA. We find no substantial evidence for such compensatory evolution in new OXPHOS proteins, which appear to be constrained by the ability to form supercomplexes. © 2018 IUBMB Life, 70(12):1240-1250, 2018.

Genome phylogeny based on gene content.

Species phylogenies derived from comparisons of single genes are rarely consistent with each other, due to horizontal gene transfer, unrecognized paralogy and highly variable rates of evolution. The advent of completely sequenced genomes allows the construction of a phylogeny that is less sensitive to such inconsistencies and more representative of whole-genomes than are single-gene trees. Here, we present a distance-based phylogeny constructed on the basis of gene content, rather than on sequence identity, of 13 completely sequenced genomes of unicellular species. The similarity between two species is defined as the number of genes that they have in common divided by their total number of genes. In this type of phylogenetic analysis, evolutionary distance can be interpreted in terms of evolutionary events such as the acquisition and loss of genes, whereas the underlying properties (the gene content) can be interpreted in terms of function. As such, it takes a position intermediate to phylogenies based on single genes and phylogenies based on phenotypic characteristics. Although our comprehensive genome phylogeny is independent of phylogenies based on the level of sequence identity of individual genes, it correlates with the standard reference of prokarytic phylogeny based on sequence similarity of 16s rRNA. Thus, shared gene content between genomes is quantitatively determined by phylogeny, rather than by phenotype, and horizontal gene transfer has only a limited role in determining the gene content of genomes.

Chemokine profiling in children and adults with symptomatic and asymptomatic respiratory viral infections.

Molecular diagnosis; Viral infection; Chemokines; Disease prognosis; CXCL10; CXCL11; CCL3; CCL4; CCL5; Random forest.

Surveillance-embedded genomic outbreak resolution of methicillin-susceptible Staphylococcus aureus in a neonatal intensive care unit.

We observed an increase in methicillin-susceptible Staphylococcus aureus (MSSA) infections at a Dutch neonatal intensive care unit. Weekly neonatal MSSA carriage surveillance and cross-sectional screenings of health care workers (HCWs) were available for outbreak tracing. Traditional clustering of MSSA isolates by spa typing and Multiple-Locus Variable number tandem repeat Analysis (MLVA) suggested that nosocomial transmission had contributed to the infections. We investigated whether whole-genome sequencing (WGS) of MSSA surveillance would provide additional evidence for transmission. MSSA isolates from neonatal infections, carriage surveillance, and HCWs were subjected to WGS and bioinformatic analysis for identification and localization of high-quality single nucleotide polymorphisms, and in-depth analysis of subsets of isolates. By measuring the genetic diversity in background surveillance, we defined transmission-level relatedness and identified isolates that had been unjustly assigned to clusters based on MLVA, while spa typing was concordant but of insufficient resolution. Detailing particular subsets of isolates provided evidence that HCWs were involved in multiple outbreaks, yet it alleviated concerns about one particular HCW. The improved resolution and accuracy of genomic outbreak analyses substantially altered the view on outbreaks, along with apposite measures. Therefore, inclusion of the circulating background population has the potential to overcome current issues in genomic outbreak inference.

Conservation of gene order: a fingerprint of proteins that physically interact.

A systematic comparison of nine bacterial and archaeal genomes reveals a low level of gene-order (and operon architecture) conservation. Nevertheless, a number of gene pairs are conserved. The proteins encoded by conserved gene pairs appear to interact physically. This observation can therefore be used to predict functions of, and interactions between, prokaryotic gene products.

Competencies on environmental health and pedagogical approaches in the nursing curriculum: A systematic review of the literature.

It has been suggested that climate change is the biggest threat to public health for the 21st Century; increased demand on health services will impact on already overstretched resources and systems will need to be able to respond. However limited attention is given to climate change and sustainability in nursing education; there is no clear guidance on curricula content for nurses or recommendations regarding the skills and competencies that will be required. Literature published in Dutch, English, German, and Spanish was searched and 32 papers met the inclusion criteria for the review. Results suggests that holistic/systems thinking is relevant to healthcare so bringing a 'sustainability lens' to nursing curricula could be seen as being consistent with wider determinants of health. The literature review has identified the educational approaches necessary to provide a broad based curriculum and a cross-disciplinary approach. The findings suggest that topics such as the use of resources, food, health promotion, globalism, disease management, and the environmental impact of delivering healthcare, if embedded in nursing education could support the nursing profession's response for this new and important aspect of healthcare.

Assessment of phylogenomic and orthology approaches for phylogenetic inference.

MOTIVATION: Phylogenomics integrates the vast amount of phylogenetic information contained in complete genome sequences, and is rapidly becoming the standard for reliably inferring species phylogenies. There are, however, fundamental differences between the ways in which phylogenomic approaches like gene content, superalignment, superdistance and supertree integrate the phylogenetic information from separate orthologous groups. Furthermore, they all depend on the method by which the orthologous groups are initially determined. Here, we systematically compare these four phylogenomic approaches, in parallel with three approaches for large-scale orthology determination: pairwise orthology, cluster orthology and tree-based orthology. RESULTS: Including various phylogenetic methods, we apply a total of 54 fully automated phylogenomic procedures to the fungi, the eukaryotic clade with the largest number of sequenced genomes, for which we retrieved a golden standard phylogeny from the literature. Phylogenomic trees based on gene content show, relative to the other methods, a bias in the tree topology that parallels convergence in lifestyle among the species compared, indicating convergence in gene content. CONCLUSIONS: Complete genomes are no guarantee for good or even consistent phylogenies. However, the large amounts of data in genomes enable us to carefully select the data most suitable for phylogenomic inference. In terms of performance, the superalignment approach, combined with restrictive orthology, is the most successful in recovering a fungal phylogeny that agrees with current taxonomic views, and allows us to obtain a high-resolution phylogeny. We provide solid support for what has grown to be a common practice in phylogenomics during its advance in recent years. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Formation of peroxisomes: present and past.

Eukaryotic cells contain functionally distinct, membrane enclosed compartments called organelles. Here we like to address two questions concerning this architectural lay out. How did this membrane complexity arise during evolution and how is this collection of organelles maintained in multiplying cells to ensure that new cells retain a complete set of them. We will try to address these questions with peroxisomes as a focal point of interest.

Inversions and the dynamics of eukaryotic gene order.

Comparisons of the gene order in closely related genomes reveal a major role for inversions in the genome shuffling process. In contrast to prokaryotes, where the inversions are predominantly large, half of the inversions between Saccharomyces cerevisiae and Candida albicans appear to be small, often encompassing only a single gene. Overall the genome rearrangement rate appears higher in eukaryotes than in prokaryotes, and the current genome data do not indicate that functional constraints on the co-expression of neighboring genes have a large role in conserving eukaryotic gene order. Nevertheless, qualitatively interesting examples of conservation of gene order in eukaryotes can be observed.

Modularity in the gain and loss of genes: applications for function prediction.

Genes that are clustered on multiple genomes and are likely to functionally interact tend to be gained or lost together during genome evolution. Here, we demonstrate that exceptions to this pattern indicate relatively distant functional interactions between the encoded proteins. Hence, this can be used to divide predicted clusters of functionally interacting proteins into sub-clusters, and as such, to refine the prediction of their function and functional interactions.

Exploitation of gene context.

Recently, a number of techniques have been proposed that use completely sequenced genomes for the function prediction of individual proteins encoded therein. They use the fusion of genes, their conserved location in operons or merely their co-occurrence in genomes to predict the existence of functional interactions between the proteins they encode. This type of information complements functional features that are predicted by classical homology-based search techniques.

Predicting disease genes using protein-protein interactions.

BACKGROUND: The responsible genes have not yet been identified for many genetically mapped disease loci. Physically interacting proteins tend to be involved in the same cellular process, and mutations in their genes may lead to similar disease phenotypes. OBJECTIVE: To investigate whether protein-protein interactions can predict genes for genetically heterogeneous diseases. METHODS: 72,940 protein-protein interactions between 10,894 human proteins were used to search 432 loci for candidate disease genes representing 383 genetically heterogeneous hereditary diseases. For each disease, the protein interaction partners of its known causative genes were compared with the disease associated loci lacking identified causative genes. Interaction partners located within such loci were considered candidate disease gene predictions. Prediction accuracy was tested using a benchmark set of known disease genes. RESULTS: Almost 300 candidate disease gene predictions were made. Some of these have since been confirmed. On average, 10% or more are expected to be genuine disease genes, representing a 10-fold enrichment compared with positional information only. Examples of interesting candidates are AKAP6 for arrythmogenic right ventricular dysplasia 3 and SYN3 for familial partial epilepsy with variable foci. CONCLUSIONS: Exploiting protein-protein interactions can greatly increase the likelihood of finding positional candidate disease genes. When applied on a large scale they can lead to novel candidate gene predictions.

POMT2 mutations cause alpha-dystroglycan hypoglycosylation and Walker-Warburg syndrome.

BACKGROUND: Walker-Warburg syndrome (WWS) is an autosomal recessive condition characterised by congenital muscular dystrophy, structural brain defects, and eye malformations. Typical brain abnormalities are hydrocephalus, lissencephaly, agenesis of the corpus callosum, fusion of the hemispheres, cerebellar hypoplasia, and neuronal overmigration, which causes a cobblestone cortex. Ocular abnormalities include cataract, microphthalmia, buphthalmos, and Peters anomaly. WWS patients show defective O-glycosylation of alpha-dystroglycan (alpha-DG), which plays a key role in bridging the cytoskeleton of muscle and CNS cells with extracellular matrix proteins, important for muscle integrity and neuronal migration. In 20% of the WWS patients, hypoglycosylation results from mutations in either the protein O-mannosyltransferase 1 (POMT1), fukutin, or fukutin related protein (FKRP) genes. The other genes for this highly heterogeneous disorder remain to be identified. OBJECTIVE: To look for mutations in POMT2 as a cause of WWS, as both POMT1 and POMT2 are required to achieve protein O-mannosyltransferase activity. METHODS: A candidate gene approach combined with homozygosity mapping. RESULTS: Homozygosity was found for the POMT2 locus at 14q24.3 in four of 11 consanguineous WWS families. Homozygous POMT2 mutations were present in two of these families as well as in one patient from another cohort of six WWS families. Immunohistochemistry in muscle showed severely reduced levels of glycosylated alpha-DG, which is consistent with the postulated role for POMT2 in the O-mannosylation pathway. CONCLUSIONS: A fourth causative gene for WWS was uncovered. These genes account for approximately one third of the WWS cases. Several more genes are anticipated, which are likely to play a role in glycosylation of alpha-DG.

STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene.

The repeated occurrence of genes in each other's neighbourhood on genomes has been shown to indicate a functional association between the proteins they encode. Here we introduce STRING (search tool for recurring instances of neighbouring genes), a tool to retrieve and display the genes a query gene repeatedly occurs with in clusters on the genome. The tool performs iterative searches and visualises the results in their genomic context. By finding the genomically associated genes for a query, it delineates a set of potentially functionally associated genes. The usefulness of STRING is illustrated with an example that suggests a functional context for an RNA methylase with unknown specificity.

Re-annotating the Mycoplasma pneumoniae genome sequence: adding value, function and reading frames.

Four years after the original sequence submission, we have re-annotated the genome of Mycoplasma pneumoniae to incorporate novel data. The total number of ORFss has been increased from 677 to 688 (10 new proteins were predicted in intergenic regions, two further were newly identified by mass spectrometry and one protein ORF was dismissed) and the number of RNAs from 39 to 42 genes. For 19 of the now 35 tRNAs and for six other functional RNAs the exact genome positions were re-annotated and two new tRNA(Leu) and a small 200 nt RNA were identified. Sixteen protein reading frames were extended and eight shortened. For each ORF a consistent annotation vocabulary has been introduced. Annotation reasoning, annotation categories and comparisons to other published data on M.pneumoniae functional assignments are given. Experimental evidence includes 2-dimensional gel electrophoresis in combination with mass spectrometry as well as gene expression data from this study. Compared to the original annotation, we increased the number of proteins with predicted functional features from 349 to 458. The increase includes 36 new predictions and 73 protein assignments confirmed by the published literature. Furthermore, there are 23 reductions and 30 additions with respect to the previous annotation. mRNA expression data support transcription of 184 of the functionally unassigned reading frames.

Variation and evolution of the citric-acid cycle: a genomic perspective.

The presence of genes encoding enzymes involved in the citric-acid cycle has been studied in 19 completely sequenced genomes. In the majority of species, the cycle appears to be incomplete or absent. Several distinct, incomplete cycles reflect adaptations to different environments. Their distribution over the phylogenetic tree hints at precursors in the evolution of the citric-acid cycle.

Assessing the reliability of RNA folding using statistical mechanics.

We have analyzed the base-pairing probability distributions of 16 S and 16 S-like, and 23 S and 23 S-like ribosomal RNAs of Archaea, Bacteria, chloroplasts, mitochondria and Eukarya, as predicted by the partition function approach for RNA folding introduced by McCaskill. A quantitative analysis of the reliability of RNA folding is done by comparing the base-pairing probability distributions with the structures predicted by comparative sequence analysis (comparative structures). We distinguish two factors that show a relationship to the reliability of RNA minimum free energy structure. The first factor is the dominance of one particular base-pair or the absence of base-pairing for a given base within the base-pairing probability distribution (BPPD). We characterize the BPPD per base, including the probability of not base-pairing, by its Shannon entropy (S). The S value indicates the uncertainty about the base-pairing of a base: low S values result from BPPDs that are strongly dominated by a single base-pair or by the absence of base-pairing. We show that bases with low S values have a relatively high probability that their minimum free energy (MFE) structure corresponds to the comparative structure. The BPPDs of prokaryotes that live at high temperatures (thermophilic Archaea and Bacteria) have, calculated at 37 degrees C, lower S values than the BPPDs of prokaryotes that live at lower temperatures (mesophilic and psychrophilic Archaea and Bacteria). This reflects an adaptation of the ribosomal RNAs to the environmental temperature. A second factor that is important to consider with regard to the reliability of MFE structure folding is a variable degree of applicability of the thermodynamic model of RNA folding for different groups of RNAs. Here we show that among the bases that show low S values, the Archaea and Bacteria have similar, high probabilities (0.96 and 0.94 in 16 S and 0.93 and 0.91 in 23 S, respectively) that the MFE structure corresponds to the comparative structure. These probabilities are lower in the chloroplasts (16 S 0.91, 23 S 0.79), mitochondria (16 S-like 0.89, 23 S-like 0.69) and Eukarya (18 S 0.81, 28 S 0.86).

Predicting function: from genes to genomes and back.

Predicting function from sequence using computational tools is a highly complicated procedure that is generally done for each gene individually. This review focuses on the added value that is provided by completely sequenced genomes in function prediction. Various levels of sequence annotation and function prediction are discussed, ranging from genomic sequence to that of complex cellular processes. Protein function is currently best described in the context of molecular interactions. In the near future it will be possible to predict protein function in the context of higher order processes such as the regulation of gene expression, metabolic pathways and signalling cascades. The analysis of such higher levels of function description uses, besides the information from completely sequenced genomes, also the additional information from proteomics and expression data. The final goal will be to elucidate the mapping between genotype and phenotype.

Homology-based fold predictions for Mycoplasma genitalium proteins.

Homology search techniques based on the iterative PSI-BLAST method in combination with various filters for low sequence complexity are applied to assign folds to all Mycoplasma genitalium proteins. The resulting procedure (implemented as a web server) is able to predict at least one domain in 37% of these proteins automatically, with an estimated accuracy higher than 98%. Taking structural features such as coiled coil or transmembrane regions aside, folds can be assigned to more than half of the globular proteins in a bacterium just by iterative sequence comparison.