Introducing the Bacterial and Viral Bioinformatics Resource Center (BV-BRC): a resource combining PATRIC, IRD and ViPR.

The National Institute of Allergy and Infectious Diseases (NIAID) established the Bioinformatics Resource Center (BRC) program to assist researchers with analyzing the growing body of genome sequence and other omics-related data. In this report, we describe the merger of the PAThosystems Resource Integration Center (PATRIC), the Influenza Research Database (IRD) and the Virus Pathogen Database and Analysis Resource (ViPR) BRCs to form the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) https://www.bv-brc.org/. The combined BV-BRC leverages the functionality of the bacterial and viral resources to provide a unified data model, enhanced web-based visualization and analysis tools, bioinformatics services, and a powerful suite of command line tools that benefit the bacterial and viral research communities.

US National Institutes of Health Prioritization of SARS-CoV-2 Variants.

Since late 2020, SARS-CoV-2 variants have regularly emerged with competitive and phenotypic differences from previously circulating strains, sometimes with the potential to escape from immunity produced by prior exposure and infection. The Early Detection group is one of the constituent groups of the US National Institutes of Health National Institute of Allergy and Infectious Diseases SARS-CoV-2 Assessment of Viral Evolution program. The group uses bioinformatic methods to monitor the emergence, spread, and potential phenotypic properties of emerging and circulating strains to identify the most relevant variants for experimental groups within the program to phenotypically characterize. Since April 2021, the group has prioritized variants monthly. Prioritization successes include rapidly identifying most major variants of SARS-CoV-2 and providing experimental groups within the National Institutes of Health program easy access to regularly updated information on the recent evolution and epidemiology of SARS-CoV-2 that can be used to guide phenotypic investigations.

Population-level mathematical modeling of antimicrobial resistance: a systematic review.

BACKGROUND: Mathematical transmission models are increasingly used to guide public health interventions for infectious diseases, particularly in the context of emerging pathogens; however, the contribution of modeling to the growing issue of antimicrobial resistance (AMR) remains unclear. Here, we systematically evaluate publications on population-level transmission models of AMR over a recent period (2006-2016) to gauge the state of research and identify gaps warranting further work. METHODS: We performed a systematic literature search of relevant databases to identify transmission studies of AMR in viral, bacterial, and parasitic disease systems. We analyzed the temporal, geographic, and subject matter trends, described the predominant medical and behavioral interventions studied, and identified central findings relating to key pathogens. RESULTS: We identified 273 modeling studies; the majority of which (> 70%) focused on 5 infectious diseases (human immunodeficiency virus (HIV), influenza virus, Plasmodium falciparum (malaria), Mycobacterium tuberculosis (TB), and methicillin-resistant Staphylococcus aureus (MRSA)). AMR studies of influenza and nosocomial pathogens were mainly set in industrialized nations, while HIV, TB, and malaria studies were heavily skewed towards developing countries. The majority of articles focused on AMR exclusively in humans (89%), either in community (58%) or healthcare (27%) settings. Model systems were largely compartmental (76%) and deterministic (66%). Only 43% of models were calibrated against epidemiological data, and few were validated against out-of-sample datasets (14%). The interventions considered were primarily the impact of different drug regimens, hygiene and infection control measures, screening, and diagnostics, while few studies addressed de novo resistance, vaccination strategies, economic, or behavioral changes to reduce antibiotic use in humans and animals. CONCLUSIONS: The AMR modeling literature concentrates on disease systems where resistance has been long-established, while few studies pro-actively address recent rise in resistance in new pathogens or explore upstream strategies to reduce overall antibiotic consumption. Notable gaps include research on emerging resistance in Enterobacteriaceae and Neisseria gonorrhoeae; AMR transmission at the animal-human interface, particularly in agricultural and veterinary settings; transmission between hospitals and the community; the role of environmental factors in AMR transmission; and the potential of vaccines to combat AMR.

The extraordinary evolutionary history of the reticuloendotheliosis viruses.

The reticuloendotheliosis viruses (REVs) comprise several closely related amphotropic retroviruses isolated from birds. These viruses exhibit several highly unusual characteristics that have not so far been adequately explained, including their extremely close relationship to mammalian retroviruses, and their presence as endogenous sequences within the genomes of certain large DNA viruses. We present evidence for an iatrogenic origin of REVs that accounts for these phenomena. Firstly, we identify endogenous retroviral fossils in mammalian genomes that share a unique recombinant structure with REVs-unequivocally demonstrating that REVs derive directly from mammalian retroviruses. Secondly, through sequencing of archived REV isolates, we confirm that contaminated Plasmodium lophurae stocks have been the source of multiple REV outbreaks in experimentally infected birds. Finally, we show that both phylogenetic and historical evidence support a scenario wherein REVs originated as mammalian retroviruses that were accidentally introduced into avian hosts in the late 1930s, during experimental studies of P. lophurae, and subsequently integrated into the fowlpox virus (FWPV) and gallid herpesvirus type 2 (GHV-2) genomes, generating recombinant DNA viruses that now circulate in wild birds and poultry. Our findings provide a novel perspective on the origin and evolution of REV, and indicate that horizontal gene transfer between virus families can expand the impact of iatrogenic transmission events.

Emergence and pandemic spread of small ruminant lentiviruses.

Small ruminant lentiviruses (SRLVs) cause chronic, persistent infections in populations of domestic sheep (Ovis aries) and goats (Capra hircus) worldwide. The vast majority of SRLV infections involve two genotypes (A and B) that spread in association with the emergence of global livestock trade. However, SRLVs have likely been present in Eurasian ruminant populations since at least the early Neolithic period. Here, we use phylogenetic and phylogeographic approaches to reconstruct the origin of pandemic SRLV strains and infer their historical pattern of global spread. We constructed an open computational resource ('Lentivirus-GLUE') via which an up-to-date database of published SRLV sequences, multiple sequence alignments (MSAs), and sequence-associated metadata can be maintained. We used data collated in Lentivirus-GLUE to perform a comprehensive phylogenetic investigation of global SRLV diversity. Phylogenies reconstructed from genome-length alignments reveal that the deep divisions in the SRLV phylogeny are consistent with an ancient split into Eastern (A-like) and Western (B-like) lineages as agricultural systems disseminated out of domestication centres during the Neolithic period. These findings are also consistent with historical and phylogeographic evidence linking the early 20th century emergence of SRLV-A to the international export of Central Asian Karakul sheep. Investigating the global diversity of SRLVs can help reveal how anthropogenic factors have impacted the ecology and evolution of livestock diseases. The open resources generated in our study can expedite these studies and can also serve more broadly to facilitate the use of genomic data in SRLV diagnostics and research.

Regulation of Hsp90 client proteins by a Cullin5-RING E3 ubiquitin ligase.

We report a link between Cullin5 (Cul5) E3 ubiquitin ligase and the heat shock protein 90 (Hsp90) chaperone complex. Hsp90 participates in the folding of its client proteins into their functional conformation. Many Hsp90 clients have been reported to be aberrantly expressed in a number of cancers. We demonstrate Cul5 interaction with members of the Hsp90 chaperone complex as well as the Hsp90 client, ErbB2. We observed recruitment of Cul5 to the site of ErbB2 at the plasma membrane and subsequent induction of polyubiquitination and proteasomal degradation. We also demonstrate Cul5 involvement in regulation of another Hsp90 client, Hif-1alpha. We observed Cul5 degradation of ErbB2 to occur independently of ElonginB-ElonginC function. The involvement of Cul5 in Hsp90 client regulation has implications in the effectiveness of Hsp90 targeted chemotherapy, which is currently undergoing clinical trials. The link between Cul5 and Hsp90 client regulation may represent an avenue for cancer drug development.

Genomic evolution of the Coronaviridae family.

The current outbreak of coronavirus disease-2019 (COVID-19) caused by SARS-CoV-2 poses unparalleled challenges to global public health. SARS-CoV-2 is a Betacoronavirus, one of four genera belonging to the Coronaviridae subfamily Orthocoronavirinae. Coronaviridae, in turn, are members of the order Nidovirales, a group of enveloped, positive-stranded RNA viruses. Here we present a systematic phylogenetic and evolutionary study based on protein domain architecture, encompassing the entire proteomes of all Orthocoronavirinae, as well as other Nidovirales. This analysis has revealed that the genomic evolution of Nidovirales is associated with extensive gains and losses of protein domains. In Orthocoronavirinae, the sections of the genomes that show the largest divergence in protein domains are found in the proteins encoded in the amino-terminal end of the polyprotein (PP1ab), the spike protein (S), and many of the accessory proteins. The diversity among the accessory proteins is particularly striking, as each subgenus possesses a set of accessory proteins that is almost entirely specific to that subgenus. The only notable exception to this is ORF3b, which is present and orthologous over all Alphacoronaviruses. In contrast, the membrane protein (M), envelope small membrane protein (E), nucleoprotein (N), as well as proteins encoded in the central and carboxy-terminal end of PP1ab (such as the 3C-like protease, RNA-dependent RNA polymerase, and Helicase) show stable domain architectures across all Orthocoronavirinae. This comprehensive analysis of the Coronaviridae domain architecture has important implication for efforts to develop broadly cross-protective coronavirus vaccines.

Early detection of emerging SARS-CoV-2 variants of interest for experimental evaluation.

Since the beginning of the COVID-19 pandemic, SARS-CoV-2 has demonstrated its ability to rapidly and continuously evolve, leading to the emergence of thousands of different sequence variants, many with distinctive phenotypic properties. Fortunately, the broad application of next generation sequencing (NGS) across the globe has produced a wealth of SARS-CoV-2 genome sequences, offering a comprehensive picture of how this virus is evolving so that accurate diagnostics, reliable therapeutics, and prophylactic vaccines against COVID-19 can be developed and maintained. The millions of SARS-CoV-2 sequences deposited into genomic sequencing databases, including GenBank, BV-BRC, and GISAID, are annotated with the dates and geographic locations of sample collection, and can be aligned to and compared with the Wuhan-Hu-1 reference genome to extract their constellation of nucleotide and amino acid substitutions. By aggregating these data into concise datasets, the spread of variants through space and time can be assessed. Variant tracking efforts have initially focused on the Spike protein due to its critical role in viral tropism and antibody neutralization. To identify emerging variants of concern as early as possible, we developed a computational pipeline to process the genomic data and assign risk scores based on both epidemiological and functional parameters. Epidemiological dynamics are used to identify variants exhibiting substantial growth over time and spread across geographical regions. Experimental data that quantify Spike protein regions targeted by adaptive immunity and critical for other virus characteristics are used to predict variants with consequential immunogenic and pathogenic impacts. The growth assessment and functional impact scores are combined to produce a Composite Score for any set of Spike substitutions detected. With this systematic method to routinely score and rank emerging variants, we have established an approach to identify threatening variants early and prioritize them for experimental evaluation.

Host restriction of HIV-1 by APOBEC3 and viral evasion through Vif.

The arms race between virus and host is a constant battle. APOBEC3 proteins are known to be potent innate cellular defenses against both endogenous retroelements and diverse retroviruses. However, retroviruses have developed their own methods to launch counter-strikes. Most primate lentiviruses encode a protein called the viral infectivity factor (Vif). Vif induces targeted destruction of APOBEC3 proteins by hijacking the cellular ubiquitin-proteasome pathway. Here we review the research that led up to the identification of A3G, the mechanisms by which APOBEC3 proteins can inhibit retroelements, and the counter-mechanisms that HIV-1 Vif has developed to evade its antiviral activities.

Characterization of a synthetic human LINE-1 retrotransposon ORFeus-Hs.

Long interspersed elements, type 1(LINE-1, L1) are the most abundant and only active autonomous retrotransposons in the human genome. Native L1 elements are inefficiently expressed because of a transcription elongation defect thought to be caused by high adenosine content in L1 sequences. Previously, we constructed a highly active synthetic mouse L1 element (ORFeus-Mm), partially by reducing the nucleotide composition bias. As a result, the transcript abundance of ORFeus-Mm was greatly increased, and its retrotransposition frequency was > 200-fold higher than its native counterpart. In this paper, we report a synthetic human L1 element (ORFeus-Hs) synthesized using a similar strategy. The adenosine content of the L1 open reading frames (ORFs) was reduced from 40% to 27% by changing 25% of the bases in the ORFs, without altering the amino acid sequence. By studying a series of native/synthetic chimeric elements, we observed increased levels of full-length L1 RNA and ORF1 protein and retrotransposition frequency, mostly proportional to increased fraction of synthetic sequence. Overall, the fully synthetic ORFeus-Hs has > 40-fold more RNA but is at most only ~threefold more active than its native counterpart (L1RP); however, its absolute retrotransposition activity is similar to ORFeus-Mm. Owing to the elevated expression of the L1 RNA/protein and its high retrotransposition ability, ORFeus-Hs and its chimeric derivatives will be useful tools for mechanistic L1 studies and mammalian genome manipulation.

Distinct determinants in HIV-1 Vif and human APOBEC3 proteins are required for the suppression of diverse host anti-viral proteins.

BACKGROUND: APOBEC3G (A3G) and related cytidine deaminases of the APOBEC3 family of proteins are potent inhibitors of many retroviruses, including HIV-1. Formation of infectious HIV-1 requires the suppression of multiple cytidine deaminases by Vif. HIV-1 Vif suppresses various APOBEC3 proteins through the common mechanism of recruiting the Cullin5-ElonginB-ElonginC E3 ubiquitin ligase to induce target protein polyubiquitination and proteasome-mediated degradation. The domains in Vif and various APOBEC3 proteins required for APOBEC3 recognition and degradation have not been fully characterized. METHODS AND FINDINGS: In the present study, we have demonstrated that the regions of APOBEC3F (A3F) that are required for its HIV-1-mediated binding and degradation are distinct from those reported for A3G. We found that the C-terminal cytidine deaminase domain (C-CDD) of A3F alone is sufficient for its interaction with HIV-1 Vif and its Vif-mediated degradation. We also observed that the domains of HIV-1 Vif that are uniquely required for its functional interaction with full-length A3F are also required for the degradation of the C-CDD of A3F; in contrast, those Vif domains that are uniquely required for functional interaction with A3G are not required for the degradation of the C-CDD of A3F. Interestingly, the HIV-1 Vif domains required for the degradation of A3F are also required for the degradation of A3C and A3DE. On the other hand, the Vif domains uniquely required for the degradation of A3G are dispensable for the degradation of cytidine deaminases A3C and A3DE. CONCLUSIONS: Our data suggest that distinct regions of A3F and A3G are targeted by HIV-1 Vif molecules. However, HIV-1 Vif suppresses A3F, A3C, and A3DE through similar recognition determinants, which are conserved among Vif molecules from diverse HIV-1 strains. Mapping these determinants may be useful for the design of novel anti-HIV inhibitors.

Differential inhibition of long interspersed element 1 by APOBEC3 does not correlate with high-molecular-mass-complex formation or P-body association.

The human cytidine deaminase APOBEC3G (A3G) and other APOBEC3 proteins exhibit differential inhibitory activities against diverse endogenous retroelements and retroviruses, including Vif-deficient human immunodeficiency virus type 1. The potential inhibitory activity of human APOBEC proteins against long interspersed element 1 (LINE-1) has not been fully evaluated. Here, we demonstrate inhibition of LINE-1 by multiple human APOBEC3 cytidine deaminases, including previously unreported activity for A3DE and A3G. More ancient members of APOBEC, cytidine deaminases AID and APOBEC2, had no detectable activity against LINE-1. A3A, which did not form high-molecular-mass (HMM) complexes and interacted poorly with P bodies, was the most potent inhibitor of LINE-1. A3A specifically recognizes LINE-1 RNA but not the other cellular RNAs tested. However, in the presence of LINE-1, A3A became associated with HMM complexes containing LINE-1 RNA. The ability of A3A to recognize LINE-1 RNA required its catalytic domain and was important for its LINE-1 suppression. Although the mechanism of LINE-1 restriction did not seem to involve DNA editing, A3A inhibited the accumulation of nascent LINE-1 DNA, suggesting interference with LINE-1 reverse transcription and/or integration or intracellular movement of LINE-1 ribonucleoprotein. Thus, association with P bodies or cellular HMM complexes could not predict the potency of APOBEC3 anti-LINE-1 activities. The catalytic domain of APOBEC3 proteins may be important for proper folding and target factors such as RNA or protein interaction in addition to cytidine deamination.