Large scale in vivo acquisition, segmentation, and 3D reconstruction of cortical vasculature using open-source functional ultrasound imaging platform The brain is composed of a dense and ramified vascular network comprising various sizes of arteries, veins, and capillaries. One way to assess the risk of cerebrovascular pathologies is to use computational models to predict the physiological effects of a reduction of blood supply and correlate these responses with observations of brain damage. Therefore, it is crucial to establish a detailed 3D organization of the brain vasculature, which could be used to develop more accurate in silico models. For this purpose, we have adapted our open-access functional ultrasound imaging platform previously designed for recording brain-wide activity that is now capable of fast and reproducible acquisition, segmentation, and reconstruction of the cortical vasculature. For the first time, it allows us to digitize the cortical vasculature in awake rodents with a [~]100 {micro}m3 spatial resolution. Contrary to most available strategies, our approach can be performed in vivo within minutes. Moreover, it is easy to implement since it neither requires exogenous contrast agents nor long post-processing time. Hence, we performed a cortex-wide reconstruction of the vasculature and its quantitative analysis, including i) classification of descending arteries versus ascending veins in more than 1500 vessels/animal, ii) quick estimation of their length. Importantly, we confirmed the relevance of our approach in a model of cortical stroke, which enables quick visualization of the ischemic lesion. This development contributes to extending the capabilities of ultrasound neuroimaging to understand better cerebrovascular pathologies such as stroke, vascular cognitive impairment, and brain tumors and is highly scalable for the clinic.
Laminar compartmentalization of attention modulation in area V4 aligns with the demands of visual processing hierarchy in the cortex Contrast is a key feature of the visual scene that aids object recognition. Attention has been shown to selectively enhance the responses to low contrast stimuli in visual area V4, a critical hub that sends projections both up and down the visual hierarchy. Veridical encoding of contrast information is a key computation in early visual areas, while later stages encode higher level features that benefit from improved sensitivity to low contrast. How area V4 meets these distinct information processing demands in the attentive state is not known. We found that attentional modulation of contrast responses in area V4 is cortical layer and cell-class specific. Putative excitatory neurons in the superficial output layers that project to higher areas show enhanced boosting of low contrast information. On the other hand, putative excitatory neurons of deep output layers that project to early visual areas exhibit contrast-independent scaling. Computational modeling revealed that such layer-wise differences may result from variations in spatial integration extent of inhibitory neurons. These findings reveal that the nature of interactions between attention and contrast in V4 is highly compartmentalized, in alignment with the demands of the visual processing hierarchy.
Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event The anaerobic ammonium oxidation (anammox) bacteria could transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we performed comprehensive phylogenomic analysis and showed a single origin of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which includes implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based approach, we estimated that anammox bacteria originated at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Ga]) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress, bioenergetics and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria. Fatal COVID-19 ARDS associated with incomplete AEC1 differentiation from the transitional state without senescence or fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. Kidney failure alters parathyroid Pin1 phosphorylation and parathyroid hormone mRNA binding proteins leading to secondary hyperparathyroidism Parathyroid hormone (PTH) regulates calcium metabolism and bone strength. Chronic kidney disease leads to secondary hyperparathyroidism (SHP) which increases morbidity and mortality. In experimental SHP, the increased PTH gene expression is due to enhanced PTH mRNA stability, mediated by changes in its interaction with stabilizing AUF1 and destabilizing KSRP. Pin1 isomerizes target proteins and leads to KSRP dephosphorylation. In SHP, Pin1 isomerase activity is decreased and phosphorylated KSRP fails to bind PTH mRNA, resulting in high PTH mRNA stability and levels. However, the up- and down-stream mechanisms by which kidney failure stimulates the parathyroid remain elusive. We now reveal a pathway where kidney failure induces parathyroid Pin1 phosphorylation, linking Pin1, KSRP and PTH mRNA stability as key players for the high PTH expression in SHP. We identified by mass-spectrometry, changes in rat parathyroid proteome and phosphoproteome profiles induced by impaired renal function, including KSRP phosphorylation at Pin1 target sites. Furthermore, both acute and chronic kidney failure led to parathyroid-specific Pin1 Ser16 and Ser71 phosphorylation, which disrupts Pin1 activity. Accordingly, pharmacologic Pin1 inhibition, that mimics the decreased Pin1 activity in SHP, increased PTH expression ex-vivo in parathyroid organ cultures and in transfected cells, through the PTH mRNA protein-interacting element and KSRP phosphorylation at potential Pin1-binding motifs. Therefore, kidney failure leads to loss of parathyroid Pin1 activity by inducing Pin1 phosphorylation. This predisposes parathyroids to increase PTH production through impaired PTH mRNA decay that is dependent on KSRP phosphorylation. Pin1 and KSRP phosphorylation and the Pin1-KSRP-PTH mRNA axis thus drive SHP. Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Dug: A Semantic Search Engine Leveraging Peer-Reviewed Knowledge to Span Biomedical Data Repositories MotivationAs the number of public data resources continues to proliferate, identifying relevant datasets across heterogenous repositories is becoming critical to answering scientific questions. To help researchers navigate this data landscape, we developed Dug: a semantic search tool for biomedical datasets utilizing evidence-based relationships from curated knowledge graphs to find relevant datasets and explain why those results are returned.

ResultsDeveloped through the National Heart, Lung, and Blood Institutes (NHLBI) BioData Catalyst ecosystem, Dug has indexed more than 15,911 study variables from public datasets. On a manually curated search dataset, Dugs total recall (total relevant results/total results) of 0.79 outperformed default Elasticsearchs total recall of 0.76. When using synonyms or related concepts as search queries, Dug (0.36) far outperformed Elasticsearch (0.14) in terms of total recall with no significant loss in the precision of its top results.

Availability and ImplementationDug is freely available at https://github.com/helxplatform/dug. An example Dug deployment is also available for use at https://search.biodatacatalyst.renci.org/.

Contactawaldrop@rti.org or scox@renci.org Artificial reefs geographical location matters more than shape, age and depth for sessile invertebrate colonization in the Gulf of Lion (NorthWestern Mediterranean Sea) AO_SCPLOWBSTRACTC_SCPLOWArtificial reefs (ARs) have been used to support fishing activities. Sessile invertebrates are essential components of trophic networks within ARs, supporting fish productivity. However, colonization by sessile invertebrates is possible only after effective larval dispersal from source populations, usually in natural habitat. While most studies focused on short term colonization by pioneer species, we propose to test the relevance of geographic location, shape, age and depth of immersion on the ARs long term colonization by species found in natural stable communities in the Gulf of Lion. We recorded the presence of five sessile invertebrates species, with contrasting life history traits and regional distribution in the natural rocky habitat, on ARs with different shapes deployed during two immersion time periods (1985 and the 2000s) and in two depth ranges (<20m and >20m). At the local level ([~]5kms), neither shape, depth nor immersion duration differentiated ARs assemblages. At the regional scale (>30kms), colonization patterns differed between species, resulting in diverse assemblages. This study highlights the primacy of geographical positioning over shape, immersion duration and depth in ARs colonization, suggesting it should be accounted for in maritime spatial planning. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Designer Scaffolds for Interfacial Bioengineering In regenerative medicine, the healing of the interfacial zone between tissues is a major challenge, yet approaches for studying the complex microenvironment of this interface remain lacking. Here, we study these complex living interfaces by manufacturing modular "blocks" of naturally porous decellularized plant-derived scaffolds of varying shapes and sizes with a computer numerical controlled mill. We demonstrate how each scaffold can be seeded with different cell types and easily assembled in a manner akin to LEGO bricks to create an engineered tissue interface (ETI). Cells were observed to migrate across the interface formed between an empty scaffold and a scaffold pre-seeded with cells. However, when both scaffolds contained cells, only a shallow cross-over zone of cell infiltration formed at the interface. As a proof-of-concept study, we utilized ETIs to investigate the interaction between lab grown bone and connective tissues. Consistent with the above, a cross-over zone of the two distinct cell types formed at the interface between scaffolds, otherwise the populations remained distinct. Finally, we demonstrate how ETIs are biocompatible in vivo, becoming vascularized and integrated into surrounding tissue after implantation. This work creates new tissue design avenues for understanding biological processes or the development of synthetic artificial tissues. Feasibility of Administering Human Pancreatic Cancer Chemotherapy in a Spontaneous Pancreatic Cancer Mouse Model BackgroundBoth modified FOLFIRINOX (mFFX) and gemcitabine/nab-paclitaxel chemotherapy regimens have been shown to improve clinical outcomes in patients with pancreatic cancer, and are often used interchangeably as the standard of care. Preclinical studies often do not use these regimens, since administering these multiagent approaches can be difficult. In this study, we assessed the feasibility of administering these two chemotherapy regimens in spontaneous pancreatic tumors using KPC mice with the ultimate goal of advancing preclinical studies.

MethodsKPC mice were created by breeding KrasLSL-G12D/+ to Trp53fl/fl;Ptf1Cre/+, resulting in KrasLSL-G12D/+;p53fl/+;Ptf1Cre/+ mice. At 14 weeks of age, mice were palpated for spontaneous tumor growth that was verified using ultrasounds. Mice with tumors under 15 mm in diameter were used. The mice were assigned to one of seven treatment regimens: 1 cycle of mFOLFIRINOX (FFX X1), 2 cycles of mFOLFIRINOX (FFX X2), 1 cycle of mFOLFIRINOX with 40Gy SBRT (FFX SBRT), 1 cycle of gemcitabine/nab-paclitaxel (GEM/AB X1), 2 cycles of gemcitabine/nab-paclitaxel (GEM/AB X2), 2 cycles of gemcitabine/nab-paclitaxel with 40Gy SBRT (GEM/AB SBRT), or no chemotherapy treatment (control).

ResultsIn total, 92 mice were included. The median OS in the FFX X2 group was slightly longer that the median OS in the FFX X1 group (15 days vs 11 days, P=0.003). Mice in the GEM/AB X2 group had longer OS when compared to mice in the GEM/AB X1 group (33.5 vs 13 days, P=0.001). Mice treated with chemotherapy survived longer than untreated control animals (median OS: 6.5 days, P <0.001). Moreover, in mice treated with chemotherapy, mice that received 2 cycles of GEM/AB X2 had the longest survival, while the FFX X1 group had the poorest OS (P <0.001). Lastly, chemotherapy followed by consolidative SBRT trended towards increased local control and survival.

ConclusionsWe demonstrate the utility and feasibility of clinicall relevant mFOLFIRINOX and gemcitabine/nab-paclitaxel in preclinical models of pancreatic cancer. Teosinte mexicana introgression modulates phosphatidylcholine levels and maize flowering time Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp.parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability, and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern US, Canada and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time. 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. Protein degradation sets the fraction of active ribosomes at vanishing growth Growing cells adopt common basic strategies to achieve optimal resource allocation under limited resource availability. Our current understanding of such "growth laws" neglects degradation, assuming that it occurs slowly compared to the cell cycle duration. Here we argue that this assumption cannot hold at slow growth, leading to important consequences. We propose a simple framework showing that at slow growth protein degradation is balanced by a fraction of "maintenance" ribosomes. Consequently, active ribosomes do not drop to zero at vanishing growth, but as growth rate diminishes, an increasing fraction of active ribosomes performs maintenance. Through a detailed analysis of compiled data, we show that the predictions of this model agree with data from E. coli and S. cerevisiae. Intriguingly, we also find that protein degradation increases at slow growth, which we interpret as a consequence of active waste management and/or recycling. Our results highlight protein turnover as an underrated factor for our understanding of growth laws across kingdoms.

Author summaryThe idea that simple quantitative relationships relate cell physiology to cellular composition dates back to the 1950s, but the recent years saw a leap in our understanding of such "growth laws", with relevant implications regarding the interdependence between growth, metabolism and biochemical networks. However, recent works on nutrient-limited growth mainly focused on laboratory conditions that are favourable to growth. Thus, our current mathematical understanding of the growth laws neglects protein degradation, under the argument that it occurs slowly compared to the timescale of the cell cycle. Instead, at slow growth the timescales of mass loss from protein degradation and dilution become comparable. In this work, we propose that protein degradation shapes the quantitative relationships between ribosome allocation and growth rate, and determines a fraction of ribosomes that do not contribute to growth and need to remain active to balance degradation. A chromosome-level assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and improves time of divergence estimates for the Cannabaceae family O_LIWe present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly resulted in fragmented assemblies. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur.
C_LIO_LIThe chromosome-level assembly of Cascade was developed by scaffolding the previously-reported Cascade assembly generated with PacBio long-read sequencing, and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations, and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assess molecular evolution in gene sequences, gene family expansion and contraction, and time divergence using Bayesian inference.
C_LIO_LIWe identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families.
C_LIO_LIThe Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.
C_LI GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. Shedding Light on the Circadian Clock of the Threespine Stickleback The circadian clock is an internal timekeeping system shared by most organisms, and knowledge about its functional importance and evolution in natural environments is still needed. Here, we investigated the circadian clock of wild-caught threespine sticklebacks (Gasterosteus aculeatus) at the behavioural and molecular levels. While their behaviour, ecology, and evolution are well studied, information on their circadian rhythms are scarce. We quantified the daily locomotor activity rhythm under a light-dark cycle (LD) and under constant darkness (DD). Under LD, all fish exhibited significant daily rhythmicity, while under DD, only 18% of individuals remained rhythmic. This interindividual variation suggests that the circadian clock controls activity only in certain individuals. Moreover, under LD, some fish were almost exclusively nocturnal, while others were active around the clock. Furthermore, the most nocturnal fish were also the least active. These results suggest that light masks activity (i.e. suppresses activity without entraining the internal clock) more strongly in some individuals than others. Finally, we quantified the expression of five clock genes in the brain of sticklebacks under DD using qPCR. We did not detect circadian rhythmicity, which could either indicate that the clock molecular oscillator is highly light-dependent, or that there was an oscillation but that we were unable to detect it. Overall, our study suggests that a strong circadian control on behavioural rhythms may not necessarily be advantageous in a natural population of sticklebacks and that the daily phase of activity varies greatly between individuals because of a differential masking effect of light.

Summary statementWe found that in wild-caught threespine sticklebacks, the circadian clock does not control locomotor activity in most, but not all, individuals. Sticklebacks are mostly nocturnal, although interindividual variation exists. Distinct SARS-CoV-2 antibody reactivity patterns elicited by natural infection and mRNA Vaccination We analyzed data from two ongoing COVID-19 longitudinal serological surveys in Orange County, CA., between April 2020 and March 2021. A total of 8,476 finger stick blood specimens were collected before and after an aggressive mRNA vaccination campaign. IgG levels were determined using a multiplex antigen microarray containing 10 SARS-CoV-2 antigens, 4 SARS, 3 MERS, 12 Common CoV, and 8 Influenza antigens. Twenty-six percent of 3,347 specimens from unvaccinated Orange County residents in December 2020 were SARS-CoV-2 seropositive. The Ab response was predominantly against nucleocapsid (NP), full length spike and the spike S2 domain. Anti-receptor binding domain (RBD) reactivity was low and there was no cross-reactivity against SARS S1 or SARS RBD. An aggressive mRNA vaccination campaign at the UCI Medical Center started on December 16, 2020 and 6,724 healthcare workers were vaccinated within 3 weeks. Seroprevalence increased from 13% in December to 79% in January, 93% in February and 99% in March. mRNA vaccination induced much higher Ab levels especially against the RBD domain and significant cross-reactivity against SARS RBD and S1 was also observed. Nucleocapsid protein Abs can be used to distinguish individuals in a population of vaccinees to classify those who have been previously infected and those who have not, because nucleocapsid is not in the vaccine. Previously infected individuals developed higher Ab titers to the vaccine than those who have not been previously exposed. These results indicate that mRNA vaccination rapidly induces a much stronger and broader Ab response than SARS-CoV-2 infection. The Skeletal Muscle Circadian Clock Regulates Titin Splicing Through RBM20 Circadian rhythms are maintained by a cell autonomous, transcriptional-translational feedback loop known as the molecular clock. While previous research suggests a role of the molecular clock in regulating skeletal muscle structure and function, no mechanisms have connected the molecular clock to sarcomere filaments. Utilizing inducible, skeletal muscle specific, Bmal1 knockout (iMSBmal1-/-) mice, we showed that knocking out skeletal muscle clock function alters titin isoform expression through both RNAseq and LC-MS. This alteration in titins spring length resulted in sarcomere length heterogeneity. We demonstrate the direct link between altered titin splicing and sarcomere length in vitro using U7 snRNPs that truncate the region of titin altered in iMSBmal1-/- muscle. Finally, we identified a mechanism whereby the skeletal muscle clock regulates titin isoform expression through transcriptional regulation of Rbm20, a potent splicing regulator of titin. Our findings demonstrate the importance of the skeletal muscle circadian clock in maintaining titin isoform through regulation of RBM20 expression. Because circadian rhythm disruption is a feature of many chronic diseases, our results highlight a novel pathway that could be targeted to maintain skeletal muscle structure and function in a range of pathologies. SARS-CoV-2 Beta Variant Substitutions Alter Spike Glycoprotein Receptor Binding Domain Structure and Stability The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) and the subsequent COVID-19 pandemic has visited a terrible cost on the world in the forms of disease, death, and economic turmoil. The rapid development and deployment of extremely effective vaccines against SARS-CoV-2 have seemingly brought within reach the end of the pandemic. However, the virus has acquired mutations; and emerging variants of concern (VOC) are more infectious and reduce the efficacy of existing vaccines. While promising efforts to combat these variants are underway, the evolutionary pressures leading to these variants are poorly understood. To that end, here we have studied the effects on the structure and function of the SARS-CoV-2 spike glycoprotein receptor-binding domain of three amino-acid substitutions found in several variants of concern, including alpha (B.1.1.7), beta (B.1.351), and gamma (P.1). We found that these substitutions alter the RBD structure, stability, and ability to bind to ACE2, in such a way as to possibly have opposing and compensatory effects. These findings provide new insights into how these VOC may have been selected for infectivity while maintaining the structure and stability of the receptor binding domain. 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. Pulled Diversification Rates, Lineages-Through-Time Plots and Modern Macroevolutionary Modelling AO_SCPLOWBSTRACTC_SCPLOWEstimating time-dependent rates of speciation and extinction from dated phylogenetic trees of extant species (timetrees), and determining how and why they vary, is key to understanding how ecological and evolutionary processes shape biodiversity. Due to an increasing availability of phylogenetic trees, a growing number of process-based methods relying on the birth-death model have been developed in the last decade to address a variety of questions in macroevolution. However, this methodological progress has regularly been criticised such that one may wonder how reliable the estimations of speciation and extinction rates are. In particular, using lineages-through-time (LTT) plots, a recent study (Louca and Pennell, 2020) has shown that there are an infinite number of equally likely diversification scenarios that can generate any timetree. This has lead to questioning whether or not diversification rates should be estimated at all. Here we summarize, clarify, and highlight technical considerations on recent findings regarding the capacity of models to disentangle diversification histories. Using simulations we demonstrate the characteristics of newly-proposed "pulled rates" and their utility. We recognize that the recent findings are a step forward in understanding the behavior of macroevolutionary modelling, but they in no way suggest we should abandon diversification modelling altogether. On the contrary, the study of macroevolution using phylogenetic trees has never been more exciting and promising than today. We still face important limitations in regard to data availability and methodological shortcomings, but by acknowledging them we can better target our joint efforts as a scientific community. Extracellular Matrix Hydrogels Promote Expression of Muscle-Tendon Junction Proteins Muscle and tendon injuries are prevalent and range from minor sprains and strains to traumatic, debilitating injuries. However, the interactions between these tissues during injury and recovery remain unclear. Three-dimensional tissue models that incorporate both tissues and a physiologically relevant junction between muscle and tendon may help understand how the two tissues interact. Here, we use tissue specific extracellular matrix (ECM) derived from muscle and tendon to determine how cells of each tissue interact with the microenvironment of the opposite tissue resulting in junction specific features. ECM materials were derived from the Achilles tendon and gastrocnemius muscle, decellularized, and processed to form tissue specific pre-hydrogel digests. ECM materials were unique in respect to protein composition and included many types of ECM proteins, not just collagens. After digestion and gelation, ECM hydrogels had similar complex viscosities which were less than type I collagen hydrogels at the same concentration. C2C12 myoblasts and tendon fibroblasts were cultured in tissuespecific ECM conditioned media or encapsulated in tissue-specific ECM hydrogels to determine cell-matrix interactions and the effects on a muscle-tendon junction marker, paxillin. ECM conditioned media had only a minor effect on upregulation of paxillin in cells cultured in monolayer. However, cells cultured within ECM hydrogels had 50-70% higher paxillin expression than cells cultured in type I collagen hydrogels. Contraction of the ECM hydrogels varied by the type of ECM used. Subsequent experiments with varying density of type I collagen (and thus contraction) showed no correlation between paxillin expression and the amount of gel contraction, suggesting that a constituent of the ECM was the driver of paxillin expression in the ECM hydrogels. In addition, the extracellular matrix protein type XXII collagen had similar expression patterns as paxillin, with smaller effect sizes. Using tissue specific ECM allowed for the de-construction of the cell-matrix interactions similar to muscletendon junctions to study the expression of MTJ specific proteins.

Impact StatementThe muscle-tendon junction is an important feature of muscle-tendon units; however, despite crosstalk between the two tissue types, it is overlooked in current research. Deconstructing the cell-matrix interactions will provide the opportunity to study significant junction specific features and markers that should be included in tissue models of the muscletendon unit, while gaining a deeper understanding of the natural junction. This research aims to inform future methods to engineer a more relevant multi-tissue platform to study the muscletendon unit. Pak1 and PP2A antagonize aPKC function to support cortical tension induced by the Crumbs-Yurt complex The Drosophila polarity protein Crumbs is essential for the establishment and growth of the apical domain in epithelial cells. The protein Yurt limits the ability of Crumbs to promote apical membrane growth, thereby defining proper apical/lateral membrane ratio that is crucial for forming and maintaining complex epithelial structures such as tubes or acini. Here, we show that Yurt also increases Myosin-dependent cortical tension downstream of Crumbs. Yurt overexpression thus induces apical constriction in epithelial cells. The kinase aPKC phosphorylates Yurt, thereby dislodging the latter from the apical domain and releasing apical tension. In contrast, the kinase Pak1 promotes Yurt dephosphorylation through activation of the phosphatase PP2A. The Pak1- PP2A module thus opposes aPKC function and supports Yurt-induced apical constriction. Hence, the complex interplay between Yurt, aPKC, Pak1 and PP2A contributes to the functional plasticity of Crumbs. Overall, our data increase our understanding of how proteins sustaining epithelial cell polarization and Myosin-dependent cell contractility interact with one another to control epithelial tissue architecture.