The Phosphatomes of the Multicellular Myxobacteria Myxococcus xanthus and Sorangium cellulosum in Comparison with Other Prokaryotic Genomes

Background

Analysis of the complete genomes from the multicellular myxobacteria Myxococcus xanthus and Sorangium cellulosum identified the highest number of eukaryotic-like protein kinases (ELKs) compared to all other genomes analyzed. High numbers of protein phosphatases (PPs) could therefore be anticipated, as reversible protein phosphorylation is a major regulation mechanism of fundamental biological processes.

Methodology

Here we report an intensive analysis of the phosphatomes of M. xanthus and S. cellulosum in which we constructed phylogenetic trees to position these sequences relative to PPs from other prokaryotic organisms.

Principal Findings

Predominant observations were: (i) M. xanthus and S. cellulosum possess predominantly Ser/Thr PPs; (ii) S. cellulosum encodes the highest number of PP2c-type phosphatases so far reported for a prokaryotic organism; (iii) in contrast to M. xanthus only S. cellulosum encodes high numbers of SpoIIE-like PPs; (iv) there is a significant lack of synteny among M. xanthus and S. cellulosum, and (v) the degree of co-organization between kinase and phosphatase genes is extremely low in these myxobacterial genomes.

Conclusions

We conclude that there has been a greater expansion of ELKs than PPs in multicellular myxobacteria.     

Population structure of mitochondrial genomes in Saccharomyces cerevisiae

Background

Rigorous study of mitochondrial functions and cell biology in the budding yeast, Saccharomyces cerevisiae has advanced our understanding of mitochondrial genetics. This yeast is now a powerful model for population genetics, owing to large genetic diversity and highly structured populations among wild isolates. Comparative mitochondrial genomic analyses between yeast species have revealed broad evolutionary changes in genome organization and architecture. A fine-scale view of recent evolutionary changes within S. cerevisiae has not been possible due to low numbers of complete mitochondrial sequences.

Results

To address challenges of sequencing AT-rich and repetitive mitochondrial DNAs (mtDNAs), we sequenced two divergent S. cerevisiae mtDNAs using a single-molecule sequencing platform (PacBio RS). Using de novo assemblies, we generated highly accurate complete mtDNA sequences. These mtDNA sequences were compared with 98 additional mtDNA sequences gathered from various published collections. Phylogenies based on mitochondrial coding sequences and intron profiles revealed that intraspecific diversity in mitochondrial genomes generally recapitulated the population structure of nuclear genomes. Analysis of intergenic sequence indicated a recent expansion of mobile elements in certain populations. Additionally, our analyses revealed that certain populations lacked introns previously believed conserved throughout the species, as well as the presence of introns never before reported in S. cerevisiae.

Conclusions

Our results revealed that the extensive variation in S. cerevisiae mtDNAs is often population specific, thus offering a window into the recent evolutionary processes shaping these genomes. In addition, we offer an effective strategy for sequencing these challenging AT-rich mitochondrial genomes for small scale projects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1664-4) contains supplementary material, which is available to authorized users.     

Genome-wide metabolic re-annotation of Ashbya gossypii: new insights into its metabolism through a comparative analysis with Saccharomyces cerevisiae and Kluyveromyces lactis

Background

Ashbya gossypii is an industrially relevant microorganism traditionally used for riboflavin production. Despite the high gene homology and gene order conservation comparatively with Saccharomyces cerevisiae, it presents a lower level of genomic complexity. Its type of growth, placing it among filamentous fungi, questions how close it really is from the budding yeast, namely in terms of metabolism, therefore raising the need for an extensive and thorough study of its entire metabolism. This work reports the first manual enzymatic genome-wide re-annotation of A. gossypii as well as the first annotation of membrane transport proteins.

Results

After applying a developed enzymatic re-annotation pipeline, 847 genes were assigned with metabolic functions. Comparatively to KEGG’s annotation, these data corrected the function for 14% of the common genes and increased the information for 52 genes, either completing existing partial EC numbers or adding new ones. Furthermore, 22 unreported enzymatic functions were found, corresponding to a significant increase in the knowledge of the metabolism of this organism. The information retrieved from the metabolic re-annotation and transport annotation was used for a comprehensive analysis of A. gossypii’s metabolism in comparison to the one of S. cerevisiae (post-WGD – whole genome duplication) and Kluyveromyces lactis (pre-WGD), suggesting some relevant differences in several parts of their metabolism, with the majority being found for the metabolism of purines, pyrimidines, nitrogen and lipids. A considerable number of enzymes were found exclusively in A. gossypii comparatively with K. lactis (90) and S. cerevisiae (13). In a similar way, 176 and 123 enzymatic functions were absent on A. gossypii comparatively to K. lactis and S. cerevisiae, respectively, confirming some of the well-known phenotypes of this organism.

Conclusions

This high quality metabolic re-annotation, together with the first membrane transporters annotation and the metabolic comparative analysis, represents a new important tool for the study and better understanding of A. gossypii’s metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-810) contains supplementary material, which is available to authorized users.     

A mathematical framework for modelling cambial surface evolution using a level set method

Background and Aims

During their lifetime, tree stems take a series of successive nested shapes. Individual tree growth models traditionally focus on apical growth and architecture. However, cambial growth, which is distributed over a surface layer wrapping the whole organism, equally contributes to plant form and function. This study aims at providing a framework to simulate how organism shape evolves as a result of a secondary growth process that occurs at the cellular scale.

Methods

The development of the vascular cambium is modelled as an expanding surface using the level set method. The surface consists of multiple compartments following distinct expansion rules. Growth behaviour can be formulated as a mathematical function of surface state variables and independent variables to describe biological processes.

Key Results

The model was coupled to an architectural model and to a forest stand model to simulate cambium dynamics and wood formation at the scale of the organism. The model is able to simulate competition between cambia, surface irregularities and local features. Predicting the shapes associated with arbitrarily complex growth functions does not add complexity to the numerical method itself.

Conclusions

Despite their slenderness, it is sometimes useful to conceive of trees as expanding surfaces. The proposed mathematical framework provides a way to integrate through time and space the biological and physical mechanisms underlying cambium activity. It can be used either to test growth hypotheses or to generate detailed maps of wood internal structure.     

Distance-based assessment of the localization of functional annotations in 3D genome reconstructions

Background

Recent studies used the contact data or three-dimensional (3D) genome reconstructions from Hi-C (chromosome conformation capture with next-generation sequencing) to assess the co-localization of functional genomic annotations in the nucleus. These analyses dichotomized data point pairs belonging to a functional annotation as “close” or “far” based on some threshold and then tested for enrichment of “close” pairs. We propose an alternative approach that avoids dichotomization of the data and instead directly estimates the significance of distances within the 3D reconstruction.

Results

We applied this approach to 3D genome reconstructions for Plasmodium falciparum, the causative agent of malaria, and Saccharomyces cerevisiae and compared the results to previous approaches. We found significant 3D co-localization of centromeres, telomeres, virulence genes, and several sets of genes with developmentally regulated expression in P. falciparum; and significant 3D co-localization of centromeres and long terminal repeats in S. cerevisiae. Additionally, we tested the experimental observation that telomeres form three to seven clusters in P. falciparum and S. cerevisiae. Applying affinity propagation clustering to telomere coordinates in the 3D reconstructions yielded six telomere clusters for both organisms.

Conclusions

Distance-based assessment replicated key findings, while avoiding dichotomization of the data (which previously yielded threshold-sensitive results).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-992) contains supplementary material, which is available to authorized users.     

Diversity arrays technology (DArT) for pan-genomic evolutionary studies of non-model organisms

Background

High-throughput tools for pan-genomic study, especially the DNA microarray platform, have sparked a remarkable increase in data production and enabled a shift in the scale at which biological investigation is possible. The use of microarrays to examine evolutionary relationships and processes, however, is predominantly restricted to model or near-model organisms.

Methodology/Principal Findings

This study explores the utility of Diversity Arrays Technology (DArT) in evolutionary studies of non-model organisms. DArT is a hybridization-based genotyping method that uses microarray technology to identify and type DNA polymorphism. Theoretically applicable to any organism (even one for which no prior genetic data are available), DArT has not yet been explored in exclusively wild sample sets, nor extensively examined in a phylogenetic framework. DArT recovered 1349 markers of largely low copy-number loci in two lineages of seed-free land plants: the diploid fern Asplenium viride and the haploid moss Garovaglia elegans. Direct sequencing of 148 of these DArT markers identified 30 putative loci including four routinely sequenced for evolutionary studies in plants. Phylogenetic analyses of DArT genotypes reveal phylogeographic and substrate specificity patterns in A. viride, a lack of phylogeographic pattern in Australian G. elegans, and additive variation in hybrid or mixed samples.

Conclusions/Significance

These results enable methodological recommendations including procedures for detecting and analysing DArT markers tailored specifically to evolutionary investigations and practical factors informing the decision to use DArT, and raise evolutionary hypotheses concerning substrate specificity and biogeographic patterns. Thus DArT is a demonstrably valuable addition to the set of existing molecular approaches used to infer biological phenomena such as adaptive radiations, population dynamics, hybridization, introgression, ecological differentiation and phylogeography.     

Saccharomyces cerevisiae: Population Divergence and Resistance to Oxidative Stress in Clinical,Domesticated and Wild Isolates

Background

Saccharomyces cerevisiae has been associated with human life for millennia in the brewery and bakery. Recently it has been recognized as an emerging opportunistic pathogen. To study the evolutionary history of S. cerevisiae, the origin of clinical isolates and the importance of a virulence-associated trait, population genetics and phenotypic assays have been applied to an ecologically diverse set of 103 strains isolated from clinics, breweries, vineyards, fruits, soil, commercial supplements and insect guts.

Methodology/Principal Findings

DNA sequence data from five nuclear DNA loci were analyzed for population structure and haplotype distribution. Additionally, all strains were tested for survival of oxidative stress, a trait associated with microbial pathogenicity. DNA sequence analyses identified three genetic subgroups within the recombining S. cerevisiae strains that are associated with ecology, geography and virulence. Shared alleles suggest that the clinical isolates contain genetic contribution from the fruit isolates. Clinical and fruit isolates exhibit high levels of recombination, unlike the genetically homogenous soil isolates in which no recombination was detected. However, clinical and soil isolates were more resistant to oxidative stress than any other population, suggesting a correlation between survival in oxidative stress and yeast pathogenicity.

Conclusions/Significance

Population genetic analyses of S. cerevisiae delineated three distinct groups, comprising primarily the (i) human-associated brewery and vineyard strains, (ii) clinical and fruit isolates (iii) and wild soil isolates from eastern U.S. The interactions between S. cerevisiae and humans potentiate yeast evolution and the development of genetically, ecologically and geographically divergent groups.     

Effect of a Standardised Dietary Restriction Protocol on Multiple Laboratory Strains of Drosophila melanogaster

Background

Outcomes of lifespan studies in model organisms are particularly susceptible to variations in technical procedures. This is especially true of dietary restriction, which is implemented in many different ways among laboratories.

Principal Findings

In this study, we have examined the effect of laboratory stock maintenance, genotype differences and microbial infection on the ability of dietary restriction (DR) to extend life in the fruit fly Drosophila melanogaster. None of these factors block the DR effect.

Conclusions

These data lend support to the idea that nutrient restriction genuinely extends lifespan in flies, and that any mechanistic discoveries made with this model are of potential relevance to the determinants of lifespan in other organisms.     

Screening the Budding Yeast Genome Reveals Unique Factors Affecting K2 Toxin Susceptibility

Background

Understanding how biotoxins kill cells is of prime importance in biomedicine and the food industry. The budding yeast (S. cerevisiae) killers serve as a convenient model to study the activity of biotoxins consistently supplying with significant insights into the basic mechanisms of virus-host cell interactions and toxin entry into eukaryotic target cells. K1 and K2 toxins are active at the cell wall, leading to the disruption of the plasma membrane and subsequent cell death by ion leakage. K28 toxin is active in the cell nucleus, blocking DNA synthesis and cell cycle progression, thereby triggering apoptosis. Genome-wide screens in the budding yeast S. cerevisiae identified several hundred effectors of K1 and K28 toxins. Surprisingly, no such screen had been performed for K2 toxin, the most frequent killer toxin among industrial budding yeasts.

Principal Findings

We conducted several concurrent genome-wide screens in S. cerevisiae and identified 332 novel K2 toxin effectors. The effectors involved in K2 resistance and hypersensitivity largely map in distinct cellular pathways, including cell wall and plasma membrane structure/biogenesis and mitochondrial function for K2 resistance, and cell wall stress signaling and ion/pH homeostasis for K2 hypersensitivity. 70% of K2 effectors are different from those involved in K1 or K28 susceptibility.

Significance

Our work demonstrates that despite the fact that K1 and K2 toxins share some aspects of their killing strategies, they largely rely on different sets of effectors. Since the vast majority of the host factors identified here is exclusively active towards K2, we conclude that cells have acquired a specific K2 toxin effectors set. Our work thus indicates that K1 and K2 have elaborated different biological pathways and provides a first step towards the detailed characterization of K2 mode of action.     

A Septin from the Filamentous Fungus A. nidulans Induces Atypical Pseudohyphae in the Budding Yeast S. cerevisiae

Background

Septins, novel cytoskeletal proteins, form rings at the bases of emerging round buds in yeasts and at the bases of emerging elongated hyphal initials in filamentous fungi.

Methodology/Principal Findings

When introduced into the yeast Saccharomyces cerevisiae, the septin AspC from the filamentous fungus Aspergillus nidulans induced highly elongated atypical pseudohyphae and spore-producing structures similar to those of hyphal fungi. AspC induced atypical pseudohyphae when S. cerevisiae pseudohyphal or haploid invasive genes were deleted, but not when the CDC10 septin gene was deleted. AspC also induced atypical pseudohyphae when S. cerevisiae genes encoding Cdc12-interacting proteins Bem4, Cla4, Gic1 and Gic2 were deleted, but not when BNI1, a Cdc12-interacting formin gene, was deleted. AspC localized to bud and pseudohypha necks, while its S. cerevisiae ortholog, Cdc12, localized only to bud necks.

Conclusions/Significance

Our results suggest that AspC competes with Cdc12 for incorporation into the yeast septin scaffold and once there alters cell shape by altering interactions with the formin Bni1. That introduction of the A. nidulans septin AspC into S. cerevisiae induces a shift from formation of buds to formation of atypical pseudohyphae suggests that septins play an important role in the morphological plasticity of fungi.     

Flavin Adenine Dinucleotide Rescues the Phenotype of Frataxin Deficiency

Background

Friedreich ataxia is a neurodegenerative disease caused by the lack of frataxin, a mitochondrial protein. We previously demonstrated that frataxin interacts with complex II subunits of the electronic transport chain (ETC) and putative electronic transfer flavoproteins, suggesting that frataxin could participate in the oxidative phosphorylation.

Methods and Findings

Here we have investigated the effect of riboflavin and its cofactors flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN) in Saccharomyces cerevisiae and Caenorhabditis elegans models of frataxin deficiency. We used a S. cerevisiae strain deleted for the yfh1 gene obtained by homologous recombination and we assessed growth in fermentable and non-fermentable cultures supplemented with either riboflavin or its derivates. Experiments with C. elegans were performed in transient knock-down worms (frh-1[RNAi]) generated by microinjection of dsRNA frh-1 into the gonads of young worms. We observed that FAD rescues the phenotype of both defective organisms. We show that cell growth and enzymatic activities of the ETC complexes and ATP production of yfh1Δ cells were improved by FAD supplementation. Moreover, FAD also improved lifespan and other physiological parameters in the C. elegans knock-down model for frataxin.

Conclusions/Significance

We propose that rescue of frataxin deficiency by FAD supplementation could be explained by an improvement in mitochondrial respiration. We suggest that riboflavin may be useful in the treatment of Friedreich ataxia.     

Random Phenotypic Variation of Yeast (Saccharomyces cerevisiae) Single-Gene Knockouts Fits a Double Pareto-Lognormal Distribution

Background

Distributed robustness is thought to influence the buffering of random phenotypic variation through the scale-free topology of gene regulatory, metabolic, and protein-protein interaction networks. If this hypothesis is true, then the phenotypic response to the perturbation of particular nodes in such a network should be proportional to the number of links those nodes make with neighboring nodes. This suggests a probability distribution approximating an inverse power-law of random phenotypic variation. Zero phenotypic variation, however, is impossible, because random molecular and cellular processes are essential to normal development. Consequently, a more realistic distribution should have a y-intercept close to zero in the lower tail, a mode greater than zero, and a long (fat) upper tail. The double Pareto-lognormal (DPLN) distribution is an ideal candidate distribution. It consists of a mixture of a lognormal body and upper and lower power-law tails.

Objective and Methods

If our assumptions are true, the DPLN distribution should provide a better fit to random phenotypic variation in a large series of single-gene knockout lines than other skewed or symmetrical distributions. We fit a large published data set of single-gene knockout lines in Saccharomyces cerevisiae to seven different probability distributions: DPLN, right Pareto-lognormal (RPLN), left Pareto-lognormal (LPLN), normal, lognormal, exponential, and Pareto. The best model was judged by the Akaike Information Criterion (AIC).

Results

Phenotypic variation among gene knockouts in S. cerevisiae fits a double Pareto-lognormal (DPLN) distribution better than any of the alternative distributions, including the right Pareto-lognormal and lognormal distributions.

Conclusions and Significance

A DPLN distribution is consistent with the hypothesis that developmental stability is mediated, in part, by distributed robustness, the resilience of gene regulatory, metabolic, and protein-protein interaction networks. Alternatively, multiplicative cell growth, and the mixing of lognormal distributions having different variances, may generate a DPLN distribution.     

Rapid Annotation of Anonymous Sequences from Genome Projects Using Semantic Similarities and a Weighting Scheme in Gene Ontology

Background

Large-scale sequencing projects have now become routine lab practice and this has led to the development of a new generation of tools involving function prediction methods, bringing the latter back to the fore. The advent of Gene Ontology, with its structured vocabulary and paradigm, has provided computational biologists with an appropriate means for this task.

Methodology

We present here a novel method called ARGOT (Annotation Retrieval of Gene Ontology Terms) that is able to process quickly thousands of sequences for functional inference. The tool exploits for the first time an integrated approach which combines clustering of GO terms, based on their semantic similarities, with a weighting scheme which assesses retrieved hits sharing a certain number of biological features with the sequence to be annotated. These hits may be obtained by different methods and in this work we have based ARGOT processing on BLAST results.

Conclusions

The extensive benchmark involved 10,000 protein sequences, the complete S. cerevisiae genome and a small subset of proteins for purposes of comparison with other available tools. The algorithm was proven to outperform existing methods and to be suitable for function prediction of single proteins due to its high degree of sensitivity, specificity and coverage.     

Waddlia chondrophila Infects and Multiplies in Ovine Trophoblast Cells Stimulating an Inflammatory Immune Response

Background

Waddlia chondrophila (W. chondrophila) is an emerging abortifacient organism which has been identified in the placentae of humans and cattle. The organism is a member of the order Chlamydiales, and shares many similarities at the genome level and in growth studies with other well-characterised zoonotic chlamydial abortifacients, such as Chlamydia abortus (C. abortus). This study investigates the growth of the organism and its effects upon pro-inflammatory cytokine expression in a ruminant placental cell line which we have previously utilised in a model of C. abortus pathogenicity.

Methodology/Principal Findings

Using qPCR, fluorescent immunocytochemistry and electron microscopy, we characterised the infection and growth of W. chondrophila within the ovine trophoblast AH-1 cell line. Inclusions were visible from 6 h post-infection (p.i.) and exponential growth of the organism could be observed over a 60 h time-course, with significant levels of host cell lysis being observed only after 36 h p.i. Expression of CXCL8, TNF-α, IL-1α and IL-1β were determined 24 h p.i. A statistically significant response in the expression of CXCL8, TNF-α and IL-1β could be observed following active infection with W. chondrophila. However a significant increase in IL-1β expression was also observed following the exposure of cells to UV-killed organisms, indicating the stimulation of multiple innate recognition pathways.

Conclusions/Significance

W. chondrophila infects and grows in the ruminant trophoblast AH-1 cell line exhibiting a complete chlamydial replicative cycle. Infection of the trophoblasts resulted in the expression of pro-inflammatory cytokines in a dose-dependent manner similar to that observed with C. abortus in previous studies, suggesting similarities in the pathogenesis of infection between the two organisms.