Systematic temporal analysis of peripheral blood transcriptomes using TrendCatcher identifies early and persistent neutrophil activation as a hallmark of severe COVID-19 Studying temporal gene expression shifts during disease progression provides important insights into the biological mechanisms that distinguish adaptive and maladaptive responses. Existing tools for the analysis of time course transcriptomic data are not designed to optimally identify distinct temporal patterns when analyzing dynamic differentially expressed genes (DDEGs). Moreover, there is a lack of methods to assess and visualize the temporal progression of biological pathways mapped from time course transcriptomic datasets. In this study, we developed an open-source R package TrendCatcher (https://github.com/jaleesr/TrendCatcher), which applies the smoothing spline ANOVA model and break point searching strategy to identify and visualize distinct dynamic transcriptional gene signatures and biological processes from longitudinal datasets. We used TrendCatcher to perform a systematic temporal analysis of COVID-19 peripheral blood transcriptomes, including bulk RNA-seq and scRNA-seq time course data. TrendCatcher uncovered the early and persistent activation of neutrophils and coagulation pathways as well as impaired type I interferon (IFN-I) signaling in circulating cells as a hallmark of patients who progressed to severe COVID-19, whereas no such patterns were identified in individuals receiving SARS- CoV-2 vaccinations or patients with mild COVID-19. These results underscore the importance of systematic temporal analysis to identify early biomarkers and possible pathogenic therapeutic targets.
Alpha-Frag: a deep neural network for fragment presence prediction improves peptide identification Identification of peptides in mass spectrometry-based proteomics typically relies on spectra matches. As MS/MS spectra record presence and intensity of fragment ions, the match should take both fragment presence similarity and intensity similarity into consideration. Fragment presence similarity can be calculated with the help of fragment presence prediction such as theoretical enumeration of all possible fragment ions or selecting non-zero intensity ions from the result of fragment intensity prediction, but neither of these two methods is accurate enough. In this work, we developed a deep neural network based model, Alpha-Frag, to predict precisely the fragment ions that should be present for a given peptide. Alpha-Frag modelled fragment presence prediction as a multi-label classification task and trained with ProteomeTools dataset. In terms of intersection over union (IoU), Alpha-Frag achieved an average of >0.7 and outperformed the benchmarks across the validation datasets. Furthermore, fragment presence similarity was calculated based on presence prediction and incorporated into the peptide statistical validation tools as an additional score to improve peptide identifications. Our preliminary experiments show that this score led to a maximum increase of 26.8% (FDR 0.1%) and 21.6% (FDR 1%) for the DDA and the DIA identification, respectively.

Significance StatementA better prediction of fragmentation for peptides in mass spectrometry (MS) is beneficial to the peptide identification. As the MS/MS spectra record two-dimensional information of fragment ions derived from precursors, mass-to-charge ratio (m/z) and their corresponding intensities, besides the fragment intensity prediction, it is necessary to study the presence prediction. Although the presence prediction can be realized by enumerating all the possible fragmentation patterns of a peptide with equal probability or by selecting non-zero intensity fragment ions from the result of fragment intensity prediction, neither of these two methods is accurate enough. In this study, deep learning is leveraged to precisely predict the fragment ions of a given peptide. Based on the fragment presence prediction, fragment presence similarity between experimental spectra and predicted spectra can be calculated which is proved to promote the peptide detections both for DDA and for DIA data.
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. Analysis of caspase-mediated suppression of the cGAS/STING pathway in Kaposi's sarcoma associated herpesvirus lytic infection reveals a dramatic cellular heterogeneity in type I interferon responses As a result of the ongoing virus-host arms race, viruses have evolved numerous immune subversion strategies, many of which are aimed at suppressing the production of type I interferons (IFNs)1. This focus on IFN evasion highlights the essential role of these cytokines in controlling viral infections. Apoptotic caspases have recently emerged as important regulators of type I IFN signaling in both non-infectious contexts and during viral infection2-5. Despite being widely considered anti-viral factors since they can trigger cell death, several apoptotic caspases promote viral replication by suppressing innate immune responses2-5. Indeed, we previously discovered that the AIDS-associated oncogenic gammaherpesvirus Kaposis sarcoma-associated herpesvirus (KSHV) exploits caspase-8 activity to suppress the antiviral type I IFN response and promote viral replication5. However, the mechanism of this novel viral immune evasion strategy is poorly understood, particularly how caspase-8 antagonizes IFN signaling during KSHV infection. Here we show that caspase activity inhibits the DNA sensor cGAS6 during KSHV lytic replication to block type I IFN induction. Furthermore, we use single-cell RNA-sequencing to reveal that the potent antiviral state conferred upon caspase inhibition is mediated by an exceptionally small percentage of IFN-{beta}-producing cells, thus uncovering further complexity of IFN regulation during viral infection. Collectively, these results provide insight into multiple levels of cellular type I IFN regulation that viruses co-opt for immune evasion. Unraveling these mechanisms can inform targeted therapeutic strategies for viral infections and reveal cellular mechanisms of regulating interferon signaling in the context of cancer and chronic inflammatory diseases. Decoding the porcine developing spatial processing system and production of human entorhinal stellate cell-like cells by a direct programming approach. Classic studies investigating how and when the entorhinal cortex (component of the memory processing system of the brain) develops have been based on traditional thymidine autoradiography and histological techniques. In this study, we take advantage of modern technologies to trace at a high resolution, the cellular complexity of the developing porcine medial entorhinal cortex by using single-cell profiling. The postnatal medial entorhinal cortex comprises 4 interneuron, 3 pyramidal neuron and 2 stellate cell populations which emerge from intermediate progenitor and immature neuron populations. We discover four MGE-derived interneurons and one CGE-derived interneuron population as well as several IN progenitors. We also identify two oligodendrocyte progenitor populations and three populations of oligodendrocytes. We perform a proof-of-concept experiment demonstrating that porcine scRNA-seq data can be used to develop novel protocols for producing human entorhinal cells in-vitro. We identified six transcription factors (RUNX1A1, SOX5, FOXP1, MEF2C, TCF4, EYA2) important in neurodevelopment and differentiation from one RELN+ stellate cell population. Using a lentiviral vector approach, we reprogrammed human induced pluripotent stem cells into stellate cell-like cells which expressed RELN, SATB2, LEF1 and BCL11B. Our findings contribute to the understanding of the formation of the brains cognitive memory and spatial processing system and provides proof-of-concept for the production of entorhinal cells from human pluripotent stem cells in-vitro. 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. Genome-scale metabolic modeling reveals SARS-CoV-2-induced metabolic changes and antiviral targets Tremendous progress has been made to control the COVID-19 pandemic caused by the SARS-CoV-2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS-CoV-2 infection using genome-scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS-CoV-2 infection. We next applied the GEM-based metabolic transformation algorithm to predict anti-SARS-CoV-2 targets that counteract the virus-induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco-2 cells. Further generating and analyzing RNA-sequencing data of remdesivir-treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti-SARS-CoV-2 drug. Our study provides clinical data-supported candidate anti-SARS-CoV-2 targets for future evaluation, demonstrating host metabolism-targeting as a promising antiviral strategy. Genome-wide analysis of the apple family 1 glycosyltransferases identified a flavonoid-modifying UGT, MdUGT83L3, which is targeted by MdMYB88 and contributes to stress adaptation Nowadays, the plant family 1 UDP-glycosyltransferases (UGTs) are more and more investigated owing to their contribution to the plant secondary metabolism and diverse biological roles. Apple (Malus domestica) is one of the most widely cultivated fruit trees with great economical importance. However, poor knowledge is known about the apple UGTs. In this study, we identified 229 members of family 1 through a genome-wide analysis of the apple UGTs, which were clustered into 18 groups, from A to R. We also performed detailed analysis to 34 apple UGTs with quantitavive RT-PCR, and discovered a number of stress-regulated UGTs. Among them, we characterized the role of MD09G1064900, also named as MdUGT83L3, which is significantly induced by salt and cold. In vivo analysis showed that it has high activity towards cyanidin, and moderate activity towards quercetin and keampferol. The transgenic callus and regenerated apple plants overexpressing MdUGT83L3 showed enhanced tolerance to salt and cold treatments. Overexpression of MdUGT83L3 also increased anthocyanin accumulation in the callus tissues and enhanced ROS clearing upon exposed to salt and cold stresses. Futhermore, via yeast-two-hybrid assay, EMSA and CHIP analysis, we also found that MdUGT83L3 could be directly regualted by MdMYB88. Our study indicated that MdUGT83L3, under the regulation of MdMYB88, plays important roles in salt and cold stress adaptation via modulating flavonoid metabolism in apple. Modeling the SARS-CoV-2 nsp1--5'-UTR complex via extended ensemble simulations Nonstructural protein 1 (nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a 180-residue protein that blocks translation of host mRNAs in SARS-CoV-2-infected cells. Although it is known that SARS-CoV-2s own RNA evades nsp1s host translation shutoff, the molecular mechanism underlying the evasion was poorly understood. We performed an extended ensemble molecular dynamics simulation to investigate the mechanism of the viral RNA evasion. Simulation results showed that the stem loop structure of the SARS-CoV-2 RNA 5-untranslated region (SL1) is recognized by both nsp1s globular region and intrinsically disordered region. The recognition presumably enables selective translation of viral RNAs. Cluster analysis of the binding mode and detailed analysis of the binding poses revealed several residues involved in the SL1 recognition mechanism. The simulation results imply that the nsp1 C-terminal helices are lifted from the 40S ribosome upon the binding of SL1 to nsp1, unblocking translation of the viral RNA. Barbary ground squirrels do not have a sentinel system but instead synchronise vigilance Coordinated behavior, such as hunting in lions and coordinated vigilance as antipredator behavior, are examples of benefits of group-living. Instead of asynchronous vigilance, some social species synchronize their vigilance bouts or take turns acting as sentinels. To increase our knowledge on the evolution of vigilance behavior, we studied whether vigilance is coordinated in Barbary ground squirrels, Atlantoxerus getulus. We show that vigilance was synchronized instead of taking turns. Multiple non-mutually exclusive hypotheses could explain synchronization: Barbary ground squirrels may perch because 1) neighbors are perched (copying effect), 2) perch synchrony may be an emergent property of the ecology as all squirrels may be satiated at the same time (collective behavior), or 3) the benefits are large in terms of evading ambush predators and scanning effectiveness (watch each others back). Particularly, in habitats where the field of view is obstructed by man-made structures and multiple individuals may be necessary to watch for terrestrial predators, synchronized vigilance may have greater fitness benefits than sentinel behavior. We conclude that it is essential to test assumptions of coordination and, thus, to analyze coordination to describe sentinel systems.

Significance StatementVigilance behavior can be vital to an animals survival. Taking turns acting as sentinels or synchronizing vigilance bouts reduces the cost of the trade-off between feeding and predation risk. A sentinel system assumes that sentinels are vigilant from raised positions, warn group members of danger, and alternate vigilance bouts. However, the assumption of alternating vigilance bouts remains poorly tested. We tested this assumption in invasive Barbary ground squirrels. We found that instead of alternating, individuals synchronized their vigilance bouts. Perch synchrony may be 1) a response to perching group members (copying effect), 2) an emergent property of the species ecology, and 3) an adaptation to anthropogenically altered habitats (watch each others back). Gymnasts' ability to modulate bioelectric sensorimotor rhythms during kinesthetic motor imagery of sports non-specific movements Previous psychological studies using questionnaires have consistently reported that athletes have superior motor imagery ability, both for sports-specific and sports non-specific movements. However, regarding motor imagery of sports non-specific movements, no physiological studies have demonstrated differences in neural activity between athletes and non-athletes. The purpose of the present study was to examine differences in bioelectric sensorimotor rhythms during kinesthetic motor imagery (KMI) of sports non-specific movements between gymnasts and non-gymnasts. We selected gymnasts as an example population because they are likely to have particularly superior motor imagery ability due to frequent usage of motor imagery including KMI as part of daily practice. Healthy young participants (16 gymnasts and 16 non-gymnasts) performed repeated motor execution and KMI of sports non-specific movements (wrist dorsiflexion and shoulder abduction of the dominant hand). Scalp electroencephalogram (EEG) was recorded over the contralateral sensorimotor cortex. During motor execution and KMI, sensorimotor EEG power is known to decrease in the - (8-15 Hz) and {beta}-bands (16-35 Hz), referred to as event-related desynchronization (ERD). We calculated the maximal peak of ERD both in the - (ERDmax) and {beta}-bands ({beta}ERDmax) as a measure of changes in corticospinal excitability. ERDmax was significantly greater in gymnasts, who subjectively evaluated their KMI as being more vivid, for both KMI tasks. On the other hand, {beta}ERDmax was greater in gymnasts only for shoulder abduction KMI. These findings suggest gymnasts signature of flexibly modulating sensorimotor rhythm with no movements, which may be the basis of their superior ability of KMI for sports non-specific movements.

New & NoteworthyKinesthetic motor imagery of sports non-specific movements was compared between gymnasts and non-gymnasts (i.e., healthy controls) from both physiological and psychological approaches. The EEG sensorimotor rhythms during kinesthetic motor imagery were more desynchronized in gymnasts who subjectively imaged their own movements as being more vivid. The work reveals novel ability in gymnasts to flexibly control their sensorimotor rhythms with no actual movements, which would be the basis of their superior ability of motor imagery. Deep learning-based quantification of arbuscular mycorrhizal fungi in plant roots Soil fungi establish mutualistic interactions with the roots of most vascular land plants. Arbuscular mycorrhizal (AM) fungi are among the most extensively characterised mycobionts to date. Current approaches to quantifying the extent of root colonisation and the abundance of hyphal structures in mutant roots rely on staining and human scoring involving simple, yet repetitive tasks prone to variations between experimenters.

We developed AMFinder which allows for automatic computer vision-based identification and quantification of AM fungal colonisation and intraradical hyphal structures on ink-stained root images using convolutional neural networks.

AMFinder delivered high-confidence predictions on image datasets of roots of multiple plant hosts (Nicotiana benthamiana, Medicago truncatula, Lotus japonicus, Oryza sativa) and captured the altered colonisation in ram1-1, str and smax1 mutants. A streamlined protocol for sample preparation and imaging allowed us to quantify mycobionts from the genera Rhizophagus, Claroideoglomus, Rhizoglomus and Funneliformis via flatbed scanning or digital microscopy including dynamic increases in colonisation in whole root systems over time.

AMFinder adapts to a wide array of experimental conditions. It enables accurate, reproducible analyses of plant root systems and will support better documentation of AM fungal colonisation analyses. AMFinder can be accessed here: https://github.com/SchornacklabSLCU/amfinder.git Stimulated oxygen transport in tissue by magnetic needle acupuncture AimTo understand physiology of magnetic needle acupuncture by investigating O2 transport in tissue during needle intervention.

MethodsO2 transport in tissue is modeled by utilizing COMSOL with magnetic needle inserted into muscle tissue in a 2D porous media. The damaged tissue has been mimicked by an extracted tissue block with 1st order O2 consumption rate. The convection-diffusion O2 transport in the damaged tissue has been further evaluated by varying magnetic flux density B0 of the needle (0-1 T), myoglobin concentration (0-1 mM), O2 tension (5-100 Torr), O2 consumption rates, tissue permeability (10-12-10-6 m2) and porosity (0.1-0.9).

Results1) Active O2 transport carried by interstitial flow is enhanced with the intervention of a magnetic needle by generating a high gradient magnetic field around the tip, which exerts a strong force (104 N/m3) on the diamagnetic interstitial fluid to accelerate the flow. 2) This interstitial flow can reach 30 m/s at B0 = 1T and strongly correlates to the needle B0 and tissue permeability. 3) The needle stirs the interstitial flow can pump O2 flux by 1-2 orders of magnitude compared to that without magnetic field. 4) The enhancement of active O2 transport by magnetic needle is site-specific to the tissue in the vicinity of the tip. This enhancement is more effective in edema condition with a high tissue permeability (>10-9 m2).

ConclusionsThe dramatic enhancement of O2 transport to restore the O2 mitochondria metabolism for dysfunctional muscle tissues is the fundamental physiological mechanism of magnetic needle acupuncture. Bioinformatic and mechanistic analysis of the palmerolide PKS-NRPS biosynthetic pathway from the microbiome of an Antarctic ascidian Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the [~]75 kb trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicohabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression. 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. Transmembrane helices are an overlooked and evolutionarily conserved source of major histocompatibility complex class I and II epitopes Cytolytic T cell responses are predicted to be biased towards membrane proteins. The peptide-binding grooves of most haplotypes of histocompatibility complex class I (MHC-I) are relatively hydrophobic, therefore peptide fragments derived from human transmembrane helices (TMHs) are predicted to be presented more often as would be expected based on their abundance in the proteome. However, the physiological reason of why membrane proteins might be over-presented is unclear. In this study, we show that the over-presentation of TMH-derived peptides is general, as it is predicted for bacteria and viruses and for both MHCI and MHC-II. Moreover, we show that TMHs are evolutionarily more conserved, because single nucleotide polymorphisms (SNPs) are present relatively less frequently in TMH-coding chromosomal regions compared to regions coding for extracellular and cytoplasmic protein regions. Thus, our findings suggest that both cytolytic and helper T cells respond more to membrane proteins, because these are evolutionary more conserved. We speculate that TMHs therefore are less prone to escape mutations that enable pathogens to evade T cell responses. Beneficial effects of novel strains of Aureobasidium pullulans produced 1,3-1,6 β-glucans on non-esterified fatty acid levels in diabetic KKAy mice Obesity, metabolic syndrome, associated lipotoxicity and its cascade of events contribute to the majority of the burden related to non-communicable diseases globally. Preventive lifestyle changes aside, several beneficial effects have been reported in type II diabetes mellitus and dyslipidaemia patients with biological response modifier glucans (BRMG) produced as an exopolysaccharide by Aureobasidium pullulans. In this study, we compared two strains (AFO-202 and N-163) that produce beta glucans in alleviating lipotoxicity. This study was performed in obese diabetic mice model of KK-Ay mice, in four groups with six subjects in each group - Group 1: sacrificed on Day 0 for baseline values; Group 2: control (drinking water); Group 3: AFO-202 beta glucan--200 mg/kg/day; Group 4: N-163 beta glucan--300 mg/kg/day. The animals in groups 2-4 had the test solutions administered by gavage once daily for 28 consecutive days. Biochemical analyses were conducted of blood glucose, triglycerides, total cholesterol, LDL cholesterol, HDL cholesterol and non-esterified fatty acids (NEFA). Group 4 (N-163) had the lowest NEFA levels, as compared to the other groups, and marginally decreased triglyceride levels. The groups had no significant differences in blood glucose, HbA1c, triglycerides, or LDL and HDL cholesterol. N-163 produced by A. pullulans decreased NEFA in a diabetic mice model in 28 days. These results, although modest, warrant further in-depth research into lipotoxicity and associated inflammatory cascades in both healthy and disease affected subjects to develop novel strategies for prevention and management. Immunomodulatory Contribution of Mast Cells to the Regenerative Biomaterial Microenvironment Bioactive immunomodulatory biomaterials have shown promise for influencing the immune response to promote tissue repair and regeneration. Macrophages and T cells have been associated with this response; however, other immune cell types have been traditionally overlooked. In this study, we investigated the role of mast cells in the regulation of the immune response to decellularized biomaterial scaffolds using a subcutaneous implant model. In mast cell-deficient mice, there was dysregulation of the expected M1 to M2 macrophage transition typically induced by the biomaterial scaffold. Polarization progression deviated in a sex specific manner with an early transition to an M2 profile in female mice, while the male response was unable to properly transition past a pro-inflammatory M1 state. Both were reversed with adoptive mast cell transfer. Further investigation of the later stage immune response in male mice determined a greater sustained pro-inflammatory gene expression profile including the IL-1 cytokine family, IL-6, alarmins, and chemokines. These results highlight mast cells as another important cell type that influences the immune response to pro-regenerative biomaterials. ModEst - Precise estimation of genome size from NGS data Precise estimates of genome sizes are important parameters for both theoretical and practical biodiversity genomics. We present here a fast, easy-to-implement and precise method to estimate genome size from the number of bases sequenced and the mean sequence coverage. To estimate the latter, we take advantage of the fact that a precise estimation of the Poisson distribution parameter lambda is possible from truncated data, restricted to the part of the coverage distribution representing the true underlying distribution. With simulations we could show that reasonable genome size estimates can be gained even from low-coverage (10X), highly discontinuous genome drafts. Comparison of estimates from a wide range of taxa and sequencing strategies with flow-cytometry estimates of the same individuals showed a very good fit and suggested that both methods yield comparable, interchangeable results. Clonal Dynamics of Cancer-Initiating Cells in a Mouse Genetic Mosaic Model of Serous Ovarian Cancer High-grade serous ovarian cancer (HGSOC) is the most lethal gynecologic malignancy. Early detection remains the most promising strategy to decrease mortality and requires a deep understanding of the cancer initiation process. Recent breakthroughs identified the fallopian tube (FT) Pax8+ cells as a cell-of-origin for HGSOC. However, the progressive trajectory from an initiated Pax8+ cell toward malignancy remains elusive. Here, we use a mouse genetic mosaic system called MADM (Mosaic Analysis with Double Markers), which generates sporadic mutant cells labeled with GFP, to trace the clonal dynamics of HGSOC development. Surprisingly, we observed dichotomous clonal expansion kinetics of individual cancer-initiating Pax8+ cells; while the vast majority stalled immediately, a few significantly expanded. Clonal analysis of wildtype cells showed concordant dichotomy, suggesting that a small fraction of Pax8+ cells intrinsically carry higher expansion potential. Furthermore, although all wildtype clones stalled shortly after development, some mutant clones continued to progress, indicating that oncogenic mutations prolong clonal expansion beyond the cell-intrinsic potential. Finally, intra-clonal fate mapping revealed that progressing mutant clones have a higher propensity to produce Pax8+ cells, linking imbalanced cell differentiation to disease progression. In summary, our studies reveal that the earliest clonal expansion during HGSOC initiation emerges when oncogenic mutations exaggerate the intrinsic expansion potential of a rare Pax8+ subset and bias its differentiation potential, warranting future investigation as a readout for early detection and cancer prevention.

SIGNIFICANCE STATEMENTEffective early detection is the most promising approach to manage cancer, but it requires a thorough understanding of cancer initiation and early progression. Serous ovarian cancer has been found to originate from Pax8+ cells in the fallopian tube. To assess the clonal dynamics of cancer-initiating cells, we established a mouse genetic mosaic model that enables clonal analysis of mutant Pax8+ cells in the fallopian tube. We uncovered that the earliest clonal event in ovarian cancer development emerges from a rare subset of stem/progenitor-like cells with intrinsic expansion potential, and that oncogenic mutations exaggerate and prolong this potential to promote cancer initiation. Our findings pinpoint a precise target for future early-detection and cancer prevention studies.

HIGHLIGHTSO_LIOvarian cancer initiation revealed by a mouse genetic mosaic model.
C_LIO_LIDichotomous expansion capacity of initiated cells at the clonal level.
C_LIO_LISynergies between cell-intrinsic properties and onco-mutations fuels clonal expansion.
C_LIO_LIImbalanced differentiation within highly expansive clones leads to Pax8 clusters.
C_LI A novel camera trap design for studying wildlife in mountain glacier ecosystems yields new insight for glacier biodiversity in the Pacific Northwest, USA ContextThe global recession of glaciers and perennial snowfields is reshaping mountain ecosystems. Beyond physical changes to the landscape and altered downstream hydrology, the implications of glacier decline for biodiversity are poorly known. Before predictions can be made about how climate change will affect wildlife in glacier-associated ecosystems, a more thorough accounting of the role that glaciers play in species life histories is needed. However, typical approaches for documenting wildlife presence and behavior--remote camera traps--are difficult to use in glaciated terrain due to limited options for securing them (e.g., no trees) and dramatic seasonal changes in snowpack.

AimsIn this study, we sought to test a novel camera trap designed for glaciated mountain ecosystems. We also aimed to use this approach to gain insight into wildlife and human usage of a mountain glacier in western North America.

MethodsWe deployed an elevational transect of uniquely designed camera traps along the western margin of the Paradise Glacier, a rapidly receding mountain glacier on the south side of Mount Rainier, WA, USA. Our simple camera trap design consisted of a wildlife camera attached to a camouflaged cylindrical cooler filled with snow and rocks.

Key resultsOur camera design proved ideal for a mountain glacier ecosystem and from June to September 2021, we detected at least 16 vertebrate species (seven birds, nine mammals) over 770 trap nights using glacier-associated habitats. Humans, primarily skiers, were the most common species detected, but we also recorded 99 observations of wildlife (birds and mammals). These included three species of conservation concern in Washington: wolverine (Gulo gulo), Cascade red fox (Vulpes vulpes cascadensis), and white-tailed ptarmigan (Lagopus leucura).

ConclusionsCollectively, our results provide proof-of-concept for a novel camera trap design that is ideal for treeless, perennially snow-covered landscapes and revealed a rich diversity of wildlife using mountain glacier habitat in the Pacific Northwest. We highlight the global need for similar studies to better understand the true scale of biodiversity that will be impacted by glacier recession in mountain ecosystems. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Using total RNA quality metrics for time since deposition estimates in degrading bloodstains 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 The Tabula Sapiens: a multiple organ single cell transcriptomic atlas of humans Molecular characterization of cell types using single cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues and tissue specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, the tissue specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type specific RNA splicing was discovered and analyzed across tissues within an individual. Supergene formation: evidence for recombination suppression among multiple functional loci within inversions. Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between set of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining an landmark-free phenotyping algorithm with multivariate genome wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We showed that the supergene controlling the striking wing-pattern polymorphism displayed by this species contains many independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are associated with colour variations only in morphs controlled by inversions, indicating that they were recruited after the formation of these inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. 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.