Innovative GenExpA software for selecting suitable reference genes for reliable normalization of gene expression in melanoma The algorithms commonly used to select the best stable reference gene in RT-qPCR data analysis have their limitations. We showed that simple selection of the reference gene or pair of genes with the lowest stability value from the pool of potential reference genes - a commonly used approach - is not sufficient to accurately and reliably normalize the target gene transcript and can lead to biologically incorrect conclusions. For reliable assessment of changes in a target gene expression level, we propose our innovative GenExpA software, which works in a manner independent of the experimental model and the normalizer used. GenExpA software selects the best reference by combining the NormFinder algorithm with progressive removal of the least stable gene from the candidate genes in a given experimental model and in the set of daughter models assigned to it. The reliability of references is validated based on the consistency of the statistical analyses of normalized target gene expression levels through all models, described by the coherence score (CS). The use of the CS value imparts a new quality to qPCR analysis because it clarifies how low the stability value of reference must be in order for biologically correct conclusions to be drawn. We tested our method on qPCR data for the B4GALT genes family in melanoma, which is characterized by a high mutation rate, and in melanocytes. GenExpA is available at https://github.com/DorotaHojaLukowicz/GenExpA or https://www.sciencemarket.pl/baza-programow-open-source#oferty.

HighlightsGenExpA - next-generation software for normalizer selection and validation

The GenExpA tool defines how low the stability value of the reference should be

The GenExpA tool determines the level of gene expression analysis robustness

The GenExpA tool increases the speed of analysis and reduces cost
Dose-dependent thresholds of dexamethasone destabilize CAR T-cell treatment efficacy Chimeric antigen receptor (CAR) T-cell therapy is potentially an effective targeted immunotherapy for glioblastoma, yet there is presently little known about the efficacy of CAR T-cell treatment when combined with the widely used anti-inflammatory and immunosuppressant glucocorticoid, dexamethasone. Here we present a mathematical model-based analysis of three patient-derived glioblastoma cell lines treated in vitro with CAR T-cells and dexamethasone. Advanced in vitro experimental cell killing assay technologies allow for highly resolved temporal dynamics of tumor cells treated with CAR T-cells and dexamethasone, making this a valuable model system for studying the rich dynamics of nonlinear biological processes with translational applications. We model the system as a nonautonomous, two-species predator-prey interaction of tumor cells and CAR T-cells, with explicit time-dependence in the clearance rate of dexamethasone. Using time as a bifurcation parameter, we show that (1) dexamethasone destabilizes coexistence equilibria between CAR T-cells and tumor cells in a dose-dependent manner and (2) as dexamethasone is cleared from the system, a stable coexistence equilibrium returns in the form of a Hopf bifurcation. With the model fit to experimental data, we demonstrate that high concentrations of dexamethasone antagonizes CAR T-cell efficacy by exhausting, or reducing the activity of CAR T-cells, and by promoting tumor cell growth. Finally, we identify a critical threshold in the ratio of CAR T-cell death to CAR T-cell proliferation rates that predicts eventual treatment success or failure that may be used to guide the dose and timing of CAR T-cell therapy in the presence of dexamethasone in patients.

Author summaryBioengineering and gene-editing technologies have paved the way for advance immunotherapies that can target patient-specific tumor cells. One of these therapies, chimeric antigen receptor (CAR) T-cell therapy has recently shown promise in treating glioblastoma, an aggressive brain cancer often with poor patient prognosis. Dexamethasone is a commonly prescribed anti-inflammatory medication due to the health complications of tumor associated swelling in the brain. However, the immunosuppressant effects of dexamethasone on the immunotherapeutic CAR T-cells are not well understood. To address this issue, we use mathematical modeling to study in vitro dynamics of dexamethasone and CAR T-cells in three patient-derived glioblastoma cell lines. We find that in each cell line studied there is a threshold of tolerable dexamethasone concentration. Below this threshold, CAR T-cells are successful at eliminating the cancer cells, while above this threshold, dexamethasone critically inhibits CAR T-cell efficacy. Our modeling suggests that in the presence of high dexamethasone reduced CAR T-cell efficacy, or increased exhaustion, can occur and result in CAR T-cell treatment failure.
Loss of histone methyltransferase SETD1B in oogenesis results in the redistribution of genomic histone 3 lysine 4 trimethylation Histone 3 lysine 4 trimethylation (H3K4me3) is an epigenetic mark found at gene promoters and CpG islands. H3K4me3 is essential for mammalian development, yet mechanisms underlying its genomic targeting are poorly understood. H3K4me3 methyltransferases SETD1B and MLL2 are essential for oogenesis. We investigated changes in H3K4me3 in Setd1b conditional knockout (cKO) oocytes using ultra-low input ChIP-seq, with comparisons to DNA methylation and gene expression analyses. H3K4me3 was redistributed in Setd1b cKO oocytes showing losses at active gene promoters associated with downregulated gene expression. Remarkably, many regions also gained H3K4me3, in particular those that were DNA hypomethylated, transcriptionally inactive and CpG-rich, which are hallmarks of MLL2 targets. Consequently, loss of SETD1B disrupts the balance between MLL2 and de novo DNA methyltransferases in determining the epigenetic landscape during oogenesis. Our work reveals two distinct, complementary mechanisms of genomic targeting of H3K4me3 in oogenesis, with SETD1B linked to gene expression and MLL2 to CpG content.

Graphical AbstractIn oogenesis, SETD1B and CXXC1 target H3K4me3 to actively transcribed gene promoters, while MLL2 targets transcriptionally inactive regions based on underlying CpG composition (upper panel). When SETD1B is ablated, H3K4me3 is lost at a subset of active promoters, resulting in downregulation of transcription (lower panel). Loss of SETD1B alters the activity of MLL2, permitting MLL2 to deposit H3K4me3 at CpG-rich regions, many of which should otherwise be DNA methylated. Thus, it is evident that MLL2 and de novo DNMTs compete for genomic occupancy late in oogenesis, and loss of SETD1B disrupts the balance of these mechanisms.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. wnt16 exerts pleiotropic effects on bone and lean mass in zebrafish Osteoporosis, a disease of bone fragility, and sarcopenia, a condition of reduced muscle mass and strength, commonly occur in the same individual--a condition termed osteosarcopenia. Genetic variants at the CPED1-WNT16 locus have pleiotropic effects on human bone mineral density (BMD) and lean tissue mass. However, our current understanding of Wnt16 and its influence on bone fails to explain effects on lean mass at this locus. In this study, we sought to better understand the biological mechanism underlying pleiotropy at the CPED1-WNT16 locus on BMD and lean mass. Using single cell analysis, microCT imaging, and genetic approaches in zebrafish, we show that wnt16 exerts dual influence on muscle and bone, and provide evidence that this occurs through a function in somite development, a period in which precursors for both of these tissues are specified. We show that mutations in wnt16 alter vertebral size and lean mass during axial skeletogenesis. We show that wnt16+ cells transcriptionally resemble cells within the amniote dermomyotome, and provide evidence they function as myogenic precursors during embryonic development. We also provide evidence that WNT16 is a gene of major effect on lean mass at the CPED-WNT16 locus. Our studies support the potential for wnt16 to exert pleiotropic effects on bone and lean mass through a function in myogenic precursors. Antibody-induced glomerulonephritis is amplified by RTEC-intrinsic IL-17 signaling and restrained by IL-17-mediated induction of the endoribonuclease Regnase-1 (Zc3h12a) Antibody-mediated glomerulonephritis (AGN) is a clinical manifestation of many autoimmune kidney diseases for which few effective treatments exist. Chronic inflammatory circuits in renal glomerular and tubular cells lead to tissue damage in AGN. These cells are targeted by the cytokine IL-17, which has recently been shown to be a central driver of the pathogenesis of AGN. However, surprisingly little is known about the regulation of pathogenic IL-17 signaling in the kidney. Here, using a well characterized mouse model of AGN, we show that IL-17 signaling in renal tubular epithelial cells (RTECs) is necessary for AGN development. We also show that Regnase-1, an RNA binding protein with endoribonuclease activity, is a negative regulator of IL-17 signaling in RTECs. Accordingly, mice with a selective Regnase-1 deficiency in RTECs exhibited exacerbated kidney dysfunction in AGN. Mechanistically, Regnase-1 inhibits IL-17-driven expression of the transcription factor I{kappa}B{xi} and consequently its downstream gene targets including Il6 and Lcn2. Moreover, deletion of Regnase-1 in human RTECs reduced inflammatory gene expression in an I{kappa}B{xi}-dependent manner. Overall, these data identify an IL-17-driven inflammatory circuit in RTECs during AGN that is constrained by Regnase-1. Monitoring fish communities through environmental DNA metabarcoding in the fish pass system of the second largest hydropower plant in the world The Itaipu Hydroelectric Power Plant is the second largest in the world in power generation. The artificial barrier created by its dam imposes an obstacle for fish migration. Thus, in 2002, a fish pass system, named Piracema Channel, was built to allow fish to access areas upstream of the reservoir. We tested the potential of environmental DNA metabarcoding to monitor the impact of both the dam and associated fish pass system in the Parana River fish communities and to compare it with traditional monitoring methods. Using a fragment of the 12S gene, we characterized richness and community composition based on amplicon sequence variants, operational taxonomic units, and zero-radius OTUs. We combined GenBank and in-house data for taxonomic assignment. We found that different bioinformatics approaches showed similar results. Also, we found a decrease in fish diversity from 2019 to 2020 probably due to the recent extreme drought experienced in southeastern Brazil. The highest alpha diversity was recorded in the mouth of the fish pass system, located in a protected valley with the highest environmental heterogeneity. Despite the clear indication that the reference databases need to be continuously improved, our results demonstrate the analytical efficiency of the metabarcoding to monitor fish species. Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae The endosymbiotic origin of plastids from cyanobacteria gave eukaryotes photosynthetic capabilities and launched the diversification of countless forms of algae. These primary plastids are found in members of the eukaryotic supergroup Archaeplastida. All known archaeplastids still retain some form of primary plastids, which are widely assumed to have a single origin. Here, we used single-cell genomics from natural samples combined with phylogenomics to infer the evolutionary origin of the phylum Picozoa, a globally distributed but seemingly rare group of marine microbial heterotrophic eukaryotes. Strikingly, the analysis of 43 single-cell genomes shows that Picozoa belong to Archaeplastida, specifically related to red algae and the phagotrophic rhodelphids. These picozoan genomes support the hypothesis that Picozoa lack a plastid, and further reveal no evidence of an early cryptic endosymbiosis with cyanobacteria. These findings change our understanding of plastid evolution as they either represent the first complete plastid loss in a free-living taxon, or indicate that red algae and rhodelphids obtained their plastids independently of other archaeplastids. Enhancers regulate 3' end processing activityto control expression of alternative 3'UTR isoforms Multi-UTR genes are widely transcribed but show cell type-specific 3'UTR isoform expression. As transcriptional enhancers regulate mRNA expression, we investigated if they also regulate mRNA isoform expression. Deletion of an endogenous enhancer of a multi-UTR gene did not impair transcript production but prevented a switch in 3'UTR isoform expression. Also, the same enhancers are able to increase transcript production in the context of single-UTR gene promoters, but they increase 3' end processing activity when paired with multi-UTR gene promoters. We show that transcription factors regulate processing activity of weak polyadenylation sites to control cell type-specific alternative 3'UTR isoform expression of widely expressed genes. The glycocalyx affects the mechanotransductive perception of the topographical microenvironment The cell/microenvironment interface is the starting point of integrin-mediated mechanotransduction, but many details of mechanotransductive signal integration remain elusive due to the complexity of the involved (extra)cellular structures, such as the glycocalyx.

We used nano-bio-interfaces reproducing the complex nanotopographical features of the extracellular matrix to analyse the glycocalyx impact on PC12 cell mechanosensing at the nanoscale (e.g., by force spectroscopy with functionalised probes). Our data demonstrates that the glycocalyx configuration affects spatio-temporal nanotopography-sensitive mechanotransductive events at the cell/microenvironment interface. Opposing effects of glycocalyx removal were observed, when comparing flat and specific nanotopographical conditions. The excessive retrograde actin flow speed and force loading are strongly reduced on certain nanotopographies upon removal of the native glycocalyx, while on the flat substrate we observe the opposite trend.

Our results highlight the importance of the glycocalyx configuration in a molecular clutch force loading-dependent cellular mechanism for mechanosensing of microenvironmental nanotopographical features. Poly(rC)-binding protein 1 limits hepatitis C virus assembly and secretion The hepatitis C virus (HCV) co-opts numerous cellular elements - including proteins, lipids, and microRNAs - to complete its viral life cycle. The cellular RNA-binding protein, poly(rC)-binding protein 1 (PCBP1), was previously reported to bind to the 5 untranslated region (UTR) of the HCV genome; however, its importance in the viral life cycle has remained unclear. Herein, we sought to clarify the role of PCBP1 in the HCV life cycle. Using the HCV cell culture (HCVcc) system, we found that knockdown of endogenous PCBP1 resulted in an overall decrease in viral RNA accumulation, yet resulted in an increase in extracellular viral titers. To dissect PCBP1s specific role in the HCV life cycle, we carried out assays for viral entry, translation, genome stability, RNA replication, as well as virion assembly and secretion. We found that PCBP1 knockdown did not directly affect viral entry, translation, RNA stability, or RNA replication, but resulted in an overall increase in infectious particle secretion. This increase in virion secretion was evident even when viral RNA synthesis was inhibited, and blocking virus secretion could partially restore the viral RNA accumulation decreased by PCBP1 knockdown. We therefore propose a model where endogenous PCBP1 normally limits virion assembly and secretion, which increases viral RNA accumulation in infected cells by preventing the departure of viral genomes packaged into virions. Overall, our findings improve our understanding of how cellular RNA-binding proteins influence viral genomic RNA utilization during the HCV life cycle. Reliance of neuronal gene expression on cohesin scales with chromatin loop length Cohesin and CTCF are major drivers of 3D genome organization, but their role in neurons is still emerging. Here we show a prominent role for cohesin in the expression of genes that facilitate neuronal maturation and homeostasis. Unexpectedly, we observed two major classes of activity-regulated genes with distinct reliance on cohesin in primary cortical neurons. Immediate early genes remained fully inducible by KCl and BDNF, and short-range enhancer-promoter contacts at the Immediate early gene Fos formed robustly in the absence of cohesin. In contrast, cohesin was required for full expression of a subset of secondary response genes characterised by long-range chromatin contacts. Cohesin-dependence of constitutive neuronal genes with key functions in synaptic transmission and neurotransmitter signaling also scaled with chromatin loop length. Our data demonstrate that key genes required for the maturation and activation of primary cortical neurons depend cohesin for their full expression, and that the degree to which these genes rely on cohesin scales with the genomic distance traversed by their chromatin contacts. 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. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. BioCAT: a novel tool to search biosynthetic gene clusters producing nonribosomal peptides with a known structure MotivationNonribosomal peptides are a class of secondary metabolites synthesized by multimodular enzymes named nonribosomal peptide synthetases and mainly produced by bacteria and fungi. It has been shown that non-ribosomal peptides have a huge structural and functional diversity including antimicrobial activity, therefore, they are of increasing interest for modern biotechnology. Methods such as NMR and LC-MS/MS allow to determine a peptide structure precisely, but it is often not a trivial task to find natural producers of them. Today, the search is usually performed manually, mostly with tools such as antiSMASH or Prism. However, there are cases when potential producers should be found among hundreds of strains, for instance, when analyzing metagenomes data. Thus, the development of automated approaches is a high-priority task for further nonribosomal peptides research.

ResultsWe developed BioCAT, a two-side approach to find biosynthetic gene clusters which may produce a given nonribosomal peptide when the structure of interesting nonribosomal peptide has already been found. Formally, BioCAT unites the antiSMASH software and the rBAN retrosynthesis tool but some improvements were added to both gene cluster and peptide chemical structure analyses. The main feature of the method is an implementation of position specific score matrix to store specificities of nonribosomal peptide synthetase modules, which has increased the alignment quality in comparison with more strict approaches developed earlier. An ensemble model was implemented to calculate the final alignment score. We tested the method on a manually curated nonribosomal peptides producers database and compared it with a competing tool called GARLIC. Finally, we showed the method applicability on several external examples.

AvailabilityBioCAT is available on the GitHub repository or via pip

Contactkonanovdmitriy@gmail.com Single molecule studies reveal branched pathways for activator-dependent assembly of RNA polymerase II pre-initiation complexes RNA polymerase II (Pol II) transcription reconstituted from purified factors suggests pre-initiation complexes (PICs) can assemble by sequential incorporation of factors at the TATA box. However, these basal transcription reactions are generally independent of activators and co-activators. To study PIC assembly under more realistic conditions, we used single-molecule microscopy to visualize factor dynamics during activator-dependent reactions in nuclear extracts. Surprisingly, Pol II, TFIIF, and TFIIE can pre-assemble on enhancer-bound activators before loading into PICs, and multiple Pol II complexes can bind simultaneously to create a localized cluster. Unlike TFIIF and TFIIE, TFIIH binding is singular and dependent on the basal promoter. Activator-tethered factors exhibit dwell times on the order of seconds. In contrast, PICs can persist on the order of minutes in the absence of nucleotide triphosphates, although TFIIE remains unexpectedly dynamic even after TFIIH incorporation. Our kinetic measurements lead to a new branched model for activator-dependent PIC assembly.

HIGHLIGHTSSingle molecule microscopy reveals unexpected dynamics of RNA Pol II and GTFs

Multiple Pol IIs cluster on UAS/enhancer-bound activators before binding the core promoter

Pol II, TFIIF, and TFIIE, but not TFIIH, can pre-assemble at the UAS/enhancer

Activators increase the rates of Pol II and GTF association with DNA

eTOC BlurbSingle-molecule microscopy experiments by Baek et al. show that RNA polymerase II and basal transcription factors TFIIF and TFIIE preassemble on UAS/enhancer-bound activators, poised for loading into initiation complexes with TFIIH at the core promoter. Transcription activators kinetically enhance factor recruitment, creating a localized cluster of polymerases at the UAS/enhancer. Fiber orientation downsampling compromises the computation of white matter tract-related deformation Incorporating neuroimaging-revealed structural details into finite element (FE) head models opens vast new opportunities to better understand brain injury mechanisms. Recently, growing efforts have been made to integrate fiber orientation from diffusion tensor imaging (DTI) into FE models to predict white matter (WM) tract-related deformation that is biomechanically characterized by tract-related strains. Commonly used approaches often downsample the spatially enriched fiber orientation to match the FE resolution with one orientation per element (i.e., element-wise orientation implementation). However, the validity of such downsampling operation and corresponding influences on the computed tract-related strains remain elusive. To address this, the current study proposed a new approach to integrate voxel-wise fiber orientation from one DTI atlas (isotropic resolution of 1 mm3) into FE models by embedding orientations from multiple voxels within one element (i.e., voxel-wise orientation implementation). By setting the responses revealed by the newly proposed voxel-wise orientation implementation as the reference, we evaluated the reliability of two previous downsampling approaches by examining the downsampled fiber orientation and the computationally predicted tract-related strains secondary to one concussive impact. Two FE models with varying element sizes (i.e., 6.37 {+/-} 1.60 mm and 1.28 {+/-} 0.55 mm, respectively) were incorporated. The results showed that, for the model with a large voxel-mesh resolution mismatch, the downsampled element-wise fiber orientation, with respect to its voxel-wise counterpart, exhibited an absolute deviation over 30{degrees} across the WM/gray matter interface and the pons regions. Accordingly, this orientation deviation compromised the computation of tract-related strains with normalized root-mean-square errors up to 30% and underestimated the peak tract-related strains up to 10%. For the other FE model with finer meshes, the downsampling-induced effects were lower, both on the fiber orientation and tract-related strains. Taken together, the voxel-wise orientation implementation is recommended in future studies as it leverages the DTI-delineated fiber orientation to a larger extent than the element-wise orientation implementation. Thus, this study yields novel insights on integrating neuroimaging-revealed fiber orientation into FE models and may better inform the computation of WM tract-related deformation, which are crucial for advancing the etiological understanding and computational predictability of brain injury.

Graphic abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=142 SRC="FIGDIR/small/471622v2_ufig1.gif" ALT="Figure 1">
View larger version (55K):
org.highwire.dtl.DTLVardef@12ebda9org.highwire.dtl.DTLVardef@1a5d624org.highwire.dtl.DTLVardef@8e068borg.highwire.dtl.DTLVardef@4f8ff5_HPS_FORMAT_FIGEXP  M_FIG C_FIG Prediction of survival after partial hepatectomy using a physiologically based pharmacokinetic model of indocyanine green liver function tests The evaluation of hepatic function and functional capacity of the liver are essential tasks in hepatology as well as in hepatobiliary surgery. Indocyanine green (ICG) is a widely applied test compound that is used in clinical routine to evaluate hepatic function. Important questions for the functional evaluation with ICG in the context of hepatectomy are how liver disease such as cirrhosis alters ICG elimination, and if postoperative survival can be predicted from preoperative ICG measurements. Within this work a physiologically based pharmacokinetic (PBPK) model of ICG was developed and applied to the prediction of the effects of a liver resection under various degrees of cirrhosis. For the parametrization of the computational model and validation of model predictions a database of ICG pharmacokinetic data was established. The model was applied (i) to study the effect of liver cirrhosis and liver resection on ICG pharmacokinetics; and (ii) to evaluate the model-based prediction of postoperative ICG-R15 as a measure for postoperative outcome. Key results are the accurate prediction of changes in ICG pharmacokinetics caused by liver cirrhosis and postoperative changes of ICG-elimination after liver resection, as validated with a wide range of data sets. Based on the PBPK model, individual survival after liver resection could be classified, demonstrating its potential value as a clinical tool. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. Express Arm Responses Appear Bilaterally on Upper-limb Muscles in An Arm Choice Reaching Task When required, humans can generate very short latency reaches towards visual targets, like catching a falling cellphone. During such rapid reaches, express arm responses are the first wave of upper limb muscle recruitment, occurring ~80-100 ms after target appearance. There is accumulating evidence that express arm responses arise from signaling along the tecto-reticulo-spinal tract, but the involvement of the reticulo-spinal tract has not been well-studied. Since the reticulospinal tract projects bilaterally, we studied whether express arm responses would be generated bilaterally. Human participants (n = 14; 7 female) performed visually guided reaches in a modified emerging target paradigm where either arm could intercept the target. We recorded electromyographic activity bilaterally from the pectoralis major muscle. Our analysis focused on target locations where participants reached with the right arm on some trials, and the left arm on others. In support of the involvement of the reticulospinal tract, express arm responses persisted bilaterally regardless of which arm reached to the target. The latency and magnitude of the express arm response did not depend on whether the arm was chosen to reach or not. However, on the reaching arm, the magnitude of the express arm response was correlated to the level of anticipatory activity. The bilateral generation of express arm responses supports the involvement of the reticulo-spinal tract. We surmise that the correlation between anticipatory activity and the magnitude of express arm responses on the reaching arm arises from convergence of cortically-derived signals with a parallel subcortical pathway mediating the express arm response.

New and NoteworthyExpress arm responses have been proposed to arise from the tecto-reticulo-spinal tract originating within the superior colliculus, but the involvement of the reticulo-spinal tract has not been well studied. Here, we show these responses appear bilaterally in a task where either arm can reach to a newly appearing stimulus. Our results suggest that the most rapid visuomotor transformations for reaching are performed by a subcortical pathway. 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. Endometriotic Organoids: A Novel In Vitro Model of Endometriotic Lesion Development. The development and progression of endometriotic lesions are poorly understood, but immune cell dysfunction and inflammation are closely associated with the pathophysiology of endometriosis. A lack of suitable 3D in vitro models permitting the study of interactions between cell types and the microenvironment is a contributing factor. To address this limitation, we developed endometriotic organoids (EO) to explore the role of epithelial-stromal interactions and model peritoneal cell invasion associated with lesion development. Using a non-adherent microwell culture system, spherical organoids were generated with endometriotic epithelial cells (12Z) combined with immortalized endometriotic stromal cells (iEc-ESC) or immortalized uterine stromal cells (iHUF). Organoids self-organized with stromal cells occupying the center and epithelial cells on the periphery of the organoid. Endometriotic organoids (EO), containing iEc-ESC, resulted in the development of stratified 12Z epithelial cells compared to those with iHUF where the 12Z cells developed as a single layered epithelium. Transcriptomic analysis found 4,522 differentially expressed genes (DEG) between EO and 12Z/iHUF organoids, and the top DEG included increased expression of interleukins and prostaglandin synthase enzymes. An overlap of the EO DEG with baboon endometriotic lesions was highly significant. Finally, to mimic invasion of endometrial tissue into the peritoneum, a model was developed using EO and extracellular matrix containing human peritoneal mesothelial cells (LP9). Invasion of EO into the extracellular matrix-LP9 layer was increased in presence of estrogen or THP1-derived proinflammatory macrophages. Taken together, our results strongly support the concept that EO are an appropriate model for dissecting mechanisms that contribute to endometriotic lesion development.

One Sentence SummaryEndometriotic organoids are an appropriate model to study epithelial-stromal interactions and model cell invasion associated with lesion development. The OPAQUE1/DISCORDIA2 myosin XI is required for phragmoplast guidance during asymmetric cell division in maize Asymmetric divisions produce cells with different fates and are critical for development. Here we show that a maize myosin XI protein, OPAQUE1 (O1), is necessary for asymmetric divisions during maize stomatal development. We analyzed stomatal precursor cells prior to and during asymmetric division to determine why o1 mutants have abnormal division planes. Cell polarization and nuclear positioning occur normally in the o1 mutant, and the future site of division is correctly specified. The defect in o1 occurs during late cytokinesis, when the plant-specific phragmoplast -- made of microtubules, actin and other proteins - forms the nascent cell plate. The phragmoplast becomes misguided and does not meet the previously established division site. Initial phragmoplast guidance is correct in o1. However, as phragmoplast expansion continues, phragmoplasts in o1 stomatal precursor cells become misguided and do not meet the cortex at the established division site. To understand how this myosin protein contributes to phragmoplast guidance, we identified O1-interacting proteins. Other myosins, specific actin-binding proteins, and maize kinesins related to the Arabidopsis thaliana division site markers PHRAGMOPLAST ORIENTING KINESINs (POKs) interact with O1. We propose that different myosins are important at multiple steps of phragmoplast expansion, and the O1 actin motor and POK-like microtubule motors work together to ensure correct late-stage phragmoplast guidance. Network Inference from Perturbation Time Course Data Networks underlie much of biology from subcellular to ecological scales. Yet, understanding what experimental data are needed and how to use them for unambiguously identifying the structure of even small networks remains a broad challenge. Here, we integrate a dynamic least squares framework into established modular response analysis (DL-MRA), that specifies sufficient experimental perturbation time course data to robustly infer arbitrary two and three node networks. DL-MRA considers important network properties that current methods often struggle to capture: (i) edge sign and directionality; (ii) cycles with feedback or feedforward loops including self-regulation; (iii) dynamic network behavior; (iv) edges external to the network; and (v) robust performance with experimental noise. We evaluate the performance of and the extent to which the approach applies to cell state transition networks, intracellular signaling networks, and gene regulatory networks. Although signaling networks are often an application of network reconstruction methods, the results suggest that only under quite restricted conditions can they be robustly inferred. For gene regulatory networks, the results suggest that incomplete knockdown is often more informative than full knockout perturbation, which may change experimental strategies for gene regulatory network reconstruction. Overall, the results give a rational basis to experimental data requirements for network reconstruction and can be applied to any such problem where perturbation time course experiments are possible. The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. Honeybee cognition as a tool for scientific engagement In addition to the crowdsourcing of experimental data, citizen science, and scientific engagement more broadly, serve as a bridge between researchers and the wider community. This serves to foster a greater understanding of the scientific method and science-based solutions generally. Apis mellifera (honeybees) are a well-established model for the study of learning and cognition and can serve as an engaging outreach system for this wider community. Here, we developed and implemented a protocol using well established honeybee conditioning protocols to safely study the effects of caffeine and dopamine on learning performance. Using this protocol, a group of high-school aged students as part of the Ryukyu Girls program demonstrated that caffeine, but not dopamine, significantly reduced the number of trials required for a successful conditioning response. This allowed these students to explore the scientific method in a relatable and engaging way.

Simple SummaryGlobal scientific literacy can be improved through widespread and effective community engagement by researchers. We propose Apis mellifera (honeybee) as an public engagement tool due to widespread awareness of colony collapse and the bees importance in food production. Moreover, their cognitive abilities make for engaging experiments. Their relative ease of cultivation means that studies can be performed cost-effectively, especially when partnering with local aperists. Using a proxy for honeybee learning, a group of non-specialist high-school-aged participants obtained data suggesting that caffeine, but not dopamine, improved learning. This hands-on experience facilitated student understanding of the scientific method, factors that shape learning and the importance of learning for hive health.