Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses.
Choice of pre-processing pipeline influences clustering quality of scRNA-seq datasets BackgroundSingle-cell RNA sequencing (scRNA-seq) has quickly become one of the most dominant techniques in modern transcriptome assessment. In particular, 10X Genomics Chromium system, with its high throughput approach, turn key and thorough user guide made this cutting-edge technique accessible to many laboratories using diverse animal models. However, standard pre-processing, including the alignment and cell filtering pipelines might not be ideal for every organism or tissue. Here we applied an alternative strategy, based on the pseudoaligner kallisto, on twenty-two publicly available single cell sequencing datasets from a wide range of tissues of eight organisms and compared the results with the standard 10X Genomics Cell Ranger pipeline.

ResultsIn most of the tested samples, kallisto produced higher sequencing read alignment rates and total gene detection rates in comparison to Cell Ranger. Although datasets processed with Cell Ranger had higher cell counts, outside of human and mouse datasets, these additional cells were routinely of low quality, containing low gene detection rates. Thorough downstream analysis of one kallisto processed dataset, obtained from the zebrafish pineal gland, revealed clearer clustering, allowing the identification of an additional photoreceptor cell type that previously went undetected. The finding of the new cluster suggests that the photoreceptive pineal gland is essentially a bi-chromatic tissue containing both green and red cone-like photoreceptors and implies that the alignment and pre-processing pipeline can affect the discovery of biologically-relevant cell types.

ConclusionWhile Cell Ranger favors higher cell numbers, using kallisto results in datasets with higher median gene detection per cell. We could demonstrate that cell type identification was not hampered by the lower cell count, but in fact improved as a result of the high gene detection rate and the more stringent filtering. Depending on the acquired dataset, it can be beneficial to favor high quality cells and accept a lower cell count, leading to an improved classification of cell types.
Structural and Functional Analyses Explain Pea KAI2 Receptor Diversity and Reveal Stereoselective Catalysis During Signal Perception KAI2 are plant /{beta} hydrolase receptors, which perceive smoke-derived butenolide signals (karrikins) and putative endogenous, yet unidentified phytohormones (KAI2-ligands, KLs). The number of functional KAI2 receptors varies among plant species. It has been suggested that KAI2 gene duplication and sub-functionalization plays an adaptative role for diverse environments or ligand diversification by altering the receptor responsiveness to specific KLs. Legumes represent one of the largest families of flowering plants and contain many essential agronomic crops. Prior to legume diversification, KAI2 underwent duplication, resulting in KAI2A and KAI2B. Integrating plant genetics, ligand perception and enzymatic assays, and protein crystallography, we demonstrate that Pisum sativum KAI2A and KAI2B act as receptors and enzymes with divergent ligand stereoselectivity. KAI2B has a stronger affinity than KAI2A towards the KAI2-ligand (-)-GR24 and remarkably hydrolyses a broader range of substrates including the strigolactone-like isomer (+)-GR24. We determine the crystal structures of PsKAI2B in apo and butenolide-bound states. The biochemical and structural analyses as well as recorded mass spectra of KAI2s reveal a transient intermediate on the catalytic serine and a stable adduct on the catalytic histidine, further illuminating the role of KAI2 not only as receptors but also as bona fide enzymes. Our work uncovers the stereoselectivity of ligand perception and catalysis by evolutionarily diverged KAI2 receptors in KAR/KL signaling pathways and proposes adaptive sensitivity to KAR/KL and strigolactone phytohormones by KAI2B. Identification of the three zinc-binding sites on Tau protein Tau protein has been extensively studied due to its key roles in microtubular cytoskeleton regulation and in the formation of aggregates found in some neurodegenerative diseases. Recently it has been shown that zinc is able to induce tau aggregation by interacting with several binding sites. However, the precise location of these sites and the molecular mechanism of zinc-induced aggregation remain unknown. Here we used Nuclear Magnetic Resonance (NMR) to identify zinc binding sites on hTau40 isoform. These experiments revealed three distinct zinc binding sites on tau, located in the N-terminal part (H14, H32, H94, and H121), the repeat region (H299, C322, H329 and H330) and the C-terminal part (H362, H374, H388 and H407). Further analysis enabled us to show that the C-terminal and the N-terminal sites are independent of each other. Using molecular simulations, we modeled the structure of each site in a complex with zinc. Given the clinical importance of zinc in tau aggregation, our findings pave the way for designing potential therapies for tauopathies.

HighlightsO_LIZinc is known to induce tau aggregation in neurodegenerative diseases
C_LIO_LIZinc binding locations and mechanism are not yet clear
C_LIO_LIUsing NMR we localized 3 zinc binding site on tau
C_LIO_LIBy molecular simulations, we proposed a modeled structure of each site
C_LIO_LIOur findings pave the way for designing potential therapies for tauopathies
C_LI Rbec: a tool for analysis of amplicon sequencing data from synthetic microbial communities SummarySynthetic microbial communities (SynComs) constitute an emerging and powerful tool in biological, biomedical, and biotechnological research. Despite recent advances in algorithms for the analysis of culture-independent amplicon sequencing data from microbial communities, there is a lack of tools specifically designed for analysing SynCom data, where reference sequences for each strain are available. Here we present Rbec, a tool designed for the analysis of SynCom data that accurately corrects PCR and sequencing errors in amplicon sequences and identifies intra-strain polymorphic variation. Extensive evaluation using mock bacterial and fungal communities show that our tool outperforms current methods for samples of varying complexity, diversity, and sequencing depth. Furthermore, Rbec also allows accurate detection of contaminants in SynCom experiments.

Availability and implementationRbec is freely available as an open-source multi-platform R package. Release versions can be obtained via Bioconductor. The developer version is maintained and can be downloaded at: https://github.com/PengfanZhang/Rbec.

Contactgarridoo@mpipz.mpg.de High expression of VRT2 increases the number of rudimentary basal spikelets in wheat Spikelets are the fundamental building blocks of Poaceae inflorescences and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat, the central spikelets produce the most and largest grains, while spikelet size gradually decreases acro- and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient correlates with the developmental age of spikelets, however the basal spikelets are developed first and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterise gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central and basal sections of a single spike than between any section belonging to consecutive developmental timepoints. We found that SVP MADS-box transcription factors, including VRT-A2, are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEP1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modelling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programmes results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain new insights into developmental genetics pathways of grass inflorescences.

One sentence summaryLarge transcriptional gradients exist within a wheat spike and are associated with rudimentary basal spikelet development, resulting in the characteristic lanceolate shape of wheat spikes. Optogenetic actuator/ERK biosensor circuits identify MAPK network nodes that shape ERK dynamics Combining single-cell measurements of ERK activity dynamics with perturbations provides insights into the MAPK network topology. We built circuits consisting of an optogenetic actuator to activate MAPK signaling and an ERK biosensor to measure single-cell ERK dynamics. This allowed us to conduct RNAi screens to investigate the role of 50 MAPK proteins in ERK dynamics. We found that the MAPK network is robust against most node perturbations. We observed that the ERK-RAF and the ERK-RSK2-SOS negative feedbacks operate simultaneously to regulate ERK dynamics. Bypassing the RSK2-mediated feedback, either by direct optogenetic activation of RAS, or by RSK2 perturbation, sensitized ERK dynamics to further perturbations. Similarly, targeting this feedback in a human ErbB2-dependent oncogenic signaling model increased the efficiency of a MEK inhibitor. The RSK2-mediated feedback is thus important for the ability of the MAPK network to produce consistent ERK outputs and its perturbation can enhance the efficiency of MAPK inhibitors. A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. Effects of substrate on shrimp growth, water quality and bacterial community in the biofloc system nursing Penaeus vannamei postlarvae This study aimed to investigate the effects of substrate on water quality, shrimp growth and bacterial community in the biofloc system with a salinity of 5{per thousand}. Two treatments, biofloc system with (sB) or without (nB) addition of elastic solid packing filler (nylon) as substrate, were set up. Penaeus vannamei postlarvae (PL, [~] stage 15) were stocked at a density of 4000 PL m-3 in each treatment with triplicates for a 28-days culture experiment, taking glucose as carbon source (C:N 15:1). Results showed that the survival rate (96.3{+/-}3.6%), FCR (0.76{+/-}0.06) and productivity (1.54{+/-}0.12 kg m-3) in sB treatment were significantly better than those in nB treatment (81.0{+/-}7.1%, 0.98{+/-}0.08 and 1.14{+/-}0.09 kg m-3, P<0.05). All water parameters were in the recommended ranges. Substrate showed significant effect on TAN, TSS, turbidity, biofloc volume, pH and carbonate alkalinity (P < 0.05). Actinobacteria (4.0-22.7%), Bacteroidetes (10.4-33.5%), Firmicutes (0.2-11.2%), Planctomycetes (4.0-14.9%) and Proteobacteria (29.4-59.0%) were the most dominant phyla for both treatments. However, the bacterial community in sB treatment showed to be significantly different from that in nB treatment (Jaccard distance 0.94{+/-}0.01, P=0.001). Substrate showed significant effects on Shannon, Heip, Pielou and Simpson index, as well as relative abundances of Actinobacteria, Bacteroidetes and Planctomycetes (P < 0.05). The results suggested that addition of substrate affected the shrimp growth, water quality and bacterial community in the biofloc system nursing P. vannamei PL with a 5{per thousand} salinity. Rec8 cohesin-mediated axis-loop chromatin architecture is required for meiotic recombination During meiotic prophase, cohesin-dependent axial structures are formed in the synaptonemal complex (SC). However, the functional correlation between these structures and cohesion remains elusive. Here, we examined the formation of cohesin-dependent axial structures in the fission yeast Schizosaccharomyces pombe. This organism forms atypical SCs composed of linear elements (LinEs) resembling the lateral elements of SC but lacking the transverse filaments. Hi-C analysis using a highly synchronous population of meiotic S. pombe cells revealed that the axis-loop chromatin structure formed in meiotic prophase was dependent on the Rec8 cohesin complex. In contrast, the Rec8-mediated formation of the axis-loop structure occurred in cells lacking components of LinEs. To dissect the functions of Rec8, we identified a rec8-F204S mutant that lost the ability to assemble the axis-loop structure without losing cohesion of sister chromatids. This mutant showed defects in the formation of the axis-loop structure and LinE assembly and thus exhibited reduced meiotic recombination. Collectively, our results demonstrate that the Rec8-dependent axis-loop structure provides a structural platform essential for LinE assembly, facilitating meiotic recombination of homologous chromosomes, independently of its role in sister chromatid cohesion. Chemomechanically voxelated niches for programmable histogenesis Stem-cell derived tissue models are commonly cultured under globally-delivered stimuli that trigger histogenesis via self-organizing activity. However, the culture of such tissue models is prone to stochastic behavior, limiting the reproducibility of cellular composition and resulting in non-physiological architectures. To overcome these shortcomings, we developed a method for printing cell niche microenvironments with microstructured cues that mediate local histogenic processes, including mechanosensing and differentiation of selected cell types. Microstructured cues include independently tunable mechano-chemical properties, with conjugated peptides, proteins, and morphogens across a range of Youngs moduli. By rationally designing niches, we mediate the structure of tissues derived from stem-cell-progenitor sources, including a bone-fat assembly from stromal mesenchyme, and embryonic tissues derived from hiPSC. We show that microstructured cues can recapitulate mechano-chemical signals resembling early embryonic histogenesis. This outcome includes a role for niche mechanics in human embryonic organization, where soft niche mechanics bias markers of mesendodermal differentiation and epithelial-to-mesenchymal-transition (EMT), as well as a demonstration of a material-mediated morphogen signaling centers able to induce foci of mesenchymal and EMT differentiation. Thus, microstructured materials can mediate local histogenic processes to enhance the structure and composition of tissue models. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. Highly efficient synthetic CRISPR RNA/Cas9-based mutagenesis for rapid cardiovascular phenotypic screening in F0 zebrafish The zebrafish is a valuable vertebrate model to study cardiovascular formation and function due to the facile visualization and rapid development of the circulatory system in its externally growing embryos. Despite having distinct advantages, zebrafish have paralogs of many important genes, making reverse genetics approaches inefficient since generating animals bearing multiple gene mutations requires substantial efforts. Here, we present a simple and robust synthetic CRISPR RNA/Cas9-based mutagenesis approach for generating biallelic F0 zebrafish knockouts. Using a dual-guide synthetic CRISPR RNA/Cas9 ribonucleoprotein (dgRNP) system, we compared the efficiency of biallelic gene disruptions following the injections of one, two, and three dgRNPs per gene into the cytoplasm or yolk. We show that simultaneous cytoplasmic injections of three distinct dgRNPs per gene into one-cell stage embryos resulted in the most efficient and consistent biallelic gene disruptions. Importantly, this triple dgRNP approach enables efficient inactivation of cell autonomous and cell non-autonomous gene function, likely due to the low mosaicism of biallelic disruptions. In support of this finding, we provide evidence that the F0 animals generated by this method fully phenocopied the endothelial and peri-vascular defects observed in corresponding stable mutant homozygotes. Moreover, this approach faithfully recapitulated the trunk vessel phenotypes resulting from the genetic interaction between two vegfr2 zebrafish paralogs. Mechanistically, investigation of genome editing and mRNA decay indicates that the combined mutagenic actions of three dgRNPs per gene lead to an increased probability of frameshift mutations, enabling efficient biallelic gene disruptions. Therefore, our approach offers a highly robust genetic platform to quickly assess novel and redundant gene function in F0 zebrafish. The renal inflammatory network of nephronophthisis STRUCTURED ABSTRACTO_ST_ABSBACKGROUNDC_ST_ABSThe majority of genetic kidney disease leading to kidney failure is caused by mutations in ciliary genes. How cilia malfunction leads to progressive kidney damage is poorly understood, but recent evidence links ciliopathy genes to CCL2 dependent macrophage recruitment in autosomal dominant polycystic kidney disease (ADPKD), the most studied renal ciliopathy. Whether or not renal inflammation is involved in other renal ciliopathies is unclear.

METHODSWe combined mice models with kidney biopsies and renal epithelial cells sampled from human urine to characterize the renal inflammatory network of nephronophthisis (NPH), the most frequent renal ciliopathy in children.

RESULTSIn human, mutations in cilia genes involved in NPH enhance urine excretion of the chemokine CCL2, causing abnormal macrophage recruitment in kidney tissues from NPH patients. Differing from ADPKD, inactivating Ccl2 specifically in mouse tubular cells does not rescue the NPH phenotype, suggesting that other inflammatory mediators are involved. Using transcriptional data from 2 NPH models, we identify a set of pro-inflammatory cytokines upregulated in this disease, independently of CCL2. The majority of detectable transcripts from this set are specifically upregulated in kidney cells from NPH patients. In line with the function of these cytokines, NPH kidneys show disproportionate neutrophils and T cells infiltrates compared to healthy subject or hypertensive and diabetic chronic kidney disease patients.

CONCLUSIONSThis study reveals that inflammation is a central aspect in human NPH and delineates a specific set of inflammatory mediators that regulates immune cell recruitment in human NPH.

SIGNIFICANCE STATEMENTMutations in genes encoding primary cilia proteins are the leading cause of genetic kidney failure. In autosomal dominant polycystic kidney disease (ADPKD), deregulated cilia signaling leads to kidney infiltration by macrophages through the chemokine CCL2. Little is known about renal inflammation in nephronophthisis (NPH), the most frequent pediatric renal ciliopathy. Using NPH mice models, tissues and cells from NPH patients, we unveil renal inflammation as preeminent feature of NPH. Remarkably, the renal inflammatory evoked by ciliary gene mutations in NPH does not overlap with ADPKD: it is CCL2 independent, involves a prominent recruitment of neutrophils and T cells and a specific cytokine signature. This unforeseen findings strengthen the link between primary cilia and renal inflammation. Source identity shapes spatial preference in primary auditory cortex during active navigation Localizing and identifying sensory objects while navigating the environment are fundamental brain functions. However, how individual objects are neuronally represented during unrestricted self-motion is mostly unexplored. We trained gerbils on a behavioral foraging paradigm that required localization and identification of sound-sources during free navigation. Chronic multi-electrode recordings in primary auditory cortex during task performance revealed previously unreported sensory object representations. Strikingly, the egocentric angle preference of the majority of spatially sensitive neurons changed significantly depending on the task-specific identity (outcome association) of the sound-source. Spatial tuning also exhibited larger temporal complexity. Moreover, we encountered egocentrically untuned neurons whose response magnitude differed between source identities. Using a neural network decoder we show that together, these neuronal response ensembles provide spatio-temporally co-existent information about both the egocentric location and the identity of individual sensory objects during self-motion, revealing a novel cortical computation principle for naturalistic sensing.

HighlightsO_LILocalization task during free navigation prompts diverse spatial tuning in gerbil A1
C_LIO_LISpatial preference of individual neurons changes with sound-source identity
C_LIO_LIEgo- and allocentric information are spatio-temporally coexistent in A1 ensembles
C_LIO_LIActive sensing reveals new cortical representations for sensory object identification
C_LI Selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. Omicron BA.1 and BA.2 Variants Increase the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE2 SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. Currently, the most common Omicron variants are the original BA.1 Omicron strain and the BA.2 variant, which became more prevalent since it first appeared. To investigate how these mutations affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ~34% and ~51% stronger binding free energies for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells. Sound localization in world and head-centered space in ferrets The location of sounds can be described in multiple coordinate systems that are defined relative to ourselves, or the world around us. Evidence from neural recordings in animals point towards the existence of both head-centered and world-centered representations of sound location in the brain; however, it is unclear whether such neural representations have perceptual correlates in the sound localization abilities of non-human listeners. Here, we establish novel behavioral tests to determine the coordinate systems in which ferrets can localize sounds. We found that ferrets could learn to discriminate between sound locations that were fixed in either world-centered or head-centered space, across wide variations in sound location in the alternative coordinate system. Using probe sounds to assess broader generalization of spatial hearing, we demonstrated that in both head and world-centered tasks, animals used continuous maps of auditory space to guide behavior. Single trial responses of individual animals were sufficiently informative that we could then model sound localization using speaker position in specific coordinate systems and accurately predict ferrets actions in held-out data. Our results indicate that an animal model in which neurons are known to be tuned to sound location in egocentric and allocentric reference frames can also localize sounds in multiple head and world-centered spaces.

Significance StatementHumans can describe the location of sounds either relative to themselves, or in the world, independent of their momentary position. These different spaces are also represented in the activity of neurons in animals, but its not clear whether non-human listeners also perceive both head and world-centered sound location. Here, we designed two behavioral tasks in which ferrets had to discriminate between two sounds using their position in the world, or relative to the head. Subjects learnt to solve both problems and showed the ability to generalize sound location in each space when presented with infrequent probe sounds. These findings reveal a perceptual correlate of neural sensitivity previously observed in the ferret brain and establish that, like humans, ferrets can access an auditory map of their local environment. Membrane contact probability: an essential and predictive character for the structural and functional studies of membrane proteins One of the unique traits of membrane proteins is that a significant fraction of their hydrophobic amino acids is exposed to the hydrophobic core of lipid bilayers rather than being embedded in the protein interior, which is often not explicitly considered in the protein structure and function predictions. Here, we propose a characteristic and predictive quantity, the membrane contact probability (MCP), to describe the likelihood of the amino acids of a given sequence being in direct contact with the acyl chains of lipid molecules. We show that MCP is complementary to solvent accessibility in characterizing the outer surface of membrane proteins, and it can be predicted for any given sequence with a machine learning-based method by utilizing a training dataset extracted from MemProtMD, a database generated from molecular dynamics simulations for the membrane proteins with a known structure. As the first of many potential applications, we demonstrate that MCP can be used to systematically improve the prediction precision of the protein contact maps and structures. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Synergy of chemotherapy and macrophage depletion leads to B cell-mediated T cell memory activation and durable triple negative breast cancer regression Immunosuppressive elements within the tumor microenvironment such as Tumor Associated Macrophages (TAMs) can present a barrier to successful anti-tumor responses by cytolytic T cells. We employed preclinical syngeneic p53 null mouse models of triple negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was used to successfully treat several highly aggressive TNBC murine mammary tumors and lung metastasis. Using this regimen and single cell RNA sequencing we characterized tumor infiltrating lymphocytes (TILs) including helper T cells and antigen-presenting B cells that were highly enriched in good responders to combination therapy. Using high dimensional imaging techniques, we identified the close spatial localization of B220+ CD86+ activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks post-treatment in one model that also exhibited long-term tumor regression post-treatment. We also characterized the transcriptional and metabolic heterogeneity of TAMS in these two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model is highly expressed and conserved in human claudin-low breast cancers, and high expression of the T12 signature correlated with reduced overall survival. This T12 tumor TAM signature may help identify human claudin-low breast cancer patients that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rationale combinations of immunotherapy.

SignificanceA treatment regimen that harnessed the immunostimulatory effects of cyclophosphamide coupled with the inhibition of CSF1R was used to successfully treat several highly aggressive claudin-low TNBC murine mammary tumors and lung metastasis. In-depth characterization of denitrifier communities across different soil ecosystems in the tundra BackgroundIn contrast to earlier assumptions, there is now mounting evidence for the role of tundra soils as important sources of the greenhouse gas nitrous oxide (N2O). However, the microorganisms involved in the cycling of N2O in this system remain largely uncharacterized. Since tundra soils are variable sources and sinks of N2O, we aimed at investigating differences in community structure across different soil ecosystems in the tundra.

ResultsWe analysed 1.4 Tb of metagenomic data from soils in northern Finland covering a range of ecosystems from dry upland soils to water-logged fens and obtained 796 manually binned and curated metagenome-assembled genomes (MAGs). We then searched for MAGs harbouring genes involved in denitrification, an important process driving N2O emissions. Communities of potential denitrifiers were dominated by microorganisms with truncated denitrification pathways (i.e., lacking one or more denitrification genes) and differed across soil ecosystems. Upland soils showed a strong N2O sink potential and were dominated by members of the Alphaproteobacteria such as Bradyrhizobium and Reyranella. Fens, which had in general net-zero N2O fluxes, had a high abundance of poorly characterized taxa affiliated with the Chloroflexota lineage Ellin6529 and the Acidobacteriota subdivision Gp23.

ConclusionsBy coupling an in-depth characterization of microbial communities with in situ measurements of N2O fluxes, our results suggest that the observed spatial patterns of N2O fluxes in the tundra are related to differences in the composition of denitrifier communities. SARS-CoV-2 infection induces soluble platelet activation markers and PAI-1 in the early moderate stage of COVID-19 IntroductionCoagulation dysfunction and thromboembolism emerge as strong comorbidity factors in severe COVID-19. However, it is unclear when particularly platelet activation markers and coagulation factors dysregulated during the pathogenesis of COVID-19. Here, we sought to assess the levels of coagulation and platelet activation markers at moderate and severe stages of COVID-19 to understand the pathogenesis.

MethodsTo understand this, hospitalized COVID-19 patients with (severe cases that required intensive care) or without pneumonia (moderate cases) were recruited. Phenotypic and molecular characterizations were performed employing basic coagulation tests including PT, APTT, D-Dimer and TFPI. The flow cytometry-based multiplex assays were performed to assess FXI, anti-thrombin, prothrombin, fibrinogen, FXIII, P-selectin, sCD40L, plasminogen, tissue-plasminogen activator (tPA), plasminogen activator inhibitor-1 (PAI-1) and D-Dimer.

ResultsThe investigations revealed induction of plasma P-selectin and CD40 ligand (sCD40L) in moderate COVID-19 cases which were significantly abolished with the progression of COVID-19 severity. Moreover, a profound reduction in plasma tissue factor pathway inhibitor (TFPI) and FXIII were identified particularly in the severe COVID-19. Further analysis revealed fibrinogen induction in both moderate and severe patients. Interestingly, an elevated PAI-1 more prominently in moderate, and tPA particularly in severe COVID-19 cases were observed. Particularly, the levels of fibrinogen and tPA directly correlated with the severity of the disease.

ConclusionsIn summary, induction of soluble P-selectin, sCD40L, fibrinogen and PAI-1 in moderate COVID-19 cases suggests the activation of platelets and coagulation system before patients require intensive care. These findings would help in designing better thromboprophylaxis to limit the COVID-19 severity. Ecological specialization, rather than the island effect, explains morphological diversification in an ancient radiation of geckos Island colonists are often assumed to experience higher levels of phenotypic diversification than continental taxa. However, empirical evidence has uncovered exceptions to this "island effect". Here, we tested this pattern using the geckos of the genus Pristurus from continental Arabia and Africa and the Socotra Archipelago. Using a recently published phylogeny and an extensive morphological dataset, we explore the differences in phenotypic evolution between Socotran and continental taxa. Moreover, we reconstructed ancestral habitat occupancy to examine if ecological specialisation is correlated with morphological change, comparing phenotypic disparity and trait evolution between habitats. We found a heterogeneous outcome of island colonisation. Namely, only one of the three colonisation events resulted in a body size increase. However, in general, Socotran species do not present higher levels or rates of morphological diversification than continental groups. Instead, habitat specialisation explains better the body size and shape evolution in Pristurus. Particularly, the colonisation of ground habitats appears as the main driver of morphological change, producing the highest disparity and evolutionary rates. Additionally, arboreal species show very similar body size and head proportions. These results reveal a determinant role of ecological mechanisms in morphological evolution and corroborate the complexity of ecomorphological dynamics in continent-island systems. 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. Genetic architecture of natural variation of cardiac performance in flies Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP)1. Genome Wide Associations Studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict Transcription Factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans.