A general theory for temperature-dependence in biology At present, there is no simple, complete, and first principles-based model for quantitatively describing the full range of observed biological temperature responses. Here, we derive a theory exhibiting these features based on the Eyring-Evans-Polanyi theory governing chemical reaction rates, and which is applicable across all scales from the micro to the macro. Assuming only that the conformational entropy of molecules changes with temperature, we derive a theory for the temperature dependence which takes the form of an exponential function modified by a power-law. Our framework leads to six deductions applicable to any biological trait that depends on temperature, and elucidates novel aspects of universal temperature responses across the tree of life, from quantum to classical scales. All predictions are well supported by data for a wide variety of biological rates and steady states, from molecular to ecological scales and across multiple taxonomic groups. In addition, we provide novel explanations of several empirical relationships including optimal values in temperature response curves.

One-Sentence SummaryWe derive a simple and universal formulae to characterize temperature responses of biological processes across the tree of life.
Joint interactions with humans may pose a higher risk of zoonotic outbreaks than interactions with conspecifics among wildlife populations at human-wildlife interfaces O_LIPandemics caused by wildlife-origin pathogens, like COVID-19, highlight the importance of understanding the ecology of zoonotic transmission and outbreaks among wildlife populations at human-wildlife interfaces. To-date, the relative effects of human-wildlife and wildlife-wildlife interactions on the likelihood of such outbreaks remain unclear.
C_LIO_LIIn this study, we used social network analysis and epidemiological Susceptible Infected Recovered (SIR) models, to track zoonotic outbreaks through wild animals joint propensities to engage in social-ecological co-interactions with humans, and their social grooming interactions with conspecifics.
C_LIO_LIWe collected behavioral and demographic data on 10 groups of macaques (Macaca spp.) living in (peri)urban environments across Asia. Outbreak sizes predicted by the SIR models were related to structural features of the social networks, and particular properties of individual animals connectivity within those networks.
C_LIO_LIOutbreak sizes were larger when the first-infected animal was highly central, in both types of networks. Across host-species, particularly for rhesus and bonnet macaques, the effects of network centrality on outbreak sizes were stronger through macaques human co-interaction networks compared to grooming networks.
C_LIO_LIOur findings, independent of pathogen-transmissibility, suggest that wildlife populations in the Anthropocene are vulnerable to zoonosis more so due to their propensities to aggregate around anthropogenic factors than their gregariousness with conspecifics. Thus, the costs of zoonotic outbreaks may outweigh the potential/perceived benefits of jointly interacting with humans to procure anthropogenic food. From One Health perspectives, animals that consistently interact with both humans and conspecifics across time and space are useful targets for disease spillover assessments and control.
C_LI
Engineered bacteria detect tumor DNA Advances in bacterial engineering have catalysed the development of living cell diagnostics and therapeutics1-3, including microbes that respond to diseases such as gut inflammation4, intestinal bleeding5, pathogens6 and hypoxic tumors7. Bacteria can easily access the entire gastrointestinal tract via oral administration8, and they can produce outputs that can be noninvasively measured in stool4 or urine7. Cellular memory, such as bistable switches4,9,10 or genomic rearrangement11, has been used to allow bacteria to store information over time. However, living biosensors have not yet been engineered to detect specific DNA sequences or mutations from outside the cell. Here, we engineer naturally competent Acinetobacter baylyi to detect donor DNA from the genomes of colorectal cancer (CRC) cells and organoids. We characterize the functionality of the biosensors in vitro with co-culture assays and then validate in vivo with sensor bacteria delivered orally or rectally into mice injected with orthotopic donor CRC organoids. We observe horizontal gene transfer from the tumor to the sensor bacteria in vivo, allowing their detection in stool. The sensor bacteria achieved 100% discrimination between mice with and without CRC using both delivery methods. Our findings establish a framework for biosensing applications that require the detection of mutations or organisms within environments that are difficult to sample. In addition, the platform can be readily expanded to include in situ production and delivery of therapeutic payloads at the detection site. Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments. A NANOG-pERK reciprocal regulatory circuit mediates Nanog autoregulation and ERK signaling dynamics. The self-renewal and differentiation potential of Embryonic stem cells (ESCs) is maintained by the regulated expression of core pluripotency factors. The expression level of core pluripotency factor Nanog is tightly regulated by a negative feedback autorepression loop. However, it remains unclear how the ESCs perceive the NANOG levels and execute autorepression. Here, we show that a dose-dependent induction of Fgfbp1 and Fgfr2 by NANOG activates an autocrine mediated ERK signaling in high-Nanog cells to trigger autorepression. pERK recruits NONO to Nanog locus to repress transcription by preventing POL2 loading. The Nanog autorepression process establishes a self-perpetuating NANOG-pERK reciprocal regulatory circuit. We further demonstrate that the reciprocal regulatory circuit induces the pERK heterogeneity and ERK signaling dynamics in pluripotent stem cells. Cellular feedback to organotelluranes displays modulation of antioxidant proteins gene expression Organotelluranes RT3 and RT4 are thiol reagents that induce mitochondrial transition pore (MTP) opening in a sensitive and insensitive manner to cyclosporin A. Although RT3 and RT4 promote glutathione depletion, paradoxically, they are also an efficient antioxidant for membrane lipids. These compounds antagonistic effects elicited the challenging question of how the gene expression of antioxidant enzymes would respond to treatment with these compounds. The influence of RT3 and RT4 on antioxidant enzyme expression was investigated in cultured aortic smooth muscle cells (ASMC). RT3 and RT4 promoted disruption of ionic calcium homeostasis, mitochondrial transmembrane potential ({Delta}{Psi}), and cell death in a dose-dependent manner. The cell death mechanisms responded qualitatively to the increase of the organotellurane concentration and changed from apoptosis to necrosis. RT3 and RT4 increased the expression of thioredoxin significantly. RT3 also increased the expression of glutaredoxin and glutathione peroxidase, slightly the catalase expression without significant effects on SOD expression. The results are consistent with GSH and protein thiol depletion and discussed based on the cell toxicity mechanism exhibited by these compounds.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=150 SRC="FIGDIR/small/425411v4_ufig1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@52bb82org.highwire.dtl.DTLVardef@19476d9org.highwire.dtl.DTLVardef@6aa99org.highwire.dtl.DTLVardef@1aa605e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Development of an exosomal gene signature to detect residual disease in dogs with osteosarcoma using a novel xenograft platform and machine learning Osteosarcoma has a guarded prognosis. A major hurdle in developing more effective osteosarcoma therapies is the lack of disease-specific biomarkers to predict risk, prognosis, or therapeutic response. Exosomes are secreted extracellular microvesicles emerging as powerful diagnostic tools. However, their clinical application is precluded by challenges in identifying disease-associated cargo from the vastly larger background of normal exosome cargo. We developed a method using canine osteosarcoma in mouse xenografts to distinguish tumor-derived from host-response exosomal mRNAs. The model allows for the identification of canine osteosarcoma-specific gene signatures by RNA sequencing and a species-differentiating bioinformatics pipeline. An osteosarcoma-associated signature consisting of five gene transcripts (SKA2, NEU1, PAF1, PSMG2, and NOB1) was validated in dogs with spontaneous osteosarcoma by qRT-PCR, while a machine learning model assigned dogs into healthy or disease groups. Serum/plasma exosomes were isolated from 53 dogs in distinct clinical groups ("healthy", "osteosarcoma", "other bone tumor", or "non-neoplastic disease"). Pre-treatment samples from osteosarcoma cases were used as the training set and a validation set from post-treatment samples was used for testing, classifying as "osteosarcoma-detected" or "osteosarcoma-NOT detected". Dogs in a validation set whose post-treatment samples were classified as "osteosarcoma-NOT detected" had longer remissions, up to 15 months after treatment. In conclusion, we identified a gene signature predictive of molecular remissions with potential applications in the early detection and minimal residual disease settings. These results provide proof-of-concept for our discovery platform and its utilization in future studies to inform cancer risk, diagnosis, prognosis, and therapeutic response.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=147 SRC="FIGDIR/small/429432v2_ufig1.gif" ALT="Figure 1">
View larger version (41K):
org.highwire.dtl.DTLVardef@1b612c4org.highwire.dtl.DTLVardef@8a1847org.highwire.dtl.DTLVardef@3ffd64org.highwire.dtl.DTLVardef@114df8c_HPS_FORMAT_FIGEXP  M_FIG C_FIG Behavioural analysis of swarming mosquitoes reveals higher hearing sensitivity than previously measured with electrophysiology methods Mosquitoes of many species mate in station-keeping swarms. Mating chases ensue as soon as a male detects the flight tones of a female with his auditory organs. Previous studies of hearing thresholds have mainly used electrophysiological methods that prevent the mosquito from flying naturally. The main aim of this study was to quantify behaviourally the sound-level threshold at which males can hear females. Free-flying male Anopheles coluzzii were released in a large arena (~2 m high x 2 m x 1 m) with a conspicuous object on the ground that stimulates swarming behaviour. Males were exposed to a range of natural and synthetic played-back sounds of female flight. We monitored the responses of males and their distance to the speaker by recording changes in their wingbeat frequency and angular speed. We show that the mean male behavioural threshold of particle-velocity hearing lies between 13-20 dB SVL (95%-CI). A conservative estimate of 20 dB SVL (i.e., < 0.5 m/s particle velocity) is already 12 to 26 dB lower than most of the published electrophysiological measurements from the Johnstons organ. In addition, we suggest that 1) the first harmonic of female flight-sound is sufficient for males to detect her presence, 2) males respond with a greater amplitude to single-female sounds than to the sound of a group of females and 3) the response of males to the playback of the flight sound of a live female is the same as that of a recorded sound of constant frequency and amplitude. Identifying and testing marker-trait associations for growth and phenology in three pine species: implications for genomic prediction In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 SNPs from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. Genetic Loci and Metabolic States Associated With Murine Epigenetic Aging Changes in DNA methylation (DNAm) are linked to aging. Here, we profile highly conserved CpGs in 339 predominantly female mice belonging to the BXD family for which we have deep longevity and genomic data. We use a  pan-mammalian microarray that provides a common platform for assaying the methylome across mammalian clades. We computed epigenetic clocks and tested associations with DNAm entropy, diet, weight, metabolic traits, and genetic variation. We describe the multifactorial variance of methylation at these CpGs, and show that high fat diet augments the age-associated changes. Entropy increases with age. The progression to disorder, particularly at CpGs that gain methylation over time, was predictive of genotype-dependent life expectancy. The longer-lived BXD strains had comparatively lower entropy at a given age. We identified two genetic loci that modulate rates of epigenetic age acceleration (EAA): one on chromosome (Chr) 11 that encompasses the Erbb2/Her2 oncogenic region, and a second on Chr19 that contains a cytochrome P450 cluster. Both loci harbor genes associated with EAA in humans including STXBP4, NKX2-3, and CUTC. Transcriptome and proteome analyses revealed associations with oxidation-reduction, metabolic, and immune response pathways. Our results highlight concordant loci for EAA in humans and mice, and demonstrate a tight coupling between the metabolic state and epigenetic aging. A bacterial extracellular vesicle-based intranasal vaccine against SARS-CoV-2 protects against disease and elicits neutralizing antibodies to wild-type virus and Delta variant Several vaccines have been introduced to combat the coronavirus infectious disease-2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Current SARS-CoV-2 vaccines include mRNA-containing lipid nanoparticles or adenoviral vectors that encode the SARS-CoV-2 Spike (S) protein of SARS-CoV-2, inactivated virus, or protein subunits. Despite growing success in worldwide vaccination efforts, additional capabilities may be needed in the future to address issues such as stability and storage requirements, need for vaccine boosters, desirability of different routes of administration, and emergence of SARS-CoV-2 variants such as the Delta variant. Here, we present a novel, well-characterized SARS-CoV-2 vaccine candidate based on extracellular vesicles (EVs) of Salmonella typhimurium that are decorated with the mammalian cell culture-derived Spike receptor-binding domain (RBD). RBD-conjugated outer membrane vesicles (RBD-OMVs) were used to immunize the golden Syrian hamster (Mesocricetus auratus) model of COVID-19. Intranasal immunization resulted in high titers of blood anti-RBD IgG as well as detectable mucosal responses. Neutralizing antibody activity against wild-type and Delta variants was evident in all vaccinated subjects. Upon challenge with live virus, hamsters immunized with RBD-OMV, but not animals immunized with unconjugated OMVs or a vehicle control, avoided body mass loss, had lower virus titers in bronchoalveolar lavage fluid, and experienced less severe lung pathology. Our results emphasize the value and versatility of OMV-based vaccine approaches. Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems. Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536. RNA degradation of blood reveals a time diagnostic window in a controlled environment Determining the time since deposition (TSD) of bloodstains would provide forensic scientists with critical information regarding the timeline of the events involving bloodshed. The physicochemical changes occurring to biomolecules as a bloodstain degrades can be used to approximate the TSD of bloodstains. Our study aims to quantify the timewise degradation trends and temperature dependence found in total RNA from bloodstains without the use of amplification, expanding the scope of the RNA TSD research which has previously targeted mRNA molecules. Whole bovine blood was stored in plastic microcentrifuge tubes at 21{degrees}C or 4{degrees}C and tested over different timepoints spanning one week. Nine RNA metrics were visually assessed and quantified using linear and mixed models; the RNA Integrity Number equivalent (RINe) and the DV200 demonstrated strong negative trends over time and statistical independence. The RINe model fit was high (R2 = 0.60), and while including the biological replicate as a random effect increased the fit for all RNA metrics, no significant differences were found between biological replicates stored at the same temperature for the RINe and DV200 metrics. Importantly, this suggests that these standardized metrics can likely be directly compared between scenarios and individuals, with DV200 having an inflection point at [~]28 hrs. This study provides a novel approach for blood TSD estimates, producing metrics that are not affected by inter-individual variation and improving our understanding of the rapid degradation occurring in bloodstains.

HIGHLIGHTSO_LIAmplification-free analysis of total RNA in degrading bloodstains.
C_LIO_LIShort-term RNA degradation assessment using high-resolution size measurements.
C_LIO_LITotal RNA quality and quantity metrics were assessed across a one-week.
C_LIO_LITotal RNA quality metrics demonstrated the strongest timewise trends.
C_LIO_LIBiological replicates produced similar results for RNA quality metrics.
C_LI Inhibition of microbial deconjugation of micellar bile acids protects against intestinal permeability and liver injury Altered host-microbe interactions and increased intestinal permeability have been implicated in the pathogenesis of a range of diseases. However, the mechanisms by which gut microbes affect epithelial barrier integrity remain unclear. Few host-produced metabolites that protect against epithelial damage have been identified, and whether microbial metabolism of host factors alters intestinal barrier function is unknown. Here, we investigate the effects of bacterial metabolism of host-produced bile acid (BA) metabolites on epithelial barrier integrity. We observe that rats fed a choline-deficient, high-fat diet (CDAHFD) exhibit reduced abundance of host-produced conjugated BAs in the intestine at early timepoints coinciding with increased permeability. We show that in vitro, conjugated BAs protect gut epithelial monolayers from damage caused by bacterially produced unconjugated BAs through micelle formation. We then demonstrate that inhibition of BA deconjugation using a small molecule inhibitor of gut bacterial bile salt hydrolase (BSH) enzymes prevents development of pathologic intestinal permeability and hepatic inflammation in CDAHFD-fed rats. Finally, we show that the predominant conjugated BAs in humans protect against epithelial barrier disruption in vitro. Our study identifies a protective role for conjugated BAs in intestinal epithelial barrier function and suggests that rational manipulation of microbial BA metabolism could be leveraged to regulate gut barrier integrity. TAP: Targeting and analysis pipeline for optimization and verification of coil placement in transcranial magnetic stimulation ObjectiveTranscranial magnetic stimulation (TMS) can modulate brain function via an electric field (E-field) induced in a region of interest (ROI) in the brain. The ROI E-field can be computationally maximized and set to match a specific reference using individualized head models to find the optimal coil placement and stimulus intensity. However, the available software lacks many practical features for prospective planning of TMS interventions and retrospective evaluation of the experimental targeting accuracy.

ApproachThe TMS targeting and analysis pipeline (TAP) software uses an MRI/fMRI-derived brain target to optimize coil placement considering experimental parameters such as the subjects hair thickness and coil placement restrictions. The coil placement optimization is implemented in SimNIBS 3.2, for which an additional graphical user interface (TargetingNavigator) is provided to visualize and adjust procedural parameters. The coil optimization process also computes the E-field at the target, allowing the selection of the TMS device intensity setting to achieve a specific E-field strength. The optimized coil placement information is prepared for neuronavigation software, which supports targeting during the TMS procedure. The neuronavigation system can record the coil placement during the experiment, and these data can be processed in TAP to quantify the accuracy of the experimental TMS coil placement and induced E-field.

Main resultsTAP was demonstrated in a study consisting of three repetitive TMS sessions in five subjects. TMS was delivered by an experienced operator under neuronavigation with the computationally optimized coil placement. Analysis of the experimental accuracy from the recorded neuronavigation data indicated coil location and orientation deviations up to about 2 mm and 2{degrees}, respectively, resulting in an 8% median decrease in the target E-field magnitude compared to the optimal placement.

SignificanceTAP supports navigated TMS with a variety of features for rigorous and reproducible stimulation delivery, including planning and evaluation of coil placement and intensity selection for E-field-based dosing. 3D Reconstructions of Mouse Skeletal Muscle Reveal a Decrease in the MICOS Complex and Altered Mitochondrial Networks Background: Skeletal muscle gradually loses mass, strength, endurance, and oxidative capacity during aging. Studies of bioenergetics and protein turnover show that mitochondria mediate this decline in function. Mitochondria cristae are essential for the production of ATP. While mitochondrial aging is associated with endoplasmic reticulum stress, fragmented mitochondria, and decreased mitochondrial capacity, the genes associated with morphological changes in mitochondria during aging still requires further elucidation Further, it is not completely understood how 3D mitochondrial networks and the specialization of mitochondria are altered during aging. Methods: We measured changes in mitochondrial morphology and mitochondrial connectivity during the aging of the mouse gastrocnemius muscle through serial block facing-scanning electron microscopy and 3D reconstruction. Nanotunnels were also measured through 3D reconstruction. CRISPR/Cas9 KO was performed to examine changes in mitochondria upon loss of the mitochondrial contact site and cristae organizing system (MICOS) complex and optic atrophy 1 (OPA-1). Metabolomics were used to find key metabolite pathways changed upon off of the MICOS complex. Results: We found changes in mitochondrial network configuration, nanotunneling, size, shape, number, contact sites, cristae organizing system (MICOS) dynamics and gene expression in skeletal muscle during aging. We also found an association of OPA-1 and the MICOS complex in the gastrocnemius with mitochondrial aging. Further, the loss of the MICOS complex was linked with decreased oxidative capacity, and altered mitochondrial metabolism. Conclusions: To our knowledge, studies of mitochondrial changes in skeletal muscle throughout aging remains limited. 3D reconstructions of nanotunnels elucidated novel patterns of mitochondrial aging in skeletal muscle. MICOS proteins decreased with age and mitochondrial morphology was similar between aged skeletal muscle and that of young mice with MICOS protein loss. In tandem, our data suggest a relationship between the MICOS complex and aging, which could be potentially linked to disease states with additional 3D reconstruction. SARS-CoV-2 spike-glycoprotein processing at S1/S2 and S2' and shedding of the ACE2 viral receptor: roles of Furin and TMPRSS2 and implications for viral infectivity and cell-to-cell fusion The spike (S)-protein of SARS-CoV-2 binds ACE2 and requires proteolytic "priming" at PRRAR685{downarrow} into S1 and S2 (cleavage at S1/S2), and "fusion-activation" at a S2 site for viral entry. In vitro, Furin cleaved peptides mimicking the S1/S2 cleavage site more efficiently than at the putative S2, whereas TMPRSS2 inefficiently cleaved both sites. In HeLa cells Furin-like enzymes mainly cleaved at S1/S2 during intracellular protein trafficking, and S2 processing by Furin at KPSKR815{downarrow} was strongly enhanced by ACE2, but not for the optimized S2 KRRKR815{downarrow} mutant (S2), whereas individual/double KR815AA mutants were retained in the endoplasmic reticulum. Pharmacological Furin-inhibitors (Boston Pharmaceuticals, BOS-inhibitors) effectively blocked endogenous S-protein processing in HeLa cells. Furthermore, we show using pseudotyped viruses that while entry by a "pH-dependent" endocytosis pathway in HEK293 cells did not require Furin processing at S1/S2, a "pH-independent" viral entry in lung-derived Calu-3 cells was sensitive to inhibitors of Furin (BOS) and TMPRSS2 (Camostat). Consistently, these inhibitors potently reduce infectious viral titer and cytopathic effects, an outcome enhanced when both compounds were combined. Quantitative analyses of cell-to-cell fusion and spike processing revealed the key importance of the Furin sites for syncytia formation. Our assays showed that TMPRSS2 enhances fusion and proteolysis at S2 in the absence of cleavage at S1/S2, an effect that is linked to ACE2 shedding by TMPRSS2. Overall, our results indicate that Furin and TMPRSS2 play synergistic roles in generating fusion-competent S-protein, and in promoting viral entry, supporting the combination of Furin and TMPRSS2 inhibitors as potent antivirals against SARS-CoV-2.

IMPORTANCESARS-CoV-2 is the etiological agent of COVID-19 that resulted in >5 million deaths. The spike protein (S) of the virus directs infection of the lungs and other tissues by binding the angiotensin-converting enzyme 2 (ACE2) receptor. For effective infection, the S-protein is cleaved at two sites: S1/S2 and S2. Cleavage at S1/S2, induces a conformational change favoring the recognition of ACE2. The S2 cleavage is critical for cell-to-cell fusion and virus entry into host cells. Our study contributes to a better understanding of the dynamics of interaction between Furin and TMPRSS2 during SARS-CoV-2 entry and suggests that the combination of a non-toxic Furin inhibitor with a TMPRSS2 inhibitor could significantly reduce viral entry in lung cells, as evidenced by an average synergistic [~]95% reduction of viral infection. This represents a powerful novel antiviral approach to reduce viral spread in individuals infected by SARS-CoV-2 or future related coronaviruses. N-VEGF, the autoregulatory arm of VEGF-A Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. Adversarial domain translation networks for fast and accurate integration of large-scale atlas-level single-cell datasets The rapid emergence of large-scale atlas-level single-cell RNA-seq datasets presents remarkable opportunities for broad and deep biological investigations through integrative analyses. However, harmonizing such datasets requires integration approaches to be not only computationally scalable, but also capable of preserving a wide range of fine-grained cell populations. We created Portal, a unified framework of adversarial domain translation to learn harmonized representations of datasets. With innovation in model and algorithm designs, Portal achieves superior performance in preserving biological variation during integration, while achieving integration of millions of cells in minutes with low memory consumption. We show that Portal is widely applicable to integrating datasets across samples, platforms and data types (including scRNA-seq, snRNA-seq and scATAC-seq). Finally, we demonstrate the power of Portal by applying it to the integration of cross-species datasets with limited shared information among them, elucidating biological insights into the similarities and divergences in the spermatogenesis process among mouse, macaque and human. GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. A time irreversible model of nucleotide substitution for the coronavirus SARS-CoV-2 SARS-CoV-2 is the cause of the worldwide epidemic of severe acute respiratory syndrome. Evolutionary studies of the virus genome will provide a predictor of the fate of COVID-19 in the near future. Recent studies of the virus genomes have shown that C to U substitutions are overrepresented in the genome sequences of SARS-CoV-2. Traditional time-reversible substitution models cannot be applied to the evolution of SARS-CoV-2 sequences. Therefore, in this study, I propose a new time-irreversible model and a new method for estimating the nucleotide substitution rate of SARS-CoV-2. Computer simulations showed that that the new method gives good estimates. I applied the new method to estimate nucleotide substitution rates of SARS-CoV-2 sequences. The result suggests that the rate of C to U substitution of SARS-Cov-2 is ten times higher than other types of substitutions. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. 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. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions.