Earliest photic zone niches probed by ancestral microbial rhodopsins For billions of years, life has continuously adapted to dynamic physical conditions near the Earths surface. Fossils and other preserved biosignatures in the paleontological record are the most direct evidence for reconstructing the broad historical contours of this adaptive interplay. However, biosignatures dating to Earths earliest history are exceedingly rare. Here, we combine phylogenetic inference of primordial rhodopsin proteins with modeled spectral features of the Precambrian Earth environment to reconstruct the paleobiological history of this essential family of photoactive transmembrane proteins. Our results suggest that ancestral microbial rhodopsins likely acted as light-driven proton pumps and were spectrally tuned toward the absorption of green light, which would have enabled their hosts to occupy depths in a water column or biofilm where UV wavelengths were attenuated. Subsequent diversification of rhodopsin functions and peak absorption frequencies was enabled by the expansion of surface ecological niches induced by the accumulation of atmospheric oxygen. Inferred ancestors retain distinct associations between extant functions and peak absorption frequencies. Our findings suggest that novel information encoded by biomolecules can be used as "paleosensors" for conditions of ancient, inhabited niches of host organisms not represented elsewhere in the paleontological record. The coupling of functional diversification and spectral tuning of this taxonomically diverse protein family underscores the utility of rhodopsins as universal testbeds for inferring remotely detectable biosignatures on inhabited planetary bodies.
Evolution of phenotypic variance provides insights into the genetic basis of adaptation Most traits are polygenic and the contributing loci can be identified by GWAS. Their adaptive architecture is, however, difficult to characterize. Here, we propose to study the adaptive architecture of traits by monitoring the evolution of their phenotypic variance during adaptation to a new environment in well-defined laboratory conditions. Extensive computer simulations show that the evolution of phenotypic variance in a replicated experimental evolution setting provides power to distinguish between oligogenic and polygenic adaptive architectures. We compared gene expression variance in male Drosophila simulans before and after 100 generations of adaptation to a novel hot environment. The variance change in gene expression was indistinguishable for genes with and without a significant change in mean expression after 100 generations of evolution. We suggest that a majority of adaptive gene expression evolution can be explained by a polygenic architecture. We propose that tracking the evolution of phenotypic variance across generations can provide an approach to characterize the adaptive architecture.
Reset Networks: Emergent Topography by Composition of Convolutional Neural Networks We introduce the Reset model, a composition of neural networks - typically several levels of convolutional neural networks - whose outputs at one level are gathered and reshaped into a spatial input for the next level. We show that units in Reset networks self-organize into clusters when trained on MNIST, Fashion MNIST, CIFAR-10 and CIFAR-100. We then show that a stronger type of self-organization, reminiscent of the topography found for numbers in parietal cortex, arises when number images are mapped onto developmentally realistic number codes. We outline the implications of this model for theories of the cortex and developmental neuroscience. Live Analysis of Position Effect Variegation (PEV) in Drosophila Reveals Different Modes of Action for HP1a and Su(var)3-9 AO_SCPLOWBSTRACTC_SCPLOWPosition Effect Variegation (PEV) results from the juxtaposition of euchromatic and heterochromatic components of eukaryotic genomes, silencing genes near the new euchromatin/heterochromatin junctions. Silencing is itself heritable through S phase, giving rise to distinctive random patterns of cell clones expressing the genes intermixed with clones in which the genes are silenced. Much of what we know about epigenetic inheritance in the soma stems from work on PEV aimed at identifying the components of the silencing machinery and its mechanism of inheritance. The roles of two central gene activities - the Su(var)3-9-encoded histone H3-Lysine-9 methyltransferase and the Su(var)205-encoded methyl-H3-Lysine-9 binding protein HP1a - have been inferred from terminal phenotypes, leaving considerable gaps in understanding how PEV behaves through development. Here, we investigate the PEV phenotypes of Su(var)3-9 and Su(var)205 mutations in live developing tissues. We discovered that mutation in Su(var)205 compromises the initial establishment of PEV in early embryogenesis. Later gains of heterochromatin-induced gene silencing are possible, but are unstable and lost rapidly. In contrast, a strain with mutation in Su(var)3-9 exhibits robust silencing early in development, but fails to maintain it through subsequent cell divisions. Our analyses show that while the terminal phenotypes of these mutations may appear identical, they have arrived at them through different developmental trajectories. We discuss how our findings expand and clarify existing models for epigenetic inheritance of heterochromatin-induced gene silencing.

SO_SCPLOWIGNIFICANCEC_SCPLOWCurrent concepts of epigenetic inheritance are exemplified by Position Effect Variegation, a phenomenon whereby heterochromatin can repress genes in clonal cell lineages. Heterochromatin is required for genome protection as it silences toxic transposable elements and prevents instability of repeat sequences. Histone H3 modified by methylation of Lysine-9 and HP1a are critical components of heterochromatin. Using live cell analysis of PEV in mutants in strains with mutations in Su(var)3-9, which encodes the histone methyltransferase, and Su(var)205, which encodes HP1a, we describe an unexpected dynamism in PEV, challenging current models of epigenetics, and revealing unexpectedly different modes of action of these two fundamental components of heterochromatin. The Hippocampus uses cryptographic principles rather than plasticity to enable animals to secretly cache and retrieve their food For animals, the ability to hide and retrieve valuable information, such as the location of food, can mean the difference between life and death. Here, we propose that to achieve this, their brain uses spatial cells similarly to how we utilize encryption for data security. Some animals are able to cache hundreds of thousands of food items annually by each individual and later retrieve most of what they themselves stashed. Rather than memorizing their cache locations as previously suggested, we propose that they use a single cryptographic-like mechanism during both caching and retrieval. The model we developed is based on hippocampal spatial cells, which respond to an animals positional attention, such as when the animal enters a specific region (place-cells) or gazes at a particular location (spatial-view-cells). We know that the region that activates each spatial cell remains consistent across subsequent visits to the same area but not between areas. This remapping, combined with the uniqueness of cognitive maps, produces a persistent crypto-hash function for both food caching and retrieval. We also show that the model stores temporal information that helps animals in food caching order preference as previously observed. This behavior, which we refer to as crypto-taxis, might also explain consistent differences in decision-making when animals are faced with a large number of alternatives such as in foraging. Modular evolution of secretion systems and virulence plasmids in a bacterial species complex BackgroundMany bacterial taxa are species complexes and uncertainties regarding the organization of their genetic diversity challenge research efforts. We utilized Agrobacterium tumefaciens, a taxon known for its phytopathogenicity and applications in transformation, as a study system and devised strategies for investigating genome diversity and evolution of species complexes.

ResultsWe utilized 35 genome assemblies to achieve a comprehensive and balanced sampling of A. tumefaciens. Our confident inference of gene content and core-genome phylogeny supported a quantitative guideline for delineating 12 species and allowed for robust investigations of genes critical in fitness and ecology. For the type VI secretion system (T6SS) involved in interbacterial competition and thought to be conserved, we detected multiple losses and one horizontal gene transfer. For the tumor-inducing plasmids (pTi) and pTi-encoded type IV secretion system (T4SS) that are essential for agrobacterial phytopathogenicity, we uncovered novel diversity and hypothesized their involvement in shaping this species complex. Intriguingly, for both T6SS and T4SS, genes encoding structural components are highly conserved, whereas extensive diversity exists for genes encoding effectors and other proteins.

ConclusionsWe demonstrated that the combination of a phylogeny-guided sampling scheme and an emphasis on high-quality assemblies provides a cost-effective approach for robust analysis in evolutionary genomics. Our strategies for multi-level investigations at scales that range from whole-genomes to intragenic domains and phylogenetic depths of between- and within-species are applicable to other bacteria. Finally, modularity observed in the molecular evolution of genes and domains is useful for inferring functional constraints and informing experimental works. Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix PurposeTo overcome some limitations of prior orientation-dependent proton transverse relaxation formalisms in white matter (WM) with a novel framework based on the generalized magic angle effect function.

MethodsA cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in human brain WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R2 and [Formula] were divided into isotropic [Formula] and anisotropic parts, [Formula], with  denoting an open angle and{varepsilon} 0 an orientation ({Phi}) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared with prior models without{varepsilon} 0 on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P [&le;] 0.05.

ResultsFit A significantly (P<0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted{varepsilon} 0 was in good agreement with calculated{varepsilon} 0 from directional diffusivities. Compared with those from healthy adult, the fitted [Formula], and  from neonates were substantially reduced but{varepsilon} 0 increased, consistent with incomplete myelination. Significant positive [Formula] and negative ( and [Formula]) correlations were found with aging (demyelination) in elderly.

ConclusionThe developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in human brain WM, shedding new light on myelin microstructural alterations at the molecular level. Genetic Impairment of Succinate Metabolism Disrupts Bioenergetic Sensing in Adrenal Neuroendocrine Cancer Metabolic dysfunction mutations can impair energy sensing and cause cancer. Loss of function of mitochondrial TCA cycle enzyme, succinate dehydrogenase B (SDHB) results in various forms of cancer typified by pheochromocytoma (PC). Here we delineate a signaling cascade where the loss of SDHB induces the Warburg effect in PC tumors, triggers dysregulation of Ca2+ homeostasis, and aberrantly activates calpain and the protein kinase Cdk5, through conversion of its cofactor from p35 to p25. Consequently, aberrant Cdk5 initiates a cascade of phospho- signaling where GSK3 inhibition inactivates energy sensing by AMP-kinase through dephosphorylation of the AMP-kinase {gamma} subunit, PRKAG2. Overexpression of p25-GFP in mouse adrenal chromaffin cells also elicits this phosphorylation signaling and causes PC tumor formation. A novel Cdk5 inhibitor, MRT3-007, reversed this phospho-cascade, invoking an anti- Warburg effect, cell cycle arrest, and senescence-like phenotype. This therapeutic approach halted tumor progression in vivo. Thus, we reveal an important novel mechanistic feature of metabolic sensing and demonstrate that its dysregulation underlies tumor progression in PC and likely other cancers.

HighlightsO_LILoss of SDHB function in pheochromocytoma causes Ca2+ dysregulation, calpain activation, and aberrant activation of the protein kinase Cdk5.
C_LIO_LIHyperactive Cdk5 deregulates a GSK3/PRKAG2/AMPK signaling cascade.
C_LIO_LIp25 overexpression and consequent aberrant Cdk5 activity in chromaffin cells causes pheochromocytoma.
C_LIO_LIInhibition of Cdk5 activates the PRKAG2/AMPK/p53 signaling to rescue cell senescence and block PC tumor progression.
C_LI

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=175 SRC="FIGDIR/small/475410v2_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@16f0a48org.highwire.dtl.DTLVardef@10c1c7aorg.highwire.dtl.DTLVardef@21b100org.highwire.dtl.DTLVardef@1eb8674_HPS_FORMAT_FIGEXP  M_FIG C_FIG Mathematical Modeling and Analysis of Mitochondrial Retrograde Signaling Dynamics Mitochondria, semi-autonomous eukaryotic organelles, participate in energy production and metabolism, making mitochondrial quality control crucial. As most mitochondrial proteins are encoded by nuclear genes, quality control depends on proper mitochondria-nucleus communication, designated mitochondrial retrograde signaling. Early studies focused on retrograde signaling participants and specific gene knockouts. However, mitochondrial signal modulation remains elusive. Using yeast, we simulated signal propagation following mitochondrial damage and proposed a mathematical model based ordinary differential equations. Mitochondrial retrograde signaling decisions were described using a Boolean model. Dynamics were analyzed through an ordinary differential equation-based model and extended to evaluate the model response to noisy damage signals. Simulation revealed localized protein concentration dynamics, including waveforms, frequency response, and robustness under noise. Retrograde signaling is bistable with three localized steady states, and increased damage compromises robustness. We elucidated mitochondrial retrograde signaling, thus providing a basis for drug design against yeast and fungi.

Author SummaryThe yeast RTG pathway regulates mitochondrial metabolism and mitochondrial quality through passing mitochondrial signal to the nucleus to modulate gene expressions. Using microscopic data of RTG proteins translocation from cytosol to nucleus, the parameters were found by fitting 16 knockout conditions simulated from the proposed differential equation-based model. Further, through dose response, ultrasensitivity, frequency response and noise interference, we demonstrated a switch-like property of RTG pathway activation, a capacity charging property of RTG protein translocation, a low pass filter property of signals in different frequencies, and noise amplification due to mitochondrial signal. The proposed mathematical model enables us to understand the dynamics and mechanisms in mitochondrial retrograde signaling in yeast and provides potential antifungal treatment strategies. Electrical Stimulated GLUT4 Signaling Attenuates Critical Illness-Associated Muscle Wasting BackgroundCritical illness myopathy (CIM) is a debilitating condition characterized by the preferential loss of the motor protein myosin. CIM is a byproduct of critical care, attributed to impaired recovery, long-term complications, and mortality. CIM pathophysiology is complex, heterogeneous and remains incompletely understood, however loss of mechanical stimuli contributes to critical illness associated muscle atrophy and weakness. Passive mechanical loading (ML) and electrical stimulation (ES) therapies augment muscle mass and function. While having beneficial outcomes, the mechanistic underpinning of these therapies is less known. Therefore, here we aimed to assess the mechanism by which chronic supramaximal ES ameliorates CIM in a unique experimental rat model of critical care.

MethodsRats were subjected to 8 days critical care conditions entailing deep sedation, controlled mechanical ventilation, and immobilization with and without direct soleus ES. Muscle size and function were assessed at the single cell level. RNAseq and Western blotting were employed to understand the mechanisms driving ES muscle outcomes in CIM.

ResultsFollowing 8 days of controlled mechanical ventilation and immobilization, soleus muscle mass, Myosin:Actin ratio and single muscle fiber maximum force normalized to cross-sectional area (specific force) were reduced by 40-50% (p< 0.0001). ES significantly reduced the loss of soleus muscle fiber cross-sectional area (CSA) and Myosin:Actin ratio by approximately 30% (p< 0.05) yet failed to effect specific force. RNAseq pathway analysis revealed downregulation of insulin signaling in the soleus muscle following critical care and GLUT4 trafficking was reduced by 55% leading to an 85% reduction of muscle glycogen content (p< 0.01). ES promoted phosphofructokinase and insulin signaling pathways to control levels (p< 0.05), consistent with the maintenance of GLUT4 translocation and glycogen levels. AMPK, but not AKT, signaling pathway was stimulated following ES, where the downstream target TBC1D4 increased 3 logFC (p= 0.029) and AMPK-specific P-TBC1D4 levels were increased approximately 2-fold (p= 0.06). Reduction of muscle protein degradation rather than protein synthesis promoted soleus CSA, as ES reduced E3 ubiquitin proteins, Atrogin-1 (p= 0.006) and MuRF1 (p= 0.08) by approximately 50%, downstream of AMPK-FoxO3.

ConclusionsES maintained GLUT4 translocation through increased AMPK-TBC1D4 signaling leading to improved muscle glucose homeostasis. Soleus CSA and myosin content was promoted through reduced protein degradation via AMPK-FoxO3 E3 ligases, Atrogin-1 and MuRF1. These results demonstrate chronic supramaximal ES reduces critical care associated muscle wasting, preserved glucose signaling and reduced muscle protein degradation in CIM. Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. AICAR-mediated selection of karyotypically normal cells in Pallister-Killian syndrome patient-derived skin fibroblasts Pallister-Killian Syndrome (PKS) is a mosaic aneuploidy syndrome typically caused by the presence of a supernumerary marker isochromosome containing two extra copies of the short arm of chromosome 12 (iso-chromosome 12p or tetrasomy 12p). This isochromosome is always present in a mosaicism state that has tissue limited distribution. PKS is characterized by severe neurodevelopmental delay, intellectual disability, multisystem involvement and congenital malformations including typical dysmorphic features and skin pigmentation anomalies.

Aneuploid cells, irrespective of the identity of the supernumerary chromosome, including cancer cells, yeast cells and mouse embryonic fibroblasts (MEFs), have been demonstrated to present a disruption of protein homeostasis and increased basal stress levels; resulting in a greater sensitivity to chemical compounds inducing cellular energy stress compared to euploid cell lines. The burden of trisomy 21 has also been recently shown to impair the proteostasis network in lymphoblastoid cell lines and fibroblasts obtained from individuals with Down syndrome.

In this study, we demonstrate that AICAR, 5-Aminoimidazole-4-carboxamide ribonucleotide, a known energy stress inducing drug with antiproliferative effects on aneuploidy cancer cells and MEFs, is also able to selectively eliminate cells carrying the isochromosme12p in PKS clones in a time and dosage dependent manner. Collectively, our results indirectly provide evidence of increased basal energy and proteotoxic stress in PKS cells carrying isochromosome 12p, and suggest a potential therapeutic drug-based strategy that, selectively acting as a stressor for aneuploid cells, may establish the euploid state in PKS and a broader spectrum of human mosaic disorders. Ultrahigh-Throughput Ambient MS: Direct Analysis at 22 Samples per Second by Infrared Matrix-Assisted Laser Desorption Electrospray Ionization Mass Spectrometry Infrared Matrix-Assisted Laser Desorption Electrospray Ionization (IR-MALDESI) mass spectrometry is an ambient-direct sampling method being developed for high-throughput, label-free, biochemical screening of large-scale compound libraries. Here, we report the development of an ultrahigh-throughput continuous motion IR-MALDESI sampling approach capable of acquiring data at rates up to 22.7 samples per second in a 384-well microtiter plate. At top speed, less than 1% analyte carryover is observed from well-to-well and signal intensity relative standard deviations (RSD) of 11.5% and 20.9% for 3 M 1-hydroxymidazolam and 12 M dextrorphan, respectively, are achieved. The ability to perform parallel kinetics studies on 384 samples with ~30s time resolution using an isocitrate dehydrogenase 1 (IDH1) enzyme assay is shown. Finally, we demonstrate the repeatability and throughput of our approach by measuring 115,200 samples from 300 microtiter plate reads consecutively over 5.54 hours with RSDs under 8.14% for each freshly introduced plate. Taken together, these results demonstrate the use of IR-MALDESI at sample acquisition rates that surpass other currently reported direct sampling mass spectrometry approaches used for high throughput compound screening.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/465730v3_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@5f916corg.highwire.dtl.DTLVardef@4473b4org.highwire.dtl.DTLVardef@82e8forg.highwire.dtl.DTLVardef@128b57b_HPS_FORMAT_FIGEXP  M_FIG C_FIG The combined treatment of Molnupiravir and Favipiravir results in a marked potentiation of efficacy in a SARS-CoV2 hamster infection model through an increased frequency of mutations in the viral genome. Favipiravir and Molnupiravir, orally available antivirals, have been reported to exert antiviral activity against SARS-CoV2. In recent days preliminary efficacy data have been reported in COVID-19 patients. We here studied the combined antiviral effect of the drugs in the SARS-CoV2 hamster infection model. We first demonstrate that Molnupiravir can reduce infectious virus titers in lungs of infected animals in a dose-dependent manner by up to 3.5 log10 which is associated with a marked improvement of virus-induced lung pathology. When animals are treated with a combination of suboptimal doses of Molnupiravir and Favipiravir (that each alone result in respectively a 1.3 log10 and 1.1 log10 reduction of infectious virus titers in the lungs), a marked combined potency is observed. Infectious virus titers in the lungs of animals treated with the combo are on average reduced by 4.5 log10 and infectious virus are no longer detected in the lungs of 60% of treated infected animals. Both drugs result in an increased mutation frequency of the remaining viral RNA recovered from the lungs. In the combo-treated hamsters an increased frequency of C-to-T and G-to-A mutations in the viral RNA is observed as compared to the single treatment groups which may explain the pronounced antiviral potency of the combination. Our findings may lay the basis for the design of clinical studies to test the efficacy of the combination of Molnupiravir and Favipiravir in the treatment of COVID-19. Relationship between shear wave velocity and muscle activation is inconsistent across different muscle types There is an increasing use of shear wave ultrasound elastography to quantify mechanical properties of muscles under various conditions such as changes muscle length and levels of activation in healthy and pathological muscle. However, little is known about the variability in shear wave velocity among muscles as most studies investigate one specific muscle. The purpose of this study was to determine if the relationship between SWV and muscle activation is consistent across muscles with different architectures: biceps brachii, tibialis anterior, and medial gastrocnemius, All measures were made at matching levels of activation and approximately at the optimal length for each muscle to control for length-dependent changes in the relationship between activation and force or stiffness. We also conducted a control experiment to determine how the passive force within a muscle alters the relationship between muscle activation and shear wave velocity. The relationship between shear wave velocity-squared and activation above 10% MVC differed across muscles, with biceps brachii and medial gastrocnemius showing a lower slope than tibialis anterior. Shear wave velocity-squared also differed between muscles at the shortest length (p<0.001) and the increase in shear wave velocity-squared with muscle lengthening also differed between muscle types (p = 0.005) Muscle-specific differences could not be explained by the architectural features such as pennation angle, fiber length, and physiological cross-sectional area. Our results demonstrate that there is not a unique relationship between muscle activation and shear wave velocity, highlighting the importance of understanding the many factors contributing to shear wave propagation in muscle before elastography can be used to make quantitative comparisons across muscle types. Honeycomb crystallography: comb formation under geometric frustrations 1As honeybees build their nests in pre-existing tree cavities, they must deal with dealing with the presence of geometric constraints, resulting in non-regular hexagons and topological defects in the comb. In this work, we study how bees adapt to their environment in order to regulate the comb structure. Specifically, we identify the irregularities in honeycomb structure in the presence of various geometric frustrations. We 3D-print experimental frames with a variety of constraints imposed on the imprinted foundations. The combs constructed by the bees show clear evidence of reoccurring patterns built by bees in response to specific geometric frustrations on these starter frames. Furthermore, using an experimental-modeling framework, we demonstrate that these patterns can be successfully modeled and replicated through a simulated annealing process, in which the minimized potential is a variation of the Lennard-Jones potential that only considers first-neighbor interactions according to a Delaunay triangulation. Our simulation results not only confirm the connection between honeycomb structures and other crystal systems such as graphene, but also show that irregularities in the honeycomb structure can be explained as the result of local interactions between honeybees and their immediate surroundings, leading to emergent global order. Additionally, our computational model can be used to describe specific strategies that bees use to effectively solve each geometric mismatch problem while minimizing cost of comb building. Conserved organ-specific microbial assemblages in different populations of a terrestrial crab Brachyuran crabs originated in the oceans and evolved specific morphological and physiological adaptations to live in freshwater, intertidal and even terrestrial habitats but the role of a selection mechanism involving symbiotic microorganisms long these colonization processes are not known. In this work we investigated the associated microbiota of three populations of a terrestrial brachyuran crab, Chiromantes haematocheir, to find evidence of a conserved crab-specific microbiome unrelated to the population of origin and dissimilar from environmental microbial assemblages. Bacterial 16S rRNA gene and fungal ITS sequences were obtained from selected crab organs and environmental matrices to profile microbial communities. In spite of the presence of truly marine larval stages and no gregarious behaviour, favouring microbiota exchanges, we found common, organ-specific microbiota, associated to the gut and the gills of the crabs (with more than 15% of the genera detected specifically enriched only in one organ). Our results suggest an early establishment of a new common, stable microbiota in the transition from water to land. The RAG1 Ubiquitin Ligase Domain Enhances Assembly and Selection of T Cell Receptor Genes to Restrain the Autoimmune Hazard of Generating T Cell Receptor Diversity RAG1/RAG2 (RAG) endonuclease-mediated assembly of diverse lymphocyte antigen receptor genes by V(D)J recombination is critical for the development and immune function of T and B cells. However, this process creates highly self-reactive cells that must be properly selected to suppress autoimmunity. The RAG1 protein contains a ubiquitin ligase domain that stabilizes RAG1 and stimulates RAG endonuclease activity in vitro. We report that mice with a mutation that inactivates the RAG1 ubiquitin ligase in vitro exhibit modestly reduced thymic cellularity, decreased assembly and altered repertoires of T cell receptor (TCR) {beta} and  genes in thymocytes, and impaired thymocyte developmental transitions that require the assembly of TCR{beta} or  genes and signaling by their proteins. These RAG1 mutant mice also exhibit less efficient positive selection and superantigen-mediated negative selection of conventional {beta} T cells, 2) impaired differentiation of iNKT lineage {beta} T cells, and 3) CD4+ {beta} T cells with elevated autoimmune potential. Our findings demonstrate that the RAG1 ubiquitin ligase domain functions in vivo to stimulate the assembly and selection of TCR{beta} and TCR genes, thereby establishing replete diversity of {beta} TCRs and {beta} T cell lineages while restraining the inherent autoimmune hazard of generating diverse antigen specificities. Metacells untangle large and complex single-cell transcriptome networks The exponential scaling of scRNA-seq data represents an important hurdle for downstream analyses. Here we develop a coarse-graining framework where highly similar cells are merged into metacells. We demonstrate that metacells not only preserve but often improve the results of downstream analyses including visualization, clustering, differential expression, cell type annotation, gene correlation, imputation, RNA velocity and data integration. By capitalizing on the redundancy inherent to scRNA-seq data, metacells significantly facilitate and accelerate the construction and interpretation of single-cell atlases, as demonstrated by the integration of 1.46 million cells from COVID-19 patients in less than two hours on a standard desktop. A multiscale lipid and cellular atlas of the human kidney The human kidney is composed of over 26 cell types that actively coordinate with each other to form higher-order structures, such as the nephron. It is not yet understood how these structures vary throughout a single organ or amongst the same organs within the human population. We have developed an extensive lipid and cellular atlas of the human kidney consisting of over 3 million cells comprising 75,000 functional tissue units (i.e., glomeruli, proximal tubules, distal tubules, and collecting ducts) from 13 human subjects. This atlas was developed using spatially registered and integrated technologies consisting of imaging mass spectrometry, multiplexed immunofluorescence, stained microscopy, and autofluorescence microscopy, to comprehensively probe large (i.e., centimeter-sized) areas of tissue. The cellular organization and lipid profiles of glomeruli, proximal tubules, distal tubules, and collecting ducts were discovered through these multimodal imaging data as well as their intra- and inter-subject variance. Relating the lipid profiles obtained from imaging mass spectrometry to distinct cell types obtained from immunofluorescence allowed us to hypothesize the functional role of specific phospholipids that have not previously been described. These hypotheses include subject characteristics, such as BMI and sex. The integrated data from the aforementioned datasets provide a valuable reference for kidney researchers, are publicly available through the NIH Human Biomolecular Atlas Program (https://portal.hubmapconsortium.org/), and discussed below. 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. Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia BackgroundIn humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals.

DataWe obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs.

ResultsIn human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data.

ConclusionsOur results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs. 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. Water levels primarily drive variation in photosynthesis and nutrient use of scrub Red Mangroves in the southeastern Florida Everglades We investigated how mangrove-island micro-elevation (i.e., habitat: center vs. edge) affects tree physiology in a scrub mangrove forest of the southeastern Everglades. We measured leaf gas exchange rates of scrub Rhizophora mangle L. trees monthly during 2019, hypothesizing that CO2 assimilation (Anet) and stomatal conductance (gsw) would decline with increasing water levels and salinity, expecting more-considerable differences at mangrove-island edges than centers, where physiological stress is greatest. Water levels varied between 0 and 60 cm from the soil surface, rising during the wet season (May-October) relative to the dry season (November-April). Porewater salinity ranged from 15 to 30 ppt, being higher at mangrove-island edges than centers. Anet maximized at 15.1 {micro}mol m-2 s-1, and gsw was typically <0.2 mol m-2 s-1, both of which were greater in the dry than the wet season and greater at island centers than edges, with seasonal variability being roughly equal to variation between habitats. After accounting for season and habitat, water level positively affected Anet in both seasons but did not affect gsw. Our findings suggest that inundation stress (i.e., water level) is the primary driver of variation in leaf gas exchange rates of scrub mangroves in the Florida Everglades, while also constraining Anet more than gsw. The interaction between inundation stress due to permanent flooding and habitat varies with season as physiological stress is alleviated at higher-elevation mangrove-island center habitats during the dry season. Freshwater inflows during the wet season, increase water levels and inundation stress at higher-elevation mangrove-island centers, but also potentially alleviate salt and sulfide stress in soils. Thus, habitat heterogeneity leads to differences in nutrient and water acquisition and use between trees growing in island centers versus edges, creating distinct physiological controls on photosynthesis, which likely affect carbon flux dynamics of scrub mangroves in the Everglades. 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. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice.