Experimental sexual selection affects the evolution of physiological and life history traits Sexual selection and sexual conflict are expected to affect all aspects of the phenotype, not only traits that are directly involved in reproduction. Here, we show coordinated evolution of multiple physiological and life history traits in response to long-term experimental manipulation of the mating system in populations of Drosophila pseudoobscura. Development time was extended under polyandry relative to monogamy in both sexes, potentially due to higher investment in traits linked to sexual selection and sexual conflict. Individuals (especially males) evolving under polyandry had higher metabolic rates and locomotor activity than those evolving under monogamy. Polyandry individuals also invested more in metabolites associated with increased endurance capacity and efficient energy metabolism and regulation, namely lipid and glycogen. Finally, polyandry males were less desiccation- and starvation-resistant than monogamy males, suggesting trade-offs between resistance and sexually selected traits. Our results provide experimental evidence that mating systems can impose selection that influences the evolution of non-sexual phenotypes such as development, activity, metabolism, and nutrient homeostasis.
Ecological correlates of gene family size in a pine-feeding sawfly genome and across Hymenoptera A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a non-eusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions and characterized their genomic distributions and evolutionary history. Our results suggest that expansions of bitter gustatory receptor (GR), clan 3 cytochrome P450 (CYP3), and antimicrobial peptide (AMP) subfamilies may have contributed to pine adaptation. By contrast, there was no evidence of recent gene family contraction via pseudogenization. Next, we compared the number of genes in these same families across insect taxa that vary in diet, dietary specialization, and social behavior. In Hymenoptera, herbivory was associated with small GR and olfactory receptor (OR) families, eusociality was associated with large OR and small AMP families, and--unlike investigations in more closely related taxa--ecological specialization was not related to gene family size. Overall, our results suggest that gene families that mediate ecological interactions may expand and contract predictably in response to particular selection pressures, however, the ecological drivers and temporal pace of gene gain and loss likely varies considerably across gene families.
Chromatin-contact atlas reveals disorder-mediated protein interactions and moonlighting chromatin-associated RBPs RNA-binding proteins (RBPs) play diverse roles in regulating co-transcriptional RNA-processing and chromatin functions, but our knowledge of the repertoire of chromatin-associated RBPs (caRBPs) and their interactions with chromatin remains limited. Here, we developed SPACE (Silica Particle Assisted Chromatin Enrichment) to isolate global and regional chromatin components with high specificity and sensitivity, and SPACEmap to identify the chromatin-contact regions in proteins. Applied to mouse embryonic stem cells, SPACE identified 1,459 chromatin-associated proteins, [~]48% of which are annotated as RBPs, indicating their dual roles in chromatin and RNA-binding. Additionally, SPACEmap stringently verified chromatin-binding of 404 RBPs and identified their chromatin-contact regions. Notably, SPACEmap showed that about half of the caRBPs bind to chromatin by intrinsically disordered regions (IDRs). Studying SPACE and total proteome dynamics from mES cells grown in 2iL and serum medium indicates significant correlation (R = 0.62). One of the most dynamic caRBPs is Dazl, which we find co-localized with PRC2 at transcription start sites of genes that are distinct from Dazl mRNA binding. Dazl and other PRC2-colocalised caRBPs are rich in intrinsically disordered regions (IDRs), which could contribute to the formation and regulation of phase-separated PRC condensates. Together, our approach provides an unprecedented insight into IDR-mediated interactions and caRBPs with moonlighting functions in native chromatin. A Histoplasma capsulatum lipid metabolic map identifies antifungal targets Lipids play a fundamental role in fungal cell biology, being essential cell membrane components and major targets of antifungal drugs. A deeper knowledge of lipid metabolism is key for developing new drugs and a better understanding of fungal pathogenesis. Here we built a comprehensive map of the Histoplasma capsulatum lipid metabolic pathway by incorporating proteomic and lipidomic analyses. We performed genetic complementation and overexpression of H. capsulatum genes in Saccharomyces cerevisiae to validate reactions identified in the map and to determine enzymes responsible for catalyzing orphan reactions. The map led to the identification of both the fatty acid desaturation and the sphingolipid biosynthesis pathways as targets for drug development. We found that the sphingolipid biosynthesis inhibitor myriocin, the fatty acid desaturase inhibitor thiocarlide and the fatty acid analog 10-thiastearic acid inhibit H. capsulatum growth in nanomolar to low micromolar concentrations. These compounds also reduced the intracellular infection in an alveolar macrophage cell line. Overall, this lipid metabolic map revealed pathways that can be targeted for drug development. Abscisic Acid Drives Stomatal Closure through Increases in Hydrogen Peroxide in Distinct Subcellular Compartments Including Mitochondria Abscisic acid (ABA) closes stomata to minimize transpiration and gas exchange in response to drought. ABA signaling increases accumulation of reactive oxygen species (ROS) in guard cells to drive stomatal closure. We used chemical dyes and genetically-encoded biosensors to identify the subcellular localizations where ABA increases hydrogen peroxide (H2O2) in guard cells. The general redox sensor, dichlorofluorescein (DCF), revealed ABA-driven ROS increases in distinct subcellular locations including chloroplasts, cytosol, and nuclei. ABA also increased fluorescence of both DCF and an H2O2-selective probe, Peroxy Orange 1, in mitochondria, which colocalize with Mitotracker Red. Whole cell and mitochondrial localized ABA-increased ROS signals were lost in an ABA-insensitive quintuple ABA receptor mutant and accentuated in an ABA-hypersensitive mutant. The ABA-induced ROS accumulation in mitochondria was lost in guard cells of mutants with defects in genes encoding H2O2 elevating respiratory burst oxidase homolog (RBOH) enzymes and reduced by treatment with the RBOH inhibitor VAS2870. ABA treatment of guard cells transformed with the genetically-encoded H2O2 reporter targeted to the cytoplasm, roGFP2-Orp1, and the mitochondrial targeted, mt-roGFP2-Orp1, revealed significant ABA-induced H2O2 increases. Furthermore, H2O2 treatment is sufficient to close stomata. These results reveal that ABA elevates H2O2 accumulation in guard cell cytosol and mitochondria to drive stomatal closure.

One sentence summaryGenetically encoded biosensors and chemical probes revealed ABA-dependent increases in hydrogen peroxide, a reactive oxygen species with signaling activity, in guard cell cytoplasm and mitochondria. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. The Backbone Network of Dynamic Functional Connectivity 1.Temporal networks have become increasingly pervasive in many real-world applications, including the functional connectivity analysis of spatially separated regions of the brain. A major challenge in analysis of such networks is the identification of noise confounds, which introduce temporal ties that are non-essential, or links that are formed by chance due to local properties of the nodes. Several approaches have been suggested in the past for static networks or temporal networks with binary weights for extracting significant ties whose likelihood cannot be reduced to the local properties of the nodes. In this work, we propose a data-driven procedure to reveal the irreducible ties in dynamic functional connectivity of resting state fRMI data with continuous weights. This framework includes a null model that estimates the latent characteristics of the distributions of temporal links through optimization, followed by a statistical test to filter the links whose formation can be reduced to the activities and local properties of their interacting nodes. We demonstrate the benefits of this approach by applying it to a resting state fMRI dataset, and provide further discussion on various aspects and advantages of it. Interpreting ciliopathy-associated missense variants of uncertain significance (VUS) in an animal model PurposeBetter methods are required to interpret the pathogenicity of disease-associated variants of uncertain significance (VUS), which cannot be actioned clinically. In this study, we explore the use of a tractable animal model (Caenorhabditis elegans) for in vivo interpretation of missense VUS alleles of TMEM67, a cilia gene associated with ciliopathies.

MethodsCRISPR/Cas9 gene editing was used to generate homozygous knock-in C. elegans worm strains carrying TMEM67 patient variants engineered into the orthologous gene (mks-3). Quantitative phenotypic assays of sensory cilia structure and function measured if the variants affect mks-3 gene function. Results from worms were validated by a genetic complementation assay in a human TMEM67 knock-out hTERT-RPE1 cell line that tests a TMEM67 signaling function.

ResultsAssays in C. elegans accurately distinguished between known benign (Asp359Glu, Thr360Ala) and known pathogenic (Glu361Ter, Gln376Pro) variants. Analysis of eight missense VUS generated evidence that three are benign (Cys173Arg, Thr176Ile, Gly979Arg) and five are pathogenic (Cys170Tyr, His782Arg, Gly786Glu, His790Arg, Ser961Tyr).

ConclusionEfficient genome editing and quantitative functional assays in C. elegans make it a tractable in vivo animal model that allows rapid, cost-effective interpretation of ciliopathy-associated missense VUS alleles. The perceptual relevance of the speech envelope-following response in normal hearing and age-related hearing impaired listeners The aim of this work was to investigate the perceptual relevance of the frequency following response to the syllable /da/ for speech intelligibility in noise based on age and hearing deficits. Recordings of the auditory evoked potential from young normal hearing (NH) and older individuals with both normal hearing and high-frequency (HF) hearing loss were analyzed. EFR metrics obtained in quiet and noise condition were calculated and correlated with speech reception. The envelope following responses were analyzed in terms of amplitude, latency and noise robustness. The response was first simulated to form predictions on the effect of cochlear synaptopathy and outer hair cell loss on the EFR. The experimental findings were in line with the computational predictions in the found observation that the EFR was reduced as a consequence of ageing and HF hearing loss. Both the audiogram and the speech EFR magnitude fell short in the individual prediction of SRT in stationary noise, but they accounted well for group performance. We also obtained within-group EFR latency with a cross covariance matrix. Validation of the method confirmed that speech EFR latency was predictive of click ABR Wave V peak latency. Moreover, statistical analysis not only showed that the robustness of the EFR obtained in the noise condition was dependent on the degree of high-frequency hearing loss in the older NH adults, but also dependent on the EFR magnitude in the NH younger adults. These findings provide evidence towards the important role of the EFR in speech-in-noise perception. How many markers are needed to robustly determine a cell's type? Our understanding of cell types has advanced considerably with the publication of single cell atlases. Marker genes play an essential role for experimental validation and computational analyses such as physiological characterization through pathway enrichment, annotation, and deconvolution. However, a framework for quantifying marker replicability and picking replicable markers is currently lacking. Here, using high quality data from the Brain Initiative Cell Census Network (BICCN), we systematically investigate marker replicability for 85 neuronal cell types. We show that, due to dataset-specific noise, we need to combine 5 datasets to obtain robust differentially expressed (DE) genes, particularly for rare populations and lowly expressed genes. We estimate that 10 to 200 meta-analytic markers provide optimal performance in downstream computational tasks. Replicable marker lists condense single cell atlases into interpretable and generalizable information about cell types, opening avenues for downstream applications, including cell type annotation, selection of gene panels and bulk data deconvolution. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. High resolution, serial imaging of early mouse and human liver bud morphogenesis in three dimensions BackgroundLiver organogenesis has thus far served as a paradigm for solid organ formation. The developing liver bud is a well established model of organogenesis, and murine genetic studies demonstrate key molecules involved in key morphogenetic changes. However, the analysis of the liver bud is typically limited to 2D tissue sections, which precludes extensive visualization, quantitation, and analysis. Further, the lack of human liver bud data has further hindered our understanding of human liver organogenesis. Therefore, new analytical and visualization approaches are needed to elicit further morphogenetic details of liver organogenesis and to elucidate differences between mouse and human liver bud growth.

ResultsTo address this need, we focused on high resolution imaging, visualization, and analysis of early liver growth by using available online databases for both mouse (EMAP, Toronto Phenogenomics center) and human (3D Atlas of Human Embryology), noninvasive multimodality imaging studies of the murine embryo, and mouse/human liver weight data. First, we performed three-dimensional (3D reconstructions) of stacked, digital tissue sections that had been initially segmented for the liver epithelium and the septum transversum mesenchyme (STM). 3D reconstruction of both mouse and human data sets enabled visualization and analysis of the dynamics of liver bud morphogenesis, including hepatic cord formation and remodeling, mechanisms of growth, and liver-epithelial STM interactions. These studies demonstrated potentially under-appreciated mechanisms of growth, including rapid exponential growth that is matched at the earliest stages by STM growth, and unique differences between mouse and human liver bud growth. To gain further insight into the exponential liver bud growth that was observed, we plotted volumetric data from 3D reconstruction together with fetal liver growth data from multimodality (optical projection tomography, magnetic resonance imaging, micro-CT) and liver weight data to compose complete growth curves during mouse (E8.5-E18) and human (day 25-300) liver development. For further analysis, we performed curve fitting and parameter estimation, using Gompertzian models, which enables the comparison between mouse and human liver bud growth, as well as comparisons to processes like liver regeneration. To demonstrate the importance of mesenchyme in rapid liver bud growth and morphogenesis in the human liver bud, we performed functional analysis in which human pluripotent stem cell (hPSC)-derived hepatic organoids were used to model collective migration that occurs in vivo, demonstrating that migration is strongly dependent upon mesenchyme.

DiscussionOur data demonstrates improved visualization with 3D images, under-appreciated and potentially new mechanisms of growth, complete liver growth curves with quantitative analysis through embryonic and fetal stages, and a new functional human stem cell-derived liver organoid assay demonstrating mesenchyme-driven collective migration. These data enhance our understanding of liver organogenesis. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. Tetrazoles as PPARγ ligands: A Structural and Computational Investigation. Diabetes is an important chronic disease affecting about 10% of the adult population in the US and over 420 million people worldwide, resulting in 1.6 million deaths every year, according to the World Health Organization. The most common type of the disease, type 2 diabetes, can be pharmacologically managed using oral hypoglycemic agents or thiazolidinediones (TZDs), such as pioglitazone, which act by activating the Peroxisome Proliferated-Activated Receptor {gamma}. Despite their beneficial effects in diabetes treatment, TZDs like rosiglitazone and troglitazone were withdrawn due to safety reasons, creating a void in the pharmacological options for the treatment of this important disease. Here, we explored a structure-based approach in the screening for new chemical probes for a deeper investigation of the effects of PPAR{gamma} activation. A class of tetrazole compounds was identified and the compounds named T1, T2 and T3 were purchased and evaluated for their ability to interact with the PPAR{gamma} ligand binding domain (LBD). The compounds were binders with micromolar range affinity, as determined by their IC50 values. A Monte Carlo simulation of the compound T2 revealed that the tetrazole ring makes favorable interaction with the polar arm of the receptor binding pocket. Finally, the crystal structure of the PPAR{gamma}-LBD-T2 complex was solved at 2.3 [A], confirming the binding mode for this compound. The structure also revealed that, when the helix H12 is mispositioned, an alternative binding conformation is observed for the ligand suggesting an H12-dependent binding conformation for the tetrazole compound. Beyond Spike: Identification of nine highly prevalent SARS-CoV-2-specific CD8 T-cell epitopes in a large Norwegian cohort T-cell epitopes with broad population coverage may form the basis for a new generation of SARS-CoV-2 vaccines. However, published studies on immunoprevalence are limited by small test cohorts, low frequencies of antigen-specific cells and lack of data correlating eluted HLA ligands with T-cell responsiveness. As the protective role of pre-existing cross-reactivity to homologous peptides is unclear, we aimed to identify SARS-CoV-2-specific minimal epitopes recognized by CD8 T-cells among 48 peptides eluted from prevalent HLA alleles, and an additional 84 predicted binders, in a large cohort of convalescents (n=83) and pre-pandemic control samples (n=19). We identified nine conserved SARS-CoV-2-specific epitopes restricted by four of the most prevalent HLA class I alleles in the Norwegian study cohort, to which responding CD8 T cells were detected in 70-100% of convalescents expressing the relevant HLA allele. Only two of these were derived from the Spike protein, included in current vaccines. We found a strong correlation between immunoprevalence and immunodominance. Thus, the CD8 T-cell response to SARS-CoV-2 is more focused than previously believed. Using a new algorithm, we predict that a vaccine including these epitopes could induce a T-cell response in 83% of Caucasians. C. elegans PEZO-1 is a Mechanosensitive Ion Channel Involved in Food Sensation PIEZO channels are force sensors essential for physiological processes, including baroreception and proprioception. The Caenorhabditis elegans genome encodes an ortholog gene of the Piezo family, pezo-1, which is expressed in several tissues, including the pharynx. This myogenic pump is an essential component of the C. elegans alimentary canal, whose contraction and relaxation are modulated by mechanical stimulation elicited by food content. Whether pezo-1 encodes a mechanosensitive ion channel and contributes to pharyngeal function remain unknown. Here, we leverage genome editing, genetics, microfluidics, and electropharyngeogram recording to establish that pezo-1 is expressed in the pharynx, including in a proprioceptive-like neuron, and regulates pharyngeal function. Knockout (KO) and gain-of-function (GOF) mutants reveal that pezo-1 is involved in fine-tuning pharyngeal pumping frequency, sensing osmolarity, and food mechanical properties. Using pressure-clamp experiments in primary C. elegans embryo cultures, we determine that pezo-1 KO cells do not display mechanosensitive currents, whereas cells expressing wild-type or GOF PEZO-1 exhibit mechanosensitivity. Moreover, infecting the Spodoptera frugiperda cell line with a baculovirus containing the G-isoform of pezo-1 (among the longest isoforms) demonstrates that pezo-1 encodes a mechanosensitive channel. Our findings reveal that pezo-1 is a mechanosensitive ion channel that regulates food sensation in worms. Directed evolution of angiotensin-converting enzyme 2 (ACE2) peptidase activity profiles for therapeutic applications Angiotensin Converting-Enzyme 2 (ACE2) is currently being investigated for its ability to beneficially modulate the Angiotensin receptor (ATR) therapeutic axis to treat multiple human diseases, but its broad substrate scope and diverse physiological roles limit its potential as a therapeutic agent. In this work we engineer ACE2 variants with enhanced hydrolytic activity and specificity toward Angiotensin-II (Ang-II) that may lead to new therapeutic candidates with improved efficacy and reduced off-target side effects. We established a yeast display-based liquid chromatography screen that enabled use of directed evolution to discover ACE2 variants with improved Ang-II activity and specificity relative to the off-target peptide substrates Apelin-13 and Angiotensin-I (Ang-I). We screened ACE2 active site libraries to reveal three substitution-tolerant positions (M360, T371 and Y510) that can be mutated to enhance ACE2s activity profile and followed up on these hits with focused double mutant libraries to further improve the enzyme. Relative to wildtype ACE2, our top variant (T371L/Y510Ile) displayed a sevenfold increase in Ang-II turnover number (kcat), a sixfold diminished catalytic efficiency (kcat/Km) on Apelin-13, and an overall decreased activity on other ACE2 substrates that were not directly assayed in directed evolution screen. At physiologically relevant substrate concentrations, T371L/Y510Ile hydrolyzes more Ang-II than wildtype ACE2 with concomitant Ang-II:Apelin-13 specificity improvements reaching 30-fold. Our efforts have delivered ATR axis-acting therapeutic candidates with relevance to both established and unexplored ACE2 therapeutic applications and provide a foundation for further ACE2 engineering efforts. 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. Social integration predicts survival in female white-faced capuchin monkeys Across multiple species of social mammals, a growing number of studies have found that that individual sociality is associated with survival. In long-lived species, like primates, lifespan is one of the main components of fitness. We used 18 years of data from the Lomas Barbudal Monkey project to quantify social integration in 11 capuchin (Cebus capucinus) groups and tested whether female survivorship was associated with females tendencies to interact with three types of partners: (1) all group members, (2) adult females, (3) adult males. We found strong evidence that females who engaged more with other females in affiliative interactions and foraged in close proximity experienced increased survivorship. We found some weak evidence that females might also benefit from engaging in more support in agonistic contexts with other females. These benefits were evident in models that account for the females rank and group size. Female interactions with all group members also increased survival, but the estimates of the effects were more uncertain. In interactions with adult males, only females who provided more grooming to males survived longer. The results presented here suggest that social integration may result in survival-related benefits. Females might enjoy these benefits through exchanging grooming for other currencies, such as coalitionary support or tolerance. α-Phenylalanyl tRNA synthetase competes with Notch signaling through its N-terminal domain The alpha subunit of the cytoplasmic Phenylalanyl tRNA synthetase (-PheRS, FARSA in humans) displays cell growth and proliferation activities and its elevated levels can induce cell fate changes and tumor-like phenotypes that do neither dependent on the canonical function of charging tRNAPhe with phenylalanine nor on stimulating general translation. In intestinal stem cells of Drosophila midguts, -PheRS levels are naturally slightly elevated and human FARSA mRNA levels are elevated in multiple cancers. In the Drosophila midgut model, elevated -PheRS levels caused the accumulation of many additional proliferating cells resembling intestinal stem cells (ISCs) and enteroblasts (EBs). This phenotype resembles partially the tumor-like phenotype described as Notch RNAi phenotype for the same cells. Genetic interactions between -PheRS and Notch suggest that their activities neutralize each other and that elevated -PheRS levels attenuate Notch signaling whether Notch induces differentiation into enterocytes, neuroblast stem cell proliferation, or transcription of a Notch reporter. These non-canonical functions all map to the N-terminal part of -PheRS which accumulates naturally in the intestine. This truncated version of -PheRS (-S) also localizes to nuclei and displays weak sequence similarity to the Notch intracellular domain (NICD), suggesting that -S might compete with the NICD for binding to a common target. Supporting this hypothesis, the tryptophan (W) residue reported to be key for the interaction between the NICD and the Su(H) BTD domain is not only conserved in -PheRS and -S, but also essential for attenuating Notch signaling.

Author SummaryAminoacyl tRNA synthetases charge tRNAs with their cognate amino acid to ensure proper decoding of the genetic code during translation. Independent of its aminoacylation function, the alpha subunit of Drosophila cytoplasmic Phenylalanyl tRNA synthetase (-PheRS, FARSA in humans) has an additional activity that promotes growth and proliferation (Ho et al., 2021). Here we describe that elevated -PheRS levels also induce cell fate changes and tumorous phenotypes in Drosophila midguts. Excessive proliferating cells with stem and progenitor cell characteristics accumulate and the composition of the terminally differentiated cells changes, too. This phenotype together with observed genetic interactions between -PheRS and Notch levels show that -PheRS counteracts Notch signaling in many different tissues and developmental stages. This novel activity of -PheRS maps to its N-terminal part, which is naturally produced. The fragment contains a DNA binding domain, translocates into nuclei, and displays essential similarities to a Notch domain that binds to the downstream transcription factor. This suggests that it might be competing with Notch for binding to a common target. Not only because Notch plays important roles in many tumors, but also because FARSA mRNA levels are considerably upregulated in many tumors, this novel activity deserves more attention for cancer research. 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. Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent global expansion A novel coronavirus, SARS-CoV-2, has caused over 190 million cases and over 4 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e. M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic one year later. By analyzing 4,013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of initial phase of global ignition over a timespan of 14-week. 11 major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted a half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implying that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261,323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e. including all the publicly available viral genomes) in order to recapitulate our findings over one-year timespan. By 25 December 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the four variants of concern (VOC) are the descendants of WE1 type. This study demonstrates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time. High performance single-cell gene regulatory network inference atscale: The Inferelator 3.0 MotivationGene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

ResultsIn this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AvailabilityThe inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/). Global structural stability and the role of cooperation in mutualistic system Dynamical systems on graphs allow to describe multiple phenomena from different areas of Science. In particular, many complex systems in Ecology are studied by this approach. In this paper we analize the mathematical framework for the study of the structural stability of each stationary point, feasible or not, introducing a generalization for this concept, defined as Global Structural Stability. This approach would fit with the proper mathematical concept of structural stability, in which we find a full description of the complex dynamics on the phase space due to nonlinear dynamics. This fact can be analyzed as an informational field grounded in a global attractor whose structure can be completely characterized. These attractors are stable under perturbation and suppose the minimal structurally stable sets. We also study in detail, mathematically and computationally, the zones characterizing different levels of biodiversity in bipartite graphs describing mutualistic antagonistic systems of population dynamics. In particular, we investigate the dependence of the region of maximal biodiversity of a system on its connectivity matrix. On the other hand, as the network topology does not completely determine the robustness of the dynamics of a complex network, we study the correlation between structural stability and several graph measures. A systematic study on synthetic and biological graphs is presented, including 10 mutualistic networks of plants and seed-dispersal and 1000 random synthetic networks. We compare the role of centrality measures and modularity, concluding the importance of just cooperation strength among nodes when describing areas of maximal biodiversity. Indeed, we show that cooperation parameters are the central role for biodiversity while other measures act as secondary supporting functions.

Author summaryWe introduce the concept of Global Structural Stability as a proper mathematical concept to fully understand biodiversity in some ecological systems. Our concept retakes the definitions in the classical works of R. Thom [1] and Andronov-Pontryagin [2]. Moreover, there exists a close relation between the structure of a complex network, described as a graph, and its associated dynamics. Mutualistic networks introduce cooperation links between two groups of species, as plant and pollinators or seed-dispersal. The understanding of organizational aspects leading to maximizing biodiversity is one of the more important research areas in Theoretical and Applied Ecology. In this work we introduce a systematic study on different graph measures in order to identify optimal organization for maximal biodiversity (defined as structural stability). Our results conclude that, for mutualistic systems, the strength in cooperation parameters are the core fact, i.e., cooperation is the real fact optimizing biodiversity among other possible structural configurations. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. PRSS2 stimulates tumor growth by remodeling the TME via repression of Tsp1 In the earliest stages of tumor development, epithelial tumors (carcinomas) are physically confined to the area of the tissue in which they form. These nascent lesions (carcinomas in situ) are sequestered from the tissue parenchyma by the basement membrane. Within the tissue parenchyma lie a myriad of cell types comprised of fibroblasts, immune and inflammatory cells and endothelial cells. Upon invasion across the basement membrane and into the tissue parenchyma, tumors must manipulate the expression of pro- and anti-tumorigenic proteins such that pro-tumorigenic factors are produced in excess to anti-tumorigenic proteins. One such anti-tumorigenic protein is Thrombospondin-1 (Tsp-1). We have previously demonstrated that stimulation of Tsp-1 in the tumor microenvironment (TME) potently inhibits tumor growth and progression and in some cases induces tumor regression. Here, we identify a novel tumor-mediated mechanism to repress the expression of Tsp-1 in the TME via secretion of the serine protease PRSS2. We demonstrate that PRSS2 represses Tsp-1, not via its enzymatic activity, but by binding to low-density lipoprotein receptor-related protein 1 (LRP1). These findings describe a novel activity for PRSS2 through binding to LRP1 and represent a potential therapeutic strategy to treat cancer by blocking the PRSS2-mediated repression of Tsp-1.