DCIS genomic signatures define biology and correlate with clinical outcome: a Human Tumor Atlas Network (HTAN) analysis of TBCRC 038 and RAHBT cohorts Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We have performed the first multiscale, integrated profiling of DCIS with clinical outcomes by analyzing 677 DCIS samples from 481 patients with 7.1 years median follow-up from the Translational Breast Cancer Research Consortium (TBCRC) 038 study and the Resource of Archival Breast Tissue (RAHBT) cohorts. We made observations on DNA, RNA, and protein expression, and generated a de novo clustering scheme for DCIS that represents a fundamental transcriptomic organization at this early stage of breast neoplasia. Distinct stromal expression patterns and immune cell compositions were identified. We found RNA expression patterns that correlate with later events. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/448585v2_ufig1.gif" ALT="Figure 1">
View larger version (44K):
org.highwire.dtl.DTLVardef@1755117org.highwire.dtl.DTLVardef@11ebdc0org.highwire.dtl.DTLVardef@13ba799org.highwire.dtl.DTLVardef@1784502_HPS_FORMAT_FIGEXP  M_FIG C_FIG HIGHLIGHTSO_LINew transcriptomic classification solution reveals 3 major subgroups in DCIS.
C_LIO_LIFour stroma-specific signatures identified.
C_LIO_LIutcome analysis identifies pathways involved in DCIS progression.
C_LIO_LICNAs characterize high risk of distant relapse IBC subtypes observed in DCIS.
C_LI
A MALDI-TOF assay identifies nilotinib as an inhibitor of inflammation in acute myeloid leukaemia Inflammatory responses are important in cancer, particularly in the context of monocyte-rich aggressive myeloid neoplasm. We developed a label-free cellular phenotypic drug discovery assay to identify anti-inflammatory drugs in human monocytes derived from acute myeloid leukaemia (AML), by tracking several biological features ionizing from only 2,500 cells using matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry. A proof-of-concept screen showed that the BCR-ABL inhibitor nilotinib, but not the structurally similar imatinib, blocks inflammatory responses. In order to identify the cellular (off-)targets of nilotinib, we performed thermal proteome profiling (TPP). Unlike imatinib, nilotinib and other later generation BCR-ABL inhibitors inhibit the p38-MK2/3 signalling axis which suppressed the expression of inflammatory cytokines, cell adhesion and innate immunity markers in activated human monocytes derived from AML. Thus, our study provides a tool for the discovery of new anti-inflammatory drugs, which could contribute to the treatment of inflammation in myeloid neoplasms and other diseases.

Key PointsLabel-free cell-based assay identifies new anti-inflammatory drugs using MALDI-TOF MS. Nilotinib reduces inflammation by inhibition of MAPK14-MK2/3 signalling axis in AML.
DENTIST - using long reads for closing assembly gaps at high accuracy BackgroundLong sequencing reads allow increasing contiguity and completeness of fragmented, short-read based genome assemblies by closing assembly gaps, ideally at high accuracy. While several gap closing methods have been developed, these methods often close an assembly gap with sequence that does not accurately represent the true sequence.

FindingsHere, we present DENTIST, a sensitive, highly accurate and automated pipeline method to close gaps in short-read assemblies with long error-prone reads. DENTIST comprehensively determines repetitive assembly regions to identify reliable and unambiguous alignments of long reads to the correct loci, integrates a consensus sequence computation step to obtain a high base accuracy for the inserted sequence, and validates the accuracy of closed gaps. Unlike previous benchmarks, we generated test assemblies that have gaps at the exact positions where real short-read assemblies have gaps. Generating such realistic benchmarks for Drosophila (134 Mb genome), Arabidopsis (119 Mb), hummingbird (1 Gb) and human (3 Gb) and using simulated or real PacBio continuous long reads, we show that DENTIST consistently achieves a substantially higher accuracy compared to previous methods, while having a similar sensitivity.

ConclusionDENTIST provides an accurate approach to improve the contiguity and completeness of fragmented assemblies with long reads. DENTISTs source code including a Snakemake workflow, conda package and Docker container is available at https://github.com/a-ludi/dentist. All test assemblies as a resource for future benchmarking are at https://bds.mpi-cbg.de/hillerlab/DENTIST/. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. 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. EMBED: Essential Microbiome Dynamics, a dimensionality reduction approach for longitudinal microbiome studies Dimensionality reduction can offer unique insights into high dimensional microbiome dynamics by leveraging collective abundance fluctuations of multiple bacteria driven by similar ecological perturbations. However, methods providing lower-dimensional representations of microbiome dynamics both at the community and individual taxa level are not currently available. To that end, we present EMBED: Essential MicroBiomE Dynamics, a probabilistic non-linear tensor factorization approach. Similar to normal mode analysis in structural biophysics, EMBED infers ecological normal modes (ECNs), which represent the unique orthogonal modes capturing the collective behavior of microbial communities. A very small number of ECNs can accurately approximate microbiome dynamics across multiple data sets. Inferred ECNs reflect specific ecological behaviors, providing natural templates along which the dynamics of individual bacteria may be partitioned. Moreover, the multi-subject treatment in EMBED systematically identifies subject-specific and universal abundance dynamics that are not detected by traditional approaches. Collectively, these results highlight the utility of EMBED as a versatile dimensionality reduction tool for studies of microbiome dynamics. Induced Cognitive Impairments Reversed by Grafts of Neural Precursors: a Longitudinal Study in a Macaque Model of Parkinson's Disease Parkinsons disease (PD) evolves over an extended and variable period in humans; several years prior to the onset of classical motor symptoms, cognitive deficits as well as sleep and biological rhythm disorders develop and worsen with disease progression, significantly impacting the quality of life of patients. The gold standard MPTP macaque model of PD recapitulates the progression of motor and non-motor symptoms over contracted periods of time.

Here, this multidisciplinary and multiparametric study follows, in five animals, the steady progression of motor and non-motor symptoms and describes their reversal following bilateral grafts of neural precursors in diverse functional domains of the basal ganglia.

Results show unprecedented recovery from cognitive symptoms in addition to a strong clinical motor recuperation. Both motor and cognitive recovery and partial circadian rhythm recovery correlate with the degree of graft integration into the host environment as well as with in-vivo levels of striatal dopaminergic innervation and function.

Given inter-individuality of disease progression and recovery the present study underlines the importance of longitudinal multidisciplinary assessments in view of clinical translation and provides empirical evidence that integration of neural precursors following transplantation efficiently restores function at multiple levels in parkinsonian non-human primates.

One Sentence SummaryEmpirical evidence that cell therapy efficiently reverts cognitive and clinical motor symptoms in the non-human primate model of Parkinsons disease. Integrative transcriptomics reveals ectopic lipid homeostasis mechanisms in non-endocrine cells of the teleost pituitary Directing both organismal homeostasis and physiological adaptation, the pituitary is a key endocrine gland in all vertebrates. It communicates the needs of the organism to different organs by secreting hormones into the bloodstream. Here, we have used the model fish medaka to investigate the developmental dynamics in the pituitary using a comprehensive RNA-seq time series. By linking developmental expression trends to single-cell RNA-seq profiles, we show how the transcriptional activities of cell types change during sexual maturation. One of the most prominent changes is the decline of the non-endocrine folliculo-stellate cell populations, and especially of rare cells expressing genes encoding secreted lipid transport proteins. As these genes are typically associated with the liver, this reveals the existence of unexpected connections between endocrine communication, lipid homeostasis, and sexual maturation. Structural insights into p300 regulation and acetylation-dependent genome organisation Histone modifications are deposited by chromatin modifying enzymes and read out by proteins that recognize the modified state. BRD4-NUT is an oncogenic fusion protein of the acetyl lysine reader BRD4 that binds to the acetylase p300 and enables formation of long-range intra- and interchromosomal interactions. We here examine how acetylation reading and writing enable formation of such interactions. We show that NUT contains an acidic transcriptional activation domain that binds to the TAZ2 domain of p300. We use NMR to investigate the structure of the complex and found that the TAZ2 domain has an autoinhibitory role for p300. NUT-TAZ2 interaction or mutations found in cancer that interfere with autoinhibition by TAZ2 allosterically activate p300. p300 activation results in a self-organizing, acetylation-dependent feed-forward reaction that enables long-range interactions by bromodomain multivalent acetyl-lysine binding. We discuss the implications for chromatin organisation, gene regulation and dysregulation in disease. Contribution Evenness: A functional redundancy metric sensitive to trait stability in microbial communities The concept of functional redundancy has received considerable attention in both the macroecology and microbial ecology literature. As a result, multiple metrics of functional redundancy have been proposed. These vary in how they weight trait levels, species abundance, functional richness, and species richness. Here we present a new functional redundancy metric tailored for community-aggregated traits, which are traits that are quantified at the community level and can be quantitatively partitioned among species. We call this metric Contribution Evenness (CE) because it measures how evenly species contribute to a community-aggregated trait. As CE is an evenness measurement, it ranges from 0 and 1, where 0 corresponds to a single species contributing to a community-aggregated trait and 1 corresponds to all species contributing equally. Using in silico simulations of species extinctions, we demonstrate that CE reflects the stability of an ecosystem function to species extinction, a hypothesized ecological consequence of functional redundancy. As a positive control and to illustrate how CE can be used with sequence data, we analyzed the functional redundancy of eight nitrogen-transforming pathways using 2,631 metagenome-assembled genomes from 47 marine sites. CE for marine nitrogen cycle marker genes was consistent with our qualitative understanding of which nitrogen pathways are most functionally redundant in the ocean. We found that, on average, the NH4+ assimilation pathway was the most functionally redundant (0.44 {+/-} 0.08) while dissimilatory nitrate reduction was the least redundant (0.005 {+/-} 0.005). As demonstrated here, CE provides a promising framework for measuring trait stability in microbiomes. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Nutrient metabolism regulates insulin granule formation in the pancreatic islet β-cell via ER redox homeostasis Defects in the pancreatic {beta}-cells secretion system are well-described in Type 2 diabetes (T2D) and include impaired proinsulin processing and a deficit in mature insulin-containing secretory granules; however, the cellular mechanisms underlying these defects and the contribution of hyperglycemia to this process remain poorly understood. Here, we used an in situ fluorescent pulse-chase strategy and proximity labeling-based quantitative proteomics analysis to study proinsulin trafficking and demonstrate a direct link to glucose metabolism via the production of redox intermediates that facilitate proinsulin export from the ER. We show that ER export of proinsulin is delayed in T2D models resulting in decreased insulin granule formation and further demonstrate this process can be regulated by NADPH and reducing equivalent availability. Together, these data highlight a critical role for nutrient metabolism and mitochondrial dysfunction in the maladaptive remodeling of the {beta}-cells secretory pathway in the decline of {beta}-cell function in T2D. Single cell profiling of functionally cured Chronic Hepatitis B patients reveals the emergence of activated innate and an altered adaptive immune response in the intra-hepatic environment Hepatitis B surface antigen (HBsAg) loss or functional cure (FC), is considered the desirable therapeutic outcome for chronic hepatitis B (CHB) patients. However, the immuno-pathological biomarkers and underlying mechanisms remain unclear. Here we present a comprehensive single cell-transcriptomic atlas together with immune-phenotyping of disease-associated cell states (DACS) isolated from intra-hepatic tissue and matched PBMCs of either CHB or FC patients. We find that the intra-hepatic environment displays specific cell identities and molecular signatures that are distinct from PBMCs. FC is associated with emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes (CD4-CTLs), and an activated innate response represented by liver-resident natural killer (LR-NK) cells. Overall, these findings provide novel insights into immuno-pathological cell states associated with FC that could serve as prognostic biomarkers. `rbioacc`: an R-package to analyse toxicokinetic data O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI Dynamic modulation of enhancer responsiveness by core promoter elements in living Drosophila embryo Regulatory interactions between enhancers and core promoters are fundamental for the temporal and spatial specificity of gene expression in development. The central role of core promoters is to initiate productive transcription in response to enhancers activation cues. However, it has not been systematically assessed how individual core promoter elements affect the induction of transcriptional bursting by enhancers. Here, we provide evidence that each core promoter element dynamically alters "enhancer responsiveness" by changing functional parameters of transcriptional bursting in developing Drosophila embryos. Quantitative live imaging analysis revealed that the timing and the continuity of burst induction are common regulatory steps on which core promoter elements impact. We also show that core promoters use the upstream TATA box as an amplifier of transcriptional bursting upon interaction with enhancers. Genome editing analysis of the pair-rule gene fushi tarazu further revealed that the endogenous TATA box and the downstream core promoter element are both essential for the correct expression and the function of key developmental genes in early embryos. We suggest that core promoter elements serve as a key regulatory module in converting enhancer activity into transcription dynamics during animal development. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. The effect of recent habitat change on genetic diversity at putatively adaptive and neutral loci in Primula veris in semi-natural grasslands Recent habitat change in semi-natural grasslands due to a lack of management has been shown to affect the genetic diversity of grassland plants. However, it is unknown how a change in local environment affects genetic diversity at adaptive loci. We applied RADseq (restriction-site associated DNA sequencing) to extract > 3,000 SNPs across 568 individuals from 32 Estonian populations of Primula veris, a plant species common to semi-natural grasslands. We evaluated the effect of recent grassland overgrowth following management abandonment on the genetic diversity at putatively neutral and adaptive loci, which we distinguished by applying three methods, i.e., linear and categorical environmental association analyses, and an FST outlier test. For validation, we randomised the genotype to sample assignments. Effects of recent habitat change on genetic diversity differed between neutral and adaptive SNP sets. Genetic diversity assessed at putatively neutral loci was similar in open and overgrown habitats but showed a significant difference between these habitat types at putatively adaptive loci: overgrown (i.e. newly established) habitats exhibited higher genetic diversity at putatively adaptive loci than open (i.e. old) habitats, likely due to the exertion of novel selection pressures imposed by new habitat conditions. This increase in genetic diversity at putatively adaptive loci in the new environment points to currently ongoing selection processes where genetic adaptation to the old habitat is potentially lost through altered allele frequencies. Our study suggests that a recent change in local habitat conditions may not be reflected in neutral loci whereas putatively adaptive loci can inform about potential selection processes. Effects of the social environment on movement-integrated habitat selection O_LIMovement links the distribution of habitats with the social environment of animals using those habitats; yet integrating movement, habitat selection, and socioecology remains an opportunity for further study.
C_LIO_LIHere, our objective was to disentangle the roles of habitat selection and social association as drivers of collective movement in a gregarious ungulate. To accomplish this objective, we (1) assessed whether socially familiar individuals form discrete social communities and whether social communities have high spatial, but not necessarily temporal, overlap; and (2) we modelled the relationship between collective movement and selection of foraging habitats using socially informed integrated step selection analysis.
C_LIO_LIWe used social network analysis to assign individuals to social communities and determine short and long-term social preference among individuals. Using integrated step selection functions (iSSF), we then modelled the effect of social processes, i.e., nearest neighbour distance and social preference, and movement behaviour on patterns of habitat selection.
C_LIO_LIBased on assignment of individuals to social communities and home range overlap analyses, individuals assorted into discrete social communities, and these communities had high spatial overlap. By unifying social network analysis with iSSF, we identified movement-dependent social association, where individuals foraged with more familiar individuals, but moved collectively with any between foraging patches.
C_LIO_LIOur study demonstrates that social behaviour and space use are inter-related based on spatial overlap of social communities and movement-dependent habitat selection. Movement, habitat selection, and social behaviour are linked in theory. Here, we put these concepts into practice to demonstrate that movement is the glue connecting individual habitat selection to the social environment.
C_LI Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. 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. 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. A time irreversible model of nucleotide substitution for the coronavirus SARS-CoV-2 SARS-CoV-2 is the cause of the worldwide epidemic of severe acute respiratory syndrome. Evolutionary studies of the virus genome will provide a predictor of the fate of COVID-19 in the near future. Recent studies of the virus genomes have shown that C to U substitutions are overrepresented in the genome sequences of SARS-CoV-2. Traditional time-reversible substitution models cannot be applied to the evolution of SARS-CoV-2 sequences. Therefore, in this study, I propose a new time-irreversible model and a new method for estimating the nucleotide substitution rate of SARS-CoV-2. Computer simulations showed that that the new method gives good estimates. I applied the new method to estimate nucleotide substitution rates of SARS-CoV-2 sequences. The result suggests that the rate of C to U substitution of SARS-Cov-2 is ten times higher than other types of substitutions. Single-dose AAV-based vaccine induces a high level of neutralizing antibodies against SARS-CoV-2 in rhesus macaques Coronavirus disease 2019 (COVID-19), which is triggered by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, continues to threaten global public health. Developing a vaccine that only requires single immunization but provides long-term protection for the prevention and control of COVID-19 is important. Here, we developed an adeno-associated virus (AAV)-based vaccine expressing a stable receptor-binding domain (SRBD) protein. The vaccine requires only a single shot but provides effective neutralizing antibodies (NAbs) over 598 days in rhesus macaques (Macaca mulatta). Importantly, our results showed that the NAbs were kept in high level and long lasting against authentic wild-type SARS-CoV-2, Beta, Delta and Omicron variants using plaque reduction neutralization test. Of note, although we detected pre-existing AAV2/9 antibodies before immunization, the vaccine still induced high and effective NAbs against COVID-19 in rhesus macaques. AAV-SRBD immune serum also efficiently inhibited the binding of ACE2 with RBD in the SARS-CoV-2 B.1.1.7 (Alpha), B.1.351 (Beta), P.1/P.2 (Gamma), B.1.617.2 (Delta), B.1.617.1/3(Kappa), and C.37 (Lambda) variants. Thus, these data suggest that the vaccine has great potential to prevent the spread of SARS-CoV-2. A biologically accurate model of directional hearing in the parasitoid fly Ormia ochracea This manuscript was compiled on October 5, 2021 Although most binaural organisms localize sound sources using neurological structures to amplify the sounds they hear, some animals use mechanically coupled hearing organs instead. One of these animals, the parasitoid fly Ormia ochracea, has astoundingly accurate sound localization abilities and can locate objects in the azimuthal plane with a precision of 2{degrees}, equal to that of humans. This is accomplished despite an intertympanal distance of only 0.5 mm, which is less than 1/100th of the wavelength of the sound emitted by the crickets that it parasitizes. In 1995, Miles et al. developed a model of hearing mechanics in O. ochracea, which works well for incoming sound angles of less than {+/-}30{degrees}, but suffers from reduced accuracy (up to 60% error) at higher angles. Even with this limitation, it has served as the basis for multiple bio-inspired microphone designs for decades. Here, we present critical improvements to the classic O. ochracea hearing model based on information from 3D reconstructions of O. ochraceas tympana. The 3D images reveal that the tympanal organ has curved lateral faces in addition to the flat front-facing prosternal membranes represented in the Miles model. To mimic these faces, we incorporated spatially-varying spring and damper coefficients that respond asymmetrically to incident sound waves, making a new quasi-two-dimensional (q2D) model. The q2D model has high accuracy (average errors of less than 10%) for the entire range of incoming sound angles. This improved biomechanical hearing model can inform the development of new technologies and may help to play a key role in developing improved hearing aids.

Significance StatementThe ability to identify the location of sound sources is critical to organismal survival and for technologies that minimize unwanted background noise, such as directional microphones for hearing aids. Because of its exceptional auditory system, the parasitoid fly Ormia ochracea has served as an important model for binaural hearing and a source of bioinspiration for building tiny directional microphones with outsized sound localization abilities. Here, we performed 3D imaging of the flys tympanal organs and used the morphological information to improve the current model for hearing in O. ochracea. This model greatly expands the range of biological accuracy from {+/-}30{degrees} to all incoming sound angles, providing a new avenue for studies of binaural hearing and further inspiration for fly-inspired technologies. A Novel Paradigm for Deep Reinforcement Learning of Biomimetic Systems Mechanisms behind neural control of movement have been an active area of research. Goal-directed movement is a common experimental setup used to understand these mechanisms and relevant neural pathways. On the one hand, optimal feedback control theory is used to model and make quantitative predictions of the coordinated activations of the effectors, such as muscles, joints or limbs. While on the other hand, evidence shows that higher centres such as Basal Ganglia and Cerebellum are involved in activities such as reinforcement learning and error correction. In this paper, we provide a framework to build a digital twin of relevant sections of the human spinal cord using our NEUROiD platform. The digital twin is anatomically and physiologically realistic model of the spinal cord at cellular, spinal networks and system level. We then build a framework to learn the supraspinal activations necessary to perform a simple goal directed movement of the upper limb. The NEUROiD model is interfaced to an Opensim model for all the musculoskeletal simulations. We use Deep Reinforcement Learning to obtain the supraspinal activations necessary to perform the goal directed movement. As per our knowledge, this is the first time an attempt is made to learn the stimulation pattern at the spinal cord level, especially by limiting the observation space to only the afferent feedback received on the Ia, II and Ib fibers. Such a setup results in a biologically realistic constrained environment for learning. Our results show that (1) Reinforcement Learning algorithm converges naturally to the triphasic response observed during goal directed movement (2) Increasing the complexity of the goal gradually helped to accelerate learning (3) Modulation of the afferent inputs were sufficient to execute tasks which were not explicitly learned, but were closely related to the learnt task. GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling.