Viral Phylodynamics

Viral phylodynamics is defined as the study of how epidemiological, immunological, and evolutionary processes act and potentially interact to shape viral phylogenies. Since the coining of the term in 2004, research on viral phylodynamics has focused on transmission dynamics in an effort to shed light on how these dynamics impact viral genetic variation. Transmission dynamics can be considered at the level of cells within an infected host, individual hosts within a population, or entire populations of hosts. Many viruses, especially RNA viruses, rapidly accumulate genetic variation because of short generation times and high mutation rates. Patterns of viral genetic variation are therefore heavily influenced by how quickly transmission occurs and by which entities transmit to one another. Patterns of viral genetic variation will also be affected by selection acting on viral phenotypes. Although viruses can differ with respect to many phenotypes, phylodynamic studies have to date tended to focus on a limited number of viral phenotypes. These include virulence phenotypes, phenotypes associated with viral transmissibility, cell or tissue tropism phenotypes, and antigenic phenotypes that can facilitate escape from host immunity. Due to the impact that transmission dynamics and selection can have on viral genetic variation, viral phylogenies can therefore be used to investigate important epidemiological, immunological, and evolutionary processes, such as epidemic spread [2], spatio-temporal dynamics including metapopulation dynamics [3], zoonotic transmission, tissue tropism [4], and antigenic drift [5]. The quantitative investigation of these processes through the consideration of viral phylogenies is the central aim of viral phylodynamics.

This is a “Topic Page” article for PLOS Computational Biology.

     

Approximate Bayesian Computation

Approximate Bayesian computation (ABC) constitutes a class of computational methods rooted in Bayesian statistics. In all model-based statistical inference, the likelihood function is of central importance, since it expresses the probability of the observed data under a particular statistical model, and thus quantifies the support data lend to particular values of parameters and to choices among different models. For simple models, an analytical formula for the likelihood function can typically be derived. However, for more complex models, an analytical formula might be elusive or the likelihood function might be computationally very costly to evaluate. ABC methods bypass the evaluation of the likelihood function. In this way, ABC methods widen the realm of models for which statistical inference can be considered. ABC methods are mathematically well-founded, but they inevitably make assumptions and approximations whose impact needs to be carefully assessed. Furthermore, the wider application domain of ABC exacerbates the challenges of parameter estimation and model selection. ABC has rapidly gained popularity over the last years and in particular for the analysis of complex problems arising in biological sciences (e.g., in population genetics, ecology, epidemiology, and systems biology).

This is a “Topic Page” article for PLOS Computational Biology.

     

Flow Cytometry Bioinformatics

Flow cytometry bioinformatics is the application of bioinformatics to flow cytometry data, which involves storing, retrieving, organizing, and analyzing flow cytometry data using extensive computational resources and tools. Flow cytometry bioinformatics requires extensive use of and contributes to the development of techniques from computational statistics and machine learning. Flow cytometry and related methods allow the quantification of multiple independent biomarkers on large numbers of single cells. The rapid growth in the multidimensionality and throughput of flow cytometry data, particularly in the 2000s, has led to the creation of a variety of computational analysis methods, data standards, and public databases for the sharing of results. Computational methods exist to assist in the preprocessing of flow cytometry data, identifying cell populations within it, matching those cell populations across samples, and performing diagnosis and discovery using the results of previous steps. For preprocessing, this includes compensating for spectral overlap, transforming data onto scales conducive to visualization and analysis, assessing data for quality, and normalizing data across samples and experiments. For population identification, tools are available to aid traditional manual identification of populations in two-dimensional scatter plots (gating), to use dimensionality reduction to aid gating, and to find populations automatically in higher dimensional space in a variety of ways. It is also possible to characterize data in more comprehensive ways, such as the density-guided binary space partitioning technique known as probability binning, or by combinatorial gating. Finally, diagnosis using flow cytometry data can be aided by supervised learning techniques, and discovery of new cell types of biological importance by high-throughput statistical methods, as part of pipelines incorporating all of the aforementioned methods. Open standards, data, and software are also key parts of flow cytometry bioinformatics. Data standards include the widely adopted Flow Cytometry Standard (FCS) defining how data from cytometers should be stored, but also several new standards under development by the International Society for Advancement of Cytometry (ISAC) to aid in storing more detailed information about experimental design and analytical steps. Open data is slowly growing with the opening of the CytoBank database in 2010 and FlowRepository in 2012, both of which allow users to freely distribute their data, and the latter of which has been recommended as the preferred repository for MIFlowCyt-compliant data by ISAC. Open software is most widely available in the form of a suite of Bioconductor packages, but is also available for web execution on the GenePattern platform.

This is a “Topic Page” article for PLOS Computational Biology.

     

Remodeling of the Rad51 DNA Strand-Exchange Protein by the Srs2 Helicase

Homologous recombination is associated with the dynamic assembly and disassembly of DNA–protein complexes. Assembly of a nucleoprotein filament comprising ssDNA and the RecA homolog, Rad51, is a key step required for homology search during recombination. The budding yeast Srs2 DNA translocase is known to dismantle Rad51 filament in vitro. However, there is limited evidence to support the dismantling activity of Srs2 in vivo. Here, we show that Srs2 indeed disrupts Rad51-containing complexes from chromosomes during meiosis. Overexpression of Srs2 during the meiotic prophase impairs meiotic recombination and removes Rad51 from meiotic chromosomes. This dismantling activity is specific for Rad51, as Srs2 Overexpression does not remove Dmc1 (a meiosis-specific Rad51 homolog), Rad52 (a Rad51 mediator), or replication protein A (RPA; a single-stranded DNA-binding protein). Rather, RPA replaces Rad51 under these conditions. A mutant Srs2 lacking helicase activity cannot remove Rad51 from meiotic chromosomes. Interestingly, the Rad51-binding domain of Srs2, which is critical for Rad51-dismantling activity in vitro, is not essential for this activity in vivo. Our results suggest that a precise level of Srs2, in the form of the Srs2 translocase, is required to appropriately regulate the Rad51 nucleoprotein filament dynamics during meiosis.     

High quality draft genome sequence of the slightly halophilic bacterium Halomonas zhanjiangensis type strain JSM 078169T (DSM 21076T) from a sea urchin in southern China

Halomonas zhanjiangensis Chen et al. 2009 is a member of the genus Halomonas, family Halomonadaceae, class Gammaproteobacteria. Representatives of the genus Halomonas are a group of halophilic bacteria often isolated from salty environments. The type strain H. zhanjiangensis JSM 078169^T was isolated from a sea urchin (Hemicentrotus pulcherrimus) collected from the South China Sea. The genome of strain JSM 078169^T is the fourteenth sequenced genome in the genus Halomonas and the fifteenth in the family Halomonadaceae. The other thirteen genomes from the genus Halomonas are H. halocynthiae, H. venusta, H. alkaliphila, H. lutea, H. anticariensis, H. jeotgali, H. titanicae, H. desiderata, H. smyrnensis, H. salifodinae, H. boliviensis, H. elongata and H stevensii. Here, we describe the features of strain JSM 078169^T, together with the complete genome sequence and annotation from a culture of DSM 21076^T. The 4,060,520 bp long draft genome consists of 17 scaffolds with the 3,659 protein-coding and 80 RNA genes and is a part of Genomic Encyclopedia of Type Strains, Phase I: the one thousand microbial genomes (KMG) project.     

Integrity and Function of the Saccharomyces cerevisiae Spindle Pole Body Depends on Connections Between the Membrane Proteins Ndc1, Rtn1, and Yop1

The nuclear envelope in Saccharomyces cerevisiae harbors two essential macromolecular protein assemblies: the nuclear pore complexes (NPCs) that enable nucleocytoplasmic transport, and the spindle pole bodies (SPBs) that mediate chromosome segregation. Previously, based on metazoan and budding yeast studies, we reported that reticulons and Yop1/DP1 play a role in the early steps of de novo NPC assembly. Here, we examined if Rtn1 and Yop1 are required for SPB function in S. cerevisiae. Electron microscopy of rtn1Δ yop1Δ cells revealed lobular abnormalities in SPB structure. Using an assay that monitors lateral expansion of the SPB central layer, we found that rtn1Δ yop1Δ SPBs had decreased connections to the NE compared to wild type, suggesting that SPBs are less stable in the NE. Furthermore, large budded rtn1Δ yop1Δ cells exhibited a high incidence of short mitotic spindles, which were frequently misoriented with respect to the mother–daughter axis. This correlated with cytoplasmic microtubule defects. We found that overexpression of the SPB insertion factors NDC1, MPS2, or BBP1 rescued the SPB defects observed in rtn1Δ yop1Δ cells. However, only overexpression of NDC1, which is also required for NPC biogenesis, rescued both the SPB and NPC associated defects. Rtn1 and Yop1 also physically interacted with Ndc1 and other NPC membrane proteins. We propose that NPC and SPB biogenesis are altered in cells lacking Rtn1 and Yop1 due to competition between these complexes for Ndc1, an essential common component of both NPCs and SPBs.     

Translational Control of the Oogenic Program by Components of OMA Ribonucleoprotein Particles in Caenorhabditis elegans

The oocytes of most sexually reproducing animals arrest in meiotic prophase I. Oocyte growth, which occurs during this period of arrest, enables oocytes to acquire the cytoplasmic components needed to produce healthy progeny and to gain competence to complete meiosis. In the nematode Caenorhabditis elegans, the major sperm protein hormone promotes meiotic resumption (also called meiotic maturation) and the cytoplasmic flows that drive oocyte growth. Prior work established that two related TIS11 zinc-finger RNA-binding proteins, OMA-1 and OMA-2, are redundantly required for normal oocyte growth and meiotic maturation. We affinity purified OMA-1 and identified associated mRNAs and proteins using genome-wide expression data and mass spectrometry, respectively. As a class, mRNAs enriched in OMA-1 ribonucleoprotein particles (OMA RNPs) have reproductive functions. Several of these mRNAs were tested and found to be targets of OMA-1/2-mediated translational repression, dependent on sequences in their 3′-untranslated regions (3′-UTRs). Consistent with a major role for OMA-1 and OMA-2 in regulating translation, OMA-1-associated proteins include translational repressors and activators, and some of these proteins bind directly to OMA-1 in yeast two-hybrid assays, including OMA-2. We show that the highly conserved TRIM-NHL protein LIN-41 is an OMA-1-associated protein, which also represses the translation of several OMA-1/2 target mRNAs. In the accompanying article in this issue, we show that LIN-41 prevents meiotic maturation and promotes oocyte growth in opposition to OMA-1/2. Taken together, these data support a model in which the conserved regulators of mRNA translation LIN-41 and OMA-1/2 coordinately control oocyte growth and the proper spatial and temporal execution of the meiotic maturation decision.     

Chemical graph generators

Chemical graph generators are software packages to generate computer representations of chemical structures adhering to certain boundary conditions. Their development is a research topic of cheminformatics. Chemical graph generators are used in areas such as virtual library generation in drug design, in molecular design with specified properties, called inverse QSAR/QSPR, as well as in organic synthesis design, retrosynthesis or in systems for computer-assisted structure elucidation (CASE). CASE systems again have regained interest for the structure elucidation of unknowns in computational metabolomics, a current area of computational biology.     

High quality draft genome sequence of Staphylococcus cohnii subsp. cohnii strain hu-01

Staphylococcus cohnii subsp. cohnii belongs to the family Staphylococcaceae in the order Bacillales, class Bacilli and phylum Firmicutes. The increasing relevance of S. cohnii to human health prompted us to determine the genomic sequence of Staphylococcus cohnii subsp. cohnii strain hu-01, a multidrug-resistant isolate from a hospital in China. Here we describe the features of S. cohnii subsp. cohnii strain hu-01, together with the genome sequence and its annotation. This is the first genome sequence of the species Staphylococcus cohnii.     

The Mechanism of Nucleotide Excision Repair-Mediated UV-Induced Mutagenesis in Nonproliferating Cells

Following the irradiation of nondividing yeast cells with ultraviolet (UV) light, most induced mutations are inherited by both daughter cells, indicating that complementary changes are introduced into both strands of duplex DNA prior to replication. Early analyses demonstrated that such two-strand mutations depend on functional nucleotide excision repair (NER), but the molecular mechanism of this unique type of mutagenesis has not been further explored. In the experiments reported here, an ade2 adeX colony-color system was used to examine the genetic control of UV-induced mutagenesis in nondividing cultures of Saccharomyces cerevisiae. We confirmed a strong suppression of two-strand mutagenesis in NER-deficient backgrounds and demonstrated that neither mismatch repair nor interstrand crosslink repair affects the production of these mutations. By contrast, proteins involved in the error-prone bypass of DNA damage (Rev3, Rev1, PCNA, Rad18, Pol32, and Rad5) and in the early steps of the DNA-damage checkpoint response (Rad17, Mec3, Ddc1, Mec1, and Rad9) were required for the production of two-strand mutations. There was no involvement, however, for the Pol η translesion synthesis DNA polymerase, the Mms2-Ubc13 postreplication repair complex, downstream DNA-damage checkpoint factors (Rad53, Chk1, and Dun1), or the Exo1 exonuclease. Our data support models in which UV-induced mutagenesis in nondividing cells occurs during the Pol ζ-dependent filling of lesion-containing, NER-generated gaps. The requirement for specific DNA-damage checkpoint proteins suggests roles in recruiting and/or activating factors required to fill such gaps.     

DNA Replication Error-Induced Extinction of Diploid Yeast

Genetic defects in DNA polymerase accuracy, proofreading, or mismatch repair (MMR) induce mutator phenotypes that accelerate adaptation of microbes and tumor cells. Certain combinations of mutator alleles synergistically increase mutation rates to levels that drive extinction of haploid cells. The maximum tolerated mutation rate of diploid cells is unknown. Here, we define the threshold for replication error-induced extinction (EEX) of diploid Saccharomyces cerevisiae. Double-mutant pol3 alleles that carry mutations for defective DNA polymerase-δ proofreading (pol3-01) and accuracy (pol3-L612M or pol3-L612G) induce strong mutator phenotypes in heterozygous diploids (POL3/pol3-01,L612M or POL3/pol3-01,L612G). Both pol3-01,L612M and pol3-01,L612G alleles are lethal in the homozygous state; cells with pol3-01,L612M divide up to 10 times before arresting at random stages in the cell cycle. Antimutator eex mutations in the pol3 alleles suppress this lethality (pol3-01,L612M,eex or pol3-01,L612G,eex). MMR defects synergize with pol3-01,L612M,eex and pol3-01,L612G,eex alleles, increasing mutation rates and impairing growth. Conversely, inactivation of the Dun1 S-phase checkpoint kinase suppresses strong pol3-01,L612M,eex and pol3-01,L612G,eex mutator phenotypes as well as the lethal pol3-01,L612M phenotype. Our results reveal that the lethal error threshold in diploids is 10 times higher than in haploids and likely determined by homozygous inactivation of essential genes. Pronounced loss of fitness occurs at mutation rates well below the lethal threshold, suggesting that mutator-driven cancers may be susceptible to drugs that exacerbate replication errors.     

Asymmetric Neuroblast Divisions Producing Apoptotic Cells Require the Cytohesin GRP-1 in Caenorhabditis elegans

Cytohesins are Arf guanine nucleotide exchange factors (GEFs) that regulate membrane trafficking and actin cytoskeletal dynamics. We report here that GRP-1, the sole Caenorhabditis elegans cytohesin, controls the asymmetric divisions of certain neuroblasts that divide to produce a larger neuronal precursor or neuron and a smaller cell fated to die. In the Q neuroblast lineage, loss of GRP-1 led to the production of daughter cells that are more similar in size and to the transformation of the normally apoptotic daughter into its sister, resulting in the production of extra neurons. Genetic interactions suggest that GRP-1 functions with the previously described Arf GAP CNT-2 and two other Arf GEFs, EFA-6 and BRIS-1, to regulate the activity of Arf GTPases. In agreement with this model, we show that GRP-1’s GEF activity, mediated by its SEC7 domain, is necessary for the posterior Q cell (Q.p) neuroblast division and that both GRP-1 and CNT-2 function in the Q.posterior Q daughter cell (Q.p) to promote its asymmetry. Although functional GFP-tagged GRP-1 proteins localized to the nucleus, the extra cell defects were rescued by targeting the Arf GEF activity of GRP-1 to the plasma membrane, suggesting that GRP-1 acts at the plasma membrane. The detection of endogenous GRP-1 protein at cytokinesis remnants, or midbodies, is consistent with GRP-1 functioning at the plasma membrane and perhaps at the cytokinetic furrow to promote the asymmetry of the divisions that require its function.     

Replisome Function During Replicative Stress Is Modulated by Histone H3 Lysine 56 Acetylation Through Ctf4

Histone H3 lysine 56 acetylation in Saccharomyces cerevisiae is required for the maintenance of genome stability under normal conditions and upon DNA replication stress. Here we show that in the absence of H3 lysine 56 acetylation replisome components become deleterious when replication forks collapse at natural replication block sites. This lethality is not a direct consequence of chromatin assembly defects during replication fork progression. Rather, our genetic analyses suggest that in the presence of replicative stress H3 lysine 56 acetylation uncouples the Cdc45–Mcm2-7–GINS DNA helicase complex and DNA polymerases through the replisome component Ctf4. In addition, we discovered that the N-terminal domain of Ctf4, necessary for the interaction of Ctf4 with Mms22, an adaptor protein of the Rtt101-Mms1 E3 ubiquitin ligase, is required for the function of the H3 lysine 56 acetylation pathway, suggesting that replicative stress promotes the interaction between Ctf4 and Mms22. Taken together, our results indicate that Ctf4 is an essential member of the H3 lysine 56 acetylation pathway and provide novel mechanistic insights into understanding the role of H3 lysine 56 acetylation in maintaining genome stability upon replication stress.     

Complete genome sequence of Coriobacterium glomerans type strain (PW2T) from the midgut of Pyrrhocoris apterus L. (red soldier bug)

Coriobacterium glomerans Haas and König 1988, is the only species of the genus Coriobacterium, family Coriobacteriaceae, order Coriobacteriales, phylum Actinobacteria. The bacterium thrives as an endosymbiont of pyrrhocorid bugs, i.e. the red fire bug Pyrrhocoris apterus L. The rationale for sequencing the genome of strain PW2^T is its endosymbiotic life style which is rare among members of Actinobacteria. Here we describe the features of this symbiont, together with the complete genome sequence and its annotation. This is the first complete genome sequence of a member of the genus Coriobacterium and the sixth member of the order Coriobacteriales for which complete genome sequences are now available. The 2,115,681 bp long single replicon genome with its 1,804 protein-coding and 54 RNA genes is part of the Genomic Encyclopedia of Bacteria and Archaea project.     

Genome sequence of the phylogenetically isolated spirochete Leptonema illini type strain (3055T)

Leptonema illini Hovind-Hougen 1979 is the type species of the genus Leptonema, family Leptospiraceae, phylum Spirochaetes. Organisms of this family have a Gram-negative-like cell envelope consisting of a cytoplasmic membrane and an outer membrane. The peptidoglycan layer is associated with the cytoplasmic rather than the outer membrane. The two flagella of members of Leptospiraceae extend from the cytoplasmic membrane at the ends of the bacteria into the periplasmic space and are necessary for their motility. Here we describe the features of the L. illini type strain, together with the complete genome sequence, and annotation. This is the first genome sequence (finished at the level of Improved High Quality Draft) to be reported from of a member of the genus Leptonema and a representative of the third genus of the family Leptospiraceae for which complete or draft genome sequences are now available. The three scaffolds of the 4,522,760 bp draft genome sequence reported here, and its 4,230 protein-coding and 47 RNA genes are part of the Genomic Encyclopedia of Bacteria and Archaea project.     

Systematic Analyses of rpm-1 Suppressors Reveal Roles for ESS-2 in mRNA Splicing in Caenorhabditis elegans

The PHR (Pam/Highwire/RPM-1) family of ubiquitin E3 ligases plays conserved roles in axon patterning and synaptic development. Genetic modifier analysis has greatly aided the discovery of the signal transduction cascades regulated by these proteins. In Caenorhabditis elegans, loss of function in rpm-1 causes axon overgrowth and aberrant presynaptic morphology, yet the mutant animals exhibit little behavioral deficits. Strikingly, rpm-1 mutations strongly synergize with loss of function in the presynaptic active zone assembly factors, syd-1 and syd-2, resulting in severe locomotor deficits. Here, we provide ultrastructural evidence that double mutants, between rpm-1 and syd-1 or syd-2, dramatically impair synapse formation. Taking advantage of the synthetic locomotor defects to select for genetic suppressors, previous studies have identified the DLK-1 MAP kinase cascade negatively regulated by RPM-1. We now report a comprehensive analysis of a large number of suppressor mutations of this screen. Our results highlight the functional specificity of the DLK-1 cascade in synaptogenesis. We also identified two previously uncharacterized genes. One encodes a novel protein, SUPR-1, that acts cell autonomously to antagonize RPM-1. The other affects a conserved protein ESS-2, the homolog of human ES2 or DGCR14. Loss of function in ess-2 suppresses rpm-1 only in the presence of a dlk-1 splice acceptor mutation. We show that ESS-2 acts to promote accurate mRNA splicing when the splice site is compromised. The human DGCR14/ES2 resides in a deleted chromosomal region implicated in DiGeorge syndrome, and its mutation has shown high probability as a risk factor for schizophrenia. Our findings provide the first functional evidence that this family of proteins regulate mRNA splicing in a context-specific manner.