Mutational profile confers increased stability of SARS-CoV-2 spike protein in Brazilian isolates

Spike (S) protein has been recognized as a promising molecular target for diagnostic, vaccines and antiviral drugs development for COVID-19. In this study, we analyzed the most predominant mutations in the S protein of Brazilian isolates and predicted the effect of these amino acid alterations to protein conformation. A total of 25,924 sequences were obtained from GISAID for five regions of Brazilian territory (Midwest, North, Northeast, South, and Southeast), according to exclusion criteria. Most of the SARS-CoV-2 isolates belongs to the G clade and showed a large occurrence of D614G, N501Y and L18F substitutions. Prediction effects of these amino acid substitutions on the structure dynamics of the spike protein indicated a positive ΔΔG values and negative ΔΔS_Vib in most cases which is associated to structural stabilization and flexibility reduction of the S protein. Mutations E484K, N501Y and K417N belong to several SARS-CoV-2 variants of concern such as Alpha, Beta, Gamma and Delta, and showed high incidence among Brazilian isolates. These mutations have been described to increase RBD affinity to ACE-2 host and abolishment of RBD affinity to potent neutralizing ant-RBD. The increase in rates of infection and reinfection requires continuous genomic surveillance studies in order to characterize emerging mutations and monitor vaccine efficacy, and thus consideration structural data and dynamics in the observed phenotypes.Communicated by Ramaswamy H. Sarma.

Parallel Evolution of Pseudomonas aeruginosa during a Prolonged ICU-Infection Outbreak

Most knowledge about Pseudomonas aeruginosa pathoadaptation is derived from studies on airway colonization in cystic fibrosis; little is known about adaptation in acute settings. P. aeruginosa frequently affects burned patients and the burn wound niche has distinct properties that likely influence pathoadaptation. This study aimed to genetically and phenotypically characterize P. aeruginosa isolates collected during an outbreak of infection in a burn intensive care unit (ICU). Sequencing reads from 58 isolates of ST1076 P. aeruginosa taken from 23 patients were independently mapped to a complete reference genome for the lineage (H25338); genetic differences were identified and were used to define the population structure. Comparative genomic analysis at single-nucleotide resolution identified pathoadaptive genes that evolved multiple, independent mutations. Three key phenotypic assays (growth performance, motility, carbapenem resistance) were performed to complement the genetic analysis for 47 unique isolates. Population structure for the ST1076 lineage revealed 11 evolutionary sublineages. Fifteen pathoadaptive genes evolved mutations in at least two sublineages. The most prominent functional classes affected were transcription/two-component regulatory systems, and chemotaxis/motility and attachment. The most frequently mutated gene was oprD, which codes for outer membrane porin involved in uptake of carbapenems. Reduced growth performance and motility were found to be adaptive phenotypic traits, as was high level of carbapenem resistance, which correlated with higher carbapenem consumption during the outbreak. Multiple prominent linages evolved each of the three traits in parallel providing evidence that they afford a fitness advantage for P. aeruginosa in the context of human burn infection. IMPORTANCE Pseudomonas aeruginosa is a Gram-negative pathogen causing infections in acutely burned patients. The precise mechanisms required for the establishment of infection in the burn setting, and adaptive traits underpinning prolonged outbreaks are not known. We have assessed genotypic data from 58 independent P. aeruginosa isolates taken from a single lineage that was responsible for an outbreak of infection in a burn ICU that lasted for almost 2.5 years and affected 23 patients. We identified a core set of 15 genes that we predict to control pathoadaptive traits in the burn infection based on the frequency with which independent mutations evolved. We combined the genotypic data with phenotypic data (growth performance, motility, antibiotic resistance) and clinical data (antibiotic consumption) to identify adaptive phenotypes that emerged in parallel. High-level carbapenem resistance evolved rapidly, and frequently, in response to high clinical demand for this antibiotic class during the outbreak.

Adaptive specialization of a unique sponge body from the Cambrian Qingjiang biota

Sponge fossils from the Cambrian black shales have attracted attention from both palaeontologists and geochemists for many years in terms of their high diversity, beautiful preservation and perplexing adaptation to inhospitable living environments. However, the body shape of these sponges, which contributes to deciphering adaptive evolution, has not been scrutinized. New complete specimens of the hexactinellid sponge Sanshapentella tentoriformis sp. nov. from the Qingjiang biota (black shale of the Cambrian Stage 3 Shuijingtuo Formation, ca 518 Ma) allow recognition of a unique dendriform body characterized by a columnar trunk with multiple conical high peaks and distinctive quadripod-shaped dermal spicules that frame each high peak. The body shape of this new sponge along with other early Cambrian hexactinellids, is classified into three morpho-groups that reflect different levels of adaptivity to the environment. The cylindrical and ovoid bodies generally adapted to a large spectrum of environments; however, the dendriform body of S. tentoriformis was restricted to the relatively deep-water, oxygen-deficient environment. From a hindsight view, the unique body shape represents a consequence of adaptation that helps maintain an effective use of oxygen and a low energy cost in hypoxic conditions.

Cell Type-Specific Differentiation Between Indica and Japonica Rice Root Tip Responses to Different Environments Based on Single-Cell RNA Sequencing

Background: As Oryza sativa ssp. indica and Oryza sativa ssp. japonica are the two major subspecies of Asian cultivated rice, the adaptative evolution of these varieties in divergent environments is an important topic in both theoretical and practical studies. However, the cell type-specific differentiation between indica and japonica rice varieties in response to divergent habitat environments, which facilitates an understanding of the genetic basis underlying differentiation and environmental adaptation between rice subspecies at the cellular level, is little known.

Methods: We analyzed a published single-cell RNA sequencing dataset to explore the differentially expressed genes between indica and japonica rice varieties in each cell type. To estimate the relationship between cell type-specific differentiation and environmental adaptation, we focused on genes in the WRKY, NAC, and BZIP transcription factor families, which are closely related to abiotic stress responses. In addition, we integrated five bulk RNA sequencing datasets obtained under conditions of abiotic stress, including cold, drought and salinity, in this study. Furthermore, we analyzed quiescent center cells in rice root tips based on orthologous markers in Arabidopsis.

Results: We found differentially expressed genes between indica and japonica rice varieties with cell type-specific patterns, which were enriched in the pathways related to root development and stress reposes. Some of these genes were members of the WRKY, NAC, and BZIP transcription factor families and were differentially expressed under cold, drought or salinity stress. In addition, LOC_Os01g16810, LOC_Os01g18670, LOC_Os04g52960, and LOC_Os08g09350 may be potential markers of quiescent center cells in rice root tips.

Conclusion: These results identified cell type-specific differentially expressed genes between indica-japonica rice varieties that were related to various environmental stresses and provided putative markers of quiescent center cells. This study provides new clues for understanding the development and physiology of plants during the process of adaptative divergence, in addition to identifying potential target genes for the improvement of stress tolerance in rice breeding applications.

Stable and Efficient Biosynthesis of 5-Aminolevulinic Acid Using Plasmid-Free Escherichia coli

5-Aminolevulinic acid (5-ALA) is a key metabolic intermediate of the heme biosynthesis pathway, which has broad application prospects in agriculture and medicine. However, segregational instability of plasmid-based expression systems and low yield have hampered large-scale manufacture of 5-ALA. In this study, two important genes of the 5-ALA C5 biosynthesis pathway, hemA and hemL, were integrated into Escherichia coli MG1655 for chemically induced chromosomal evolution (CIChE). The highest hemA and hemL copy-number, 98 per genome, was obtained in CIChE strain MG136. The 5-ALA titer of this strain reached 2724 mg/L in optimized condition. Then, after undergoing adaptative evolution and the deletion of recA, strain MG136a ΔrecA::FRT could stably produce 4550 mg/L 5-ALA from glucose, 450 times the amount produced by hemA- hemL single copy strain MG1655-hemAL. This study constructed a plasmid-free E. coli strain for 5-ALA production, which will provide the basis for further manipulation of metabolic regulation and optimization of fermentation.