Global diversity and distribution of prophages are lineage-specific within the Ralstonia solanacearum plant pathogenic bacterium species complex Ralstonia solanacearum is a destructive plant pathogenic bacterium and the causative agent of bacterial wilt disease, infecting over 200 plant species worldwide. In addition to chromosomal genes, its virulence is mediated by mobile genetic elements including integrated DNA of bacteriophages, i.e., prophages, which may carry fitness-associated auxiliary genes or modulate host gene expression. Although experimental studies have characterised several prophages that shape R. solanacearum virulence, the global diversity, distribution, and wider functional gene content of R. solanacearum prophages is unknown. In this study, prophages were identified in a diverse collection of 192 R. solanacearum draft genome assemblies originating from six continents. Prophages were identified bioinformatically and their diversity investigated using genetic distance measures, gene content, GC, and total length. Prophage distribution was characterised using metadata on R. solanacearum geographic origin and lineage classification (phylotypes), and their functional gene content was assessed by identifying putative prophage-encoded auxiliary genes. In total, 343 intact prophages were identified, forming ten genetically distinct clusters. These included five prophage clusters belonging to the Inoviridae, Myoviridae, and Siphoviridae phage families, and five uncharacterised clusters, possibly representing novel, previously undescribed phages. The prophages had broad geographical distribution being present across multiple continents. However, they were generally host phylogenetic lineage-specific, and overall, prophage diversity was proportional to the genetic diversity of their hosts. The prophages contained a myriad of auxiliary genes involved in metabolism and virulence of both phage and bacteria. Our results show that while R. solanacearum prophages are highly diverse globally, they make lineage-specific contributions to the R. solanacearum accessory genome, which could have resulted from shared coevolutionary history.
Loci for prediction of penicillin and tetracycline susceptibility in Neisseria gonorrhoeae: a genome wide association study BackgroundWhile Neisseria gonorrhoeae poses an urgent public health threat because of increasing antimicrobial resistance, much of the circulating population remains susceptible to historical treatment regimens. Point-of-care diagnostics that report susceptibility could allow for reintroduction of these regimens, but development of such diagnostics has been limited to ciprofloxacin, for which susceptibility can be predicted from a single locus.

MethodsWe assembled a dataset of 12,045 N. gonorrhoeae genomes with phenotypic resistance data for tetracycline (n = 3,611) and penicillin (n = 6,935). Using conditional genome wide association studies (GWAS), we sought to define genetic variants associated with susceptibility to penicillin and tetracycline. We evaluated the sensitivity and specificity of these variants for predicting susceptibility and non-resistance in our collection of gonococcal genomes.

FindingsIn our conditional penicillin GWAS, the presence of a genetic variant defined by a non-mosaic penA allele without an insertion at codon 345 was significantly associated with penicillin susceptibility and had the highest negative effect size of significant variants (p = 5.0 x 10-14, {beta} = -2.5). In combination with the absence of blaTEM, this variant predicted penicillin susceptibility with high specificity (99.8%) and modest sensitivity (36.7%). For tetracycline, the wild type allele at rpsJ codon 57, encoding valine, was significantly associated with tetracycline susceptibility (p = 5.6 x 10-16, {beta} = -1.6) after conditioning on the presence of tetM. The combination of rpsJ codon 57 allele and tetM absence predicted tetracycline susceptibility with high specificity (97.2%) and sensitivity (88.7%).

InterpretationAs few as two genetic loci can predict susceptibility to penicillin and tetracycline in N. gonorrhoeae with high specificity. Molecular point-of-care diagnostics targeting these loci have the potential to increase available treatments for gonorrhea.

FundingNational Institute of Allergy and Infectious Diseases, the National Science Foundation, and the Smith Family Foundation

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed with the terms "Neisseria gonorrhoeae" and "diagnostic" or "assay" plus "penicillin" or "tetracycline" for reports in any language published up to July 1, 2021. We additionally searched for "Neisseria gonorrhoeae" and "genome wide association study". We found that previously proposed molecular diagnostics for penicillin and tetracycline susceptibility either exclusively focused on plasmid-mediated resistance (i.e., targeting blaTEM or tetM) or did not include variants in genes encoding antibiotic targets (e.g., did not include penA or rpsJ). Targets for molecular surveillance have focused on resistance-associated alleles rather than susceptibility-associated alleles. We did not find any previous penicillin or tetracycline conditional genome wide association studies (GWAS) in N. gonorrhoeae.

Added value of this studyTo identify targets for molecular diagnostics that predict penicillin and tetracycline susceptibility, we conducted GWAS conditioning on the presence of plasmid-mediated resistance determinants to detect chromosomal loci with the highest association with susceptibility. We discovered a sequence (penA_01) that differentiates susceptible isolates from those with a resistance-associated insertion at codon 345 and from those with mosaic penA alleles, which is associated with penicillin susceptibility. We also found that rpsJ codon 57 was the chromosomal locus contributing the most to tetracycline susceptibility. The combination of these chromosomal loci and the absence of plasmid encoded determinants predicts penicillin and tetracycline susceptibility with high specificity in both a large global collection of N. gonorrhoeae and a validation dataset consisting of recently published genomes from CDCs Gonococcal Isolate Surveillance Program (GISP) surveillance collected in 2018.

Implications of all the available evidenceThe chromosomal loci penA_01 and rpsJ codon 57 in combination with plasmid loci blaTEM and tetM are candidates for the development of point-of-care molecular diagnostics for penicillin and tetracycline susceptibility. The loci may be combined with the currently available ciprofloxacin susceptibility diagnostics to predict susceptibility to multiple antibiotics. Additionally, our study suggests that conditional GWAS focused on variants associated with susceptibility may be a promising approach to identify minimal sets of loci for molecular diagnostics and surveillance.
Spatio-temporal patterns of multi-trophic biodiversity and food-web characteristics uncovered across a river catchment using environmental DNA Accurate characterisation of ecological communities with respect to their biodiversity and food-web structure is essential for conservation. However, combined empirical study of biodiversity and multi-trophic food-webs at a large spatial and temporal resolution has been prohibited by the lack of appropriate access to such data from natural systems. Here, we assessed biodiversity and food-web characteristics across a 700 km2 riverine network over seasons using environmental DNA. We found contrasting biodiversity patterns between major taxonomic groups; local richness showed seasonally dependent and statistically significant increases and decreases towards downstream positions within the catchment for fish and bacteria respectively, while invertebrate richness remained unchanged with increased downstream position. The local food-webs, formed by these taxonomic groups, also showed a variation in their structure, such as link density and nestedness, to both space and time, yet these patterns did not necessarily mirror those of biodiversity and functional feeding characteristics. In order to conserve species diversity as well as their functional trophic integrity of communities, patterns of biodiversity and food-web characteristics must thus be jointly studied, as our results suggest that they are not directly scalable to each other even at the same spatial and temporal scales. Kinetics and Mechanism of Fentanyl Dissociation from the mu-Opioid Receptor Driven by illicit fentanyl, opioid related deaths have reached the highest level in 2020. Currently, an opioid over-dose is resuscitated by the use of naloxone, which competitively binds and antagonizes the {micro}-opioid receptor (mOR). Thus, knowledge of the residence times of opioids at mOR and the unbinding mechanisms is valuable for assessing the effectiveness of naloxone. In the present study we calculate the fentanyl-mOR dissociation time and elucidate the mechanism by applying an enhanced sampling molecular dynamics (MD) technique. Two sets of metadynamics simulations with different initial structures were performed while accounting for the protonation state of the conserved H2976.52, which has been suggested to modulate the ligand-mOR affinity and binding mode. Surprisingly, with the N{delta}-protonated H2976.52, fentanyl can descend as much as 10 [A] below the level of the conserved D1473.32 before escaping the receptor, and has a calculated residence time{tau} of 38 s. In contrast, with the N - and doubly protonated H2976.52, the calculated{tau} are 2.6 s and 0.9 s, respectively. Analysis suggests that formation of the piperidine-Hid297 hydrogen bond strengthens the hydrophobic contacts with the transmembrane helix (TM) 6, allowing fentanyl to explore a deep pocket. Considering the experimental{tau} of [~]4 min for fentanyl and the role of TM6 in mOR activation, we suggest that the deep insertion mechanism may be biologically relevant. The work paves the way for large-scale computational predictions of opioid dissociation rates to inform the regulatory decision. The profound role of the histidine protonation state found here may shift the paradigm in computational studies of ligand-receptor kinetics. An early signaling GARP transcription factor regulates the entry point for differentiation in Giardia Transcriptional regulation of differentiation is critical for parasitic pathogens to adapt to environmental changes and regulate transmission. In response to encystation stimuli, Giardia lamblia cells shift from G1+S phase of the cell cycle to G2+M. By 2-4 hours, cyst wall proteins are upregulated, indicating that key regulatory steps occur within the first four hours of encystation. However, all characterized transcription factors (TFs) in Giardia have only been investigated at later time points of encystation. How TFs initiate encystation and link it to the cell cycle remains enigmatic. Here, we systematically screened six putative early upregulated TFs for nuclear localization, established their dynamic expression profiles, and determined their functional role in regulating encystation. We found a critical repressor, GLP4 that increases rapidly after 30 min of encystation stimuli and downregulates encystation specific markers, including Cyst Wall Proteins (CWPs) and enzymes in the cyst N-acetyl galactosamine (GalNAc) pathway. Depletion of GLP4 increases cyst production. Importantly, we observe that G2+M cells exhibit higher levels of CWP1 resulting from the activation of MYB2, a TF previously linked to encystation in Giardia. GLP4 upregulation occurs in G1+S cells, suggesting a role in repressing MYB2 and encystation specific genes in the G1+S phase of the cell cycle. Furthermore, we demonstrate that depletion of GLP4 upregulates MYB2 and promotes encystation while overexpression of GLP4 downregulates MYB2 and represses encystation. Together, these results suggest that Giardia employs a dose-dependent transcriptional response that involves the cell cycle regulated repressor GLP4 to orchestrate MYB2 and entry into the encystation pathway.

ImportanceTransition between life cycle stages is a common feature among parasitic pathogens and its regulation must be optimized to balance persistence of infection with transmission. The early transcription factors (TFs) regulating commitment to differentiate are totally unknown in Giardia. In this work, we identified GLP4, a previously uncharacterized GARP-like TF, as an early-acting transcriptional repressor that inhibits G1+S cells from entering the encystation pathway. GLP4 is therefore a key regulator controlling the balance between proliferative growth and terminal differentiation into infective cysts. Regulatory perturbations of ribosome allocation reshape the growth proteome with a trade-off in adaptation capacity Bacteria regulate their cellular resource allocation to enable fast growth-adaptation to a variety of environmental niches. We studied the ribosomal allocation, growth and expression profiles of two sets of fast-growing mutants of Escherichia coli K-12 MG1655 in glucose minimal medium. Mutants with only 3 of the seven copies of ribosomal RNA operons grew faster than the wild-type strain in minimal media and show similar phenotype to previously studied rpoB mutants. Higher growth rates due to increased ribosome content affected resource allocation. Expression profiles of fast-growing mutants shared downregulation of hedging functions and upregulated growth functions. Mutants showed longer diauxic shifts and reduced activity of gluconeogenic promoters during glucose-acetate shifts, suggesting reduced availability of the RNA Polymerase for expressing hedging proteome. These results show that the regulation of ribosomal allocation underlies the growth/hedging phenotypes obtained from laboratory evolution experiments. We show how two different regulatory perturbations (rRNA promoters or rpoB mutations) reshape the proteome for growth with a concomitant fitness cost

HighlightsMutants with only 3 ribosomal operons grow faster than wild-type in minimal medium

 {Delta}4 rrn and rpoB mutants share phenotypic traits

Faster growth of mutants is achieved by increased ribosome content

Fast-growing mutants display reduced hedging expression and adaptation trade-offs

Despite similar ribosomal content in rich medium the mutants present growth defects Imbalanced segregation of recombinant haplotypes in hybrid populations reveals inter- and intrachromosomal Dobzhansky-Muller incompatibilities Dobzhansky-Muller incompatibilities (DMIs) are a major component of reproductive isolation between species. DMIs imply negative epistasis, and are exposed when two diverged populations hybridize. Mapping the locations of DMIs has largely relied on classical genetic mapping. Approaches to date are hampered by low power and the challenge of identifying DMI loci on the same chromosome, because strong initial linkage of parental haplotypes weakens statistical tests. Here, we propose new statistics to infer negative epistasis from haplotype frequencies in hybrid populations. When two divergent populations hybridize, the variance in heterozygosity at two loci decreases faster with time at DMI loci than at random pairs of loci. When two populations hybridize at near-even admixture proportions, the deviation of the observed variance from its expectation becomes negative for the DMI pair. This negative deviation enables us to detect intermediate to strong negative epistasis both within and between chromosomes. In practice, the detection window in hybrid populations depends on the demographic scenario, the recombination rate, and the strength of epistasis. When the initial proportion of the two parental populations is uneven, only strong DMIs can be detected with our method unless migration prevents parental haplotypes from being lost. We use the new statistics to infer candidate DMIs from three hybrid populations of swordtail fish. We identify numerous new DMI candidates, some of which are inferred to interact with several loci within and between chromosomes. Moreover, we discuss our results in the context of an expected enrichment in intrachromosomal over interchromosomal DMIs.

Author SummaryGenetic incompatibility in the form of (Bateson-)Dobzhansky-Muller incompatibilities (DMIs) is an important component of reproductive isolation between species. However, the evolutionary role of DMIs during the process of speciation is contentious. DMIs occur when two or more genetic variants interact to reduce their carriers fitness. Once recombination combines incompatible variants in hybrids, selection acts to remove these variants from the population. One step towards addressing the evolutionary role of DMIs is to quantify the prevalence of DMIs in incipient and hybridizing species. Here, we present statistics that are sensitive to the resulting recombinant imbalance and that can indicate the location of DMIs in hybrid genomes in various demographic scenarios. We use simulations to show that the time window during which a DMI is detectable depends on its genomic location, the severity of the DMI, and the populations demography. Importantly, our statistic distinguishes genetic associations arising due to physical linkage from those arising due to gene interactions, which allows for the inference of both inter- and intrachromosomal DMIs. Applying our statistics to three hybrid populations of swordtail fish, we confirm previously known DMIs and identify new candidate incompatibilities. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Force redistribution during clathrin-mediated endocytosis revealed by new phase-separating force sensors Forces are central to countless cellular processes, yet in vivo force measurement at the molecular scale remains difficult if not impossible. During clathrin-mediated endocytosis, forces produced by the actin cytoskeleton are transmitted to the plasma membrane by a multi-protein coat for membrane deformation. However, the magnitudes of these forces remain unknown. Here, we present new in vivo force sensors that phase separate under force. We measured the forces on the fission yeast HIP1R homologue End4p, a protein that links the membrane to the actin cytoskeleton. End4p is under [~]15 pN force near the actin cytoskeleton, [~] 11 pN near the clathrin lattice, and [~]9 pN near the plasma membrane. Our results demonstrate that forces are collected and redistributed across the endocytic machinery.

One-Sentence SummaryPhase-separating coiled-coil force sensors reveal an unexpected mode of force transduction during endocytosis in vivo. Systematic monitoring of 2-Cys peroxiredoxin-derived redox signals unveiled its role in attenuating carbon assimilation rate Transmission of reductive cues from the photosynthetic electron transport chain to redox-regulated proteins plays a crucial role in activating chloroplast metabolism. However, deciphering the role of their counterbalanced oxidative signals is challenging due to monitoring difficulties. Here, we demonstrate the light-depended redox modification of chloroplast-targeted 2-Cys peroxiredoxins and introduce peroxiredoxin-based biosensors to monitor photosynthetically-derived oxidative signals. By employing a set of genetically encoded biosensors, we show the induction of oxidative signals under habitual light intensities and their inverse relationship with NADPH levels, unraveling the combined activity of reducing and oxidizing signals in fine-tuning chloroplast metabolism. A faster increase in carbon assimilation rates during photosynthesis induction phase was measured in plants deficient in 2-Cys peroxiredoxins compared to wild-type, suggesting the involvement of oxidative signals in attenuating photosynthesis under variable light environments. We suggest that oxidative signals measured by peroxiredoxin-based biosensors reflect the limitation to photosynthesis imposed by the redox regulatory system.

One-Sentence SummaryA genetically encoded biosensor unmasked the dominant role of photosynthetically-derived oxidative signals under habitual conditions. Relationship between shear wave velocity and muscle activation is inconsistent across different muscle types There is an increasing use of shear wave ultrasound elastography to quantify mechanical properties of muscles under various conditions such as changes muscle length and levels of activation in healthy and pathological muscle. However, little is known about the variability in shear wave velocity among muscles as most studies investigate one specific muscle. The purpose of this study was to determine if the relationship between SWV and muscle activation is consistent across muscles with different architectures: biceps brachii, tibialis anterior, and medial gastrocnemius, All measures were made at matching levels of activation and approximately at the optimal length for each muscle to control for length-dependent changes in the relationship between activation and force or stiffness. We also conducted a control experiment to determine how the passive force within a muscle alters the relationship between muscle activation and shear wave velocity. The relationship between shear wave velocity-squared and activation above 10% MVC differed across muscles, with biceps brachii and medial gastrocnemius showing a lower slope than tibialis anterior. Shear wave velocity-squared also differed between muscles at the shortest length (p<0.001) and the increase in shear wave velocity-squared with muscle lengthening also differed between muscle types (p = 0.005) Muscle-specific differences could not be explained by the architectural features such as pennation angle, fiber length, and physiological cross-sectional area. Our results demonstrate that there is not a unique relationship between muscle activation and shear wave velocity, highlighting the importance of understanding the many factors contributing to shear wave propagation in muscle before elastography can be used to make quantitative comparisons across muscle types. Testing for phylogenetic signal in single-cell RNA-seq data AO_SCPLOWBSTRACTC_SCPLOWPhylogenetic methods are emerging as a useful tool to understand cancer evolutionary dynamics, including tumor structure, heterogeneity, and progression. Most currently used approaches utilize either bulk whole genome sequencing (WGS) or single-cell DNA sequencing (scDNA-seq) and are based on calling copy number alterations and single nucleotide variants (SNVs). Here we explore the potential of single-cell RNA sequencing (scRNA-seq) to reconstruct cancer evolutionary dynamics. scRNA-seq is commonly applied to explore differential gene expression of cancer cells throughout tumor progression. The method exacerbates the single-cell sequencing problem of low yield per cell with uneven expression levels. This accounts for low and uneven sequencing coverage and makes SNV detection and phylogenetic analysis challenging. In this paper, we demonstrate for the first time that scRNA-seq data contains sufficient evolutionary signal and can be utilized in phylogenetic analyses. We explore and compare results of such analyses based on both expression levels and SNVs called from scRNA-seq data. Both techniques are shown to be useful for reconstructing phylogenetic relationships between cells, reflecting the clonal composition of a tumor. Both standardized expression values and SNVs appear to be equally capable of reconstructing a similar pattern of phylogenetic relationship. This pattern is stable even when phylogenetic uncertainty is taken in account. Our results open up a new direction of somatic phylogenetics based on scRNA-seq data. Further research is required to refine and improve these approaches to capture the full picture of somatic evolutionary dynamics in cancer. Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Darwin's small and medium ground finches might have taste preferences, but not for human foods Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwins finches consume human-introduced foods preferentially; however, it remains unclear why. Here we presented pastry with flavour profiles typical of human foods (oily, salty, sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive selection of human foods. If human-food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site-differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human-food flavours were preferred compared to control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwins finches may have acquired a tolerance to the flavours of human food, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization. Meiotic DSB-independent role of protein phosphatase 4 in Hop1 assembly to promote meiotic chromosome axis formation in budding yeast Dynamic changes in chromosomal structure that occur during meiotic prophase play an important role in the progression of meiosis. Among these structures, meiosis-specific chromosomal axis-loop structures are important as a scaffold for integrated control between the meiotic recombination reaction and the associated checkpoint system to ensure accurate chromosome segregation. However, the molecular mechanism of the initial step of chromosome axis-loop construction is not well understood. Here, we showed that in budding yeast, protein phosphatase 4 (PP4) that primarily counteracts Mec1/Tel1 phosphorylation is required to promote the assembly of chromosomal axis components Hop1 and Red1 onto meiotic chromatin via interaction with Hop1. PP4, in contrast, did not affect Rec8 assembly. Notably, unlike the previously known function of PP4, this novel function of PP4 was independent of meiotic DSB-dependent Tel1/Mec1 kinase activities. The defect in Hop1/Red1 assembly in the absence of PP4 function was not suppressed by dysfunction of Pch2, which removes Hop1 protein from the chromosome axis. This suggests that PP4 is required for the initial step of chromatin loading of Hop1 rather than stabilization of Hop1 on axes. These results indicate that phosphorylation/dephosphorylation-mediated regulation of Hop1 recruitment onto chromatin occurs during chromosome axis construction before meiotic double-strand break formation. GUIdEStaR (G-quadruplex, uORF, IRES, Epigenetics, Small RNA, Repeats), the integrated metadatabase in conjunction with neural network methods GUIdEStaR integrates existing databases of various types of G-quadruplex, upstream Open Reading Frame (uORF), Internal Ribosome Entry Site (IRES), methylation to RNA and histone protein, small RNA, and repeats. GUIdEStaR consists of approx. 40,000 genes and 320,000 transcripts. An mRNA transcript is divided into 5 regions (5UTR, 3UTR, exon, intron, and biological region) where each region contains presence-absence data of 169 different types of elements. Recently, artificial intelligence (AI) based analysis of sequencing data has been gaining popularity in the area of bioinformatics. GUIdEStaR generates datasets that can be used as inputs to AI methods. At the GUIdEStaR homepage, users submit gene symbols by clicking a "Send" button, and shortly result files in CSV format are available for download at the result website. Users have an option to send the result files to their email addresses. Additionally, the entire database and the example Java codes are also freely available for download. Here, we demonstrate the database usage with three neural network classification studies-1) small RNA study for classifying transcription factor (TF) genes into either one of TF mediated by small RNA originated from SARS-CoV-2 or by human microRNA (miRNA), 2) cell membrane receptor study for classifying receptor genes as either with virus interaction or without one, and 3) nonsense mediated mRNA decay (NMD) study for classifying cell membrane and nuclear receptors as either NMD target or non-target. GUIdEStaR is available for access to the easy-to-use web-based database at www.guidestar.kr and for download at https://sourceforge.net/projects/guidestar. Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and High-Frequency mutation N501Y in different levels of control measures Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1 058 020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3 % of global COVID-19 cases and determined three types of mutation dynamics: High-Frequency, Medium-Frequency, and Low-Frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that Medium-Frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of High-Frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with High-Frequency or not-High-Frequency and both follows similar trends in different levels of stringency. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Human 14-3-3 proteins site-selectively bind the mutational hotspot region of SARS-CoV-2 nucleoprotein modulating its phosphoregulation Phosphorylated within its Ser/Arg-rich region, the SARS-CoV-2 nucleoprotein (N) recruits the phosphopeptide-binding human 14-3-3 proteins that play a well-recognized role in replication of many viruses. Here we use genetic code expansion to demonstrate that phosphorylation of SARS-CoV-2 N at either of two pseudo-repeats centered at Ser197 and Thr205 is sufficient for 14-3-3 binding. According to fluorescence anisotropy, the pT205-motif, present in SARS-CoV-2 but not in SARS-CoV, is preferred over the pS197-motif by all seven human 14-3-3 isoforms, which display unforeseen pT205/pS197 binding selectivity hierarchy. Crystal structures demonstrate that pS197 and pT205 are mutually exclusive 14-3-3-binding sites, whereas SAXS and biochemical data indicate 14-3-3 binding occludes the Ser/Arg-rich region, inhibiting its dephosphorylation. This Ser/Arg-rich region of N is highly prone to mutations, as exemplified by the Omicron and Delta variants, with our data suggesting how the strength of its 14-3-3 binding can be linked with the replicative fitness of the virus. HSP90 modulates T2R bitter taste receptor nitric oxide production and innate immune responses in human airway epithelial cells and macrophages Bitter taste receptors (T2Rs) are G protein-coupled receptors (GPCRs) expressed in various cell types including ciliated airway epithelial cells and macrophages. T2Rs in these two airway innate immune cell types are activated by bitter products, including some secreted by Pseudomonas aeruginosa, leading to Ca2+-dependent activation of endothelial nitric oxide (NO) synthase (eNOS). NO enhances mucociliary clearance and has direct antibacterial effects in ciliated epithelial cells and increases phagocytosis by macrophages. Using biochemistry and live cell imaging, we explored the role of heat shock protein 90 (HSP90) in regulating T2R-dependent NO pathways in primary sinonasal epithelial cells, primary monocyte-derived macrophages, and a human bronchiolar cell line (H441). We used immunofluorescence to show that H441 cells express eNOS and certain T2Rs and that the bitterant denatonium benzoate activates NO production in an HSP90-dependent manner in cells grown either as submerged cultures and at air liquid interface. In primary sinonasal epithelial cells, we determined that HSP-90 inhibition reduces T2R-stimulated NO production and ciliary beating which are crucial for pathogen clearance. In primary monocyte-derived macrophages, we found that HSP-90 is integral to T2R-stimulated NO production and phagocytosis of FITC-labeled Escherichia coli and pHrodo-Staphylococcus aureus. Our study demonstrates that HSP90 serves an innate immune role by regulating NO production downstream of T2R signaling by augmenting eNOS activation without impairing upstream calcium signaling. These findings suggest that HSP90 plays an important role in airway antibacterial innate immunity and may be an important target in airway diseases like chronic rhinosinusitis, asthma, or cystic fibrosis. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. A Bayesian model selection approach to mediation analysis Mediation analysis is a powerful tool for discovery of causal relationships. We describe a Bayesian model selection approach to mediation analysis that is implemented in our bmediatR software. Using simulations, we show that bmediatR performs as well or better than established methods including the Sobel test, while allowing greater flexibility in both model specification and in the types of inference that are possible. We applied bmediatR to genetic data from mice and human cell lines to demonstrate its ability to derive biologically meaningful findings. The Bayesian model selection framework is extensible to support a wide variety of mediation models. ScRNA-seq Expression of IFI27 and APOC2 Identifies Four Alveolar Macrophage Superclusters in Healthy BALF Alveolar macrophages (AMs) reside on the luminal surface of the airways and alveoli, ensuring proper gas exchange by ingesting cellular debris and pathogens, and regulating inflammatory responses. Therefore, understanding the heterogeneity and diverse roles played by AMs and recruited monocytes is critical for treating airway diseases. We performed single-cell RNA sequencing on 113,213 bronchoalveolar lavage cells from four healthy and three uninflamed cystic fibrosis subjects and identified four AM superclusters (families) based on the expression of IFI27 and APOC2 genes. These 4 AM families have at least eight distinct functional members (subclusters) named after their differentially expressed gene(s): IGF1, CCL18, CXCL5, Cholesterol, Chemokine, Metallothionein, Interferon and small-cluster AMs. Interestingly, the Chemokine cluster further divides with each subcluster selectively expressing a unique combination of chemokines. One of the most striking observations, besides the heterogeneity, is the conservation of AM family members in relatively equal ratio across all AM superclusters and individuals. Finally, transcriptional data and TotalSeq technology were used to investigate cell surface markers that distinguish resident AMs from recruited monocytes. Overall, functional specializations of the AM subclusters suggest that there are AM niches with defined programming states, highlighting a clear division of labor. Reading Modality Shapes Reading Network in Proficient Blind Readers of Braille The neural basis of reading is highly consistent across many languages and scripts. Are there alternative neural routes to reading? How does the sensory modality of symbols (tactile vs. visual) influence their neural representations? We examined these questions by comparing reading of visual print (sighted group, n=19) and tactile Braille (congenitally blind group, n=19). Blind and sighted readers were presented with written (words, consonant strings, non-letter shapes) and spoken stimuli (words, backward speech) that varied in word-likeness. Consistent with prior work, the ventral occipitotemporal cortex (vOTC) was active during Braille and visual reading. A posterior/anterior vOTC word-form gradient was observed only in sighted readers with more anterior regions preferring larger orthographic units (words). No such gradient was observed in blind readers. Consistent with connectivity predictions, in blind compared to sighted readers, posterior parietal cortices were recruited to a greater degree and contained word-preferring patches. Lateralization of Braille in blind readers was predicted by laterality of spoken language and reading hand. The effect of spoken language increased along a cortical hierarchy, whereas effect of reading-hand waned. These results suggested that the neural basis of reading is influenced by symbol modality and spoken language and support connectivity-based views of cortical function.