GEM-DeCan: Improved tumor immune microenvironment profiling through novel gene expression and DNA methylation signatures predicts immunotherapy response Quantifying the proportion of the different cell types present in tumor biopsies remains a priority in cancer research. So far, a number of deconvolution methods have emerged for estimating cell composition using reference signatures, either based on gene expression or on DNA methylation from purified cells. These two deconvolution approaches could be complementary to each other, leading to even more performant signatures, in cases where both data types are available. However, the potential relationship between signatures based on gene expression and those based on DNA methylation remains underexplored.

Here we present five new deconvolution signature matrices, based on DNA methylation or RNAseq data, which can estimate the proportion of immune cells and cancer cells in a tumour sample. We test these signature matrices on available datasets for in-silico and in-vitro mixtures, peripheral blood, cancer samples from TCGA, bone marrow from multiple myeloma patients and a single-cell melanoma dataset. Cell proportions estimates based on deconvolution performed using our signature matrices, implemented within the EpiDISH framework, show comparable or better correlation with FACS measurements of immune cell-type abundance and with various estimates of cancer sample purity and composition than existing methods.

Using publicly available data of 3D chromatin structure in haematopoietic cells, we expanded the list of genes to be included in the RNAseq signature matrices by considering the presence of methylated CpGs in gene promoters or in genomic regions which are in 3D contact with these promoters. Our expanded signature matrices have improved performance compared to our initial RNAseq signature matrix. Finally, we show the value of our signatures in predicting patient response to immune checkpoint inhibitors in three melanoma and one bladder cancer cohorts, based on bulk tumour sample gene expression.

We also provide GEM-DeCan: a snakemake pipeline, able to run an analysis from raw sequencing data to deconvolution based on various gene expression signature matrices, both for bulk RNASeq and DNA methylation data. The code for producing the signature matrices and reproducing all the figures of this paper is available on the GEM-DeCan repository.
Tourmaline: a containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake BackgroundAmplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results.

FindingsWe developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files--a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences--it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm, assigns taxonomy to the resulting representative sequences, performs analyses of taxonomic, alpha, and beta diversity, and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S rRNA amplicon dataset from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks.

ConclusionsAutomated workflows like Tourmaline enable rapid analysis of environmental and biomedical amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline.
DNA supercoiling-mediated collective behavior of co-transcribing RNA polymerases Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or "push" between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate via the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the torsional response of supercoiled DNA. Using stochastic simulations, we demonstrate long-range cooperation between co-transcribing RNAPs mediated by DNA supercoiling. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern a key cellular process. The "Alu-ome" shapes the epigenetic environment of regulatory elements controlling cellular defense Promoters and enhancers are sites of transcription initiation (TSSs) and carry active histone modifications, including H3K4me1, H3K4me3, and H3K27ac. Yet, the principles governing the boundaries of such regulatory elements are still poorly characterized. Alu elements are good candidates for a boundary function, being highly abundant in gene-rich regions, while essentially excluded from regulatory elements. Here, we show that the interval from the TSS to the first upstream Alu accommodates all H3K4me3 and most H3K27ac marks, while excluding DNA methylation. Remarkably, the average length of this intervals greatly varies in-between tissues, being longer in stem-and shorter in immune-cells. The very shortest TSS-to-Alu intervals were observed at promoters active in T cells, particularly at immune genes, correlating with RNA polymerase II transcription through the first Alu and accumulation of H3K4me1 signal on this first Alu. Finally, DNA methylation at first-Alus was found to evolved with age, regressing from young to middle-aged, then recovering later in life. Thus, the first Alus upstream of TSSs appear as dynamic boundaries marking the transition from DNA methylation to active histone modifications at regulatory elements, while also participating in the recording of immune gene transcriptional events by positioning H3K4me1-modified nucleosomes. Sufu regulation of Hedgehog signaling in P19 cells is required for proper astrocyte differentiation Hedgehog signaling is essential for vertebrate development, however, less is known about the negative regulators that influence this pathway during the differentiation of cell fates. Using the mouse P19 embryonal carcinoma cell model, Suppressor of Fused (SUFU), a negative regulator of the Hedgehog pathway, was investigated during retinoic acid-induced neural differentiation. We found Hedgehog signaling was activated in the early phase of neural differentiation but was inactive during terminal differentiation of neurons and astrocytes. This early activation was required for neural differentiation however, it alone was not sufficient to induce neural lineages. SUFU, which regulates signaling at the level of GLI, remained relatively unchanged during the differentiation process, but its loss through CRISPR-Cas9 gene editing resulted in ectopic expression of Hedgehog target genes. Interestingly, these SUFU-deficient cells were unable to differentiate in the absence of retinoic acid, but when differentiated using retinoic acid, they showed delayed and decreased astrocyte differentiation; neuron differentiation was unaffected. Ectopic activation of Hh target genes in SUFU-deficient cells remained throughout retinoic acid-induced differentiation and this phenotype was accompanied by the loss of GLI3, despite the presence of the Gli3 message. Thus, the study indicates the proper timing and proportion of astrocyte differentiation requires SUFU, and its normal regulation of GLI3 to maintain Hh signaling in an inactive state. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Missing the forest for the trees: Do seaweed ecosystems mitigate atmospheric CO2 emissions? Global seaweed carbon sequestration estimates are currently taken as the fraction of the net primary production (NPP) exported to the deep ocean. However, this perspective does not account for CO2 from the consumption of external subsidies. Here we clarify: i) the role of export relative to seaweed net ecosystem production (NEP) for a closed system and one more likely open to subsidies; ii) the importance of subsidies by compiling published estimates of NEP from seaweed-dominated ecosystems; and iii) discuss their impact on the global seaweed net carbon balance and other sequestration constraints as a mitigation service. Examples of seaweed NEP (n = 18) were sparse and variable. Nevertheless, the average NEP (-9.2mmol C m-2 day-1 SE {+/-} 11.6) suggested that seaweed ecosystems are a C source, becoming increasingly heterotrophic as their export is consumed. Critically, mitigation of greenhouse gas emissions was mixed relative to their replacement or baseline states, and where CO2 is supplied independently of organic metabolism and atmospheric exchange we caution a sole reliance on NEP or NPP. This will ensure a more accurate seaweed mitigation assessment, one that does exceed their capacity and is effective within a compliance and carbon trading scheme. Climate change and chill accumulation: implications for tree fruit production in cold-winter regions Winter chill accumulation is critical for the productivity and profitability of perennial tree fruit systems. Several studies have quantified the impacts of global warming on chill accumulation in the warmer production regions of the world, where insufficient chill events occur and their frequency is increasing. In contrast, we focus on a region with relatively cold winters-the Pacific Northwest United States (PNW)-where insufficient chill events are currently absent, and quantify the potential for introduction of these risks under climate change. Results identified spatial variation within the PNW, with chill accumulation projected to increase in some areas but decrease in others. There was also spatiotemporal variation in the driving factors of changes to chill accumulation. Even with decreases in chill accumulations, there are likely minimal issues with insufficient chill accumulation. However, delayed chill accumulation in combination with advances in the onset of heat accumulation can potentially shift the region from one where spring phenology is primarily forcing-driven to one where the dynamic interplay between chilling and forcing processes become important. These interactions might create production risks for varieties with high chill requirements, post mid-21st-century under high emissions scenarios. Future work should focus on understanding, modeling, and projecting responses across these overlapping chilling and forcing processes. Additionally, given significant spatial differences across a relatively small geographic range, it is also critical to understand and model these dynamics at a local landscape resolution for regions such as the PNW. Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin BackgroundEvolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation at many enhancers, distal regulatory elements that control spatio-temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific, code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

ResultsWe present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that is able to predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Then, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity and liver enhancers associated with liver regeneration have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

ConclusionThe framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements. Understanding the Spin Crossover Dynamical Effects of the Dioxygen Binding and Activation on HOD Enzyme For the cofactor-free 1-H-3-hydroxy-4-oxoquinaldine-2,4-dioxygenase (HOD), the dioxygen (O2) dependent steps are rate-limiting along with a spin state crossover to the singlet spin state. Here, the primary triplet O2 molecule activation on the 2-methyl-3-hydroxy-4(1H)-quinolone (MHQ) is investigated, and the catalytic role of the intersystem crossing effects is highlighted by directly comparing results from the Born-Oppenheimer dynamics and non-adiabatic surface hopping dynamics. This work confirms non-adiabatic dynamical effects are essential to modulate the O2 activation on the substrate MHQ. The time scale of the equilibration and conversion from triplet to singlet state should be in the range of a few hundreds of femtoseconds. We hope this work provides us a fresh look at the underlying physics of dioxygen activation reactions involving more than one spin state. An Essential Role of UBXN3B for B Lymphopoiesis Hematopoiesis is finely regulated to enable timely production of the right numbers and types of mature immune cells to maintain tissue homeostasis. Dysregulated hematopoiesis may compromise antiviral immunity and/or exacerbate immunopathogenesis. Herein, we report an essential role of UBXN3B in maintenance of hematopoietic homeostasis and restriction of immunopathogenesis during respiratory viral infection. Ubxn3b deficient (Ubxn3b-/-) mice are highly vulnerable to SARS-CoV-2 and influenza A infection, characterized by more severe lung immunopathology, lower virus-specific IgG, significantly fewer B cells, but more myeloid cells than Ubxn3b+/+ littermates. This aberrant immune compartmentalization is recapitulated in uninfected Ubxn3b-/- mice. Mechanistically, UBXN3B controls precursor B-I (pre-BI) transition to pre-BII and subsequent proliferation in a cell-intrinsic manner, by maintaining BLNK protein stability and pre-BCR signaling. These results reveal an essential role of UBXN3B for the early stage of B cell development. Increasing triacylglycerol formation and lipid storage by unsaturated lipids protects renal proximal tubules in diabetes In diabetic patients, dyslipidemia frequently contributes to organ damage such as diabetic kidney disease (DKD). DKD is associated with excessive renal deposition of triacylglycerol (TAG) in lipid droplets (LD). Yet, it is unclear whether LDs play a protective or damaging role and how this might be influenced by dietary patterns. By using a diabetes mouse model, we find here that high fat diet enriched in the unsaturated oleic acid (OA) caused more lipid storage in LDs in renal proximal tubular cells (PTC) but less tubular damage than a corresponding butter diet with the saturated palmitic acid (PA). Mechanistically, we identify endoplasmic reticulum (ER) stress as the main cause of PA-induced PTC injury. ER stress is caused by elevated cellular levels of saturated TAG precursors and to higher membrane order in the ER. The resulting cell death is preceded by a transcriptional rewiring of phospholipid metabolism. Simultaneous addition of OA rescues the cytotoxic effects by normalizing membrane order and by increasing the total TAG amount. The latter also stimulates the formation of LDs that in turn can release unsaturated lipids upon demand by lipolysis. Our study thus clarifies mechanisms underlying PA-induced cell stress in PTCs and emphasizes the importance of olive oil for the prevention of DKD. MODEL PREDICTIVE CONTROL OF CANCER CELLULAR DYNAMICS: A NEW STRATEGY FOR THERAPY DESIGN Recent advancements in Cybergenetics have led to the development of new computational and experimental platforms that enable to robustly steer cellular dynamics by applying external feedback control. Such technologies have never been applied to regulate intracellular dynamics of cancer cells. Here, we show in silico that adaptive model predictive control (MPC) can effectively be used to steer signalling dynamics in Non-Small Cell Lung Cancer (NSCLC) cells to resemble those of wild-type cells, and to support the design of combination therapies. Our optimisation-based control algorithm enables tailoring the cost function to force the controller to alternate different drugs and/or reduce drug exposure, minimising both drug-induced toxicity and resistance to treatment. Our results pave the way for new cybergenetics experiments in cancer cells for drug combination therapy design. 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. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. Integrated analysis of cervical squamous cell carcinoma cohorts from three continents reveals conserved subtypes of prognostic significance. Human papillomavirus (HPV)-associated cervical cancer represents one of the leading causes of cancer death worldwide. Although low-middle income countries are disproportionately affected, our knowledge of the disease predominantly originates from populations in high-income countries. Using the largest multi-omic analysis of cervical squamous cell carcinoma (CSCC) to date, totalling 643 tumours and representing patient populations from the USA, Europe and Sub-Saharan Africa, we identify two CSCC subtypes (C1 and C2) with differing prognosis. C1 tumours are largely HPV16-driven, display increased cytotoxic T-lymphocyte infiltration and frequently harbour PIK3CA and EP300 mutations. C2 tumours are associated with shorter overall survival, are frequently driven by HPVs from the HPV18-containing alpha-7 clade, harbour alterations in the Hippo signalling pathway and increased expression of immune checkpoint genes, B7-H3 (also known as CD276) and NT5E (also known as CD73) and PD-L2 (also known as PDCD1LG2). In conclusion, we identify two novel, therapy-relevant CSCC subtypes that share the same defining characteristics across three geographically diverse cohorts. Ineffectual AEC1 differentiation from KRT8hi transitional state without fibrosis is associated with fatal COVID-19 ARDS 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. Sniffing Out New Friends: Similarity in Body-Odor Predicts the Quality of Same-Sex Non-Romantic Dyadic Interactions Most are familiar with the notion of socially "clicking" with someone, namely sensing an immediate bond that can lead to strong and often long-lasting friendships. The mechanisms underlying such rapid bonding remain unclear. Given that body-odor similarity is a critical cue for social interaction in non-human mammals, we tested the hypothesis that body-odor similarly contributes to bonding in same-sex non-romantic human dyads. We observed that objective ratings obtained with an electronic nose, and subjective ratings obtained from human smellers, converged to suggest that click-friends smell more similar to each other than random dyads. Remarkably, we then found that we could use the electronic nose to predict which strangers would later form better dyadic interactions. Thus, humans may literally sniff-out new friends based on similarities in body-odor. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. Evaluation of methods for estimating coalescence times using ancestral recombination graphs The ancestral recombination graph (ARG) is a structure that describes the joint genealogies of sampled DNA sequences along the genome. Recent computational methods have made impressive progress towards scalably estimating whole-genome genealogies. In addition to inferring the ARG, some of these methods can also provide ARGs sampled from a defined posterior distribution. Obtaining good samples of ARGs is crucial for quantifying statistical uncertainty and for estimating population genetic parameters such as effective population size, mutation rate, and allele age. Here, we use standard neutral coalescent simulations to benchmark the estimates of pairwise coalescence times from three popular ARG inference programs: ARGweaver, Relate, and tsinfer+tsdate. We compare 1) the true coalescence times to the inferred times at each locus; 2) the distribution of coalescence times across all loci to the expected exponential distribution; 3) whether the sampled coalescence times have the properties expected of a valid posterior distribution. We find that inferred coalescence times at each locus are most accurate in ARGweaver, and often more accurate in Relate than in tsinfer+tsdate. However, all three methods tend to overestimate small coalescence times and underestimate large ones. Lastly, the posterior distribution of ARGweaver is closer to the expected posterior distribution than Relates, but this higher accuracy comes at a substantial trade-off in scalability. The best choice of method will depend on the number and length of input sequences and on the goal of downstream analyses, and we provide guidelines for the best practices. Trypanosomal variant surface glycoprotein expression in human African trypanosomiasis patients Trypanosoma brucei gambiense, an extracellular protozoan parasite, is the primary causative agent of human African trypanosomiasis. T. b. gambiense is endemic to West and Central Africa, where it is transmitted by the bite of infected tsetse flies. In the bloodstream of an infected host, the parasite evades antibody recognition by altering the Variant Surface Glycoprotein (VSG) that forms a dense coat on its cell surface through a process known as antigenic variation. Each VSG has a variable N-terminal domain exposed to the host and a less variable C-terminal domain that is partially hidden from host antibodies. Our lab developed VSG-seq, a targeted RNA-seq method, to study VSG gene expression in T. brucei. Studies using VSG-seq to characterize antigenic variation in a mouse model have revealed marked diversity in VSG expression within parasite populations, but this finding has not yet been validated in a natural human infection. Here, we used VSG-seq to analyze VSGs expressed in the blood of 12 patients infected with T. b. gambiense. The number of VSGs identified per patient ranged from one to 14; notably, two VSGs were shared by more than one patient. Analysis of expressed VSG N-terminal domain types revealed that 82% of expressed VSGs encoded a type B N-terminus, a bias not seen in datasets from other T. brucei subspecies. C-terminal types in T. b. gambiense infection were also restricted. These results demonstrate a bias either in the underlying VSG repertoire of T. b. gambiense or in this subspecies expression of VSGs during infection. This work demonstrates the feasibility of using VSG-seq to study antigenic variation in human infections and highlights the importance of understanding VSG repertoires in the field.

Author SummaryHuman African Trypanosomiasis is a neglected tropical disease primarily caused by the extracellular parasite Trypanosoma brucei gambiense. To avoid elimination by the host, these parasites repeatedly replace their dense surface coat of Variant Surface Glycoprotein (VSG). Despite the important role of VSGs in prolonging infection, VSG expression during natural human infections is poorly understood. A better understanding of natural VSG gene expression dynamics can clarify the mechanisms that T. brucei uses to alter its VSG coat and improve how trypanosomiasis is diagnosed in humans. We analyzed the expressed VSGs detected in the blood of patients with trypanosomiasis. Our findings indicate that parasites in both natural and experimental infections express diverse sets of VSGs. Juvenile Gouldian Finches (Erythrura gouldiae) form sibling sub-groups during social integration The formation of social relationships in complex groups is critical in shaping patterns of social organization and behavioral development. In many birds, young individuals remain dependent on their parents for extended periods but must abruptly transition to navigating interactions in the wider group after independence. While lack of social relationships during this period is detrimental in the development of later social skills, little is known about the social relationships juveniles form after independence in many bird species. In this study we describe patterns of social interactions in juvenile Gouldian Finches (Erythrura gouldiae) during transitions from family groups into flocks of unfamiliar individuals. A group of 20 juveniles from 4 families were introduced into two flocks. After introductions juveniles showed a gradient of approach rates with most approaches directed towards siblings, followed by juvenile peers, adult females, and lastly adult males. Significant preferences for siblings resulted in the emergence of sibling subgroups within the larger social network. This active self-assortment of siblings suggests that sibling sub-groups are an important bridge linking social connections within the family and the wider group. Such findings suggests that sibling relationships have a critical role in the socialization after independence, as well as structuring the social organization of Gouldian finch flocks. Tackling polypharmacology of kinase inhibitors and selecting chemical kinase probes by transcription factor activity profiling. Protein kinase inhibitors (PKI) evolved as promising drugs for multiple diseases. However, PKIs promiscuity causes polypharmacological effects challenging therapeutic development. We show that the polypharmacology of PKIs can be inferred from their impact on transcription factors (TF) lying at the apexes of signaling pathways. Using a high-content reporter assay, we evaluated TF activity profiles (TFAP) in cells treated with Akt, CDK, Aurora, Raf, MEK, and ERK inhibitors. For each kinase family, we found kinase-specific consensus TFAP signatures signifying the on-target PKI activity. Remarkably, proximal kinases had the highest similarity of their PKIs consensus signatures. Furthermore, we show that individual PKIs exhibited the consensus, "on-target signatures" within certain "specificity windows" and distinct "off-target signatures" at other concentrations. We also show that the TFAP signatures permit pinpointing the off-target PKI activity at irrelevant kinases and non-kinase effectors. Therefore, the TFAP approach provides clear-cut quantitative metrics for assessing PKI activity and concentration ranges where the on-target activity dominates cell response. By illuminating PKI biology invisible to target-based techniques, this effect-based approach streamlines the selection of kinase chemical probes and PKI drug prioritization. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. Phylogenomic discordance suggests polytomies along the backbone of the large genus Solanum Premise of the studyEvolutionary studies require solid phylogenetic frameworks, but increased volumes of phylogenomic data have revealed incongruent topologies among gene trees in many organisms both between and within genomes. Some of these incongruences indicate polytomies that may remain impossible to resolve. Here we investigate the degree of gene-tree discordance in Solanum, one of the largest flowering plant genera that includes the cultivated potato, tomato, and eggplant, as well as 24 minor crop plants.

MethodsA densely sampled species-level phylogeny of Solanum is built using unpublished and publicly available Sanger sequences comprising 60% of all accepted species (742 spp.) and nine regions (ITS, waxy, and seven plastid markers). The robustness of this topology is tested by examining a full plastome dataset with 140 species and a nuclear target-capture dataset with 39 species of Solanum (Angiosperms353 probe set).

Key resultsWhile the taxonomic framework of Solanum remained stable, gene tree conflicts and discordance between phylogenetic trees generated from the target-capture and plastome datasets were observed. The latter correspond to regions with short internodal branches, and network analysis and polytomy tests suggest the backbone is composed of three polytomies found at different evolutionary depths. The strongest area of discordance, near the crown node of Solanum, could potentially represent a hard polytomy.

ConclusionsWe argue that incomplete lineage sorting due to rapid diversification is the most likely cause for these polytomies, and that embracing the uncertainty that underlies them is crucial to understand the evolution of large and rapidly radiating lineages. 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.