Cell Rearrangement Generates Pattern Emergence as a Function of Temporal Morphogen Exposure. As tissues elongate, cell rearrangement alters positional information in manner that must be accounted for to generate gene expression pattern. How this is achieved during paraxial mesoderm elongation is unknown. By reverse-engineering gene regulatory networks that predict single cell expression trajectories across the tissue, we find a network capable of recapitulating the full range of dynamic differentiation profiles observed both in vivo and in vitro. Simulating gene expression profiles on in toto cell tracking data sets reveal that temporal exposure to Wnt and FGF is generated by cell movement. The absence of reversal in gene expression towards a more premature gene expression state predicts the generation of aberrant tbx6 expression in the posterior progenitor zone that we then confirm by quantitative single cell imaging. Taken together, these results demonstrate cell rearrangement tunes the dynamics of mesoderm progenitor differentiation to generate pattern emergence as a function of temporal Wnt and FGF exposure.
R7 photoreceptor axon targeting depends on the relative levels of lost and found expression in R7 and its synaptic partners As neural circuits form, growing processes select the correct synaptic partners through interactions between cell surface proteins. The presence of such proteins on two neuronal processes may lead to either adhesion or repulsion; however, the consequences of mismatched expression have rarely been explored. Here we show that the Drosophila CUB-LDL protein Lost and found (Loaf) is required in the UV-sensitive R7 photoreceptor for normal axon targeting only when Loaf is also present in its synaptic partners. Although targeting occurs normally in loaf mutant animals, removing loaf from photoreceptors or expressing it in their postsynaptic neurons Tm5a/b or Dm9 in a loaf mutant causes mistargeting of R7 axons. Loaf localizes primarily to intracellular vesicles including endosomes. We propose that Loaf regulates the trafficking or function of one or more cell surface proteins, and an excess of these proteins on the synaptic partners of R7 prevents the formation of stable connections.
Alien chromatin other than the GST-encoding Fhb7 candidate confers Fusarium head blight resistance in wheat breeding The lack of resistance resources is a major bottleneck for wheat Fusarium head blight (FHB) resistance breeding. Three wheat-Th. elongatum FHB resistant translocation lines have been developed and used for wheat breeding without yield penalty. Transcriptomic analysis identified a derivative glutathione S-transferase transcript T26102, which was homologous to Fhb7 and induced dramatically by Fusarium graminearum. Unlike other studies, Fhb7 homologs were detected not only in Thinopyrum but also in Elymus, Leymus, Pseudoroegeria and Roegeria. We also found that several wheat-Th. ponticum derivatives carrying Fhb7 and its homologs were highly susceptible to FHB. Moreover, the transgenic plants expressing Fhb7 and its homolog on different backgrounds did not improve the FHB resistance.

One Sentence SummaryThe GST-encoding Fhb7 candidate cannot improve Fusarium head blight resistance in wheat breeding. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS. Rescue of stalled clathrin-mediated endocytosis by asymmetric Arp2/3-mediated actin assembly Actin assembly facilitates vesicle formation in several trafficking pathways. Clathrin-mediated endocytosis (CME) shows elevated actin assembly dependence under high membrane tension. Why actin assembly at CME sites occurs heterogeneously even within the same cell, and how assembly forces are harnessed, are not fully understood. Here, endocytic dynamics, actin presence, and geometry of CME proteins from three different functional modules, were analyzed using three-dimensional (3D) super-resolution microscopy, live-cell imaging, and machine-learning-based computation. When hundreds of CME events were compared, sites with actin assembly showed a distinct signature, a delay between completion of coat expansion and vesicle scission, indicating that actin assembly occurs preferentially at stalled CME sites. N-WASP is recruited to one side of CME sites where it is positioned to stimulate asymmetric actin assembly. We propose that asymmetric actin assembly rescues stalled CME sites by pulling vesicles into the cell much like a bottle opener pulls off a bottle cap. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Multi-omics-based label-free metabolic flux inference reveals obesity-associated dysregulatory mechanisms in liver glucose metabolism Glucose homeostasis is maintained by modulation of metabolic flux. Enzymes and metabolites regulate the involved metabolic pathways. Dysregulation of glucose homeostasis is a pathological event in obesity. Analyzing metabolic pathways and the mechanisms contributing to obesity-associated dysregulation in vivo is challenging. Here, we introduce OMELET: Omics-Based Metabolic Flux Estimation without Labeling for Extended Trans-omic Analysis. OMELET uses metabolomic, proteomic, and transcriptomic data to identify changes in metabolic flux, and to quantify contributions of metabolites, enzymes, and transcripts to the changes in metabolic flux. By evaluating the livers of fasting ob/ob mice, we found that increased metabolic flux through gluconeogenesis resulted primarily from increased transcripts, whereas that through the pyruvate cycle resulted from both increased transcripts and changes in substrates of metabolic enzymes. With OMELET, we identified mechanisms underlying the obesity-associated dysregulation of metabolic flux in liver.

HighlightsO_LIWe created OMELET to infer metabolic flux and its regulation from multi-omic data.
C_LIO_LIGluconeogenic and pyruvate cycle fluxes increased in fasting ob/ob mice.
C_LIO_LITranscripts increases mediated the increase in gluconeogenic fluxes in ob/ob mice.
C_LIO_LIIncreases in transcripts and substrates enhanced pyruvate cycle flux in ob/ob mice.
C_LI Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. Temporally robust occupancy frequency distributions in riverine ecosystems explained by local biodiversity regulation The distribution of species site occupancy--the occupancy frequency distribution (OFD)--is a fundamental pattern in spatial ecology. Despite decades of research, the mechanisms responsible for the shape of the distribution remain incompletely understood. Here, we simultaneously examine both spatial and temporal patterns of site occupancy in communities of macroinverterbrates, macrophytes and diatoms. We find that the shape of the OFD at the catchment scale is typically time invariant despite substantial turnover in composition and species proportional occupancy over the 30 year time series. This can been seen as evidence in favour of an ecological mechanism for the shape of the OFD driven by local regulation of species richness and other biodiversity properties. We show that a simple, locally saturated patch occupancy model explains both the shape of the OFD and the temporal turnover rate in English rivers. This model is mathematically analogous to Hubbells neutral theory however concerns patches occupied by a given species rather its population size. Our analysis supports the view that metacommunities exist in a dynamic steady state arising from local ecological constraints and that this intrinsic regulation can drives the emergence of a wide variety of macroecological patterns. We also show that interspecific differences are unlikely to play a dominant role in driving large scale ecosystem structure. T606-phosphorylation deprives the function of Kaiso as a transcription and oncogenic factor It is well known that Kaiso protein encoded by ZBTB33 gene is a transcription repressor and that Kaiso-P120ctn interaction increases the shift of Kaiso from the nucleus into the cytoplasm. However, the regulatory mechanisms of Kaiso compartmentalization are far from clear. Here, we reported that AKT1 could phosphorylate 606-threonine residue (T606) within the RSSTIP motif of Kaiso in the cytoplasm. The T606-phosphorylated Kaiso (pT606-Kaiso) could directly bind to 14-3-3 family proteins and the depletion of T606 phosphorylation by T606A mutation abolished most of the Kaiso-14-3-3 binding. In addition, the Kaiso-P120ctn interaction was essential for the pT606-Kaiso accumulation in the cytoplasm. Notably, enforced 14-3-3{sigma} (SFN) overexpression could increase the pT606-Kaiso accumulation in the cytoplasm and de-repress the transcription of Kaiso target gene CDH1. Decreased amounts of both pT606-Kaiso and CDH1 proteins were frequently observed in human gastric cancer tissues relative to paired normal controls. The mRNA levels of 14-3-3{sigma} and Kaiso target gene CDH1 were positively and significantly correlated with each other in bioinformatics analyses using publicly available RNA-seq datasets for human normal tissues (n=11688, r=0.60, p<0.001) in the GTEx project and for cancer cell lines (n=1156, r=0.41, p<0.001) in the CCLE project. Furthermore, Kaiso T606A mutant (unable to be phosphorylated) significantly increased the migration and invasion of cancer cells in vitro as well as boosted the growth of these cells in vivo. In conclusion, Kaiso could be phosphorylated by AKT1 at the T606 and the pT606-Kaiso accumulates in the cytoplasm through binding to 14-3-3/P120ctn that de-represses the expression of Kaiso target gene CDH1 in normal tissues. Decreased Kaiso phosphorylation may contribute to the development of gastrointestinal cancer. The status of Kaiso phosphorylation is a determinant factor for the role of Kaiso in the development of cancer. Chromatin accessibility landscape and active transcription factors in primary human invasive lobular and ductal breast carcinomas Background: Invasive lobular breast carcinoma (ILC), the second most prevalent histological subtype of breast cancer, exhibits unique molecular features compared with the more common invasive ductal carcinoma (IDC). While genomic and transcriptomic features of ILC and IDC have been characterized, genome-wide chromatin accessibility pattern differences between ILC and IDC remain largely unexplored. Methods: Here, we characterized tumor-intrinsic chromatin accessibility differences between ILC and IDC using primary tumors from The Cancer Genome Atlas (TCGA) breast cancer assay for transposase-accessible chromatin with sequencing (ATAC-seq) dataset. Results: We identified distinct patterns of genome-wide chromatin accessibility in ILC and IDC. Inferred patient-specific transcription factor (TF) motif activities revealed regulatory differences between and within ILC and IDC tumors. EGR1, RUNX3, TP63, STAT6, SOX family, and TEAD family TFs were higher in ILC, while ATF4, PBX3, SPDEF, PITX family, and FOX family TFs were higher in IDC. Conclusions: This study reveals the distinct epigenomic features of ILC and IDC and the active TFs driving cancer progression that may provide valuable information on patient prognosis. Neuronal Activation of the Gastrointestinal Tract Shapes the Gut Environment in Mice The enteric nervous system (ENS) coordinates responses that shape the intestinal milieu, and participates in bi-directional communication with various organs including the brain. Tools to study, and knowledge of, ENS development and function lag behind brain research. Herein, we deploy recombinant adeno-associated viral (rAAV) vectors with enhanced tropism for the gut to map and activate gut associated neurons in mice with spatial and temporal resolution. rAAV-mediated fluorescent labelling coupled with whole-tissue clearing methods in the small intestine and colon reveal novel neuronal architectures in the gut of mice. Viral delivery of genetically-encoded calcium indicators enable in vivo measurement of ENS activity. Further, we employed chemogentics to specifically activate gut neurons that express choline acetyltransferase (ChAT+) or tyrosine hydroxylase (TH+). Targeted activation of ChAT+ or TH+ neuronal populations associated with the gastrointestinal (GI) tract alters the intestinal transcriptome, mouse and microbial proteomes, metagenome, and fecal metabolome. We reveal previously unknown and broad roles for enteric neurons in modulating GI physiology, mucosal immunity, and microbiome structure, and provide evidence for novel functions by the ENS such as regulating fungal colonization and shaping of bile acid profiles in the gut. These advanced experimental tools and rich datasets with multi-parameter characterization of the gut ecosystem may facilitate expanded understanding and study of the ENS. Targeting CISH enhances natural cytotoxicity receptor signaling and reduces NK cell exhaustion to improve solid tumor immunity Cytokine inducible SH2-containing protein (CISH) is a natural killer (NK) cell negative regulator of cytokine signaling pathway. To further understand CISH functions in NK cells, we developed a conditional Cish-deficient mouse model in NK cells (Cishfl/flNcr1Ki/+). We detected no developmental or homeostatic difference in NK cells.

However, global gene expression of Cishfl/flNcr1Ki/+ NK cells compared to Cish+/+Ncr1Ki/+ NK cells revealed upregulation of pathways and genes associated with NK cell cycling and activation. We show that CISH does not only regulate interleukin-15 (IL-15) signaling pathways but also natural cytotoxicity receptors (NCR) pathways. Indeed, CISH protein expression level increases upon NCR triggering. Primed Cishfl/flNcr1Ki/+ NK cells display increased activation upon NCR stimulation. Cishfl/flNcr1Ki/+ NK cells display lower activation thresholds and Cishfl/flNcr1Ki/+ mice are more resistant to tumor metastasis. Remarkably, we found that Cishfl/flNcr1Ki/+ mice were also more resistant to primary breast cancer growth in addition to superior control of spontaneous tumor metastasis. CISH deletion favors NK cell accumulation to the primary tumor, optimizes NK cell killing properties and decreases TIGIT immune checkpoint receptor expression, limiting NK cell exhaustion. Finally, we argue that specifically enhancing NK cell function is sufficient to boost anti-tumor response to both primary and secondary tumor models. Using CRISPRi, we then targeted CISH in human NK-92 or primary NK cells. According to the results in our mouse model, CISH deletion favors NCR signaling and anti-tumor functions in human NK cells. Our results validate CISH as an emerging therapeutic target to enhance NK cell immunotherapy. TAP: Targeting and analysis pipeline for optimization and verification of coil placement in transcranial magnetic stimulation ObjectiveTranscranial magnetic stimulation (TMS) can modulate brain function via an electric field (E-field) induced in a region of interest (ROI) in the brain. The ROI E-field can be computationally maximized and set to match a specific reference using individualized head models to find the optimal coil placement and stimulus intensity. However, the available software lacks many practical features for prospective planning of TMS interventions and retrospective evaluation of the experimental targeting accuracy.

ApproachThe TMS targeting and analysis pipeline (TAP) software uses an MRI/fMRI-derived brain target to optimize coil placement considering experimental parameters such as the subjects hair thickness and coil placement restrictions. The coil placement optimization is implemented in SimNIBS 3.2, for which an additional graphical user interface (TargetingNavigator) is provided to visualize and adjust procedural parameters. The coil optimization process also computes the E-field at the target, allowing the selection of the TMS device intensity setting to achieve a specific E-field strength. The optimized coil placement information is prepared for neuronavigation software, which supports targeting during the TMS procedure. The neuronavigation system can record the coil placement during the experiment, and these data can be processed in TAP to quantify the accuracy of the experimental TMS coil placement and induced E-field.

Main resultsTAP was demonstrated in a study consisting of three repetitive TMS sessions in five subjects. TMS was delivered by an experienced operator under neuronavigation with the computationally optimized coil placement. Analysis of the experimental accuracy from the recorded neuronavigation data indicated coil location and orientation deviations up to about 2 mm and 2{degrees}, respectively, resulting in an 8% median decrease in the target E-field magnitude compared to the optimal placement.

SignificanceTAP supports navigated TMS with a variety of features for rigorous and reproducible stimulation delivery, including planning and evaluation of coil placement and intensity selection for E-field-based dosing. Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd. The activity and functions of soil microbial communities in the Finnish sub-Arctic vary across vegetation types Increased microbial activity in high-latitude soils due to climate change might lead to higher greenhouse gas (GHG) emissions. However, mechanisms of microbial GHG production and consumption in tundra soils are not thoroughly understood. We analyzed 116 soil metatranscriptomes from 73 sites in the Finnish sub-Arctic to investigate how the diversity and functional potential of bacterial and archaeal communities vary across vegetation types and soil layers. Soils differed in physicochemical conditions, with meadow soils being characterized by higher pH and low soil organic matter (SOM) and carbon/nitrogen ratio whereas dwarf shrub-dominated ecosystems with high SOM and low pH. Actinobacteria, Acidobacteria, Alphaproteobacteria, and Planctomycetes predominated all communities but there were significant differences on genus level between vegetation types, as plant polymer degrading groups were more active in shrub-dominated soils compared to meadows. Given that climate change scenarios predict expansion in dwarf shrubs at high latitudes, our results indicate that the rate of carbon turnover in tundra soils may increase in the future. Additionally, transcripts of methanotrophs were detected in the mineral layer of all soils, potentially moderating methane fluxes from deeper layers. In all, this study provides new insights into possible shifts in tundra microbial diversity and activity with climate change. Generalization of Cortical MOSTest Genome-Wide Associations Within and Across Samples Genome-Wide Association studies have typically been limited to single phenotypes, given that high dimensional phenotypes incur a large multiple comparisons burden: ~1 million tests across the genome times the number of phenotypes. Recent work demonstrates that a Multivariate Omnibus Statistic Test (MOSTest) is well powered to discover genomic effects distributed across multiple phenotypes. Applied to cortical brain MRI morphology measures, MOSTest has resulted in a drastic improvement in power to discover loci when compared to established approaches (min-P). One question that arises is how well these discovered loci replicate in independent data. Here we perform 10 times cross validation within 35,644 individuals from UK Biobank for imaging measures of cortical area, thickness and sulcal depth (>1,000 dimensionality for each). By deploying a replication method that aggregates discovered effects distributed across multiple phenotypes, termed PolyVertex Score (PVS), we demonstrate a higher replication yield and comparable replication rate of discovered loci for MOSTest (# replicated loci: 348-845, replication rate: 94-95%) in independent data when compared with the established min-P approach (# replicated loci: 31-68, replication rate: 65-80%). An out-of-sample replication of discovered loci was conducted with a sample of 8,336 individuals from the Adolescent Brain Cognitive Development(R) (ABCD) study, who are on average 50 years younger than UK Biobank individuals. We observe a higher replication yield and comparable replication rate of MOSTest compared to min-P. This finding underscores the importance of using well-powered multivariate techniques for both discovery and replication of high dimensional phenotypes in Genome-Wide Association studies. Spatial phylogenetics of Japanese ferns: Patterns, processes, and implications for conservation PremiseBiodiversity is often only measured with species richness. However, this ignores evolutionary history and is not sufficient for making conservation decisions. Here, we characterize multiple facets and drivers of biodiversity to understand how these relate to bioregions and conservation status in the ferns of Japan.

MethodsWe compiled a community dataset of 1,239 20 km x 20 km grid-cells including 672 taxa based on > 300,000 specimen records. We combined this with a phylogeny and functional traits to analyze taxonomic, phylogenetic, and functional diversity, and modeled biodiversity metrics in response to environmental factors and reproductive mode. Hierarchical clustering was used to delimit bioregions. Conservation status and threats were assessed by comparing the overlap of significantly diverse grid-cells with conservation zones and range maps of native Japanese deer.

ResultsTaxonomic richness was highest at mid-latitudes. Phylogenetic and functional diversity and phylogenetic endemism were highest in small southern islands. Relative phylogenetic and functional diversity were high at high and low latitudes, and low at mid-latitudes. Grid-cells were grouped into three (phylogenetic) or four (taxonomic) major bioregions. Temperature and apomixis were identified as drivers of biodiversity patterns. Conservation status was generally high for grid-cells with significantly high biodiversity, but the threat due to herbivory by deer was greater for taxonomic richness than other metrics.

ConclusionsOur integrative approach reveals previously undetected patterns and drivers of biodiversity in the ferns of Japan. Future conservation efforts should recognize that threats can vary by biodiversity metric and consider multiple multiple metrics when establishing conservation priorities. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Protein contact map prediction using multiple sequence alignment dropout and consistency learning for sequences with fewer homologs The prediction of protein contact map needs enough normalized number of effective sequence (Nf) in multiple sequence alignment (MSA). When Nf is small, the predicted contact maps are often not satisfactory. To solve this problem, we randomly selected a small part of sequence homologs for proteins with large Nf to generate MSAs with small Nf. From these MSAs, input features were generated and were passed through a consistency learning network, aiming to get the same results when using the features generated from the MSA with large Nf. The results showed that this method effectively improves the prediction accuracy of protein contact maps with small Nf. Protective host-dependent antagonism among 	Pseudomonas in the Arabidopsis phyllosphere The plant microbiome is a rich biotic environment, comprising numerous taxa. The community structure of these colonizers is constrained by multiple factors, including host-microbe and microbe-microbe interactions, as well as the interplay between the two. While much can be learned from pairwise relationships between individual hosts and microbes, or individual microbes with themselves, the ensemble of interrelations between the host and microbial consortia may lead to different outcomes that are not easily predicted from the individual interactions. Their study can thus provide new insights into the complex relationship between plants and microbes. Of particular importance is how strain-specific such plant-microbe-microbe interactions are, and how they eventually affect plant health. Here, we test strain-level interactions in the phyllosphere between groups of co-existing commensal and pathogenic Pseudomonas among each other and with A. thaliana, by employing synthetic communities of genome-barcoded isolates. We found that commensal Pseudomonas prompted a host response leading to a selective inhibition of a specific pathogenic lineage, resulting in plant protection. The extent of plant protection, however, was dependent on plant genotype, indicating that these effects are host-mediated. There were similar genotype-specific effects on the microbe side, as we could pinpoint an individual Pseudomonas isolate as the predominant cause for this differential interaction. Collectively, our work highlights how within-species genetic differences on both the host and microbe side can have profound effects on host-microbe-microbe dynamics. The paradigm that we have established provides a platform for the study of host-dependent microbe-microbe competition and cooperation in the A. thaliana-Pseudomonas system. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life. Collective polymerase dynamics emerge from DNA supercoiling during transcription All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts and transcription stops. DNA supercoiling can also mediate non-local interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton) whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.

PACS numbers: A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics.