The genome sequence of the avian vampire fly (Philornis downsi), an invasive nest parasite of Darwin's finches in Galapagos The invasive avian vampire fly (Philornis downsi) is considered one of the greatest threats to the unique and endemic avifauna of the Galapagos Islands, Ecuador. The fly parasitizes nearly every passerine species, including Darwins finches, in the Galapagos. The fly is thought to have been introduced from mainland Ecuador, although the full pathway of invasion is not yet known. The majority of research to date has focused on the effects of the fly on the fitness of avian hosts and explorations of mitigation methods. A lag in research related to the genetics of this invasion demonstrates, in part, a need to develop full-scale genomic resources with which to address further questions within this system. In this study, an adult P. downsi collected from San Cristobal Island within the Galapagos archipelago was sequenced to generate a high-quality genome assembly. We examined various features of the genome (e.g., coding regions, non-coding transposable elements) and carried out comparative genomics analysis against other dipteran genomes. We identified lists of gene families that are significantly expanding/contracting in P. downsi that are related to insecticide resistance, detoxification, and potential feeding ecology and counter defense against host immune responses. The P. downsi genome assembly provides an important foundational resource for studying the molecular basis of its successful invasion in the Galapagos and the dynamics of its population across multiple islands. The findings of significantly changing gene families associated with insecticide resistance and immune responses highlight the need for further investigations into the role of different gene families in aiding the flys successful invasion. Furthermore, this genomic resource will also better help inform future research studies and mitigation strategies aimed at minimizing the flys impact on the birds of the Galapagos.
Network centrality is the most important variable determining gene sequence evolution rates in Schizosaccharomyces pombe Which variables determine the rate of gene and protein sequence evolution is a central question in evolutionary genomics. In the model organism fission yeast (Schizosaccharomyces pombe), the determinants of the rate of sequence evolution have yet to be determined. Previous studies in other organisms have typically found gene expression levels to be most significant, with numerous other variables identified as having a smaller impact. Here, partial least squares regression (PLS) and partial correlation analysis are used to model sequence evolution rates in the fission yeast genome by a range of variables. Variable importance in projection (VIP) scores as well as partial correlation coefficients are calculated for each variable, and used as estimates of the influence of each independent variable on sequence evolution rate. Unlike many previous studies in other organisms, centrality in the PPI network is shown to be the most important variable, and gene expression found to be less influential. Considerable heterogeneity is found in the influence of different gene ontology terms as well as amino acid composition. However, the majority of variance in constraint in fission yeast remains unexplained by this study, indicating that variables not yet considered and stochastics probably have considerable impact on the rate of molecular evolution.
Microbial loop of a Proterozoic ocean analogue Meromictic Lake Cadagno, an ancient ocean analogue, is known for its permanent stratification and persistent anoxygenic microbial bloom within the chemocline. Although the anaerobic microbial ecology of the lake has been extensively studied for at least 25 years, a comprehensive picture of the microbial food web linking the bacterial layer to phytoplankton and viruses, with explicit measures of primary and secondary production, is still missing. This study sought to better understand the abundances and productivity of microbes in the context of nutrient biogeochemical cycling across the stratified zones of Lake Cadagno. Photosynthetic pigments and chloroplast 16S rRNA gene phylogenies suggested the presence of eukaryotic phytoplankton through the water column. Evidence supported high abundances of Ankyra judayi, a high-alpine adapted chlorophyte, in the oxic mixolimnion where oxygenic-primary production peaked. Through the no- and low-oxygen chemocline and monimolimnion, chlorophytes related to Closteriopsis acicularis, a known genus of meromictic lakes, and Parachlorella kessleri were observed. Chromatium, anoxygenic phototrophic sulfur bacteria, dominated the chemocline along with Lentimicrobium, a genus of known fermenters whose abundance was newly reported in Lake Cadagno. Secondary production peaked in the chemocline suggesting primary producers depend on hetrotrophs for nutrient remineralization. As previously observed, sulfur reducing bacteria (SRBs), especially Desulfocapsa and Desulfobulbus, were present in the chemocline and anoxic monimolimnion. Virus-to-microbe ratios (VMR) peaked in the zone of phytoplankton, yet were at a minimum at the peak of Chromatium. These dynamic trends suggest viruses may play a role in the modulation of oxygenic and anoxygenic photo- and chemosynthesis in Lake Cadagno and other permanently stratified systems.

ImportanceAs a window to the past, the study offers insights into the role of microbial guilds of Proterozoic ocean chemoclines in the production and recycling of organic matter of sulfur- and ammonia-containing ancient oceans. The new observations described here suggest that eukaryotic algae were persistent in the low oxygen upper-chemocline in association with purple and green sulfur bacteria in the lower half of the chemocline. Further, this study provides the first insights into Lake Cadagno viral ecology. High viral abundances suggested viruses may be essential components of the chemocline where their activity may result in the release and recycling of organic matter. The framework developed in this study through the integration of diverse geochemical and biological data types lays the foundation for future studies to quantitatively resolve the processes performed by discrete populations comprising the microbial loop in this early anoxic ocean analogue. Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain. Computational Simulations of the 4-D Micro-Circulatory Network in Zebrafish Tail Amputation and Regeneration Wall shear stress (WSS) in the micro-vasculature contributes to biomechanical cues that regulate mechanotransduction underlying vascular development, regeneration, and homeostasis. We hereby elucidate the interplay between hemodynamic shear forces and luminal remodeling in response to vascular injury and regeneration in the zebrafish model of tail amputation. Using the transgenic Tg(fli1:eGFP; Gata1:ds-red) line, we were able to track the enhanced green-fluorescent protein (eGFP)-labeled endothelial lining of the 3-D microvasculature for post-image segmentation and reconstruction of fluid domain for computational fluid dynamics (CFD) simulation. At 1 day post amputation (dpa), dorsal aorta (DA) and posterior cardinal vein (PCV) were severed, and vasoconstriction developed in the dorsal longitudinal anastomotic vessel (DLAV) with a concomitant increase in WSS in the segmental vessels (SV) proximal to the amputation site and a decrease in WSS in SVs distal to amputation. Simultaneously, we observed angiogenesis commencing at the tips of the amputated DLAV and PCV where WSS was minimal in the absence of blood flow. At 2 dpa, vasodilation occurred in a pair of SVs proximal to amputation, resulting in increased flow rate and WSS, whereas in the SVs distal to amputation, WSS normalized to the baseline. At 3 dpa, the flow rate in the arterial SV proximal to amputation continued to rise and merged with DLAV that formed a new loop with PCV. Thus, our CFD modeling uncovered a well-coordinated micro-vascular adaptation process following tail amputation, accompanied by the rise and fall of WSS and dynamic changes in flow rate during vascular regeneration. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. 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. Development of an exosomal gene signature to detect residual disease in dogs with osteosarcoma using a novel xenograft platform and machine learning Osteosarcoma has a guarded prognosis. A major hurdle in developing more effective osteosarcoma therapies is the lack of disease-specific biomarkers to predict risk, prognosis, or therapeutic response. Exosomes are secreted extracellular microvesicles emerging as powerful diagnostic tools. However, their clinical application is precluded by challenges in identifying disease-associated cargo from the vastly larger background of normal exosome cargo. We developed a method using canine osteosarcoma in mouse xenografts to distinguish tumor-derived from host-response exosomal mRNAs. The model allows for the identification of canine osteosarcoma-specific gene signatures by RNA sequencing and a species-differentiating bioinformatics pipeline. An osteosarcoma-associated signature consisting of five gene transcripts (SKA2, NEU1, PAF1, PSMG2, and NOB1) was validated in dogs with spontaneous osteosarcoma by qRT-PCR, while a machine learning model assigned dogs into healthy or disease groups. Serum/plasma exosomes were isolated from 53 dogs in distinct clinical groups ("healthy", "osteosarcoma", "other bone tumor", or "non-neoplastic disease"). Pre-treatment samples from osteosarcoma cases were used as the training set and a validation set from post-treatment samples was used for testing, classifying as "osteosarcoma-detected" or "osteosarcoma-NOT detected". Dogs in a validation set whose post-treatment samples were classified as "osteosarcoma-NOT detected" had longer remissions, up to 15 months after treatment. In conclusion, we identified a gene signature predictive of molecular remissions with potential applications in the early detection and minimal residual disease settings. These results provide proof-of-concept for our discovery platform and its utilization in future studies to inform cancer risk, diagnosis, prognosis, and therapeutic response.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=147 SRC="FIGDIR/small/429432v2_ufig1.gif" ALT="Figure 1">
View larger version (41K):
org.highwire.dtl.DTLVardef@1b612c4org.highwire.dtl.DTLVardef@8a1847org.highwire.dtl.DTLVardef@3ffd64org.highwire.dtl.DTLVardef@114df8c_HPS_FORMAT_FIGEXP  M_FIG C_FIG Long-range cortical synchronization supports abrupt visual learning Visual plasticity declines sharply after the critical period, yet we easily learn to recognize new faces and places even as adults. Such learning is often characterized by a "moment of insight", an abrupt and dramatic improvement in recognition. The mechanisms that support abrupt learning are unknown, but one hypothesis is that they involve changes in synchronization between brain regions. To test this hypothesis, we used a behavioral task in which non-human primates rapidly learned to recognize novel images and to associate them with specific responses. Simultaneous recordings from inferotemporal and prefrontal cortices revealed a transient synchronization of neural activity between these areas that peaked around the moment of insight. Synchronization was strongest between inferotemporal sites that encoded images and reward-sensitive prefrontal sites. Moreover, its magnitude intensified gradually over image exposures, suggesting that abrupt learning culminates from the search for informative signals within a circuit linking sensory information to task demands. Methylome dynamic upon proteasome inhibition by the Pseudomonas syringae virulence factor Syringolin A DNA methylation is an important epigenetic mark required for proper gene expression and silencing of transposable elements. DNA methylation patterns can be modified by environmental factors such as pathogen infection, where modification of DNA methylation can be associated with plant resistance. To counter the plant defense pathways, pathogens produce effectors molecule, several of which act as proteasome inhibitors. Here we investigated the effect of proteasome inhibition by the bacterial virulence factor Syringolin A on genome-wide DNA methylation. We show that Syringolin A treatment results in an increase of DNA methylation at centromeric and pericentromeric regions of Arabidopsis chromosomes. We identify several CHH DMRs that are enriched in the proximity of transcriptional start sites. Syringolin A treatment does not result in significant changes in small RNA composition. However, significant changes in genome transcriptional activity can be observed, including a strong upregulation of resistance genes that are located on chromosomal arms. We hypothesize that DNA methylation changes could be linked to the upregulation of some atypical members of the de novo DNA methylation pathway: AGO3, AGO9 and DRM1. Our data suggests that modification of genome-wide DNA methylation resulting from an inhibition of the proteasome by bacterial effectors could be part of a epi-genomic arms race against pathogens. Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd. SARS-CoV-2 variants resist antibody neutralization and broaden host ACE2 usage New SARS-CoV-2 variants continue to emerge from the current global pandemic, some of which can replicate faster and with greater transmissibility and pathogenicity. In particular, UK501Y.V1 identified in UK, SA501Y.V2 in South Africa, and BR501Y.V3 in Brazil are raising serious concerns as they spread quickly and contain spike protein mutations that may facilitate escape from current antibody therapies and vaccine protection. Here, we constructed a panel of 28 SARS-CoV-2 pseudoviruses bearing single or combined mutations found in the spike protein of these three variants, as well as additional nine mutations that within or close by the major antigenic sites in the spike protein identified in the GISAID database. These pseudoviruses were tested against a panel of monoclonal antibodies (mAbs), including some approved for emergency use to treat SARS-CoV-2 infection, and convalescent patient plasma collected early in the pandemic. SA501Y.V2 pseudovirus was the most resistant, in magnitude and breadth, against mAbs and convalescent plasma, followed by BR501Y.V3, and then UK501Y.V1. This resistance hierarchy corresponds with Y144del and 242-244del mutations in the N-terminal domain as well as K417N/T, E484K and N501Y mutations in the receptor binding domain (RBD). Crystal structural analysis of RBD carrying triple K417N-E484K-N501Y mutations found in SA501Y.V2 bound with mAb P2C-1F11 revealed a molecular basis for antibody neutralization and escape. SA501Y.V2 and BR501Y.V3 also acquired substantial ability to use mouse and mink ACE2 for entry. Taken together, our results clearly demonstrate major antigenic shifts and potentially broadening the host range of SA501Y.V2 and BR501Y.V3, which pose serious challenges to our current antibody therapies and vaccine protection. Natural variations in reproductive timing of brood size and X-chromosome nondisjunction in Caenorhabditis elegans Reproductive timing underlies self-sperm production and X-chromosome nondisjunction rate in Caenorhabditis elegans. These traits can be measured as brood size and male-production rate, and brood size decreases and male-production rate increases as worms age. These phenotypic changes occur simultaneously; however, whether these changes share common underlying genetic architectures still remains unclear. To enhance our understanding of reproductive timing in C. elegans, we measured and compared these two reproductive timing-related traits of 96 wild strains during early, late and total reproductive periods. We found that the two traits exhibited natural phenotypic variation with few outliers and a similar reproductive timing pattern as previous reports. The brood size and the male-production rate were not genetically correlated, implying that they might have different genetic architectures and that meiosis quality can be maintained despite increased progeny-production in the late reproductive period caused by more self-sperm. We also identified loci and candidate genetic variants significantly associated with male-production rate in the late and total reproductive periods. Our results provide an insight into life history traits in wild C. elegans strains. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. Standardized genome-wide function prediction enables comparative functional genomics: a new application area for Gene Ontologies in plants BackgroundGenome-wide gene function annotations are useful for hypothesis generation and for prioritizing candidate genes potentially responsible for phenotypes of interest. We functionally annotated the genes of 18 crop plant genomes across 14 species using the GOMAP pipeline.

ResultsBy comparison to existing GO annotation datasets, GOMAP-generated datasets cover more genes, contain more GO terms, and produce datasets similar in quality (based on precision and recall metrics using existing gold standards as the basis for comparison). From there, we sought to determine whether the datasets across multiple species could be used together to carry out comparative functional genomics analyses in plants. To test the idea and as a proof of concept, we created dendrograms of functional relatedness based on terms assigned for all 18 genomes. These dendrograms were compared to well-established species-level evolutionary phylogenies to determine whether trees derived were in agreement with known evolutionary relationships, which they largely are. Where discrepancies were observed, we determined branch support based on jack-knifing then removed individual annotation sets by genome to identify the annotation sets causing unexpected relationships.

ConclusionsGOMAP-derived functional annotations used together across multiple species generally retain sufficient biological signal to recover known phylogenetic relationships based on genome-wide functional similarities, indicating that comparative functional genomics across species based on GO data hold promise for generating novel hypotheses about comparative gene function and traits. A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. A vesicle microrheometer for high-throughput viscosity measurements of lipid and polymer membranes Viscosity is a key property of cell membranes that controls mobility of embedded proteins and membrane remodeling. Measuring it is challenging because existing approaches involve complex experimental designs and/or models, and the applicability of some is limited to specific systems and membrane compositions. As a result there is scarcity of data and the reported values for membrane viscosity vary by orders of magnitude for the same system. Here, we show how viscosity of bilayer membranes can be obtained from the transient deformation of giant unilamellar vesicles. The approach enables a non-invasive, probe-independent and high-throughput measurement of the viscosity of bilayers made of lipids or polymers with a wide range of compositions and phase state. Pure lipid and single-phase mixed bilayers are found to behave as Newtonian fluids with strain-rate independent viscosity, while phase-separated and diblock-copolymers systems exhibit shear-thinning in the explored range of strain rates 1-2000 s-1. The results also reveal that electrically polarized bilayers can be significantly more viscous than charge-neutral bilayers. These findings suggest that biomembrane viscosity is a dynamic property that can be actively modulated not only by composition but also by membrane polarization, e.g., as in action potentials. A thermodynamic model for water activity and redox potential in evolution and development Reactions involving water and oxygen are basic features of geological and biological processes. To understand how life interacts with its environment requires monitoring interactions with H2O and O2 at timescales relevant to not only organismal growth but also billions of years of geobiological evolution. Here, chemical transformations intrinsic to evolution and development were characterized by analyzing data from recent phylostratigraphic and proteomic studies. This two-stage analysis involves obtaining chemical metrics (carbon oxidation state and stoichiometric hydration state) from the elemental compositions of proteins followed by modeling the relative stabilities of target proteins against a proteomic background to infer thermodynamic parameters (oxygen fugacity, water activity, and virtual redox potential (Eh)). The main results of this study are a rise in carbon oxidation state of proteins spanning the time of the Great Oxidation Event, a rise in virtual redox potential that coincides with the likely emergence of aerobic metabolism, a rise in carbon oxidation state of proteins inferred from the transcriptome in late stages of Bacillus subtilis biofilm growth, and a drop in stoichiometric hydration state of the fruit fly developmental proteome at the same time as a drop in organismal water content. Stoichiometric hydration state also decreases for proteins with more recent gene ages and through stages of biofilm development, leading to predictions of higher hydration potentials at earlier time points. By building chemical representations of protein evolution and developmental proteomes, exciting new types of geobiochemical proxies can be developed. A constraints-based theory of the primary cause of senescence: imbalance of epigenetic and non-epigenetic information in histone crosstalk Cellular aging has been progressively elucidated by science. However, the fundamental cause of senescence--i.e., why organisms age at the multicellular-individual level--remains unclear. A recent theory of individuated multicellularity describes the emergence and growth of crucial information content for cell differentiation. This information is mostly conveyed in the non-epigenetic (i.e., transcription uncorrelated) histone crosstalk near transcription start sites. According to this theory, the non-epigenetic content emerges and grows at the expense of the information capacity for epigenetic content. If this "reassignment" of information capacity continues after adulthood, it may explain the senescence phenomenon. Here, I present a novel, falsifiable theory describing an uninterrupted growth of capacity for non-epigenetic information at the expense of that for epigenetic information not only during ontogeny but also throughout adulthood. As a byproduct, this continuous "reassignment" of capacity effectively creates an information imbalance in histone crosstalk, which in turn overregulates transcriptional levels. This overregulation is to be understood as transcriptional levels becoming more and more accurate but also less and less precise with respect to the needs of the multicellular individual--up to the point of dysfunctionality. This epigenetic/non-epigenetic information imbalance is proposed to be the primary reason why individuated multicellular organisms senesce. 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. A small molecule that promotes cellular senescence prevents fibrogenesis and tumorigenesis in vitro BackgroundUncontrolled proliferative diseases such as fibrosis or cancer can be fatal. We previously found that a compound containing a chromone scaffold, ONG41008, had potent anti-fibrotic effects in diseased human lung myofibroblasts but not normal human lung fibroblasts.

MethodsWe investigated the effects of ONG41008 on tumor cells, and compared these effects with those in pathogenic myofibrotic cells and normal fibroblasts cells.

FindingsStimulation of A549 lung carcinoma epithelial cells with ONG41008 resulted in cellular senescence, indicating that dysregulated cell proliferation is common to fibrotic cells and tumor cells. Replicative senescence of A549 cells resulted in multinucleation, which was followed by oncogene-induced senescence. There was significant upregulation of expression and nuclear translocation of p-TP53 and p16 in ONG41008-treated A549 cells, and all cells died after 72 hr. Similar effects occurred after ONG41008 treatment in several human aggressive cancer cell lines such as PANC1, MCF7, PC3, or primary non-small cell lung carcinoma cells. Unlike cisplatin, ONG41008 was not toxic to normal human lung fibroblasts or primary prostate epithelial cells, suggesting ONG41008 can distinguish the intracellular microenvironment between normal cells and aged or diseased cells. This effect might occur as a result of the increased NAD/NADH ratio or increased lactate dehydrogenase levels in aged or diseased cells.

InterpretationTo our best knowledge, this is the first study to show that a small molecule can arrest uncontrolled proliferation during fibrogenesis or tumorigenesis in vitro. ONG41008 could be a potential drug for a broad range of fibrotic or tumorigenic diseases.

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSThe notion that aging is a disease and that diseases occur as a consequence of aging was first put forward by David Sinclair and colleagues (Aging Cell;14(4):497-510). There are biological systems that provide evidence for this notion; for example, somatic cells can revert to embryonic cells, producing younger somatic cells. This phenomenon underlies induced pluripotent stem cells. Another example is that some types of jellyfish can live forever. These examples suggest that a counter-aging program exists in animals. Human diseases are the manifestations of cell aging generated by the accumulation of somatic mutations. Aged and pathogenic cells are senescent, so a drug that specifically targeted senescent cells might initiate a cellular program that could ameliorate age-associated disease. Indeed, the kinase inhibitor dasatinib induces cellular senescence (Clin Ther 2007 29:2289-2308). In 2017, two drugs that target senescent cells were identified: quercetin and fisetin. These drugs selectively kill senescent cells, and are referred to as senolytic drugs or senotherapeutics (Aging 2017 8;9(3):955-963). Although it is established that senescent cells accumulate in cancer and idiopathic pulmonary fibrosis (Nat Commun 2017 23;8 :14532), the effect of senolytic drugs in these diseases is largely unknown.

Added value of this studyThis study characterized a novel drug, termed ONG41008, which was found to have both senogenic and senolytic effects in cell-based assays. ONG41008 induced senescence in myofibroblasts and several cancer cell lines representative of aggressive human cancers, which was followed by cell death. Importantly, ONG41008 exhibited essentially no toxicity on normal human lung fibroblasts or primary prostate epithelial cells.

Implication of all the evidenceBased on our results, we believe that ONG41008 is a potent inducer of cellular senescence (replicative senescence and oncogene-induced senescence) and causes arrest of uncontrolled, pathogenic proliferation of myofibroblasts or cancer cells. Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems. Polycomb-mediated repression of paternal chromosomes maintains haploid dosage in diploid embryos of Marchantia Complex mechanisms regulate gene dosage throughout eukaryotic life cycles. Mechanisms controlling gene dosage have been extensively studied in animals, however it is unknown how generalizable these mechanisms are to diverse eukaryotes. Here, we use the haploid plant Marchantia polymorpha to assess gene dosage control in its short-lived diploid embryo. We show that throughout embryogenesis, paternal chromosomes are repressed resulting in functional haploidy. The paternal genome is targeted for genomic imprinting by the Polycomb mark H3K27me3 starting at fertilization, rendering the maternal genome in control of embryogenesis. Maintaining haploid gene dosage by this new form of imprinting is essential for embryonic development. Our findings illustrate how haploid-dominant species can regulate gene dosage through paternal chromosome inactivation and initiates the exploration of the link between life cycle history and gene dosage in a broader range of organisms. Environment-dependence of the expression of mutational load and species' range limits Theoretical and empirical research on the causes of species range limits suggests the contribution of several intrinsic and extrinsic factors, with potentially complex interactions among them. An intrinsic factor proposed by recent theory is mutational load increasing towards range edges because of genetic drift. Furthermore, environmental quality may decline toward range edges and enhance the expression of load. Here we tested whether the expression of mutational load associated with range limits in the North American plant Arabidopsis lyrata is enhanced under stressful environmental conditions by comparing the performance of within- versus between-population crosses at common garden sites across the species distribution and beyond. Heterosis, reflecting the expression of load, increased with heightened estimates of genomic load and with environmental stress caused by warming, but the interaction was not significant. We conclude that range-edge populations suffer from a twofold genetic Allee effect caused by increased mutational load and stress-dependent load linked to general heterozygote deficiency, but no synergistic effect between them. Community composition exceeds area as a predictor of long-term conservation value Conserving biodiversity often requires deciding which sites to prioritise for protection. Predicting the impact of habitat loss is a major challenge, however, since impacts can be distant from the perturbation in both space and time. Here we study the long-term impacts of habitat loss in a mechanistic metacommunity model in terms of both immediate extinctions and secondary species losses. We find that biomass-at-site, closely related to site area, is a poor predictor of long-term regional species losses following site removal. Knowledge of the compositional distinctness (average between-site Bray-Curtis dissimilarity) of the removed community can markedly improve the prediction of impacts at the regional scale, even when biotic responses play out at substantial spatial or temporal distance from the removed site. Fitting our model to empirical species-by-site tables describing Andean diatoms and Brazilian lichen-fungi, we show that compositional distinctness surpasses area as a predictor of long-term species losses in the empirically relevant parameter range. Our results robustly demonstrate that site area alone is not sufficient to gauge conservation priorities; analysis of compositional distinctness permits improved prioritisation at low cost. 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.