Standardised method for cardiomyocyte isolation and purification from individual murine neonatal, infant, and adult hearts Primary cardiomyocytes are invaluable for understanding postnatal heart development. However, a universal method to obtain freshly purified cardiomyocytes without using different agedependent isolation procedures and cell culture, is lacking. Here, we report the development of a standardised method that allows rapid isolation and purification of high-quality cardiomyocytes from individual neonatal through to adult C57BL/6J murine hearts. Langendorff retrograde perfusion, which is currently limited to adult hearts, was adapted for use in neonatal and infant hearts by developing an easier in situ aortic cannulation technique. Tissue digestion conditions were optimised to achieve efficient digestion of hearts of all ages in a comparable timeframe (<14 min). This resulted in a high yield (1.56-2.2x106 cells/heart) and viability (~70-100%) of cardiomyocytes post-isolation. An immunomagnetic cell separation step was then applied to yield highly purified cardiomyocytes (~95%) as confirmed by immunocytochemistry, flow cytometry, and qRT-PCR. For cell-type specific studies, cardiomyocyte DNA, RNA, and protein could be extracted in sufficient yields to conduct molecular experiments. We generated transcriptomic datasets for neonatal cardiomyocytes from individual hearts, for the first time, which revealed nine sex-specific genes (FDR<0.05) encoded on the sex chromosomes. Finally, we also developed an in situ fixation protocol that preserved the native cytoarchitecture of cardiomyocytes (~94% rodshaped post-isolation), and used it to evaluate cell morphology during cardiomyocyte maturation, as well as capture spindle-shaped neonatal cells undergoing cytokinesis. Together, these procedures allow molecular and morphological profiling of high-quality cardiomyocytes from individual hearts of any postnatal age.
Surface-associated antigen induces B-cell permeabilization and lysosome exocytosis facilitating antigen uptake and presentation to T-cells B-cell receptor (BCR)-mediated antigen internalization and presentation are essential for humoral memory immune responses. Antigen encountered by B-cells is often tightly associated with the surface of pathogens and/or antigen-presenting cells. Internalization of such antigens requires myosin-mediated traction forces and extracellular release of lysosomal enzymes, but the mechanism triggering lysosomal exocytosis is unknown. Here we show that BCR-mediated recognition of antigen tethered to beads, to planar lipid-bilayers or expressed on cell surfaces causes localized plasma membrane (PM) permeabilization, a process that requires BCR signaling and non-muscle myosin II activity. B-cell permeabilization triggers PM repair responses involving lysosomal exocytosis, and B-cells permeabilized by surface-associated antigen internalize more antigen than cells that remain intact. Higher affinity antigens cause more B-cell permeabilization and lysosomal exocytosis and are more efficiently presented to T-cells. Thus, PM permeabilization by surface-associated antigen triggers a lysosome-mediated B-cell resealing response, providing the extracellular hydrolases that facilitate antigen internalization and presentation.
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. R7 photoreceptor axon targeting depends on the relative levels of lost and found expression in R7 and its synaptic partners As neural circuits form, growing processes select the correct synaptic partners through interactions between cell surface proteins. The presence of such proteins on two neuronal processes may lead to either adhesion or repulsion; however, the consequences of mismatched expression have rarely been explored. Here we show that the Drosophila CUB-LDL protein Lost and found (Loaf) is required in the UV-sensitive R7 photoreceptor for normal axon targeting only when Loaf is also present in its synaptic partners. Although targeting occurs normally in loaf mutant animals, removing loaf from photoreceptors or expressing it in their postsynaptic neurons Tm5a/b or Dm9 in a loaf mutant causes mistargeting of R7 axons. Loaf localizes primarily to intracellular vesicles including endosomes. We propose that Loaf regulates the trafficking or function of one or more cell surface proteins, and an excess of these proteins on the synaptic partners of R7 prevents the formation of stable connections. swCAM: estimation of subtype-specific expressions in individual samples with unsupervised sample-wise deconvolution MotivationComplex biological tissues are often a heterogeneous mixture of several molecularly distinct cell or tissue subtypes. Both subtype compositions and expressions in individual samples can vary across different biological states or conditions. Computational deconvolution aims to dissect patterns of bulk gene expression data into subtype compositions and subtype-specific expressions. Typically, existing deconvolution methods can only estimate averaged subtype-specific expressions in a population, while detecting differential expressions or co-expression networks in particular subtypes requires unique subtype expression estimates in individual samples. Different from population-level deconvolution, however, individual-level deconvolution is mathematically an underdetermined problem because there are more variables than observations.

ResultsWe report a sample-wise Convex Analysis of Mixtures (swCAM) method that can estimate subtype proportions and subtype-specific expressions in individual samples from bulk tissue transcriptomes. We extend our previous CAM framework to include a new term accounting for between-sample variations and formulate swCAM as a nuclear-norm and{ell} 2,1-norm regularized matrix factorization problem. We determine hyperparameter values using a cross-validation scheme with random entry exclusion and obtain a swCAM solution using an efficient alternating direction method of multipliers. The swCAM is implemented in open-source R scripts. Experimental results on realistic simulation data show that swCAM can accurately estimate subtype-specific expressions in individual samples and successfully extract co-expression networks in particular subtypes that are otherwise unobtainable using bulk expression data. Application of swCAM to bulk-tissue data of 320 samples from bipolar disorder patients and controls identified changes in cell proportions, expression and coexpression modules in patient neurons. Mitochondria related genes showed significant changes suggesting an important role of energy dysregulation in bipolar disorder.

Availability and implementationThe R Scripts of swCAM is freely available at https://github.com/Lululuella/swCAM. A users guide and a vignette are provided.

Contactyuewang@vt.edu

Supplementary informationSupplementary data are available at Bioinformatics online. Liver ductal organoids reconstruct intrahepatic biliary trees in decellularized liver grafts Three-dimensional scaffolds decellularized from native organs are a promising technique to establish engineered liver grafts and overcome the current shortage of donor organs. However, limited sources of bile duct cells and inappropriate cell distribution in bioengineered liver grafts have hindered their practical application. Organoid technology is anticipated to be an excellent tool for the advancement of regenerative medicine. In the present study, we reconstructed intrahepatic bile ducts in a rat decellularized liver graft by recellularization with liver ductal organoids. Using an ex vivo perfusion culture system, we demonstrated the biliary characteristics of repopulated mouse liver organoids, which maintained bile duct markers and reconstructed biliary tree-like networks with luminal structures. We also established a method for the co-recellularization with engineered bile ducts and primary hepatocytes, revealing the appropriate cell distribution to mimic the native liver. We then utilized this model in human organoids to demonstrate the reconstructed bile ducts. Our results show that liver ductal organoids are a potential cell source for bile ducts from bioengineered liver grafts using three-dimensional scaffolds. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible The members of the tribe Brassiceae share a whole genome triplication (WGT), and one proposed model for its formation is a "two-step" pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using POInT (the Polyploidy Orthology Inference Tool). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring prior to the arrival of the third. We also, for the 90,000 individual genes in our study, make parental "subgenome" assignments, inferring, with measured uncertainty, which of the progenitor genomes of the allohexaploidy each gene derives from. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, where subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty. 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. Scaling effects of temperature on parasitism from individuals to host-parasite systems Parasitism is expected to change in a warmer future, but whether warming leads to substantial increases in parasitism remains unclear. Understanding how warming effects on parasitism in individual hosts (e.g., parasite load) translate to effects on population-level parasitism (e.g., prevalence, R0) remains a major knowledge gap. We analyzed the temperature dependence of parasitism at both host and population levels in thirteen empirical vector-borne host-parasite systems and found a strong, significant positive correlation between the thermal optima of individual- and population-level parasitism. We also found a significant, positive correlation in eleven environmentally-transmitted parasite systems, though several of these systems exhibited thermal optima >5{o}C apart between individual and population levels. Similarly, parasitism thermal optima were close to host performance thermal optima in vector-borne systems but not in environmentally-transmitted systems. We then adapted and simulated simple models for both transmission modes and found a similar pattern to the empirical systems: thermal optima in vector-borne systems were more strongly correlated across scales compared to environmentally-transmitted systems. Generally, our results suggest that information on the temperature-dependence, and specifically the thermal optimum, at either the individual- or population-level should provide a useful--though not quantitatively exact--baseline for predicting temperature dependence at the other level, especially in vector-borne parasite systems. Environmentally-transmitted parasitism may operate by a different set of rules, in which temperature-dependence is decoupled in some systems, requiring the need for trait-based studies of temperature-dependence at individual and population levels. Mapping the genetic and environmental aetiology of autistic traits in Sweden and the United Kingdom BackgroundAutistic traits are influenced by both genetic and environmental factors, and are known to vary geographically in prevalence. But to what extent does their aetiology also vary from place to place?

MethodsWe applied a novel spatial approach to data on autistic traits from two large twin studies, the Child and Adolescent Twin Study in Sweden (CATSS; N=16,677, including 8,307 twin pairs) and the Twins Early Development Study in the UK (TEDS; N=11,594, including 5,796 twin pairs), to explore how the influence of nature and nurture on autistic traits varies from place to place.

ResultsWe present maps of gene- and environment-by geography interactions in Sweden and the United Kingdom (UK), showing geographical variation in both genetic and environmental influences across the two countries. In Sweden genetic influences appear higher in the far south and in a band running across the centre of the country. Environmental influences appear greatest in the south and north, with reduced environmental influence across the central band. In the UK genetic influences appear greater in the south, particularly in more central southern areas and the southeast, the Midlands and the north of England. Environmental influences appear greatest in the south and east of the UK, with less influence in the north and the west.

ConclusionsWe hope this systematic approach to identifying aetiological interactions will inspire research to examine a wider range of previously unknown environmental influences on the aetiology of autistic traits. By doing so, we will gain greater understanding of how these environments draw out or mask genetic predisposition and interact with other environmental influences in the development of autistic traits. Understanding the Spin Crossover Dynamical Effects of the Dioxygen Binding and Activation on HOD Enzyme For the cofactor-free 1-H-3-hydroxy-4-oxoquinaldine-2,4-dioxygenase (HOD), the dioxygen (O2) dependent steps are rate-limiting along with a spin state crossover to the singlet spin state. Here, the primary triplet O2 molecule activation on the 2-methyl-3-hydroxy-4(1H)-quinolone (MHQ) is investigated, and the catalytic role of the intersystem crossing effects is highlighted by directly comparing results from the Born-Oppenheimer dynamics and non-adiabatic surface hopping dynamics. This work confirms non-adiabatic dynamical effects are essential to modulate the O2 activation on the substrate MHQ. The time scale of the equilibration and conversion from triplet to singlet state should be in the range of a few hundreds of femtoseconds. We hope this work provides us a fresh look at the underlying physics of dioxygen activation reactions involving more than one spin state. Multi-omics-based label-free metabolic flux inference reveals obesity-associated dysregulatory mechanisms in liver glucose metabolism Glucose homeostasis is maintained by modulation of metabolic flux. Enzymes and metabolites regulate the involved metabolic pathways. Dysregulation of glucose homeostasis is a pathological event in obesity. Analyzing metabolic pathways and the mechanisms contributing to obesity-associated dysregulation in vivo is challenging. Here, we introduce OMELET: Omics-Based Metabolic Flux Estimation without Labeling for Extended Trans-omic Analysis. OMELET uses metabolomic, proteomic, and transcriptomic data to identify changes in metabolic flux, and to quantify contributions of metabolites, enzymes, and transcripts to the changes in metabolic flux. By evaluating the livers of fasting ob/ob mice, we found that increased metabolic flux through gluconeogenesis resulted primarily from increased transcripts, whereas that through the pyruvate cycle resulted from both increased transcripts and changes in substrates of metabolic enzymes. With OMELET, we identified mechanisms underlying the obesity-associated dysregulation of metabolic flux in liver.

HighlightsO_LIWe created OMELET to infer metabolic flux and its regulation from multi-omic data.
C_LIO_LIGluconeogenic and pyruvate cycle fluxes increased in fasting ob/ob mice.
C_LIO_LITranscripts increases mediated the increase in gluconeogenic fluxes in ob/ob mice.
C_LIO_LIIncreases in transcripts and substrates enhanced pyruvate cycle flux in ob/ob mice.
C_LI Differential metabolism of arsenicals regulates Fps1-mediated arsenite transport Arsenic is an environmental toxin that exists mainly as pentavalent arsenate and trivalent arsenite. Both forms activate the yeast SAPK Hog1, but with different consequences. We describe a mechanism by which cells distinguish between these arsenicals through one-step metabolism to differentially regulate the bidirectional glycerol channel Fps1, an adventitious port for arsenite. Cells exposed to arsenate reduce it to thiol-reactive arsenite, which modifies a set of cysteine residues in target proteins; whereas cells exposed to arsenite metabolize it to methylarsenite, which modifies an additional set of cysteine residues. Hog1 becomes arsenylated, which prevents it from closing Fps1. However, this block is overcome in cells exposed to arsenite through methylarsenylation of Acr3, an arsenite efflux pump that we found also regulates Fps1 directly. This adaptation allows cells to restrict arsenite entry through Fps1, but also allows its exit when produced from arsenate exposure. These results have broad implications for understanding how SAPKs activated by diverse stressors can drive stress-specific outputs.

Summary (for JCB table of contents)

Jongmin Lee, David E. Levin

Lee and Levin investigate the question of how different stressors can drive divergent outputs from an active SAPK. The work describes the mechanism by which two forms of arsenic that both activate the yeast SAPK Hog1 differentially regulate the glycerol channel Fps1. Beta-cell specific Insr deletion promotes insulin hypersecretion and improves glucose tolerance prior to global insulin resistance Abstract

Insulin receptor (Insr) protein can be found at higher levels in pancreatic {beta}-cells than in most other tissues, but the consequences of {beta}-cell insulin resistance remain enigmatic. Ins1cre allele was used to delete Insr specifically in {beta}-cells of both female and male mice. Experimental mice were compared to Ins1cre-containing littermate controls at multiple ages and on multiple diets. RNA-seq of purified recombined {beta}-cells revealed transcriptomic consequences of Insr loss, which differed between female and male mice. Action potential and calcium oscillation frequencies were increased in Insr knockout {beta}- cells from female, but not male mice, whereas only male {beta}InsrKO mice had reduced ATP-coupled oxygen consumption rate and reduced expression of genes involved in ATP synthesis. Female {beta}InsrKO and {beta}InsrHET mice exhibited elevated insulin release in perifusion experiments, during hyperglycemic clamps, and following i.p. glucose challenge. Deletion of Insr did not alter {beta}-cell area up to 9 months of age, nor did it impair hyperglycemia-induced proliferation. Based on our data, we adapted a mathematical model to include {beta}-cell insulin resistance, which predicted that {beta}-cell Insr knockout would improve glucose tolerance depending on the degree of whole-body insulin resistance. Indeed, glucose tolerance was significantly improved in female {beta}InsrKO and {beta}InsrHET mice when compared to controls at 9, 21 and 39 weeks, and also in insulin-sensitive 4-week old males. We did not observe improved glucose tolerance in older male mice or in high fat diet-fed mice, corroborating the prediction that global insulin resistance obscures the effects of {beta}-cell specific insulin resistance. The propensity for hyperinsulinemia was associated with mildly reduced fasting glucose and increased body weight. We further validated our main in vivo findings using the Ins1-CreERT transgenic line and found that male mice had improved glucose tolerance 4 weeks after tamoxifen-mediated Insr deletion. Collectively, our data show that loss of {beta}-cell Insr contributes to glucose-induced hyperinsulinemia, thereby improving glucose homeostasis in otherwise insulin sensitive dietary and age contexts. Nuclear and plastid phylogenomic analyses provide insights into the biogeography and reticulate evolution of the Sino-Japanese disjunctive Diabelia (Caprifoliaceae) Understanding biological diversity and the mechanisms of the Sino-Japanese disjunctions are major challenges in eastern Asia biogeography. The Sino-Japanese flora has been broadly studied as an ideal model for plant phylogeography. Diabelia (Caprifoliaceae) is an East Asian genus, with a disjunctive distribution across the Sino-Japanese region. However, relationships within Diabelia remain elusive. In this study, we reconstructed the phylogeny of Diabelia and inferred historical biogeography and evolutionary patterns based on nuclear and plastid sequences from target enrichment and genome skimming approaches, respectively. We found that the main clades within Diabelia were discordant between nuclear and plastid trees. Both nuclear and plastid phylogenetic analyses supported five main clades: D. serrata, D. tetrasepala, D. sanguinea, D. spathulata var. stenophylla and D. spathulata var. spathulata. Species network analyses revealed that Diabelia tetrasepala is likely the result of a hybridization event. Divergence time estimation and ancestral area reconstructions showed that Diabelia originated in Japan during the early Miocene, with subsequent vicariance and dispersal events between Japan and Korea, and between Japan and China. Overall, our results support the division of Diabelia into five main clades and the recognition of five species in the genus. This research provides new insights in the species delimitation and speciation processes of taxonomically complex lineages such as Diabelia. 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. The geometry of cortical representations of touch in rodents Neural responses are often highly heterogeneous non-linear functions of multiple task variables, a signature of a high-dimensional geometry of the neural representations. We studied the representational geometry in the somatosensory cortex of mice trained to report the curvature of objects using their whiskers. High-speed videos of the whisker movements revealed that the task can be solved by linearly integrating multiple whisker contacts over time. However, the neural activity in somatosensory cortex reflects a process of non-linear integration of spatio-temporal features of the sensory inputs. Although the responses at first appear disorganized, we could identify an interesting structure in the representational geometry: different whisker contacts are disentangled variables represented in approximately, but not fully, orthogonal subspaces of the neural activity space. The observed geometry allows linear readouts to perform a broad class of tasks of different complexities without compromising the ability to generalize to novel situations. Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event The anaerobic ammonium oxidation (anammox) bacteria could transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we performed comprehensive phylogenomic analysis and showed a single origin of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which includes implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based approach, we estimated that anammox bacteria originated at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Ga]) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress, bioenergetics and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria. Effect of Adolescent Stress on Adult Morphine-Induced Behavioral Sensitization is Dependent Upon Genetic Background 1Deaths related to opioid use have skyrocketed in the United States, leading to a public health epidemic. Research has shown that both biological (genes) and environmental (stress) precursors are linked to opioid use. In particular, stress during adolescence - a critical period of frontal lobe development - influences the likelihood of abusing drugs. However, little is known about the biological mechanisms through which adolescent stress leads to long-term risk of opioid use, or whether genetic background moderates this response. Male and female C57BL/6J and BALB/cJ mice were exposed to chronic variable social stress (CVSS) or control conditions throughout adolescence and then tested for morphine locomotor sensitization or morphine consumption in adulthood. To examine possible mechanisms that underlie stress-induced changes in morphine behaviors, we assessed physiological changes in response to acute stress exposure and prefrontal cortex miRNA gene expression. Adolescent stress did not influence morphine sensitization or consumption in BALB/cJ animals, and there was limited evidence of stress effects in female C57BL/6J mice. In contrast, male C57BL/6J mice exposed to adolescent CVSS had blunted morphine sensitization compared to control animals; no differences were observed in the acute locomotor response to morphine administration or morphine consumption. Physiologically, C57BL/6J mice exposed to CVSS had an attenuated corticosterone recovery following an acute stressor and downregulation of twelve miRNA in the prefrontal cortex compared to control mice. The specificity of the effects for C57BL/6J versus BALB/cJ mice provides evidence of a gene by environmental interaction influencing opioid behaviors. Long-term differences in stress reactivity or miRNA expression suggests two possible biological mechanisms to evaluate in future research. 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. Plant-environment microscopy tracks interactions of Bacillus subtilis with plant roots across the entire rhizosphere Our understanding of plant-microbe interactions in soil is limited by the difficulty of observing processes at the microscopic scale throughout plants large volume of influence. Here, we present the development of 3D live microscopy for resolving plant-microbe interactions across the environment of an entire seedling growing in a transparent soil in tailor-made mesocosms, maintaining physical conditions for the culture of both plants and microorganisms. A tailor made dual-illumination light-sheet system acquired scattering signals from the plant whilst fluorescence signals were captured from transparent soil particles and labelled microorganisms, allowing the generation of quantitative data on samples approximately 3600 mm3 in size with as good as 5 m resolution at a rate of up to one scan every 30 minutes. The system tracked the movement of Bacillus subtilis populations in the rhizosphere of lettuce plants in real time, revealing previously unseen patterns of activity. Motile bacteria favoured small pore spaces over the surface of soil particles, colonising the root in a pulsatile manner. Migrations appeared to be directed towards the root cap, the point "first contact", before subsequent colonisation of mature epidermis cells. Our findings show that microscopes dedicated to live environmental studies present an invaluable tool to understand plant-microbe interactions. Vitamin C Binds to SARS Coronavirus-2 Main Protease Essential for Viral Replication There is an urgent need for anti-viral agents that treat and/or prevent Covid-19 caused by SARS-Coronavirus (CoV-2) infections. The replication of the SARS CoV-2 is dependent on the activity of two cysteine proteases, a papain-like protease, PL-pro, and the 3C-like protease known as main protease Mpro or 3CLpro. The shortest and the safest path to clinical use is the repurposing of drugs with binding affinity to PLpro or 3CLpro that have an established safety profile in humans. Several studies have reported crystal structures of SARS-CoV-2 main protease in complex with FDA approved drugs such as those used in treatment of hepatitis C. Here, we report the crystal structure of 3CLpro in complex Vitamin C (L-ascorbate) bound to the proteins active site at 2.5 [A]ngstrom resolution. The crystal structure of the Vitamin C 3CLpro complex may aid future studies on the effect of Vitamin C not only on the coronavirus main protease but on related proteases of other infectious viruses. Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. Sex biases in cancer and autoimmune disease incidence are strongly positively correlated with mitochondrial gene expression across human tissues Cancer occurs more frequently in men while autoimmune diseases (AIDs) occur more frequently in women. To explore whether these sex biases have a common basis, we collected 170 AID incidence studies from many countries for tissues that have both a cancer type and an AID that arise from that tissue. Analyzing a total of 182 country-specific, tissue-matched cancer-AID incidence rate sex bias data pairs, we find that, indeed, the sex biases observed in the incidence of AIDs and cancers that occur in the same tissue are positively correlated across human tissues. The common key factor whose levels across human tissues are most strongly associated with these incidence rate sex biases is the sex bias in the expression of the 37 genes encoded in the mitochondrial genome.