Rapid Emergence of SARS-CoV-2 in the Greater New York Metropolitan Area: Geolocation, Demographics, Positivity Rates, and Hospitalization for 46 793 Persons Tested by Northwell Health.

BACKGROUND: In March 2020, the greater New York metropolitan area became an epicenter for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. The initial evolution of case incidence has not been well characterized. METHODS: Northwell Health Laboratories tested 46 793 persons for SARS-CoV-2 from 4 March through 10 April. The primary outcome measure was a positive reverse transcription-polymerase chain reaction test for SARS-CoV-2. The secondary outcomes included patient age, sex, and race, if stated; dates the specimen was obtained and the test result; clinical practice site sources; geolocation of patient residence; and hospitalization. RESULTS: From 8 March through 10 April, a total of 26 735 of 46 793 persons (57.1%) tested positive for SARS-CoV-2. Males of each race were disproportionally more affected than females above age 25, with a progressive male predominance as age increased. Of the positive persons, 7292 were hospitalized directly upon presentation; an additional 882 persons tested positive in an ambulatory setting before subsequent hospitalization, a median of 4.8 days later. Total hospitalization rate was thus 8174 persons (30.6% of positive persons). There was a broad range (>10-fold) in the cumulative number of positive cases across individual zip codes following documented first caseincidence. Test positivity was greater for persons living in zip codes with lower annual household income. CONCLUSIONS: Our data reveal that SARS-CoV-2 incidence emerged rapidly and almost simultaneously across a broad demographic population in the region. These findings support the premise that SARS-CoV-2 infection was widely distributed prior to virus testing availability.

BpWrapper: BioPerl-based sequence and tree utilities for rapid prototyping of bioinformatics pipelines.

BACKGROUND: Automated bioinformatics workflows are more robust, easier to maintain, and results more reproducible when built with command-line utilities than with custom-coded scripts. Command-line utilities further benefit by relieving bioinformatics developers to learn the use of, or to interact directly with, biological software libraries. There is however a lack of command-line utilities that leverage popular Open Source biological software toolkits such as BioPerl ( http://bioperl.org ) to make many of the well-designed, robust, and routinely used biological classes available for a wider base of end users. RESULTS: Designed as standard utilities for UNIX-family operating systems, BpWrapper makes functionality of some of the most popular BioPerl modules readily accessible on the command line to novice as well as to experienced bioinformatics practitioners. The initial release of BpWrapper includes four utilities with concise command-line user interfaces, bioseq, bioaln, biotree, and biopop, specialized for manipulation of molecular sequences, sequence alignments, phylogenetic trees, and DNA polymorphisms, respectively. Over a hundred methods are currently available as command-line options and new methods are easily incorporated. Performance of BpWrapper utilities lags that of precompiled utilities while equivalent to that of other utilities based on BioPerl. BpWrapper has been tested on BioPerl Release 1.6, Perl versions 5.10.1 to 5.25.10, and operating systems including Apple macOS, Microsoft Windows, and GNU/Linux. Release code is available from the Comprehensive Perl Archive Network (CPAN) at https://metacpan.org/pod/Bio::BPWrapper . Source code is available on GitHub at https://github.com/bioperl/p5-bpwrapper . CONCLUSIONS: BpWrapper improves on existing sequence utilities by following the design principles of Unix text utilities such including a concise user interface, extensive command-line options, and standard input/output for serialized operations. Further, dozens of novel methods for manipulation of sequences, alignments, and phylogenetic trees, unavailable in existing utilities (e.g., EMBOSS, Newick Utilities, and FAST), are provided. Bioinformaticians should find BpWrapper useful for rapid prototyping of workflows on the command-line without creating custom scripts for comparative genomics and other bioinformatics applications.

Evaluation of a human adenovirus viral load assay using the Altona RealStar® PCR test.

This study evaluated the performance of the Altona Diagnostics RealStar® Adenovirus Research Use Only (RUO) real-time PCR reagents for HAdV quantitation in plasma samples from immunodeficient patients. The assay was linear from 2.30-9.17 log10 copies/mL (coefficient of determination; R2=0.998) with limits of detection and quantification of 2.19 log10 and 2.30 log10 copies/mL (>95% positivity rate), respectively. Assay precision was highly reproducible with coefficients of variance ranging from 0% to 4.7%. A comparison of 66 matched samples showed good agreement (R2=0.845) between the Altona and the reference laboratory assay, with an average negative bias (-0.24 log10 copies/mL). Genotyping analysis demonstrated that HAdV species B and C accounted for 77% of the positive samples. A significant (≥0.9 log10) difference in quantitation between both tests was found for three HAdV types (HAdV types A12, B14 and F41). In conclusion, the Altona RealStar® test is a reliable and sensitive assay for HAdV DNA quantitation.

Gastric Carcinomas With Lymphoid Stroma: An Evaluation of the Histopathologic and Molecular Features.

Gastric carcinoma with lymphoid stroma is an uncommon variant enriched for mutually exclusive Epstein-Barr virus (EBV) positivity and mismatch repair (MMR) deficiency. We performed this study to evaluate molecular alterations in this morphologically homogeneous subtype and compare them with 295 conventional gastric cancers analyzed in The Cancer Genome Atlas study. We identified 31 study cases and subjected them to in situ hybridization for EBV-encoded RNAs and assessment for MMR status. Immunostains for PD-L1, ß-catenin, and HER2 were performed; extracted DNA was sequenced with a Comprehensive Cancer Panel. Most study patients were older adult men with stage I or II disease (76%). Tumors were classified as EBV/MMR-proficient (MMR-P) (n=7), EBV/MMR deficient (n=12), and EBV/MMR-P (n=12). EBV/MMR-P tumors were usually located in the proximal stomach (83%) and showed heterogenous growth patterns with glandular differentiation (83%). Tumors in all groups showed numerous tumor infiltrating lymphocytes and PD-L1 expression, infrequent nuclear ß-catenin accumulation (10%), and lacked both membranous HER2 staining and HER2 amplification. EBV/MMR-deficient tumors showed significantly higher tumor mutation burden (P=0.001) and KRAS alterations (56%) compared with EBV/MMR-P tumors (9%, P=0.05). TP53 variants were more common among EBV/MMR-P tumors (82%) compared with EBV/MMR proficient (0%, P=0.01) and EBV/MMR-deficient (11%, P<0.01) tumors. Alterations in KRAS, ARID1A, PIK3CA, and TP53 followed similar patterns of distribution compared with The Cancer Genome Atlas dataset. We conclude that gastric carcinomas with lymphoid stroma show a spectrum of molecular changes and frequent PD-L1 expression, raising the possibility that this subgroup of tumors may be susceptible to checkpoint inhibitors and/or agents that target receptor tyrosine kinase-mediated signaling.

Novel Organelles with Elements of Bacterial and Eukaryotic Secretion Systems Weaponize Parasites of Drosophila.

The evolutionary success of parasitoid wasps, a highly diverse group of insects widely used in biocontrol, depends on a variety of life history strategies in conflict with those of their hosts [1]. Drosophila melanogaster is a natural host of parasitic wasps of the genus Leptopilina. Attack by L. boulardi (Lb), a specialist wasp to flies of the melanogaster group, activates NF-κB-mediated humoral and cellular immunity. Inflammatory blood cells mobilize and encapsulate Lb eggs and embryos [2-5]. L. heterotoma (Lh), a generalist wasp, kills larval blood cells and actively suppresses immune responses. Spiked virus-like particles (VLPs) in wasp venom have clearly been linked to wasps' successful parasitism of Drosophila [6], but the composition of VLPs and their biotic nature have remained mysterious. Our proteomics studies reveal that VLPs lack viral coat proteins but possess a pharmacopoeia of (1) the eukaryotic vesicular transport system, (2) immunity, and (3) previously unknown proteins. These novel proteins distinguish Lh from Lb VLPs; notably, some proteins specific to Lh VLPs possess sequence similarities with bacterial secretion system proteins. Structure-informed analyses of an abundant Lh VLP surface and spike-tip protein, p40, reveal similarities to the needle-tip invasin proteins SipD and IpaD of Gram-negative bacterial type-3 secretion systems that breach immune barriers and deliver virulence factors into mammalian cells. Our studies suggest that Lh VLPs represent a new class of extracellular organelles and share pathways for protein delivery with both eukaryotic microvesicles and bacterial surface secretion systems. Given their mixed prokaryotic and eukaryotic properties, we propose the term mixed-strategy extracellular vesicle (MSEV) to replace VLP.

RNA-DNA sequence differences in Saccharomyces cerevisiae.

Alterations of RNA sequences and structures, such as those from editing and alternative splicing, result in two or more RNA transcripts from a DNA template. It was thought that in yeast, RNA editing only occurs in tRNAs. Here, we found that Saccharomyces cerevisiae have all 12 types of RNA-DNA sequence differences (RDDs) in the mRNA. We showed these sequence differences are propagated to proteins, as we identified peptides encoded by the RNA sequences in addition to those by the DNA sequences at RDD sites. RDDs are significantly enriched at regions with R-loops. A screen of yeast mutants showed that RDD formation is affected by mutations in genes regulating R-loops. Loss-of-function mutations in ribonuclease H, senataxin, and topoisomerase I that resolve RNA-DNA hybrids lead to increases in RDD frequency. Our results demonstrate that RDD is a conserved process that diversifies transcriptomes and proteomes and provide a mechanistic link between R-loops and RDDs.

From Mollusks to Medicine: A Venomics Approach for the Discovery and Characterization of Therapeutics from Terebridae Peptide Toxins.

Animal venoms comprise a diversity of peptide toxins that manipulate molecular targets such as ion channels and receptors, making venom peptides attractive candidates for the development of therapeutics to benefit human health. However, identifying bioactive venom peptides remains a significant challenge. In this review we describe our particular venomics strategy for the discovery, characterization, and optimization of Terebridae venom peptides, teretoxins. Our strategy reflects the scientific path from mollusks to medicine in an integrative sequential approach with the following steps: (1) delimitation of venomous Terebridae lineages through taxonomic and phylogenetic analyses; (2) identification and classification of putative teretoxins through omics methodologies, including genomics, transcriptomics, and proteomics; (3) chemical and recombinant synthesis of promising peptide toxins; (4) structural characterization through experimental and computational methods; (5) determination of teretoxin bioactivity and molecular function through biological assays and computational modeling; (6) optimization of peptide toxin affinity and selectivity to molecular target; and (7) development of strategies for effective delivery of venom peptide therapeutics. While our research focuses on terebrids, the venomics approach outlined here can be applied to the discovery and characterization of peptide toxins from any venomous taxa.

Genetic variation in insulin-induced kinase signaling.

Individual differences in sensitivity to insulin contribute to disease susceptibility including diabetes and metabolic syndrome. Cellular responses to insulin are well studied. However, which steps in these response pathways differ across individuals remains largely unknown. Such knowledge is needed to guide more precise therapeutic interventions. Here, we studied insulin response and found extensive individual variation in the activation of key signaling factors, including ERK whose induction differs by more than 20-fold among our subjects. This variation in kinase activity is propagated to differences in downstream gene expression response to insulin. By genetic analysis, we identified cis-acting DNA variants that influence signaling response, which in turn affects downstream changes in gene expression and cellular phenotypes, such as protein translation and cell proliferation. These findings show that polymorphic differences in signal transduction contribute to individual variation in insulin response, and suggest kinase modulators as promising therapeutics for diseases characterized by insulin resistance.

Molecular Diversity and Gene Evolution of the Venom Arsenal of Terebridae Predatory Marine Snails.

Venom peptides from predatory organisms are a resource for investigating evolutionary processes such as adaptive radiation or diversification, and exemplify promising targets for biomedical drug development. Terebridae are an understudied lineage of conoidean snails, which also includes cone snails and turrids. Characterization of cone snail venom peptides, conotoxins, has revealed a cocktail of bioactive compounds used to investigate physiological cellular function, predator-prey interactions, and to develop novel therapeutics. However, venom diversity of other conoidean snails remains poorly understood. The present research applies a venomics approach to characterize novel terebrid venom peptides, teretoxins, from the venom gland transcriptomes of Triplostephanus anilis and Terebra subulata. Next-generation sequencing and de novo assembly identified 139 putative teretoxins that were analyzed for the presence of canonical peptide features as identified in conotoxins. To meet the challenges of de novo assembly, multiple approaches for cross validation of findings were performed to achieve reliable assemblies of venom duct transcriptomes and to obtain a robust portrait of Terebridae venom. Phylogenetic methodology was used to identify 14 teretoxin gene superfamilies for the first time, 13 of which are unique to the Terebridae. Additionally, basic local algorithm search tool homology-based searches to venom-related genes and posttranslational modification enzymes identified a convergence of certain venom proteins, such as actinoporin, commonly found in venoms. This research provides novel insights into venom evolution and recruitment in Conoidean predatory marine snails and identifies a plethora of terebrid venom peptides that can be used to investigate fundamental questions pertaining to gene evolution.

Phylogenomic identification of regulatory sequences in bacteria: an analysis of statistical power and an application to Borrelia burgdorferi sensu lato.

UNLABELLED: Phylogenomic footprinting is an approach for ab initio identification of genome-wide regulatory elements in bacterial species based on sequence conservation. The statistical power of the phylogenomic approach depends on the degree of sequence conservation, the length of regulatory elements, and the level of phylogenetic divergence among genomes. Building on an earlier model, we propose a binomial model that uses synonymous tree lengths as neutral expectations for determining the statistical significance of conserved intergenic spacer (IGS) sequences. Simulations show that the binomial model is robust to variations in the value of evolutionary parameters, including base frequencies and the transition-to-transversion ratio. We used the model to search for regulatory sequences in the Lyme disease species group (Borrelia burgdorferi sensu lato) using 23 genomes. The model indicates that the currently available set of Borrelia genomes would not yield regulatory sequences shorter than five bases, suggesting that genome sequences of additional B. burgdorferi sensu lato species are needed. Nevertheless, we show that previously known regulatory elements are indeed strongly conserved in sequence or structure across these Borrelia species. Further, we predict with sufficient confidence two new RpoS binding sites, 39 promoters, 19 transcription terminators, 28 noncoding RNAs, and four sets of coregulated genes. These putative cis- and trans-regulatory elements suggest novel, Borrelia-specific mechanisms regulating the transition between the tick and host environments, a key adaptation and virulence mechanism of B. burgdorferi. Alignments of IGS sequences are available on BorreliaBase.org, an online database of orthologous open reading frame (ORF) and IGS sequences in Borrelia. IMPORTANCE: While bacterial genomes contain mostly protein-coding genes, they also house DNA sequences regulating the expression of these genes. Gene regulatory sequences tend to be conserved during evolution. By sequencing and comparing related genomes, one can therefore identify regulatory sequences in bacteria based on sequence conservation. Here, we describe a statistical framework by which one may determine how many genomes need to be sequenced and at what level of evolutionary relatedness in order to achieve a high level of statistical significance. We applied the framework to Borrelia burgdorferi, the Lyme disease agent, and identified a large number of candidate regulatory sequences, many of which are known to be involved in regulating the phase transition between the tick vector and mammalian hosts.

BorreliaBase: a phylogeny-centered browser of Borrelia genomes.

BACKGROUND: The bacterial genus Borrelia (phylum Spirochaetes) consists of two groups of pathogens represented respectively by B. burgdorferi, the agent of Lyme borreliosis, and B. hermsii, the agent of tick-borne relapsing fever. The number of publicly available Borrelia genomic sequences is growing rapidly with the discovery and sequencing of Borrelia strains worldwide. There is however a lack of dedicated online databases to facilitate comparative analyses of Borrelia genomes. DESCRIPTION: We have developed BorreliaBase, an online database for comparative browsing of Borrelia genomes. The database is currently populated with sequences from 35 genomes of eight Lyme-borreliosis (LB) group Borrelia species and 7 Relapsing-fever (RF) group Borrelia species. Distinct from genome repositories and aggregator databases, BorreliaBase serves manually curated comparative-genomic data including genome-based phylogeny, genome synteny, and sequence alignments of orthologous genes and intergenic spacers. CONCLUSIONS: With a genome phylogeny at its center, BorreliaBase allows online identification of hypervariable lipoprotein genes, potential regulatory elements, and recombination footprints by providing evolution-based expectations of sequence variability at each genomic locus. The phylo-centric design of BorreliaBase (http://borreliabase.org) is a novel model for interactive browsing and comparative analysis of bacterial genomes online.