title: Influence of Nesting Material Composition on Tick Tube Use by Peromyscus leucopus

abstract: Host-targeted acaricides are a valuable tool for the reduction of ticks and tick-borne disease. Tick tubes (also known as tick control tubes) are commercially available products containing permethrin-treated nesting materials. Through superficial acaricide application to Peromyscus mice, tick tubes reduce populations of the blacklegged tick, Ixodes scapularis Say. Results of prior field trials have varied, suggesting that mouse behavior as well as the scale of the intervention and the composition of the local host community are important determinants of efficacy. Here we evaluated behaviors related to nest material collection by P. leucopus. Two forms of nest materials used in commercial tick tube products (cotton batting and balls) were assessed through side-by-side comparisons over a four-week period. We quantified cotton uptake by monitoring weekly changes in material weight and used video surveillance to categorize and assess mouse behaviors. The odds of cotton batting being taken from tubes was 2.14 times greater than cotton balls but the process was less efficient; mice removed 0.35 g less cotton batting for each removal event and required 2.17 times longer to complete the removal. While cotton balls were readily carried in the jaws of mice, batting required separating smaller fragments from the mass before placement in the oral cavity. Video surveillance suggested that a small number of mice were super users and responsible for 22% of the 119 visits in which material was removed. Combined, material weight loss and video-captured removal events improve our understanding of host usage of nest materials but also raise questions about dissemination of the material in nests of the local mouse community.

category: ecology


title: Fixation dynamics of beneficial alleles in prokaryotic polyploid chromosomes and plasmids

abstract: Theoretical population genetics has been mostly developed for sexually reproducing diploid and for monoploid (haploid) organisms, focusing on eukaryotes. The evolution of bacteria and archaea is often studied by models for the allele dynamics in monoploid populations. However, many prokaryotic organisms harbor multicopy replicons - chromosomes and plasmids - and theory for the allele dynamics in populations of polyploid prokaryotes remains lacking. Here we present a population genetics model for replicons with multiple copies in the cell. Using this model, we characterize the fixation process of a dominant beneficial mutation at two levels: the phenotype and the genotype. Our results show that, depending on the mode of replication and segregation, the fixation time of mutant phenotypes may precede the genotypic fixation time by many generations; we term this time interval the heterozygosity window. We furthermore derive concise analytical expressions for the occurrence and length of the heterozygosity window, showing that it emerges if the copy number is high and selection strong. Replicon ploidy thus allows for the maintenance of genetic variation following phenotypic adaptation and consequently for reversibility in adaptation to fluctuating environmental conditions.

category: evolutionary biology


title: Cancer-driving mutations are enriched in genic regions intolerant to germline variation

abstract: Large reference datasets of protein-coding variation in human populations have allowed us to determine which genes and genic sub-regions are intolerant to germline genetic variation. There is also a growing number of genes implicated in severe Mendelian diseases that overlap with genes implicated in cancer. Here, we hypothesized that mitotically mutable genic sub-regions that are intolerant to germline variation are enriched for cancer-driving mutations. We introduce a new metric, OncMTR, which uses 125,748 exomes in the gnomAD database to identify genic sub-regions intolerant to germline variation but enriched for hematologic somatic variants. We demonstrate that OncMTR can significantly predict driver mutations implicated in hematologic malignancies. Divergent OncMTR regions were enriched for cancer-relevant protein domains, and overlaying OncMTR scores on protein structures identified functionally important protein residues. Finally, we performed a rare variant, gene-based collapsing analysis on an independent set of 394,694 exomes from the UK Biobank and find that OncMTR dramatically improves genetic signals for hematologic malignancies. Our web app enables easy visualization of OncMTR scores for each protein-coding gene (https://astrazeneca-cgr-publications.github.io/OncMTR-Viewer/).

category: genomics


title: Chromosome-level genome assembly and annotation of two lineages of the ant Cataglyphis hispanica: steppingstones towards genomic studies of hybridogenesis and thermal adaptation in desert ants

abstract: Cataglyphis are thermophilic ants that forage during the day when temperatures are highest and sometimes close to their critical thermal limit. Several Cataglyphis species have evolved unusual reproductive systems such as facultative queen parthenogenesis or social hybridogenesis, which have not yet been investigated in detail at the molecular level. We generated high-quality genome assemblies for two hybridogenetic lineages of the Iberian ant Cataglyphis hispanica using long-read Nanopore sequencing and exploited chromosome conformation capture (3C) sequencing to assemble contigs into 26 and 27 chromosomes, respectively. Males of one lineage were karyotyped to confirm the number of chromosomes inferred from 3C data. We obtained transcriptomic data to assist gene annotation and built custom repeat libraries for each of the two assemblies. Comparative analyses with 19 other published ant genomes were also conducted. These new genomic resources pave the way for exploring the genetic mechanisms underlying the remarkable thermal adaptation and the molecular mechanisms associated with transitions between different genetic systems characteristics of the ant genus Cataglyphis.


title: RNA binding to human METTL3-METTL14 restricts N6-deoxyadenosine methylation of DNA in vitro

abstract: Methyltransferase like-3 (METTL3) and METTL14 complex transfers a methyl group from S-adenosyl-L-methionine to N6 amino group of adenosine bases in RNA (m6A) and DNA (m6dA). Emerging evidence highlights a role of METTL3-METTL14 in the chromatin context, especially in processes where DNA and RNA are held in close proximity. However, a mechanistic framework about specificity for substrate RNA/DNA and their interrelationship remain unclear. By systematically studying methylation activity and binding affinity to a number of DNA and RNA oligos with different propensities to form inter- or intra-molecular duplexes or single-stranded molecules in vitro, we uncover an inverse relationship for substrate binding and methylation and show that METTL3-METTL14 preferentially catalyzes the formation of m6dA in single-stranded DNA (ssDNA), despite weaker binding affinity to DNA. In contrast, it binds structured RNAs with high affinity, but methylates the target adenosine in RNA (m6A) much less efficiently than it does in ssDNA. We also show that METTL3-METTL14-mediated methylation of DNA is largely regulated by structured RNA elements prevalent in long noncoding and other cellular RNAs.

category: biochemistry


title: Identification of a phage display-derived peptide interacting with the N-terminal region of Factor VII activating protease (FSAP) enables characterization of zymogen activation.

abstract:
Increased Factor VII activating protease (FSAP) activity has a protective effect in diverse disease conditions as inferred from studies in FSAP-/- mice and humans deficient in FSAP activity due to a single nucleotide polymorphism. The activation of FSAP zymogen in plasma is mediated by extracellular histones that are released during tissue injury or inflammation or by positively charged surfaces. However, it is not clear if this activation mechanism is specific and amenable to manipulation. Using a phage display approach we have identified a peptide, NNKC9/41, that activates pro-FSAP in plasma. Other commonly found zymogens in the plasma were not activated. Binding studies with FSAP domain deletion mutants indicate that the N-terminus of FSAP is the key interaction site of this peptide. Blocking the contact pathway of coagulation did not influence pro-FSAP activation by the peptide. In a monoclonal antibody screen, we identified MA-FSAP-38C7 that prevented the activation of pro-FSAP by the peptide. This antibody bound to the LESLDP sequence (amino acids 30-35) in the N-terminus of FSAP. The plasma clotting time was shortened by NNKC9/41 and this was reversed by MA-FSAP-38C7 demonstrating the utility of this peptide. Identification of this peptide, and the corresponding interaction site, provides proof of principle that it is possible to activate a single protease zymogen in blood in a specific manner. Peptide NNKC/41 will be useful as a tool to delineate the molecular mechanism of activation of pro-FSAP in more detail, elucidate its biological role.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=105 SRC="FIGDIR/small/475526v1_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@112eea6org.highwire.dtl.DTLVardef@1f1d7eaorg.highwire.dtl.DTLVardef@967b03org.highwire.dtl.DTLVardef@1662c51_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Critical-like bistable dynamics in the resting-state human brain

abstract: Brain activity exhibits scale-free avalanche dynamics and power-law long-range temporal correlations (LRTCs) across the nervous system. This has been thought to reflect "brain criticality", i.e., brains operating near a critical phase transition between disorder and excessive order. Neuronal activity is, however, metabolically costly and may be constrained by activity-limiting mechanisms and resource depletion, which could make the phase transition discontinuous and bistable. Observations of bistability in awake human brain activity have nonetheless remained scarce and its functional significance unclear. First, using computational modelling where bistable synchronization dynamics emerged through local positive feedback, we found bistability to occur exclusively in a regime of critical-like dynamics. We then assessed bistability in vivo with resting-state magnetoencephalography and stereo-encephalography. Bistability was a robust characteristic of cortical oscillations throughout frequency bands from {delta} (3-7 Hz) to high-{gamma} (100-225 Hz). As predicted by modelling, bistability and LRTCs were positively correlated. Importantly, while moderate levels of bistability were positively correlated with executive functioning, excessive bistability was associated with epileptic pathophysiology and predictive of local epileptogenicity. Critical bistability is thus a salient feature of spontaneous human brain dynamics in awake resting-state and is both functionally and clinically significant. These findings expand the framework of brain criticality and show that critical-like neuronal dynamics in vivo involves both continuous and discontinuous phase transitions in a frequency-, neuroanatomy-, and state-dependent manner.

category: neuroscience


title: Robust Working Memory through Short-Term Synaptic Plasticity

abstract: Working memory has long been thought to arise from sustained spiking/attractor dynamics. However, recent work has suggested that short-term synaptic plasticity (STSP) may help maintain attractor states over gaps in time with little or no spiking. To determine if STSP endows additional functional advantages, we trained artificial recurrent neural networks (RNNs) with and without STSP to perform an object working memory task. We found that RNNs with and without STSP were both able to maintain memories over distractors presented in the middle of the memory delay. However, RNNs with STSP showed activity that was similar to that seen in the cortex of monkeys performing the same task. By contrast, RNNs without STSP showed activity that was less brain-like. Further, RNNs with STSP were more robust to noise and network degradation than RNNs without STSP. These results show that STSP can not only help maintain working memories, it also makes neural networks more robust.


title: Assessing Comparative Microbiome Performance in Plant Cell Wall Deconstruction Using Multi-omics-Informed Network Analysis

abstract: Plant cell walls are interwoven structures recalcitrant to degradation. Both native and adapted microbiomes are particularly effective at plant cell wall deconstruction. Studying these deconstructive microbiomes provides an opportunity to assess microbiome performance and relate it to specific microbial populations and enzymes. To establish a system assessing comparative microbiome performance, parallel microbiomes were cultivated on sorghum (Sorghum bicolor L. Moench) from compost inocula. Biomass loss and biochemical assays indicated that these microbiomes diverged in their ability to deconstruct biomass. Network reconstructions from time-dependent gene expression identified key deconstructive groups within the adapted sorghum-degrading communities, including Actinotalea, Filomicrobium, and Gemmanimonadetes populations. Functional analysis of gene expression demonstrated that the microbiomes proceeded through successional stages that are linked to enzymes that deconstruct plant cell wall polymers. This combination of network and functional analysis highlighted the importance of celluloseactive Actinobacteria in differentiating the performance of these microbiomes.

category: microbiology


title: TRF2 rescues pathogenic phenotypes in Duchenne muscular dystrophy cardiomyocytes derived from human iPSCs

abstract: Duchenne muscular dystrophy (DMD) is a severe muscle wasting disease caused by the lack of dystrophin. Heart failure, driven by cardiomyocyte death, fibrosis, and the development of dilated cardiomyopathy, is the leading cause of death in DMD patients. Current treatments decrease the mechanical load on the heart, but do not address the root cause of dilated cardiomyopathy: cardiomyocyte death. Previously, we showed that telomere shortening is a hallmark of cardiomyocytes in DMD patient cardiac tissues and in a mouse model of DMD we generated that manifests dilated cardiomyopathy. We also found that this telomere shortening is recapitulated in DMD cardiomyocytes differentiated from patient-derived induced pluripotent stem cells (iPSCs). Here we further characterize that DMD cardiomyocytes exhibit reduced cell size, nuclear size, and sarcomere density compared to healthy isogenic controls. The telomere-binding protein, TRF2, is a core component of the shelterin complex, which protects chromosome ends and is expressed at lower levels in DMD cardiomyocytes compared to controls. We investigated whether preservation of telomere length could be achieved with upregulation of TRF2. We found that TRF2 upregulation in DMD cardiomyocytes increased telomere lengths, cell size, nuclear size, sarcomere density, and cell survival. These data suggest TRF2 gene therapy has the potential to delay the onset of dilated cardiomyopathy.

category: cell biology


title: TRF2 rescues telomere attrition and prolongs cell survival in Duchenne muscular dystrophy cardiomyocytes derived from human iPSCs


title: SOX17 Is Not Required for the Derivation and Maintenance of Mouse Extraembryonic Endoderm Stem Cell Lines

abstract: Extraembryonic endoderm stem (XEN) cell lines can be derived and maintained in vitro and reflect the primitive endoderm cell lineage. SOX17 is thought to be required for the derivation and maintenance of mouse XEN cell lines. Here we have re-evaluated this requirement for SOX17. We derived multiple SOX17-deficient XEN cell lines from preimplantation embryos of a SOX17-Cre knockout strain and chemically converted multiple SOX17-deficient embryonic stem cell lines into XEN cell lines by transient culturing with retinoic acid and Activin A. We confirmed the XEN profile of SOX17-deficient cell lines by immunofluorescence with various markers, by NanoString gene expression analyses, and by their contribution to the extraembryonic endoderm of chimeric embryos produced by injecting these cells into blastocysts. Thus, SOX17 is not required for the derivation and maintenance of XEN cell lines.

category: developmental biology


title: Quantifying biological carbon pump pathways with a data-constrained mechanistic model ensemble approach

abstract: The ability to constrain the mechanisms that transport organic carbon into the deep ocean is complicated by the multiple physical, chemical, and ecological processes that intersect to create, transform, and transport particles in the ocean. In this manuscript we develop and parameterize a data-assimilative model of the multiple pathways of the biological carbon pump (NEMUROBCP). The mechanistic model is designed to represent sinking particle flux, active transport by vertically migrating zooplankton, and passive transport by subduction and vertical mixing, while also explicitly representing multiple biological and chemical properties measured directly in the field (including nutrients, phytoplankton and zooplankton taxa, carbon dioxide and oxygen, nitrogen isotopes, and 234Thorium). Using 30 different data types (including standing stock and rate measurements related to nutrients, phytoplankton, zooplankton, and non-living organic matter) from Lagrangian experiments conducted on 11 cruises from four ocean regions, we conduct an objective statistical parameterization of the model and generate one million different potential parameter sets that are used for ensemble model simulations. The model simulates in situ parameters that were assimilated (net primary production and gravitational particle flux) and parameters that were withheld (234Thorium and nitrogen isotopes) with reasonable accuracy. Model results show that gravitational flux of sinking particles and vertical mixing of organic matter from the surface ocean are more important biological pump pathways than active transport by vertically-migrating zooplankton. However, these processes are regionally variable, with sinking particles most important in oligotrophic areas of the Gulf of Mexico and California, sinking particles and vertical mixing roughly equivalent in productive regions of the CCE and the subtropical front in the Southern Ocean, and active transport an important contributor in the Eastern Tropical Pacific. We further find that mortality at depth is an important component of active transport when mesozooplankton biomasses are high, but that it is negligible in regions with low mesozooplankton biomass. Our results also highlight the high degree of uncertainty, particularly amongst mesozooplankton functional groups, that is derived from uncertainty in model parameters, with important implications from results that rely on non-ensemble model outputs. We also discuss the implications of our results for other data assimilation approaches.


title: Seasonal variability in carbon:234thorium ratios of suspended and sinking particles in coastal Antarctic waters: Field data and modeling synthesis

abstract:
238U-234Th disequilibrium is a powerful tool for investigating particle cycling and carbon export associated with the oceans biological carbon pump. However, the interpretation of this method is complicated by multiple processes that can modify carbon:thorium ratios over small spatial scales. We investigated seasonal variability in the thorium and carbon cycles at a coastal site in the Western Antarctic Peninsula. Throughout the ice-free summer season, we quantified carbon and 234Th vertical flux, total water column 234Th, particulate 234Th, and the C:234Th ratios of sinking material and bulk suspended material. Simultaneous identification and separation of fecal pellets from sinking material showed that fecal pellets (primarily from krill) contributed 56% of carbon flux and that as a result of lower C:234Th ratios than suspended particles, these fecal pellets were primary drivers of variability in the C:234Th ratios of sinking material. Bulk suspended particles had highly variable C:234Th ratios and were consistently elevated in the euphotic zone relative to deeper waters. The fraction of 234Th adsorbed onto particles was positively correlated with chlorophyll and particulate organic carbon (POC) concentrations. The C:234Th ratios of suspended particles were positively correlated with POC, although during the spring diatom bloom C:234Th ratios were lower than would have been predicted based on POC concentrations alone. We hypothesize that diatom production of transparent exopolymers may have led to enhanced rates of thorium adsorption during the bloom, thus decreasing the C:234Th ratios. We used a Bayesian model selection approach to develop and parameterize mechanistic models to simulate thorium sorption dynamics. The best model incorporated one slowly-sinking POC pool and rapidly-sinking fecal pellets, with second-order sorption kinetics. The model accurately simulated temporal patterns in the C:234Th ratios of sinking and suspended particles and the fraction of 234Th adsorbed to particles. However, it slightly over-estimated C:234Th ratios during the spring (diatom-dominated) bloom and underestimated C:234Th ratios during the fall (mixed-assemblage) bloom. Optimized model parameters for thorium sorption and desorption were 0.0047 {+/-} 0.0002 m3 mmol C-1 d-1 and 0.017 {+/-} 0.008 d-1, respectively. Our results highlight the important role that specific taxa can play in modifying the C:234Th ratio of sinking and suspended particles and provide guidance for future studies that use 234Th measurements to investigate the functional relationships driving the efficiency of the biological pump.

HIGHLIGHTSInvestigated thorium and carbon cycling over full ice-free season

C:234Th ratios of sinking particles were controlled by low C:234Th of fecal pellets

C:234Th ratios of suspended particles were correlated with chlorophyll and POC

Diatom abundance may have led to high particulate thorium during spring bloom

Second-order thorium sorption kinetics model accurately simulates C:234Th ratios


title: Reexamining Waddington: Canalization and new mutations are not required for the evolution of genetic assimilation

abstract: Over 65 years ago, Waddington demonstrated ancestrally phenotypically plastic traits can evolve to become constitutive, a process he termed genetic assimilation. Genetic assimilation evolves rapidly, assumed to be in large part due to segregating genetic variation only expressed in rare/novel environments, but otherwise phenotypically cryptic. Despite previous work suggesting a substantial role of cryptic genetic variation contributing to the evolution of genetic assimilation, some have argued for a prominent role for new mutations of large effect concurrent with selection. Interestingly, Waddington was less concerned by the relative contribution of CGV or new variants, but aimed to test the role of canalization, an evolved form of robustness. While canalization has been extensively studied, its role in the evolution of genetic assimilation is disputed, in part because explicit tests of evolved robustness are lacking. To address these questions, we recreated Waddingtons selection experiments on an environmentally sensitive change in Drosophila wing morphology (crossvein development), using many independently evolved replicate lineages. Using these, we show that 1) a polygenic CGV, but not new variants of large effect are largely responsible for the evolved response demonstrated using both genomic and genetic approaches. 2) Using both environmental manipulations and mutagenesis of the evolved lineages that there is no evidence for evolved changes in canalization contributing to genetic assimilation. Finally, we demonstrate that 3) CGV has potentially pleiotropic and fitness consequences in natural populations and may not be entirely "cryptic".


title: Sex-specific regulation of metabolic health and vertebrate lifespan by AMP biosynthesis

abstract: Energy homeostasis is disrupted with age, which then fuels multiple age-related pathologies. The AMP-activated protein kinase (AMPK) is the primary sensor of cellular energy in eukaryotes. However, the genetic regulation of vertebrate aging by AMPK remains poorly understood. Here, we manipulate energy levels in the turquoise killifish by mutating APRT, a key enzyme in AMP biosynthesis. These manipulations produced a male-specific lifespan extension and restored metabolic plasticity. Exploring the observed sex differences using an integrated omics approach implicated the mitochondria as an important player. Mechanistically, APRT regulated mitochondrial functions and AMPK activity, mimicking energy starvation in heterozygous cells. A fasting-like state was also detected, particularly in heterozygous males, which leads to resistance to high-fat diet. Finally, life-long intermittent fasting eliminated the male-specific longevity benefits mediated by the APRT mutation. These observations identify the AMP/AMPK axis as a sex-specific regulator of vertebrate longevity and metabolic health.

category: genetics


title: A pleiotropic chemoreceptor facilitates the functional coupling of pheromone production and perception

abstract: Optimal mating decisions depend on the robust coupling of signal production and perception because independent changes in either could carry a fitness cost. However, since the perception and production of mating signals are often mediated by different tissues and cell types, the mechanisms that drive and maintain their coupling remain unknown for most animal species. Here, we show that in Drosophila, sensory perception and production of an inhibitory mating pheromone are co-regulated by Gr8a, a member of the Gustatory receptor gene family. Specifically, we found that the pleiotropic action of Gr8a independently regulates the perception of pheromones by the chemosensory systems of males and females, as well as their production in the fat body and oenocytes of males. These findings provide a relatively simple molecular explanation for how pleiotropic receptors maintain robust mating signaling systems at the population and species levels.


title: The genomic landscape of contemporary western Remote Oceanians

abstract:
The Vanuatu archipelago served as a gateway to Remote Oceania during one of the most extensive human migrations to uninhabited lands ~3,200 years ago. Ancient DNA studies suggest an initial settlement by East Asian-related peoples that was quickly followed by the arrival of Papuan-related populations, leading to a major population turnover. Yet, there is uncertainty over the population processes and the sociocultural factors that have shaped the genomic diversity of ni-Vanuatu, who present nowadays among the worlds highest linguistic and cultural diversity. Here, we report genome-wide data for 1,433 contemporary ni-Vanuatu from 29 different islands, including 287 couples. We find that ni-Vanuatu derive their East Asian- and Papuan-related ancestry from the same source populations and descend from relatively synchronous admixture events that occurred ~1,700-2,300 years ago, indicating a peopling history common to all the archipelago. However, our analyses reveal that the Papuan population turnover was geographically uneven, and that the genetic contribution of Papuan-related peoples was male-biased. Furthermore, we detect Polynesian ancestry arriving -600-1,000 years ago to South Vanuatu, and map its distribution to both Polynesian- and non-Polynesian-speaking islands. Lastly, we provide evidence for a tendency of spouses to carry similar genetic ancestry, when accounting for relatedness avoidance. The signal is not driven by strong genetic effects of specific loci or trait-associated variants, suggesting that it results instead from social assortative mating. Altogether, our findings provide insight into both the genetic history of ni-Vanuatu populations and how sociocultural processes have shaped the diversity of their genomes.

Significance statementVanuatu served as a gateway for the peopling of Remote Oceania but how demography and culture have shaped the genomic diversity of its peoples remains partially understood. Here, we report genomic data for 1,433 contemporary ni-Vanuatu covering 29 islands. We demonstrate that their genomic diversity derives from relatively synchronous, geographically uneven and sex-biased admixture events between the Austronesian-related first settlers and the Papuan-related subsequent incomers. We show that the more recent arrival of Polynesian ancestry did not always trigger shifts to Polynesian languages. Lastly, we find evidence for non-random mating, probably owing to social assortative mating. Our findings reveal the impact of multiple migrations, admixture and sociocultural practices on the genomic landscape of a key region of the South Pacific.


title: Evolutionarily stable gene clusters shed light on the common grounds of pathogenicity in the Acinetobacter calcoaceticus-baumanniicomplex

abstract:
Nosocomial pathogens of the Acinetobacter calcoaceticus-baumannii (ACB) complex are a cautionary example for the world-wide spread of multi- and pan-drug resistant bacteria. Aiding the urgent demand for novel therapeutic targets, comparative genomics studies between pathogens and their apathogenic relatives shed light on the genetic basis of human-pathogen interaction. Yet, existing studies are limited in taxonomic scope, sensing of the phylogenetic signal, and resolution by largely analyzing genes isolated from their functional contexts. Here, we explored more than 3,000 Acinetobacter genomes in a phylogenomic framework integrating orthology-based phylogenetic profiling and microsynteny conservation analyses. This allowed to delineate gene clusters in the type strain A. baumannii ATCC 19606 whose evolutionary conservation indicates a functional integration of the subsumed genes. These evolutionarily stable gene clusters (ESGCs) reveal metabolic pathways, transcriptional regulators residing next to their targets but also tie together sub-clusters with distinct functions to form higher-order functional modules. We shortlisted 150 ESGCs that either co-emerged with, or are found preferentially in, the pathogenic ACB clade. They unveil, at an unprecedented resolution, the genetic makeup that coincides with the manifestation of the pathogenic phenotype in the last common ancestor of the ACB clade. Key innovations are the remodeling of the regulatory-effector cascade connecting LuxR/LuxI quorum sensing via an intermediate messenger to biofilm formation, the extension of micronutrient scavenging systems, and the increase of metabolic flexibility by exploiting carbon sources that are provided by the human host. Specifically, we could show that only members of the ACB clade use kynurenine as a sole carbon and energy source, a substance produced by humans to fine-tune the antimicrobial innate immune response. In summary, this study provides a rich and unbiased set of novel testable hypotheses on how pathogenic Acinetobacter interact with and ultimately infect their human host. They disclose promising routes for future therapeutic strategies.

Author SummaryThe spread of multi- and pan-drug resistant bacterial pathogens is a world-wide threat to human health. Understanding the genetics of host colonization and infection is essential for devising novel ways of treatment. Acinetobacter baumannii, a nosocomial pathogen ranked top by the World Health Organization in the list of bacteria for which novel therapeutic approaches are needed, is a prime example. Here, we have carved out the genetic make-up that distinguishes A. baumannii and its pathogenic next relatives comprising the A. baumannii-calcoaceticus complex from other and mostly apathogenic Acinetobacter species. We found a rich spectrum of pathways and regulatory modules that reveals how the pathogens have modified biofilm formation, iron scavenging, and their carbohydrate metabolism to adapt to their human host. Among these, the capability to metabolize kynurenine is particularly intriguing. Humans produce this substance to contain bacterial invaders and to fine-tune the innate immune response. But A. baumannii and its allies have turned the table. They feed on kynurenine, and thus most likely also dysregulate human immune response. In summary, our study substantially deepens the understanding of how a highly critical pathogen interacts with its host, which substantially eases the identification of novel targets for innovative therapeutic strategies.


title: Genomic features underlie the co-option of SVA transposons as cis-regulatory elements in human pluripotent stem cells

abstract: Domestication of transposable elements (TEs) into functional cis-regulatory elements is a widespread phenomenon. However, the mechanisms behind why some TEs are co-opted as functional enhancers while others are not are underappreciated. SINE-VNTR-Alus (SVAs) are the youngest group of transposons in the human genome, where [~]3,700 copies are annotated, nearly half of which are human-specific. Many studies indicate that SVAs are among the most frequently co-opted TEs in human gene regulation, but the mechanisms underlying such processes have not yet been thoroughly investigated. Here, we leveraged CRISPR-interference (CRISPRi), computational and functional genomics to elucidate the genomic features that underlie SVA domestication into human stem-cell gene regulation. We found that [~]750 SVAs are co-opted as functional cis-regulatory elements in human induced pluripotent stem cells. These SVAs are significantly closer to genes and harbor more transcription factor binding sites than non-co-opted SVAs. We show that a long DNA motif composed of flanking YY1/2 and OCT4 binding sites is enriched in the co-opted SVAs and that these two transcription factors bind consecutively on the TE sequence. We used CRISPRi to epigenetically repress active SVAs in stem cell-like NCCIT cells. Epigenetic perturbation of active SVAs strongly attenuated YY1/OCT4 binding and influenced neighboring gene expression. Ultimately, SVA repression resulted in [~]3,000 differentially expressed genes, 131 of which were the nearest gene to an annotated SVA. In summary, we demonstrated that SVAs modulate human gene expression, and uncovered that location and sequence composition contribute to SVA domestication into gene regulatory networks.


title: Evaluation of an optimized protocol and Illumina ARTIC V4 primer pool for sequencing of SARS-CoV-2 using COVIDSeq™ and DRAGEN™ COVID Lineage App workflow

abstract: Next-Generation Sequencing based genomic surveillance has been widely implemented for identification and tracking of emerging SARS-CoV-2 variants to guide the Public Health response to the COVID-19 pandemic. Amplicon-based assays, such as the Illumina(R) COVIDSeq Test (RUO) and COVIDSeq Assay (RUO), enable scalable sequencing of SARS-CoV-2, leveraging V3 and V4 primer designs from the ARTIC community and DRAGEN COVID Lineage App analysis available on Illumina BaseSpace. We report here a comparison of COVIDSeq performance for SARS-CoV-2 genome reporting using the ARTIC V3 based primer pool (including primers for human control genes) that is provided with the COVIDSeq kit versus the ARTIC V4 based Illumina COVIDSeq V4 primer pool, using an optimized protocol and DRAGEN COVID Lineage App analysis. The data indicates that both primer pools enable robust reporting of SARS-CoV-2 variants. The Illumina COVIDSeq V4 primer pool has superior performance for SARS-CoV-2 genome reporting, particularly in samples with low virus load, and is therefore the recommended primer pool for genomic surveillance of SARS-CoV-2 for research use using COVIDSeq.


title: A simple bypass assay for DNA polymerases shows hypermutating variants associated with cancer show mechanistic differences in vitro

abstract: Errors made by DNA polymerases contribute to both natural variation and, in extreme cases, to genome instability and its associated diseases. Recently the importance of polymerase misincorporation in disease has been highlighted by the identification of cancer-associated polymerase variants and the recognition that a subgroup of these variants have a hypermutation phenotype in tumours. We have developed a bypass assay to rapidly determine the tendency of a polymerase to misincorporate in vitro. We have used the assay to compare misincorporation by wild-type, exonuclease defective and two hypermutating DNA polymerase e variants, P286R and V411L. The assay clearly distinguished between the misincorporation rates of wild type, exonuclease dead and P286R polymerases. However, the V411L polymerase showed different misincorporation characteristics to P286R, suggesting that these variants cause hypermutation by different mechanisms. Using this assay misincorporation opposite a templated C nucleotide was consistently higher than for other nucleotides, and this caused predominantly C to T transitions. This is consistent with the observation that C to T transitions are commonly seen in POLE mutant tumours.


title: Fluorescence visualization of deep-buried hollow organs

abstract: High-definition fluorescence imaging of deep-buried organs is still challenging. Here, we develop bright fluorophores emitting to 1700 nm by enhancing electron donating ability and reducing donor-acceptor distance. In parallel, the heavy water functions as the solvent of the delicately designed fluorophores, effectively reducing the fluorescent signal loss caused by the absorption by water. The near-infrared-II (NIR-II, 900-1880 nm) emission is eventually recovered and extended beyond 1400 nm. Compared with the spectral range beyond 1500 nm, the one beyond 1400 nm gives a more accurate fluorescence visualization of the hollow organs, owing to the absorption-induced scattering suppression. In addition, the intraluminal lesions containing much water are simultaneously negatively stained, leading to a stark contrast for precise diagnosis. Eventually, the intraluminally perfused fluorescent probes are excreted from mice and thus no obvious side effects emerge. This general method can provide new avenues for future biomedical imaging of deep and highly scattering tissues.

category: bioengineering


title: Application of sequential cyclic compression on cancer cells in a flexible microdevice

abstract: Mechanical forces shape physiological structure and function within cell and tissue microenvironments, during which cells strive to restore their shape or develop an adaptive mechanism to maintain cell integrity depending on strength and type of the mechanical loading. While some cells are shown to experience permanent plastic deformation after a repetitive mechanical tensile loading and unloading, the impact of such repetitive compression on plastic deformation of cells is yet to be understood. As such, the ability to apply cyclic compression is crucial for any experimental setup aimed at the study of mechanical compression taking place in cell and tissue microenvironments. Here, we demonstrate such cyclic compression using a microfluidic compression platform on live cell actin in SKOV-3 ovarian cancer cells. Live imaging of the actin cytoskeleton dynamics of the compressed cells was performed for varying pressures applied sequentially in ascending order during cell compression. Additionally, recovery of the compressed cells was investigated by capturing actin cytoskeleton and nuclei profiles of the cells at zero time and 24 h-recovery after compression in end point assays. This was performed for a range of mild pressures within the physiological range. Results showed that the phenotypical response of compressed cells during recovery after compression with 20.8 kPa differed observably from that for 15.6 kPa. This demonstrated the ability of the platform to aid in the capture of minute differences in cell behaviour as a result of being compressed at various pressures in physiologically relevant manner. Differences observed between compressed cells fixed at zero time or after 24 h-recovery suggest that SKOV-3 cells exhibit deformations at the time of the compression, a proposed mechanism cells use to prevent mechanical damage. The extent of recovery of compressed cells, inferred by imaging cell membrane bulges and actin cytoskeleton and measuring the shape descriptors of cell nuclei, can give insights into the plasticity of cancer cells. Thus, biomechanical responses of SKOV-3 ovarian cancer cells to sequential cyclic compression and during recovery after compression could be revealed in a flexible microdevice. As demonstrated in this work, the microfluidic platform could control the strength and duration of cyclic compression, while enabling the observation of morphological, cytoskeletal and nuclear changes and recovery in compressed cells, thus providing a powerful new tool for the study of mechanobiological processes in cancer and cell biology.


title: Single-Cell Proteomics Study of Immune Cell Diversity by Quantitating 465 Proteins

abstract: The identification and characterization of T cell subpopulations is critical to reveal cell development throughout life and immune responses to environmental factors. Next-generation sequencing technologies have dramatically advanced the single-cell genomics and transcriptomics for T cell classification. However, gene expression is often not correlated with protein expression, and immunotyping is mostly accepted in the protein format. Current single-cell proteomics technologies are either limited in multiplex capacity or not sensitive enough to detect the critical functional proteins. Herein we present a cyclic multiplex in situ tagging (Cyclic MIST) technology to simultaneously measure 465 proteins, a scale of >10 times than similar technologies, in single cells. Such a high multiplexity is achieved by reiterative staining of the single cells coupled with a MIST array for detection. This technology has been thoroughly validated through comparison with flow cytometry and fluorescence immunostaining techniques. Both THP1 and CD4+ T cells are analyzed by the Cyclic MIST technology, and over 300 surface markers have been profiled to classify the subpopulations. This represents the most comprehensive mapping of the diversity of immune cells at the protein level. With additional information from intracellular proteins of the same single cells, our technology can potentially facilitate mechanistic studies of immune responses, particularly cytokine storm that results in sepsis.


title: Retinal axial motion analysis and implications for real-time correction in human retinal imaging

abstract: High-resolution ophthalmic imaging devices including spectral-domain and full-field optical coherence tomography (SDOCT and FFOCT) are adversely affected by the presence of continuous involuntary retinal axial motion. Here, we thoroughly quantify and characterize retinal axial motion with both high temporal resolution (200,000 A-scans/s) and high axial resolution (4.5 {micro}m), recorded over a typical data acquisition duration of 3 s with an SDOCT device over 14 subjects. We demonstrate that although breath-holding can help decrease large-and-slow drifts, it increases small-and-fast fluctuations, which is not ideal when motion compensation is desired. Finally, by simulating the action of an axial motion stabilization control loop, we show that a loop rate of 1.2 kHz is ideal to achieve 100% robust clinical in-vivo retinal imaging.


title: Putative host-derived insertions in the genome of circulating SARS-CoV-2 variants

abstract:
Insertions in the SARS-CoV-2 genome have the potential to drive viral evolution, but the source of the insertions is often unknown. Recent proposals have suggested that human RNAs could be a source of some insertions, but the small size of many insertions makes this difficult to confirm. Through an analysis of available direct RNA sequencing data from SARS-CoV-2 infected cells, we show that viral-host chimeric RNAs are formed through what are likely stochastic RNA-dependent RNA polymerase template switching events. Through an analysis of the publicly available GISAID SARS-CoV-2 genome collection, we identified two genomic insertions in circulating SARS-CoV-2 variants that are identical to regions of the human 18S and 28S rRNAs. These results provide direct evidence of the formation of viral-host chimeric sequences and the integration of host genetic material into the SARS-CoV-2 genome, highlighting the potential importance of host-derived insertions in viral evolution.

IMPORTANCEThroughout the COVID-19 pandemic, the sequencing of SARS-CoV-2 genomes has revealed the presence of insertions in multiple globally circulating lineages of SARS-CoV-2, including the Omicron variant. The human genome has been suggested to be the source of some of the larger insertions, but evidence for this kind of event occurring is still lacking. Here, we leverage direct RNA sequencing data and SARS-CoV-2 genomes to show host-viral chimeric RNAs are generated in infected cells and two large genomic insertions have likely been formed through the incorporation of host rRNA fragments into the SARS-CoV-2 genome. These host-derived insertions may increase the genetic diversity of SARS-CoV-2 and expand its strategies to acquire genetic materials, potentially enhancing its adaptability, virulence, and spread.

category: bioinformatics


title: Systematic comparison of ranking aggregation methods for gene lists in experimental results

abstract: A common experimental output in biomedical science is a list of genes implicated in a given biological process or disease. The results of a group of studies answering the same, or similar, questions can be combined by meta-analysis to find a consensus or a more reliable answer. Ranking aggregation methods can be used to combine gene lists from various sources in meta-analyses. Evaluating a ranking aggregation method on a specific type of dataset before using it is required to support the reliability of the result since the property of a dataset can influence the performance of an algorithm. Evaluation of aggregation methods is usually based on a simulated database especially for the algorithms designed for gene lists because of the lack of a known truth for real data. However, simulated datasets tend to be too small compared to experimental data and neglect key features, including heterogeneity of quality, relevance and the inclusion of unranked lists. In this study, a group of existing methods and their variations which are suitable for meta-analysis of gene lists are compared using simulated and real data. Simulated data was used to explore the performance of the aggregation methods as a function of emulating the common scenarios of real genomics data, with various heterogeneity of quality, noise level, and a mix of unranked and ranked data using 20000 possible entities. In addition to the evaluation with simulated data, a comparison using real genomic data on the SARS-CoV-2 virus, cancer (NSCLC), and bacteria (macrophage apoptosis) was performed. We summarise our evaluation results in terms of a simple flowchart to select a ranking aggregation method for genomics data.


title: A visual atlas of genes tissue-specific pathological roles

abstract:
Dysregulation of a genes function, either due to mutations or impairments in regulatory networks, often triggers pathological states in the affected tissue. Comprehensive mapping of these apparent gene-pathology relationships is an ever daunting task, primarily due to genetic pleiotropy and lack of suitable computational approaches. With the advent of high throughput genomics platforms and community scale initiatives such as the Human Cell Landscape (HCL) project [1], researchers have been able to create gene expression portraits of healthy tissues resolved at the level of single cells. However, a similar wealth of knowledge is currently not at our finger-tip when it comes to diseases. This is because the genetic manifestation of a disease is often quite heterogeneous and is confounded by several clinical and demographic covariates. To circumvent this, we mined ~18 million PubMed abstracts published till May 2019 and selected ~6.1 million of them that describe the pathological role of genes in different diseases. Further, we employed a word embedding technique from the domain of Natural Language Processing (NLP) to learn vector representation of entities such as genes, diseases, tissues, etc., in a way such that their relationship is preserved in a vector space. Notably, Pathomap, by the virtue of its underpinning theory, also learns transitive relationships. Pathomap provided a vector representation of words indicating a possible association between DNMT3A/BCOR with CYLD cutaneous syndrome (CCS). The first manuscript reporting this finding was not part of our training data.

Key pointsO_LIWe mined ~18 million PubMed abstracts to extract latent knowledge pertaining to tissue specific pathological roles of genes.
C_LIO_LIWe found well-defined gene modules implicated in disease pathogenesis in anatomically proximal tissues.
C_LIO_LIWe demonstrated an ahead of time discovery of the association between DNMT3A/BCOR with CYLD cutaneous syndrome (CCS), as a knowledge synthesis use-case.
C_LI


title: A graph-embedded topic model enables characterization of diverse pain phenotypes among UK Biobank individuals

abstract: Large biobank repositories of clinical conditions and medications data open opportunities to investigate the phenotypic disease network. To enable systematic investigation of entire structured phenomes, we present graph embedded topic model (GETM). Our contributions are two folds in terms of method and applications. On the methodology side, we offer two main contributions in GETM. First, to aid topic inference, we integrate existing biomedical knowledge graph information in the form of pre-trained graph embedding into the embedded topic model. Second, leveraging deep learning techniques, we developed a variational autoencoder framework to infer patient phenotypic mixture by modeling multi-modal discrete patient medical records. In particular, for interpretability, we use a linear decoder to simultaneously infer the bi-modal distributions of the disease conditions and medications. On the application side, we applied GETM to UK Biobank (UKB) self-reported clinical phenotype data, which contains 443 self-reported medical conditions and 802 self-reported medications for 457,461 individuals. Compared to existing methods, GETM demonstrates overall superior performance in imputing missing conditions and medications. Here, we focused on characterizing pain phenotypes recorded in the questionnaire of the UKB individuals. GETM accurately predicts the status of chronic musculoskeletal (CMK) pain, chronic pain by body-site, and non-specific chronic pain using past conditions and medications. Our analyses revealed not only the known pain-related topics but also the surprising predominance of medications and conditions in the cardiovascular category among the most predictive topics across chronic pain phenotypes.


title: Uncovering One-Dimensional Reaction Coordinate that Underlies Structure-Function Relationship of Proteins

abstract: Understanding the mechanism of functional protein dynamics is critical to understanding protein functions. Reaction coordinates is a central topic in protein dynamics and the grail is to find the one-dimensional reaction coordinate that can fully determine the value of committor (i.e. the reaction probability in configuration space) for any protein configuration. We present a powerful new method that can, for the first time, identify the rigorous one-dimensional reaction coordinate in complex molecules. This one-dimensional reaction coordinate is determined by a fundamental mechanical operator--the generalized work functional. This method only requires modest computational cost and can be readily applied to large molecules. Most importantly, the generalized work functional is the physical origin of the collectivity in functional protein dynamics and provides a tentative roadmap that connects the structure of a protein to its function.

category: biophysics


title: Quantification of salt stress in wheat leaves by Raman spectroscopy and machine learning

abstract: The salinity level of the growing medium has diverse effects on the development of plants, including both physical and biochemical changes. To determine the salt stress level of a plant endures, one can measure these structural and chemical changes. Raman spectroscopy and biochemical analysis are some of the most common techniques in the literature. Here, we present a combination of machine learning and Raman spectroscopy with which we can both find out the biochemical change that occurs while the medium salt concentration changes and predict the level of salt stress a wheat sample experiences accurately using our trained regression models. In addition, by applying different machine learning algorithms, we compare the level of success for different algorithms and determine the best method to use in this application. Production units can take actions based on the quantitative information they get from the trained machine learning models related to salt stress, which can potentially increase efficiency and avoid the loss of crops.


title: Deviation from the Residue-Based Linear Free Energy Relationship Reveals a Non-Native Structure in Protein Folding

abstract:
Multiprobe measurements, such as NMR and hydrogen exchange study, can provide the equilibrium constant K and kinetic rate constant k of the structural changes of a polypeptide on a per-residue basis. We previously found a linear relationship between residue-specific log K values and residuespecific log k values for the two-state topological isomerization of a 27-residue peptide. To test the general applicability of the residue-based linear free energy relationship (rbLEFR), we performed a literature search to collect residue-specific equilibrium and kinetic constants in various exchange processes, including protein folding, coupled folding and binding of intrinsically disordered peptides, and structural fluctuations of folded proteins. The good linearity in a substantial number of log-log plots proved that the rbLFER holds for the structural changes in a wide variety of protein-related phenomena. Protein molecules quickly fold into their native structures and change their conformations smoothly. Theoretical studies and molecular simulations advocate that the physicochemical basis is the consistency principle and the minimal frustration principle: Non-native structures/interactions are absent or minimized along the folding pathway. The linearity of the residue-based free energy relationship demonstrates experimentally the absence of non-native structures in transition states. In this context, the hydrogen exchange study of apomyoglobin folding intermediates is particularly interesting. We found that the residues that deviated from the linear relationship corresponded to the non-native structure, which had been identified by other experiments. The rbLFER provides a unique and practical method to probe the dynamic aspects of the transition states of protein molecules.

HighlightsO_LIA collection of equilibrium and kinetic constants of structural changes of proteins
C_LIO_LIResidue-based linear free energy relationship widely holds between the two constants
C_LIO_LIrbLFER indicates the absence of non-ground state structures in transition states
C_LIO_LIrbLFER is an experiment proof of the consistency principle of protein folding
C_LIO_LIDeviations from the linear relation suggest special structures in transition states
C_LI

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY


title: Characterizing protein protonation microstates using Monte Carlo sampling

abstract: Proteins are polyelectrolytes with acidic or basic amino acids making up {approx}25% of the residues. The protonation state of all Asp, Glu, Arg, Lys, His and other protonatable residues, cofactors and ligands define each protonation microstate. As all of these residues will not be fully ionized or neutral, proteins exist in a mixture of microstates. The microstate distribution changes with pH. As the protein environment modifies the proton affinity of each site the distribution may also change in different reaction intermediates or as ligands are bound. Particular protonation microstates may be required for function, while others exist simply because there are many states with similar energy. Here, the protonation microstates generated in Monte Carlo sampling in MCCE are characterized in HEW lysozyme as a function of pH and bacterial photosynthetic reaction centers (RCs) in different reaction intermediates. The lowest energy and highest probability microstates are compared. The {Delta}G, {Delta}H and {Delta}S between the four protonation states of Glu35 and Asp52 in lysozyme are shown to be calculated with reasonable precision. A weighted Pearson correlation analysis identifies coupling between residue protonation states in RCs and how they change when the quinone in the QB site is reduced.


title: Radiofrequency ablation remodels the tumor microenvironment and promotes systemic immunomodulation in pancreatic cancer

abstract:
Background and AimsPancreatic ductal adenocarcinoma (PDAC) is characterized by resistance to therapy. A major contributing factor to therapeutic failure is profound desmoplasia and a well-documented hypoxic tumor microenvironment (TME). In PDAC, several therapeutic approaches, including chemotherapy and radiation alone or combined with immune checkpoint inhibitors, have shown minimal therapeutic success, placing an imperative need for the discovery and application of innovative treatments. Endoscopic ultrasound guided radiofrequency ablation (EUS-RFA) is a promising immunomodulator therapy for PDAC. In this work, we hypothesized RFA promotes local and systemic stromal and immunomodulating effects that can be identified for new combination therapeutic strategies.

MethodsTo test our hypothesis, a syngeneic PDAC mouse model was performed by symmetrically injecting 100k murine KPC cells in bilateral flanks of C57BL/6 female mice. RFA treatment initiated when tumors reached 200-500 mm3 and was performed only in the right flank. The left flank tumor (non-RFA contralateral side) was used as a paired control for further analysis.

ResultsRFA promoted a significant reduction in tumor growth rate 4 days after treatment in RFA treated and non-RFA side contralateral tumors from treated mice when compared to controls. Histological analysis revealed a significant increase in expression of cleaved Caspase3 in RFA treated tumors. In addition, collagen deposition and CD31+ cells were significantly elevated in RFA side and non-RFA contralateral tumors from RFA treated mice. Proteome profiling showed changes in C5a and IL-23 in RFA responsive tumors, indicating a role of RFA in modulating intratumoral inflammatory responses.

ConclusionsThese data indicate RFA promotes local and systemic anti-tumor responses in a syngeneic mouse model of PDAC implicating RFA treatment for local tumors as well as metastatic disease.

Graphical Abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

category: cancer biology


title: Peptide-centric analyses of human plasma enable increased resolution of biological insights into non-small cell lung cancer relative to protein-centric analysis

abstract: Comprehensive assessment of the human proteome remains challenging due to multiple forms of a protein, or proteoforms, arising from alternative splicing, allelic variation, and protein modifications. As proteoforms can serve distinct functions and act as functional links between genotype and phenotype, proteoform-level knowledge is critical in understanding the molecular mechanisms underlying health and disease. However, identification of proteoforms requires unbiased protein coverage at amino acid resolution. Scalable, deep, and unbiased proteomics studies have been impractical due to cumbersome and lengthy workflows required for complex samples, like blood plasma. Here, we demonstrate the power of the Proteograph Product Suite in enabling unbiased, deep, and rapid proteomics at scale in a proof-of-concept proteoform analysis to dissect differences between protein isoforms in plasma samples from 80 healthy controls and 61 patients with early-stage non-small-cell lung cancer (NSCLC). Processing the 141 plasma samples with Proteograph yielded 22,993 peptides corresponding to 2,569 protein groups at a confidence of 1% false discovery rate. We extracted four proteins with peptides with significant abundance differences (p < 0.05; Benjamini-Hochberg corrected) in healthy control and cancer plasma samples. For one, the abundance variation can be explained by underlying annotated protein isoforms. For a second, we find evidence for differentially transcribed isoforms in the broader sequence data, but not in the known annotated protein isoforms. The others may be explained by novel isoforms or post-translational modifications. In addition, we sought to identify protein variants arising from allelic variation. To this end, we first performed whole exome sequencing on buffy coat samples from 29 individuals in the NSCLC study. Then, we created personalized mass spectrometry search databases for each individual subject from the exome sequences. From these libraries, we identified 422 protein variants, one of which has previously been shown to relate to lung cancer. In conclusion, our results demonstrate that Proteograph can generate unbiased and deep plasma proteome profiles that enable identification of proteoforms present in plasma at a scale sufficient to enable population-scale proteomic studies powered to reveal novel mechanistic and biomedical insights.


title: Diurnal stomatal apertures and density ratios affect whole-canopy stomatal conductance, water-use efficiency and yield

abstract: Key physiological traits of plants, such as transpiration and stomatal conductance, are usually studied under steady-state conditions or modeled using only a few measured data points. Those measurements do not reflect the dynamic behavior of the plant in response to field conditions. To overcome this bottleneck, we used a gravimetric functional-phenotyping platform and a reverse-phenotyping method to examine the dynamic whole-plant water-regulation responses of tomato introgression lines and compared those responses with several years of yield performance in commercial fields. Ideotype lines had highly plastic stomatal conductance and high abaxial-adaxial stomatal density ratios and the size of their stomatal apertures peaked early in the day under water-deficit conditions. These traits resulted in dynamic daily water-use efficiency, which allowed for the rapid recovery of transpiration when irrigation was resumed after a period of imposed drought. We found that stomatal density, the abaxial-adaxial stomatal density ratio and the time of maximum stomatal apertures are crucial for plant adaptation and productivity under drought-stress conditions. Abaxial stomatal density was also found to be strongly correlated with the expression of the stomatal-development genes SPCH and ZEP. This study demonstrates how a reverse functional phenotyping approach based on field yield data, continuous and simultaneous whole-plant water-balance measurements and anatomical examination of individual leaves can help us to understand and identify dynamic and complex yield-related physiological traits.

category: plant biology


title: ARF small GTPases in the developmental function mediated by ARF regulators GNOM and VAN3

abstract: ARF small GTPases are molecular switches acting in intracellular trafficking. Their cycles of activity are controlled by regulators, ARF Guanine nucleotide Exchange Factors (ARF-GEFs) and ARF GTPase Activating Proteins (ARF-GAPs). The ARF-GEF GNOM (GN) and the ARF-GAP VAN3 share a prominent function in auxin-mediated developmental patterning, but the ARFs which they might control were not identified. We conducted a loss-of-function and localization-based screening of the ARF/ARF-LIKE gene family in Arabidopsis thaliana with the primary aim of identifying functional partners of GN and VAN3, while extending the limited understanding of this gene group as a whole. We identified a function of ARLA1 in branching angle control. Mutants lacking the variably localized ARLB1, ARFB1, ARFC1, ARFD1, and ARF3, even in high order combinations, do not exhibit any evident phenotypes. Loss of function arfa1 phenotypes support a major role of ARFA1 in growth and development overall, but patterning defects typical to gn loss of function are not found. ARFA1 are not localized at the plasma membrane, where GN and VAN3 carry out developmental patterning function according to current models. Taken together, putative ARF partners of GN and VAN3 in developmental patterning cannot be conclusively identified.


title: The non-neutral, multi-level, phenotypic impact of synonymous mutations revealed through heterologous gene expression in human cells.

abstract:
Redundancy in the genetic code allows for differences in transcription and/or translation efficiency between sequences carrying synonymous polymorphisms, potentially leading to phenotypic changes. It is commonly admitted that the evolution codon usage bias ("CUB", the over-representation of certain codons in a genome, a gene or in positions along a gene) are driven by a combination of neutral and selective processes, but their relative contribution is a matter of debate, especially in mammals. Particularly, integrative studies quantifying the phenotypic impact of CUB at different molecular and cellular levels are lacking. Here we report a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells.

Six synonymous versions of the shble antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was assessed by means of: i) mRNA-to-DNA and protein-to-mRNA ratios for each shble version; ii) cellular fluorescence, using flow cytometry, as a proxy for single-cell level expression; and iii) real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness.

We show that differences in CUB strongly impact the molecular and the cellular phenotype: i) they result in large differences in mRNA and in protein levels, as well in mRNA-to-protein ratio; ii) they introduce splicing events not predicted by current algorithms; iii) they lead to reproducible phenotypic heterogeneity; iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression.

We interpret that CUB modulate mRNA availability and suitability for translation in human cells, leading to differences in protein levels and eventually eliciting phenotypic differences.

AUTHOR SUMMARYThe genetic code is redundant, with several codons encoding for the same amino acid. These synonymous codons are not used with equal frequencies. Instead, codon usage bias (CUB) varies between species, genes, and even positions along a gene. At each of these levels, CUB are shaped by the overall balance between mutational biases and selection forces. To shed light on the molecular mechanisms underlying molecular and organism phenotypes, integrative studies quantifying the phenotypic impact of CUB at different levels of biological integration are necessary. Here, we monitored the multiscale changes induced by six synonymous versions of an antibiotic resistance gene independently expressed in a human cell line. We show that: 1. both mRNA levels, protein levels, and their ratios are affected by CUB; 2. potential effects on mRNA (splicing and availability for translation) seem to constitute the main level of action; 3. cell fitness is severely impacted by a trade-off between the burden of heterologous expression and the benefits of antibiotic resistance. This integrative study provides new insights on translation regulation and the associated phenotypic impact in human cells, associated to CU.

category: systems biology


title: The multi-level phenotypic impact of synonymous substitutions: heterologous gene expression in human cells.


title: Diversity of cave Phlebotomines (Diptera: Psychodidae) from a Colombian cave

abstract: Sandflies are vector species of Leishmania, among many other pathogens, with a global distribution and a variety of ecological niches. Previous samplings have found that karstic formations (i.e., caves and folds formed by the erosion of limestone) serve as a natural habitat to sandfly species. The majority of samplings of cave sandfly diversity have occurred in Brazil and to date none have studied the species composition in a cave in the Northern Andes. We collected sandflies in the Cave "Los Guacharos", in the state of Antioquia, Colombia. The sampling was carried out during two consecutive nights in September 2019. CDC-type light traps were installed inside the cavern and in other surrounding karst systems (caves and folds). In total, we identified 18 species of sandfly from the cave and surrounding karst systems, including three new records for Colombia (Bichromomyia olmeca nociva, Brumptomyia brumpti, and Warileya leponti), and provide the first karstic reports for four other species (Lutzomyia gomezi, Lutzomyia hartmanni, Pintomyia ovallesi, and Psychodopygus panamensis). We then used the results of our survey and published literature to test two hypotheses. First, that sandfly diversity in Neotropical caves is richest nearer to the equator and second that there is a phylogenetic signal of karstic habitat use in sandflies. Counter to our predictions, we found no evidence that diversity follows a latitudinal gradient. Further, we find no evidence of a phylogenetic signal of karstic habitat use, instead finding that the use of caves likely evolved multiple times across several genera. Our results highlight the importance of a wide sampling to understand the natural habitat of sandflies and other disease vectors.

category: zoology


title: Strong Aversive Conditioning Triggers a Long-Lasting Generalized Aversion

abstract: Generalization is an adaptive mnemonic process in which an animal can leverage past learning experiences to navigate future scenarios, but overgeneralization is a hallmark feature of anxiety disorders. Therefore, understanding the synaptic plasticity mechanisms that govern memory generalization and its persistence is an important goal. Here, we demonstrate that strong CTA conditioning results in a long-lasting generalized aversion that persists for at least two weeks. Using brain slice electrophysiology and activity-dependent labeling of the conditioning-active neuronal ensemble within the gustatory cortex, we find that strong CTA conditioning induces a long-lasting increase in synaptic strengths that occurs uniformly across superficial and deep layers of GC. Repeated exposure to salt, the generalized tastant, causes a rapid attenuation of the generalized aversion that correlates with a reversal of the CTA-induced increases in synaptic strength. Unlike the uniform strengthening that happens across layers, reversal of the generalized aversion results in a more pronounced depression of synaptic strengths in superficial layers. Finally, the generalized aversion and its reversal do not impact the acquisition and maintenance of the aversion to the conditioned tastant (saccharin). The strong correlation between the generalized aversion and synaptic strengthening, and the reversal of both in superficial layers by repeated salt exposure, strongly suggests that the synaptic changes in superficial layers contribute to the formation and reversal of the generalized aversion. In contrast, the persistence of synaptic strengthening in deep layers correlates with the persistence of CTA. Taken together, our data suggest that layer-specific synaptic plasticity mechanisms separately govern the persistence and generalization of CTA memory.


title: Two-photon patterned photostimulation with low-power, high-efficiency and reliable single-cell optogenetic control

abstract: Two-photon optogenetics enables selectively stimulating individual cells for manipulating neuronal ensembles. As the general photostimulation strategy, the patterned two-photon excitation has enabled millisecond-timescale activation for single or multiple neurons, but its activation efficiency is suffered from high laser power due to low beam-modulation efficiency. Here, we develop a high- efficiency beam-shaping method based on the Gerchberg-Saxton (GS) algorithm with spherical-distribution initial phase (GSSIP) to reduce the patterned two-photon excitation speckles and intensity. It can well control the phase of shaped beams to attain speckle-free accurate patterned illumination with an improvement of 44.21% in the modulation efficiency compared with that of the traditional GS algorithm. A combination of temporal focusing and the GSSIP algorithm (TF-GSSIP) achieves patterned focusing through 500-m-thickness mouse brain slices, which is 2.5 times deeper than the penetration depth of TF-GS with the same signal-to-noise ratio (SNR). With our method, the laser power can be reduced to only 55.56% of that with traditional method (the temporal focusing with GS, TF-GS) to reliably evoke GCaMP6s response in C1V1-expressing cultured neurons with single-cell resolution. Besides, the photostimulation efficiency is remarkably increased by 80.19% at the same excitation density of 0.27 mW/m2. This two-photon stimulation method with low-power, reliable and patterned illumination may pave the way for analyzing neural circuits and neural coding and decoding mechanism.


title: Sex differences in the alcohol-mediated modulation of BLA network states

abstract:
About 85% of adults in the United States report drinking alcohol in their lifetime. Mood disorders, like generalized anxiety disorder and major depression, are highly comorbid with alcohol use. The basolateral amygdala (BLA) is an area of the brain that is heavily implicated in both mood disorders and alcohol use disorder. Importantly, modulation of BLA network/oscillatory states via parvalbumin-positive (PV) GABAergic interneurons has been shown to control the behavioral expression of fear and anxiety. Further, PV interneurons express a high density of {delta}-subunit-containing GABAA receptors (GABAARs), which are sensitive to low concentrations of alcohol. Our lab previously demonstrated that {delta}-subunit-containing GABAARs on PV interneurons in the BLA influence voluntary ethanol intake and anxiety-like behavior in withdrawal. Therefore, we hypothesized that the effects of alcohol may modulate BLA network states that have been associated with fear and anxiety behaviors via {delta}-GABAARs on PV interneurons in the BLA. Given the impact of ovarian hormones on the expression of {delta}-GABAARs, we examined the ability of alcohol to modulate local field potentials (LFPs) in the BLA from male and female C57BL/6J and Gabrd-/- mice after acute and repeated exposure to alcohol. Here, we demonstrate that acute and repeated alcohol can differentially modulate oscillatory states in male and female C57BL/6J mice, a process which involves {delta}-GABAARs. This is the first study to demonstrate that alcohol is capable of altering network states implicated in both anxiety and alcohol use disorders.

Significance StatementAlcohol use disorder and mood disorders are highly comorbid. The basolateral amygdala (BLA) is implicated in both of these disorders, but the mechanisms contributing to their shared pathophysiology remain uncertain. Here we demonstrate that acute and repeated alcohol exposure can alter network oscillations in the BLA which control the behavioral expression of fear and anxiety. These data suggest that alcohol may directly influence network states associated with mood. Further, we demonstrate sex differences in alcohols ability to modulate BLA network states, an effect involving {delta}-GABAA receptors, which may contribute to sex differences in alcohol intake and comorbid mood disorders. These data potentially point to a novel mechanism mediating the effects of alcohol on affective states.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY


title: Anterior cingulate metabolite levels, memory, and inhibitory control in abstinent men and women with Alcohol Use Disorder

abstract:
AimsAlcohol use disorder (AUD), has been shown to have harmful cognitive and physiological effects, including altered brain chemistry. Further, although men and women may differ in vulnerability to the neurobiological effects of AUD, results of existing studies have been conflicting. Brain metabolite levels and cognitive functions were examined in a cross section of men with AUD (AUDm) and women with AUD (AUDw) to determine degree of abnormalities after extended periods of abstinence (mean, six years), and to evaluate gender differences in cognitive and metabolite measures.

MethodsParticipants were 40 abstinent individuals with AUD (22 AUDw, 18 AUDm) and 50 age-equivalent non-AUD comparison participants (26 NCw, 24 NCm). Proton magnetic resonance spectroscopy (MRS) was employed at 3 Tesla to acquire metabolite spectra from the dorsal anterior cingulate cortex (dACC). Brain metabolites N-acetylaspartate (NAA), choline (Cho), myo-Inositol (mI), and glutamate & glutamine (Glx) were examined relative to measures of memory and inhibitory control.

ResultsMetabolite levels in the AUD group showed no significant differences from the NC group. Memory and inhibitory-control impairments were observed in the AUD group. There also were significant group-specific associations between metabolite ratios and measures of inhibitory control. There were no Group-by-Gender interactions for the four metabolite ratios.

ConclusionsThese findings demonstrate that brain metabolite levels in men and women with AUD, following long-term abstinence, do not differ from individuals without AUD. The data also provide evidence of associations between metabolite levels and measures of inhibitory control, a functional domain important for curtailing harmful drinking.


title: Molecular Structure Study on the Polyelectrolyte Properties of Actin Filaments

abstract: An accurate characterization of the polyelectrolyte properties of actin filaments might provide a deeper understanding of the fundamental mechanisms governing the intracellular ionic wave packet propagation in neurons. Infinitely long cylindrical models for actin filaments and approximate electrochemical theories for the electrolyte solutions were recently used to characterize these properties in in vitro and intracellular conditions. This article uses a molecular structure model for actin filaments to investigate the impact of roughness and finite size on the mean electrical potential, ionic density distributions, currents, and conductivities. We solved the electrochemical theories numerically without further approximations. Our findings bring new insights into the electrochemical interactions between a filaments irregular surface charge density and the surrounding medium. The irregular shape of the filament structure model generated pockets, or hot spots, where the current density reached higher or lower magnitudes than those in neighboring areas throughout the filament surface. It also revealed the formation of a well-defined asymmetric electrical double layer with a thickness larger than that commonly used for symmetric models.


title: Nrf2 attenuates the innate immune response after experimental myocardial infarction

abstract:
ObjectivesWe aimed to investigate the contribution of the transcription factor nuclear factor erythroid-derived 2-like 2 (Nrf2) to the inflammatory response after experimental myocardial infarction (MI).

BackgroundThere is compelling evidence implicating dysregulated inflammation in the mechanism of ventricular remodeling and heart failure (HF) after MI. The transcription factor Nrf2 (encoded by Nfe2l2) is a promising target in this context. It impedes transcriptional upregulation of pro-inflammatory cytokines and is anti-inflammatory in various murine models.

MethodsWe subjected Nrf2-/- mice and wild type (WT) controls to permanent left coronary artery (LCA) ligation. The inflammatory response was investigated with fluorescence-activated cell sorting (FACS) analysis of peripheral blood and heart cell suspensions, together with qRT-PCR of infarcted tissue for chemokines and their receptors. To investigate whether Nrf2-mediated transcription is a dedicated function of leukocytes, we interrogated publicly available RNA-sequencing (RNA-seq) data from mouse hearts after permanent LCA ligation for Nrf2-regulated gene (NRG) expression.

ResultsFACS analysis demonstrated a profoundly inflamed phenotype in the hearts of global Nrf2-/- mice as compared to WT mice after MI. Moreover, infarcted tissue from Nrf2-/- mice displayed higher expression of inflammatory cytokines, chemokines, and their receptors, including IL6, Ccl2, and Cxcr4. RNA-seq analysis showed upregulated NRG expression in WT mice after MI compared to untreated mice, which was significantly higher in bioinformatically isolated CCR2+ cells.

ConclusionsTaken together, the results suggest that Nrf2 signalling in leukocytes, and possibly CCR2+ monocyte-derived cardiac resident macrophages, may be potential targets to prevent post-MI ventricular remodeling.

category: immunology


title: RelCoVax(R), a two antigen subunit protein vaccine candidate against SARS-CoV-2 induces strong immune responses in mice

abstract: The COVID-19 pandemic has spurred an unprecedented movement to develop safe and effective vaccines against the SARS-CoV-2 virus to immunize the global population. The first set of vaccine candidates that received emergency use authorization targeted the spike (S) glycoprotein of the SARS-CoV-2 virus that enables virus entry into cells via the receptor binding domain (RBD). Recently, multiple variants of SARS-CoV-2 have emerged with mutations in S protein and the ability to evade neutralizing antibodies in vaccinated individuals. We have developed a dual RBD and nucleocapsid (N) subunit protein vaccine candidate named RelCoVax(R) through heterologous expression in mammalian cells (RBD) and E. coli (N). The RelCoVax(R) formulation containing a combination of aluminum hydroxide (alum) and a synthetic CpG oligonucleotide as adjuvants elicited high antibody titers against RBD and N proteins in mice after a prime and boost dose regimen administered 2 weeks apart. The vaccine also stimulated cellular immune responses with a potential Th1 bias as evidenced by increased IFN-{gamma} release by splenocytes from immunized mice upon antigen exposure particularly N protein. Finally, the serum of mice immunized with RelCoVax(R) demonstrated the ability to neutralize two different SARS-CoV-2 viral strains in vitro including the Delta strain that has become dominant in many regions of the world and can evade vaccine induced neutralizing antibodies. These results warrant further evaluation of RelCoVax(R) through advanced studies and contribute towards enhancing our understanding of multicomponent subunit vaccine candidates against SARS-CoV-2.


title: Membrane interactions and uncoating of Aichi virus, a picornavirus that lacks a VP4

abstract:
Kobuviruses are an unusual and poorly characterised genus within the picornavirus family, and can cause gastrointestinal enteric disease in humans, livestock and pets. The human Kobuvirus, Aichi virus (AiV) can cause severe gastroenteritis and deaths in children below the age of five years, however this is a very rare occurrence. During the assembly of most picornaviruses (e.g. poliovirus, rhinovirus and foot-and-mouth disease virus), the capsid precursor protein VP0 is cleaved into VP4 and VP2. However, Kobuviruses retain an uncleaved VP0. From studies with other picornaviruses, it is known that VP4 performs the essential function of pore formation in membranes, which facilitates transfer of the viral genome across the endosomal membrane and into the cytoplasm for replication. Here, we employ genome exposure and membrane interaction assays to demonstrate that pH plays a critical role in AiV uncoating and membrane interactions. We demonstrate that incubation at low pH alters the exposure of hydrophobic residues within the capsid, enhances genome exposure and enhances permeabilisation of model membranes. Furthermore, using peptides we demonstrate that the N-terminus of VP0 mediates membrane pore formation in model membranes, indicating that this plays an analogous function to VP4.

ImportanceTo initiate infection, viruses must enter a host cell and deliver their genome into the appropriate location. The picornavirus family of small non-enveloped RNA viruses includes significant human and animal pathogens and are also models to understand the process of cell entry. Most picornavirus capsids contain the internal protein VP4, generated from cleavage of a VP0 precursor. During entry, VP4 is released from the capsid. In enteroviruses this forms a membrane pore, which facilitates genome release into the cytoplasm. Due to high levels of sequence similarity, it is expected to play the same role for other picornaviruses. Some picornaviruses, such as Aichi virus, retain an intact VP0, and it is unknown how these viruses re-arrange their capsids and induce membrane permeability in the absence of VP4. Here we have used Aichi virus as a model VP0 virus to test for conservation of function between VP0 and VP4. This could enhance understanding of pore function and lead to development of novel therapeutic agents that block entry.


title: Live imaging of Yersinia translocon formation and immune recognition in host cells

abstract: Yersinia enterocolitica employs a type three secretion system (T3SS) to translocate immunosuppressive effector proteins into host cells. To this end, the T3SS assembles a translocon/pore complex composed of the translocator proteins YopB and YopD in host cell membranes serving as an entry port for the effectors. The translocon is formed in a Yersinia-containing pre-phagosomal compartment that is connected to the extracellular space. As the phagosome matures, the translocon and the membrane damage it causes are recognized by the cell-autonomous immune system. We infected cells in the presence of fluorophore-labeled ALFA-tag-binding nanobodies with a Y. enterocolitica strain expressing YopD labeled with an ALFA-tag. Thereby we could record the integration of YopD into translocons and its intracellular fate in living host cells. YopD was integrated into translocons around 2 min after uptake of the bacteria into a phosphatidylinositol-4,5-bisphosphate enriched pre-phagosomal compartment and remained there for 27 min on average. Damaging of the phagosomal membrane as visualized with recruitment of GFP-tagged galectin-3 occurred in the mean around 14 min after translocon formation. Shortly after recruitment of galectin-3, guanylate-binding protein 1 (GBP-1) was recruited to phagosomes, which was accompanied by a decrease in the signal intensity of translocons, suggesting their degradation. In sum, we were able for the first time to film the spatiotemporal dynamics of Yersinia T3SS translocon formation and degradation and its sensing by components of the cell-autonomous immune system.


title: Pilot scale study on UV-C inactivation of bacterial endospores and virus particles in whole milk: evaluation of system efficiency and product quality

abstract: UV-C processing of whole milk (WM) using a designed pilot scale Dean flow system was conducted at flow rates (11.88, 23.77, and 47.55 gph), Reynolds number ranges from 2890-11562 and the Dean number (at curved region) calculated as (648-2595) to inactivate bacterial endospores and virus particles. Biodosimetry studies were conducted to quantify the reduction equivalent fluence at selected experimental conditions. Results revealed that the fluence distribution improved as flow rate increases, attributed to increase in Dean effects and turbulence intensity. Microbial inactivation studies conducted at 47.55 gph showed 0.91 {+/-} 0.15 and 2.14 {+/-} 0.19 log reduction/ pass for B. cereus endospores and T1UV phage. Linear inactivation trend was observed against number of passes which clearly demonstrates equivalent fluence delivery during each pass. Lipid peroxidation value and volatiles profile did not change significantly at UV fluence of 60 mJ/cm2. Lower EEO value signifies the higher electrical efficiency of the system.


title: Visualization of SpoVAEa protein dynamics in dormant spores of Bacillus cereus and dynamic changes in their germinosomes and SpoVAEa during germination

abstract:
Bacillus cereus spores, like most Bacillus spores, can survive for years, germinate when their surroundings become suitable, and spore germination proteins play an important role in the initiation of germination. Because germinated spores lose dormant spores extreme resistance, information on the function of germination proteins could be useful in developing new strategies to control B. cereus spores. Prior work has shown that: i) the channel protein SpoVAEa exhibits high frequency movement in the outer leaflet of the inner membrane (IM) in dormant spores of B. subtilis; ii) the formation of the foci termed the germinosome between two germination proteins, the germinant receptor GerR and the scaffold protein GerD, in developing spores of B. cereus is slower than foci formation by GerR and GerD individually. However, the dynamics of movement of SpoVAEa in B. cereus spores, and the behaviour of the germinosome in B. cereus spore germination are unclear. In this study, we found that SpoVAEa fluorescent foci in dormant spores of B. cereus move on the IM, but slower than in B. subtilis spores, and likely colocalize transiently with GerD-mScarlet-I in the germinosome. Our results further indicate that: i) expression of GerR-SGFP2 and SpoVAEa-SGFP2 with GerD-mScarlet-I from a plasmid leads to more heterogeneity and lower efficiency of spore germination in B. cereus; and ii) germinosome foci observed by Fluorescence Resonance Energy Transfer (FRET) between GerR-SGFP2 and GerD-mScarlet-I b can be lost soon after the spore phase transition. However this is not always the case, as some GerR-SGFP2 foci and GerD-mScarlet-I foci continued to exist, colocalize, and even show a weak FRET signal. These data highlight the heterogeneous behaviour of spore germination protein complexes and indicate that some complexes may persist well beyond the initiation of germination.

IMPORTANCEBacillus cereus is commonly present in soil and does harm to humans via contaminated food. In this study, we used B. cereus spores to investigate the movement of the spore-specific channel protein SpoVAEa, the interaction between SpoVAEa and the germinosome scaffold protein GerD, as well as the dynamics of a number of germination proteins in spore germination. Our results expand upon observations of the interactions between specific B. cereus spore germination proteins, in particular the GerR germinant receptor A, B and C subunits and GerD, as well as between SpoVAEa and GerD. The approaches used in this work could also be used to examine the interactions between GerD and SpoVAEa and other germination proteins in spores of other Bacillus species.


title: The Salmonella transmembrane effector SteD hijacks AP1-mediated vesicular trafficking for delivery to antigen-loading MHCII compartments

abstract: SteD is a transmembrane effector of the Salmonella SPI-2 type III secretion system that inhibits T cell activation by reducing the amounts of at least three proteins - major histocompatibility complex II (MHCII), CD86 and CD97 - from the surface of antigen-presenting cells. SteD specifically localises at the trans-Golgi network (TGN) and MHCII compartments; however, the targeting, membrane integration and trafficking of SteD are not understood. Using systematic mutagenesis, we identify distinct regions of SteD that are required for these processes. We show that SteD integrates into membranes of the ER/Golgi through a two-step mechanism of membrane recruitment from the cytoplasm followed by integration. SteD then migrates to and accumulates within the TGN. From here it hijacks the host adaptor protein (AP)1-mediated trafficking pathway from the TGN to MHCII compartments. AP1 binding and post-TGN trafficking require a short sequence in the N-terminal cytoplasmic tail of SteD that resembles the AP1-interacting dileucine sorting signal, but in inverted orientation, suggesting convergent evolution.


title: Establishment of a transparent soil system to study Bacillus subtilis chemical ecology

abstract: Bacterial secondary metabolites are structurally diverse molecules that drive microbial interaction by altering growth, cell differentiation, and signaling. Bacillus subtilis, a Gram-positive soil-dwelling bacterium, produces a wealth of secondary metabolites, among them, lipopeptides have been vastly studied by their antimicrobial, antitumor, and surfactant activities. However, the natural functions of secondary metabolites in the lifestyles of the producing organism remain less explored under natural conditions, i.e. in soil. Here, we describe a hydrogel-based transparent soil system to investigate B. subtilis chemical ecology under controllable soil-like conditions. The transparent soil matrix allows the growth of B. subtilis and other isolates gnotobiotically and under nutrient-controlled conditions. Additionally, we show that transparent soil allows the detection of lipopeptides production and dynamics by HPLC-MS and MALDI-MS imaging, along with fluorescence imaging of 3-dimensional bacterial assemblages. We anticipate that this affordable and highly controllable system will promote bacterial chemical ecology research and help to elucidate microbial interactions driven by secondary metabolites.


title: Differential response of Candida species morphologies and isolates to fluconazole and boric acid

abstract: Candida albicans is the most prevalent cause of vulvovaginal candidiasis ( yeast infection) and recurrent vulvovaginal candidiasis, though the incidence of non-albicans yeast species is increasing. The azole fluconazole is the primary antifungal drug used to treat R/VVC yet isolates from some species have intrinsic resistance to fluconazole, and recurrent infection can occur even with fluconazole-susceptible populations. The second-line broad-spectrum antimicrobial drug, boric acid, is an alternative treatment that has been found to successfully treat complicated VVC infections. Far less is known about how boric acid inhibits growth of yeast isolates in different morphologies compared to fluconazole. We found significant differences in drug resistance and drug tolerance (the ability of a subpopulation to grow slowly in high levels of drug) between C. albicans, C. glabrata, and C. parapsilosis isolates, with the specific relationships dependent on both drug and phenotype. Population-level variation for both susceptibility and tolerance was broader for fluconazole than boric acid in all species. Unlike fluconazole, which neither prevented hyphal formation nor disrupted mature biofilms, boric acid inhibited C. albicans hyphal formation and reduced mature biofilm biomass and metabolic activity in all isolates in a dose-dependent manner. Variation in planktonic response did not generally predict biofilm phenotypes. Overall, our findings illustrate that boric acid is broadly effective at inhibiting growth across many isolates and morphologies, which could explain why it is an effective treatment for R/VVC.


title: Molnupiravir combined with different repurposed drugs further inhibits SARS-CoV-2 infection in human nasal epithelium in vitro

abstract: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in late 2019, has caused a worldwide pandemic with unprecedented economic and societal impact. Currently, several vaccines are available, and multitudes of antiviral treatments have been proposed and tested. Although many of the vaccines show high clinical efficacy, they are not equally accessible worldwide. Additionally, due to the continuous emergence of new virus variants, and generally short duration of immunity, the development of safe and effective antiviral treatments remains of the utmost importance. Since the emergence of SARS-CoV-2, substantial efforts have been undertaken to repurpose existing and approved drugs for accelerated clinical testing and potential emergency use authorizations. However, drug-repurposing using high throughput screenings in cellular assays, often identify hits that later prove ineffective in clinical studies. Our approach was to evaluate the activity of compounds that have either been tested clinically or already undergone extensive preclinical profiling, using a standardized in vitro model of human nasal epithelium. Secondly, we evaluated drug combinations using sub-maximal doses of each active single compound. Here, we report the antiviral effects of 95 single compounds and 30 combinations. The data show that selected drug combinations including 10 M of molnupiravir, a viral RNA-dependent RNA polymerase (RdRp) inhibitor, effectively inhibit SARS-CoV-2 replication. This indicates that such combinations are worthy of further evaluation as potential treatment strategies against coronavirus disease 2019 (COVID-19).


title: Fusogenicity and neutralization sensitivity of the SARS-CoV-2 Delta sublineage AY.4.2

abstract: SARS-CoV-2 lineages are continuously evolving. As of December 2021, the AY.4.2 Delta sub-lineage represented 20 % of sequenced strains in UK and has been detected in dozens of countries. It has since then been supplanted by the Omicron variant. AY.4.2 displays three additional mutations (T95I, Y145H and A222V) in the N-terminal domain (NTD) of the spike when compared to the original Delta variant (B.1.617.2) and remains poorly characterized. Here, we analyzed the fusogenicity of the AY.4.2 spike and the sensitivity of an authentic AY.4.2 isolate to neutralizing antibodies. The AY.4.2 spike exhibited similar fusogenicity and binding to ACE2 than Delta. The sensitivity of infectious AY.4.2 to a panel of monoclonal neutralizing antibodies was similar to Delta, except for the anti-RBD Imdevimab, which showed incomplete neutralization. Sensitivity of AY.4.2 to sera from individuals having received two or three doses of Pfizer or two doses of AstraZeneca vaccines was reduced by 1.7 to 2.1 fold, when compared to Delta. Our results suggest that mutations in the NTD remotely impair the efficacy of anti-RBD antibodies. The temporary spread of AY.4.2 was not associated with major changes in spike function but rather to a partially reduced neutralization sensitivity.