Phylogeny-aware linear B-cell epitope predictor detects targets associated with immune response to orthopoxviruses.

We introduce a phylogeny-aware framework for predicting linear B-cell epitope (LBCE)-containing regions within proteins. Our approach leverages evolutionary information by using a taxonomic scaffold to build models trained on hierarchically structured data. The resulting models present performance equivalent or superior to generalist methods, despite using simpler features and a fraction of the data volume required by current state-of-the-art predictors. This allows the utilization of available data for major pathogen lineages to facilitate the prediction of LBCEs for emerging infectious agents. We demonstrate the efficacy of our approach by predicting new LBCEs in the monkeypox (MPXV) and vaccinia viruses. Experimental validation of selected targets using sera from infected patients confirms the presence of LBCEs, including candidates for the differential serodiagnosis of recent MPXV infections. These results point to the use of phylogeny-aware predictors as a useful strategy to facilitate the targeted development of immunodiagnostic tools.

Donor-derived acute myeloid leukemia in solid organ transplantation.

We report the transmission of acute myeloid leukemia (AML) undetected at donation from a deceased organ donor to two kidneys and one liver recipients. We reviewed the medical records, and performed molecular analyses and whole exome sequencing (WES) to ascertain AML donor origin and its molecular evolution. The liver recipient was diagnosed 11 months after transplantation and died from complications 2 months later. The two kidney recipients (R1 and R2) were diagnosed 19 and 20 months after transplantation and both received treatment for leukemia. R1 died of complications 11 months after diagnosis, while R2 went into complete remission for 44 months, before relapsing. R2 died 10 months later of complications from allogenic bone marrow transplantation. Microsatellite analysis demonstrated donor chimerism in circulating cells from both kidney recipients. Targeted molecular analyses and medical records revealed NPM1 mutation present in the donor and recipients, while FLT3 was mutated only in R1. These findings were confirmed by WES, which revealed additional founder and clonal mutations, and HLA genomic loss in R2. In conclusion, we report the first in-depth genomic analysis of AML transmission following solid organ transplantation, revealing distinct clonal evolution, and providing a potential molecular explanation for tumor escape.

Coevolved Positions Represent Key Functional Properties in the Trypsin-Like Serine Proteases Protein Family.

Trypsin-like serine proteases are a group of homologous enzymes which exert multiple roles in both vertebrate and invertebrate organisms. Key properties of these enzymes include their activation from an inactive zymogen form to their active form by cleavage of residues in their N-terminus, the presence of a conserved catalytic triad of residues, and the existence of different patterns of substrate selectivity for residue cleavage between the various members of this protein family. In this article, we apply the decomposition of residue coevolution networks computational method to find sets of residues related to some of these key properties, especially to zymogen activation. Positive selection detection, normal modes analysis, and the calculation of thermal couplings between the bovine trypsinogen and bovine trypsin structures residues yielded further information for understanding the zymogen activation process and highlighted the importance of some of the coevolved set residues during these transitions.

Molecular characterization and gene expression analysis of the pro-inflammatory cytokines IL-1ß and IL-8 in the South American fish Piaractus mesopotamicus challenged with Aeromonas dhakensis.

In the present study, the complete characterization of cDNA and genomic sequences of IL-1ß and IL-8, as well as the expression profile of these genes in the South American fish pacu (Piaractus mesopotamicus) is provided. The full-length pmIL-1ß cDNA was composed of 1208 nucleotides that would produce a precursor peptide with 273 amino acid residues. A putative caspase-1 cleavage site, similar to what is found in mammalian IL-1ß, was identified producing a mature peptide with a theoretical molecular weight of 17.21 kDa. The pmIL-8 cDNA sequence consisted of 1019 nucleotides which encoded a 95-amino acid protein with a theoretical molecular weight of 10.43 kDa that showed all typical CXC chemokine features, including a 20-residue signal peptide and four conserved cysteine residues. Constitutive mRNA expression was detected for both genes in the liver, head kidney, gill, intestine, skin and spleen. After a bacterial challenge, up-regulation was detected for both pmIL-1ß and pmIL-8 in the spleen and head kidney at 12 h post-infection. At 24 h post-infection there was a decrease in the expression of both genes, with pmIL-8 showing a significant down-regulation in the liver and head kidney when compared to the control groups.

Characterization and gene expression analysis of pacu (Piaractus mesopotamicus) inducible nitric oxide synthase (iNOS) following Aeromonas dhakensis infection.

Nitric oxide (NO) is an important effector molecule which is involved in a myriad of biological processes, including immune responses against pathogens such as parasites, virus and bacteria. During the inflammatory processes in vertebrates, NO is produced by the inducible nitric oxide synthase (iNOS) enzyme in practically all nucleated cells to suppress or kill intracellular pathogens. The aim of the present study was to characterize the full coding region of the iNOS gene of pacu (Piaractus mesopotamicus), an economically and ecologically important South American fish species, and to analyze mRNA expression levels following intraperitoneal infection with the pathogenic bacterium Aeromonas dhakensis by means of quantitative real time PCR (qPCR). The results showed that the pacu iNOS transcript is 3237 bp in length, encoding a putative protein composed of 1078 amino acid residues. The amino acid sequence showed similarities ranging from 69.03% to 94.34% with other teleost fish and 57.70% with the human iNOS, with all characteristic domains and cofactor binding sites of the enzyme detected. Phylogenetic analysis showed that the iNOS from the red-bellied piranha, another South American characiform, was the closest related sequence to the pacu iNOS. iNOS transcripts were constitutively detected in the liver, spleen and head kidney, and there was a significant upregulation in the liver and spleen at 12, 24 and 48 h after infection with A. dhakensis. No significant variations were observed in the head kidney during the periods analyzed. These results show that iNOS expression was induced by A. dhakensis infection and suggest that this enzyme may be involved in the response to this bacterium in pacu.

Sequence-independent characterization of viruses based on the pattern of viral small RNAs produced by the host.

Virus surveillance in vector insects is potentially of great benefit to public health. Large-scale sequencing of small and long RNAs has previously been used to detect viruses, but without any formal comparison of different strategies. Furthermore, the identification of viral sequences largely depends on similarity searches against reference databases. Here, we developed a sequence-independent strategy based on virus-derived small RNAs produced by the host response, such as the RNA interference pathway. In insects, we compared sequences of small and long RNAs, demonstrating that viral sequences are enriched in the small RNA fraction. We also noted that the small RNA size profile is a unique signature for each virus and can be used to identify novel viral sequences without known relatives in reference databases. Using this strategy, we characterized six novel viruses in the viromes of laboratory fruit flies and wild populations of two insect vectors: mosquitoes and sandflies. We also show that the small RNA profile could be used to infer viral tropism for ovaries among other aspects of virus biology. Additionally, our results suggest that virus detection utilizing small RNAs can also be applied to vertebrates, although not as efficiently as to plants and insects.

Genome-Wide Survey of Genes Under Positive Selection in Avian Pathogenic Escherichia coli Strains.

The ability to obtain bacterial genomes from the same host has allowed for comparative studies that help in the understanding of the molecular evolution of specific pathotypes. Avian pathogenic Escherichia coli (APEC) is a group of extraintestinal strains responsible for causing colibacillosis in birds. APEC is also suggested to possess a role as a zoonotic agent. Despite its importance, APEC pathogenesis still has several cryptic pathogenic processes that need to be better understood. In this work, a genome-wide survey of eight APEC strains for genes with evidence of recombination revealed that ∼14% of the homologous groups evaluated present signs of recombination. Enrichment analyses revealed that nine Gene Ontology (GO) terms were significantly more represented in recombinant genes. Among these GO terms, several were noted to be ATP-related categories. The search for positive selection in these APEC genomes revealed 32 groups of homologous genes with evidence of positive selection. Among these groups, we found several related to cell metabolism, as well as several uncharacterized genes, beyond the well-known virulence factors ompC, lamB, waaW, waaL, and fliC. A GO term enrichment test showed a prevalence of terms related to bacterial cell contact with the external environment (e.g., viral entry into host cell, detection of virus, pore complex, bacterial-type flagellum filament C, and porin activity). Finally, the genes with evidence of positive selection were retrieved from genomes of non-APEC strains and tested as were done for APEC strains. The result revealed that none of the groups of genes presented evidence of positive selection, confirming that the analysis was effective in inferring positive selection for APEC and not for E. coli in general, which means that the study of the genes with evidence of positive selection identified in this study can contribute for the better understanding of APEC pathogenesis processes.

Daytime soybean transcriptome fluctuations during water deficit stress.

BACKGROUND: Since drought can seriously affect plant growth and development and little is known about how the oscillations of gene expression during the drought stress-acclimation response in soybean is affected, we applied Illumina technology to sequence 36 cDNA libraries synthesized from control and drought-stressed soybean plants to verify the dynamic changes in gene expression during a 24-h time course. Cycling variables were measured from the expression data to determine the putative circadian rhythm regulation of gene expression. RESULTS: We identified 4866 genes differentially expressed in soybean plants in response to water deficit. Of these genes, 3715 were differentially expressed during the light period, from which approximately 9.55% were observed in both light and darkness. We found 887 genes that were either up- or down-regulated in different periods of the day. Of 54,175 predicted soybean genes, 35.52% exhibited expression oscillations in a 24 h period. This number increased to 39.23% when plants were submitted to water deficit. Major differences in gene expression were observed in the control plants from late day (ZT16) until predawn (ZT20) periods, indicating that gene expression oscillates during the course of 24 h in normal development. Under water deficit, dissimilarity increased in all time-periods, indicating that the applied stress influenced gene expression. Such differences in plants under stress were primarily observed in ZT0 (early morning) to ZT8 (late day) and also from ZT4 to ZT12. Stress-related pathways were triggered in response to water deficit primarily during midday, when more genes were up-regulated compared to early morning. Additionally, genes known to be involved in secondary metabolism and hormone signaling were also expressed in the dark period. CONCLUSIONS: Gene expression networks can be dynamically shaped to acclimate plant metabolism under environmental stressful conditions. We have identified putative cycling genes that are expressed in soybean leaves under normal developmental conditions and genes whose expression oscillates under conditions of water deficit. These results suggest that time of day, as well as light and temperature oscillations that occur considerably affect the regulation of water deficit stress response in soybean plants.

Divergent immune microenvironments in two tumor nodules from a patient with mismatch repair-deficient prostate cancer.

Patients with prostate cancer (PC) generally do not respond favorably to immune checkpoint inhibitors, which may be due to a low abundance of tumor-infiltrating lymphocytes even when mutational load is high. Here, we identified a patient who presented with high-grade primary prostate cancer with two adjacent tumor nodules. While both nodules were mismatch repair-deficient (MMRd), exhibited pathogenic MSH2 and MSH6 alterations, had a high tumor mutational burden (TMB), and demonstrated high microsatellite instability (MSI), they had markedly distinct immune phenotypes. The first displayed a dense infiltrate of lymphocytes ("hot nodule"), while the second displayed significantly fewer infiltrating lymphocytes ("cold nodule"). Whole-exome DNA analysis found that both nodules shared many identical mutations, indicating that they were derived from a single clone. However, the cold nodule appeared to be sub-clonal relative to the hot nodule, suggesting divergent evolution of the cold nodule from the hot nodule. Whole-transcriptome RNA analysis found that the cold nodule demonstrated lower expression of genes related to antigen presentation (HLA) and, paradoxically, classical tumor immune tolerance markers such as PD-L1 (CD274) and CTLA-4. Immune cell deconvolution suggested that the hot nodule was enriched not only in CD8+ and CD4 + T lymphocytes, but also in M1 macrophages, activated NK cells, and γδ T cells compared to the cold nodule. This case highlights that MMRd/TMB-high PC can evolve to minimize an anti-tumor immune response, and nominates downregulation of antigen presentation machinery (HLA loss) as a potential mechanism of adaptive immune evasion in PC.

A directed approach for the identification of transcripts harbouring the spliced leader sequence and the effect of trans-splicing knockdown in Schistosoma mansoni.

Schistosomiasis is a major neglected tropical disease caused by trematodes from the genus Schistosoma. Because schistosomes exhibit a complex life cycle and numerous mechanisms for regulating gene expression, it is believed that spliced leader (SL) trans-splicing could play an important role in the biology of these parasites. The purpose of this study was to investigate the function of trans-splicing in Schistosoma mansoni through analysis of genes that may be regulated by this mechanism and via silencing SL-containing transcripts through RNA interference. Here, we report our analysis of SL transcript-enriched cDNA libraries from different S. mansoni life stages. Our results show that the trans-splicing mechanism is apparently not associated with specific genes, subcellular localisations or life stages. In cross-species comparisons, even though the sets of genes that are subject to SL trans-splicing regulation appear to differ between organisms, several commonly shared orthologues were observed. Knockdown of trans-spliced transcripts in sporocysts resulted in a systemic reduction of the expression levels of all tested trans-spliced transcripts; however, the only phenotypic effect observed was diminished larval size. Further studies involving the findings from this work will provide new insights into the role of trans-splicing in the biology of S. mansoni and other organisms. All Expressed Sequence Tags generated in this study were submitted to dbEST as five different libraries. The accessions for each library and for the individual sequences are as follows: (i) adult worms of mixed sexes (LIBEST_027999: JZ139310 - JZ139779), (ii) female adult worms (LIBEST_028000: JZ139780 - JZ140379), (iii) male adult worms (LIBEST_028001: JZ140380 - JZ141002), (iv) eggs (LIBEST_028002: JZ141003 - JZ141497) and (v) schistosomula (LIBEST_028003: JZ141498 - JZ141974).

Why are some plants taller? Researchers on the unveiling of genetic variation associated with complex quantitative phenotypes.

Francisco Pereira Lobo, Giovanni Marques de Castro, and Felipe Campelo are part of an international team of collaborators that developed CALANGO, a comparative genomics tool to investigate quantitative genotype-phenotype relationships. Their Patterns article highlights how the tool integrates species-centric data to perform genome-wide search and detect genes potentially involved in the emergence of complex quantitative traits across species. Here, they talk about their view of data science, their experience with interdisciplinary research, and the potential applications of their tool.

Chronic ethanol exposure impairs alveolar leukocyte infiltration during pneumococcal pneumonia, leading to an increased bacterial burden despite increased CXCL1 and nitric oxide levels.

Ethanol abuse is a risk factor for the development of pneumonia caused by Streptococcus pneumoniae, a critical pathogen for public health. The aim of this article was to investigate the inflammatory mechanisms involved in pneumococcal pneumonia that may be associated with chronic ethanol exposure. Male C57BL6/J-Unib mice were exposed to 20% (v/v) ethanol for twelve weeks and intranasally infected with 5x104 CFU of S. pneumoniae. Twenty-four hours after infection, lungs, bronchoalveolar lavage and blood samples were obtained to assess the consequences of chronic ethanol exposure during infection. Alcohol-fed mice showed increased production of nitric oxide and CXCL1 in alveoli and plasma during pneumococcal pneumonia. Beside this, ethanol-treated mice exhibited a decrease in leukocyte infiltration into the alveoli and reduced frequency of severe lung inflammation, which was associated with an increase in bacterial load. Curiously, no changes were observed in survival after infection. Taken together, these results demonstrate that chronic ethanol exposure alters the inflammatory response during S. pneumoniae lung infection in mice with a reduction in the inflammatory infiltrate even in the presence of higher levels of the chemoattractant CXCL1.

Mpox vaccine and infection-driven human immune signatures.

Background: Mpox (formerly known as monkeypox) outbreaks outside endemic areas peaked in July 2022, infecting > 85,000 people and raising concerns about our preparedness against this emerging viral pathogen. Licensed and approved for mpox, the JYNNEOS vaccine has fewer side effects than previous smallpox vaccines and demonstrated efficacy against mpox infection in humans. Comparing JYNNEOS vaccine- and mpox-induced immunity is imperative to evaluate JYNNEOS' immunogenicity and inform vaccine administration and design. Methods: We examined the polyclonal serum (ELISA) and single B cell (heavy chain gene and transcriptome data) antibody repertoires and T cells (AIM and ICS assays) induced by the JYNNEOS vaccine as well as mpox infection. Findings: Gene-level plasmablast and antibody responses were negligible and JYNNEOS vaccinee sera displayed minimal binding to recombinant mpox proteins and native proteins from the 2022 outbreak strain. In contrast, recent mpox infection (within 20-102 days) induced robust serum antibody responses to A29L, A35R, A33R, B18R, and A30L, and to native mpox proteins, compared to vaccinees. JYNNEOS vaccine recipients presented comparable CD4 and CD8 T cell responses against orthopox peptides to those observed after mpox infection. Interpretation: JYNNEOS immunization does not elicit a robust B cell response, and its immunogenicity may be mediated by T cells. Funding: Research reported in this publication was supported, in part, by the National Cancer Institute of the National Institutes of Health under Award Number U54CA267776, U19AI168631(VS), as well as institutional funds from the Icahn School of Medicine.

CALANGO: A phylogeny-aware comparative genomics tool for discovering quantitative genotype-phenotype associations across species.

Living species vary significantly in phenotype and genomic content. Sophisticated statistical methods linking genes with phenotypes within a species have led to breakthroughs in complex genetic diseases and genetic breeding. Despite the abundance of genomic and phenotypic data available for thousands of species, finding genotype-phenotype associations across species is challenging due to the non-independence of species data resulting from common ancestry. To address this, we present CALANGO (comparative analysis with annotation-based genomic components), a phylogeny-aware comparative genomics tool to find homologous regions and biological roles associated with quantitative phenotypes across species. In two case studies, CALANGO identified both known and previously unidentified genotype-phenotype associations. The first study revealed unknown aspects of the ecological interaction between Escherichia coli, its integrated bacteriophages, and the pathogenicity phenotype. The second identified an association between maximum height in angiosperms and the expansion of a reproductive mechanism that prevents inbreeding and increases genetic diversity, with implications for conservation biology and agriculture.

Mpox vaccine and infection-driven human immune signatures: an immunological analysis of an observational study.

BACKGROUND: Monkeypox virus has recently infected more than 88 000 people, raising concerns about our preparedness against this emerging viral pathogen. Licensed and approved for mpox, the JYNNEOS vaccine has fewer side-effects than previous smallpox vaccines and has shown immunogenicity against monkeypox in animal models. This study aims to elucidate human immune responses to JYNNEOS vaccination compared with mpox-induced immunity. METHODS: Peripheral blood mononuclear cells and sera were obtained from ten individuals vaccinated with one or two doses of JYNNEOS and six individuals diagnosed with monkeypox virus infection. Samples were obtained from seven individuals before vaccination to serve as a baseline. We examined the polyclonal serum (ELISA) and single B-cell (heavy chain gene and transcriptome data) antibody repertoires and T-cell responses (activation-induced marker and intracellular cytokine staining assays) induced by the JYNNEOS vaccine versus monkeypox virus infection. FINDINGS: All participants were men between the ages of 21 and 60 years, except for one woman in the group of mpox-convalescent individuals, and none had previous orthopoxvirus exposure. All mpox cases were mild. Vaccinee samples were collected 6-33 days after the first dose and 5-40 days after the second dose. Mpox-convalescent samples were collected 20-102 days after infection. In vaccine recipients, gene-level plasmablast and antibody responses were negligible and sera displayed moderate binding to recombinant orthopoxviral proteins (A29L, A35R, E8L, A30L, A27L, A33R, B18R, and L1R) and native proteins from the 2022 monkeypox outbreak strain. By contrast, recent monkeypox virus infection (within 20-102 days) induced robust serum antibody responses to monkeypox virus proteins and to native monkeypox virus proteins from a viral isolate obtained during the 2022 outbreak. JYNNEOS vaccine recipients presented robust orthopoxviral CD4+ and CD8+ T-cell responses. INTERPRETATION: Infection with monkeypox virus resulted in robust B-cell and T-cell responses, whereas immunisation with JYNNEOS elicited more robust T-cell responses. These data can help to inform vaccine design and policies for preventing mpox in humans. FUNDING: National Cancer Institute (National Institutes of Health), National Institute of Allergy and Infectious Diseases (National Institutes of Health), and Icahn School of Medicine.

Spatiotemporal dynamics of the resistome and virulome of riverine microbiomes disturbed by a mining mud tsunami.

Aquatic ecosystems are highly vulnerable to anthropogenic activities. However, it remains unclear how the microbiome responds to press disturbance events in these ecosystems. We examined the impact of the world's largest mining disaster (Brazil, 2015) on sediment microbiomes in two disturbed rivers compared to an undisturbed river during 390 days post-disturbance. The diversity and structure of the virulome and microbiome, and of antibiotic and metal resistomes, consistently differed between the disturbed and undisturbed rivers, particularly at day 7 post-disturbance. 684 different ARGs were predicted, 38% were exclusive to the disturbed rivers. Critical antibiotic resistance genes (ARGs), e.g., mcr and ereA2, were significantly more common in the disturbed microbiomes. 401 different ARGs were associated with mobile genetic elements (MGEs), 95% occurred in the disturbed rivers. While plasmids were the most common MGEs with a broad spectrum of ARGs, spanning 16 antibiotic classes, integrative conjugative elements (ICEs) and integrons disseminated ARGs associated with aminoglycoside and tetracycline, and aminoglycoside and beta-lactam, respectively. A significant increase in the relative abundance of class 1 integrons, ICEs, and pathogens was identified at day 7 in the disturbed microbiomes, 72-, 14- and 3- fold higher, respectively, compared with the undisturbed river. Mobile ARGs associated with ESKAPEE group pathogens, while metal resistance genes and virulence factor genes in nonpathogenic hosts predominated in all microbiomes. Network analysis showed highly interconnected ARGs in the disturbed communities, including genes targeting antibiotics of last resort. Interactions between copper and beta-lactam/aminoglycoside/macrolide resistance genes, mostly mobile and critical, were also uncovered. We conclude that the mud tsunami resulted in resistome expansion, enrichment of pathogens, and increases in promiscuous and mobile ARGs. From a One Health perspective, mining companies need to move toward more environmentally friendly and socially responsible mining practices to reduce risks associated with pathogens and critical and mobile ARGs.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families.

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

Exploring antibiotic resistance in environmental integron-cassettes through intI-attC amplicons deep sequencing.

INTRODUCTION: Freshwater ecosystems provide propitious conditions for the acquisition and spread of antibiotic resistance genes (ARGs), and integrons play an important role in this process. MATERIAL AND METHODS: In the present study, the diversity of putative environmental integron-cassettes, as well as their potential bacterial hosts in the Velhas River (Brazil), was explored through intI-attC and 16S rRNA amplicons deep sequencing. RESULTS AND DISCUSSION: ORFs related to different biological processes were observed, from DNA integration to oxidation-reduction. ARGs-cassettes were mainly associated with class 1 mobile integrons carried by pathogenic Gammaproteobacteria, and possibly sedentary chromosomal integrons hosted by Proteobacteria and Actinobacteria. Two putative novel ARG-cassettes homologs to fosB3 and novA were detected. Regarding 16SrRNA gene analysis, taxonomic and functional profiles unveiled Proteobacteria, Bacteroidetes, Firmicutes, and Actinobacteria as dominant phyla. Betaproteobacteria, Alphaproteobacteria, and Actinobacteria classes were the main contributors for KEGG orthologs associated with resistance. CONCLUSIONS: Overall, these results provide new information about environmental integrons as a source of resistance determinants outside clinical settings and the bacterial community in the Velhas River.

Analysis of a Marseillevirus Transcriptome Reveals Temporal Gene Expression Profile and Host Transcriptional Shift.

Marseilleviruses comprise a family of large double-stranded DNA viruses belonging to the proposed order "Megavirales." These viruses have a circular genome of ∼370 kbp, coding hundreds of genes. Over a half of their genes are associated with AT-rich putative promoter motifs, which have been demonstrated to be important for gene regulation. However, the transcriptional profile of Marseilleviruses is currently unknown. Here we used RNA sequencing technology to get a general transcriptional profile of Marseilleviruses. Eight million 75-bp-long nucleotide sequences were robustly mapped to all 457 genes initially predicted for Marseillevirus isolate T19, the prototype strain of the family, and we were able to assemble 359 viral contigs using a genome-guided approach with stringent parameters. These reads were differentially mapped to the genes according to the replicative cycle time point from which they were obtained. Cluster analysis indicated the existence of three main temporal categories of gene expression, early, intermediate and late, which were validated by quantitative reverse transcription polymerase chain reaction assays targeting several genes. Genes belonging to different functional groups exhibited distinct expression levels throughout the infection cycle. We observed that the previously predicted promoter motif, AAATATTT, as well as new predicted motifs, were not specifically related to any of the temporal or functional classes of genes, suggesting that other components are involved in temporally regulating virus transcription. Moreover, the host transcription machinery is heavily altered, and many genes are down regulated, including those related to translation process. This study provides an overview of the transcriptional landscape of Marseilleviruses.