Tatajuba -- Exploring the distribution of homopolymer tracts Length variation of homopolymeric tracts, which induces phase variation, is known to regulate gene expression leading to phenotypic variation in a wide range of bacterial species. There is no specialised bioinformatics software which can, at scale, exhaustively explore and describe these features from sequencing data. Identifying these is non-trivial as sequencing and bioinformatics methods are prone to introducing artefacts when presented with homopolymeric tracts due to the decreased base diversity. We present tatajuba, which can automatically identify potential homopolymeric tracts and their putative phenotypic impact, allowing for rapid investigation. We use it to detect all tracts in two separate datasets, one of Campylobacter jejuni and one of three Bordetella species, and to highlight those tracts that are polymorphic across samples. With this we confirm homopolymer tract variation with phenotypic impact found in previous studies and additionally find many more with potential variability. The software is written in C and is available under the open source license GNU GPL version 3 from https://github.com/quadram-institute-bioscience/tatajuba.
LazySampling and LinearSampling: Fast Stochastic Sampling of RNA Secondary Structure with Applications to SARS-CoV-2 Many RNAs fold into multiple structures at equilibrium. The classical stochastic sampling algorithm can sample secondary structures according to their probabilities in the Boltzmann ensemble, and is widely used. However, this algorithm, consisting of a bottom-up partition function phase followed by a top-down sampling phase, suffers from three limitations: (a) the formulation and implementation of the sampling phase are unnecessarily complicated; (b) the sampling phase repeatedly recalculates many redundant recursions already done during the partition function phase; (c) the partition function runtime scales cubically with the sequence length. These issues prevent stochastic sampling from being used for very long RNAs such as the full genomes of SARS-CoV-2. To address these problems, we first adopt a hypergraph framework under which the sampling algorithm can be greatly simplified. We then present three sampling algorithms under this framework, among which the LazySampling algorithm is the fastest by eliminating redundant work in the sampling phase via on-demand caching. Based on LazySampling, we further replace the cubic-time partition function by a linear-time approximate one, and derive LinearSampling, an end-to-end linear-time sampling algorithm that is orders of magnitude faster than the standard one. For instance, LinearSampling is 176x faster (38.9s vs. 1.9h) than Vienna RNAsubopt on the full genome of Ebola virus (18,959 nt). More importantly, LinearSampling is the first RNA structure sampling algorithm to scale up to the full-genome of SARS-CoV-2 without local window constraints, taking only 69.2 seconds on its reference sequence (29,903 nt). The resulting sample correlates well with the experimentally-guided structures. On the SARS-CoV-2 genome, LinearSampling finds 23 regions of 15 nt with high accessibilities, which are potential targets for COVID-19 diagnostics and drug design.

See code: https://github.com/LinearFold/LinearSampling
An integrative approach leads to the discovery of a novel anti-leukemic peptide from human milk Chemotherapy in childhood leukemia is associated with late morbidity in leukemic survivors, while certain patient subsets are relatively resistant to standard chemotherapy. It is therefore important to identify new agents with sensitivity and selectivity towards leukemic cells, while having less systemic toxicity. Peptide-based therapeutics has gained much attention during the last few years. Here, we used an integrative workflow combining mass spectrometric peptide library construction, in silico anticancer peptide screening, and in vitro leukemic cell studies to discover a novel anti-leukemic peptide having 3+charges and alpha-helical structure, namely HMP-S7, from human breast milk. HMP-S7 showed cytotoxic activity against four distinct leukemic cell lines in a dose-dependent manner but had no effect on solid malignancies or representative normal cells. HMP-S7 induced leukemic cell death by penetrating the plasma membrane to enter the cytoplasm and cause leakage of lactate dehydrogenase, thus acting in a membranolytic manner. Importantly, HMP-S7 exhibited anti-leukemic effect against patient-derived leukemic cells ex vivo. In conclusion, HMP-S7 is a selective anti-leukemic peptide with promise which requires further validation in preclinical and clinical studies.

TeaserIn silico screening of naturally occurring human milk peptides discovers a new anticancer peptide that kills leukemic cells in vitro and ex vivo. Genomics-based deconvolution of multiplexed transcriptional responses to wood smoke particles defines rapid AHR signaling dynamics. Heterogeneity of respirable particulates and compounds complicates understanding transcriptional responses to air pollution. Here, we address this problem through applying precision nuclear run-on sequencing (PRO-seq) to measure nascent transcription and the assay for transposase-accessible chromatin using sequencing (ATAC-seq) to airway epithelial cells after wood smoke particle (WSP) exposure. We used transcription factor enrichment analysis to agnostically identify temporally distinct roles for the TCF/SRF family, the aryl hydrocarbon receptor (AHR), and NFkB in regulating transcriptional changes induced by WSP. Transcription of canonical targets of the AHR, such as CYP1A1 and AHRR, was robustly increased after just 30 minutes of WSP exposure, and we discovered novel AHR-regulated pathways and targets including the DNA methyltransferase, DNMT3L. Transcription of these genes and associated enhancers rapidly returned to near baseline by 120 minutes. The kinetics of AHR- and NFkB-regulated responses to WSP were distinguishable based on the timing of both transcriptional responses and chromatin remodeling, with induction of several cytokines implicated in maintaining NFkB responses through 120 minutes of exposure. In aggregate, our data establish a direct and primary role for AHR in mediating airway epithelial responses to WSP and identify crosstalk between AHR and NFkB signaling in controlling pro-inflammatory gene expression. This work also defines an integrated genomics-based strategy for deconvoluting multiplexed transcriptional responses to heterogeneous environmental exposures. 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. LabGym: quantification of user-defined animal behaviors using learning-based holistic assessment Quantifying animal behavior is important for many branches of biological research. Current computational tools for quantifying behaviors typically rely on a few pre-defined, simplified features to identify a behavior. However, such an approach underuses the behavioral information and restricts the tools applicability to limited behavior types or animal species. Here we report a new tool, LabGym, for quantifying animal behaviors without such limitations. Combining a novel approach for effective evaluation of animal motion with customizable convolutional recurrent networks for capturing spatiotemporal details, LabGym provides holistic behavioral assessment and accurately identify user-defined animal behaviors without restrictions on behavior types or animal species. It then calculates various parameters to provide quantitative measurements of the intensity and dynamics of each behavior. LabGym requires neither any intermediate step that causes information loss nor programming knowledge from users for post-hoc analysis. It tracks multiple animals simultaneously in various experimental settings for high-throughput and versatile analysis, and provides users a way to generate visualizable behavioral datasets that are valuable resources for research in neuroscience, behavioral science, and computer science. We demonstrate its efficacy in capturing subtle behavioral changes in animals ranging from soft-bodied invertebrates to mammals. Fixed or random? On the reliability of mixed-effects models for a small number of levels in grouping variables O_LIBiological data are often intrinsically hierarchical. Due to their ability to account for such dependencies, mixed-effects models have become a common analysis technique in ecology and evolution. While many questions around their theoretical foundations and practical applications are solved, one fundamental question is still highly debated: When facing a low number of levels should we model a grouping (blocking, clustering) variable as a random or fixed effect? In such situation, the variance of the random effect is imprecise, but whether this affects the statistical properties of the population effect is unclear.
C_LIO_LIHere, we analyzed the consequences of including a grouping variable as fixed or random effect in the correctly specified and other possible miss-specified models (too simple or too complex models) for data with small number of levels (2 - 8). For all these options, we calculated type I error rates and power. Moreover, we show how these statistical properties change with the study design.
C_LIO_LIWe found that the model choice does not influence the statistical properties of the population effect when the effect is the same at all levels However, if an ecological effect differs among levels, using a random slope and intercept model, and switching to a fixed-effect model only in case of a singular fit, avoids overconfidence in the results. Additionally, power and type I error are strongly influenced by the number of and difference between levels.
C_LIO_LIWe conclude that inferring the correct random effect structure is of high importance to get correct statistical properties. When in doubt, we recommend starting with the simpler model and using model diagnostics to identify missing components. When having identified the correct structure, we encourage to start with a mixed-effects model independent of the number of levels and switch to a fixed-effect model only in case of a singular fit. With these recommendations, we allow for more informative choices about study design and data analysis and thus make ecological inference with mixed-effects models more robust for small number of levels.
C_LI Differential fate of acellular vascular scaffolds in vivo involves macrophages and T-cell subsets Biological scaffold or implant is a popular choice for the preparation of tissue-engineered organs and has the potential to address donor shortage in the clinics. However, biological scaffolds prepared by physical or chemical agents cause damage to the extracellular matrix by potentially inducing immune responses after implantation. The current study explores an alternative route for the preparation of acellular scaffolds and explores the fate of the prepared scaffolds in a milieu of immune cells following implantation without using immunosuppressant. Using the syngeneic (Lewis male-Lewis female) and allogeneic (Brown Norway male-Lewis female) models and different tissue routes (subcutaneous vs omentum) for transplantation, normal blood vascular scaffolds were implanted which was converted to acellular vascular scaffolds by in vivo natural decellularization at the end of 2 months of observation. We also prepared chemically decellularized acellular scaffolds from normal untreated blood vascular scaffolds using a cocktail of chemicals which was also similarly placed in subcutaneous and omentum sites. Here, we applied in-depth quantitative proteomics along with histology and image analysis to comprehensively describe and compare the proteome of the natural and chemically decellularized scaffold. Our data confirm that site-specific advantages exist in modulating the ECM and regulating the immune responses (macrophage and T cells) following implantation, which possibly led to the production of an acellular scaffold (natural decellularization) under in vivo conditions. The current approach opens up the possibility to create tailor-made acellular scaffolds to build functional blood vessels. In addition, the identification of different tissue sites facilitates differential immune response against the scaffolds. This study provides a rich resource aimed toward an enhanced mechanistic understanding to study immune responses under similar settings in the field of transplantation and regenerative medicine.

Impact statementThe development of a scaffold helps in the preparation of a functional organ in the clinics. In the current study, we prepared an acellular vascular scaffold by utilizing site specific tissue changes and vis-a-vis compared with a conventionally chemically prepared biological scaffold at genomic and protein level, which helped us to identify immunological trigger following implantation. The current study which was carried out without any immunosuppressive agents could help to establish (a) alternative strategies for preparing biological scaffolds as well as (b) implantable sites as potential bioreactors to circumvent any adverse immune reactions for acceptance of the scaffold/implant post implantation. Post-translational flavinylation is associated with diverse extracytosolic redox functionalities throughout bacterial life Disparate redox activities that take place beyond the bounds of the prokaryotic cell cytosol must connect to membrane or cytosolic electron pools. Proteins post-translationally flavinylated by the enzyme ApbE mediate electron transfer in several characterized extracytosolic redox systems but the breadth of functions of this modification remains unknown. Here we present a comprehensive bioinformatic analysis of 31,910 prokaryotic genomes that provides evidence of extracytosolic ApbEs within ~50% of bacteria and the involvement of flavinylation in numerous uncharacterized biochemical processes. By mining flavinylation-associated gene clusters, we identify five protein classes responsible for transmembrane electron transfer and two domains of unknown function (DUF2271 and DUF3570) that are flavinylated by ApbE. We observe flavinylation/iron transporter gene colocalization patterns that implicate functions in iron reduction and assimilation. We find associations with characterized and uncharacterized respiratory oxidoreductases that highlight roles of flavinylation in respiratory electron transport chains. Finally, we identify interspecies gene cluster variability consistent with flavinylation/cytochrome functional redundancies and discover a class of "multi-flavinylated proteins'' that may resemble multiheme cytochromes in facilitating longer distance electron transfer. These findings provide key mechanistic insight into an important facet of bacterial physiology and establish flavinylation as a functionally diverse mediator of extracytosolic electron transfer.

Impact StatementBioinformatic and biochemical studies find that covalently bound flavins are common and participate in wide-ranging extracytosolic redox activities throughout bacterial life. Bio-hybrid micro-swimmers propelled by flagella isolated from C. reinhardtii Bio-hybrid micro-swimmers, composed of biological entities integrated with synthetic constructs, actively transport cargo by converting chemical energy into mechanical work in a fluid at low Reynolds number, where viscous drag dominates over inertia. Here, using isolated and demembranated flagella from green algae Chlamydomonas reinhardtii (C. reinhardtii), we build efficient axonemally-driven micro-swimmers that consume ATP to propel micron-sized beads. Depending on the calcium concentration, we observed two main classes of motion: Whereas beads move along curved trajectories at calcium concentrations below 0.03 mM, they are propelled along straight paths when the calcium concentration increases. In this regime, they reached velocities of approximately 20 m/sec, comparable to human sperm velocity in vivo. We relate this transition to the properties of beating axonemes, in particular the reduced static curvature with increasing calcium concentration. To quantify the motion, we used mode decomposition of the flagellar waveform, and we studied both analytically and numerically the propulsion of the bead as a function of the axonemal waveform and bead-axoneme attachment geometry. While our analysis semi-quantitatively describes the experimental results, it also reveals the existence of a counter-intuitive propulsion regime where the speed of the axonemally-driven bead increases with the size of the bead. Moreover, we demonstrated that asymmetric, sideways attachment of the axoneme to the bead can also contribute to the rotational velocity of the micro-swimmer. The uncovered mechanism has potential applications in the fabrication of synthetic micro-swimmers, and in particular, bio-actuated medical micro-robots for targeted drug delivery. State-dependent network interactions differentially gate sensory input at the motor and command neuron level in Caenorhabditis elegans Neural responses are influenced by both external stimuli and internal network states. While network states have been linked to behavioral and stimulus states, little is known about how sensory inputs are filtered by whole-brain activity to affect motor and command neurons. Here, we recorded whole-brain activity of Caenorhabditis elegans experiencing bacterial food stimuli, and modeled how sensory inputs affect motor and command neurons in a network state-dependent manner. First, we classified active neurons into six functional clusters: two sensory neuron clusters (ON, OFF), and four motor/command neuron clusters (AVA, RME, SMDD, SMDV). Using encoding models, we found that ON and OFF sensory neurons that respond to onset and removal of bacteria, respectively, employ different adaptation strategies. Next, we used decoding models to show that bacterial onset and removal differentially drive AVA and RME cluster activity. To explore state-dependent effects on AVA and RME clusters, we developed a model that identified network states and fitted submodels for each state to predict how each of the six functional clusters drive AVA and RME cluster activity. We also identified network states in which AVA and RME clusters were either largely unperturbed by or receptive to bacterial sensory input. Furthermore, this model allowed us to disentangle the state-dependent contributions of stimulus timescales and bacterial content to neural activity. Collectively, we present an interpretable approach for modeling network dynamics that goes beyond implication of neurons in particular states, and moves toward explicitly dissecting how neural populations work together to produce state dependence.

Significance StatementA major function of the brain is to transform sensory information into behavior. As the first receiver of sensory input, sensory neuron activity is often most correlated with stimulus features. However, this high-fidelity representation of sensory input becomes diluted as it travels to downstream neurons, where sensory information is integrated with network activity. By the time sensory information reaches motor neurons, it is often difficult to dissociate the influence of sensory input from the influence of network activity. Here, we describe a method that is fully interpretable such that we can show how neural populations on a whole-brain scale interact to produce network states. From there, we can attribute motor neuron activity to network history and sensory input. The Ensemble of Gene Regulatory Networks at Mutation-Selection Balance The evolution of diverse phenotypes both involves and is constrained by molecular interaction networks. When these networks influence patterns of expression, we refer to them as gene regulatory networks (GRNs). Here, we develop a quasi-species model of GRN evolution. With this model, we prove that-across a broad spectrum of viability and mutation functions-the dynamics converge to a stationary distribution over GRNs. Next, we show from first principles how the frequency of GRNs at equilibrium will be proportional to each GRNs eigenvector centrality in the genotype network. Finally, we determine the structural characteristics of GRNs that are favored in response to a range of selective environments and mutational constraints. Our work connects GRN evolution to quasi-species models, and thus can provide a mechanistic explanation for the topology of GRNs experiencing various evolutionary forces. Coordination of genome replication and anaphase entry by rDNA copy number in S. cerevisiae Eukaryotes maintain hundreds of copies of ribosomal DNA (rDNA), many more than required for ribosome biogenesis, suggesting a yet undefined role for large rDNA arrays outside of ribosomal RNA synthesis. We demonstrate that reducing the Saccharomyces cerevisiae rDNA array to 35 copies, which is sufficient for ribosome function, shifts rDNA from being the latest replicating region in the genome to one of the earliest. This change in replication timing results in delayed genome-wide replication and classic replication defects. We present evidence that the requirement for rDNA to replicate late, which is conserved among eukaryotes, also coordinates the completion of genome replication with anaphase entry through the proper sequestration of the mitotic exit regulator Cdc14p in the rDNA-containing nucleolus. Our findings suggest that, instead of being a passive repetitive element, the large late-replicating rDNA array plays an active role in genome replication and cell cycle control. Human papillomavirus integration transforms chromatin to drive oncogenesis Human papillomavirus (HPV) drives almost all cervical cancers and up to [~]70% of head and neck cancers. Frequent integration into the host genome occurs only for tumourigenic strains of HPV. We hypothesized that changes in the epigenome and transcriptome contribute to the tumourigenicity of HPV. We found that viral integration events often occurred along with changes in chromatin state and expression of genes near the integration site. We investigated whether introduction of new transcription factor binding sites due to HPV integration could invoke these changes. Some regions within the HPV genome, particularly the position of a conserved CTCF binding site, showed enriched chromatin accessibility signal. ChIP-seq revealed that the conserved CTCF binding site within the HPV genome bound CTCF in 4 HPV+ cancer cell lines. Significant changes in CTCF binding pattern and increases in chromatin accessibility occurred exclusively within 100 kbp of HPV integration sites. The chromatin changes co-occurred with out-sized changes in transcription and alternative splicing of local genes. We analyzed the essentiality of genes upregulated around HPV integration sites of The Cancer Genome Atlas (TCGA) HPV+ tumours. HPV integration upregulated genes which had significantly higher essentiality scores compared to randomly selected upregulated genes from the same tumours. Our results suggest that introduction of a new CTCF binding site due to HPV integration reorganizes chromatin and upregulates genes essential for tumour viability in some HPV+ tumours. These findings emphasize a newly recognized role of HPV integration in oncogenesis. Reconstructing SNP Allele and Genotype Frequencies from GWAS Summary Statistics The emergence of genomewide association studies (GWAS) has led to the creation of large repositories of human genetic variation, creating enormous opportunities for genetic research and worldwide collaboration. Methods that are based on GWAS summary statistics seek to leverage such records, overcoming barriers that often exist in individual-level data access while also offering significant computational savings. Here, we propose a novel framework that can reconstruct allelic and genotypic counts/frequencies for each SNP from case-control GWAS summary statistics. Our framework is simple and efficient without the need of any complicated underlying assumptions. Illustrating the great potential of this framework we also propose three summary-statistics-based applications implemented in a new software package (ReACt): GWAS meta-analysis (with and without sample overlap), case-case GWAS, and, for the first time, group polygenic risk score (PRS) estimation. We evaluate our methods against the current state-of-the-art on both synthetic data and real genotype data and show high performance in power and error control. Our novel group PRS method based on summary statistics could not be achieved prior to our proposed framework. We demonstrate here the potential applications and advantages of this approach. Our work further highlights the great potential of summary-statistics-based methodologies towards elucidating the genetic background of complex disease and opens up new avenues for research. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life. Adaptive anti-tumor immunity is orchestrated by a population of CCL5-producing tissue-resident NK cells T cell-directed cancer immunotherapy often fails to generate lasting tumor control. Harnessing additional effectors of the immune response against tumors may strengthen the clinical benefit of immunotherapies. Here, we demonstrate that therapeutic targeting of the IFN{gamma}-IL-12 pathway relies on the ability of a population of tissue-resident NK (trNK) cells to orchestrate an anti-tumor microenvironment. Particularly, utilizing an engineered adenoviral platform, we show that paracrine IL-12 enhances functional DC-CD8 T cell interactions to generate adaptive anti-tumor immunity. This effect depends on the abundance of trNK cells and specifically their capacity to produce the cDC1-chemoattractant CCL5. Failure to respond to IL-12 and other IFN{gamma}-inducing therapies such as immune checkpoint blockade in tumors with low trNK cell infiltration could be overcome by intra-tumoral delivery of CCL5. Our findings reveal a novel barrier for T cell-focused therapies and offer mechanistic insights into how T cell-NK cell-DC crosstalk can be enhanced to promote anti-tumor immunity and overcome resistance.

SignificanceWe identified the lack of CCL5-producing, tissue-resident NK (trNK) cells as a barrier to T cell-focused therapies. While IL-12 induces anti-tumoral DC-T cell crosstalk in trNK cellrich tumors, resistance to IL-12 or anti-PD-1 in trNK cellpoor tumors can be overcome by the additional delivery of CCL5. Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Fatal COVID-19 ARDS associated with incomplete AEC1 differentiation from the transitional state without senescence or fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Robust and Remote Measurement of Heart Rate Based on A Surveillance Camera Arrhythmia is a marked symptom of many cardiovascular diseases (CVDs). The accurate and timely detection of heart rate can greatly reduce the harm of arrhythmia to people. However, it is still a challenge to robustly and remotely measure heart rate in daily life due to the changing environmental conditions during measurement, such as the varying light intensity, the movement of people, and the uncertain distance between the sensor to people. In this study, we propose a method to accurately measure human heart rate within a distance of 4.5 meters under different light intensities by simply using a surveillance camera. After a 20-second color video of a persons hand is captured by the camera, a method based on the Fast Fourier Transform (FFT) algorithm is designed to extract the blood volume pulse wave to calculate the heart rate. According to the comparison between the real heart rate and results measured by electrocardiography (ECG), the proposed method achieves an accuracy of 95.8% when the measurement is performed within a distance of 4.5 meters and 90% when within 5.0 meters. Our experiments show that when the illuminance varies between 100-1000 lux (lighting level indoor), we still get the correct results. Our experiments also demonstrate that the proposed method accurately obtain heart rate even when the light intensity is below 32 lux (300-500 lux in a workplace environment). The methods strong adaptability to changing environmental conditions makes it applicable to many scenarios, such as homes of the elderly, classrooms, and other public spaces. The Temporal Statistics of Musical Rhythm across Western Genres: An Amplitude Modulation Phase Hierarchy Model Statistical learning by the human brain plays a core role in the development of cognitive systems like language and music. Both music and speech have structured inherent rhythms, however the acoustic sources of these rhythms are debated. Theoretically, rhythm structures in both systems may be related to a novel set of acoustic statistics embedded in the amplitude envelope, statistics originally revealed by modelling childrens nursery rhymes. Here we apply similar modelling to explore whether the amplitude modulation (AM) timescales underlying rhythm in music match those in child-directed speech (CDS). Utilising AM-driven phase hierarchy modelling previously applied to infant-directed speech (IDS), adult-directed speech (ADS) and CDS, we test whether the physical stimulus characteristics that yield speech rhythm in IDS and CDS describe rhythm in music. Two models were applied. One utilized a low-dimensional representation of the auditory signal adjusted for known mechanisms of the human cochlear, and the second utilized probabilistic amplitude demodulation, estimating the modulator (envelope) and carriers using Bayesian inference. Both models revealed a similar hierarchically-nested temporal modulation structure across Western musical genres and instruments. Core bands of AM and spectral patterning matched prior analyses of IDS and CDS, and music showed strong phase dependence between slower bands of AMs, again matching IDS and CDS. This phase dependence is critical to the perception of rhythm. Control analyses modelling other natural sounds (wind, rain, storms, rivers) did not show similar temporal modulation structures and phase dependencies. We conclude that acoustic rhythm in language and music has a shared statistical basis. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Fertility costs of cryptic viral infections in a model social insect Declining insect populations emphasize the importance of understanding the drivers underlying reductions in insect fitness. Here, we investigated viruses as a threat to social insect reproduction, using honey bees as a model species. We report that in a sample of N = 93 honey bee (Apis mellifera) queens from nine beekeeping operations across a wide geographic range, high levels of natural viral infection are associated with decreased ovary mass. We confirmed this finding in an independent sample of N = 54 queens. Failed (poor quality) queens displayed higher levels of viral infection, reduced sperm viability, smaller ovaries, and altered ovary protein composition compared to healthy queens. We experimentally infected queens with Israeli acute paralysis virus (IAPV) and found that the ovary masses of IAPV-injected queens were significantly smaller than control queens, demonstrating a causal relationship between viral infection and ovary size. Queens injected with IAPV also had significantly lower expression of vitellogenin, the main source of nutrition deposited into developing oocytes, and higher levels of heat-shock proteins (HSPs), which are part of the honey bees antiviral response. This work together shows that viral infections occurring naturally in the field are compromising queen reproductive success. Neprosin belongs to a new family of glutamic peptidase family based on in silico evidence Neprosin was first discovered in the insectivorous tropical pitcher plants of Nepenthes species as a novel protease with prolyl endopeptidase (PEP) activity. Neprosin has two uncharacterized domains of neprosin activation peptide and neprosin. A previous study has shown neprosin activity in hydrolyzing proline-rich gliadin, a gluten component that triggers celiac disease. In this study, we performed in silico structure-function analysis to investigate the catalytic mechanism of neprosin. Neprosin sequences lack the catalytic triad and motifs of PEP family S9. Protein structures of neprosins from Nepenthes x ventrata (NvNpr) and N. rafflesiana (NrNpr1) were generated by ab initio methods and comparatively assessed to obtain high-quality models. Structural alignment of models to experimental structures in the Protein Data Bank (PDB) found a high structural similarity to glutamic peptidases. Further investigations reveal other resemblances to the glutamic peptidases with low optimum pH that activates the enzyme via autoproteolysis for maturation. Two highly conserved glutamic acid residues, which are stable according to the molecular dynamics simulation, can be found at the active site of the substrate cleft. Protein docking demonstrated that mature neprosins bind well with potent antigen I-gliadin at the putative active site. Taken together, neprosins represent a new glutamic peptidase family, with a putative catalytic dyad of two glutamic acids. This study illustrates a hypothetical enzymatic mechanism of the neprosin family and demonstrates the useful application of an accurate ab initio protein structure prediction in the structure-function study of a novel protein family. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. Disruption of the microphysiological niche alters matrix deposition and causes loss of mitochondrial Ca2+ control in skeletal muscle fibers Cells rapidly lose their physiological phenotype upon disruption of their extracellular matrix (ECM)-intracellular cytoskeleton interactions. Here, we investigated acute effects of ECM disruption on cellular and mitochondrial morphology, transcriptomic signatures, and Ca2+ handling in adult mouse skeletal muscle fibers. Adult skeletal muscle fibers were isolated from mouse toe muscle either by collagenase-induced dissociation of the ECM or by mechanical dissection that leaves the contiguous ECM intact. Experiments were generally performed four hours after cell isolation. At this time, there were striking differences in the gene expression patterns between fibers isolated with the two methods; 24h after cell isolation, enzymatically dissociated fibers had transcriptomic signatures resembling dystrophic phenotypes. Mitochondrial appearance was grossly similar in the two groups, but 3D electron microscopy revealed shorter and less branched mitochondria in enzymatically dissociated than in mechanically dissected fibers. Similar increases in free cytosolic [Ca2+] during repeated tetanic stimulation were accompanied by marked mitochondrial Ca2+ uptake only in enzymatically dissociated muscle fibers. The aberrant mitochondrial Ca2+ uptake was partially prevented by the mitochondrial Ca2+ uniporter inhibitor Ru360 and by cyclosporine A and NV556, which inhibit the mitochondrial protein Ppif (also called cyclophilin D). Importantly, inhibition of Ppif with NV556 significantly improved survival of mice with mitochondrial myopathy in which muscle mitochondria take up excessive amounts of Ca2+ even with an intact ECM. In conclusion, skeletal muscle fibers isolated by collagenase-induced dissociation of the ECM display aberrant mitochondrial Ca2+ uptake, which involves a Ppif-dependent mitochondrial Ca2+ influx resembling that observed in mitochondrial myopathies.