Intense P.1 (Gamma) diversification followed by rapid Delta substitution in Southern Brazil: a SARS-CoV-2 genomic epidemiology study.

The analyses of genetic traits, dispersion patterns and phylogenomics are essential for understanding the evolutionary forces driving SARS-CoV-2 viruses in these three years of COVID-19 pandemics. Brazil is one of the most affected countries in the world and not sufficient genomic studies have been performed. The emergence of P.1 lineage led to one of the most serious public health crises on record. Our study presents the genomic sequencing and characterization of 412 samples from Rio Grande do Sul state, in the Brazilian Southern region, during Gamma and Delta epidemic waves, in 2021. Additionally, molecular evolution tests were performed to identify positively selected sites in Brazil between 2020 and 2022, as well as offer some evolutionary perspective about the maintenance of multiple spike mutations in Omicron lineages. Genomic epidemiology analysis has indicated an intense P.1 (Gamma) diversification followed by rapid Delta substitution in Southern Brazil.

Genomic Surveillance of SARS-CoV-2 in Healthcare Workers: A Critical Sentinel Group for Monitoring the SARS-CoV-2 Variant Shift.

SARS-CoV-2 genome surveillance is important for monitoring risk groups and health workers as well as data on new cases and mortality rate due to COVID-19. We characterized the circulation of SARS-CoV-2 variants from May 2021 to April 2022 in the state of Santa Catarina, southern Brazil, and evaluated the similarity between variants present in the population and healthcare workers (HCW). A total of 5291 sequenced genomes demonstrated the circulation of 55 strains and four variants of concern (Alpha, Delta, Gamma and Omicron-sublineages BA.1 and BA.2). The number of cases was relatively low in May 2021, but the number of deaths was higher with the Gamma variant. There was a significant increase in both numbers between December 2021 and February 2022, peaking in mid-January 2022, when the Omicron variant dominated. After May 2021, two distinct variant groups (Delta and Omicron) were observed, equally distributed among the five Santa Catarina mesoregions. Moreover, from November 2021 to February 2022, similar variant profiles between HCW and the general population were observed, and a quicker shift from Delta to Omicron in HCW than in the general population. This demonstrates the importance of HCW as a sentinel group for monitoring disease trends in the general population.

Genomic characterization and molecular evolution of SARS-CoV-2 in Rio Grande do Sul State, Brazil.

SARS-CoV-2 is the virus responsible for the COVID-19 and has afflicted the world since the end of 2019. Different lineages have been discovered and the Gamma lineage, which started the second wave of infections, was first described in Brazil, one of the most affected countries by pandemic. Therefore, this study analyzed SARS-CoV-2 sequenced genomes from Esteio city in Rio Grande do Sul, Southern Brazil. We also comparatively analyzed genomes of the two first years of the pandemic from Rio Grande do Sul state for understanding their genomic and evolutionary patterns. The phylogenomic analysis showed monophyletic groups for Alpha, Gamma, Delta and Omicron, as well as for other circulating lineages in the state. Molecular evolutionary analysis identified several sites under adaptive selection in membrane and nucleocapsid proteins which could be related to a prevalent stabilizing effect on membrane protein structure, as well as majoritarily destabilizing effects on C-terminal nucleocapsid domain.

Detection of Intestinal Dysbiosis in Post-COVID-19 Patients One to Eight Months after Acute Disease Resolution.

The intestinal microbiota plays an important role in the immune response against viral infections, modulating both innate and adaptive immune responses. The cytokine storm is associated with COVID-19 severity, and the patient's immune status is influenced by the intestinal microbiota in a gut-lung bidirectional interaction. In this study, we evaluate the intestinal microbiota of Brazilian patients in different post-COVID-19 periods, and correlate this with clinical data and the antibiotic therapy used during the acute phase. DNA extracted from stool samples was sequenced and total anti-SARS-CoV-2 antibodies and C-reactive protein were quantified. Compared with controls, there were significant differences in the microbiota diversity in post-COVID-19 patients, suggesting an intestinal dysbiosis even several months after acute disease resolution. Additionally, we detected some genera possibly associated with the post-COVID-19 dysbiosis, including Desulfovibrio, Haemophillus, Dialister, and Prevotella, in addition to decreased beneficial microbes, associated with antibiotic-induced dysbiosis, such as Bifidobacterium and Akkermansia. Therefore, our hypothesis is that dysbiosis and the indiscriminate use of antibiotics during the pandemic may be associated with post-COVID-19 clinical manifestations. In our study, 39% (n = 58) of patients reported symptoms, including fatigue, dyspnea, myalgia, alopecia, anxiety, memory loss, and depression. These data suggest that microbiota modulation may represent a target for recovery from acute COVID-19 and a therapeutic approach for post-COVID-19 sequelae.

Emergence of Two Distinct SARS-CoV-2 Gamma Variants and the Rapid Spread of P.1-like-II SARS-CoV-2 during the Second Wave of COVID-19 in Santa Catarina, Southern Brazil.

The western mesoregion of the state of Santa Catarina (SC), Southern Brazil, was heavily affected as a whole by the COVID-19 pandemic in early 2021. This study aimed to evaluate the dynamics of the SARS-CoV-2 virus spreading patterns in the SC state from March 2020 to April 2021 using genomic surveillance. During this period, there were 23 distinct variants, including Beta and Gamma, among which the Gamma and related lineages were predominant in the second pandemic wave within SC. A regionalization of P.1-like-II in the Western SC region was observed, concomitant to the increase in cases, mortality, and the case fatality rate (CFR) index. This is the first evidence of the regionalization of the SARS-CoV-2 transmission in SC and it highlights the importance of tracking the variants, dispersion, and impact of SARS-CoV-2 on the public health systems.

Healthcare-Associated Infections-Related Bacteriome and Antimicrobial Resistance Profiling: Assessing Contamination Hotspots in a Developing Country Public Hospital.

Hospital-built environment colonization by healthcare-associated infections-related bacteria (HAIrB) and the interaction with their occupants have been studied to support more effective tools for HAI control. To investigate HAIrB dynamics and antimicrobial resistance (AMR) profile we carried out a 6-month surveillance program in a developing country public hospital, targeting patients, hospital environment, and healthcare workers, using culture-dependent and culture-independent 16S rRNA gene sequencing methods. The bacterial abundance in both approaches shows that the HAIrB group has important representativeness, with the taxa Enterobacteriaceae, Pseudomonas, Staphylococcus, E. coli, and A. baumannii widely dispersed and abundant over the time at the five different hospital units included in the survey. We observed a high abundance of HAIrB in the patient rectum, hands, and nasal sites. In the healthcare workers, the HAIrB distribution was similar for the hands, protective clothing, and mobile phones. In the hospital environment, the healthcare workers resting areas, bathrooms, and bed equipment presented a wide distribution of HAIrB and AMR, being classified as contamination hotspots. AMR is highest in patients, followed by the environment and healthcare workers. The most frequently detected beta-lactamases genes were, bla SHV-like, bla OXA- 23 -like, bla OXA- 51 -like, bla KPC-like, bla CTX-M- 1, bla CTX-M- 8, and bla CTX-M- 9 groups. Our results demonstrate that there is a wide spread of antimicrobial resistance due to HAIrB in the hospital environment, circulating among patients and healthcare workers. The contamination hotspots identified proved to be constant over time. In the fight for patient safety, these findings can reorient practices and help to set up new guidelines for HAI control.

Equivolumetric Protocol Generates Library Sizes Proportional to Total Microbial Load in 16S Amplicon Sequencing.

High-throughput sequencing of 16S rRNA amplicon has been extensively employed to perform microbiome characterization worldwide. As a culture-independent methodology, it has allowed high-level profiling of sample bacterial composition directly from samples. However, most studies are limited to information regarding relative bacterial abundances (sample proportions), ignoring scenarios in which sample microbe biomass can vary widely. Here, we use an equivolumetric protocol for 16S rRNA amplicon library preparation capable of generating Illumina sequencing data responsive to input DNA, recovering proportionality between observed read counts and absolute bacterial abundances within each sample. Under specified conditions, we show that the estimation of colony-forming units (CFU), the most common unit of bacterial abundance in classical microbiology, is challenged mostly by resolution and taxon-to-taxon variation. We propose Bayesian cumulative probability models to address such issues. Our results indicate that predictive errors vary consistently below one order of magnitude for total microbial load and abundance of observed bacteria. We also demonstrate our approach has the potential to generalize to previously unseen bacteria, but predictive performance is hampered by specific taxa of uncommon profile. Finally, it remains clear that high-throughput sequencing data are not inherently restricted to sample proportions only, and such technologies bear the potential to meet the working scales of traditional microbiology.

The presence of SARS-CoV-2 RNA in human sewage in Santa Catarina, Brazil, November 2019.

Human sewage from Florianopolis (Santa Catarina, Brazil) was analyzed for severe acute respiratory syndrome coronavirus-2 (SARS-CoV2) from October 2019 until March 2020. Twenty five ml of sewage samples were clarified and viruses concentrated using a glycine buffer method coupled with polyethylene glycol precipitation, and viral RNA extracted using a commercial kit. SARS-CoV-2 RNA was detected by RT-qPCR using oligonucleotides targeting N1, S and two RdRp regions. The results of all positive samples were further confirmed by a different RT-qPCR system in an independent laboratory. S and RdRp amplicons were sequenced to confirm identity with SARS-CoV-2. Genome sequencing was performed using two strategies; a sequence-independent single-primer amplification (SISPA) approach, and by direct metagenomics using Illumina's NGS. SARS-CoV-2 RNA was detected on 27th November 2019 (5.49 ± 0.02 log10 SARS-CoV-2 genome copies (GC) L-1), detection being confirmed by an independent laboratory and genome sequencing analysis. The samples in the subsequent three events were positive by all RT-qPCR assays; these positive results were also confirmed by an independent laboratory. The average load was 5.83 ± 0.12 log10 SARS-CoV-2 GC L-1, ranging from 5.49 ± 0.02 log10 GC L-1 (27th November 2019) to 6.68 ± 0.02 log10 GC L-1 (4th March 2020). Our findings demonstrate that SARS-CoV-2 was likely circulating undetected in the community in Brazil since November 2019, earlier than the first reported case in the Americas (21st January 2020).

Swab pooling: A new method for large-scale RT-qPCR screening of SARS-CoV-2 avoiding sample dilution.

To minimize sample dilution effect on SARS-CoV-2 pool testing, we assessed analytical and diagnostic performance of a new methodology, namely swab pooling. In this method, swabs are pooled at the time of collection, as opposed to pooling of equal volumes from individually collected samples. Paired analysis of pooled and individual samples from 613 patients revealed 94 positive individuals. Having individual testing as reference, no false-positives or false-negatives were observed for swab pooling. In additional 18,922 patients screened with swab pooling (1,344 pools), mean Cq differences between individual and pool samples ranged from 0.1 (Cr.I. -0.98 to 1.17) to 2.09 (Cr.I. 1.24 to 2.94). Overall, 19,535 asymptomatic patients were screened using 4,400 RT-qPCR assays. This corresponds to an increase of 4.4 times in laboratory capacity and a reduction of 77% in required tests. Therefore, swab pooling represents a major alternative for reliable and large-scale screening of SARS-CoV-2 in low prevalence populations.

Unusual RNA plant virus integration in the soybean genome leads to the production of small RNAs.

Horizontal gene transfer (HGT) is known to be a major force in genome evolution. The acquisition of genes from viruses by eukaryotic genomes is a well-studied example of HGT, including rare cases of non-retroviral RNA virus integration. The present study describes the integration of cucumber mosaic virus RNA-1 into soybean genome. After an initial metatranscriptomic analysis of small RNAs derived from soybean, the de novo assembly resulted a 3029-nt contig homologous to RNA-1. The integration of this sequence in the soybean genome was confirmed by DNA deep sequencing. The locus where the integration occurred harbors the full RNA-1 sequence followed by the partial sequence of an endogenous mRNA and another sequence of RNA-1 as an inverted repeat and allowing the formation of a hairpin structure. This region recombined into a retrotransposon located inside an exon of a soybean gene. The nucleotide similarity of the integrated sequence compared to other Cucumber mosaic virus sequences indicates that the integration event occurred recently. We described a rare event of non-retroviral RNA virus integration in soybean that leads to the production of a double-stranded RNA in a similar fashion to virus resistance RNAi plants.

ASR5 is involved in the regulation of miRNA expression in rice.

KEY MESSAGE: The work describes an ASR knockdown transcriptomic analysis by deep sequencing of rice root seedlings and the transactivation of ASR cis-acting elements in the upstream region of a MIR gene. MicroRNAs are key regulators of gene expression that guide post-transcriptional control of plant development and responses to environmental stresses. ASR (ABA, Stress and Ripening) proteins are plant-specific transcription factors with key roles in different biological processes. In rice, ASR proteins have been suggested to participate in the regulation of stress response genes. This work describes the transcriptomic analysis by deep sequencing two libraries, comparing miRNA abundance from the roots of transgenic ASR5 knockdown rice seedlings with that of the roots of wild-type non-transformed rice seedlings. Members of 59 miRNA families were detected, and 276 mature miRNAs were identified. Our analysis detected 112 miRNAs that were differentially expressed between the two libraries. A predicted inverse correlation between miR167abc and its target gene (LOC_Os07g29820) was confirmed using RT-qPCR. Protoplast transactivation assays showed that ASR5 is able to recognize binding sites upstream of the MIR167a gene and drive its expression in vivo. Together, our data establish a comparative study of miRNAome profiles and is the first study to suggest the involvement of ASR proteins in miRNA gene regulation.

The source of the river as a nursery for microbial diversity.

Bacteria are highly diverse and ubiquitous organisms that play a key role as drivers for ecosystem processes. The application of NGS (next-generation sequencing technologies) for 16S analysis has been broadly used for understanding bacterioplankton composition and structure. Most of studies conducted on aquatic ecosystems with 16S NGS have been in seawater and lakes. A few studies using NGS have been conducted in river environments and have suggested the presence of a bacterial seed-bank. We performed 16S highly variable V4 region high-throughput analysis in the Sinos River, which is located in one of most important Brazilian industrial centers. This region has several contrasts in its environmental characteristics, presenting a longitudinal gradient of eutrophication and making it a remarkable study site for observing the dynamics of bacterioplankton. We demonstrated consistent evidence for the existence of a bacterial seed-bank and its longitudinal persistence. Seasonal shifts reinforce the importance of the source of the river in maintaining the bacterial seed-bank that spreads throughout the river. Therefore, the preservation of the source of the river is important not only for hydrologic reasons but also to maintain the microbial composition and the ecological integrity of the river.

The Wall-associated Kinase gene family in rice genomes.

The environment is a dynamic system in which life forms adapt. Wall-Associated Kinases (WAK) are a subfamily of receptor-like kinases associated with the cell wall. These genes have been suggested as sensors of the extracellular environment and triggers of intracellular signals. They belong to the ePK superfamily with or without a conserved arginine before the catalytic subdomain VIB, which characterizes RD and non-RD WAKs. WAK is a large subfamily in rice. We performed an extensive comparison of WAK genes from A. thaliana (AtWAK), O. sativa japonica and indica subspecies (OsWAK). Phylogenetic studies and WAK domain characterization allowed for the identification of two distinct groups of WAK genes in Arabidopsis and rice. One group corresponds to a cluster containing only OsWAKs that most likely expanded after the monocot-dicot separation, which evolved into a non-RD kinase class. The other group comprises classical RD-kinases with both AtWAK and OsWAK representatives. Clusterization analysis using extracellular and kinase domains demonstrated putative functional redundancy for some genes, but also highlighted genes that could recognize similar extracellular stimuli and activate different cascades. The gene expression pattern of WAKs in response to cold suggests differences in the regulation of the OsWAK genes in the indica and japonica subspecies. Our results also confirm the hypothesis of functional diversification between A. thaliana and O. sativa WAK genes. Furthermore, we propose that plant WAKs constitute two evolutionarily related but independent subfamilies: WAK-RD and WAK-nonRD. Recognition of this structural division will further provide insights to understanding WAK functions and regulations.

New insights into aluminum tolerance in rice: the ASR5 protein binds the STAR1 promoter and other aluminum-responsive genes.

Aluminum (Al) toxicity in plants is one of the primary constraints in crop production. Al³âº, the most toxic form of Al, is released into soil under acidic conditions and causes extensive damage to plants, especially in the roots. In rice, Al tolerance requires the ASR5 gene, but the molecular function of ASR5 has remained unknown. Here, we perform genome-wide analyses to identify ASR5-dependent Al-responsive genes in rice. Based on ASR5_RNAi silencing in plants, a global transcriptome analysis identified a total of 961 genes that were responsive to Al treatment in wild-type rice roots. Of these genes, 909 did not respond to Al in the ASR5_RNAi plants, indicating a central role for ASR5 in Al-responsive gene expression. Under normal conditions, without Al treatment, the ASR5_RNAi plants expressed 1.756 genes differentially compared to the wild-type plants, and 446 of these genes responded to Al treatment in the wild-type plants. Chromatin immunoprecipitation followed by deep sequencing identified 104 putative target genes that were directly regulated by ASR5 binding to their promoters, including the STAR1 gene, which encodes an ABC transporter required for Al tolerance. Motif analysis of the binding peak sequences revealed the binding motif for ASR5, which was confirmed via in vitro DNA-binding assays using the STAR1 promoter. These results demonstrate that ASR5 acts as a key transcription factor that is essential for Al-responsive gene expression and Al tolerance in rice.

isomiRID: a framework to identify microRNA isoforms.

SUMMARY: MicroRNAs (miRNAs) have been extensively studied owing to their important regulatory roles in genic expression. An increasingly number of reports are performing extensive data mining in small RNA sequencing libraries to detect miRNAs isoforms and also 5' and 3' post-transcriptional nucleotide additions, as well as edited miRNAs sequences. A ready to use pipeline, isomiRID, was developed to standardize and automatize the search for miRNAs isoforms in high-throughput small RNA sequencing libraries. AVAILABILITY: isomiRID is a command line Python script available at http://www.ufrgs.br/RNAi/isomiRID/.

Heavy metal-associated isoprenylated plant protein (HIPP): characterization of a family of proteins exclusive to plants.

Metallochaperones are key proteins for the safe transport of metallic ions inside the cell. HIPPs (heavy metal-associated isoprenylated plant proteins) are metallochaperones that contain a metal binding domain (HMA) and a C-terminal isoprenylation motif. In this study, we provide evidence that proteins of this family are found only in vascular plants and may be separated into five distinct clusters. HIPPs may be involved in (a) heavy metal homeostasis and detoxification mechanisms, especially those involved in cadmium tolerance, (b) transcriptional responses to cold and drought, and (c) plant-pathogen interactions. In particular, our results show that the rice (Oryza sativa) HIPP OsHIPP41 gene is highly expressed in response to cold and drought stresses, and its product is localized in the cytosol and the nucleus. The results suggest that HIPPs play an important role in the development of vascular plants and in plant responses to environmental changes.

Identifying conserved and novel microRNAs in developing seeds of Brassica napus using deep sequencing.

MicroRNAs (miRNAs) are important post-transcriptional regulators of plant development and seed formation. In Brassica napus, an important edible oil crop, valuable lipids are synthesized and stored in specific seed tissues during embryogenesis. The miRNA transcriptome of B. napus is currently poorly characterized, especially at different seed developmental stages. This work aims to describe the miRNAome of developing seeds of B. napus by identifying plant-conserved and novel miRNAs and comparing miRNA abundance in mature versus developing seeds. Members of 59 miRNA families were detected through a computational analysis of a large number of reads obtained from deep sequencing two small RNA and two RNA-seq libraries of (i) pooled immature developing stages and (ii) mature B. napus seeds. Among these miRNA families, 17 families are currently known to exist in B. napus; additionally 29 families not reported in B. napus but conserved in other plant species were identified by alignment with known plant mature miRNAs. Assembled mRNA-seq contigs allowed for a search of putative new precursors and led to the identification of 13 novel miRNA families. Analysis of miRNA population between libraries reveals that several miRNAs and isomiRNAs have different abundance in developing stages compared to mature seeds. The predicted miRNA target genes encode a broad range of proteins related to seed development and energy storage. This work presents a comparative study of the miRNA transcriptome of mature and developing B. napus seeds and provides a basis for future research on individual miRNAs and their functions in embryogenesis, seed maturation and lipid accumulation in B. napus.

Metatranscriptomic analysis of small RNAs present in soybean deep sequencing libraries.

A large number of small RNAs unrelated to the soybean genome were identified after deep sequencing of soybean small RNA libraries. A metatranscriptomic analysis was carried out to identify the origin of these sequences. Comparative analyses of small interference RNAs (siRNAs) present in samples collected in open areas corresponding to soybean field plantations and samples from soybean cultivated in greenhouses under a controlled environment were made. Different pathogenic, symbiotic and free-living organisms were identified from samples of both growth systems. They included viruses, bacteria and different groups of fungi. This approach can be useful not only to identify potentially unknown pathogens and pests, but also to understand the relations that soybean plants establish with microorganisms that may affect, directly or indirectly, plant health and crop production.