Computing folding pathways between RNA secondary structures

Given an RNA sequence and two designated secondary structures A, B, we describe a new algorithm that computes a nearly optimal folding pathway from A to B. The algorithm, RNAtabupath, employs a tabu semi-greedy heuristic, known to be an effective search strategy in combinatorial optimization. Folding pathways, sometimes called routes or trajectories, are computed by RNAtabupath in a fraction of the time required by the barriers program of Vienna RNA Package. We benchmark RNAtabupath with other algorithms to compute low energy folding pathways between experimentally known structures of several conformational switches. The RNApathfinder web server, source code for algorithms to compute and analyze pathways and supplementary data are available at http://bioinformatics.bc.edu/clotelab/RNApathfinder.     

SWALO: scaffolding with assembly likelihood optimization

Scaffolding, i.e. ordering and orienting contigs is an important step in genome assembly. We present a method for scaffolding using second generation sequencing reads based on likelihoods of genome assemblies. A generative model for sequencing is used to obtain maximum likelihood estimates of gaps between contigs and to estimate whether linking contigs into scaffolds would lead to an increase in the likelihood of the assembly. We then link contigs if they can be unambiguously joined or if the corresponding increase in likelihood is substantially greater than that of other possible joins of those contigs. The method is implemented in a tool called Swalo with approximations to make it efficient and applicable to large datasets. Analysis on real and simulated datasets reveals that it consistently makes more or similar number of correct joins as other scaffolders while linking very few contigs incorrectly, thus outperforming other scaffolders and demonstrating that substantial improvement in genome assembly may be achieved through the use of statistical models. Swalo is freely available for download at https://atifrahman.github.io/SWALO/.     

Accurate,Fully-Automated NMR Spectral Profiling for Metabolomics

Many diseases cause significant changes to the concentrations of small molecules (a.k.a. metabolites) that appear in a person’s biofluids, which means such diseases can often be readily detected from a person’s “metabolic profile"—i.e., the list of concentrations of those metabolites. This information can be extracted from a biofluids Nuclear Magnetic Resonance (NMR) spectrum. However, due to its complexity, NMR spectral profiling has remained manual, resulting in slow, expensive and error-prone procedures that have hindered clinical and industrial adoption of metabolomics via NMR. This paper presents a system, BAYESIL, which can quickly, accurately, and autonomously produce a person’s metabolic profile. Given a 1D ¹H NMR spectrum of a complex biofluid (specifically serum or cerebrospinal fluid), BAYESIL can automatically determine the metabolic profile. This requires first performing several spectral processing steps, then matching the resulting spectrum against a reference compound library, which contains the “signatures” of each relevant metabolite. BAYESIL views spectral matching as an inference problem within a probabilistic graphical model that rapidly approximates the most probable metabolic profile. Our extensive studies on a diverse set of complex mixtures including real biological samples (serum and CSF), defined mixtures and realistic computer generated spectra; involving > 50 compounds, show that BAYESIL can autonomously find the concentration of NMR-detectable metabolites accurately (~ 90% correct identification and ~ 10% quantification error), in less than 5 minutes on a single CPU. These results demonstrate that BAYESIL is the first fully-automatic publicly-accessible system that provides quantitative NMR spectral profiling effectively—with an accuracy on these biofluids that meets or exceeds the performance of trained experts. We anticipate this tool will usher in high-throughput metabolomics and enable a wealth of new applications of NMR in clinical settings. BAYESIL is accessible at http://www.bayesil.ca.     

A Real-Time All-Atom Structural Search Engine for Proteins

Protein designers use a wide variety of software tools for de novo design, yet their repertoire still lacks a fast and interactive all-atom search engine. To solve this, we have built the Suns program: a real-time, atomic search engine integrated into the PyMOL molecular visualization system. Users build atomic-level structural search queries within PyMOL and receive a stream of search results aligned to their query within a few seconds. This instant feedback cycle enables a new “designability”-inspired approach to protein design where the designer searches for and interactively incorporates native-like fragments from proven protein structures. We demonstrate the use of Suns to interactively build protein motifs, tertiary interactions, and to identify scaffolds compatible with hot-spot residues. The official web site and installer are located at http://www.degradolab.org/suns/ and the source code is hosted at https://github.com/godotgildor/Suns (PyMOL plugin, BSD license), https://github.com/Gabriel439/suns-cmd (command line client, BSD license), and https://github.com/Gabriel439/suns-search (search engine server, GPLv2 license).

This is a PLOS Computational Biology Software Article

     

Caenorhabditis elegans Fibroblast Growth Factor Receptor Signaling Can Occur Independently of the Multi-Substrate Adaptor FRS2

The components of receptor tyrosine kinase signaling complexes help to define the specificity of the effects of their activation. The Caenorhabditis elegans fibroblast growth factor receptor (FGFR), EGL-15, regulates a number of processes, including sex myoblast (SM) migration guidance and fluid homeostasis, both of which require a Grb2/Sos/Ras cassette of signaling components. Here we show that SEM-5/Grb2 can bind directly to EGL-15 to mediate SM chemoattraction. A yeast two-hybrid screen identified SEM-5 as able to interact with the carboxy-terminal domain (CTD) of EGL-15, a domain that is specifically required for SM chemoattraction. This interaction requires the SEM-5 SH2-binding motifs present in the CTD (Y¹⁰⁰⁹ and Y¹⁰⁸⁷), and these sites are required for the CTD role of EGL-15 in SM chemoattraction. SEM-5, but not the SEM-5 binding sites located in the CTD, is required for the fluid homeostasis function of EGL-15, indicating that SEM-5 can link to EGL-15 through an alternative mechanism. The multi-substrate adaptor protein FRS2 serves to link vertebrate FGFRs to Grb2. In C. elegans, an FRS2-like gene, rog-1, functions upstream of a Ras/MAPK pathway for oocyte maturation but is not required for EGL-15 function. Thus, unlike the vertebrate FGFRs, which require the multi-substrate adaptor FRS2 to recruit Grb2, EGL-15 can recruit SEM-5/Grb2 directly.FIBROBLAST growth factors (FGFs) play important roles in many developmental and physiological processes, including cell migration, angiogenesis, proliferation, differentiation, and survival (Ornitz and Itoh 2001; Polanska et al. 2009). Mammals have a battery of both FGF ligands and high-affinity receptors to carry out this diverse set of important functions. These ligands and their receptors are generated from a set of 18 genes encoding FGFs and 4 genes encoding their receptors (Eswarakumar et al. 2005). Upon ligand binding, fibroblast growth factor receptors (FGFRs) dimerize, activating their intrinsic tyrosine kinase activity, which causes both autophosphorylation on intracellular tyrosine residues and phosphorylation of additional substrates (Eswarakumar et al. 2005). These phosphorylation events lead to the assembly of a signaling complex around the activated receptor, ultimately promoting various downstream signaling pathways (Eswarakumar et al. 2005).A large portion of mammalian FGFR signaling is mediated by the multi-substrate adaptor protein FRS2/snt-1 (Kouhara et al. 1997; Hadari et al. 1998, 2001; Lax et al. 2002; Gotoh et al. 2005). FRS2 constitutively associates with the juxtamembrane region of the FGFR via its amino-terminal PTB domain (Xu et al. 1998; Ong et al. 2000). Upon FGFR activation, FRS2 becomes heavily phosphorylated, allowing it to recruit Grb2 and Shp2 via their SH2 domains (Kouhara et al. 1997; Eswarakumar et al. 2005). Since these components cannot associate with the receptor in the absence of FRS2 (Hadari et al. 2001), FRS2 serves as an essential link between the activated receptor and many downstream signal transduction pathways.The understanding of FGF-stimulated signal transduction pathways has been aided by the study of FGF signaling in model organisms. Powerful genetic screens and the reduced complexity of the set of FGFs and their receptors in both Drosophila melanogaster and Caenorhabditis elegans have helped promote an understanding of the conserved aspects of FGF signaling pathways (Huang and Stern 2005; Polanska et al. 2009). In C. elegans, FGF signaling is mediated by two FGF ligands, EGL-17 and LET-756, and a single FGF receptor, EGL-15 (DeVore et al. 1995; Burdine et al. 1997; Roubin et al. 1999). The EGL-15 FGFR is structurally very similar to mammalian FGF receptors, with the highest level of sequence conservation found within the intracellular tyrosine kinase domain and the three extracellular immunoglobulin (IG) domains.Similar to mammalian FGFRs, alternative splicing also generates functionally distinct EGL-15 isoforms. A major structural difference between EGL-15 and other FGFRs lies in an additional domain located between the first IG domain and the acid box of EGL-15. This EGL-15-specific insert is encoded by a pair of mutually exclusive fifth exons, generating two EGL-15 isoforms, 5A and 5B, with different functions (Goodman et al. 2003). Alternative splicing also affects the sequence at the very end of the carboxy-terminal domain (CTD) of EGL-15 (Goodman et al. 2003), giving rise to four distinct C-terminal isoform types, referred to as types I–IV (see Figure 1A).Open in a separate windowFigure 1.—The SEM-5∷EGL-15 interaction is dependent upon the SEM-5 SH2 binding sites in the CTD of EGL-15. (A) Important residues in the EGL-15 CTD isoforms. The C-terminal portion of the kinase domain and the CTD are shown. The gray box is common to all C-terminal isoforms. Sequences of the CTD isoforms can be found in Figure S5. (B) SEM-5 binds directly to EGL-15 Y1009 and Y1087. A full-length SEM-5 prey was mated to eight EGL-15 baits: empty vector control (pBTM116) and derivatives containing variants of either the type I or the type IV EGL-15(Intra). Growth is an indication of an interaction between SEM-5 and the EGL-15 bait. Dilutions of the culture mixture are indicated above. hFGFR1, a bait containing the corresponding portion of the intracellular domain of the human FGFR1. Similar results were obtained using the human SEM-5 ortholog, Grb2.EGL-15, like its mammalian orthologs, is also involved in a large variety of functions, including cell migration guidance affecting the sex myoblast (SM) (Stern and Horvitz 1991; DeVore et al. 1995; Goodman et al. 2003) and CAN cells (Fleming et al. 2005), muscle arm extension (Dixon et al. 2006), a number of processes controlling terminal axon morphology (Bulow et al. 2004), muscle protein degradation (Szewczyk and Jacobson 2003), and fluid homeostasis (Huang and Stern 2004). A conserved FGFR signaling pathway in C. elegans was established by identifying the genes necessary for the role of EGL-15 in fluid homeostasis, and much of this same pathway is utilized in other functions of EGL-15 (DeVore et al. 1995; Borland et al. 2001; Huang and Stern 2004, 2005). The identification of these genes was facilitated by a temperature-sensitive mutation affecting a receptor tyrosine phosphatase, CLR-1 (CLeaR), which functions to negatively regulate EGL-15 signaling (Kokel et al. 1998). Mutations in clr-1 abolish this regulatory constraint on EGL-15, resulting in fluid accumulation in the pseudocoelomic cavity due to hyperactive EGL-15 signaling. The buildup of this clear fluid, resulting in the Clr phenotype, is easily scored and can be used to identify suppressors (soc, suppressor of Clr) that reduce EGL-15 signaling efficiency. These suppressors define an EGL-15 signaling pathway necessary for fluid homeostasis, which involves the activation of the Ras/MAPK cascade via the SEM-5/Grb2 adaptor protein, the let-341/sos-1 Sos-like guanine nucleotide exchange factor, and a PTP-2-SOC-1/Shp2-Gab1 cassette (Borland et al. 2001).Mutations in egl-15 have been identified on the basis of their effects on either fluid homeostasis or the guidance of the migrating SMs (DeVore et al. 1995; Goodman et al. 2003). Most of these alleles fall into an allelic series and affect the general aspects of EGL-15 signaling (Goodman et al. 2003). The strongest of these alleles confers an early larval arrest phenotype, whereas weaker alleles confer either a scrawny body morphology (Scr) or just the ability to suppress the Clr phenotype (Soc). While the weakest of these alleles does not affect egg laying, more highly compromised mutants show an egg-laying defect due to the mispositioning of the SMs. Four egl-15 alleles specifically affect SM migration; when homozygous, these mutations cause dramatic mispositioning of the SMs, but do not cause a Soc phenotype (Goodman et al. 2003). Three of these are nonsense mutations in exon 5A and eliminate the 5A EGL-15 isoform. The phenotype of these mutants highlights the specific requirement of the 5A isoform for SM migration guidance. The fourth mutation in this class, egl-15(n1457), is a nonsense mutation that truncates the carboxy-terminal domain, specifically implicating the CTD in SM migration guidance.Immediately following their birth at the end of the first larval stage (L1), the two bilaterally symmetric SMs undergo anteriorly directed migrations to final positions that flank the precise center of the gonad (Sulston and Horvitz 1977). In the middle of the third larval stage, the SMs divide to generate 16 cells that differentiate into the egg-laying muscles. Multiple mechanisms help guide the migrations of the SMs (Chen and Stern 1998; Branda and Stern 2000), including a chemoattraction mediated by EGL-15 that guides the SMs to their precise final positions (Burdine et al. 1998). The EGL-17 FGF serves as the chemoattractive cue, emanating from central gonadal cells (Branda and Stern 2000). In the absence of this chemoattraction, SMs are posteriorly displaced (Stern and Horvitz 1991). While mispositioned SMs still generate sex muscles, these muscles end up too far posterior to attach properly, causing the animal to be defective in egg laying (Egl).The signal transduction pathway downstream of EGL-15 that mediates SM chemoattraction is not well established. Several lines of evidence implicate SEM-5/Grb2, LET-341/Sos, and LET-60/Ras in SM chemoattraction. The roles of components in the Ras-MAPK cascade in this event are less clear (Sundaram et al. 1996; Chen et al. 1997; Chen and Stern 1998). A crucial gap in our understanding lies in the link between activated EGL-15 and the downstream signaling components. Here we show that SEM-5, the C. elegans GRB2 ortholog, appears to bind directly to SH2 binding sites within the carboxy terminal tail of EGL-15. These interactions are required for SM chemoattraction, but not for the essential function of EGL-15.     

c-Type Cytochrome Assembly in Saccharomyces cerevisiae: A Key Residue for Apocytochrome c1/Lyase Interaction

The electron transport chains in the membranes of bacteria and organelles generate proton-motive force essential for ATP production. The c-type cytochromes, defined by the covalent attachment of heme to a CXXCH motif, are key electron carriers in these energy-transducing membranes. In mitochondria, cytochromes c and c₁ are assembled by the cytochrome c heme lyases (CCHL and CC₁HL) and by Cyc2p, a putative redox protein. A cytochrome c₁ mutant with a CAPCH heme-binding site instead of the wild-type CAACH is strictly dependent upon Cyc2p for assembly. In this context, we found that overexpression of CC₁HL, as well as mutations of the proline in the CAPCH site to H, L, S, or T residues, can bypass the absence of Cyc2p. The P mutation was postulated to shift the CXXCH motif to an oxidized form, which must be reduced in a Cyc2p-dependent reaction before heme ligation. However, measurement of the redox midpoint potential of apocytochrome c₁ indicates that neither the P nor the T residues impact the thermodynamic propensity of the CXXCH motif to occur in a disulfide vs. dithiol form. We show instead that the identity of the second intervening residue in the CXXCH motif is key in determining the CCHL-dependent vs. CC₁HL-dependent assembly of holocytochrome c₁. We also provide evidence that Cyc2p is dedicated to the CCHL pathway and is not required for the CC₁HL-dependent assembly of cytochrome c₁.THE c-type cytochromes, also referred to as cytochrome c, represent a universal class of heme-containing proteins that function as electron carriers in the energy-transducing pathways of bacteria, plastids, and mitochondria (Thöny-Meyer 1997; Nakamoto et al. 2000; Bonnard et al. 2010). Because cytochromes c carry a heme covalently attached to a CXXCH motif, they constitute an attractive object of study to address the question of cofactor protein assembly. The biochemical requirements for cytochrome c assembly were deduced from in vivo and in vitro studies, and the conclusion is that both apocytochromes c and heme are transported independently across at least one biological membrane and maintained as reduced prior to catalysis of the heme attachment reaction (Allen et al. 2003; Hamel et al. 2009; Kranz et al. 2009; Sanders et al. 2010). Bacterial cytochromes c are assembled in the periplasmic space, a compartment where cysteine pairs in proteins form disulfide bonds in reactions catalyzed by dedicated enzymes (Inaba 2009; Kadokura and Beckwith 2010). The current thinking holds that a c-type apocytochrome is a substrate of the disulfide bond-forming pathway, which introduces an intramolecular disulfide between the two cysteines of the CXXCH sequence (Allen et al. 2003; Sanders et al. 2010). This disulfide needs to be reduced to a dithiol to provide free sulfhydryls for the heme ligation. Consistent with this view is the fact that groups of specific oxido-reductases that constitute a transmembrane dithiol-disulfide relay from the cytosol to the periplasmic space have been shown to function as c-type cytochrome assembly factors (Allen et al. 2003; Kadokura et al. 2003; Mapller and Hederstedt 2006; Sanders et al. 2010). The proposal that the components of this pathway control the in vivo redox status of the CXXCH sulfhydryls has been inferred from the presence of motifs in their protein sequences that are consistent with a function in redox chemistry and also from the demonstration that their recombinant forms participate in dithiol–disulfide exchange reactions (Monika et al. 1997; Setterdahl et al. 2000). Moreover, the ability of exogenous thiol compounds to bypass the lack of these factors in vivo substantiates the view that the redox components have a disulfide-reducing activity in the pathway (e.g., Sambongi and Ferguson 1994; Fabianek et al. 1998; Beckett et al. 2000; Deshmukh et al. 2000; Bardischewsky and Friedrich 2001; Erlendsson and Hederstedt 2002; Erlendsson et al. 2003; Feissner et al. 2005; Turkarslan et al. 2008).While the role of these pathways is well established in bacteria, much less is known about the components that catalyze thiol/disulfide chemistry in the mitochondrial intermembrane space (IMS), which is topologically equivalent to the bacterial periplasm. By analogy with the bacterial pathways, the participation of redox-active factors that catalyze thiol formation is expected, as the mitochondrial IMS houses two c-type cytochromes, the soluble cytochrome c and the membrane-bound cytochrome c₁, both of which function in respiration. In fungi, heme attachment to apocytochromes c and c₁ is dependent upon the IMS resident cytochrome c and c₁ heme lyases, CCHL and CC₁HL, although the exact role of these lyases in the assembly process is still unclear (Dumont et al. 1987; Zollner et al. 1992). Conversion of apocytochrome to holocytochrome c depends only on CCHL, while apocytochrome c₁ can be acted upon by both CCHL and CC₁HL (Matner and Sherman 1982; Dumont et al. 1987; Stuart et al. 1990; Zollner et al. 1992; Bernard et al. 2003). In animals, apoforms of cytochromes c and c₁ are assembled by a unique heme lyase, HCCS, which carries both the CCHL and CC₁HL activities (Prakash et al. 2002; Schwarz and Cox 2002; Bernard et al. 2003).Cyc2p, a component first described as a mitochondrial biogenesis factor in yeast (Matner and Sherman 1982; Dumont et al. 1993; Pearce et al. 1998; Sanchez et al. 2001), was recently rediscovered in the context of cytochrome c₁ maturation (Bernard et al. 2003). Cyc2p is located at the mitochondrial inner membrane with its C-terminal domain containing a non-covalently bound FAD exposed to the IMS (Bernard et al. 2005). A redox function for Cyc2p is likely based on the finding that a recombinant form of the molecule exhibits a NAD(P)H-dependent reductase activity (Bernard et al. 2005). However, as Cyc2p activity is not essential for the maturation process, a functional redundancy was postulated based on the fact that a cyc2-null mutant still assembles holoforms of cytochromes c and c₁ (Bernard et al. 2005). The absolute requirement of Cyc2p was revealed via genetic analysis of the cyc2-null cyt1-34 combination that displays a synthetic respiratory-deficient phenotype with loss of holocytochrome c₁ assembly (Bernard et al. 2005). The cyt1-34 mutation maps to the gene encoding cytochrome c₁ and results in a CAPCH heme-binding site replacing the wild-type CAACH site (Bernard et al. 2005). The synthetic interaction is specific for the cyt1-34 allele carrying the A-to-P mutation and is not observed in a cyc2-null cyt1-48 strain carrying an A-to-D mutation at the heme-binding site of apocytochrome c₁ (Bernard et al. 2005). The fact that Cyc2p becomes essential when the cytochrome c₁ heme-binding site carries an A-to-P mutation suggests that the CXXCH motif could be the target of Cyc2p action in vivo. One possible interpretation for this observation is that the P residue alters the reactivity of the cysteinyl thiols to redox chemistry so that the apocytochrome c₁ CAPCH heme-binding site occurs in an oxidized (disulfide) form, which must be reduced in a Cyc2p-dependent reaction before heme attachment can occur.In this article, we have undertaken a genetic approach to elucidate this pathway and searched for suppressors that alleviate the respiratory deficiency of the cyc2-null cyt1-34 strain. Either overexpression of CC₁HL or replacement of the P mutation in the heme-binding site by H, L, S, or T residues restore the assembly of holocytochrome c₁. In vitro measurement of redox potential of apoforms of CA(A/P/T)CH cytochrome c₁ indicates that there is no change in the thermodynamic stability of the disulfide at the CXXCH motif that could account for the Cyc2p-dependent assembly of cytochrome c₁. Genetic studies reveal that the replacement of the second A residue at the CAACH motif by H, L, P, S, and T residues is key in determining the conversion of apocytochrome c₁ to its corresponding holoform via the CCHL and/or CC₁HL-dependent pathway. We also demonstrate that Cyc2p is a component dedicated to the CCHL pathway and is not required for the CC₁HL-dependent assembly of cytochrome c₁. We propose that the CAPCH cytochrome c₁ is strictly dependent upon CCHL and Cyc2p for its assembly but becomes a substrate of CC₁HL upon overexpression of CC₁HL or in the presence of H, L, S, or T mutations.     

Changes in the Nuclear Envelope Environment Affect Spindle Pole Body Duplication in Saccharomyces cerevisiae

The Saccharomyces cerevisiae nuclear membrane is part of a complex nuclear envelope environment also containing chromatin, integral and peripheral membrane proteins, and large structures such as nuclear pore complexes (NPCs) and the spindle pole body. To study how properties of the nuclear membrane affect nuclear envelope processes, we altered the nuclear membrane by deleting the SPO7 gene. We found that spo7Δ cells were sickened by the mutation of genes coding for spindle pole body components and that spo7Δ was synthetically lethal with mutations in the SUN domain gene MPS3. Mps3p is required for spindle pole body duplication and for a variety of other nuclear envelope processes. In spo7Δ cells, the spindle pole body defect of mps3 mutants was exacerbated, suggesting that nuclear membrane composition affects spindle pole body function. The synthetic lethality between spo7Δ and mps3 mutants was suppressed by deletion of specific nucleoporin genes. In fact, these gene deletions bypassed the requirement for Mps3p entirely, suggesting that under certain conditions spindle pole body duplication can occur via an Mps3p-independent pathway. These data point to an antagonistic relationship between nuclear pore complexes and the spindle pole body. We propose a model whereby nuclear pore complexes either compete with the spindle pole body for insertion into the nuclear membrane or affect spindle pole body duplication by altering the nuclear envelope environment.THE nuclear envelope is composed of distinct outer and inner nuclear membranes. The outer nuclear membrane is continuous with the endoplasmic reticulum. The inner nuclear membrane is associated with a unique set of proteins, some of which mediate interactions between the nuclear envelope and chromatin (reviewed in Zhao et al. 2009). Nuclear pore complexes traverse both membranes and allow transport of proteins and solutes between the cytoplasm and the nucleus. The inner and outer nuclear membranes fuse in the region surrounding each nuclear pore complex.In animal cells, the nuclear envelope disassembles as cells enter mitosis and reassembles upon mitotic exit. Nuclear envelope breakdown allows the association of chromosomes with spindle microtubules, which are nucleated from centrosomes that reside in the cytoplasm. In contrast, certain types of fungi, such as the budding yeast Saccharomyces cerevisiae, undergo closed mitosis, where the nuclear envelope remains intact throughout the entire cell cycle. Closed mitosis is possible because the yeast centrosome-equivalent, the spindle pole body (SPB), is embedded in the nuclear envelope, allowing the SPB to nucleate both cytoplasmic and nuclear microtubules.SPB duplication requires a mechanism for inserting the new SPB into the nuclear envelope (reviewed in Jaspersen and Winey 2004). The new SPB begins to form in late G₁/early S phase as satellite material deposited on the cytoplasmic face of an electron-dense region of the nuclear envelope, called the half-bridge. The satellite material matures into a duplication plaque, which is then inserted into the nuclear membrane and becomes the daughter SPB. Many genes are known to be required for SPB duplication, and this process has been carefully examined cytologically (Rose and Fink 1987; Winey et al. 1991, 1993; Spang et al. 1995; Bullitt et al. 1997; Adams and Kilmartin 1999; Elliott et al. 1999; Schramm et al. 2000; Jaspersen et al. 2002; Nishikawa et al. 2003; Araki et al. 2006). However, the exact mechanisms by which SPB duplication and insertion occur remain a mystery.Equally unclear is how nuclear pore complexes are inserted into an intact nuclear envelope (reviewed in Hetzer and Wente 2009). For both the SPB and nuclear pore complexes, the inner and outer nuclear membranes must fuse to allow insertion into the nuclear envelope. Yeast and vertebrate nuclear pore complexes each have four pore membrane (POM) nucleoporins containing transmembrane domains. Other nucleoporins have motifs with potential for bending membranes or sensing membrane curvature. Thus, certain nuclear pore complex components may have the ability to alter the nuclear membrane or stabilize particular membrane conformations (Devos et al. 2004, 2006; Alber et al. 2007; Drin et al. 2007). It is interesting to note that, in S. cerevisiae, nuclear pore complexes are enriched in the vicinity of the SPB (Heath et al. 1995; Winey et al. 1997; Adams and Kilmartin 1999), but the significance of this phenomenon is not known. The SPB and nuclear pore complexes share at least two common components, the integral membrane protein Ndc1p and the small calcium-binding protein Cdc31p (Chial et al. 1998; Fischer et al. 2004). Ndc1p is thought to play a role in insertion of both SPBs and nuclear pore complexes into the nuclear membrane.SUN domain proteins are a conserved family of inner nuclear membrane proteins that interact with specific outer nuclear membrane proteins to form a physical bridge across the nuclear envelope (reviewed in Hiraoka and Dernburg 2009; Razafsky and Hodzic 2009). One of the components of the S. cerevisiae SPB is the SUN domain protein Mps3p. The N terminus of Mps3p is in the nucleoplasm, while the C terminus, containing the SUN domain, is found in the space between the inner and outer nuclear membranes. In addition to the SPB, Mps3p localizes to multiple foci at the nuclear periphery, and these two pools of Mps3p have distinct functions (Jaspersen et al. 2002, 2006; Nishikawa et al. 2003). At the SPB, Mps3p is required for half-bridge formation and early steps of SPB duplication, and cells compromised for Mps3p function accumulate in mitosis with a single SPB and a monopolar spindle (Jaspersen et al. 2002; Nishikawa et al. 2003). At the nuclear periphery, Mps3p is involved in tethering telomeres to the nuclear envelope in mitosis and meiosis, sequestering DNA double-strand breaks away from recombination factors, and associating with soluble chromatin proteins (Antoniacci et al. 2004, 2007; Bupp et al. 2007; Conrad et al. 2007, 2008; Oza et al. 2009; Schober et al. 2009).While many structural features of the yeast nucleus have been identified, little is known about how the physical properties of the nuclear membrane contribute to processes that occur at the nuclear envelope. As noted above, resident proteins of the nuclear envelope may affect nuclear membrane properties. In addition, the nuclear membrane is affected by altering lipid biosynthesis, for example, by inactivating the phosphatidic acid (PA) phosphohydrolase Pah1p or by inactivating the phosphates complex, made of Spo7p and Nem1p, which activates Pah1p. In the absence of Spo7p, Nem1p, or Pah1p, cells exhibit nuclear envelope extensions and extensive ER membrane sheets, and they also have altered membrane lipid composition, including a decrease in phosphatidylcholine and an increase in PA, phosphatidylethanolamine, and phosphatidylinositol (Siniossoglou et al. 1998; Santos-Rosa et al. 2005; Campbell et al. 2006; Han et al. 2006). These three proteins are unique among phospholipid biosynthesis proteins in their ability to affect nuclear morphology upon gene disruption (Han et al. 2008). A similar phenotype was seen upon overexpression of DGK1, which counteracts the activity of Pah1p by converting diacylglycerol to PA, leading to an increase in PA levels at the nuclear envelope (Han et al. 2008). Consistent with a conserved role for Pah1p in regulating nuclear envelope processes, deletion of either NEM1 or SPO7 is synthetically lethal with deletions of certain nucleoporin genes (Siniossoglou et al. 1998), and inactivation of the PAH1 homolog in Caenorhabditis elegans, LPIN-1, results in defects in nuclear envelope disassembly and reassembly (Golden et al. 2009; Gorjanacz and Mattaj 2009).To identify processes that are affected by altered nuclear membrane properties, we screened for pathways that are compromised in spo7Δ cells. We found that SPO7 inactivation strongly influences the SPB. By screening for proteins that could alleviate spo7Δ-induced SPB defects, we uncovered an unexpected inhibitory role for nucleoporins in SPB function, revealing that nuclear pore complexes, or components thereof, act antagonistically to the SPB in the nuclear envelope. Taken together, our findings indicate that the nuclear envelope environment is important for the function of protein complexes and biological processes occurring at the nuclear periphery.     

Telomerase Recruitment in Saccharomyces cerevisiae Is Not Dependent on Tel1-Mediated Phosphorylation of Cdc13

In Saccharomyces cerevisiae, association between the Est1 telomerase subunit and the telomere-binding protein Cdc13 is essential for telomerase to be recruited to its site of action. A current model proposes that Tel1 binding to telomeres marks them for elongation, as the result of phosphorylation of a proposed S/TQ cluster in the telomerase recruitment domain of Cdc13. However, three observations presented here argue against one key aspect of this model. First, the pattern of Cdc13 phosphatase-sensitive isoforms is not altered by loss of Tel1 function or by mutations introduced into two conserved serines (S249 and S255) in the Cdc13 recruitment domain. Second, an interaction between Cdc13 and Est1, as monitored by a two-hybrid assay, is dependent on S255 but Tel1-independent. Finally, a derivative of Cdc13, cdc13–(S/TQ)₁₁→(S/TA)₁₁, in which every potential consensus phosphorylation site for Tel1 has been eliminated, confers nearly wild-type telomere length. These results are inconsistent with a model in which the Cdc13–Est1 interaction is regulated by Tel1-mediated phosphorylation of the Cdc13 telomerase recruitment domain. We propose an alternative model for the role of Tel1 in telomere homeostasis, which is based on the assumption that Tel1 performs the same molecular task at double-strand breaks (DSBs) and chromosome termini.TELOMERE length homeostasis is a highly regulated process that must balance telomere loss (as the result of incomplete replication and/or nucleolytic degradation) with telomeric repeat addition (through the action of telomerase and/or recombination). In the budding yeast Saccharomyces cerevisiae, a key regulatory event is recruitment of telomerase to chromosome ends by the telomere end-binding protein Cdc13 (Nugent et al. 1996; Evans and Lundblad 1999; Pennock et al. 2001; Bianchi et al. 2004; Chan et al. 2008). Recruitment relies on a direct interaction between Cdc13 and the Est1 subunit of telomerase (Pennock et al. 2001), which brings the catalytic core of the enzyme to its site of action. Disruption of this interaction, due to mutations in either CDC13 (cdc13-2) or EST1 (est1-60), results in an Est (ever-shorter-telomere) phenotype, as manifested by progressive telomere shortening and an eventual senescence phenotype. The recruitment activity of Cdc13, which resides in a 15-kDa N-terminal domain (Pennock et al. 2001), is sufficient to direct telomerase even to nontelomeric sites (Bianchi et al. 2004). As predicted by the recruitment model, association of telomerase with telomeres is greatly reduced in strains expressing the recruitment-defective cdc13-2 allele (Chan et al. 2008).Telomerase action at individual telomeres is highly regulated. Using an assay that monitors telomere addition at single nucleotide resolution (single telomere extension, STEX), Lingner and colleagues showed that only ∼7% of telomeres with wild-type (i.e., 300 bp) length are elongated by telomerase during a single cell cycle (Teixeira et al. 2004). However, as telomere length declines, the extension frequency increases: ∼20% of telomeres 200 bp in length and >40% of 100-bp-long telomeres are elongated (Teixeira et al. 2004; Arneric and Lingner 2007). The mechanism by which telomerase distinguishes short from long telomeres has been the subject of intense investigation. Work from a number of laboratories has led to the proposal that Tel1-dependent phosphorylation of Cdc13 at underelongated telomeres mediates the interaction between Cdc13 and the telomerase-associated Est1 protein, thus ensuring that telomerase is directed to the shortest telomeres in a population. In support of this model, the Est1 and Est2 telomerase subunits exhibit enhanced association with telomeres that have been artificially shortened, whereas Cdc13 displays length-independent association with telomeres (Bianchi and Shore 2007; Sabourin et al. 2007). This suggests that the preferential elongation of shorter telomeres is controlled at the level of recruitment of the telomerase holoenzyme by Cdc13. Furthermore, efficient association of Est1 and Est2 with chromosome ends requires Tel1 and Mre11 (which acts in the same pathway as Tel1 for telomere length regulation; Nugent et al. 1998; Ritchie and Petes 2000) but not Mec1 (Takata et al. 2005; Goudsouzian et al. 2006). Tel1 itself is also telomere bound (Takata et al. 2004), with enhanced binding to shorter telomeres (Bianchi and Shore 2007; Hector et al. 2007; Sabourin et al. 2007; Abdallah et al. 2009), although there is considerable controversy over the degree and timing of Tel1 association with chromosome termini during the cell cycle. As expected for a key regulator of the interaction between Cdc13 and a telomerase subunit, a tel1-Δ strain has short telomeres (Lustig and Petes 1986), although telomere length is not impaired enough to confer the Est phenotype displayed by cdc13-2 and est1-60 strains.Implicit in the above proposal is that Cdc13 must be a direct substrate of Tel1. In support of this, Teng and colleagues reported several years ago that the recruitment domain of Cdc13 has a cluster of potential Tel1 (and/or Mec1) phosphorylation sites (Tseng et al. 2006). Substrates of the DNA damage kinases often contain several closely spaced phosphorylation sites, termed S/TQ cluster domains (SCDs), which are considered a structural hallmark of many DNA damage-response proteins (Traven and Heierhorst 2005). On the basis of in vitro kinase assays with GST fusions to 75- to 90-amino-acid portions of the Cdc13 recruitment domain, Tseng et al. 2006 concluded that four SQ sites in the recruitment domain of Cdc13 are overlapping substrates for both Tel1 and Mec1, leading to the proposal that telomerase recruitment in S. cerevisiae is regulated by Tel1-dependent phosphorylation of Cdc13.The above model makes a key prediction: in a tel1-Δ strain, telomerase should no longer exhibit a length-dependent pattern of elongation. However, preferential elongation of short telomeres still occurs at native chromosome ends in the absence of Tel1 (Arneric and Lingner 2007). In addition, Petes and colleagues have argued, on the basis of epistasis data, that Tel1 performs an indirect role in the telomerase pathway, rather than directly targeting a telomerase regulator (Ritchie et al. 1999; Ritchie and Petes 2000). These observations are not easily explained, if preferential recognition of short telomeres by telomerase is mediated by Tel1-dependent phosphorylation of Cdc13. In this current study, we have re-examined the evidence for phosphorylation of Cdc13 as a regulatory mechanism for telomere length homeostasis. We report on a series of observations that indicate that Tel1 contributes to telomere length control through a mechanism other than phosphorylation of the Cdc13 S/TQ cluster.     

HeurAA: Accurate and Fast Detection of Genetic Variations with a Novel Heuristic Amplicon Aligner Program for Next Generation Sequencing

Next generation sequencing (NGS) of PCR amplicons is a standard approach to detect genetic variations in personalized medicine such as cancer diagnostics. Computer programs used in the NGS community often miss insertions and deletions (indels) that constitute a large part of known human mutations. We have developed HeurAA, an open source, heuristic amplicon aligner program. We tested the program on simulated datasets as well as experimental data from multiplex sequencing of 40 amplicons in 12 oncogenes collected on a 454 Genome Sequencer from lung cancer cell lines. We found that HeurAA can accurately detect all indels, and is more than an order of magnitude faster than previous programs. HeurAA can compare reads and reference sequences up to several thousand base pairs in length, and it can evaluate data from complex mixtures containing reads of different gene-segments from different samples. HeurAA is written in C and Perl for Linux operating systems, the code and the documentation are available for research applications at http://sourceforge.net/projects/heuraa/     

Wham: Identifying Structural Variants of Biological Consequence

Existing methods for identifying structural variants (SVs) from short read datasets are inaccurate. This complicates disease-gene identification and efforts to understand the consequences of genetic variation. In response, we have created Wham (Whole-genome Alignment Metrics) to provide a single, integrated framework for both structural variant calling and association testing, thereby bypassing many of the difficulties that currently frustrate attempts to employ SVs in association testing. Here we describe Wham, benchmark it against three other widely used SV identification tools–Lumpy, Delly and SoftSearch–and demonstrate Wham’s ability to identify and associate SVs with phenotypes using data from humans, domestic pigeons, and vaccinia virus. Wham and all associated software are covered under the MIT License and can be freely downloaded from github (https://github.com/zeeev/wham), with documentation on a wiki (http://zeeev.github.io/wham/). For community support please post questions to https://www.biostars.org/.

This is PLOS Computational Biology software paper.

     

MetAMOS: a modular and open source metagenomic assembly and analysis pipeline

We describe MetAMOS, an open source and modular metagenomic assembly and analysis pipeline. MetAMOS represents an important step towards fully automated metagenomic analysis, starting with next-generation sequencing reads and producing genomic scaffolds, open-reading frames and taxonomic or functional annotations. MetAMOS can aid in reducing assembly errors, commonly encountered when assembling metagenomic samples, and improves taxonomic assignment accuracy while also reducing computational cost. MetAMOS can be downloaded from: https://github.com/treangen/MetAMOS.     

A Genetic Survey of Fluoxetine Action on Synaptic Transmission in Caenorhabditis elegans

Fluoxetine is one of the most commonly prescribed medications for many behavioral and neurological disorders. Fluoxetine acts primarily as an inhibitor of the serotonin reuptake transporter (SERT) to block the removal of serotonin from the synaptic cleft, thereby enhancing serotonin signals. While the effects of fluoxetine on behavior are firmly established, debate is ongoing whether inhibition of serotonin reuptake is a sufficient explanation for its therapeutic action. Here, we provide evidence of two additional aspects of fluoxetine action through genetic analyses in Caenorhabditis elegans. We show that fluoxetine treatment and null mutation in the sole SERT gene mod-5 eliminate serotonin in specific neurons. These neurons do not synthesize serotonin but import extracellular serotonin via MOD-5/SERT. Furthermore, we show that fluoxetine acts independently of MOD-5/SERT to regulate discrete properties of acetylcholine (Ach), gamma-aminobutyric acid (GABA), and glutamate neurotransmission in the locomotory circuit. We identified that two G-protein–coupled 5-HT receptors, SER-7 and SER-5, antagonistically regulate the effects of fluoxetine and that fluoxetine binds to SER-7. Epistatic analyses suggest that SER-7 and SER-5 act upstream of AMPA receptor GLR-1 signaling. Our work provides genetic evidence that fluoxetine may influence neuronal functions and behavior by directly targeting serotonin receptors.FLUOXETINE is a selective serotonin reuptake inhibitor (SSRI) and has made a major impact on the treatment of many behavioral disorders. The empirical action of SSRIs is blocking the serotonin reuptake transporter (SERT). SERT is localized in the plasma membrane and transports extracellular serotonin (5-HT) into the cytoplasm (Blakely et al. 1991; Hoffman et al. 1991), this being the major mechanism of terminating 5-HT signaling. Consequently, SSRIs are thought to exert therapeutic effects by blocking SERT from removal of 5-HT in the synaptic clef, thereby increasing the level of 5-HT signals (Schatzberg and Nemeroff 2004). However, several observations point to additional actions of SSRIs on the 5-HT system and neuronal functions. First, knockout of SERT in mouse caused a marked reduction of 5-HT in the brain (Bengel et al. 1998). Second, a variety of studies with cultured mammalian cells and mouse brain slices showed that SSRIs and tricyclic antidepressant agents (TCAs) have high affinities to many 5-HT receptor subtypes and act as agonists or antagonists depending on particular receptors being tested (Ni and Miledi 1997; Kroeze and Roth 1998; Eisensamer et al. 2003). Third, genetic analyses of the nematode Caenorhabditis elegans in our laboratory and others showed that fluoxetine and the TCA imipramine (Tofrani) could influence behavior independent of SERT function (Weinshenker et al. 1995; Ranganathan et al. 2001; Dempsey et al. 2005). In this study, we carried out a systematic survey of SSRIs treatment in C. elegans to gain new insights into actions of SSRIs on the 5-HT system and other neurotransmitter systems.In both vertebrates and invertebrates, 5-HT functions as a neuromodulator to either facilitate or inhibit synaptic transmission of other neurotransmitters (Fink and Gothert 2007). Modulation of synaptic activity by 5-HT signaling underscores the synaptic plasticity involved in stress responses, learning, adaptation, and memory (Kandel 2001; Zhang et al. 2005). The role of 5-HT in C. elegans was initially identified through pharmacological experiments showing that exogenous 5-HT can promptly induce changes in a variety of behaviors, including feeding, egg laying, and locomotion (Avery and Horvitz 1990; Weinshenker et al. 1995; Nurrish et al. 1999). The relevance of these behaviors to endogenous 5-HT has since been validated through studies of mutants of 5-HT signaling. Importantly, multiple 5-HT receptors may function in distinct cells synergistically or antagonistically to regulate a specific behavior (Carnell et al. 2005; Dernovici et al. 2007; Murakami and Murakami 2007; Hapiak et al. 2009). In nearly all tested paradigms, fluoxetine and imipramine induce behavioral changes similarly to exogenous 5-HT (Weinshenker et al. 1995; Nurrish et al. 1999), implying that fluoxetine regulates 5-HT inputs to these neural circuits. However, the tryptophan hydroxylase gene tph-1 is required for 5-HT biosynthesis in C. elegans (Sze et al. 2000), mod-5 encodes its sole SERT (Ranganathan et al. 2001), and yet fluoxetine could stimulate egg laying and inhibit locomotion in mod-5 and tph-1 mutants (Weinshenker et al. 1995; Choy and Thomas 1999; Ranganathan et al. 2001; Dempsey et al. 2005). These findings provided a basis for further investigation into genes and synaptic functions regulated by 5-HT and the impact of fluoxetine on 5-HT signaling.Here we present genetic evidence of multifaceted effects of fluoxetine on the 5-HT system and its downstream targets in C. elegans. We show that fluoxetine treatment and loss of MOD-5/SERT function do not simply increase presynaptic 5-HT signals. Rather, they may eliminate 5-HT in specific neurons. Furthermore, fluoxetine acts independently of SERT to regulate 5-HT serotonin receptors and their downstream targets involved in acetylcholine (ACh), gamma-aminobutyric acid (GABA), and glutamate neurotransmission.