SARS-CoV-2 clade dynamics and their associations with hospitalisations during the first two years of the COVID-19 pandemic.

BACKGROUND: The COVID-19 pandemic was characterised by rapid waves of disease, carried by the emergence of new and more infectious SARS-CoV-2 virus variants. How the pandemic unfolded in various locations during its first two years has yet to be sufficiently covered. To this end, here we are looking at the circulating SARS-CoV-2 variants, their diversity, and hospitalisation rates in Estonia in the period from March 2000 to March 2022. METHODS: We sequenced a total of 27,550 SARS-CoV-2 samples in Estonia between March 2020 and March 2022. High-quality sequences were genotyped and assigned to Nextstrain clades and Pango lineages. We used regression analysis to determine the dynamics of lineage diversity and the probability of clade-specific hospitalisation stratified by age and sex. RESULTS: We successfully sequenced a total of 25,375 SARS-CoV-2 genomes (or 92%), identifying 19 Nextstrain clades and 199 Pango lineages. In 2020 the most prevalent clades were 20B and 20A. The various subsequent waves of infection were driven by 20I (Alpha), 21J (Delta) and Omicron clades 21K and 21L. Lineage diversity via the Shannon index was at its highest during the Delta wave. About 3% of sequenced SARS-CoV-2 samples came from hospitalised individuals. Hospitalisation increased markedly with age in the over-forties, and was negligible in the under-forties. Vaccination decreased the odds of hospitalisation in over-forties. The effect of vaccination on hospitalisation rates was strongly dependent upon age but was clade-independent. People who were infected with Omicron clades had a lower hospitalisation likelihood in age groups of forty and over than was the case with pre-Omicron clades regardless of vaccination status. CONCLUSIONS: COVID-19 disease waves in Estonia were driven by the Alpha, Delta, and Omicron clades. Omicron clades were associated with a substantially lower hospitalisation probability than pre-Omicron clades. The protective effect of vaccination in reducing hospitalisation likelihood was independent of the involved clade.

eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs.

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

The interplay of DNA repair context with target sequence predictably biasses Cas9-generated mutations.

The genome engineering capability of the CRISPR/Cas system depends on the DNA repair machinery to generate the final outcome. Several genes can have an impact on mutations created, but their exact function and contribution to the result of the repair are not completely characterised. This lack of knowledge has limited the ability to comprehend and regulate the editing outcomes. Here, we measure how the absence of 21 repair genes changes the mutation outcomes of Cas9-generated cuts at 2,812 synthetic target sequences in mouse embryonic stem cells. Absence of key non-homologous end joining genes Lig4, Xrcc4, and Xlf abolished small insertions and deletions, while disabling key microhomology-mediated repair genes Nbn and Polq reduced frequency of longer deletions. Complex alleles of combined insertion and deletions were preferentially generated in the absence of Xrcc6. We further discover finer structure in the outcome frequency changes for single nucleotide insertions and deletions between large microhomologies that are differentially modulated by the knockouts. We use the knowledge of the reproducible variation across repair milieus to build predictive models of Cas9 editing results that outperform the current standards. This work improves our understanding of DNA repair gene function, and provides avenues for more precise modulation of CRISPR/Cas9-generated mutations.

Rabies in a postpandemic world: resilient reservoirs, redoubtable riposte, recurrent roadblocks, and resolute recidivism.

Rabies is an ancient disease. Two centuries since Pasteur, fundamental progress occurred in virology, vaccinology, and diagnostics-and an understanding of pathobiology and epizootiology of rabies in testament to One Health-before common terminological coinage. Prevention, control, selective elimination, and even the unthinkable-occasional treatment-of this zoonosis dawned by the twenty-first century. However, in contrast to smallpox and rinderpest, eradication is a wishful misnomer applied to rabies, particularly post-COVID-19 pandemic. Reasons are minion. Polyhostality encompasses bats and mesocarnivores, but other mammals represent a diverse spectrum of potential hosts. While rabies virus is the classical member of the genus, other species of lyssaviruses also cause the disease. Some reservoirs remain cryptic. Although global, this viral encephalitis is untreatable and often ignored. As with other neglected diseases, laboratory-based surveillance falls short of the notifiable ideal, especially in lower- and middle-income countries. Calculation of actual burden defaults to a flux within broad health economic models. Competing priorities, lack of defined, long-term international donors, and shrinking local champions challenge human prophylaxis and mass dog vaccination toward targets of 2030 for even canine rabies impacts. For prevention, all licensed vaccines are delivered to the individual, whether parenteral or oral-essentially 'one and done'. Exploiting mammalian social behaviors, future 'spreadable vaccines' might increase the proportion of immunized hosts per unit effort. However, the release of replication-competent, genetically modified organisms selectively engineered to spread intentionally throughout a population raises significant biological, ethical, and regulatory issues in need of broader, transdisciplinary discourse. How this rather curious idea will evolve toward actual unconventional prevention, control, or elimination in the near term remains debatable. In the interim, more precise terminology and realistic expectations serve as the norm for diverse, collective constituents to maintain progress in the field.

g:Profiler-interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update).

g:Profiler is a reliable and up-to-date functional enrichment analysis tool that supports various evidence types, identifier types and organisms. The toolset integrates many databases, including Gene Ontology, KEGG and TRANSFAC, to provide a comprehensive and in-depth analysis of gene lists. It also provides interactive and intuitive user interfaces and supports ordered queries and custom statistical backgrounds, among other settings. g:Profiler provides multiple programmatic interfaces to access its functionality. These can be easily integrated into custom workflows and external tools, making them valuable resources for researchers who want to develop their own solutions. g:Profiler has been available since 2007 and is used to analyse millions of queries. Research reproducibility and transparency are achieved by maintaining working versions of all past database releases since 2015. g:Profiler supports 849 species, including vertebrates, plants, fungi, insects and parasites, and can analyse any organism through user-uploaded custom annotation files. In this update article, we introduce a novel filtering method highlighting Gene Ontology driver terms, accompanied by new graph visualizations providing a broader context for significant Gene Ontology terms. As a leading enrichment analysis and gene list interoperability service, g:Profiler offers a valuable resource for genetics, biology and medical researchers. It is freely accessible at https://biit.cs.ut.ee/gprofiler.

Systematic visualisation of molecular QTLs reveals variant mechanisms at GWAS loci.

Splicing quantitative trait loci (QTLs) have been implicated as a common mechanism underlying complex trait associations. However, utilising splicing QTLs in target discovery and prioritisation has been challenging due to extensive data normalisation which often renders the direction of the genetic effect as well as its magnitude difficult to interpret. This is further complicated by the fact that strong expression QTLs often manifest as weak splicing QTLs and vice versa, making it difficult to uniquely identify the underlying molecular mechanism at each locus. We find that these ambiguities can be mitigated by visualising the association between the genotype and average RNA sequencing read coverage in the region. Here, we generate these QTL coverage plots for 1.7 million molecular QTL associations in the eQTL Catalogue identified with five quantification methods. We illustrate the utility of these QTL coverage plots by performing colocalisation between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. We find that while visually confirmed splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases. All our association summary statistics and QTL coverage plots are freely available at https://www.ebi.ac.uk/eqtl/.

Neuroanatomy of the crocodylomorph Portugalosuchus azenhae from the late cretaceous of Portugal.

We present the first detailed braincase anatomical description and neuroanatomical study of Portugalosuchus azenhae, from the Cenomanian (Late Cretaceous) of Portugal. This eusuchian crocodylomorph was originally described as a putative Crocodylia and one of the oldest representatives of this clade; however, its phylogenetic position remains controversial. Based on new data obtained from high resolution Computed Tomography images (by micro-CT scan), this study aims to improve the original description of this taxon and also update the scarce neuroanatomical knowledge of Eusuchia and Crocodylia from this time interval, a key period to understand the origin and evolution of these clades. The resulting three-dimensional models from the CT data allowed a detailed description of its well-preserved neurocranium and internal cavities. Therefore, it was possible to reconstruct the cavities of the olfactory region, nasopharyngeal ducts, brain, nerves, carotid arteries, blood vessels, paratympanic sinus system and inner ear, which allowed to estimate some neurosensorial capabilities. By comparison with other crocodylomorphs, these analyses showed that Portugalosuchus, back in the Cenomanian, already displayed an olfactive acuity, sight, hearing and cognitive skills within the range of that observed in other basal eusuchians and crocodylians, including extant species. In addition, and in order to test its disputed phylogenetic position, these new anatomical data, which helped to correct and complete some of the original observations, were included in one of the most recent morphology-based phylogenies. The position of Portugalosuchus differs slightly from the original publication since it is now located as a "thoracosaurid" within Gavialoidea, but still as a crocodylian. Despite all this, to better contrast these results, additional phylogenetic analyses including this new morphological character coding together with DNA data should be performed.

Prediction of prime editing insertion efficiencies using sequence features and DNA repair determinants.

Most short sequences can be precisely written into a selected genomic target using prime editing; however, it remains unclear what factors govern insertion. We design a library of 3,604 sequences of various lengths and measure the frequency of their insertion into four genomic sites in three human cell lines, using different prime editor systems in varying DNA repair contexts. We find that length, nucleotide composition and secondary structure of the insertion sequence all affect insertion rates. We also discover that the 3' flap nucleases TREX1 and TREX2 suppress the insertion of longer sequences. Combining the sequence and repair features into a machine learning model, we can predict relative frequency of insertions into a site with R = 0.70. Finally, we demonstrate how our accurate prediction and user-friendly software help choose codon variants of common fusion tags that insert at high efficiency, and provide a catalog of empirically determined insertion rates for over a hundred useful sequences.

ICTV Virus Taxonomy Profile: Rhabdoviridae 2022.

The family Rhabdoviridae comprises viruses with negative-sense (-) RNA genomes of 10-16 kb. Virions are typically enveloped with bullet-shaped or bacilliform morphology but can also be non-enveloped filaments. Rhabdoviruses infect plants or animals, including mammals, birds, reptiles, amphibians or fish, as well as arthropods, which serve as single hosts or act as biological vectors for transmission to animals or plants. Rhabdoviruses include important pathogens of humans, livestock, fish or agricultural crops. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Rhabdoviridae, which is available at ictv.global/report/rhabdoviridae.

A single intranasal dose of human parainfluenza virus type 3-vectored vaccine induces effective antibody and memory T cell response in the lungs and protects hamsters against SARS-CoV-2.

Respiratory tract vaccination has an advantage of needle-free delivery and induction of mucosal immune response in the portal of SARS-CoV-2 entry. We utilized human parainfluenza virus type 3 vector to generate constructs expressing the full spike (S) protein of SARS-CoV-2, its S1 subunit, or the receptor-binding domain, and tested them in hamsters as single-dose intranasal vaccines. The construct bearing full-length S induced high titers of neutralizing antibodies specific to S protein domains critical to the protein functions. Robust memory T cell responses in the lungs were also induced, which represent an additional barrier to infection and should be less sensitive than the antibody responses to mutations present in SARS-CoV-2 variants. Following SARS-CoV-2 challenge, animals were protected from the disease and detectable viral replication. Vaccination prevented induction of gene pathways associated with inflammation. These results indicate advantages of respiratory vaccination against COVID-19 and inform the design of mucosal SARS-CoV-2 vaccines.

A compendium of uniformly processed human gene expression and splicing quantitative trait loci.

Many gene expression quantitative trait locus (eQTL) studies have published their summary statistics, which can be used to gain insight into complex human traits by downstream analyses, such as fine mapping and co-localization. However, technical differences between these datasets are a barrier to their widespread use. Consequently, target genes for most genome-wide association study (GWAS) signals have still not been identified. In the present study, we present the eQTL Catalogue ( https://www.ebi.ac.uk/eqtl ), a resource of quality-controlled, uniformly re-computed gene expression and splicing QTLs from 21 studies. We find that, for matching cell types and tissues, the eQTL effect sizes are highly reproducible between studies. Although most QTLs were shared between most bulk tissues, we identified a greater diversity of cell-type-specific QTLs from purified cell types, a subset of which also manifested as new disease co-localizations. Our summary statistics are freely available to enable the systematic interpretation of human GWAS associations across many cell types and tissues.

Braincase anatomy of extant Crocodylia, with new insights into the development and evolution of the neurocranium in crocodylomorphs.

Present-day crocodylians exhibit a remarkably akinetic skull with a highly modified braincase. We present a comprehensive description of the neurocranial osteology of extant crocodylians, with notes on the development of individual skeletal elements and a discussion of the terminology used for this project. The quadrate is rigidly fixed by multiple contacts with most braincase elements. The parabasisphenoid is sutured to the pterygoids (palate) and the quadrate (suspensorium); as a result, the basipterygoid joint is completely immobilized. The prootic is reduced and externally concealed by the quadrate. It has a verticalized buttress that participates in the canal for the temporal vasculature. The ventrolateral processes of the otoccipitals completely cover the posteroventral region of the braincase, enclose the occipital nerves and blood vessels in narrow bony canals and also provide additional sutural contacts between the braincase elements and further consolidate the posterior portion of the crocodylian skull. The otic capsule of crocodylians has a characteristic cochlear prominence that corresponds to the lateral route of the perilymphatic sac. Complex internal structures of the otoccipital (extracapsular buttress) additionally arrange the neurovascular structures of the periotic space of the cranium. Most of the braincase elements of crocodylians are excavated by the paratympanic pneumatic sinuses. The braincase in various extant crocodylians has an overall similar structure with some consistent variation between taxa. Several newly observed features of the braincase are present in Gavialis gangeticus and extant members of Crocodylidae to the exclusion of alligatorids: the reduced exposure of the prootic buttress on the floor of the temporal canal, the sagittal nuchal crest of the supraoccipital projecting posteriorly beyond the postoccipital processes and the reduced paratympanic pneumaticity. The most distinctive features of the crocodylian braincase (fixed quadrate and basipterygoid joint, consolidated occiput) evolved relatively rapidly at the base of Crocodylomorpha and accompanied the initial diversification of this clade during the Late Triassic and Early Jurassic. We hypothesize that profound rearrangements in the individual development of the braincases of basal crocodylomorphs underlie these rapid evolutionary modifications. These rearrangements are likely reflected in the embryonic development of extant crocodylians and include the involvement of neomorphic dermal anlagen in different portions of the developing chondrocranium, the extensive ossification of the palatoquadrate cartilage as a single expanded quadrate and the anteromedial inclination of the quadrate.

Wear patterns and dental functioning in an Early Cretaceous stegosaur from Yakutia, Eastern Russia.

Isolated stegosaurian teeth from the Early Cretaceous high-latitude (palaeolatitude estimate of N 62°- 66.5°) Teete locality in Yakutia (Eastern Siberia, Russia) are characterized by a labiolingually compressed, slightly asymmetrical and mesiodistally denticulated (9-14 denticles) crown, a pronounced ring-like cingulum, as well as a "complex network of secondary ridges". The 63 teeth (found during on-site excavation in 2012, 2017-2019 and screen-washing in 2017-2019) most likely belong to one species of a derived (stegosaurine) stegosaur. Most of the teeth exhibit a high degree of wear and up to three wear facets has been observed on a single tooth. The prevalence of worn teeth with up to three wear facets and the presence of different types of facets (including steeply inclined and groove-like) indicate the tooth-tooth contact and precise dental occlusion in the Teete stegosaur. The microwear pattern (mesiodistally or slightly obliquely oriented scratches; differently oriented straight and curved scratches on some wear facets) suggest a complex jaw mechanism with palinal jaw motion. Histological analysis revealed that the Teete stegosaur is characterized by relatively short tooth formation time (95 days) and the presence of a "wavy enamel pattern". Discoveries of a "wavy enamel pattern" in the Teete stegosaur, in a Middle Jurassic stegosaur from Western Siberia, and in the basal ceratopsian Psittacosaurus, suggest that this histological feature is common for different ornithischian clades, including ornithopods, marginocephalians, and thyreophorans. A juvenile tooth in the Teete sample indicates that stegosaurs were year-round residents and reproduced in high latitudes. The combination of high degree of tooth wear with formation of multiple wear facets, complex jaw motions, relatively short tooth formation time and possibly high tooth replacement rates is interpreted as a special adaptation for a life in high-latitude conditions or, alternatively, as a common stegosaurian adaptation making stegosaurs a successful group of herbivorous dinosaurs in the Middle Jurassic-Early Cretaceous and enabeling them to live in both low- and high-latitude ecosystems.

Effects of Overexpression of the Egyptian Fruit Bat Innate Immune Genes on Filovirus Infections in the Host Cells.

Bats constitute a large and diverse group of mammals with unique characteristics. One of these is the ability of bats to maintain various pathogens, particularly viruses, without evidence of disease. The innate immune system has been implicated as one of the important components involved in this process. However, in contrast to the human innate immune system, little data is available for bats. In the present study we generated 23 fusion constructs of innate immune genes of Egyptian fruit bat (Rousettus aegyptiacus) with mCherry as a fluorescent reporter. We evaluated the effects of overexpressing these genes on the replication of Marburg and Ebola viruses in the Egyptian fruit bat cell line R06EJ. Both viruses were substantially inhibited by overexpression of type I, II and III interferons, as well as by DDX58 (RIG-I), IFIH1, and IRF1. Our observations suggest that the broad antiviral activity of these genes reported previously in human cells is conserved in Egyptian fruit bats and these possess anti-filovirus activities that may contribute to the efficient virus clearance.

PAWER: protein array web exploreR.

BACKGROUND: Protein microarray is a well-established approach for characterizing activity levels of thousands of proteins in a parallel manner. Analysis of protein microarray data is complex and time-consuming, while existing solutions are either outdated or challenging to use without programming skills. The typical data analysis pipeline consists of a data preprocessing step, followed by differential expression analysis, which is then put into context via functional enrichment. Normally, biologists would need to assemble their own workflow by combining a set of unrelated tools to analyze experimental data. Provided that most of these tools are developed independently by various bioinformatics groups, making them work together could be a real challenge. RESULTS: Here we present PAWER, the online web tool dedicated solely to protein microarray analysis. PAWER enables biologists to carry out all the necessary analysis steps in one go. PAWER provides access to state-of-the-art computational methods through the user-friendly interface, resulting in publication-ready illustrations. We also provide an R package for more advanced use cases, such as bespoke analysis workflows. CONCLUSIONS: PAWER is freely available at https://biit.cs.ut.ee/pawer .

Species-Specific Evolution of Ebola Virus during Replication in Human and Bat Cells.

Ebola virus (EBOV) causes a severe, often fatal disease in humans and nonhuman primates. Within the past decade, EBOV has caused two large and difficult-to-control outbreaks, one of which recently ended in the Democratic Republic of the Congo. Bats are the likely reservoir of EBOV, but little is known of their relationship with the virus. We perform serial passages of EBOV in human and bat cells and use circular sequencing to compare the short-term evolution of the virus. Virus populations passaged in bat cells have sequence markers indicative of host RNA editing enzyme activity, including evidence for ADAR editing of the EBOV glycoprotein. Multiple regions in the EBOV genome appear to have undergone adaptive evolution when passaged in bat and human cells. Individual mutated viruses are rescued and characterized. Our results provide insight into the host species-specific evolution of EBOV and highlight the adaptive flexibility of the virus.

gprofiler2 -- an R package for gene list functional enrichment analysis and namespace conversion toolset g:Profiler.

g:Profiler ( https://biit.cs.ut.ee/gprofiler) is a widely used gene list functional profiling and namespace conversion toolset that has been contributing to reproducible biological data analysis already since 2007. Here we introduce the accompanying R package, gprofiler2, developed to facilitate programmatic access to g:Profiler computations and databases via REST API. The gprofiler2 package provides an easy-to-use functionality that enables researchers to incorporate functional enrichment analysis into automated analysis pipelines written in R. The package also implements interactive visualisation methods to help to interpret the enrichment results and to illustrate them for publications. In addition, gprofiler2 gives access to the versatile gene/protein identifier conversion functionality in g:Profiler enabling to map between hundreds of different identifier types or orthologous species. The gprofiler2 package is freely available at the CRAN repository.

Taxonomy of the order Mononegavirales: update 2019.

In February 2019, following the annual taxon ratification vote, the order Mononegavirales was amended by the addition of four new subfamilies and 12 new genera and the creation of 28 novel species. This article presents the updated taxonomy of the order Mononegavirales as now accepted by the International Committee on Taxonomy of Viruses (ICTV).

g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update).

Biological data analysis often deals with lists of genes arising from various studies. The g:Profiler toolset is widely used for finding biological categories enriched in gene lists, conversions between gene identifiers and mappings to their orthologs. The mission of g:Profiler is to provide a reliable service based on up-to-date high quality data in a convenient manner across many evidence types, identifier spaces and organisms. g:Profiler relies on Ensembl as a primary data source and follows their quarterly release cycle while updating the other data sources simultaneously. The current update provides a better user experience due to a modern responsive web interface, standardised API and libraries. The results are delivered through an interactive and configurable web design. Results can be downloaded as publication ready visualisations or delimited text files. In the current update we have extended the support to 467 species and strains, including vertebrates, plants, fungi, insects and parasites. By supporting user uploaded custom GMT files, g:Profiler is now capable of analysing data from any organism. All past releases are maintained for reproducibility and transparency. The 2019 update introduces an extensive technical rewrite making the services faster and more flexible. g:Profiler is freely available at https://biit.cs.ut.ee/gprofiler.

Priorities in applied research to ensure programmatic success in the global elimination of canine rabies.

The elimination of human rabies mediated by dogs is attainable in concept, based upon current sensitive and specific diagnostic methods, existing safe and effective human and veterinary vaccines and a sound virological, pathological and epidemiological understanding of the disease. Globally, all developed countries achieved this goal. Regionally, major progress occurred throughout the Americas. However, less advancement is evident in Africa and Asia. Our objective was to concentrate upon those salient improvements to extant tools and methods over the next five years which could assist and simplify the task for both those developing countries that have already begun the process, as well as other localities in the earlier stages of consideration. We considered several categories of applied research which could be accomplished in the short term, based upon the available scientific evidence and recent recommendations from subject matter experts and key opinion leaders, focused upon perceived major limitations to prior program success. Areas of concentration included: laboratory-based surveillance, pathogen detection and characterization; human rabies prophylaxis; veterinary biologics; implementation of canine vaccination; and oral vaccination of free-ranging community dogs. Further real-time application in these core areas with proven techniques and technology would simplify attaining not only the global goal focused subtly upon human mortality, but the actual elimination of canine rabies as well.