Integrative analysis reveals aged clonal B cells, microenvironment and c-Myc activation in the origin of age-related lymphoma While cancer is an age-related disease, most studies focus on genetically engineered younger mouse models. Here, we uncover how cancer develops as a consequence of the naturally aged immune system in mice. B-cell lymphoma frequently occurs in aged mice and is associated with increased cell size, splenomegaly, and a novel clonal B-cell population. Age-emergent B cells clonally expand outside of germinal centers driven by somatic mutations, activated c-Myc and hypermethylated promoters, and both genetically and epigenetically recapitulate human follicular and diffuse-large B-cell lymphomas. Mechanistically, mouse cancerous B cells originate from age-associated B cells, which are atypical memory B cells. Age-associated B cells secrete a spectrum of proinflammatory cytokines and activate paracrinally the expression of c-Myc in surrounding B cells. Although clonal B cells are a product of an aging microenvironment, they evolve being self-sufficient and support malignancy when transferred into young mice. Inhibition of mTOR and c-Myc attenuates premalignant changes in B cells during aging and emerges as a therapeutic strategy to delay the onset of age-related lymphoma. Together, we uncover how aging generates cancerous B cells, characterize a model that captures the origin of spontaneous cancer during aging and identify interventions that may postpone age-associated lymphoma.
Transcriptional regulation of HSCs in Aging and MDS reveals DDIT3 as a Potential Driver of Dyserythropoiesis Myelodysplastic syndromes (MDS) are hematopoietic stem cell (HSC) malignancies characterized by ineffective hematopoiesis, with increased incidence in elderly individuals. In this work, we analyzed the transcriptome of human HSCs purified from young and elderly healthy donors, as well as MDS patients, identifying transcriptional alterations following eight different patterns of expression. While aging-associated lesions seemed to predispose HSCs to myeloid transformation, disease-specific alterations may trigger MDS development. Among MDS-specific lesions, we detected the upregulation of the transcription factor DDIT3. Overexpression of DDIT3 in human healthy HSCs induced an MDS-like transcriptional state, and a delay in erythropoiesis. Such effect was associated with downregulation of transcription factors required for normal erythropoiesis, and with a failure in the activation of their transcriptional programs. Moreover, DDIT3 knockdown in CD34+ cells from MDS patients with anemia was able to restore erythropoiesis. These results identify DDIT3 as a driver of dyserythropoiesis, and a potential therapeutic target to restore the inefficient erythropoiesis characterizing MDS patients.

STATEMENT OF SIGNIFICANCEThis study defines how human aging and MDS development are characterized by transcriptional alterations in HSCs that follow different patterns, some of which may contribute to myeloid transformation. Among them, we demonstrate how MDS-specific upregulation of DDIT3 in HSCs induces dyserythropoiesis, while its knockdown in HSPCs from MDS patients restores proper erythroid differentiation.
The Cif proteins from Wolbachia prophage WO modify sperm genome integrity to establish cytoplasmic incompatibility Inherited microorganisms can selfishly manipulate host reproduction to drive through populations. In Drosophila melanogaster, germline expression of the native Wolbachia prophage WO proteins CifA and CifB cause cytoplasmic incompatibility (CI) in which embryos from infected males and uninfected females suffer catastrophic mitotic defects and lethality; however, in infected females, CifA expression rescues the embryonic lethality and thus imparts a fitness advantage to the maternally-transmitted Wolbachia. Despite widespread relevance to sex determination, evolution, and vector control, the mechanisms underlying when and how CI impairs male reproduction remain unknown and a topic of debate. Here we use cytochemical, microscopic, and transgenic assays in D. melanogaster to demonstrate that CifA and CifB proteins of wMel localize to nuclear DNA throughout the process of spermatogenesis. Cif proteins cause abnormal histone retention in elongating spermatids and protamine deficiency in mature sperms that travel to the female reproductive tract with Cif proteins. Notably, protamine gene knockouts enhance wild type CI. In ovaries, CifA localizes to germ cell nuclei and cytoplasm of early-stage egg chambers, however Cifs are absent in late-stage oocytes and subsequently in fertilized embryos. Finally, CI and rescue are contingent upon a newly annotated CifA bipartite nuclear localization sequence. Together, our results strongly support the Host Modification model of CI in which Cifs initially modify the paternal and maternal gametes to bestow CI-defining embryonic lethality and rescue. Delta-range coupling between prefrontal cortex and hippocampus supported by respiratory rhythmic input from the olfactory bulb in freely behaving rats An explosion of recent findings firmly demonstrated that brain activity and cognitive function in rodents and humans are modulated synchronously with nasal respiration. Rhythmic respiratory (RR) coupling of wide-spread forebrain activity was confirmed using advanced techniques, including current source density analysis, single unit firing, and phase modulation of local gamma activity, creating solid premise for investigating how higher networks use this mechanism in their communication. Here we show essential differences in the way prefrontal cortex (PFC) and hippocampus (HC) process the RR signal from the olfactory bulb (OB) allowing dynamic PFC-HC coupling utilizing this input. We used inter-regional coherences and their correlations in rats, breathing at low rate ([~]2 Hz) at rest, outside of the short sniffing bouts. We found strong and stable OB-PFC coherence, contrasting OB-HC coherence which was low but highly variable. PFC-HC coupling, however, primarily correlated with the latter, indicating that HC access to the PFC output is dynamically regulated by the responsiveness of HC to the common rhythmic drive. This pattern was present in both theta and non-theta states of waking, whereas PFC-HC communication appeared protected from RR synchronization in sleep states. The findings help to understand the mechanism of rhythmic modulation of non-olfactory cognitive processes by the on-going regular respiration, reported in rodents as well as humans. These mechanisms may be impaired when nasal breathing is limited or in OB-pathology, including malfunctions of the OB epithelium due to infections, such as in COVID-19. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. The adapted Activity-By-Contact model for enhancer-gene assignment and its application to single-cell data Identifying regulatory regions in the genome is of great interest for understanding the epigenomic landscape in cells. One fundamental challenge in this context is to find the target genes whose expression the regulatory regions affect. A recent successful method is the Activity-By-Contact (ABC) model [10] which scores enhancer-gene interactions based on enhancer activity and the contact frequency of an enhancer to its target gene. It requires two types of assays to measure enhancer activity, which limits the applicability. The ABC score describes regulatory interactions entirely from a genes perspective, and does not account for all the candidate target genes of an enhancer. Moreover, there is no implementation available which would allow for an integration with transcription factor (TF) binding information or an efficient analysis of single-cell data. We show that the ABC-model is comparably accurate with only one assay to measure enhancer activity. Further, we demonstrate that it can yield a higher accuracy by adapting the enhancer activity according to the number of contacts the enhancer has to its candidate target genes. We combined our adapted ABC-model with TF binding information and illustrate an analysis of a single-cell ATAC-seq data set of the human heart. We were able to characterise cell type-specific regulatory interactions as well as to prioritise candidate TFs that drive cell type-specific expression. All executed processing steps are incorporated into our new computational pipeline STARE. The software is available at https://github.com/schulzlab/STARE. The maize gene maternal derepression of r1 (mdr1) encodes a DNA glycosylase that demethylates DNA and reduces siRNA expression in endosperm Demethylation of transposons can activate expression of nearby genes and cause imprinted gene expression in endosperm, and it is hypothesized to lead to expression of transposon siRNAs that reinforce silencing in the next generation through transfer either into egg or embryo. Here we describe maternal derepression of r1 (mdr1), which encodes a DNA glycosylase with homology to Arabidopsis DEMETER and which is partially responsible for demethylation of thousands of regions in endosperm. Instead of promoting siRNA expression in endosperm, MDR1 activity inhibits it. Methylation of most repetitive DNA elements in endosperm is not significantly affected by MDR1, with an exception of Helitrons. While maternally-expressed imprinted genes preferentially overlap with MDR1 demethylated regions, the majority of genes that overlap demethylated regions are not imprinted. Double mutant megagametophytes lacking both MDR1 and its close homolog DNG102 result in early seed failure, and double mutant microgametophytes fail pre-fertilization. These data establish DNA demethylation by glycosylases as essential in maize endosperm and pollen and suggest that neither transposon repression nor genomic imprinting are its main function in endosperm. Immune complex solubility and size govern Fc-gamma receptor responses: A scalable cell-based reporter system Fc-gamma receptor (Fc{gamma}R) activation by soluble IgG immune complexes (sICs) represents a major mechanism of inflammation in certain autoimmune diseases such as systemic lupus erythematosus (SLE). A robust and scalable test system allowing for the detection and quantification of sIC bioactivity is missing. Previously described Fc{gamma}R interaction assays are limited to certain Fc{gamma}Rs, lack scalability and flexibility, are not indicative of receptor activation or lack sensitivity towards sIC size. We developed a comprehensive reporter cell panel detecting individual activation of Fc{gamma}Rs from humans and the mouse. The reporter cell lines were integrated into an assay format that provides flexible read-outs enabling the quantification of sIC reactivity via ELISA or a fast detection using flow cytometry. This identified Fc{gamma}RIIA(H) and Fc{gamma}RIIIA as the most sIC-sensitive Fc{gamma}Rs in our test system. Applying the assay we demonstrate that sICs versus immobilized ICs are fundamentally different Fc{gamma}R-ligands with regard to Fc{gamma}R preference and signal strength. Reaching a detection limit in the very low nanomolar range, the assay proved also to be sensitive to sIC stoichiometry and size enabling for the first time a complete reproduction of the Heidelberger-Kendall precipitation curve in terms of immune receptor activation. Analyzing sera from SLE patients and mouse models of lupus and arthritis proved that sIC-dependent Fc{gamma}R activation has predictive capabilities regarding severity of SLE disease. The new methodology provides a sensitive, scalable and comprehensive tool to evaluate the size, amount and bioactivity of sICs in all settings.

One Sentence SummaryIn this study we established a comprehensive Fc{gamma}R reporter cell assay enabling the detection and quantification of soluble immune complexes generated in experimental and clinical settings. Structure-function analysis of the nsp14 N7-guanine methyltransferase reveals an essential role in Betacoronavirus replication As coronaviruses (CoVs) replicate in the host cell cytoplasm, they rely on their own capping machinery to ensure the efficient translation of their mRNAs, protect them from degradation by cellular 5 exoribonucleases, and escape innate immune sensing. The CoV nonstructural protein 14 (nsp14) is a bi-functional replicase subunit harboring an N-terminal 3'-to-5' exoribonuclease (ExoN) domain and a C-terminal (N7-guanine)-methyltransferase (N7-MTase) domain that is presumably involved in viral mRNA capping. Here, we aimed to integrate structural, biochemical, and virological data to assess the importance of conserved N7-MTase residues for nsp14s enzymatic activities and virus viability. We revisited the crystal structure of severe acute respiratory syndrome (SARS)-CoV nsp14 to perform an in silico comparative analysis between betacoronaviruses. We identified several residues likely involved in the formation of the N7-MTase catalytic pocket, which presents a fold distinct from the Rossmann fold observed in most known MTases. Next, for SARS-CoV and Middle East respiratory syndrome-CoV, site-directed mutagenesis of selected residues was used to assess their importance for in vitro enzymatic activity. Most of the engineered mutations abolished N7-MTase activity, while not affecting nsp14-ExoN activity. Upon reverse engineering of these mutations into different betacoronavirus genomes, we identified two substitutions (R310A and F426A in SARS-CoV nsp14) abrogating virus viability and one mutation (H424A) yielding a crippled phenotype across all viruses tested. Our results identify the N7-MTase as a critical enzyme for betacoronavirus replication and define key residues of its catalytic pocket that can be targeted to design inhibitors with a potential pan-coronaviral activity spectrum.

Significance StatementThe ongoing SARS-CoV-2 pandemic emphasizes the urgent need to develop efficient broad-spectrum anti-CoV drugs. The structure-function characterization of conserved CoV replicative enzymes is key to identifying the most suitable drug targets. Using a multidisciplinary comparative approach and different betacoronaviruses, we characterized the key conserved residues of the nsp14 (N7-guanine)-methyltransferase, a poorly defined subunit of the CoV mRNA-synthesizing machinery. Our study highlights the unique structural features of this enzyme and establishes its essential role in betacoronavirus replication, while identifying two residues that are critical for the replication of the four betacoronaviruses tested, including SARS-CoV-2. Microsecond simulations and CD spectroscopy reveals the intrinsically disordered nature of SARS-CoV-2 Spike-C-terminal cytoplasmic tail (residues 1242-1273) in isolation All available SARS-CoV-2 spike protein crystal and cryo-EM structures have shown missing electron densities for cytosolic C-terminal regions. Generally, the missing electron densities point towards the intrinsically disordered nature of the protein region. This curiosity has led us to investigate the C terminal cytosolic region of the spike glycoprotein of SARS-CoV-2 in isolation. The cytosolic regions is supposed to be from 1235-1273 residues or 1242-1273 residues depending on what prediction tool we use. Therefore, we have demonstrated the structural conformation of cytosolic region and its dynamics through computer simulations up to microsecond timescale using OPLS and CHARMM forcefields. The simulations have revealed the unstructured conformation of cytosolic region. Also, in temperature dependent replica-exchange molecular dynamics simulations it has shown its conformational dynamics. Further, we have validated our computational observations with circular dichroism (CD) spectroscopy-based experiments and found its signature spectra at 198 nm which is also adding the analysis as its intrinsically disordered nature. We believe that our findings will surely help us understand the structure-function relationship of the spike proteins cytosolic region.

Graphical Abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY ATP-binding cassette protein ABCF1 couples gene transcription with maintenance of genome integrity in embryonic stem cells OCT4 and SOX2 confer pluripotency by recruiting coactivators to activate stem cell-specific gene expression programs. However, the composition of coactivator complexes and their roles in maintaining stem cell fidelity remain unclear. Here we report the identification of ATP-binding cassette subfamily F member 1 (ABCF1) as a critical coactivator for OCT4/SOX2. ABCF1 is required for pluripotency gene expression and stem cell self-renewal. ABCF1 binds co-dependent coactivators XPC and DKC1 via its intrinsically disordered region and stimulates transcription by linking SOX2 to the transcription machinery. Furthermore, in response to pathogen infection and DNA damage, ABCF1 binds intracellular DNAs accumulated in cells, concomitant with loss of SOX2 interaction and pluripotency gene transcription. This results in spontaneous differentiation of compromised stem cells and elimination from the self-renewing population. Thus, ABCF1 directly couples pluripotency gene transcription with sensing aberrant DNAs and acts as a checkpoint for self-renewal to safeguard stem cell fidelity and genome integrity. Visualization of lake succession as indicated by macrophyte species composition and trophic levels in lakes subdivided by volcanic eruption at Akan Caldera, Hokkaido, Japan It is not feasible to continuously observe ecological succession in lakes because of the long time-scales generally involved. Thus, the process has been inductively deduced by comparing many lakes with different succession states, or indirectly simulated by tracking studies of smaller water bodies, experiments using microcosms or mesocosms, and reconstruction of lake history by sediment analysis. However, the reality of succession processes in large lakes with slow succession is not well understood, and new approaches are needed. Theoretically, in a group of large and small lakes of similar ages and with similar initial and watershed environments, the rate of nutrient accumulation in each lake depends on the ratio of watershed area to lake size, and the lakes are predicted to evolve to different trophic levels over time. Here, we tested this hypothesis on the 10 lakes of varying sizes in Akan Caldera, Japan, which were formed thousands of years ago by fragmentation due to volcanic eruptions within the caldera. Topographic and water quality studies showed that the ratio of accumulated watershed area to lake size (area and volume), expressed logarithmically, had a positive linear regression with the total phosphorus concentration, an indicator of trophic level. The trophic levels of the lakes were diverse, including oligotrophic, mesotrophic, and eutrophic types in the traditional "lake type" classification based on total phosphorus and chlorophyll-a concentrations. Furthermore, 21 species of aquatic macrophytes were observed by a diving survey, and the plant species composition was classified into five groups corresponding to the trophic status of the lakes, indicating a conventional "hydrarch succession". The diversity of water quality and aquatic vegetation in a group of lakes with similar origins paves the way for new comparative studies of lakes, including large lakes. Allopurinol Blocks The Formation And Progression Of Aortic Aneurysm In A Mouse Model Of Marfan Syndrome BackgroundThe pathogenesis and progression of aortic aneurysm in Marfan syndrome (MFS) involves dysregulated TGF-{beta} and nitric oxide signaling, altered hemodynamics, and biomechanical forces. Increasing evidence indicates that redox stress participates in MFS aortopathy, though its contribution is not well established. We reported elevated reactive oxygen species (ROS) formation and NADPH oxidase NOX4 upregulation in MFS patients and mice aortic samples. Here we address the contribution of xanthine oxidoreductase (XOR), which is transcribed as dehydrogenase (XDH) and then post-translationally modified to oxidase (XO), catabolizing mainly purines into uric acid and ROS.

Methods and ResultsWe report that in MFS patients, XOR protein levels increased in aortic samples. In MFS mice (Fbn1C1041G/+), XOR mRNA transcripts augmented in the aorta of young but not older MFS mice. Moreover, the enzymatic activity of XO form increased with respect to XDH. The administration of the XOR inhibitor allopurinol (ALO) blocked the progression of aortic root aneurysm in MFS mice. ALO was also protective when administrated before the appearance of aneurysm. MFS-induced elastic fiber fragmentation, fibrotic remodeling, nuclear translocation of pNRF2, and increased 3-nitrotyrosine levels in the aortic tunica media as well as endothelial dysfunction were all prevented by allopurinol. The large overproduction of H2O2 in MFS aorta did not occur in allopurinol-treated mice whose uric acid plasma levels remained unaltered.

ConclusionsThis study strengthens the concept that redox stress is an important determinant of aortic aneurysm formation and progression in MFS and supports a clinical trial for ALO in the pharmacological treatment of MFS aortopathy.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/464182v3_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@eab02aorg.highwire.dtl.DTLVardef@1d3b39eorg.highwire.dtl.DTLVardef@1b54d0forg.highwire.dtl.DTLVardef@1d61ad0_HPS_FORMAT_FIGEXP  M_FIG C_FIG N-glycan chitobiose core biosynthesis by Agl24 strengthens the hypothesis of an archaeal origin of the eukaryal N-glycosylation Protein N-glycosylation is the most common posttranslational modifications found in all three domains of life. The crenarchaeal N-glycosylation begins with the synthesis of a lipid-linked chitobiose core structure, identical to that in eukaryotes. Here, we report the identification of a thermostable archaeal beta-1,4-N-acetylglucosaminyltransferase, named archaeal glycosylation enzyme 24 (Agl24), responsible for the synthesis of the N-glycan chitobiose core. Biochemical characterization confirmed the function as an inverting {beta}-D-GlcNAc-(1[-&gt;]4)--D-GlcNAc-diphosphodolichol glycosyltransferase. Substitution of a conserved histidine residue, found also in the eukaryotic and bacterial homologs, demonstrated its functional importance for Agl24. Furthermore, bioinformatics and structural modeling revealed strong similarities between Agl24 and both the eukaryotic Alg14/13 and a distant relation to the bacterial MurG, which catalyze the identical or a similar process, respectively. Our data, complemented by phylogenetic analysis of Alg13 and Alg14, revealed similar sequences in Asgardarchaeota, further supporting the hypothesis that the Alg13/14 homologs in eukaryotes have been acquired during eukaryogenesis.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=95 SRC="FIGDIR/small/427365v2_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@13d1ffdorg.highwire.dtl.DTLVardef@1ff48d0org.highwire.dtl.DTLVardef@1056e75org.highwire.dtl.DTLVardef@196ff58_HPS_FORMAT_FIGEXP  M_FIG C_FIG HighlightsO_LIFirst identification and characterization of a thermostable {beta}-D-GlcNAc-(1[-&gt;]4)--D-GlcNAc-diphosphodolichol glycosyltransferase (GT family 28) in Archaea.
C_LIO_LIA highly conserved histidine, within a GGH motif in Agl24, Alg14, and MurG, is essential for function of Agl24.
C_LIO_LIAgl24-like homologs are broadly distributed among Archaea.
C_LIO_LIThe eukaryotic Alg13 and Alg14 are closely related to the Asgard homologs, suggesting their acquisition during eukaryogenesis.
C_LI 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. Maternal microbiota Bifidobacterium promotes placental morphogenesis, nutrient transport and fetal growth in mice The gut microbiota plays a central role in regulating host metabolism. However, while substantial progress has been made in discerning how the microbiota influences host functions post birth and beyond, little is known about how key members of the maternal gut microbiota can influence feto-placental growth. Here, using germ-free and specific-pathogen-free mice, we demonstrate that the bacterium Bifidobacterium breve UCC2003 modulates maternal body adaptations, placental structure and nutrient transporter capacity, with implications for fetal metabolism and growth. The effects of B. breve UCC2003 on feto-placental growth are mediated, in part, by changes in the maternal and placental metabolome (i.e. acetate and carnitine). Histological analysis of the placenta confirmed that Bifidobacterium modifies placental structure via changes in Igf2P0, Dlk1, Mapk1 and Mapk14 expression. Additionally, B. breve UCC2003, acting through Slc2a1 and Fatp1-4 transporters, was shown to restore fetal glycaemia and fetal growth in association with changes in the fetal hepatic transcriptome. This study provides knowledge towards a novel and safe therapeutic strategy for treating pregnancy disorders via modulation of the maternal gut microbiota. Statistical Analysis of Three-Dimensional Chromatin Packing Domains Determined by Chromatin Scanning Transmission Electron Microscopy (ChromSTEM) Chromatin organization over a wide range of length scales plays a critical role in the regulation of transcription and deciphering the interplay of these processes requires high-resolution, three-dimensional, quantitative imaging of chromatin structure in vitro. Herein, we introduce ChromSTEM, a method that utilizes high angle annular dark-field imaging and tomography in scanning transmission electron microscopy combined with DNA-specific staining for electron microscopy. We utilized ChromSTEM to quantify chromatin structure in cultured cells and the statistical packing behavior of the chromatin polymer. Using chromatin mass and density analysis, we observed that chromatin forms spatially well-defined higher-order domains which are around 100 nm in radius, with a radially decreasing mass-density from the center to the periphery. Although the morphological properties of the domains vary within the same cell line, they seem to exhibit greater heterogeneity across cell lines, underlying a potential role of statistical chromatin packing in regulating cell-type-specific gene expression. HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. Exploring the natural origins of SARS-CoV-2 in the light of recombination The lack of an identifiable intermediate host species for the proximal animal ancestor of SARS-CoV-2, and the large geographical distance between Wuhan and where the closest evolutionary related coronaviruses circulating in horseshoe bats (Sarbecoviruses) have been identified, is fuelling speculation on the natural origins of SARS-CoV-2. We have comprehensively analysed phylogenetic relations between SARS-CoV-2, and the related bat and pangolin Sarbecoviruses sampled so far. Determining the likely recombination events reveals a highly reticulate evolutionary history within this group of coronaviruses. Clustering of the inferred recombination events is non-random with evidence that Spike, the main target for humoral immunity, is beside a recombination hotspot likely driving antigenic shift in the ancestry of bat Sarbecoviruses. Coupled with the geographic ranges of their hosts and the sampling locations, across southern China, and into Southeast Asia, we confirm horseshoe bats, Rhinolophus, are the likely SARS-CoV-2 progenitor reservoir species. By tracing the recombinant sequence patterns, we conclude that there has been relatively recent geographic movement and co-circulation of these viruses ancestors, extending across their bat host ranges in China and Southeast Asia over the last 100 years or so. We confirm that a direct proximal ancestor to SARS-CoV-2 is yet to be sampled, since the closest relative shared a common ancestor with SARS-CoV-2 approximately 40 years ago. Our analysis highlights the need for more wildlife sampling to (i) pinpoint the exact origins of SARS-CoV-2s animal progenitor, and (ii) survey the extent of the diversity in the related Sarbecoviruses phylogeny that present high risk for future spillover.

HighlightsO_LIThe origin of SARS-CoV-2 can be traced to horseshoe bats, genus Rhinolophus, with ranges in both China and Southeast Asia.
C_LIO_LIThe closest known relatives of SARS-CoV-2 exhibit frequent transmission among their Rhinolophus host species.
C_LIO_LISarbecoviruses have undergone extensive recombination throughout their evolutionary history.
C_LIO_LIAccounting for the mosaic patterns of these recombinants is important when inferring relatedness to SARS-CoV-2.
C_LIO_LIBreakpoint patterns are consistent with recombination hotspots in the coronavirus genome, particularly upstream of the pike open reading frame with a coldspot in S1.
C_LI Outcome-selective reinstatement is predominantly context-independent, and associated with c-Fos activation in the posterior dorsomedial striatum Research from human and animal studies has found that responding that has been successfully reduced following treatment can return upon exposure to certain contexts. An individual in recovery from alcohol use disorder, for example, might relapse to drinking upon visiting their favourite bar. However, most of these data have been derived from experiments involving a single (active) response, and the context-dependence of returned responding in situations involving choice between multiple actions and outcomes is less well-understood. We thus investigated how outcome-selective reinstatement - a procedure involving choice between multiple actions - was affected by altering the physical context in rats. In Experiment 1, rats were trained over 6 days to press a left lever for one food outcome (pellets or sucrose) and a right lever for the other outcome. Then, rats received an extinction session in either the same context (A) as lever press training, or in a different context (B). Rats were tested immediately (5 minutes) after extinction in Context A or B such that there were four groups in total: AAA, ABB, ABA, and AAB. Reinstatement testing consisted of one food outcome being delivered  freely (i.e. unearned by lever pressing and unsignalled by cues) to the food magazine every 4 minutes in the following order: Sucrose, Pellet, Pellet, Sucrose. Selective reinstatement was considered intact if pellet delivery increased pressing selectively on the pellet lever, and sucrose delivery selectively increased pressing on the sucrose lever. This result (Reinstated > Nonreinstated) was observed for rats in group AAA and ABB, but not rats in groups ABA and AAB. Experiment 2 was conducted identically, except that rats received two extinction sessions over two days and tested one day later. This time, all groups demonstrated intact outcome-selective reinstatement regardless of context. Analysis of c-Fos expression in several brain regions revealed that only c-Fos expression in the posterior dorsomedial striatum (pDMS) was related to intact reinstatement performance. Overall, these results suggest that outcome-selective reinstatement is predominantly context-independent, and that intact reinstatement is related to neuronal activity in the pDMS.

HighlightsO_LIOutcome-selective reinstatement is predominantly context-independent
C_LIO_LIOutcome-selective reinstatement is entirely context-independent after multiple extinction sessions
C_LIO_LIOutcome-selective reinstatement increases c-Fos expression in dorsomedial striatum
C_LIO_LIc-Fos expression in orbitofrontal cortex and dorsal hippocampus is unaffected by selective reinstatement.
C_LI Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. 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. Slowing down an analyte in a nanopore An unresolved problem in nanopore sensing is the high translocation speed ([~]10-100 monomers/s) of an analyte (nucleotide, DNA, amino acid (AA), peptide) through the pore. Here a method based on reversing the pore voltage and changing the solution pH is described. A simplified Fokker-Planck model shows mean translocation times of 1-10 ms in a nanopore of length 10 nm. Simulations show that a positive-negative voltage profile can trap an analyte for [~]1 ms. This method can be used for free nucleotides, single AAs, oligonucleotides, and oligopeptides. Its applicability to existing nanopore sensing and sequencing techniques and implementation issues are discussed. Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.) Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Significance StatementUsing traditional and automated/high-throughput (using a novel deep learning approach) phenotyping methods we studied leaf anatomical variation in sunflower. Genome-wide association (GWA) analyses identified numerous genomic regions underlying individual trait variation and regions underlying major multivariate axes of phenotypic variation. These results illustrate the value of employing a multivariate approach to GWA analyses and shed light on the extent to which leaf trait (co-)variation can be genetically decoupled to explore novel phenotypic space. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development.