Measuring transcription factor binding and gene expression using barcoded self-reporting transposon calling cards and transcriptomes Calling cards technology using self-reporting transposons enables the identification of DNA-protein interactions through RNA sequencing. Although immensely powerful, current implementations of calling cards in bulk experiments on populations of cells are technically cumbersome and require many replicates to identify independent insertions into the same genomic locus. Here, we have drastically reduced the cost and labor requirements of calling card experiments in bulk populations of cells by introducing a DNA barcode into the calling card itself. An additional barcode incorporated during reverse transcription enables simultaneous transcriptome measurement in a facile and affordable protocol. We demonstrate that barcoded self-reporting transposons recover in vitro binding sites for four basic helix-loop-helix transcription factors with important roles in cell fate specification: ASCL1, MYOD1, NEUROD2, and NGN1. Further, simultaneous calling cards and transcriptional profiling during transcription factor overexpression identified both binding sites and gene expression changes for two of these factors. Lastly, we demonstrated barcoded calling cards can record binding in vivo in the mouse brain. In sum, RNA-based identification of transcription factor binding sites and gene expression through barcoded self-reporting transposon calling cards and transcriptomes is an efficient and powerful method to infer gene regulatory networks in a population of cells.
Detecting cell-of-origin and cancer-specific methylation features of cell-free DNA from Nanopore sequencing The Oxford Nanopore (ONT) platform provides portable and rapid genome sequencing, and its ability to natively profile DNA methylation without complex sample processing is attractive for clinical sequencing. We recently demonstrated ONT shallow whole-genome sequencing to detect copy number alterations (CNA) from the circulating tumor DNA (ctDNA) of cancer patients. Here, we show that cell-type and cancer-specific methylation changes can also be detected, as well as cancer-associated fragmentation signatures. This feasibility study suggests that ONT shallow WGS could be a powerful tool for liquid biopsy, especially real-time medical applications.
Common coupling map advances GPCR-G protein selectivity Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein  couplome. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on GPCR-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling towards the exploitation of G protein signaling bias in design of safer drugs. The cell growth cycle: resource allocation and subcellular pattern formation in budding yeast. Metabolic oscillation percolate throughout cellular physiology. The transcriptome oscillates between gene cohorts encoding for biosynthesis and growth, and for catabolism and stress response. However, long protein half-lives dampen effects on protein abundance, and the function of periodic transcription remains elusive. To test three prevailing hypotheses, we performed RNA-seq analysis during a dynamic state of the system, just preceding a bifurcation of system dynamics to a longer period. A few enzymes with short protein half-live and high transcript abundance amplitude could play a direct role in metabolic dynamics (hypothesis H1), e.g., carbonic anhydrase, sulfate uptake and the glyoxylate cycle. Substantial changes in the relative duration of expression of the antagonistic cohorts precede the bifurcation, supporting the idea of a function in cellular resource allocation (H2). The pulse-width modulation model, a mathematical formulation of this idea, explains the dependence of the cycle period on the culture growth rate. This pulse-like model of cell growth provides a first theoretical framework, where the phenomenon is understood as a mechanism of cellular resource allocation and protein homeostasis. The universal temporal program of transcription encodes the spatial structure of the cell, and we further suggest a function in pattern formation of the growing cell (H3). Alpha satellite insertion close to an ancestral centromeric region Human centromeres are mainly composed of alpha satellite DNA hierarchically organized as higher-order repeats (HORs). Alpha satellite dynamics is shown by sequence homogenization in centromeric arrays and by its transfer to other centromeric locations, for example during the maturation of new centromeres. We identified during prenatal aneuploidy diagnosis by FISH a de novo insertion of alpha satellite DNA from the centromere of chromosome 18 (D18Z1) into cytoband 15q26. Although bound by CENP-B, this locus did not acquire centromeric functionality as demonstrated by lack of constriction and absence of CENP-A binding. The insertion was associated with a 2.8 kbp deletion and likely occurred in the paternal germline. The site was enriched in long terminal repeats (LTRs) and located ~10 Mbp from the location where a centromere was ancestrally seeded and became inactive in the common ancestor of humans and apes 20-25 million years ago. Long read mapping to the T2T-CHM13 human genome assembly revealed that the insertion derives from a specific region of chromosome 18 centromeric 12-mer HOR array in which the monomer size follows a regular pattern. The rearrangement did not directly disrupt any gene or predicted regulatory element and did not alter the methylation status of the surrounding region, consistent with the absence of phenotypic consequences in the carrier. This case demonstrates a likely rare but new class of structural variation that we name  alpha satellite insertion. It also expands our knowledge on alphoid DNA dynamics and conveys the possibility that alphoid arrays can relocate near vestigial centromeric sites. Dynamic interaction of BcsD subunit in type I bacterial cellulose synthase Cellulose is a promising biological material for supporting sustainable human life. This natural polymer is synthesized by cellulose synthase, a protein complex in the cell membrane. Cellulose synthase in bacteria is a hetero-subunit complex, and its subunit organization varies widely depending on the species. In the type I bacterial cellulose synthase complex, the BcsD (bacterial cellulose synthase D) protein is believed to play an important role in producing cellulose with long slender fiber morphology and high crystallinity, given the phenotype of the bcsD-deficient mutant and the specific existence of the type I operon in bacterial species synthesizing crystalline cellulose microfibrils such as Acetobacter. In this study, we successfully established a heterogeneously expressed Bcs protein in Escherichia coli as an experimental model and conducted biochemical studies for the BcsD protein and the other three major subunits of bacterial cellulose synthase, BcsA, BcsB, and BcsC. It has been shown that the BcsD protein interacts with the functionally required minimal subunits of the BcsAB complex, as well as the BcsC protein. Furthermore, it was shown that BcsD interacts with the BcsAB complex in two modes: direct protein-protein interactions and indirect interactions through the product cellulose. The former and latter modes represent the basal and active states of the type I bacterial cellulose synthase, respectively. This dynamic behavior of the BcsD protein is important for the type I bacterial cellulose synthase complex to regulate the crystallization process of cellulose. Modular biosynthesis of plant hemicellulose and its impact on yeast cells BackgroundThe carbohydrate polymers that encapsulate plants cells have benefited humans for centuries and have valuable biotechnological uses. In the past five years, exciting possibilities have emerged in the engineering of polysaccharide-based biomaterials. Despite impressive advances on bacterial cellulose-based hydrogels, comparatively little is known about how plant hemicelluloses can be reconstituted and modulated in cells suitable for biotechnological purposes.

ResultsHere, we assembled cellulose synthase-like A (CSLA) enzymes using an optimized Pichia pastoris platform to produce tunable heteromannan (HM) polysaccharides in yeast. By swapping the domains of plant mannan and glucomannan synthases, we engineered chimeric CSLA proteins that made {beta}-1,4-linked mannan in quantities surpassing those of the native enzymes while minimizing the burden on yeast growth. Prolonged expression of a glucomannan synthase from Amorphophallus konjac was toxic to yeast cells: reducing biomass accumulation and ultimately leading to compromised cell viability. However, an engineered glucomannan synthase as well as CSLA pure mannan synthases and a CSLC glucan synthase did not inhibit growth. Interestingly, Pichia cell size could be increased or decreased depending on the composition of the CSLA protein sequence. HM yield and glucose incorporation could be further increased by co-expressing chimeric CSLA proteins with a MANNAN-SYNTHESIS-RELATED (MSR) co-factor from Arabidopsis thaliana.

ConclusionThe results provide novel routes for the engineering of polysaccharide-based biomaterials that are needed for a sustainable bioeconomy. The characterization of chimeric cellulose synthase-like enzymes in yeast offers an exciting avenue to produce plant polysaccharides in a tunable manner. Furthermore, cells modified with non-toxic plant polysaccharides such as {beta}-mannan offer a modular chassis to produce and encapsulate sensitive cargo such as therapeutic proteins. Patterned apoptosis has an instructive role for local growth and tissue shape in a fast-growing epithelium What regulates organ size and shape remains one of the fundamental mysteries of modern biology. So far, research in this area has primarily focused on deciphering the regulation in time and space of growth and cell division, while the contribution of cell death has been much more neglected. This includes studies of the Drosophila wing imaginal disc, the prospective fly wing which undergoes massive growth during larval stage, and represents one of the best characterised systems for the study of growth and patterning. So far, it has been assumed that cell death was relatively neglectable in this tissue and as a result the pattern of growth was usually attributed to the distribution of cell division. Here, using systematic mapping and registration combined with quantitative assessment of clone size and disappearance, we show for the first time that cell death is not neglectable, and outline a persistent pattern of cell death and clone elimination in the disc. Local variation of cell death is associated with local variation of clone size, pointing to an impact of cell death on local growth which is not fully compensated by proliferation. Using morphometric analyses of adult wing shape and genetic perturbations, we provide evidence that patterned death affects locally and globally adult wing shape and size. This study describes a roadmap for accurate assessment of the contribution of cell death to tissue shape, and outlines for the first time an important instructive role of cell death in modulating quantitatively local growth and the morphogenesis of a fast-growing tissue. Genome annotation of Caenorhabditis briggsae by TEC-RED identifies new exons, paralogs, and conserved and novel operons The nematode Caenorhabditis briggsae is routinely used in comparative and evolutionary studies involving its well-known cousin C. elegans. The C. briggsae genome sequence has accelerated research by facilitating the generation of new resources, tools, and functional studies of genes. While substantial progress has been made in predicting genes and start sites, experimental evidence is still lacking in many cases. Here, we report an improved annotation of the C. briggsae genome using the Trans-spliced Exon Coupled RNA End Determination (TEC-RED) technique. In addition to identifying the 5 ends of expressed genes, we have discovered operons and paralogs. In summary, our analysis yielded 10,243 unique 5 end sequence tags with matches in the C. briggsae genome. Of these, 6,395 were found to represent 4,252 unique genes along with 362 paralogs and 52 previously unknown exons. These genes included 14 that are exclusively trans-spliced in C. briggsae when compared with C. elegans orthologs. A major contribution of this study is the identification of 493 operons, of which two-thirds are fully supported by tags. In addition, two SL1-type operons were discovered. Interestingly, comparisons with C. elegans showed that only 40% of operons are conserved. Of the remaining operons, 73 are novel, including 12 that entirely lack orthologs in C. elegans. Further analysis revealed that four of the 12 novel operons are conserved in C. nigoni. Altogether, the work described here has significantly advanced our understanding of the C. briggsae system and serves as a rich resource to aid biological studies involving this species. Ineffectual AEC1 differentiation from KRT8hi transitional state without fibrosis is associated with fatal COVID-19 ARDS 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. Deep Transfer Learning of Drug Responses by Integrating Bulk and Single-cell RNA-seq data Massively bulk RNA sequencing databases incorporating drug screening have opened up an avenue to inform the optimal clinical application of cancer drugs. Meanwhile, the growing single-cell RNA sequencing (scRNA-seq) data contributes to improving therapeutic effectiveness by studying the heterogeneity of drug responses for cancer cell subpopulations. There is a clear significance in developing computational biology approaches to predict and interpret cancer drug response in single cell data from clinical samples. Here, we introduce scDEAL, a deep transfer learning framework for cancer drug response prediction at single-cell level by integrating large-scale bulk cell line data. The true innovation of scDEAL is to translate cancer cell line drug responses into predicting clinical drug responses via learning relations of gene expressions and drug responses at bulk-level and transfer to predict drug responses in scRNA-seq. Another innovation is the integrated gradient feature interpretation to infer a comprehensive set of signature genes to reveal potential drug resistance mechanisms. We benchmarked scDEAL on six scRNA-seq datasets and indicate its model interpretability through these case studies. We believe that this work may help study cell reprogramming, drug selection, and repurposing for improving therapeutic efficacy. Brain deletion of Serf2 shifts amyloid conformation in a mouse model Loss of protein homeostasis accelerates ageing and contributes to age-related brain diseases like Parkinsons and Alzheimers diseases. Disease-specific proteins become aggregation-prone and form amyloid deposits that characterize disease, but the biological mechanisms driving their transition in cells and tissues are poorly understood. MOAG-4/SERF accelerates aggregation and toxicity of disease-related proteins in worm- and cell models. Whether MOAG-4/SERF affects amyloid pathology in brain, however, has remained unknown. Here we show that brain-specific knockout of SERF2 is viable and shifts aggregation in a mouse model for amyloid-beta aggregation. Without affecting amyloid-beta levels, knockout of SERF2 redirects its assembly into structurally and morphologically distinct amyloid fibrils and deposits. Our results show that SERF modifies the structural outcome of the aggregation process in mammalian brain. With structural polymorphisms recently being associated with disease severity and progression, our results may help to understand the molecular basis for structural diversity and find possibilities for structure-based interventions.

HighlightsO_LILoss of SERF2 slows embryonic development and causes perinatal lethality
C_LIO_LISERF2 affects proliferation in a cell-autonomous fashion
C_LIO_LIBrain-specific Serf2 knockout does not affect viability or A{beta} production
C_LIO_LIBrain deletion of Serf2 shifts the amyloid conformation of A{beta}
C_LI The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. 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. Modeling predicts differences in CAR T cell signaling due to biological variability In recent decades, chimeric antigen receptors (CARs) have been successfully used to generate engineered T cells capable of recognizing and eliminating cancer cells. The structure of CARs frequently includes costimulatory domains, which enhance the T cell response upon antigen encounter. However, it is not fully known how the CAR co-stimulatory domains influence T cell activation in the presence of biological variability. In this work, we used mathematical modeling to elucidate how the inclusion of one such co-stimulatory molecule, CD28, impacts the response of a population of engineered T cells under different sources of variability. Particularly, our simulations demonstrate that CD28-bearing CARs mediate a faster and more consistent population response under both target antigen variability and kinetic rate variability. We identify kinetic parameters that have the most impact on mediating cell activation. Finally, based on our findings, we propose that enhancing the catalytic activity of lymphocyte-specific protein tyrosine kinase (LCK) can result in drastically reduced and more consistent response times among heterogeneous CAR T cell populations. Blebbistatin protects iPSC-CMs from hypercontraction and facilitates automated patch-clamp based electrophysiological study Recently, there have been great advances in cardiovascular channelopathy modeling and drug safety pharmacology using human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). The automated patch-clamp (APC) technique overcomes the disadvantages of the manual patch-clamp (MPC) technique, which is labor intensive and gives low output. However, the application of the APC platform is still limited in iPSC-CM based research, due to the difficulty in maintaining the high quality of single iPSC-CMs during dissociation and recording. In this study, we improved the method for single iPSC-CM preparation by applying 2.5 M blebbistatin (BB, an excitation-contraction coupling uncoupler) throughout APC procedures (dissociation, filtration, storage, and recording). Under non-BB buffered condition, iPSC-CMs in suspension showed a severe bleb-like morphology. However, BB-supplement led to significant improvements in morphology and INa recording, and we even obtained several CMs that showed spontaneous action potentials with typical morphology. Furthermore, APC faithfully recapitulated the single-cell electrophysiological phenotypes of iPSC-CMs derived from Brugada syndrome patients, as detected with MPC. Our study indicates that APC is capable of replacing MPC in the modeling of cardiac channelopathies using human iPSC-CMs by providing high-quality data with higher throughput. Point Mutations in SARS-CoV-2 Variants Induce Long-Range Dynamical Perturbations in Neutralizing Antibodies Monoclonal antibodies are emerging as a viable treatment for the coronavirus disease 19 (COVID-19). However, newly evolved variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can reduce the efficacy of currently available antibodies and can diminish vaccine-induced immunity. Here, we demonstrate that the microscopic dynamics of neutralizing monoclonal antibodies can be profoundly modified by the mutations present in the spike proteins of the SARS-COV-2 variants currently circulating in the world population. The dynamical perturbations within the antibody structure, which alter the thermodynamics of antigen recognition, are diverse and can depend both on the nature of the antibody and on the spatial location of the spike mutation. The correlation between the motion of the antibody and that of the spike receptor binding domain (RBD) can also be changed, modulating binding affinity. Using protein-graph-connectivity networks, we delineated the mutant-induced modifications in the information-flow along allosteric pathway throughout the antibody. Changes in the collective dynamics were spatially distributed both locally and across long-range distances within the antibody. On the receptor side, we identified an anchor-like structural element that prevents the detachment of the antibodies; individual mutations there can significantly affect the antibody binding propensity. Our study provides insight into how virus neutralization by monoclonal antibodies can be impacted by local mutations in the epitope via a change in dynamics. This realization adds a new layer of sophistication to the efforts for rational design of monoclonal antibodies against new variants of SARS-CoV2, taking the allostery in the antibody into consideration. Elevation of major constitutive heat shock proteins is heat shock factor independent and essential for establishment and growth of Lgl loss and Yorkie gain-mediated tumors in Drosophila Cancer cells overexpress heat shock proteins (Hsps), the major components of cellular stress response, to overcome the diversity of stresses and to survive. However, the specific roles of Hsps in initiation and establishment of cancers remain unclear. Using loss of Lgl function-mediated epithelial tumorigenesis in Drosophila, we induced tumorigenic MARCM clones of different genetic backgrounds in wing imaginal discs to examine the temporal and spatial expression and roles of major heat shock proteins in tumor growth. The constitutively expressed Hsp83, Hsc70 (heat shock cognate), Hsp60 and Hsp27 show elevated levels in all cells of the tumorigenic clone since early stages, which persists till their transformation. However, the stress-inducible Hsp70 is expressd only in a few cells at later stage of established tumorous clones that show high F-actin aggregation. Intriguingly, levels of heat shock factor (HSF), the master regulator of Hsps, remain unaltered in these tumorous cells and its down-regulation does not affect tumorigenic growth of lgl- clones overexpressing Yorkie, although down-regulation of Hsp83 prevents their survival and growth. Interestingly, overexpression of HSF or Hsp83 in lgl- cells makes them competitively successful in establishing tumorous clones. These results show that the major constitutively expressed Hsps, but not the stress-inducible Hsp70, are involved in early as well as late stages of epithelial tumors and their elevated expression in lgl- clones co-overexpressing Yorkie is independent of HSF. In Vivo Multi-Day Calcium Imaging of CA1 Hippocampus in Freely Moving Rats Reveals a High Preponderance of Place Cells with Consistent Place Fields Calcium imaging using GCaMP indicators and miniature microscopes has been used to image cellular populations during long timescales and in different task phases, as well as to determine neuronal circuit topology and organization. Because the hippocampus (HPC) is essential for tasks of memory, spatial navigation, and learning, calcium imaging of large populations of HPC neurons can provide new insight on cell changes over time during these tasks. All reported HPC in vivo calcium imaging experiments have been done in mouse. However, rats have many behavioral and physiological experimental advantages over mice. In this paper, we present the first (to our knowledge) in vivo calcium imaging from CA1 hippocampus in freely moving male rats. Using the UCLA Miniscope, we demonstrate that, in rat, hundreds of cells can be visualized and held across weeks. We show that calcium events in these cells are highly correlated with periods of movement, with few calcium events occurring during periods without movement. We additionally show that an extremely large percent of cells recorded during a navigational task are place cells (77.3{+/-}5.0%, surpassing the percent seen during mouse calcium imaging), and that these cells enable accurate decoding of animal position and can be held over days with consistent place fields in a consistent spatial map. A detailed protocol is included, and implications of these advancements on in vivo imaging and place field literature are discussed.

Significance statementIn vivo calcium imaging in freely moving animals allows the visualization of cellular activity across days. In this paper, we present the first in vivo Ca2+ recording from CA1 hippocampus in freely moving rats. We demonstrate that hundreds of cells can be visualized and held across weeks, and that calcium activity corresponds to periods of movement. We show that a high percentage (77.3{+/-}5.0%) of imaged cells are place cells, and that these place cells enable accurate decoding and can be held stably over days with little change in field location. Because the hippocampus is essential for many tasks involving memory, navigation, and learning, imaging of large populations of HPC neurons can shed new insight on cellular activity changes and organization. Multi-omic analysis of the Arabidopsis clock activator mutant rve 4 6 8 reveals connections to carbohydrate metabolism and proteasome regulation Plants are able to sense changes in their light environments, such as the onset of day and night, as well as anticipate these changes in order to adapt and survive. Central to this ability is the plant circadian clock, a molecular circuit that precisely orchestrates plant cell processes over the course of a day. REVEILLE proteins (RVEs) are recently discovered members of the plant circadian circuitry that activate the evening complex and PRR genes to maintain regular circadian oscillation. The RVE 8 protein and its two homologs, RVE 4 and 6, have been shown to limit the length of the circadian period, with rve 4 6 8 triple-knockout plants possessing an elongated period along with increased leaf surface area, biomass, cell size and delayed flowering relative to wild-type Col-0 plants. Here, using a multi-omics approach consisting of phenomics, transcriptomics, proteomics, and metabolomics we draw novel connections between RVE8-like proteins and a number of core plant cell processes. In particular, we reveal that loss of RVE8-like proteins results in altered carbohydrate, organic acid and lipid metabolism, including a starch excess phenotype at dawn. We further demonstrate that rve 4 6 8 plants have lower levels of 20S proteasome subunits and possess significantly reduced proteasome activity, potentially explaining the increase in cell-size observed in RVE8-like mutants. Overall, this robust, multi-omic dataset, provides substantial new insights into the far reaching impact RVE8-like proteins have on the diel plant cell environment. Barbary ground squirrels do not have a sentinel system but instead synchronise vigilance Coordinated behavior, such as hunting in lions and coordinated vigilance as antipredator behavior, are examples of benefits of group-living. Instead of asynchronous vigilance, some social species synchronize their vigilance bouts or take turns acting as sentinels. To increase our knowledge on the evolution of vigilance behavior, we studied whether vigilance is coordinated in Barbary ground squirrels, Atlantoxerus getulus. We show that vigilance was synchronized instead of taking turns. Multiple non-mutually exclusive hypotheses could explain synchronization: Barbary ground squirrels may perch because 1) neighbors are perched (copying effect), 2) perch synchrony may be an emergent property of the ecology as all squirrels may be satiated at the same time (collective behavior), or 3) the benefits are large in terms of evading ambush predators and scanning effectiveness (watch each others back). Particularly, in habitats where the field of view is obstructed by man-made structures and multiple individuals may be necessary to watch for terrestrial predators, synchronized vigilance may have greater fitness benefits than sentinel behavior. We conclude that it is essential to test assumptions of coordination and, thus, to analyze coordination to describe sentinel systems.

Significance StatementVigilance behavior can be vital to an animals survival. Taking turns acting as sentinels or synchronizing vigilance bouts reduces the cost of the trade-off between feeding and predation risk. A sentinel system assumes that sentinels are vigilant from raised positions, warn group members of danger, and alternate vigilance bouts. However, the assumption of alternating vigilance bouts remains poorly tested. We tested this assumption in invasive Barbary ground squirrels. We found that instead of alternating, individuals synchronized their vigilance bouts. Perch synchrony may be 1) a response to perching group members (copying effect), 2) an emergent property of the species ecology, and 3) an adaptation to anthropogenically altered habitats (watch each others back). Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Comparative infection and pathogenesis of SARS-CoV-2 Omicron and Delta variants in aged and young Syrian hamsters Coronavirus disease 2019 continues to batter the world with the unceasing introduction of new variants of the causative virus, SARS-CoV-2. In order to understand differences in disease caused by variants of concern and to develop variant-specific vaccines, suitable small animal models are required that mimic disease progression in humans at various stages of life. In this study, we compared the dynamics of infection with two SARS-CoV-2 variants of concern (Delta and Omicron) in aged (>1 year 3 months old) and young (<5 weeks old) Syrian hamsters (Mesocricetus auratus). We show that no weight loss occurred in Omicron infected groups regardless of age, while infection with the Delta variant caused weight loss of up to 10% by day 7 post-infection with slower and incomplete recovery in the aged group. Omicron replicated to similar levels as Delta in the lungs, trachea and nasal turbinates, with no significant differences in the tissue viral loads of aged versus young animals for either variant. In contrast to rare necrosis observed in Omicron-infected animals regardless of age, severe necrosis was observed in the olfactory epithelium in Delta-infected animals. Omicron infection also resulted in mild pulmonary disease in both young and aged animals compared to the moderate acute necrotizing bronchointerstitial pneumonia seen in Delta-infected animals. These results suggest that Omicron infection results in an attenuated clinical disease outlook in Syrian hamsters compared to infection with the Delta variant irrespective of age. Implicating Gene and Cell Networks Responsible for Differential COVID-19 Host Responses via an Interactive Single Cell Web Portal Numerous studies have provided single-cell transcriptome profiles of host responses to SARS-CoV-2 infection. Critically lacking however is a datamine that allows users to compare and explore cell profiles to gain insights and develop new hypotheses. To accomplish this, we harmonized datasets from COVID-19 and other control condition blood, bronchoalveolar lavage, and tissue samples, and derived a compendium of gene signature modules per cell type, subtype, clinical condition, and compartment. We demonstrate approaches to probe these via a new interactive web portal (http://toppcell.cchmc.org/ COVID-19). As examples, we develop three hypotheses: (1) a multicellular signaling cascade among alternatively differentiated monocyte-derived macrophages whose tasks include T cell recruitment and activation; (2) novel platelet subtypes with drastically modulated expression of genes responsible for adhesion, coagulation and thrombosis; and (3) a multilineage cell activator network able to drive extrafollicular B maturation via an ensemble of genes strongly associated with risk for developing post-viral autoimmunity. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. A northern range extension of a Canadian species of Special Concern, Dielis pilipes (Hymenoptera: Scoliidae), in the Okanagan Valley of British Columbia The only known Canadian records of the yellow scarab hunter wasp, Dielis pilipes (Saussure) (Hymenoptera: Scoliidae), are from the southern Okanagan and Similkameen valleys of British Columbia. We report a 25-km northern range extension of the species into the Ponderosa Pine Biogeoclimatic Ecosystem Classification zone, collected in an unmanaged agricultural field in Summerland, British Columbia. This finding is of conservation importance and has implications for natural biological control of ten-lined June beetles, Polyphylla decemlineata (Say) and P. crinita LeConte (Coleoptera: Scarabaeidae), incidental agricultural pests in the Okanagan.