The Egyptian Rousette Genome Reveals Unexpected Features of Bat Antiviral Immunity.

Bats harbor many viruses asymptomatically, including several notorious for causing extreme virulence in humans. To identify differences between antiviral mechanisms in humans and bats, we sequenced, assembled, and analyzed the genome of Rousettus aegyptiacus, a natural reservoir of Marburg virus and the only known reservoir for any filovirus. We found an expanded and diversified KLRC/KLRD family of natural killer cell receptors, MHC class I genes, and type I interferons, which dramatically differ from their functional counterparts in other mammals. Such concerted evolution of key components of bat immunity is strongly suggestive of novel modes of antiviral defense. An evaluation of the theoretical function of these genes suggests that an inhibitory immune state may exist in bats. Based on our findings, we hypothesize that tolerance of viral infection, rather than enhanced potency of antiviral defenses, may be a key mechanism by which bats asymptomatically host viruses that are pathogenic in humans.

"Super-Spreaders" and Person-to-Person Transmission of Andes Virus in Argentina.

BACKGROUND: From November 2018 through February 2019, person-to-person transmission of Andes virus (ANDV) hantavirus pulmonary syndrome occurred in Chubut Province, Argentina, and resulted in 34 confirmed infections and 11 deaths. Understanding the genomic, epidemiologic, and clinical characteristics of person-to-person transmission of ANDV is crucial to designing effective interventions. METHODS: Clinical and epidemiologic information was obtained by means of patient report and from public health centers. Serologic testing, contact-tracing, and next-generation sequencing were used to identify ANDV infection as the cause of this outbreak of hantavirus pulmonary syndrome and to reconstruct person-to-person transmission events. RESULTS: After a single introduction of ANDV from a rodent reservoir into the human population, transmission was driven by 3 symptomatic persons who attended crowded social events. After 18 cases were confirmed, public health officials enforced isolation of persons with confirmed cases and self-quarantine of possible contacts; these measures most likely curtailed further spread. The median reproductive number (the number of secondary cases caused by an infected person during the infectious period) was 2.12 before the control measures were enforced and decreased to 0.96 after the measures were implemented. Full genome sequencing of the ANDV strain involved in this outbreak was performed with specimens from 27 patients and showed that the strain that was present (Epuyén/18-19) was similar to the causative strain (Epilink/96) in the first known person-to-person transmission of hantavirus pulmonary syndrome caused by ANDV, which occurred in El Bolsón, Argentina, in 1996. Clinical investigations involving patients with ANDV hantavirus pulmonary syndrome in this outbreak revealed that patients with a high viral load and liver injury were more likely than other patients to spread infection. Disease severity, genomic diversity, age, and time spent in the hospital had no clear association with secondary transmission. CONCLUSIONS: Among patients with ANDV hantavirus pulmonary syndrome, high viral titers in combination with attendance at massive social gatherings or extensive contact among persons were associated with a higher likelihood of transmission. (Funded by the Ministerio de Salud y Desarrollo Social de la Nación Argentina and others.).

Virus genomes reveal factors that spread and sustained the Ebola epidemic.

The 2013-2016 West African epidemic caused by the Ebola virus was of unprecedented magnitude, duration and impact. Here we reconstruct the dispersal, proliferation and decline of Ebola virus throughout the region by analysing 1,610 Ebola virus genomes, which represent over 5% of the known cases. We test the association of geography, climate and demography with viral movement among administrative regions, inferring a classic 'gravity' model, with intense dispersal between larger and closer populations. Despite attenuation of international dispersal after border closures, cross-border transmission had already sown the seeds for an international epidemic, rendering these measures ineffective at curbing the epidemic. We address why the epidemic did not spread into neighbouring countries, showing that these countries were susceptible to substantial outbreaks but at lower risk of introductions. Finally, we reveal that this large epidemic was a heterogeneous and spatially dissociated collection of transmission clusters of varying size, duration and connectivity. These insights will help to inform interventions in future epidemics.

Molecular Evidence of Sexual Transmission of Ebola Virus.

A suspected case of sexual transmission from a male survivor of Ebola virus disease (EVD) to his female partner (the patient in this report) occurred in Liberia in March 2015. Ebola virus (EBOV) genomes assembled from blood samples from the patient and a semen sample from the survivor were consistent with direct transmission. The genomes shared three substitutions that were absent from all other Western African EBOV sequences and that were distinct from the last documented transmission chain in Liberia before this case. Combined with epidemiologic data, the genomic analysis provides evidence of sexual transmission of EBOV and evidence of the persistence of infective EBOV in semen for 179 days or more after the onset of EVD. (Funded by the Defense Threat Reduction Agency and others.).

Cynomolgus macaque (Macaca fascicularis) immunoglobulin heavy chain locus description.

Cynomolgus macaques (Macaca fascicularis) have become an important animal model for biomedical research. In particular, it is the animal model of choice for the development of vaccine candidates associated with emerging dangerous pathogens. Despite their increasing importance as animal models, the cynomolgus macaque genome is not fully characterized, hindering molecular studies for this model. More importantly, the lack of knowledge about the immunoglobulin (IG) locus organization directly impacts the analysis of the humoral response in cynomolgus macaques. Recent advances in next generation sequencing (NGS) technologies to analyze IG repertoires open the opportunity to deeply characterize the humoral immune response. However, the IG locus organization for the animal is required to completely dissect IG repertoires. Here, we describe the localization and organization of the rearranging IG heavy (IGH) genes on chromosome 7 of the cynomolgus macaque draft genome. Our annotation comprises 108 functional genes which include 63 variable (IGHV), 38 diversity (IGHD), and 7 joining (IGHJ) genes. For validation, we provide RNA transcript data for most of the IGHV genes and all of the annotated IGHJ genes, as well as proteomic data to validate IGH constant genes. The description and annotation of the rearranging IGH genes for the cynomolgus macaques will significantly facilitate scientific research. This is particularly relevant to dissect the immune response during vaccination or infection with dangerous pathogens such as Ebola, Marburg and other emerging pathogens where non-human primate models play a significant role for countermeasure development.

De novo transcriptome reconstruction and annotation of the Egyptian rousette bat.

BACKGROUND: The Egyptian Rousette bat (Rousettus aegyptiacus), a common fruit bat species found throughout Africa and the Middle East, was recently identified as a natural reservoir host of Marburg virus. With Ebola virus, Marburg virus is a member of the family Filoviridae that causes severe hemorrhagic fever disease in humans and nonhuman primates, but results in little to no pathological consequences in bats. Understanding host-pathogen interactions within reservoir host species and how it differs from hosts that experience severe disease is an important aspect of evaluating viral pathogenesis and developing novel therapeutics and methods of prevention. RESULTS: Progress in studying bat reservoir host responses to virus infection is hampered by the lack of host-specific reagents required for immunological studies. In order to establish a basis for the design of reagents, we sequenced, assembled, and annotated the R. aegyptiacus transcriptome. We performed de novo transcriptome assembly using deep RNA sequencing data from 11 distinct tissues from one male and one female bat. We observed high similarity between this transcriptome and those available from other bat species. Gene expression analysis demonstrated clustering of expression profiles by tissue, where we also identified enrichment of tissue-specific gene ontology terms. In addition, we identified and experimentally validated the expression of novel coding transcripts that may be specific to this species. CONCLUSION: We comprehensively characterized the R. aegyptiacus transcriptome de novo. This transcriptome will be an important resource for understanding bat immunology, physiology, disease pathogenesis, and virus transmission.

Monitoring of Ebola Virus Makona Evolution through Establishment of Advanced Genomic Capability in Liberia.

To support Liberia's response to the ongoing Ebola virus (EBOV) disease epidemic in Western Africa, we established in-country advanced genomic capabilities to monitor EBOV evolution. Twenty-five EBOV genomes were sequenced at the Liberian Institute for Biomedical Research, which provided an in-depth view of EBOV diversity in Liberia during September 2014-February 2015. These sequences were consistent with a single virus introduction to Liberia; however, shared ancestry with isolates from Mali indicated at least 1 additional instance of movement into or out of Liberia. The pace of change is generally consistent with previous estimates of mutation rate. We observed 23 nonsynonymous mutations and 1 nonsense mutation. Six of these changes are within known binding sites for sequence-based EBOV medical countermeasures; however, the diagnostic and therapeutic impact of EBOV evolution within Liberia appears to be low.

Orthohantavirus diversity in Central-East Argentina: Insights from complete genomic sequencing on phylogenetics, Geographic patterns and transmission scenarios.

Hantavirus Pulmonary Syndrome (HPS), characterized by its high fatality rate, poses a significant public health concern in Argentina due to the increasing evidence of person-to-person transmission of Andes virus. Several orthohantaviruses were described in the country, but their phylogenetic relationships were inferred from partial genomic sequences. The objectives of this work were to assess the viral diversity of the most prevalent orthohantaviruses associated with HPS cases in the Central-East (CE) region of Argentina, elucidate the geographic patterns of distribution of each variant and reconstruct comprehensive phylogenetic relationships utilizing complete genomic sequencing. To accomplish this, a detailed analysis was conducted of the geographic distribution of reported cases within the most impacted province of the region. A representative sample of cases was then selected to generate a geographic map illustrating the distribution of viral variants. Complete viral genomes were obtained from HPS cases reported in the region, including some from epidemiologically linked cases. The phylogenetic analysis based on complete genomes defined two separate clades in Argentina: Andes virus in the Southwestern region and Andes-like viruses in other parts of the country. In the CE region, Buenos Aires virus and Lechiguanas virus clearly segregate in two subclades. Complete genomes were useful to distinguish person-to-person transmission from environmental co-exposure to rodent population. This study enhances the understanding of the genetic diversity, geographical spread, and transmission dynamics of orthohantaviruses in Central Argentina and prompt to consider the inclusion of Buenos Aires virus and Lechiguanas virus in the species Orthohantavirus andesense, as named viruses.

Monkeypox virus genomic accordion strategies.

The 2023 monkeypox (mpox) epidemic was caused by a subclade IIb descendant of a monkeypox virus (MPXV) lineage traced back to Nigeria in 1971. Person-to-person transmission appears higher than for clade I or subclade IIa MPXV, possibly caused by genomic changes in subclade IIb MPXV. Key genomic changes could occur in the genome's low-complexity regions (LCRs), which are challenging to sequence and are often dismissed as uninformative. Here, using a combination of highly sensitive techniques, we determine a high-quality MPXV genome sequence of a representative of the current epidemic with LCRs resolved at unprecedented accuracy. This reveals significant variation in short tandem repeats within LCRs. We demonstrate that LCR entropy in the MPXV genome is significantly higher than that of single-nucleotide polymorphisms (SNPs) and that LCRs are not randomly distributed. In silico analyses indicate that expression, translation, stability, or function of MPXV orthologous poxvirus genes (OPGs), including OPG153, OPG204, and OPG208, could be affected in a manner consistent with the established "genomic accordion" evolutionary strategies of orthopoxviruses. We posit that genomic studies focusing on phenotypic MPXV differences should consider LCR variability.

Cutting edge: A role for inside-out signaling in TCR regulation of CD28 ligand binding.

Efficient T cell activation depends on the engagement of both TCR and CD28, although the molecular mechanisms that control this signal integration are not fully understood. Using fluorescence resonance energy transfer, we show that T cell activation can drive a reorientation of the cytosolic tails of the CD28 dimer. However, this is not mediated through CD28 ligand binding. Rather, TCR signaling itself mediates this conformation change in CD28. We also show that TCR signaling can induce CD28-ligand interactions. Although the CD28 dimer appears to bind ligand monovalently in solution, we show that both ligand binding sites are required to efficiently recruit CD28 to the immunological synapse. These results suggest, that analogous to the cross-talk from TCR that regulates integrin activation, TCR-initiated inside-out signaling may induce a conformational change to the extracellular domains of CD28, enabling ligand binding and initiating CD28 signaling.

Induced protection from a CCHFV-M DNA vaccine requires CD8+ T cells.

Crimean-Congo hemorrhagic fever (CCHF) is a World Health Organization prioritized disease because its broad distribution and severity of disease make it a global health threat. Despite advancements in preclinical vaccine development for CCHF virus (CCHFV), including multiple platforms targeting multiple antigens, a clear definition of the adaptive immune correlates of protection is lacking. Levels of neutralizing antibodies in vaccinated animal models do not necessarily correlate with protection, suggesting that cellular immunity, such as CD8+ T cells, might have an important role in protection in this model. Using a well-established IFN-I antibody blockade mouse model (IS) and a DNA-based vaccine encoding the CCHFV M-segment glycoprotein precursor, we investigated the role of humoral and T cell immunity in vaccine-mediated protection in mice genetically devoid of these immune compartments. We found that in the absence of the B-cell compartment (µMT knockout mice), protection provided by the vaccine was not reduced. In contrast, in the absence of CD8+ T cells (CD8+ knockout mice) the vaccine-mediated protection was significantly diminished. Importantly, humoral responses to the vaccine in CD8+ T-cell knockout mice were equivalent to wild-type mice. These findings indicated that CD8+ T-cell responses are necessary and sufficient to promote protection in mice vaccinated with the M-segment DNA vaccine. Identifying a crucial role of the cellular immunity to protect against CCHFV should help guide the development of CCHFV-targeting vaccines.

Andes Virus Genome Mutations That Are Likely Associated with Animal Model Attenuation and Human Person-to-Person Transmission.

We performed whole-genome sequencing with bait enrichment techniques to analyze Andes virus (ANDV), a cause of human hantavirus pulmonary syndrome. We used cryopreserved lung tissues from a naturally infected long-tailed colilargo, including early, intermediate, and late cell culture, passages of an ANDV isolate from that animal, and lung tissues from golden hamsters experimentally exposed to that ANDV isolate. The resulting complete genome sequences were subjected to detailed comparative genomic analysis against American orthohantaviruses. We identified four amino acid substitutions related to cell culture adaptation that resulted in attenuation of ANDV in the typically lethal golden hamster animal model of hantavirus pulmonary syndrome. Changes in the ANDV nucleocapsid protein, glycoprotein, and small nonstructural protein open reading frames correlated with mutations typical for ANDV strains associated with increased virulence in the small-animal model. Finally, we identified three amino acid substitutions, two in the small nonstructural protein and one in the glycoprotein, that were only present in the clade of viruses associated with efficient person-to-person transmission. Our results indicate that there are single-nucleotide polymorphisms that could be used to predict strain-specific ANDV virulence and/or transmissibility. IMPORTANCE Several orthohantaviruses cause the zoonotic disease hantavirus pulmonary syndrome (HPS) in the Americas. Among them, HPS caused by Andes virus (ANDV) is of great public health concern because it is associated with the highest case fatality rate (up to 50%). ANDV is also the only orthohantavirus associated with relatively robust evidence of person-to-person transmission. This work reveals nucleotide changes in the ANDV genome that are associated with virulence attenuation in an animal model and increased transmissibility in humans. These findings may pave the way to early severity predictions in future ANDV-caused HPS outbreaks.

Peripheral immune responses to filoviruses in a reservoir versus spillover hosts reveal transcriptional correlates of disease.

Several filoviruses, including Marburg virus (MARV), cause severe disease in humans and nonhuman primates (NHPs). However, the Egyptian rousette bat (ERB, Rousettus aegyptiacus), the only known MARV reservoir, shows no overt illness upon natural or experimental infection, which, like other bat hosts of zoonoses, is due to well-adapted, likely species-specific immune features. Despite advances in understanding reservoir immune responses to filoviruses, ERB peripheral blood responses to MARV and how they compare to those of diseased filovirus-infected spillover hosts remain ill-defined. We thus conducted a longitudinal analysis of ERB blood gene responses during acute MARV infection. These data were then contrasted with a compilation of published primate blood response studies to elucidate gene correlates of filovirus protection versus disease. Our work expands on previous findings in MARV-infected ERBs by supporting both host resistance and disease tolerance mechanisms, offers insight into the peripheral immunocellular repertoire during infection, and provides the most direct known cross-examination between reservoir and spillover hosts of the most prevalently-regulated response genes, pathways and activities associated with differences in filovirus pathogenesis and pathogenicity.

Delayed viral clearance despite high number of activated T cells during the acute phase in Argentinean patients with hantavirus pulmonary syndrome.

BACKGROUND: The hallmarks of HPS are increase of vascular permeability and endothelial dysfunction. Although an exacerbated immune response is thought to be implicated in pathogenesis, clear evidence is still elusive. As orthohantaviruses are not cytopathic CD8+ T cells are believed to be the central players involved in pathogenesis. METHODS: Serum and blood samples from Argentinean HPS patients were collected from 2014 to 2019. Routine white blood cell analyses, quantification and characterization of T-cell phenotypic profile, viral load, neutralizing antibody response and quantification of inflammatory mediators were performed. FINDINGS: High numbers of activated CD4+ and CD8+ T cells were found in all HPS cases independently of disease severity. We found increased levels of some proinflammatory mediators during the acute phase of illness. Nonetheless, viral RNA remained high, showing a delay in clearance from blood up to late convalescence, when titers of neutralizing antibodies reached a high level. INTERPRETATION: The high activated phenotypic profile of T cells seems to be unable to resolve infection during the acute and early convalescent phases, and it was not associated with the severity of the disease. Thus, at least part of the activated T cells could be induced by the dysregulated inflammatory response in an unspecific manner. Viral clearance seems to have been more related to high titers of neutralizing antibodies than to the T-cell response. FUNDING: This work was supported mainly by the Administración Nacional de Laboratorios e Institutos de Salud (ANLIS) "Dr. Carlos Malbrán". Further details of fundings sources is included in the appendix.

Genomic features of humoral immunity support tolerance model in Egyptian rousette bats.

Bats asymptomatically harbor many viruses that can cause severe human diseases. The Egyptian rousette bat (ERB) is the only known reservoir for Marburgviruses and Sosuga virus, making it an exceptional animal model to study antiviral mechanisms in an asymptomatic host. With this goal in mind, we constructed and annotated the immunoglobulin heavy chain locus, finding an expansion on immunoglobulin variable genes associated with protective human antibodies to different viruses. We also annotated two functional and distinct immunoglobulin epsilon genes and four distinctive functional immunoglobulin gamma genes. We described the Fc receptor repertoire in ERBs, including features that may affect activation potential, and discovered the lack of evolutionary conserved short pentraxins. These findings reinforce the hypothesis that a differential threshold of regulation and/or absence of key immune mediators may promote tolerance and decrease inflammation in ERBs.

Molecular analysis of the 2012 Bundibugyo virus disease outbreak.

Bundibugyo virus (BDBV) is one of four ebolaviruses known to cause disease in humans. Bundibugyo virus disease (BVD) outbreaks occurred in 2007-2008 in Bundibugyo District, Uganda, and in 2012 in Isiro, Province Orientale, Democratic Republic of the Congo. The 2012 BVD outbreak resulted in 38 laboratory-confirmed cases of human infection, 13 of whom died. However, only 4 BDBV specimens from the 2012 outbreak have been sequenced. Here, we provide BDBV sequences from seven additional patients. Analysis of the molecular epidemiology and evolutionary dynamics of the 2012 outbreak with these additional isolates challenges the current hypothesis that the outbreak was the result of a single spillover event. In addition, one patient record indicates that BDBV's initial emergence in Isiro occurred 50 days earlier than previously accepted. Collectively, this work demonstrates how retrospective sequencing can be used to elucidate outbreak origins and provide epidemiological contexts to a medically relevant pathogen.

Novel viruses in hard ticks collected in the Republic of Korea unveiled by metagenomic high-throughput sequencing analysis.

Ticks are vectors of a wide range of zoonotic viruses of medical and veterinary importance. Recently, metagenomics studies demonstrated that they are also the source of potentially pathogenic novel viruses. During the period from 2015 to 2017, questing ticks were collected by dragging the vegetation from geographically distant locations in the Republic of Korea (ROK) and a target-independent high-throughput sequencing method was utilized to study their virome. A total of seven viruses, including six putative novel viral entities, were identified. Genomic analysis showed that the novel viruses were most closely related to members in the orders Jingchuvirales and Bunyavirales. Phylogenetic reconstruction showed that the Bunyavirales-like viruses grouped in the same clade with other viruses within the Nairovirus and Phlebovirus genera, while the novel Jingchuvirales-like virus grouped together with other viruses within the family Chuviridae. Real-time RT-PCR was used to determine the geographic distribution and prevalence of these viruses in adult ticks. These novel viruses have a wide geographic distribution in the ROK with prevalences ranging from 2% to 18%. Our study expands the knowledge about the composition of the tick virome and highlights the wide diversity of viruses they harbor in the ROK. The discovery of novel viruses associated with ticks in the ROK highlights the need for an active tick-borne disease surveillance program to identify possible reservoirs of putative novel human pathogens.

Asymptomatic Infection of Marburg Virus Reservoir Bats Is Explained by a Strategy of Immunoprotective Disease Tolerance.

Marburg virus (MARV) is among the most virulent pathogens of primates, including humans. Contributors to severe MARV disease include immune response suppression and inflammatory gene dysregulation ("cytokine storm"), leading to systemic damage and often death. Conversely, MARV causes little to no clinical disease in its reservoir host, the Egyptian rousette bat (ERB). Previous genomic and in vitro data suggest that a tolerant ERB immune response may underlie MARV avirulence, but no significant examination of this response in vivo yet exists. Here, using colony-bred ERBs inoculated with a bat isolate of MARV, we use species-specific antibodies and an immune gene probe array (NanoString) to temporally characterize the transcriptional host response at sites of MARV replication relevant to primate pathogenesis and immunity, including CD14+ monocytes/macrophages, critical immune response mediators, primary MARV targets, and skin at the inoculation site, where highest viral loads and initial engagement of antiviral defenses are expected. Our analysis shows that ERBs upregulate canonical antiviral genes typical of mammalian systems, such as ISG15, IFIT1, and OAS3, yet demonstrate a remarkable lack of significant induction of proinflammatory genes classically implicated in primate filoviral pathogenesis, including CCL8, FAS, and IL6. Together, these findings offer the first in vivo functional evidence for disease tolerance as an immunological mechanism by which the bat reservoir asymptomatically hosts MARV. More broadly, these data highlight factors determining disparate outcomes between reservoir and spillover hosts and defensive strategies likely utilized by bat hosts of other emerging pathogens, knowledge that may guide development of effective antiviral therapies.

On-Demand Patient-Specific Phenotype-to-Genotype Ebola Virus Characterization.

Biosafety, biosecurity, logistical, political, and technical considerations can delay or prevent the wide dissemination of source material containing viable virus from the geographic origin of an outbreak to laboratories involved in developing medical countermeasures (MCMs). However, once virus genome sequence information is available from clinical samples, reverse-genetics systems can be used to generate virus stocks de novo to initiate MCM development. In this study, we developed a reverse-genetics system for natural isolates of Ebola virus (EBOV) variants Makona, Tumba, and Ituri, which have been challenging to obtain. These systems were generated starting solely with in silico genome sequence information and have been used successfully to produce recombinant stocks of each of the viruses for use in MCM testing. The antiviral activity of MCMs targeting viral entry varied depending on the recombinant virus isolate used. Collectively, selecting and synthetically engineering emerging EBOV variants and demonstrating their efficacy against available MCMs will be crucial for answering pressing public health and biosecurity concerns during Ebola disease (EBOD) outbreaks.

Signals and sequences that control CD28 localization to the central region of the immunological synapse.

During T cell interaction with APC, CD28 is recruited to the central region (cSMAC) of the immunological synapse. CD28-mediated signaling through PI3K results in the recruitment of protein kinase C-theta (PKCtheta) to the cSMAC, activation of NF-kappaB, and up-regulation of IL-2 transcription. However, the mechanism that mediates CD28 localization to the cSMAC and the functional consequences of CD28 localization to the cSMAC are not understood. In this report, we show that CD28 recruitment and persistence at the immunological synapse requires TCR signals and CD80 engagement. Addition of mAb to either MHC class II or CD80 results in the rapid displacement of CD28 from the immunological synapse. Ligand binding is not sufficient for CD28 localization to the immunological synapse, as truncation of the cytosolic tail of CD28 disrupts synapse localization without effecting the ability of CD28 to bind CD80. Furthermore, a single point mutation in the CD28 cytosolic tail (tyrosine 188) interferes with the ability of CD28 to preferentially accumulate at the cSMAC. PKCtheta distribution at the immunological synapse mirrors the distribution of tyrosine 188-mutated CD28, indicating that CD28 drives the localization of PKCtheta even when CD28 is not localized to the cSMAC. Mutation of tyrosine 188 also results in diminished activation of NF-kappaB, suggesting that CD28-mediated localization of PKCtheta to the cSMAC is important for efficient signal transduction. These data reinforce the importance of the interplay of signals between TCR and CD28 and suggest that CD28 signaling through PCKtheta may be mediated through localization to the cSMAC region of the immunological synapse.