Transmission from vaccinated individuals in a large SARS-CoV-2 Delta variant outbreak.

An outbreak of over 1,000 COVID-19 cases in Provincetown, Massachusetts (MA), in July 2021-the first large outbreak mostly in vaccinated individuals in the US-prompted a comprehensive public health response, motivating changes to national masking recommendations and raising questions about infection and transmission among vaccinated individuals. To address these questions, we combined viral genomic and epidemiological data from 467 individuals, including 40% of outbreak-associated cases. The Delta variant accounted for 99% of cases in this dataset; it was introduced from at least 40 sources, but 83% of cases derived from a single source, likely through transmission across multiple settings over a short time rather than a single event. Genomic and epidemiological data supported multiple transmissions of Delta from and between fully vaccinated individuals. However, despite its magnitude, the outbreak had limited onward impact in MA and the US overall, likely due to high vaccination rates and a robust public health response.

Benchmarking Metagenomics Tools for Taxonomic Classification.

Metagenomic sequencing is revolutionizing the detection and characterization of microbial species, and a wide variety of software tools are available to perform taxonomic classification of these data. The fast pace of development of these tools and the complexity of metagenomic data make it important that researchers are able to benchmark their performance. Here, we review current approaches for metagenomic analysis and evaluate the performance of 20 metagenomic classifiers using simulated and experimental datasets. We describe the key metrics used to assess performance, offer a framework for the comparison of additional classifiers, and discuss the future of metagenomic data analysis.

COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets.

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63+ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

Lassa Virus Genetics.

In a pattern repeated across a range of ecological niches, arenaviruses have evolved a compact four-gene genome to orchestrate a complex life cycle in a narrow range of susceptible hosts. A number of mammalian arenaviruses cross-infect humans, often causing a life-threatening viral hemorrhagic fever. Among this group of geographically bound zoonoses, Lassa virus has evolved a unique niche that leads to significant and sustained human morbidity and mortality. As a biosafety level 4 pathogen, direct study of the pathogenesis of Lassa virus is limited by the sparse availability, high operating costs, and technical restrictions of the high-level biocontainment laboratories required for safe experimentation. In this chapter, we introduce the relationship between genome structure and the life cycle of Lassa virus and outline reverse genetic approaches used to probe and describe functional elements of the Lassa virus genome. We then review the tools used to obtain viral genomic sequences used for phylogeny and molecular diagnostics, before shifting to a population perspective to assess the contributions of phylogenetic analysis in understanding the evolution and ecology of Lassa virus in West Africa. We finally consider the future outlook and clinical applications for genetic study of Lassa virus.

Distinguishing Severe Acute Respiratory Syndrome Coronavirus 2 Persistence and Reinfection: A Retrospective Cohort Study.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reinfection is poorly understood, partly because few studies have systematically applied genomic analysis to distinguish reinfection from persistent RNA detection related to initial infection. We aimed to evaluate the characteristics of SARS-CoV-2 reinfection and persistent RNA detection using independent genomic, clinical, and laboratory assessments. METHODS: All individuals at a large academic medical center who underwent a SARS-CoV-2 nucleic acid amplification test (NAAT) ≥45 days after an initial positive test, with both tests between 14 March and 30 December 2020, were analyzed for potential reinfection. Inclusion criteria required having ≥2 positive NAATs collected ≥45 days apart with a cycle threshold (Ct) value <35 at repeat testing. For each included subject, likelihood of reinfection was assessed by viral genomic analysis of all available specimens with a Ct value <35, structured Ct trajectory criteria, and case-by-case review by infectious diseases physicians. RESULTS: Among 1569 individuals with repeat SARS-CoV-2 testing ≥45 days after an initial positive NAAT, 65 (4%) met cohort inclusion criteria. Viral genomic analysis characterized mutations present and was successful for 14/65 (22%) subjects. Six subjects had genomically supported reinfection, and 8 subjects had genomically supported persistent RNA detection. Compared to viral genomic analysis, clinical and laboratory assessments correctly distinguished reinfection from persistent RNA detection in 12/14 (86%) subjects but missed 2/6 (33%) genomically supported reinfections. CONCLUSIONS: Despite good overall concordance with viral genomic analysis, clinical and Ct value-based assessments failed to identify 33% of genomically supported reinfections. Scaling-up genomic analysis for clinical use would improve detection of SARS-CoV-2 reinfections.

Combining genomics and epidemiology to track mumps virus transmission in the United States.

Unusually large outbreaks of mumps across the United States in 2016 and 2017 raised questions about the extent of mumps circulation and the relationship between these and prior outbreaks. We paired epidemiological data from public health investigations with analysis of mumps virus whole genome sequences from 201 infected individuals, focusing on Massachusetts university communities. Our analysis suggests continuous, undetected circulation of mumps locally and nationally, including multiple independent introductions into Massachusetts and into individual communities. Despite the presence of these multiple mumps virus lineages, the genomic data show that one lineage has dominated in the US since at least 2006. Widespread transmission was surprising given high vaccination rates, but we found no genetic evidence that variants arising during this outbreak contributed to vaccine escape. Viral genomic data allowed us to reconstruct mumps transmission links not evident from epidemiological data or standard single-gene surveillance efforts and also revealed connections between apparently unrelated mumps outbreaks.

Determining the Incidence of Asymptomatic SARS-CoV-2 Among Early Recipients of COVID-19 Vaccines (DISCOVER-COVID-19): A Prospective Cohort Study of Healthcare Workers Before, During and After Vaccination.

The impact of coronavirus disease 2019 vaccination on viral characteristics of breakthrough infections is unknown. In this prospective cohort study, incidence of severe acute respiratory syndrome coronavirus 2 infection decreased following vaccination. Although asymptomatic positive tests were observed following vaccination, the higher cycle thresholds, repeat negative tests, and inability to culture virus raise questions about their clinical significance.

Genomic Analysis of Lassa Virus during an Increase in Cases in Nigeria in 2018.

During 2018, an unusual increase in Lassa fever cases occurred in Nigeria, raising concern among national and international public health agencies. We analyzed 220 Lassa virus genomes from infected patients, including 129 from the 2017-2018 transmission season, to understand the viral populations underpinning the increase. A total of 14 initial genomes from 2018 samples were generated at Redeemer's University in Nigeria, and the findings were shared with the Nigerian Center for Disease Control in real time. We found that the increase in cases was not attributable to a particular Lassa virus strain or sustained by human-to-human transmission. Instead, the data were consistent with ongoing cross-species transmission from local rodent populations. Phylogenetic analysis also revealed extensive viral diversity that was structured according to geography, with major rivers appearing to act as barriers to migration of the rodent reservoir.

Outbreak of SARS-CoV-2 Infections, Including COVID-19 Vaccine Breakthrough Infections, Associated with Large Public Gatherings - Barnstable County, Massachusetts, July 2021.

During July 2021, 469 cases of COVID-19 associated with multiple summer events and large public gatherings in a town in Barnstable County, Massachusetts, were identified among Massachusetts residents; vaccination coverage among eligible Massachusetts residents was 69%. Approximately three quarters (346; 74%) of cases occurred in fully vaccinated persons (those who had completed a 2-dose course of mRNA vaccine [Pfizer-BioNTech or Moderna] or had received a single dose of Janssen [Johnson & Johnson] vaccine ≥14 days before exposure). Genomic sequencing of specimens from 133 patients identified the B.1.617.2 (Delta) variant of SARS-CoV-2, the virus that causes COVID-19, in 119 (89%) and the Delta AY.3 sublineage in one (1%). Overall, 274 (79%) vaccinated patients with breakthrough infection were symptomatic. Among five COVID-19 patients who were hospitalized, four were fully vaccinated; no deaths were reported. Real-time reverse transcription-polymerase chain reaction (RT-PCR) cycle threshold (Ct) values in specimens from 127 vaccinated persons with breakthrough cases were similar to those from 84 persons who were unvaccinated, not fully vaccinated, or whose vaccination status was unknown (median = 22.77 and 21.54, respectively). The Delta variant of SARS-CoV-2 is highly transmissible (1); vaccination is the most important strategy to prevent severe illness and death. On July 27, CDC recommended that all persons, including those who are fully vaccinated, should wear masks in indoor public settings in areas where COVID-19 transmission is high or substantial.* Findings from this investigation suggest that even jurisdictions without substantial or high COVID-19 transmission might consider expanding prevention strategies, including masking in indoor public settings regardless of vaccination status, given the potential risk of infection during attendance at large public gatherings that include travelers from many areas with differing levels of transmission.

Genomic investigation of atypical malaria cases in Kanel, northern Senegal.

BACKGROUND: The diagnosis of malaria cases in regions where the malaria burden has decreased significantly and prevalence is very low is more challenging, in part because of reduced clinical presumption of malaria. The appearance of a cluster of malaria cases with atypical symptoms in Mbounguiel, a village in northern Senegal where malaria transmission is low, in September 2018 exemplifies this scenario. The collaboration between the National Malaria Control Programme (NMCP) at the Senegal Ministry of Health and the Laboratory of Parasitology and Mycology at Cheikh Anta Diop University worked together to evaluate this cluster of malaria cases using molecular and serological tools. METHODS: Malaria cases were diagnosed primarily by rapid diagnostic test (RDT), and confirmed by photo-induced electron transfer-polymerase chain reaction (PET-PCR). 24 single nucleotide polymorphisms (SNPs) barcoding was used for Plasmodium falciparum genotyping. Unbiased metagenomic sequencing and Luminex-based multi-pathogen antibody and antigen profiling were used to assess exposure to other pathogens. RESULTS: Nine patients, of 15 suspected cases, were evaluated, and all nine samples were found to be positive for P. falciparum only. The 24 SNPs molecular barcode showed the predominance of polygenomic infections, with identifiable strains being different from one another. All patients tested positive for the P. falciparum antigens. No other pathogenic infection was detected by either the serological panel or metagenomic sequencing. CONCLUSIONS: This work, undertaken locally within Senegal as a collaboration between the NMCP and a research laboratory at University of Cheikh Anta Diop (UCAD) revealed that a cluster of malaria cases were caused by different strains of P. falciparum. The public health response in real time demonstrates the value of local molecular and genomics capacity in affected countries for disease control and elimination.

miR-222 isoforms are differentially regulated by type-I interferon.

Endogenous microRNAs (miRNAs) often exist as multiple isoforms (known as "isomiRs") with predominant variation around their 3'-end. Increasing evidence suggests that different isomiRs of the same family can have diverse functional roles, as recently demonstrated with the example of miR-222-3p 3'-end variants. While isomiR levels from a same miRNA family can vary between tissues and cell types, change of templated isomiR stoichiometry to stimulation has not been reported to date. Relying on small RNA-sequencing analyses, we demonstrate here that miR-222-3p 3'-end variants >23 nt are specifically decreased upon interferon (IFN) ß stimulation of human fibroblasts, while shorter isoforms are spared. This length-dependent dynamic regulation of long miR-222-3p 3'-isoforms and >40 other miRNA families was confirmed in human monocyte-derived dendritic cells following infection with Salmonella Typhimurium, underlining the breadth of 3'-length regulation by infection, beyond the example of miR-222-3p. We further show that stem-loop miRNA Taqman RT-qPCR exhibits selectivity between 3'-isoforms, according to their length, and that this can lead to misinterpretation of results when these isoforms are differentially regulated. Collectively, and to our knowledge, this work constitutes the first demonstration that the stoichiometry of highly abundant templated 3'-isoforms of a same miRNA family can be dynamically regulated by a stimulus. Given that such 3'-isomiRs can have different functions, our study underlines the need to consider isomiRs when investigating miRNA-based regulation.

Molecular epidemiology of Plasmodium falciparum by multiplexed amplicon deep sequencing in Senegal.

BACKGROUND: Molecular epidemiology can provide important information regarding the genetic diversity and transmission of Plasmodium falciparum, which can assist in designing and monitoring elimination efforts. However, malaria molecular epidemiology including understanding the genetic diversity of the parasite and performing molecular surveillance of transmission has been poorly documented in Senegal. Next Generation Sequencing (NGS) offers a practical, fast and high-throughput approach to understand malaria population genetics. This study aims to unravel the population structure of P. falciparum and to estimate the allelic diversity, multiplicity of infection (MOI), and evolutionary patterns of the malaria parasite using the NGS platform. METHODS: Multiplex amplicon deep sequencing of merozoite surface protein 1 (PfMSP1) and merozoite surface protein 2 (PfMSP2) in fifty-three P. falciparum isolates from two epidemiologically different areas in the South and North of Senegal, was carried out. RESULTS: A total of 76 Pfmsp1 and 116 Pfmsp2 clones were identified and 135 different alleles were found, 56 and 79 belonged to the pfmsp1 and pfmsp2 genes, respectively. K1 and IC3D7 allelic families were most predominant in both sites. The local haplotype diversity (Hd) and nucleotide diversity (π) were higher in the South than in the North for both genes. For pfmsp1, a high positive Tajima's D (TD) value was observed in the South (D = 2.0453) while negative TD value was recorded in the North (D = - 1.46045) and F-Statistic (Fst) was 0.19505. For pfmsp2, non-directional selection was found with a highly positive TD test in both areas and Fst was 0.02111. The mean MOI for both genes was 3.07 and 1.76 for the South and the North, respectively, with a statistically significant difference between areas (p = 0.001). CONCLUSION: This study revealed a high genetic diversity of pfmsp1 and pfmsp2 genes and low genetic differentiation in P. falciparum population in Senegal. The MOI means were significantly different between the Southern and Northern areas. Findings also showed that multiplexed amplicon deep sequencing is a useful technique to investigate genetic diversity and molecular epidemiology of P. falciparum infections.

Amplicon deep sequencing of kelch13 in Plasmodium falciparum isolates from Senegal.

BACKGROUND: In 2006, the Senegalese National Malaria Control Programme recommended artemisinin-based combination therapy (ACT) with artemether-lumefantrine as the first-line treatment for uncomplicated Plasmodium falciparum malaria. To date, multiple mutations associated with artemisinin delayed parasite clearance have been described in Southeast Asia in the Pfk13 gene, such as Y493H, R539T, I543T and C580Y. Even though ACT remains clinically and parasitologically efficacious in Senegal, the spread of resistance is possible as shown by the earlier emergence of resistance to chloroquine in Southeast Asia that subsequently spread to Africa. Therefore, surveillance of artemisinin resistance in malaria endemic regions is crucial and requires the implementation of sensitive tools, such as next-generation sequencing (NGS) which can detect novel mutations at low frequency. METHODS: Here, an amplicon sequencing approach was used to identify mutations in the Pfk13 gene in eighty-one P. falciparum isolates collected from three different regions of Senegal. RESULTS: In total, 10 SNPs around the propeller domain were identified; one synonymous SNP and nine non-synonymous SNPs, and two insertions. Three of these SNPs (T478T, A578S and V637I) were located in the propeller domain. A578S, is the most frequent mutation observed in Africa, but has not previously been reported in Senegal. A previous study has suggested that A578S could disrupt the function of the Pfk13 propeller region. CONCLUSION: As the genetic basis of possible artemisinin resistance may be distinct in Africa and Southeast Asia, further studies are necessary to assess the new SNPs reported in this study.

Widespread Shortening of 3' Untranslated Regions and Increased Exon Inclusion Are Evolutionarily Conserved Features of Innate Immune Responses to Infection.

The contribution of pre-mRNA processing mechanisms to the regulation of immune responses remains poorly studied despite emerging examples of their role as regulators of immune defenses. We sought to investigate the role of mRNA processing in the cellular responses of human macrophages to live bacterial infections. Here, we used mRNA sequencing to quantify gene expression and isoform abundances in primary macrophages from 60 individuals, before and after infection with Listeria monocytogenes and Salmonella typhimurium. In response to both bacteria we identified thousands of genes that significantly change isoform usage in response to infection, characterized by an overall increase in isoform diversity after infection. In response to both bacteria, we found global shifts towards (i) the inclusion of cassette exons and (ii) shorter 3' UTRs, with near-universal shifts towards usage of more upstream polyadenylation sites. Using complementary data collected in non-human primates, we show that these features are evolutionarily conserved among primates. Following infection, we identify candidate RNA processing factors whose expression is associated with individual-specific variation in isoform abundance. Finally, by profiling microRNA levels, we show that 3' UTRs with reduced abundance after infection are significantly enriched for target sites for particular miRNAs. These results suggest that the pervasive usage of shorter 3' UTRs is a mechanism for particular genes to evade repression by immune-activated miRNAs. Collectively, our results suggest that dynamic changes in RNA processing may play key roles in the regulation of innate immune responses.

Negative Selection by Spiral Inertial Microfluidics Improves Viral Recovery and Sequencing from Blood.

In blood samples from patients with viral infection, it is often important to separate viral particles from human cells, for example, to minimize background in performing viral whole genome sequencing. Here, we present a microfluidic device that uses spiral inertial microfluidics with continuous circulation to separate host cells from viral particles and free nucleic acid. We demonstrate that this device effectively reduces white blood cells, red blood cells, and platelets from both whole blood and plasma samples with excellent recovery of viral nucleic acid. Furthermore, microfluidic separation leads to greater viral genome coverage and depth, highlighting an important application of this device in processing clinical samples for viral genome sequencing.

bacterial infection drives the expression dynamics of microRNAs and their isomiRs.

The optimal coordination of the transcriptional response of host cells to infection is essential for establishing appropriate immunological outcomes. In this context, the role of microRNAs (miRNAs)--important epigenetic regulators of gene expression--in regulating mammalian immune systems is increasingly well recognised. However, the expression dynamics of miRNAs, and that of their isoforms, in response to infection remains largely unexplored. Here, we characterized the genome-wide miRNA transcriptional responses of human dendritic cells, over time, to various mycobacteria differing in their virulence as well as to other bacteria outside the genus Mycobacterium, using small RNA-sequencing. We detected the presence of a core temporal response to infection, shared across bacteria, comprising 49 miRNAs, highlighting a set of miRNAs that may play an essential role in the regulation of basic cellular responses to stress. Despite such broadly shared expression dynamics, we identified specific elements of variation in the miRNA response to infection across bacteria, including a virulence-dependent induction of the miR-132/212 family in response to mycobacterial infections. We also found that infection has a strong impact on both the relative abundance of the miRNA hairpin arms and the expression dynamics of miRNA isoforms. That we observed broadly consistent changes in relative arm expression and isomiR distribution across bacteria suggests that this additional, internal layer of variability in miRNA responses represents an additional source of subtle miRNA-mediated regulation upon infection. Collectively, this study increases our understanding of the dynamism and role of miRNAs in response to bacterial infection, revealing novel features of their internal variability and identifying candidate miRNAs that may contribute to differences in the pathogenicity of mycobacterial infections.

Evidence of Ebola Virus Replication and High Concentration in Semen of a Patient During Recovery.

In one patient over time, we found that concentration of Ebola virus RNA in semen during recovery is remarkably higher than blood at peak illness. Virus in semen is replication-competent with no change in viral genome over time. Presence of sense RNA suggests replication in cells present in semen.

A genomic portrait of the genetic architecture and regulatory impact of microRNA expression in response to infection.

MicroRNAs (miRNAs) are critical regulators of gene expression, and their role in a wide variety of biological processes, including host antimicrobial defense, is increasingly well described. Consistent with their diverse functional effects, miRNA expression is highly context dependent and shows marked changes upon cellular activation. However, the genetic control of miRNA expression in response to external stimuli and the impact of such perturbations on miRNA-mediated regulatory networks at the population level remain to be determined. Here we assessed changes in miRNA expression upon Mycobacterium tuberculosis infection and mapped expression quantitative trait loci (eQTL) in dendritic cells from a panel of healthy individuals. Genome-wide expression profiling revealed that ∼40% of miRNAs are differentially expressed upon infection. We find that the expression of 3% of miRNAs is controlled by proximate genetic factors, which are enriched in a promoter-specific histone modification associated with active transcription. Notably, we identify two infection-specific response eQTLs, for miR-326 and miR-1260, providing an initial assessment of the impact of genotype-environment interactions on miRNA molecular phenotypes. Furthermore, we show that infection coincides with a marked remodeling of the genome-wide relationships between miRNA and mRNA expression levels. This observation, supplemented by experimental data using the model of miR-29a, sheds light on the role of a set of miRNAs in cellular responses to infection. Collectively, this study increases our understanding of the genetic architecture of miRNA expression in response to infection, and highlights the wide-reaching impact of altering miRNA expression on the transcriptional landscape of a cell.