Search

BackgroundMost organisms rely on a molecular circadian clock to orchestrate a wide range of physiological processes to match the 24-hour day. These molecular clocks are typically based on a negative feedback loop among a small set of proteins that govern the circadian output. Light or other environmental conditions can reset the circadian clock, but true circadian behaviors continue to cycle even in constant darkness, with an intrinsic period called the free-running period (FRP). Spiders have unusual FRPs, with some species having extremely short FRP (e.g. 18 hours for trashline orb weaver), and many having highly variable FRPs (intraspecific variation of up to 10 hours). In the absence of any genetic model of circadian rhythms in spiders, we developed a mathematical model to optimize experimental conditions for identifying circadian genes that also respond to light cues. ResultsOur mathematical model involved a single gene that encodes a protein that inhibits its own transcription. In our model, light degrades the circadian transcript, which allows a broad range of FRPs to be entrained to a 24-hour day. Our model predicted that exposing spiders to a pulse of light in the middle of the night would cause a pattern of expression between two later time points that was opposite the pattern exhibited by spiders who did not receive a pulse of light. RNA-sequencing of four groups of adult female orb weaving spiders, Metazygia wittfeldae, under these experimental conditions resulted in 528 significantly differentially expressed (DE) transcripts between the two collection times or between the light pulse and no light pulse. Consistent with our model, we found a cluster of transcripts with the flipped pattern of expression between the two collection times, dependent on the application of light.ConclusionsOur DE transcripts represent the first genetic evidence for circadian output in spiders. Furthermore, those transcripts with a flipped pattern of expression represent prime candidates for light-sensitive circadian genes, which may be involved in entraining the circadian clock to light. Functions of these genes varied from growth and development to reproduction to gene regulation, consistent with other circadian systems.

Background Spiders have evolved two types of sticky capture threads: one with wet adhesive spun by ecribellate orb-weavers and another with dry adhesive spun by cribellate spiders. The evolutionary history of cribellate capture threads is especially poorly understood. Here, we use genomic approaches to catalog the spider-specific silk gene family (spidroins) for the cribellate orb-weaver Uloborus diversus. Results We show that the cribellar spidroin, which forms the puffy fibrils of cribellate threads, has three distinct repeat units, one of which is conserved across cribellate taxa separated by ~ 250 Mya. We also propose candidates for a new silk type, paracribellar spidroins, which connect the puffy fibrils to pseudoflagelliform support lines. Moreover, we describe the complete repeat architecture for the pseudoflagelliform spidroin (Pflag), which contributes to extensibility of pseudoflagelliform axial fibers. Conclusions Our finding that Pflag is closely related to Flag, supports homology of the support lines of cribellate and ecribellate capture threads. It further suggests an evolutionary phase following gene duplication, in which both Flag and Pflag were incorporated into the axial lines, with subsequent loss of Flag in uloborids, and increase in expression of Flag in ecribellate orb-weavers, explaining the distinct mechanical properties of the axial lines of these two groups.

The origin of aggregate silk glands and their production of wet adhesive silks is considered a key innovation of the Araneoidea, a superfamily of spiders that build orb-webs and cobwebs. Orb-web weavers place aggregate glue on an extensible capture spiral, whereas cobweb weavers add it to the ends of strong, stiff fibers, called gumfoot lines. Here we describe the material behavior and quantitative proteomics of the aggregate glues of two cobweb weaving species, the Western black widow, Latrodectus hesperus, and the common house spider, Parasteatoda tepidariorum. For each species respectively, we identified 48 and 33 proteins that were significantly more abundant in the portion of the gumfoot line with glue than in its fibers. These proteins were more highly glycosylated and phosphorylated than proteins found in silk fibers without glue, which likely explains aggregate glue stickiness. Most glue-enriched proteins were of anterior aggregate gland origin, supporting the hypothesis that cobweb weavers’ posterior aggregate glue is specialized for another function. We found that cobweb weaver glue droplets are stiffer and tougher than the adhesive of most orb-web weaving species. Attributes of gumfoot glue protein composition that likely contribute to this stiffness include the presence of multiple protein families with conserved cysteine residues, a bimodal distribution of isoelectric points, and families with conserved functions in protein aggregation, all of which should contribute to cohesive protein-protein interactions. House spider aggregate droplets were more responsive to humidity changes than black widow droplets, which could be mediated by differences in protein sequence, post-translational modifications, the non-protein components of the glue droplets, and/or the larger amount of aqueous material that surrounds the adhesive cores of their glue droplets.

Molecular studies of the secretory glands involved in spider silk production have revealed candidate genes for silk synthesis and a complicated history of spider silk gene evolution. However, differential gene expression profiles of the multiple silk gland types within an individual orb-web weaving spider are lacking. Each of these gland types produces a functionally distinct silk type. Comparison of gene expression among spider silk gland types would provide insight into the genes that define silk glands generally from non-silk gland tissues, and the genes that define silk glands from each other. Here, we perform 3' tag digital gene expression profiling of the seven silk gland types of the silver garden orb weaver Argiope argentata. Five of these gland types produce silks that are non-adhesive fibers, one silk includes both fibers and glue-like adhesives, and one silk is exclusively glue-like. We identify 1275 highly expressed, significantly upregulated, and tissue specific silk gland specific transcripts (SSTs). These SSTs include seven types of spider silk protein encoding genes known as spidroin genes. We find that the fiber-producing major ampullate and minor ampullate silk glands have more similar expression profiles than any other pair of glands. We also find that a subset of the SSTs is enriched for transmembrane transport and oxidoreductases, and that these transcripts highlight differences and similarities among the major ampullate, minor ampullate, and aggregate silk glands. Furthermore, we show that the wet glue-producing aggregate glands have the most unique SSTs, but still share some SSTs with fiber producing glands. Aciniform glands were the only gland type to share a majority of SSTs with other silk gland types, supporting previous hypotheses that duplication of aciniform glands and subsequent divergence of the duplicates gave rise to the multiple silk gland types within an individual spider.

The human microbiome has been proposed as a tool to investigate different forensic questions, including for the identification of multiple personal information. However, the fragmented state of the publicly available data has retarded the development of analysis techniques and, therefore, the implementation of microbiomes as a forensic tool. To address this, we introduce the forensic microbiome database (FMD), which is a collection of 16S rRNA data and associated metadata generated from publicly available data. The raw data was further normalized and processed using a pipeline to create a standardized data set for downstream analysis. We present a website allowing for the exploration of geolocation signals in the FMD. The website allows users to investigate the taxonomic differences between microbiomes harvested from different locations and to predict the geolocation of their data based on the FMD sequences. All the results are presented in dynamic graphics to allow for a rapid and intuitive investigation of the taxonomic distributions underpinning the geolocation signals and prediction between locations. Apart from the forensic aspect, the database also allows exploration and comparison of microbiome samples from different geolocation and between different body sites. The goal of the FMD is to provide the scientific and non-scientific communities with data and tools to explore the possibilities of microbiomes to answer forensic questions and serve as a model for any future such databases.1

To survive in terrestrial and aquatic environments, spiders often rely heavily on their silk. The vast majority of silks that have been studied are from orb-web or cob-web weaving species, leaving the silks of water-associated spiders largely undescribed. We characterize transcripts, proteins, and silk fibres from the semi-aquatic spider Dolomedes triton. From silk gland RNAseq libraries, we report 18 silk transcripts representing four categories of known silk protein types: aciniform, ampullate, pyriform, and tubuliform. Proteomic and structural analyses (scanning electron microscopy, energy dispersive X-ray spectrometry, contact angle) of the D. triton submersible egg sac reveal similarities to silks from aquatic caddisfly larvae. We identified two layers in D. triton egg sacs, notably a highly hydrophobic outer layer with a different elemental composition compared to egg sacs of terrestrial spiders. These features may provide D. triton egg sacs with their water repellent properties.

Chronic liver disease due to alcohol-use disorder contributes markedly to the global burden of disease and mortality1–3. Alcoholic hepatitis is a severe and life-threatening form of alcohol-associated liver disease. The gut microbiota promotes ethanol-induced liver disease in mice4, but little is known about the microbial factors that are responsible for this process. Here we identify cytolysin—a two-subunit exotoxin that is secreted by Enterococcus faecalis5,6—as a cause of hepatocyte death and liver injury. Compared with non-alcoholic individuals or patients with alcohol-use disorder, patients with alcoholic hepatitis have increased faecal numbers of E. faecalis. The presence of cytolysin-positive (cytolytic) E. faecalis correlated with the severity of liver disease and with mortality in patients with alcoholic hepatitis. Using humanized mice that were colonized with bacteria from the faeces of patients with alcoholic hepatitis, we investigated the therapeutic effects of bacteriophages that target cytolytic E. faecalis. We found that these bacteriophages decrease cytolysin in the liver and abolish ethanol-induced liver disease in humanized mice. Our findings link cytolytic E. faecalis with more severe clinical outcomes and increased mortality in patients with alcoholic hepatitis. We show that bacteriophages can specifically target cytolytic E. faecalis, which provides a method for precisely editing the intestinal microbiota. A clinical trial with a larger cohort is required to validate the relevance of our findings in humans, and to test whether this therapeutic approach is effective for patients with alcoholic hepatitis. In patients with alcoholic hepatitis, cytolysin-positive Enterococcus faecalis strains are correlated with liver disease severity and increased mortality, and in mouse models these strains can be specifically targeted by bacteriophages.

Background Bacterial pan-genomes, comprised of conserved and variable genes across multiple sequenced bacterial genomes, allow for identification of genomic regions that are phylogenetically discriminating or functionally important. Pan-genomes consist of large amounts of data, which can restrict researchers ability to locate and analyze these regions. Multiple software packages are available to visualize pan-genomes, but currently their ability to address these concerns are limited by using only pre-computed data sets, prioritizing core over variable gene clusters, or by not accounting for pan-chromosome positioning in the viewer. Results We introduce PanACEA (Pan-genome Atlas with Chromosome Explorer and Analyzer), which utilizes locally-computed interactive web-pages to view ordered pan-genome data. It consists of multi-tiered, hierarchical display pages that extend from pan-chromosomes to both core and variable regions to single genes. Regions and genes are functionally annotated to allow for rapid searching and visual identification of regions of interest with the option that user-supplied genomic phylogenies and metadata can be incorporated. PanACEA’s memory and time requirements are within the capacities of standard laptops. The capability of PanACEA as a research tool is demonstrated by highlighting a variable region important in differentiating strains of Enterobacter hormaechei. Conclusions PanACEA can rapidly translate the results of pan-chromosome programs into an intuitive and interactive visual representation. It will empower researchers to visually explore and identify regions of the pan-chromosome that are most biologically interesting, and to obtain publication quality images of these regions. Electronic supplementary material The online version of this article (10.1186/s12859-018-2250-y) contains supplementary material, which is available to authorized users.

The skin is a complex living ecosystem harboring diverse microbial communities. Its highly variable properties and influence of intrinsic and extrinsic factors creates unique microenvironments where niche-specific microbes thrive. As part of the skin, hair supports its own microbial habitat that is also intra and inter-personal variable. This little explored substrate has significant potential in forensics microbiome research due to the unique signatures that are available on an individual. To further investigate this, we explored the hair microbiota from scalp and pubic regions in healthy adults to investigate how the hair shaft microenvironment varies microbially. Our results suggest that there are distinct differences between the microbial communities identified on hair shafts originating from different parts of the body. The taxonomic composition of the communities from different hair sources are most reminiscent of those identified from their associated cutaneous region. We further demonstrate that the hair microbiota varies by geographical origin and has the potential to be used to predict the source location of the hair.

Motivation The vast number of available sequenced bacterial genomes occasionally exceeds the facilities of comparative genomic methods or is dominated by a single outbreak strain, and thus a diverse and representative subset is required. Generation of the reduced subset currently requires a priori supervised clustering and sequence-only selection of medoid genomic sequences, independent of any additional genome metrics or strain attributes. Results The Gaussian Genome Representative Selector with Prioritization (GGRaSP) R-package described below generates a reduced subset of genomes that prioritizes maintaining genomes of interest to the user as well as minimizing the loss of genetic variation. The package also allows for unsupervised clustering by modeling the genomic relationships using a Gaussian mixture model to select an appropriate cluster threshold. We demonstrate the capabilities of GGRaSP by generating a reduced list of 315 genomes from a genomic dataset of 4600 Escherichia coli genomes, prioritizing selection by type strain and by genome completeness. Availability and implementaion GGRaSP is available at https://github.com/JCVenterInstitute/ggrasp/. Supplementary information Supplementary data are available at Bioinformatics online.

Objectives To investigate the genomic context of a novel resistance island (RI) in multiply antibiotic-resistant Acinetobacter baumannii clinical isolates and global isolates. Methods Using a combination of long and short reads generated from the Oxford Nanopore and Illumina platforms, contiguous chromosomes and plasmid sequences were determined. BLAST-based analysis was used to identify the RI insertion target. Results Genomes of four multiply antibiotic-resistant A. baumannii clinical strains, from a US hospital system, belonging to prevalent MLST ST2 (Pasteur scheme) and ST281 (Oxford scheme) clade F isolates were sequenced to completion. A class 1 integron carrying aadB (tobramycin resistance) and aadA2 (streptomycin/spectinomycin resistance) was identified. The class 1 integron was 6.8 kb, bounded by IS26 at both ends, and embedded in a new target location between an α/β-hydrolase and a reductase. Due to its novel insertion site and unique RI composition, we suggest naming this novel RI AbGRI4. Molecular analysis of global A. baumannii isolates identified multiple AbGRI4 RI variants in non-ST2 clonal lineages, including variations in the resistance gene cassettes, integron backbone and insertion breakpoints at the hydrolase gene. Conclusions A novel RI insertion target harbouring a class 1 integron was identified in a subgroup of ST2/ST281 clinical isolates. Variants of the RI suggested evolution and horizontal transfer of the RI across clonal lineages. Long- and short-read hybrid assembly technology completely resolved the genomic context of IS-bounded RIs, which was not possible using short reads alone.

Spiders are commonly found in terrestrial environments and many rely heavily on their silks for fitness related tasks such as reproduction and dispersal. Although rare, a few species occupy aquatic or semi-aquatic habitats and for them, silk-related specializations are also essential to survive in aquatic environments. Most spider silks studied to date are from cob-web and orb-web weaving species, leaving the silks from many other terrestrial spiders as well as water-associated spiders largely undescribed. Here, we characterize silks from three Dictynoidea species: the aquatic spiders Argyroneta aquatica and Desis marina as well as the terrestrial Badumna longinqua. From silk gland RNA-Seq libraries, we report a total of 47 different homologs of the spidroin (spider fibroin) gene family. Some of these 47 spidroins correspond to known spidroin types (aciniform, ampullate, cribellar, pyriform, and tubuliform), while other spidroins represent novel branches of the spidroin gene family. We also report a hydrophobic amino acid motif (GV) that, to date, is found only in the spidroins of aquatic and semi-aquatic spiders. Comparison of spider silk sequences to the silks from other water-associated arthropods, shows that there is a diversity of strategies to function in aquatic environments.

Successful treatment of Acinetobacter baumannii infections require early and appropriate antimicrobial therapy. One of the first steps in this process is understanding which β-lactamase (bla) alleles are present and in what combinations. Thus, we performed WGS on 98 carbapenem-resistant A. baumannii (CR Ab). In most isolates, an acquired blaOXA carbapenemase was found in addition to the intrinsic blaOXA allele. The most commonly found allele was blaOXA-23 (n = 78/98). In some isolates, blaOXA-23 was found in addition to other carbapenemase alleles: blaOXA-82 (n = 12/78), blaOXA-72 (n = 2/78) and blaOXA-24/40 (n = 1/78). Surprisingly, 20% of isolates carried carbapenemases not routinely assayed for by rapid molecular diagnostic platforms, i.e., blaOXA-82 and blaOXA-172; all had ISAba1 elements. In 8 CR Ab, blaOXA-82 or blaOXA-172 was the only carbapenemase. Both blaOXA-24/40 and its variant blaOXA-72 were each found in 6/98 isolates. The most prevalent ADC variants were blaADC-30 (21%), blaADC-162 (21%), and blaADC-212 (26%). Complete combinations are reported.

Gene duplication and positive selection can be important determinants of the evolution of venom, a protein-rich secretion used in prey capture and defense. In a typical model of venom evolution, gene duplicates switch to venom gland expression and change function under the action of positive selection, which together with further duplication produces large gene families encoding diverse toxins. Although these processes have been demonstrated for individual toxin families, high-throughput multitissue sequencing of closely related venomous species can provide insights into evolutionary dynamics at the scale of the entire venom gland transcriptome. By assembling and analyzing multitissue transcriptomes from the Western black widow spider and two closely related species with distinct venom toxicity phenotypes, we do not find that gene duplication and duplicate retention is greater in gene families with venom gland biased expression in comparison with broadly expressed families. Positive selection has acted on some venom toxin families, but does not appear to be in excess for families with venom gland biased expression. Moreover, we find 309 distinct gene families that have single transcripts with venom gland biased expression, suggesting that the switching of genes to venom gland expression in numerous unrelated gene families has been a dominant mode of evolution. We also find ample variation in protein sequences of venom gland–specific transcripts, lineage-specific family sizes, and ortholog expression among species. This variation might contribute to the variable venom toxicity of these species.

Background: There are a large and growing number of non-tuberculous mycobacteria (NTM) species that have been isolated, identified, and described in the literature, yet there are many clinical isolates which are not assignable to known species even when the genome has been sequenced. Additionally, a recent manuscript has proposed the reclassification of the Mycobacterium genus into five distinct genera. Methods: We describe using a community standard fast average nucleotide identity (ANI) approximation method, MASH, for classifying NTM genomes by comparison to a resource of type strain genomes and proxy genomes. We evaluate the genus reclassification proposal in light of our ANI, MLST, and pan-genome work. Results: We describe here a sequencing study of hundreds of clinical NTM isolates. To aid in characterizing these isolates we defined a multi-locus sequence typing (MLST) schema for NTMs which can differentiate strains at the species and subspecies level using eight ribosomal protein genes. We determined and deposited the allele profiles for 2,802 NTM and Mycobacterium tuberculosis complex strains in PubMLST. Conclusions: The MLST schema and our pan-genome analysis of Mycobacteria can help inform the design of marker-gene diagnostics. The ANI comparisons likewise can assist in the classification of unknown genomes, even from previously unknown species.

A powerful system for studying protein aggregation, particularly rapid self-assembly, is spider silk. Spider silks are proteinaceous and silk proteins are synthesized and stored within silk glands as liquid dope. As needed, liquid dope is near-instantaneously transformed into solid fibers or viscous adhesives. The dominant constituents of silks are spidroins (spider fibroins) and their terminal domains are vital for the tight control of silk self-assembly. To better understand spidroin termini, we used target capture and deep sequencing to identify spidroin gene sequences from six species representing the araneoid families of Araneidae, Nephilidae, and Theridiidae. We obtained 145 terminal regions, of which 103 are newly annotated here, as well as novel variants within nine diverse spidroin types. Our comparative analyses demonstrated the conservation of acidic, basic, and cysteine amino acid residues across spidroin types that had been proposed to be important for monomer stability, dimer formation, and self-assembly from a limited sampling of spidroins. Computational, protein homology modeling revealed areas of spidroin terminal regions that are highly conserved in three-dimensions despite sequence divergence across spidroin types. Analyses of our dense sampling of terminal regions suggest that most spidroins share stabilization mechanisms, dimer formation, and tertiary structure, despite producing functionally distinct materials.

The evolution of specialized tissues with novel functions, such as the silk synthesizing glands in spiders, is likely an influential driver of adaptive success. Large-scale gene duplication events and subsequent paralog divergence are thought to be required for generating evolutionary novelty. Such an event has been proposed for spiders, but not tested. We de novo assembled transcriptomes from three cobweb weaving spider species. Based on phylogenetic analyses of gene families with representatives from each of the three species, we found numerous duplication events indicative of a whole genome or segmental duplication. We estimated the age of the gene duplications relative to several speciation events within spiders and arachnids and found that the duplications likely occurred after the divergence of scorpions (order Scorpionida) and spiders (order Araneae), but before the divergence of the spider suborders Mygalomorphae and Araneomorphae, near the evolutionary origin of spider silk glands. Transcripts that are expressed exclusively or primarily within black widow silk glands are more likely to have a paralog descended from the ancient duplication event and have elevated amino acid replacement rates compared with other transcripts. Thus, an ancient large-scale gene duplication event within the spider lineage was likely an important source of molecular novelty during the evolution of silk gland-specific expression. This duplication event may have provided genetic material for subsequent silk gland diversification in the true spiders (Araneomorphae).

Spider silk synthesis is an emerging model for the evolution of tissue-specific gene expression and the role of gene duplication in functional novelty, but its potential has not been fully realized. Accordingly, we quantified transcript (mRNA) abundance in seven silk gland types and three non-silk gland tissues for three cobweb-weaving spider species. Evolutionary analyses based on expression levels of thousands of homologous transcripts and phylogenetic reconstruction of 605 gene families demonstrated conservation of expression for each gland type among species. Despite serial homology of all silk glands, the expression profiles of the glue-forming aggregate glands were divergent from fiber-forming glands. Also surprising was our finding that shifts in gene expression among silk gland types were not necessarily coupled with gene duplication, even though silk-specific genes belong to multi-paralog gene families. Our results challenge widely accepted models of tissue specialization and significantly advance efforts to replicate silk-based high-performance biomaterials.

et al., [Background]: The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum. [Results]: We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication. [Conclusions]: Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes., This work was supported by NIH grant NHGRI U54 HG003273 to RAG, the National Science Foundation (IOS-0951886 to NAA and DEB-1257053 to JHW), a Leverhulme visiting fellowship for EES (VF-2012-016), funding and PhD studentships (DJL, LG and AS) from Oxford Brookes University, and a Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) scholarship to CLBP. N-MP was funded by the Deutsche Forschungsgemeinschaft (grant numbers PR 1109/ 4-1, PR 1109/7-1 and PR 1109/6-1 to N-MP). Additional financial backing has been received from the Göttingen Graduate School for Neurosciences, Biophysics and Molecular Biosciences (GGNB), the Göttingen Center for Molecular Biosciences (GZMB), and the University of Göttingen (GAU). NT is supported by a Christiane Nüsslein-Volhard-Foundation fellowship and a “Women in Science” Award by L'Oréal Deutschland and the Deutsche UNESCO-Kommission. NP has been funded by the Volkswagen Foundation (project number: 85 983) and the Emmy Noether Programme of the Deutsche Forschungsgemeinschaft (grant number: PO 1648/3-1). Funding to RJ was provided by the Swedish Research Council VR grant 621-2011-4703.

The introduction of DNA fingerprinting to forensic science rapidly expanded the available evidence that could be garnered from a crime scene and used in court cases. Next generation sequencing technologies increased available genetic data that could be used as evidence by orders of magnitude, and as such, significant additional genetic information is now available for use in forensic science. This includes DNA from the bacteria that live in and on humans, known as the human microbiome. Next generation sequencing of the human microbiome demonstrates that its bacterial DNA can be used to uniquely identify an individual, provide information about their life and behavioral patterns, determine the body site where a sample came from, and estimate postmortem intervals. Bacterial samples from the environment and objects can also be leveraged to address similar questions about the individual(s) who interacted with them. However, the applications of this new field in forensic sciences raises concerns on current methods used in sample processing, including sample collection, storage, and the statistical power of published studies. These areas of human microbiome research need to be fully addressed before microbiome data can become a regularly incorporated evidence type and routine procedure of the forensic toolkit. Here, we summarize information on the current status of microbiome research as applies to the forensic field, the mathematical models used to make predictions, and the possible legal and practical difficulties that can limit the application of microbiomes in forensic science.

Summary The JCVI pan-genome pipeline is a collection of programs to run PanOCT and tools that support and extend the capabilities of PanOCT. PanOCT (pan-genome ortholog clustering tool) is a tool for pan-genome analysis of closely related prokaryotic species or strains. The JCVI Pan-Genome Pipeline wrapper invokes command-line utilities that prepare input genomes, invoke third-party tools such as NCBI Blast+, run PanOCT, generate a consensus pan-genome, annotate features of the pan-genome, detect sets of genes of interest such as antimicrobial resistance (AMR) genes and generate figures, tables and html pages to visualize the results. The pipeline can run in a hierarchical mode, lowering the RAM and compute resources used. Availability and implementation Source code, demo data, and detailed documentation are freely available at https://github.com/JCVenterInstitute/PanGenomePipeline.

Spiders (order Araneae) rely on their silks for essential tasks, such as dispersal, prey capture, and reproduction. Spider silks are largely composed of spidroins, members of a protein family that are synthesized in silk glands. As needed, silk stored in silk glands is extruded through spigots on the spinnerets. Nearly all studies of spider silks have been conducted on females; thus, little is known about male silk biology. To shed light on silk use by males, we compared silk gene expression profiles of mature males to those of females from three cob-web weaving species (Theridiidae). We de novo assembled species-specific male transcriptomes from Latrodectus hesperus, Latrodectus geometricus, and Steatoda grossa followed by differential gene expression analyses. Consistent with their complement of silk spigots, male theridiid spiders express appreciable amounts of aciniform, major ampullate, minor ampullate, and pyriform spidroin genes but not tubuliform spidroin genes. The relative expression levels of particular spidroin genes varied between sexes and species. Because mature males desert their prey-capture webs and become cursorial in their search for mates, we anticipated that major ampullate (dragline) spidroin genes would be the silk genes most highly expressed by males. Indeed, major ampullate spidroin genes had the highest expression in S. grossa males. However, minor ampullate spidroin genes were the most highly expressed spidroin genes in L. geometricus and L. hesperus males. Our expression profiling results suggest species-specific adaptive divergence of silk use by male theridiids.

Background Orb-web weaving spiders and their relatives use multiple types of task-specific silks. The majority of spider silk studies have focused on the ultra-tough dragline silk synthesized in major ampullate glands, but other silk types have impressive material properties. For instance, minor ampullate silks of orb-web weaving spiders are as tough as draglines, due to their higher extensibility despite lower strength. Differences in material properties between silk types result from differences in their component proteins, particularly members of the spidroin (spider fibroin) gene family. However, the extent to which variation in material properties within a single silk type can be explained by variation in spidroin sequences is unknown. Here, we compare the minor ampullate spidroins (MiSp) of orb-weavers and cobweb weavers. Orb-web weavers use minor ampullate silk to form the auxiliary spiral of the orb-web while cobweb weavers use it to wrap prey, suggesting that selection pressures on minor ampullate spidroins (MiSp) may differ between the two groups. Results We report complete or nearly complete MiSp sequences from five cobweb weaving spider species and measure material properties of minor ampullate silks in a subset of these species. We also compare MiSp sequences and silk properties of our cobweb weavers to published data for orb-web weavers. We demonstrate that all our cobweb weavers possess multiple MiSp loci and that one locus is more highly expressed in at least two species. We also find that the proportion of β-spiral-forming amino acid motifs in MiSp positively correlates with minor ampullate silk extensibility across orb-web and cobweb weavers. Conclusions MiSp sequences vary dramatically within and among spider species, and have likely been subject to multiple rounds of gene duplication and concerted evolution, which have contributed to the diverse material properties of minor ampullate silks. Our sequences also provide templates for recombinant silk proteins with tailored properties. Electronic supplementary material The online version of this article (doi:10.1186/s12862-017-0927-x) contains supplementary material, which is available to authorized users.

Orb-web weaving spiders produce a variety of task-specific silks from specialized silk glands. The genetics underlying the synthesis of specific silk types are largely unknown, and transcriptome analysis could be a powerful approach for identifying candidate genes. However, de novo assembly and expression profiling of silk glands with RNA-sequencing (RNAseq) are problematic because the few known gene transcripts for silk proteins are extremely long and highly repetitive. To identify candidate genes for tubuliform (egg case) silk synthesis by the orb-weaver Argiope argentata (Araneidae), we estimated transcript abundance using two sequencing methods: RNAseq reads from throughout the length of mRNA molecules, and 3′ digital gene expression reads from the 3′ region of mRNA molecules. Both analyses identified similar sets of genes as differentially expressed when comparing tubuliform and nonsilk gland tissue. However, incompletely assembled silk gene transcripts were identified as differentially expressed because of RNAseq read alignments to highly repetitive regions, confounding interpretation of RNAseq results. Homologues of egg case silk protein (ECP) genes were upregulated in tubuliform glands. This discovery is the first description of ECP homologues in an araneid. We also propose additional candidate genes involved in synthesis of tubuliform or other silk types.

Spiders in the superfamily Araneoidea produce viscous glue from aggregate silk glands. Aggregate glue coats prey-capture threads and hampers the escape of prey from webs, thereby increasing the foraging success of spiders. cDNAs for Aggregate Spider Glue 1 (ASG1) and 2 (ASG2) have been previously described from the golden orb-weaver, Nephila clavipes and Western black widow, Latrodectus hesperus. To further investigate aggregate glues, we assembled ASG1 and ASG2 from genomic target capture libraries constructed from three species of cob-web weavers and three species of orb-web weavers, all araneoids. We show that ASG1 is unlikely to be a glue, but rather is part of a widespread arthropod gene family, the peritrophic matrix proteins. For ASG2, we demonstrate its remarkable architectural and sequence similarities to spider silk fibroins, indicating that ASG2 is a member of the spidroin gene family. Thus, spidroins have diversified into glues in addition to task-specific, high performance fibers.

Spider silk research has largely focused on spidroins, proteins that are the primary components of spider silk fibers. Although a number of spidroins have been characterized, other types of proteins associated with silk synthesis are virtually unknown. Previous analyses of tissue-specific RNA-seq libraries identified 647 predicted genes that were differentially expressed in silk glands of the Western black widow, Latrodectus hesperus. Only ∼5% of these silk-gland specific transcripts (SSTs) encode spidroins; although the remaining predicted genes presumably encode other proteins associated with silk production, this is mostly unverified. Here, we used proteomic analysis of multiple silk glands and dragline silk fiber to investigate the translation of the differentially expressed genes. We find 48 proteins encoded by the differentially expressed transcripts in L. hesperus major ampullate, minor ampullate, and tubuliform silk glands and detect 17 SST encoded proteins in major ampullate silk fibers. The observed proteins include known silk-related proteins, but most are uncharacterized, with no annotation. These unannotated proteins likely include novel silk-associated proteins. Major and minor ampullate glands have the highest overlap of identified proteins, consistent with their shared, distinctive ampullate shape and the overlapping functions of major and minor ampullate silks. Our study substantiates and prioritizes predictions from differential expression analysis of spider silk gland transcriptomes.

Background Animal venoms attract enormous interest given their potential for pharmacological discovery and understanding the evolution of natural chemistries. Next-generation transcriptomics and proteomics provide unparalleled, but underexploited, capabilities for venom characterization. We combined multi-tissue RNA-Seq with mass spectrometry and bioinformatic analyses to determine venom gland specific transcripts and venom proteins from the Western black widow spider (Latrodectus hesperus) and investigated their evolution. Results We estimated expression of 97,217 L. hesperus transcripts in venom glands relative to silk and cephalothorax tissues. We identified 695 venom gland specific transcripts (VSTs), many of which BLAST and GO term analyses indicate may function as toxins or their delivery agents. ~38% of VSTs had BLAST hits, including latrotoxins, inhibitor cystine knot toxins, CRISPs, hyaluronidases, chitinase, and proteases, and 59% of VSTs had predicted protein domains. Latrotoxins are venom toxins that cause massive neurotransmitter release from vertebrate or invertebrate neurons. We discovered ≥ 20 divergent latrotoxin paralogs expressed in L. hesperus venom glands, significantly increasing this biomedically important family. Mass spectrometry of L. hesperus venom identified 49 proteins from VSTs, 24 of which BLAST to toxins. Phylogenetic analyses showed venom gland specific gene family expansions and shifts in tissue expression. Conclusions Quantitative expression analyses comparing multiple tissues are necessary to identify venom gland specific transcripts. We present a black widow venom specific exome that uncovers a trove of diverse toxins and associated proteins, suggesting a dynamic evolutionary history. This justifies a reevaluation of the functional activities of black widow venom in light of its emerging complexity. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-366) contains supplementary material, which is available to authorized users.

Summary Investigations into the role of tomato ARF6 and ARF8 reveal that they are critical components in floral and gynoecium development before anthesis., Auxin regulates the expression of diverse genes that affect plant growth and development. This regulation requires AUXIN RESPONSE FACTORS (ARFs) that bind to the promoter regions of these genes. ARF6 and ARF8 in Arabidopsis thaliana are required to promote inflorescence stem elongation and late stages of petal, stamen, and gynoecium development. All seed plants studied thus far have ARF6 and ARF8 orthologues as well as the microRNA miR167, which targets ARF6 and ARF8. Whether these genes have broadly conserved roles in flower development is not known. To address this question, the effects of down-regulation of ARF6 and ARF8 were investigated through transgenic expression of Arabidopsis MIR167a in tomato, which diverged from Arabidopsis before the radiation of dicotyledonous plants approximately 90–112 million years ago. The transgenic tomato plants overexpressing MIR167a exhibited reductions in leaf size and internode length as well as shortened petals, stamens, and styles. More significantly, the transgenic plants were female-sterile as a result of failure of wild-type pollen to germinate on the stigma surface and/or to grow through the style. RNA-Seq analysis identified many genes with significantly altered expression patterns, including those encoding products with functions in ‘transcription regulation’, ‘cell wall’ and ‘lipid metabolism’ categories. Putative orthologues of a subset of these genes were also differentially expressed in Arabidopsis arf6 arf8 mutant flowers. These results thus suggest that ARF6 and ARF8 have conserved roles in controlling growth and development of vegetative and flower organs in dicots.

Ciliated protozoa have two kinds of nuclei: Macronuclei (MAC) and Micronuclei (MIC). In some ciliate classes, such as spirotrichs, most genes undergo several layers of DNA rearrangement during macronuclear development. Because of such processes, these organisms provide ideal systems for studying mechanisms of recombination and gene rearrangement. Here, we describe a database that contains all spirotrich genes for which both MAC and MIC versions are sequenced, with consistent annotation and easy access to all the features. An interface to query the database is available at http://oxytricha.princeton.edu/dimorphism/database.htm.

Background Spiders (Order Araneae) are essential predators in every terrestrial ecosystem largely because they have evolved potent arsenals of silk and venom. Spider silks are high performance materials made almost entirely of proteins, and thus represent an ideal system for investigating genome level evolution of novel protein functions. However, genomic level resources remain limited for spiders. Results We de novo assembled a transcriptome for the Western black widow (Latrodectus hesperus) from deeply sequenced cDNAs of three tissue types. Our multi-tissue assembly contained ~100,000 unique transcripts, of which > 27,000 were annotated by homology. Comparing transcript abundance among the different tissues, we identified 647 silk gland-specific transcripts, including the few known silk fiber components (e.g. six spider fibroins, spidroins). Silk gland specific transcripts are enriched compared to the entire transcriptome in several functions, including protein degradation, inhibition of protein degradation, and oxidation-reduction. Phylogenetic analyses of 37 gene families containing silk gland specific transcripts demonstrated novel gene expansions within silk glands, and multiple co-options of silk specific expression from paralogs expressed in other tissues. Conclusions We propose a transcriptional program for the silk glands that involves regulating gland specific synthesis of silk fiber and glue components followed by protecting and processing these components into functional fibers and glues. Our black widow silk gland gene repertoire provides extensive expansion of resources for biomimetic applications of silk in industry and medicine. Furthermore, our multi-tissue transcriptome facilitates evolutionary analysis of arachnid genomes and adaptive protein systems. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-365) contains supplementary material, which is available to authorized users.