The recent Zika virus (ZIKV) outbreak has linked ZIKV with microcephaly and other central nervous system pathologies in humans. Astrocytes are among the first cells to respond to ZIKV infection in the brain and are also targets for virus infection. In this study, we investigated the interaction between ZIKV and primary human brain cortical astrocytes (HBCA). HBCAs were highly sensitive to representatives of both Asian and African ZIKV lineages and produced high viral yields. The infection was associated with limited immune cytokine/chemokine response activation; the highest increase of expression, following infection, was seen in CXCL-10 (IP-10), interleukin-6, 8, 12, and CCL5 (RANTES). Ultrastructural changes in the ZIKV-infected HBCA were characterized by electron tomography (ET). ET reconstructions elucidated high-resolution 3D images of the proliferating and extensively rearranged endoplasmic reticulum (ER) containing viral particles and virus-induced vesicles, tightly juxtaposed to collapsed ER cisternae. The results confirm that human astrocytes are sensitive to ZIKV infection and could be a source of proinflammatory cytokines in the ZIKV-infected brain tissue.

We followed the presence of Zika virus (ZIKV) in four healthy adults (two men and two women), for periods ranging from 78 to 298 days post symptom onset. The patients were evaluated regarding the presence of the virus in different body fluids (blood, saliva, urine and semen), development of immune responses (including antibodies, cytokines and chemokines), and virus genetic variation within samples collected from semen and urine during the infection course. The analysis was focused primarily on the two male patients who shed the virus for up to 158 days after the initial symptoms. ZIKV particles were detected in the spermatozoa cytoplasm and flagella, in immature sperm cells and could also be isolated from semen in cell culture, confirming that the virus is able to preserve integrity and infectivity during replication in the male reproductive system (MRS). Despite the damage caused by ZIKV infection within the MRS, our data showed that ZIKV infection did not result in infertility at least in one of the male patients. This patient was able to conceive a child after the infection. We also detected alterations in the male genital cytokine milieu, which could play an important role in the replication and transmission of the virus which could considerably increase the risk of ZIKV sexual spread. In addition, full genome ZIKV sequences were obtained from several samples (mainly semen), which allowed us to monitor the evolution of the virus within a patient during the infection course. We observed genetic changes over time in consensus sequences and lower frequency intra-host single nucleotide variants (iSNV), that suggested independent compartmentalization of ZIKV populations in the reproductive and urinary systems. Altogether, the present observations confirm the risks associated with the long-term replication and shedding of ZIKV in the MRS and help to elucidate patterns of intra-host genetic evolution during long term replication of the virus.

The causative agent of COVID-19 pandemic, SARS-CoV-2, has a 29,903 bases positive-sense single-stranded RNA genome. RNAs exhibit about 150 modified bases that are essential for proper function. Among internal modified bases, the N⁶-methyladenosine, or m6A, is the most frequent, and is implicated in SARS-CoV-2 immune response evasion. Although the SARS-CoV-2 genome is RNA, almost all genomes sequenced thus far are, in fact, reverse transcribed complementary DNAs. This process reduces the true complexity of these viral genomes because the incorporation of dNTPs hides RNA base modifications. Here, we present an initial exploration of Nanopore direct RNA sequencing to assess the m6A residues in the SARS-CoV-2 sequences of ORF3a, E, M, ORF6, ORF7a, ORF7b, ORF8, N, ORF10 and the 3′-untranslated region. We identified fifteen m6A methylated positions, of which, six are in ORF N. Additionally, because m6A is associated with the DRACH motif, we compared its distribution in major SARS-CoV-2 variants. Although DRACH is highly conserved among variants, we show that variants Beta and Eta have a fourth position C > U change in DRACH at 28,884b that could affect methylation. This is the first report of direct RNA sequencing of a Brazilian SARS-CoV-2 sample coupled with the identification of modified bases.

The current COVID-19 pandemic has emphasized the vulnerability of communities living in the urban outskirts and informal settlements. The lack of reliable COVID-19 case data highlights the importance and application of wastewater-based epidemiology. This study aimed to monitor the COVID-19 trends in four vulnerable urban communities (slums and low-income neighborhoods) in metropolitan São Paulo by assessing the SARS-CoV-2 RNA viral load in wastewater. We analyzed 160 samples from May 2020 to June 2021 with weekly or fortnightly samplings. The samples were ultracentrifuged with glycine elution and quantified by N1/N2 SARS-CoV-2 RT-qPCR. The results of positivity were 100% (Paraisópolis, Heliópolis and Cidade Tiradentes) and 76.9% (Vila Brasilândia). The new case numbers of COVID-19, counted from the onset of symptoms, positively correlated with SARS-CoV-2 N1 viral loads from the two largest communities (p<0.001). SARS-CoV-2 infectivity was tested in Vero E6 cells after concentration with the two techniques, ultrafiltration (Centricon® Plus-70 10 kDa) and sucrose cushion ultracentrifugation, but none of the evaluated samples presented positive results. Next-generation sequencing (NGS) analysis from samples collected in March and August 2021 revealed the presence of the clade 20 J (lineage P.1) belonging to the most prevalent circulating variant in the country. Our results showed that wastewater surveillance data can be used as complementary indicators to monitor the dynamics and temporal trends of COVID-19. The infectivity test results strengthened the evidence of low risk of infection associated with SARS-CoV-2 in wastewater.

Nearly two decades after the last epidemic caused by a severe acute respiratory syndrome coronavirus (SARS-CoV), newly emerged SARS-CoV-2 quickly spread in 2020 and precipitated an ongoing global public health crisis. Both the continuous accumulation of point mutations, owed to the naturally imposed genomic plasticity of SARS-CoV-2 evolutionary processes, as well as viral spread over time, allow this RNA virus to gain new genetic identities, spawn novel variants and enhance its potential for immune evasion. Here, through an in-depth phylogenetic clustering analysis of upwards of 200,000 whole-genome sequences, we reveal the presence of previously unreported and hitherto unidentified mutations and recombination breakpoints in Variants of Concern (VOC) and Variants of Interest (VOI) from Brazil, India (Beta, Eta and Kappa) and the USA (Beta, Eta and Lambda). Additionally, we identify sites with shared mutations under directional evolution in the SARS-CoV-2 Spike-encoding protein of VOC and VOI, tracing a heretofore-undescribed correlation with viral spread in South America, India and the USA. Our evidence-based analysis provides well-supported evidence of similar pathways of evolution for such mutations in all SARS-CoV-2 variants and sub-lineages. This raises two pivotal points: (i) the co-circulation of variants and sub-lineages in close evolutionary environments, which sheds light onto their trajectories into convergent and directional evolution, and (ii) a linear perspective into the prospective vaccine efficacy against different SARS-CoV-2 strains.

The causative agent of COVID-19 pandemic, the SARS-CoV-2 coronavirus, has a 29,903 bases positive-sense single-stranded RNA genome. RNAs exhibit about 100 modified bases that are essential for proper function. Among internal modified bases, the N ⁶ -methyladenosine, or m6A, is the most frequent, and is implicated in SARS-CoV-2 immune response evasion. Although the SARS-CoV-2 genome is RNA, almost all genomes sequenced so far are in fact, reverse transcribed complementary DNAs. This process reduces the true complexity of these viral genomes because incorporation of dNTPs hides RNA base modifications. Here, in this perspective paper, we present an initial exploration of the Nanopore direct RNA sequencing to assess the m6A residues in the SARS-CoV-2 sequences of ORF3a, E, M, ORF6, ORF7a, ORF7b, ORF8, N, ORF10 and the 3’-untranslated region. We identified 15 m6A methylated positions, of which, 6 are in ORF N. Also, because m6A is associated with the DRACH motif, we compared its distribution in major SARS-CoV-2 variants. Although DRACH is highly conserved among variants we show that variants Beta and Eta have a fourth position C>U change in DRACH at 28,884b that could affect methylation. The Nanopore technology offers a unique opportunity for the study of viral epitranscriptomics. This technique is PCR-free and is not sequencing-by-synthesis, therefore, no PCR bias and synthesis errors are introduced. The modified bases are preserved and assessed directly with no need for chemical treatments or antibodies. This is the first report of direct RNA sequencing of a Brazilian SARS-CoV-2 sample coupled with the identification of modified bases.