NOMA: a high-throughput microchip for robust, sequential measurements of secretions from the same single-cells Despite recent advances in single-cell analysis technologies, lacking simple methods to keep the live single-cells traceable for longitudinal detection reliably poses a significant obstacle in single-cell secretion analysis. Here we reported the high-density NOMA (narrow-opening microwell array) microchip that realized the retention of[&ge;] 97% of trapped single cells in dedicated spatial locations during repetitive detection procedures, verified with both adherent and suspension cells by two researchers independently. We applied it in monitoring single-cell protein secretions sequentially from the same single cells, and we found the digital protein secretion patterns dominate the protein secretion. We also demonstrated the microchip for longitudinally tracking IL-8 and the CD81+EV secretions from the same single-cells over days, which revealed the presence of "super secretors" within the cell population be more persistent to secrete protein or extracellular vesicle for an extended period. The NOMA platform reported here is simple, robust, and easy to operate for tracking sequential measurements from the same single cells, representing a novel and informative tool to inspire new observations in biomedical research.

Table of Contents Graphic

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=66 SRC="FIGDIR/small/427049v4_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@aa69c7org.highwire.dtl.DTLVardef@6bc79org.highwire.dtl.DTLVardef@3fcbbforg.highwire.dtl.DTLVardef@597f4b_HPS_FORMAT_FIGEXP  M_FIG C_FIG
Genomically Mined Acoustic Reporter Genes Enable In Vivo Monitoring of Tumors and Tumor-Homing Bacteria A major outstanding challenge in the fields of biological research, synthetic biology and cell-based medicine is the difficulty of visualizing the function of natural and engineered cells noninvasively inside opaque organisms. Ultrasound imaging has the potential to address this challenge as a widely available technique with a tissue penetration of several centimeters and spatial resolution below 100 {micro}m. Recently, the first genetically encoded acoustic reporters were developed based on bacterial gas vesicles to link ultrasound signals to molecular and cellular function. However, the properties of these first-generation acoustic reporter genes (ARGs) resulted in limited sensitivity and specificity for imaging in the in vivo context. Here, we describe second-generation ARGs with greatly improved acoustic properties and expression characteristics, identified through a phylogenetic screen of candidate gene clusters from diverse bacteria and archaea. The resulting constructs offer major qualitative and quantitative improvements, including much stronger ultrasound contrast, the ability to produce nonlinear signals distinguishable from background tissue, and stable long-term expression. We demonstrate the utility of these next-generation ARGs by imaging in situ gene expression in mouse models of breast cancer and tumor-homing probiotic bacteria, revealing the unique spatial distributions of tumor growth and colonization by therapeutic cells noninvasively in living subjects.
Combined epigenetic and metabolic inhibition blocks platinum-induced ovarian cancer stem cell enrichment High grade serous ovarian cancer (HGSOC) is the most common and aggressive type of ovarian cancer. Platinum resistance is a common occurrence in HGSOC and a main cause of tumor relapse resulting in high patient mortality rates. Recurrent OC is enriched in aldehyde dehydrogenase (ALDH)+ ovarian cancer stem cells (OCSCs), which are resistant to platinum agents. We demonstrated that acute platinum treatment induced a DNA damage-dependent decrease in BRCA1 levels. In a parallel response associated with G2/M arrest, platinum treatment also induced an increase in expression of NAMPT, the rate limiting regulator of NAD+ production from the salvage pathway, and levels of NAD+, the cofactor required for ALDH activity. Concurrent inhibition of DNA methyltransferases (DNMTs) and NAMPT synergistically abrogated the platinum-induced increase in OCSCs. Combining pharmacological inhibitors of DNMT and NAMPT with carboplatin reduced tumorigenesis and OCSC percentage in vivo. We conclude that both epigenetic and metabolic alterations lead to platinum induced OCSC enrichment, providing preclinical evidence that in the neoadjuvant setting, combining DNMT and NAMPT inhibitors with platinum has the potential to reduce OC recurrence and avert the development of platinum resistance. Integrative Modelling of Signalling Network Dynamics Identifies Cell Type-selective Therapeutic Strategies for FGFR4-driven Cancers Oncogenic FGFR4 signalling represents a potential therapeutic target in many cancer types, including triple negative breast cancer (TNBC) and hepatocellular carcinoma (HCC). However, resistance to single-agent therapy directed at FGFR4 remains a major challenge, prompting the need to identify more effective combinatorial therapeutic strategies. Here, we integrated computational network modelling and experimental validation to characterise dynamic reprogramming of the FGFR4 signalling network in TNBC following FGFR4 kinase inhibition. We found that AKT, which signals downstream of FGFR4, displayed a rapid and potent reactivation following FGFR4 targeting. Through model-based simulation and systematic prediction of the effect of co-targeting specific network nodes, we predicted, and validated experimentally, strong synergism of co-targeting FGFR4 and particular ErbB kinases or AKT, but not the upstream kinase PI3K. Further, incorporation of protein expression data from hundreds of cancer cell types enabled us to adapt our model to other diverse cellular contexts, leading to the prediction that while AKT rebound occurs frequently, it is not a general phenomenon. Instead, ERK is reactivated in a subset of cell types, including the FGFR4-driven HCC cell line Hep3B. This was subsequently corroborated, and moreover, co-targeting FGFR4 and MEK in Hep3B cells markedly enhanced inhibition of cell proliferation. Overall, these findings provide novel insights into the dynamics of drug-induced network remodelling in cancer cells, highlight the impact of protein expression heterogeneity on network response to targeted therapy and identify candidate cell type-selective combination treatments for FGFR4-driven cancer. Network Analysis Reveals Different Cellulose Degradation Strategies across Trichoderma harzianum Strains Associated with XYR1 and CRE1 Trichoderma harzianum, whose gene expression is tightly controlled by the transcription factors (TFs) XYR1 and CRE1, is a potential candidate for hydrolytic enzyme production. Here, we performed a network analysis of T. harzianum IOC-3844 and T. harzianum CBMAI-0179 to explore how the regulation of these TFs varies between these strains. In addition, we explored the evolutionary relationships of XYR1 and CRE1 protein sequences among Trichoderma spp. The results of the T. harzianum strains were compared with those of Trichoderma atroviride CBMAI-0020, a mycoparasitic species. Although transcripts encoding carbohydrate-active enzymes (CAZymes), TFs, transporters, and proteins with unknown functions were coexpressed with cre1 or xyr1, other proteins indirectly related to cellulose degradation were identified. The enriched GO terms describing the transcripts of these groups differed across all strains, and several metabolic pathways with high similarity between both regulators but strain-specific differences were identified. In addition, the CRE1 and XYR1 subnetworks presented different topology profiles in each strain, likely indicating differences in the influences of these regulators according to the fungi. The hubs of the cre1 and xyr1 groups included transcripts not yet characterized or described as being related to cellulose degradation. The first-neighbor analyses confirmed the results of the profile of the coexpressed transcripts in cre1 and xyr1. The analyses of the shortest paths revealed that CAZymes upregulated under cellulose degradation conditions are most closely related to both regulators, and new targets between such signaling pathways were discovered. Although the evaluated T. harzianum strains are phylogenetically close and their amino acid sequences related to XYR1 and CRE1 are very similar, the set of transcripts related to xyr1 and cre1 differed, suggesting that each T. harzianum strain used a specific regulation strategy for cellulose degradation. More interestingly, our findings may suggest that XYR1 and CRE1 indirectly regulate genes encoding proteins related to cellulose degradation in the evaluated T. harzianum strains. An improved understanding of the basic biology of fungi during the cellulose degradation process can contribute to the use of their enzymes in several biotechnological applications and pave the way for further studies on the differences across strains of the same species. Topaz1, an essential gene for murine spermatogenesis, down-regulates the expression of numerous testis-specific long non-coding RNAs. Spermatogenesis involves coordinated processes, including meiosis, to produce functional gametes. We previously reported Topaz1 as a germ cell-specific gene highly conserved in vertebrates. Topaz1 knockout males are sterile with testes that lack haploid germ cells because of meiotic arrest after prophase I. To better characterize Topaz1-/- testes, we used RNA-sequencing analyses at two different developmental stages (P16 and P18). The absence of TOPAZ1 disturbed the expression of genes involved in microtubule and/or cilium mobility, which was consistent with testicular histology showing the disruption of microtubules and centrosomes. Moreover, a quarter of P18 dysregulated genes are long non-coding RNAs (lncRNAs), and three of them are testis-specific and located in spermatocytes, their expression starting between P11 and P15. The suppression of one of them, 4939463O16Rik, did not alter fertility although sperm parameters were disturbed and sperm concentration fell. The transcriptome of P18-4939463O16Rik-/- testes was altered and the molecular pathways affected included microtubule-based processes, the regulation of cilium movement and spermatogenesis. The absence of TOPAZ1 protein or 4930463O16Rik produced the same enrichment clusters in mutant testes despite a contrasted phenotype on male fertility. In conclusion, TOPAZ1 appeared to stabilize the expression of numerous lncRNAs. Its suppression is not essential for fertility but required during the terminal differentiation of male gametes.

Author SummaryThe Topaz1 gene was initially characterized during the initiation of meiosis in the sheep fetal ovary. In order to determine its function, a KO of the murine gene was performed. In this species, only males were sterile and spermatogenesis was blocked before the first meiotic division. Here, we show that cytoskeletal elements are markedly disturbed in mutant testes, indicating that these elements play an important function in spermatogenesis. While the mitotic spindle of spermatogonia was normal, the meiotic spindle of spermatocytes was hemi-spindle-shaped and the homologous chromosome pairs could position themselves on the equatorial plate. In addition, lncRNAs account for 25% of genes whose expression in testes varies significantly in the absence of Topaz1. This suggests a key role for these factors in spermatogenesis. Largely testis-specific, they may be involved in spermatogenesis and play a more or less critical role in mouse fertility, which probably also depends on their redundancies. Mouse oocytes do not contain a Balbiani body Oocytes spend the majority of their lifetime in a primordial, dormant state. Unlike many somatic cell types and mature oocytes, the cellular and molecular biology of primordial oocytes is largely unexplored. Yet, studying their cellular biology is necessary to understand the mechanisms through which oocytes maintain cellular fitness for decades, and why they eventually fail with age.

A hallmark of primordial oocytes in many species is the Balbiani body, a non-membrane bound compartment that contains the majority of mitochondria in the oocyte cytoplasm. The Balbiani body has been proposed to be essential for maintaining mitochondria in a healthy state during long-lasting dormancy, however, the architecture and function of the mammalian Balbiani body remains unknown.

Here, we develop enabling methods for live-imaging based comparative characterization of Xenopus, mouse and human primordial oocytes. We show that primordial oocytes in all three vertebrate species contain active mitochondria, Golgi apparatus and lysosomes. We further demonstrate that human and Xenopus oocytes have a Balbiani body characterized by a dense accumulation of mitochondria in their cytoplasm. However, despite previous reports, we did not find a Balbiani body in mouse oocytes. Instead, we demonstrate what was previously used as a marker for the Balbiani body in mouse primordial oocytes is in fact a ring-shaped Golgi apparatus that is not functionally associated with oocyte dormancy. Our work provides the first insights into the organisation of the cytoplasm in mammalian primordial oocytes, and clarifies relative advantages and limitations of choosing different model organisms for studying oocyte dormancy. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Mapping densely packed GPIIb/IIIa receptors in murine blood platelets with expansion microscopy Interrogating small platelets and their densely packed, highly abundant receptor landscape is key to understand platelet clotting, a process that can save lives when stopping blood loss after an injury, but also kill when causing heart attack, stroke or pulmonary embolism. The underlying key receptor distributions and interactions, in particular the relevance of integrin clustering, are not fully understood is because of highly abundant and densely distributed GPIIb/IIIa receptors. This makes receptor distributions difficult to assess even by super-resolution fluorescence microscopy. Here, we combine dual-color expansion and confocal microscopy with colocalization analysis to assess platelet receptor organization without the need of a super-resolution microscope. We show that 4x expansion is highly straight-forward for super-resolution microscopy of platelets, while 10x expansion provides higher precision at the price of increased efforts in sample preparation and imaging. Quantifying various receptor colocalization scenarios we demonstrate that expansion microscopy can pinpoint receptor distributions and interactions in resting and activated platelets being superior to conventional methods that fail in such dense 3D scenarios with highly abundant receptors. We reveal the presence of GPIIb/IIIa clusters in resting platelets, which are not affected by platelet activation indicating that they contribute to the rapid platelet response during platelet clotting.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=80 SRC="FIGDIR/small/431449v3_ufig1.gif" ALT="Figure 1">
View larger version (29K):
org.highwire.dtl.DTLVardef@c4e38forg.highwire.dtl.DTLVardef@37c1adorg.highwire.dtl.DTLVardef@8c8597org.highwire.dtl.DTLVardef@196bd69_HPS_FORMAT_FIGEXP  M_FIG C_FIG Gene expression underlying parenting and being parented reflect canalization, not plasticity, in response to environmental stress Parenting buffers offspring from hostile environments, but it is not clear how or if the genes that underlie parenting change their expression under environmental stress. We recently demonstrated that for the subsocial carrion beetle, Nicrophorus orbicollis, temperature during parenting does not affect parenting phenotypes. Here, we ask if transcriptional changes associated with parenting are likewise robust to environmental stress. The absence of a transcriptional response for parenting under stress would suggest that the genetic  programs for parenting and being parented are canalized. Conversely, a robust transcriptional response would suggest that plasticity of underlying gene expression is critical for maintaining behavioral stability, and that these mechanisms provide a potential target for selection in the face of environmental change. We test these alternatives by characterizing gene expression of parents and offspring with and without parent-offspring interactions under a benign and a stressful temperature. We found that parent-offspring interactions elicit distinct transcriptional responses of parents and larvae irrespective of temperature. We further detected robust changes of gene expression in beetles breeding at 24{degrees} C compared to 20{degrees} C irrespective of family interaction. However, no strong interaction between parent-offspring interaction and temperature was detected for either parents or larvae. We therefore conclude that canalization, not plasticity of gene expression, most likely explains the absence of behavioral plasticity under thermal stress. This result suggests that species may not have the genetic variation needed to respond to all environmental change, especially for complex phenotypes.

IMPACT STATEMENTParenting is thought to evolve to buffer offspring from harsh or inhospitable conditions. However, once evolved, does parenting continue to buffer against new environmental challenges such as might occur under global warming? One way organisms achieve different phenotypes without evolving different genes for each trait is by changing timing or extent of gene expression. Changes in gene expression influencing physiological traits such as hypoxia have been shown to facilitate adaptation to changing environments. It is less clear if changes in gene expression are also important for behavioral plasticity in complex social traits, such as parental care. In this study, we ask whether families of the subsocial burying beetle, Nicrophorus orbicollis, show altered expression of parent or offspring genes that may help compensate for reduced offspring success under global warming. To test the extent that gene expression underlies behavioral plasticity in parenting and offspring responses to parenting, we characterize the gene expression profiles of whole families (males, females, and larvae) with and without parent-offspring interactions under two thermal conditions: 20{degrees} C and 24{degrees} C. We found that heat stress predominantly induced downregulation of genes, consistent with a constrained acclimation response. However, when temperature was considered in conjunction with parent-offspring interactions, genetic plasticity was extremely limited. That is, the genes underlying parent-offspring interactions did not show strong signatures of differential expression under more stressful thermal conditions. Given that we previously found no effects of temperature on parental behavior, we suggest that constrained, not plastic gene expression, underpins stability of parent-offspring interactions and limits the buffering capacity for family interactions under environmental stress. Insights from this study deepen our understanding of the interplay between genetic and phenotypic plasticity in complex behavioral traits and how traits such as parenting may limit responses to novel environmental challenges. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. 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. Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.) Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Significance StatementUsing traditional and automated/high-throughput (using a novel deep learning approach) phenotyping methods we studied leaf anatomical variation in sunflower. Genome-wide association (GWA) analyses identified numerous genomic regions underlying individual trait variation and regions underlying major multivariate axes of phenotypic variation. These results illustrate the value of employing a multivariate approach to GWA analyses and shed light on the extent to which leaf trait (co-)variation can be genetically decoupled to explore novel phenotypic space. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. CDK9 and PP2A regulate RNA polymerase II transcription termination and coupled RNA maturation. CDK9 is a critical kinase required for the productive transcription of protein-coding genes by RNA polymerase II (pol II). As part of P-TEFb, CDK9 phosphorylates the carboxyl-terminal domain (CTD) of pol II and elongation factors, including SPT5, which allows pol II to elongate past the early elongation checkpoint (EEC) encountered soon after initiation. We show that, in addition to halting pol II at the EEC, loss of CDK9 activity causes premature termination of transcription across the last exon, loss of polyadenylation factors from chromatin, and loss of polyadenylation of nascent transcripts. Inhibition of the phosphatase PP2A abrogates the premature termination and loss of polyadenylation caused by CDK9 inhibition, indicating that this kinase/phosphatase pair regulates transcription elongation and RNA processing at the end of protein-coding genes. Our phosphoproteomic analyses after CDK9 inhibition, using either DRB or an analog-sensitive CDK9 cell line, confirm the splicing factor SF3B1 as an additional key target of this kinase. These results emphasize the important roles that CDK9 plays in coupling transcription elongation and termination to RNA maturation downstream of the EEC. The maize gene maternal derepression of r1 (mdr1) encodes a DNA glycosylase that demethylates DNA and reduces siRNA expression in endosperm Demethylation of transposons can activate expression of nearby genes and cause imprinted gene expression in endosperm, and it is hypothesized to lead to expression of transposon siRNAs that reinforce silencing in the next generation through transfer either into egg or embryo. Here we describe maternal derepression of r1 (mdr1), which encodes a DNA glycosylase with homology to Arabidopsis DEMETER and which is partially responsible for demethylation of thousands of regions in endosperm. Instead of promoting siRNA expression in endosperm, MDR1 activity inhibits it. Methylation of most repetitive DNA elements in endosperm is not significantly affected by MDR1, with an exception of Helitrons. While maternally-expressed imprinted genes preferentially overlap with MDR1 demethylated regions, the majority of genes that overlap demethylated regions are not imprinted. Double mutant megagametophytes lacking both MDR1 and its close homolog DNG102 result in early seed failure, and double mutant microgametophytes fail pre-fertilization. These data establish DNA demethylation by glycosylases as essential in maize endosperm and pollen and suggest that neither transposon repression nor genomic imprinting are its main function in endosperm. Synergistic block of SARS-CoV-2 infection by combined drug inhibition of the host entry factors PIKfyve kinase and TMPRSS2 protease Repurposing FDA-approved inhibitors able to prevent infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) could provide a rapid path to establish new therapeutic options to mitigate the effects of coronavirus disease 2019 (COVID-19). Proteolytic cleavages of the spike S protein of SARS-CoV-2, mediated by the host cell proteases cathepsin and TMPRSS2, alone or in combination, are key early activation steps required for efficient infection. The PIKfyve kinase inhibitor apilimod interferes with late endosomal viral traffic, and through an ill-defined mechanism prevents in vitro infection through late endosomes mediated by cathepsin. Similarly, inhibition of TMPRSS2 protease activity by camostat mesylate or nafamostat mesylate prevents infection mediated by the TMPRSS2-dependent and cathepsin-independent pathway. Here, we combined the use of apilimod with camostat mesylate or nafamostat mesylate and found an unexpected ~5-10-fold increase in their effectiveness to prevent SARS-CoV-2 infection in different cell types. Comparable synergism was observed using both, a chimeric vesicular stomatitis virus (VSV) containing S of SARS-CoV-2 (VSV-SARS-CoV-2) and SARS-CoV-2 virus. The substantial ~5-fold or more decrease of half maximal effective concentrations (EC50 values) suggests a plausible treatment strategy based on the combined use of these inhibitors.

IMPORTANCEInfection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is causing the coronavirus disease 2019 (COVID-2019) global pandemic. There are ongoing efforts to uncover effective antiviral agents that could mitigate the severity of the disease by controlling the ensuing viral replication. Promising candidates include small molecules that inhibit the enzymatic activities of host proteins, thus preventing SARS-CoV-2 entry and infection. They include Apilimod, an inhibitor of PIKfyve kinase and camostat mesylate and nafamostat mesylate, inhibitors of TMPRSS2 protease. Our research is significant for having uncovered an unexpected synergism in the effective inhibitory activity of apilimod used together with camostat mesylate or with nafamostat mesylate. 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. Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. Automated methodology for optimal selection of the minimum electrode subset for accurate EEG source estimation based on Genetic Algorithm optimization High-density Electroencephalography (HD-EEG) has been proven to be the most accurate option to estimate the neural activity inside the brain. Multiple studies report the effect of electrode number on source localization for specific sources and specific electrode configurations. The electrodes for these configurations have been manually selected to uniformly cover the entire head, going from 32 to 128 electrodes, but electrodes were not selected according to their contribution to accuracy. In this work, an optimization-based study is proposed to determine the minimum number of electrodes and identify optimal combinations of electrodes that can retain the localization accuracy of HD-EEG reconstructions. This optimization approach incorporates scalp landmark positions of widely used EEG montages. In this way, a systematic search for the minimum electrode subset is performed for single and multiple source localization problems. The Non-dominated Sorting Genetic Algorithm II (NSGA-II) combined with source reconstruction methods is used to formulate a multi-objective optimization problem that concurrently minimizes 1) the localization error for each source and 2) the number of required EEG electrodes. The method can be used for evaluating the source localization quality of low-density EEG systems (e.g. consumer-grade wearable EEG). We performed an evaluation over synthetic and real EEG dataset with known ground-truth. The experimental results show that optimal subsets with 6 electrodes can obtain an equal or better accuracy than HD-EEG (with more than 200 channels) for a single source case. This happened when reconstructing a particular brain activity in more than 88% of the cases (in synthetic signals) and 63% (in real signals), and in more than 88% and 73% of the cases when considering optimal combinations with 8 channels. For a multiple-source case of three sources (only with synthetic signals), it was found that optimized combinations of 8, 12 and 16 electrodes attained an equal or better accuracy than HD-EEG with 231 electrodes in at least 58%, 76%, and 82% of the cases respectively. Additionally, for such electrode numbers, a lower mean error and standard deviation than with 231 electrodes were obtained.

HighlightsO_LIThe number of EEG electrodes and their locations can be optimized for reconstructing the brain source activity.
C_LIO_LIOptimally selected EEG electrodes can retain the accuracy of high density montages (256, 128 chs) for brain source estimation, when electrodes are selected according to their contribution to accuracy.
C_LIO_LIWith optimization, selected combinations of EEG electrodes will flexibilize the estimation of the source activity.
C_LI Regulation of the PKD2 channel function by TACAN Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in membrane receptor PKD1 or cation channel PKD2. TACAN (also named TMEM120A), recently reported as an ion channel in neuron cells for mechano and pain sensing, is also distributed in diverse non-neuronal tissues such as kidney, heart and intestine, suggesting its involvement in other functions. In this study, we found that TACAN is in complex with PKD2 in native renal cell lines. Using the two-electrode voltage clamp in Xenopus oocytes we found that TACAN inhibited the channel activity of PKD2 gain-of-function mutant F604P. The first and last transmembrane domains of TACAN were found to interact with the PKD2 C-and N-terminal portions, respectively. We showed that the TACAN N-terminus acted as a blocking peptide and that TACAN inhibits the PKD2 function through the PKD2/TACAN binding. By patch clamping in mammalian cells, we found that TACAN inhibits both the single channel conductance and open probability of PKD2 and mutant F604P. PKD2 co-expressed with TACAN, but not PKD2 alone, exhibited pressure sensitivity. Furthermore, we also found that TACAN aggravates PKD2-dependent tail curvature and pronephric cysts in larval zebrafish, in support of the in vitro inhibitory effects of TACAN. In summary, this study revealed that TACAN acts as a PKD2 inhibitor and mediates mechano sensitivity of the PKD2/TACAN channel complex. Extending Bayesian Elo-rating to quantify the steepness of dominance hierarchies The steepness of dominance hierarchies provides information about the degree of competition within animal social groups and is thus an important concept in socioecology. The currently most widely-used metrics to quantify steepness are based on Davids scores (DS) derived from dominance interaction networks. One serious drawback of these DS-based metrics is that they are biased, i.e., network density systematically decreases steepness values. Here, we provide a novel approach to estimate steepness based on Elo-ratings, implemented in a Bayesian framework (STEER: Steepness estimation with Elo-rating). Our new metric has two key advantages. First, STEER is unbiased, precise and more robust to data density than DS-based steepness. Second, it provides explicit probability distributions of the estimated steepness coefficient, which allows uncertainty assessment. In addition, it relies on the same underlying concept and is on the same scale as the original measure, and thus allows comparison to existing published results. We evaluate and validate performance of STEER by means of experimentation on empirical and artificial data sets and compare its performance to that of several other steepness estimators. Our results suggest that STEER provides a considerable improvement over existing methods to estimate dominance hierarchy steepness. We provide an R package EloSteepness to calculate the new steepness measure, and also show an example of using steepness in a comparative analysis. Severe COVID-19 associated variants linked to chemokine receptor gene control in monocytes and macrophages Genome-wide association studies have identified 3p21.31 as the main risk locus for severe disease in COVID-19 patients, although underlying biological mechanisms remain elusive. We performed a comprehensive epigenomic dissection of the 3p21.31 locus, identifying a CTCF-dependent tissue-specific 3D regulatory chromatin hub that controls the activity of several tissue-homing chemokine receptor (CCR) genes in monocytes and macrophages. Risk SNPs colocalized with regulatory elements and were linked to increased expression of CCR1, CCR2 and CCR5 in monocytes and macrophages. As excessive organ infiltration of inflammatory monocytes and macrophages is a hallmark of severe COVID-19, our findings provide a rationale for the genetic association of 3p21.31 variants with elevated risk of hospitalization upon SARS-CoV-2 infection. Confronting population models with experimental microcosm data: from trajectory matching to state-space models O_LIPopulation and community ecology traditionally has a very strong theoretical foundation with well-known models, such as the logistic and its many variations, and many modification of the classical Lotka-Volterra predator-prey and interspecific competition models. More and more, these classical models are confronted to data via fitting to empirical time-series for purposes of projections or for estimating model parameters of interest. However, the interface between mathematical population or community models and data, provided by a statistical model, is far from trivial.
C_LIO_LIIn order to help empiricists, especially researchers working with experimental laboratory populations in micro- and mesocosms, make informed decisions, we here ask which error structure one should use when fitting classical deterministic ODE models to empirical data, from single species to community dynamics and trophic interactions. We use both realistically simulated data and empirical data from microcosms to answer this question in a Bayesian framework.
C_LIO_LIWe find that many model parameters can be estimated precisely with an appropriate choice of error structure using pure observation error or state-space models, if observation errors are not too high. However, Allee effect models are typically hard to identify and state-space models should be preferred with when model complexity increases.
C_LIO_LIOur work shows that, at least in the context of experimental laboratory populations, deterministic models can be used to describe stochastic population dynamics that include process variability and observation error. Also, many scenarios do not require a complex state-space model formulation and simpler trajectory matching is sufficient for accurate parameter estimates. Finally, we provide a comprehensive tutorial for fitting these models in R.
C_LI