PRDM9 losses in vertebrates are coupled to those of paralogs ZCWPW1 and ZCWPW2 In most mammals and likely throughout vertebrates, the gene PRDM9 specifies the locations of meiotic double strand breaks; in mice and humans at least, it also aids in their repair. For both roles, many of the molecular partners remain unknown. Here, we take a phylogenetic approach to identify genes that may be interacting with PRDM9, by leveraging the fact that PRDM9 arose before the origin of vertebrates, but was lost many times, either partially or entirely--and with it, its role in recombination. As a first step, we characterize PRDM9 domain composition across 446 vertebrate species, inferring at least thirteen independent losses. We then use the interdigitation of PRDM9 orthologs across vertebrates to test whether it co-evolved with any of 241 candidate genes co-expressed with PRDM9 in mice or associated with recombination phenotypes in mammals. Accounting for the phylogenetic relationship among species, we find two genes whose presence and absence is unexpectedly coincident with that of PRDM9: ZCWPW1, which was recently shown to facilitate double strand break repair, and its paralog ZCWPW2, as well as more tentative evidence for TEX15 and FBXO47. ZCWPW2 is expected to be recruited to sites of PRDM9 binding; its tight coevolution with PRDM9 across vertebrates suggests that it is a key interactor, with a role either in recruiting the recombination machinery or in double strand break repair.

Author SummaryOur understanding of meiotic recombination in mammals has seen great progress over the past 15 years, spurred in part by the convergence of lines of evidence from molecular biology, statistical genetics and evolutionary biology. We now know that in most mammals and likely in many vertebrates, the gene PRDM9 specifies the location of meiotic double strand breaks and that in mice and humans at least, it also aids in their repair. For both roles, however, many of the molecular partners remain unknown. To search for these, we take a phylogenetic approach, leveraging the fact that the complete PRDM9 has been lost at least thirteen times in vertebrates and thus that its presence is interdigitated across species. By this approach, we identify two genes whose presence or absence across vertebrates is coupled to the presence or absence of PRDM9, ZCWPW1 and ZCWPW2, as well as two genes for which the evidence is weaker, TEX15 and FBXO47. ZCWPW1 was recently shown to be recruited to sites of PRDM9 binding and to aid in the repair of double strand breaks. ZCWPW2 is likely recruited to sites of PRDM9 binding as well; its tight coevolution with PRDM9 across vertebrates suggests that it plays an important role either in double strand break formation, potentially as the missing link that recruits the recombination machinery to sites of PRDM9 binding, or in double strand break repair.
Factors behind intraspecific competition during post-glacial recolonization: a Norway spruce case study O_LITo a large extent, the current distribution and population structure of many species were shaped by the cycles of isolation in glacial refugia and subsequent population expansions. Isolation in, and post-glacial expansion through, heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a complex interplay between evolutionary forces.
C_LIO_LIWe investigated population differentiation and adaptations of Norway spruce for juvenile growth, diameter of the trunk at breast height and its components: wood density and tracheid traits. Data from 4,461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of QST vs FST, and to carry out multiple Genome-Wide Association Studies.
C_LIO_LIWe show that (i) the macroscopic resultant trait (trunk diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection, and (ii) adaptation predated the Last Glaciation Maximum.
C_LIO_LIAltogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly pre-adapted mainly through growth-related traits.
C_LI
Uncovering non-random binary patterns within sequences of intrinsically disordered proteins Sequence-ensemble relationships of intrinsically disordered proteins (IDPs) are governed by binary patterns such as the linear clustering or mixing of specific residues or residue types with respect to one another. To enable the discovery of potentially important, shared patterns across sequence families, we describe a computational method referred to as NARDINI for Non-random Arrangement of Residues in Disordered Regions Inferred using Numerical Intermixing. This work was partially motivated by the observation that parameters that are currently in use for describing different binary patterns are not interoperable across IDPs of different amino acid compositions and lengths. In NARDINI, we generate an ensemble of scrambled sequences to set up a composition-specific null model for the patterning parameters of interest. We then compute a series of pattern-specific z-scores to quantify how each pattern deviates from a null model for the IDP of interest. The z-scores help in identifying putative non-random linear sequence patterns within an IDP. We demonstrate the use of NARDINI derived z-scores by identifying sequence patterns in three well-studied IDP systems. We also demonstrate how NARDINI can be deployed to study archetypal IDPs across homologs and orthologs. Overall, NARDINI is likely to aid in designing novel IDPs with a view toward engineering new sequence-function relationships or uncovering cryptic ones. We further propose that the z-scores introduced here are likely to be useful for theoretical and computational descriptions of sequence-ensemble relationships across IDPs of different compositions and lengths. Hypertrophic Chondrocytes Serve as a Reservoir for Marrow Associated Skeletal Stem and Progenitor Cells, Osteoblasts, and Adipocytes During Skeletal Development Hypertrophic chondrocytes give rise to osteoblasts during skeletal development; however, the process by which these non-mitotic cells make this transition is not well understood. Prior studies have also suggested that skeletal stem and progenitor cells (SSPCs) localize to the surrounding periosteum and serve as a major source of marrow associated SSPCs, osteoblasts, osteocytes, and adipocytes during skeletal development. To further understand the cell transition process by which hypertrophic chondrocytes contribute to osteoblasts or other marrow associated cells, we utilized inducible and constitutive hypertrophic chondrocyte lineage tracing and reporter mouse models (Col10a1CreERT2; R26-tdTomatof/+ and Col10a1Cre; R26-tdTomatof/+) in combination with a PDGFRa-H2B-GFP transgenic line, single cell RNA-sequencing, bulk RNA-sequencing, immunofluorescence staining, and cell transplantation assays. Our data demonstrate that hypertrophic chondrocytes undergo a process of dedifferentiation to generate marrow associated SSPCs that serve as a primary source of osteoblasts during skeletal development. These hypertrophic chondrocyte derived SSPCs commit to a CXCL12-abundant reticular (CAR) cell phenotype during skeletal development and demonstrate unique abilities to recruit vasculature and promote bone marrow establishment, while also contributing to the adipogenic lineage. 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. Radiation-reprogrammed glioma stem cells generate vascular-like cells to build a trophic niche driving tumor recurrence Treatment-refractory glioma stem and tumor cells exhibit phenotypic plasticity driving recurrence, but the underlying molecular mechanisms remain to be elucidated. Here, we employed single-cell and whole transcriptomic analyses to discover that radiation induces a dynamic shift in functional states of glioma cells allowing for acquisition of vascular endothelial-like and pericyte-like cell phenotypes. These vascular-like cells provide a trophic niche to promote proliferation of irradiated glioma cells, and their selective depletion results in reduced tumor growth post-treatment in vivo. Mechanistically, the acquisition of vascular-like phenotype is driven by increased chromatin accessibility and H3K27 acetylation in specific vascular gene regions post-treatment. Blocking P300 histone acetyltransferase activity reverses the epigenetic changes induced by radiation, and inhibits the phenotypic transition and tumor growth. Our findings highlight an important role for P300 histone acetyltransferase in treatment-induced plasticity and opens a new therapeutic avenue for preventing glioma recurrence.

SignificanceOur study demonstrates that radiation therapy promotes glioma resistance by inducing vascular-like phenotypes in GSC that, in turn, aid in proliferation of the remaining tumor cells. This phenotype switch is mediated by P300 HAT, and inhibition of this enzyme is a potential therapeutic target for preventing glioma recurrence following radiation. 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. Spatio-temporal, optogenetic control of gene expression in organoids Organoids derived from stem cells become increasingly important to study human development and to model disease. However, methods are needed to control and study spatio-temporal patterns of gene expression in organoids. To this aim, we combined optogenetics and gene perturbation technologies to activate or knock-down RNA of target genes, at single-cell resolution and in programmable spatio-temporal patterns. To illustrate the usefulness of our approach, we locally activated Sonic Hedgehog (SHH) signaling in an organoid model for human neurodevelopment. High-resolution spatial transcriptomic and single-cell analyses showed that this local induction was sufficient to generate stereotypically patterned organoids in three dimensions and revealed new insights into SHHs contribution to gene regulation in neurodevelopment.

With this study, we propose optogenetic perturbations in combination with spatial transcriptomics as a powerful technology to reprogram and study cell fates and tissue patterning in organoids. Development of a Novel PROTAC using the Nucleic Acid Aptamer as a Targeting Ligand for Tumor Selective Degradation of Nucleolin PROteolysis TArgeting Chimeras (PROTACs) induce targeted protein degradation by hijacking the intracellular ubiquitin proteasome system, thus emerging as a new strategy for drug development. However, most PROTACs generated lack cell-type selectivity and have poorly soluble in water. To address this drawback, we developed a novel PROTAC ZL216 using aptamer AS1411 as a targeting ligand of nucleolin to conjugate with a small molecule ligand of E3 ligase VHL, which shows high aqueous solubility and serum stability. Based on the differential expression of nucleolin on the cell surface, ZL216 could bind to and internalize into breast cancer cells, but not normal breast cells. Furthermore, we revealed that ZL216 promoted the formation of a nucleolin-ZL216-VHL ternary complex in breast cancer cells and potently induced nucleolin degradation in vitro and in vivo. As a result, ZL216 inhibited the proliferation and migration of breast cancer cells. These studies demonstrate that in addition to peptides and small molecule compounds, nuclei acid aptamers can also be used to generate PROTACs, which broadens the toolbox constructing PROTACs and provides a promising strategy for development of tumor-selective PROTACs. Maf/ham1-like pyrophosphatases of non-canonical nucleotides are host-specific partners of viral RNA-dependent RNA polymerases Cassava brown streak disease (CBSD), dubbed the "Ebola of plants", is a serious threat to food security in Africa caused by two viruses of the family Potyviridae: cassava brown streak virus (CBSV) and Ugandan (U)CBSV. Intriguingly, U/CBSV, along with another member of this family and one secoviridae, are the only known RNA viruses encoding a protein of the Maf/ham1-like family, a group of widespread pyrophosphatase of non-canonical nucleotides (ITPase) expressed by all living organisms. Despite the socio-economic impact of CDSD, the relevance and role of this atypical viral factor has not been yet established. Here, using an infectious cDNA clone and reverse genetics, we demonstrate that UCBSV requires the ITPase activity for infectivity in cassava, but not in the model plant Nicotiana benthamiana. HPLC MS/MS experiments showed that, quite likely, this host-specific constraint is due to an unexpected high concentration of non-canonical nucleotides in cassava. Finally, protein analyses and experimental evolution of mutant viruses indicated that keeping a fraction of the yielded UCBSV ITPase covalently bound to the viral RNA-dependent RNA polymerase (RdRP) optimizes viral fitness, and this seems to be a feature shared by the other members of the Potyviridae family expressing Maf/ham1-like proteins. All in all, our work (i) reveals that the over-accumulation of non-canonical nucleotides in the host might have a key role in antiviral defense, and (ii) provides the first example of an RdRP-ITPase partnership, reinforcing the idea that RNA viruses are incredibly versatile at adaptation to different host setups. Naming of human diseases on the wrong side of history BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. In Mycobacterium abscessus, the stringent factor Rel regulates metabolism, but is not the only (p)ppGpp synthase The stringent response is a broadly conserved stress response system that exhibits functional variability across bacterial clades. Here, we characterize the role of the stringent factor Rel in the non-tuberculous mycobacterial pathogen, Mycobacterium abscessus (Mab). We found that deletion of rel does not ablate (p)ppGpp synthesis, and that rel does not provide a survival advantage in several stress conditions, or in antibiotic treatment. Transcriptional data show that RelMab is involved in regulating expression of anabolism and growth genes in stationary phase. However, it does not activate transcription of stress response or antibiotic resistance genes, and actually represses transcription of many antibiotic resistance genes. This work shows that there is an unannotated (p)ppGpp synthetase in Mab.

ImportanceIn this study, we examined the functional roles of the stringent factor Rel in Mycobacterium abscessus (Mab). In most species, stringent factors synthesize the alarmone (p)ppGpp, which globally alters transcription to promote growth arrest and survival under stress and in antibiotic treatment. Our work shows that in Mab, an emerging pathogen which is resistant to many antibiotics, the stringent factor Rel is not solely responsible for synthesizing (p)ppGpp. We find that RelMab downregulates many metabolic genes under stress, but does not upregulate stress response genes and does not promote antibiotic tolerance. This study implies that there is another critical but unannotated (p)ppGpp synthetase in Mab, and suggests that RelMab inhibitors are unlikely to sensitize Mab infections to antibiotic treatment. Genomic evidence of paternal genome elimination in the globular springtail Allacma fusca Paternal genome elimination (PGE) - a type of reproduction in which males inherit but fail to pass on their fathers genome - evolved independently in six to eight arthropod clades. Thousands of species, including several important for agriculture, reproduce via this mode of reproduction. While PGE is well established in some of the clades, the evidence in globular springtails (Symphypleona) remains elusive, even though they represent the oldest and most species rich clade putatively reproducing via PGE. We sequenced genomic DNA from whole bodies of Allacma fusca males with high fractions (>27.5%) of sperm to conclusively confirm that all the sperm carry one parental haplotype only. Although it is suggestive that the single haplotype present in sperm is maternally inherited, definitive genetic proof of the parent of origin is still needed. The genomic approach we developed allows for detection of genotypic differences between germline and soma in all species with sufficiently high fraction of germline in their bodies. This opens new opportunities for scans of reproductive modes in small organisms. Differed cell division angle, position of cell proliferative area, and localization of ANGUSTIFOLIA3 in lateral organs Leaf meristem is a cell proliferative zone present in the lateral organ primordia. In this study, we investigated how the proliferative zone affects the final morphology of the lateral organs. We examined how cell proliferative zones differ in the primordia of planar floral organs and polar auxin transport inhibitor (PATI)-treated leaves from normal foliage leaf primordia of Arabidopsis thaliana with a focus on the spatial accumulation pattern of ANGUSTIFOLIA3 (AN3), a key element for leaf meristem positioning. We found that organ shape changes by PATI treatment were correlated to the angle of the cell division plane relative to the leaf primordia axis in the leaf meristem (cell division angle), but not with leaf-meristem positioning, size of the leaf meristem, or the localization pattern of AN3 protein. In contrast, different shapes between sepals and petals compared with foliage leaves were associated with both altered meristem position associated with altered AN3 expression patterns and different distributions of cell division angles. These results suggest that lateral organ shapes are regulated via two aspects: position of meristem and cell division angles

Summary statementDifferent lateral organs with different morphology possess different properties of meristems; cell division angles, position of cell proliferative area and AN3 localization patterns. Integrative small and long RNA-omics analysis of human healing and non-healing wounds discovers cooperating microRNAs as therapeutic targets MicroRNAs (miR) are important posttranscriptional regulators and exhibit a high potential to be utilized in diagnosis and therapy. However, our insufficient knowledge of the miR-mediated gene regulation in human skin wound healing severely hinders the identification of clinically relevant miRs. Here, we performed paired small RNA and long RNA sequencing in human tissue samples, including matched skin and acute wounds collected at each healing stage and chronic non-healing venous ulcers (VU). With integrative small and long RNA-omics analysis, we developed a compendium (https://www.xulandenlab.com/humanwounds-mirna-mrna), which will be an open, comprehensive resource to broadly aid wound healing research. With this first clinical, wound-centric resource of miRs and mRNAs, we identified 17 pathologically relevant miRs that exhibited abnormal VU expression and displayed their targets enriched explicitly in the VU gene signature. Intermeshing regulatory networks controlled by these miRs revealed their high cooperativity in contributing to chronic wound pathology characterized by persistent inflammation and proliferative phase initiation failure. Furthermore, we demonstrated that miR-34a, miR-424, and miR-516, upregulated in VU, cooperatively suppressed keratinocyte growth while promoting inflammatory response. Collectively, our study opens the possibility of developing innovative wound treatment that targets pathologically relevant cooperating miRs to attain higher therapeutic efficacy and specificity. How grower decision-making is modelled impacts the effect of plant disease control. While the spread of plant disease depends strongly on biological factors driving transmission, it also has a human dimension. Disease control depends on decisions made by individual growers, who are in turn influenced by a broad range of factors. Despite this, human behaviour has rarely been included in plant epidemic models. Considering Cassava Brown Streak Disease, we model how the perceived increase in profit due to disease management influences participation in clean seed systems (CSS). Our models are rooted in game theory, with growers making strategic decisions based on the expected profitability of different control strategies. We find that both the information used by growers to assess profitability and the perception of economic and epidemiological parameters influence long-term participation in the CSS. Over-estimation of infection risk leads to lower participation in the CSS, as growers perceive that paying for the CSS will be futile. Additionally, even though good disease management can be achieved through the implementation of CSS, and a scenario where all controllers use the CSS is achievable when growers base their decision on the average of their entire strategy, CBSD is rarely eliminated from the system. These results are robust to stochastic and spatial effects. Our work highlights the importance of including human behaviour in plant disease models, but also the significance of how that behaviour is included. The RAG1 Ubiquitin Ligase Domain Enhances Assembly and Selection of T Cell Receptor Genes to Restrain the Autoimmune Hazard of Generating T Cell Receptor Diversity RAG1/RAG2 (RAG) endonuclease-mediated assembly of diverse lymphocyte antigen receptor genes by V(D)J recombination is critical for the development and immune function of T and B cells. However, this process creates highly self-reactive cells that must be properly selected to suppress autoimmunity. The RAG1 protein contains a ubiquitin ligase domain that stabilizes RAG1 and stimulates RAG endonuclease activity in vitro. We report that mice with a mutation that inactivates the RAG1 ubiquitin ligase in vitro exhibit modestly reduced thymic cellularity, decreased assembly and altered repertoires of T cell receptor (TCR) {beta} and  genes in thymocytes, and impaired thymocyte developmental transitions that require the assembly of TCR{beta} or  genes and signaling by their proteins. These RAG1 mutant mice also exhibit less efficient positive selection and superantigen-mediated negative selection of conventional {beta} T cells, 2) impaired differentiation of iNKT lineage {beta} T cells, and 3) CD4+ {beta} T cells with elevated autoimmune potential. Our findings demonstrate that the RAG1 ubiquitin ligase domain functions in vivo to stimulate the assembly and selection of TCR{beta} and TCR genes, thereby establishing replete diversity of {beta} TCRs and {beta} T cell lineages while restraining the inherent autoimmune hazard of generating diverse antigen specificities. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. Exploratory analysis of multiple traits co-adaptations in the population history During the history of range expansion, the populations encounter with variety of environments. They respond to the local environments by modifying the mutually interacting traits. Therefore, to understand the whole life history of the populations, it is ideal to capture the history of their range expansion with reference to the series of surrounding environments and to infer the coadaptation of the multiple traits. Toward this end, we provide an exploratory analysis based on the features of populations: site frequency spectra of populations, population-specific FST, association between genes and environments, positive selections on traits mapped on the admixture graph, and GWAS results. Correspondence analysis of genes, environments, and traits provides a birds-eye view of the history of population differentiation and range expansion and various types of environmental selections at the times. Principal component analysis of the estimated trait-specific polygenic adaptations mapped on the admixture graph enables to understand the coadaptation of multiple traits. The potential usefulness was confirmed by analyzing a public dataset of wild poplar in northwestern America. In response to the northern cold temperature and longer daylength, the populations increased the photosynthetic activity and nutrient use efficiency at the expense of the risk of pathogen invasion, and in response to warm temperature, they increased the growth. At higher altitude, they shifted the maximum activity to earlier period in spring to reduce the activity in dry summer. The R codes for our representation method and simulations of population colonization used in this study are available as supplementary script. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Yolk sac-derived bipotent endothelial-macrophage progenitor cells are a self-renewing tissue-resident source of inflammation and neovascularization postnatally Converging evidence indicates that extra-embryonic yolk sac is the source of both macrophages and endothelial cells in adult mouse tissues. Prevailing views are that these yolk sac-derived cells are maintained after birth by proliferative self-renewal in their differentiated states. Here we identify clonogenic, self-renewing endothelial-macrophage (EndoMac) progenitor cells in postnatal mouse aorta, heart and lung, that are independent of definitive hematopoiesis and derive from a CX3CR1+ and CSF1R+ yolk sac source. These bipotent progenitors are highly proliferative and vasculogenic, contributing to adventitial neovascularization in the aortic wall and forming perfused blood vessels after adoptive transfer into ischemic tissue. We establish a regulatory role for angiotensin II, which enhances their clonogenic, self-renewal and differentiation properties. Our findings demonstrate that tissue-resident EndoMac progenitors drive local inflammatory and vasculogenic responses by contributing to the renewal and expansion of yolk sac-derived macrophages and endothelial cells postnatally. 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. Reconsider phase reconstruction in signals with dynamic periodicity from the modern signal processing perspective Phase is the most fundamental physical quantity when we study an oscillatory time series. There are many tools aiming to estimate phase, most of them are developed based on the analytic function model. Unfortunately, this approach might not be suitable for modern signals with intrinsic nonstartionary structure, including multiple oscillatory components, each with time-varying frequency, amplitude, and non-sinusoidal oscillation, e.g., biomedical signals. Specifically, due to the lack of consensus of model and algorithm, phases estimated from signals simultaneously recorded from different sensors for the same physiological system from the same subject might be different. This fact might challenge reproducibility, communication, and scientific interpretation and thus we need a standardized approach with theoretical support over a unified model. In this paper, after summarizing existing models for phase and discussing the main challenge caused by the above-mentioned intrinsic nonstartionary structure, we introduce the adaptive non-harmonic model (ANHM), provide a definition of phase called fundamental phase, which is a vector-valued function describing the dynamics of all oscillatory components in the signal, and suggest a time-varying bandpass filter (tvBPF) scheme based on time-frequency analysis tools to estimate the fundamental phase. The proposed approach is validated with a simulated database and a real-world database with experts labels, and it is applied to two real-world databases, each of which has biomedical signals recorded from different sensors, to show how to standardize the definition of phase in the real-world experimental environment. Specifically, we report that the phase describing a physiological system, if properly modeled and extracted, is immune to the selected sensor for that system, while other approaches might fail. In conclusion, the proposed approach resolves the above-mentioned scientific challenge. We expect its scientific impact on a broad range of applications. 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. Robust and efficient online auditory psychophysics Most human auditory psychophysics research has historically been conducted in carefully controlled environments with calibrated audio-equipment, and over potentially hours of repetitive testing with expert listeners. Here, we operationally define such conditions as having high  auditory hygiene. Hereof, conducting auditory psychophysical paradigms online presents a serious challenge, in that results may hinge on absolute sound presentation level, reliably estimated perceptual thresholds, and sustained motivation and attention. We introduce a set of procedures that address these challenges and facilitate auditory hygiene for online auditory psychophysics. First, we establish a simple means of setting sound presentation levels. Across a set of four level-setting conditions conducted in person, we demonstrate the stability and robustness of this amplitude setting procedure in open air and controlled settings. Second, we test participants tone-in-noise thresholds using widely adopted online experiment platforms and demonstrate that reliable threshold estimates can be derived online in approximately one minute of testing. Third, using these amplitude and threshold setting procedures to establish participant-specific stimulus conditions, we show that an online implementation of the classic probe-signal paradigm can be used to demonstrate frequency-selective attention on an individual participant basis, using a third of the trials used in recent in-lab experiments. Finally, we show how threshold and attentional measures relate to well-validated assays of online participants in-task motivation, fatigue, and confidence. This demonstrates the promise of online auditory psychophysics for asking new auditory neuroscience questions quickly, efficiently, and with more diverse samples. Code for the tests is publicly available through Pavlovia and Gorilla.