Decoding the porcine developing spatial processing system and production of human entorhinal stellate cell-like cells by a direct programming approach. Classic studies investigating how and when the entorhinal cortex (component of the memory processing system of the brain) develops have been based on traditional thymidine autoradiography and histological techniques. In this study, we take advantage of modern technologies to trace at a high resolution, the cellular complexity of the developing porcine medial entorhinal cortex by using single-cell profiling. The postnatal medial entorhinal cortex comprises 4 interneuron, 3 pyramidal neuron and 2 stellate cell populations which emerge from intermediate progenitor and immature neuron populations. We discover four MGE-derived interneurons and one CGE-derived interneuron population as well as several IN progenitors. We also identify two oligodendrocyte progenitor populations and three populations of oligodendrocytes. We perform a proof-of-concept experiment demonstrating that porcine scRNA-seq data can be used to develop novel protocols for producing human entorhinal cells in-vitro. We identified six transcription factors (RUNX1A1, SOX5, FOXP1, MEF2C, TCF4, EYA2) important in neurodevelopment and differentiation from one RELN+ stellate cell population. Using a lentiviral vector approach, we reprogrammed human induced pluripotent stem cells into stellate cell-like cells which expressed RELN, SATB2, LEF1 and BCL11B. Our findings contribute to the understanding of the formation of the brains cognitive memory and spatial processing system and provides proof-of-concept for the production of entorhinal cells from human pluripotent stem cells in-vitro.
Tracheal tube fusion in Drosophila involves release of exosomes from multivesicular bodies Extracellular vesicles (EVs) comprise diverse types of cell-released membranous structures that are thought to play important roles in intercellular communication. Despite extensive work on the formation and functions of EVs in cultured cells, studies of EVs in vivo have remained scarce. We report here that EVs are present in the developing lumen of tracheal tubes in Drosophila embryos. We defined two distinct EV subpopulations, one of which contains the Munc13-4 homologue Staccato (Stac) and is spatially and temporally associated with tracheal tube fusion events. The formation of Stac-positive luminal EVs depends on the tip-cell-specific GTPase Arl3, which is also required for the formation of Stac-positive multivesicular bodies, suggesting that Stac-EVs derive from fusion of Stac-MVBs with the luminal membrane in tip cells during anastomosis formation. The GTPases Rab27 and Rab35 cooperate downstream of Arl3 to promote Stac-MVB formation and tube fusion. We propose that Stac-MVBs act as membrane reservoirs that facilitate tracheal lumen fusion in a process regulated by Arl3, Rab27, Rab35, and Stac/Munc13-4.
Potential geographic distribution and ecological niche of New World dobsonflies (Megaloptera, Corydalidae): the case of the Neartic-Neotropical transition zone The dobsonflies (Megaloptera, Corydalidae) are an interesting group of insects, but among the New World dobsonflies, it is not known the effect of the Neartic-Neotropical transition zone on their biogeographic distribution. Here we studied at the species level, the records on the geographic range of the dobsonflies that occurred in and near the transition zone. We presented potential geographic distributions based on ecological niche models for several species of dobsonflies. Results suggested that the geographic range of dobsonflies in the transition zone is associated to mountainous formations and that most species favor warm climates with higher precipitation rates. Climate types tend to be important for species that show narrow geographic ranges, but precipitation tends to be the most important variable to explain species dispersion. Overall, our models support the dispersion of dobsonflies from the Neotropics to North America and explain the two endemic species in Mexico as the result of the formation of the transition zone. Plausible Mechanistic Insights in Biofilm Eradication Potential of against Candida spp. using In Situ Synthesized Tyrosol Functionalized Chitosan Gold Nanoparticles as a versatile Antifouling Coating on Implant Surfaces In the present study, tyrosol functionalized chitosan gold nanoparticles (Chi-TY-AuNPs) were prepared as an alternative treatment strategy to combat fungal infections. Various biophysical techniques were used to characterize the synthesized Chi-TY-AuNPs. The antifungal and antibiofilm activities of Chi-TY-AuNPs were evaluated against C. albicans and C. glabrata and efforts have been made to elucidate the possible mechanism of action. Chi-TY-AuNPs showed a high fungicidal effect against both sessile and planktonic cells of Candida spp. Additionally, Chi-TY-AuNPs completely eradicated (100%) the mature biofilms of both the Candida spp. FESEM analysis highlighted the morphological alterations in Chi-TY-AuNPs treated Candida biofilm cells. Effect of Chi-TY-AuNPs on the ECM components showed significant reduction in protein content in C. glabrata biofilm and substantial decrease in extracellular DNA (eDNA) content of both the Candida spp. ROS generation analysis using DCFDA-PI staining showed high ROS levels in both the Candida spp., whereas pronounced ROS production was observed in Chi-TY-AuNPs treated C. glabrata biofilm. Biochemical analysis revealed decreased ergosterol content in Chi-TY-AuNPs treated C. glabrata cells, while inconsequential changes were observed in C. albicans. Furthermore, the transcriptional expression of selected genes (ergosterol biosynthesis, efflux, sterol importer, and glucan biogenesis) was reduced in C. glabrata in response to Chi-TY-AuNPs except ERG11 and CDR1. Conclusively the result showed the biofilm inhibition and biofilm eradication efficacy of Chi-TY-AuNPs in both the Candida spp. Findings of the present study manifest Chi-TY-AuNPs as a potential therapeutic solution to Candida biofilm-related chronic infections and overcome biofilm antifungal resistance. To knot or not to knot: Multiple conformations of the SARS-CoV-2 frameshifting RNA element The SARS-CoV-2 frameshifting RNA element (FSE) is an excellent target for therapeutic intervention against Covid-19. This small gene element employs a shifting mechanism to pause and backtrack the ribosome during translation between Open Reading Frames 1a and 1b, which code for viral polyproteins. Any interference with this process has profound effect on viral replication and propagation. Pinpointing the structures adapted by the FSE and associated structural transformations involved in frameshifting has been a challenge. Using our graph-theory-based modeling tools for representing RNA secondary structures, "RAG" (RNA-As-Graphs), and chemical structure probing experiments, we show that the 3-stem H-type pseudoknot (3_6 dual graph), long assumed to be the dominant structure has a viable alternative, an HL-type 3-stem pseudoknot (3_3) for longer constructs. In addition, an unknotted 3-way junction RNA (3_5) emerges as a minor conformation. These three conformations share Stems 1 and 3, while the different Stem 2 may be involved in a conformational switch and possibly associations with the ribo-some during translation. For full-length genomes, a stem-loop motif (2_2) may compete with these forms. These structural and mechanistic insights advance our understanding of the SARS-CoV-2 frameshifting process and concomitant virus life cycle, and point to three avenues of therapeutic intervention. 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. Leaf Aqueous Extract of Manilkara hexandra influenced Glucose Metabolism in fish: An observation with Labeo rohita Fingerlings as model organism Manilkara hexandra (Roxb.) Dubard (common name Khirni, family: Sapotaceae) is an evergreen traditionally used medicinal plant. The present study aimed to see the effects of leaf aqueous extract of M. hexandra (LAEMH) on digestive and glucose metabolic enzyme action and serum glucose level in rohu, Labeo rohita, fingerlings. Experimental fish were fed a basal diet (Group I), and diets supplemented with LAEMH at 300 mg (Group II) and 600 mg (Group III) kg-1body weight. A significant reduction in serum glucose was noticed in the treated groups when compared with the control group. The reduced amylase activity was noticed in the treated fingerlings in a dose-dependent manner. However, lipase and protease activities didnt differ significantly among the experimental groups. Reduced serum glucose level might be correlated with the decline in the activity of digestive amylase in fish. Group II (3.35{+/-}0.19 U) and Group III (3.49{+/-}0.13 U) were recorded with reduced glucose-6-phosphatase activities than the control group (4.4{+/-}0.39 U). Moreover, the study revealed a decline in fructose-1,6-phosphatase activities in the treated groups in comparison to the control group. The decline in the activities of the metabolic enzymes might be associated with the non-availability of glucose owing to reduced activity of the digestive amylase in the treated groups. In conclusion, the present study established the hypoglycaemic effect of the leaf aqueous extract of M. hexandra in a fish model. Haplotype-resolved powdery mildew resistance loci reveal the impact of heterozygous structural variation on NLR genes in Muscadinia rotundifolia Muscadinia rotundifolia cv. Trayshed is a valuable source of resistance to grape powdery mildew. It carries two powdery mildew resistance-associated genetic loci, Run1.2 on chromosome 12 and Run2.2 on chromosome 18. In this study, we identified, phased, and reconstructed the two haplotypes of each resistance-associated locus. Haplotype phasing allowed the identification of several structural variation events between haplotypes of both loci. Combined with manual refinement of the gene models, we found that the heterozygous structural variants affected the gene content, with some resulting in duplicated or hemizygous nucleotide-binding leucine-rich repeat (NLR) genes. Structural variations also impacted the domain composition of some NLRs. These findings emphasize the need of generating haplotype-resolved sequences instead of using consensus sequences for identifying haplotype-specific candidate genes. Comparison of the NLRs in the Run1.2 and Run2.2 loci indicated that the two loci are composed of a different number and classes of NLR genes. We provide a list of candidate NLR genes from the Run1.2b and Run2.2 loci, whose expression suggests a role in powdery mildew resistance in Trayshed. These first complete and haplotype-resolved resistance-associated loci, and their candidate NLR genes, represent new resources to develop powdery mildew-resistant grape cultivars. Tau assemblies enter the cytosol in a cholesterol sensitive process essential to seeded aggregation Accumulating evidence supports a prion-like mechanism in the spread of assembled tau in neurodegenerative diseases. Prion-like spread is proposed to require the transit of tau assemblies to the interior of neurons in order to seed aggregation of native, cytosolic tau. This process is poorly understood and remains largely hypothetical. Here, we develop sensitive techniques to quantify the cytosolic entry of tau in real-time. We find that tau does not promote its own entry but, rather, is wholly dependent on cellular machinery. We find that entry to the widely used reporter cell line HEK293 requires clathrin whereas entry to neurons does not. Cholesterol depletion or knockdown of cholesterol transport protein Niemann-Pick type C1 in neurons renders cells highly vulnerable to cytosolic entry and seeded aggregation. Our findings establish entry as the rate-limiting step in seeded aggregation and demonstrate that dysregulated cholesterol, a feature of several neurodegenerative diseases, potentiates tau aggregation.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=170 SRC="FIGDIR/small/449238v3_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@1849023org.highwire.dtl.DTLVardef@13129eorg.highwire.dtl.DTLVardef@2b4736org.highwire.dtl.DTLVardef@15038ac_HPS_FORMAT_FIGEXP  M_FIG Graphical Abstract C_FIG Principles of mRNA targeting via the Arabidopsism6A-binding protein ECT2 Specific recognition of N6-methyladenosine (m6A) in mRNA by RNA-binding proteins containing a YT521-B homology (YTH) domain is important in eukaryotic gene regulation. The Arabidopsis YTH-domain protein ECT2 is thought to bind to mRNA at URU(m6A)Y sites, yet RR(m6A)CH is the canonical m6A consensus site in all eukaryotes and ECT2 functions require m6A binding activity. Here, we apply iCLIP (individual-nucleotide resolution cross-linking and immunoprecipitation) and HyperTRIBE (targets of RNA-binding proteins identified by editing) to define high-quality target sets of ECT2, and analyze the patterns of enriched sequence motifs around ECT2 crosslink sites. Our analyses show that ECT2 does in fact bind to RR(m6A)CH. Pyrimidine-rich motifs are enriched around, but not at m6A-sites, reflecting a preference for N6-adenosine methylation of RRACH/GGAU islands in pyrimidine-rich regions. Such motifs, particularly oligo-U and UNUNU upstream of m6A sites, are also implicated in ECT2 binding via its intrinsically disordered region (IDR). Finally, URUAY-type motifs are enriched at ECT2 crosslink sites, but their distinct properties suggest function as sites of competition between binding of ECT2 and as yet unidentified RNA-binding proteins. Our study provides coherence between genetic and molecular studies of m6A-YTH function in plants, and reveals new insight into the mode of RNA recognition by YTH-domain-containing proteins. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. 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. 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. Human genetic analyses of organelles highlight the nucleus, but not the mitochondrion, in age-related trait heritability Most age-related human diseases are accompanied by a decline in cellular organelle integrity, including impaired lysosomal proteostasis and defective mitochondrial oxidative phosphorylation. An open question, however, is the degree to which inherited variation in or near genes encoding each organelle contributes to age-related disease pathogenesis. Here, we evaluate if genetic loci encoding organelle proteomes confer greater-than-expected age-related disease risk. As mitochondrial dysfunction is a "hallmark" of aging, we begin by assessing nuclear and mitochondrial DNA loci near genes encoding the mitochondrial proteome and surprisingly observe a lack of enrichment across 24 age-related traits. Within nine other organelles, we find no enrichment with one exception: the nucleus, where enrichment emanates from nuclear transcription factors. In agreement, we find that genes encoding several organelles tend to be "haplosufficient," while we observe strong purifying selection against heterozygous protein-truncating variants impacting the nucleus. Our work identifies common variation near transcription factors as having outsize influence on age-related trait risk, motivating future efforts to determine if and how this inherited variation then contributes to observed age-related organelle deterioration. 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. The first genotypic fitness network in a vertebrate reveals that hybridization increases access to novel fitness peaks Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 SNPs that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation. 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. Genetic variation associated with plastic and homeostatic growth responses to drought in Arabidopsis O_LINatural genetic variation influences plant responses to environmental stressors. However, the extent to which such variation underlies plastic versus homeostatic response phenotypes has been little studied.
C_LIO_LIWe quantified the extent of drought-induced changes in leaf area in a set of Iberian Arabidopsis accessions and then performed association studies correlating variation in this phenotype with genomic and transcriptomic variation.
C_LIO_LIDrought-induced plastic reductions in relative leaf area typified accessions originating from productive environments. In contrast, homeostasis in relative leaf area typified accessions originating from unproductive environments. Genome-Wide Association Studies (GWAS), Transcriptome Wide Association Studies (TWAS), and expression GWAS (eGWAS) highlighted the importance of auxin-related processes and, in particular, the potential role of the SMALL AUXIN UP RNA 26 (SAUR26) gene in conferring leaf area plasticity. Homeostatic responses in relative leaf area were associated with a diverse gene set and positively associated with a higher intrinsic water use efficiency (WUEi) as confirmed in a TWAS metanalysis of previously published {delta}13C measurements.
C_LIO_LIWe have identified not only candidate "plasticity genes" but also candidate "homeostasis genes" controlling leaf area. Our results demonstrate the value of a combined GWAS, TWAS, and eGWAS approach to identify mechanisms underlying phenotypic responses to stress.
C_LI

One-sentence summaryInformation on phenotype, genotype, and transcript abundance is integrated to identify both plasticity and homeostasis genes and processes associated with local adaptation to drought stress in Arabidopsis accessions of the Iberian Peninsula. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. Identification of a stereotypic molecular arrangement of endogenous glycine receptors at spinal cord synapses Precise quantitative information about the molecular architecture of synapses is essential to understanding the functional specificity and downstream signaling processes at specific populations of synapses. Glycine receptors (GlyRs) are the primary fast inhibitory neurotransmitter receptors in the spinal cord and brainstem. These inhibitory glycinergic networks crucially regulate motor and sensory processes. Thus far the nanoscale organization of GlyRs underlying the different network specificities has not been defined. Here, we have quantitatively characterized the molecular arrangement and ultra-structure of glycinergic synapses in spinal cord tissue using quantitative super-resolution correlative light and electron microscopy (SR-CLEM). We show that endogenous GlyRs exhibit equal receptor-scaffold occupancy and constant packing densities of about 2000 GlyRs m-2 at synapses across the spinal cord and throughout adulthood, even though ventral horn synapses have twice the total copy numbers, larger postsynaptic domains and more convoluted morphologies than dorsal horn synapses. We demonstrate that this stereotypic molecular arrangement is maintained at glycinergic synapses in the oscillator mouse model of the neuromotor disease hyperekplexia despite a decrease in synapse size, indicating that the molecular organization of GlyRs is preserved in this hypomorph. We thus conclude that the morphology and size of inhibitory postsynaptic specializations rather than differences in GlyR packing determine the postsynaptic strength of glycinergic neurotransmission in motor and sensory spinal cord networks. The Landscape of Receptor-Mediated Precision Cancer Combination Therapy: A Single-Cell Perspective The availability of single-cell transcriptomics data opens new opportunities for rational design of combination cancer treatments. Mining such data, we employed combinatorial optimization techniques to explore the landscape of optimal combination therapies in solid tumors including brain, head and neck, melanoma, lung, breast and colon cancers. We assume that each individual therapy can target any one of 1269 genes encoding cell surface receptors, which may be targets of CAR-T, conjugated antibodies or coated nanoparticle therapies. As a baseline case, we studied the killing of at least 80% of the tumor cells while sparing more than 90% of the non-tumor cells in each patient, as a putative regimen. We find that in most cancer types, personalized combinations composed of at most four targets are then sufficient. However, the number of distinct targets that one would need to assemble to treat all patients in a cohort accordingly would be around 10 in most cases. Further requiring that the target genes be also lowly expressed in healthy tissues uncovers qualitatively similar trends. However, as one asks for more stringent and selective killing beyond the baseline regimen we focused on, we find that the number of targets needed rises rapidly. Emerging individual promising receptor targets include PTPRZ1, which is frequently found in the optimal combinations for brain and head and neck cancers, and EGFR, a recurring target in multiple tumor types. In sum, this systematic single-cell based characterization of the landscape of combinatorial receptor-mediated cancer treatments establishes first of their kind estimates on the number of targets needed, identifying promising ones for future development. Targeting CISH enhances natural cytotoxicity receptor signaling and reduces NK cell exhaustion to improve solid tumor immunity Cytokine inducible SH2-containing protein (CISH) is a natural killer (NK) cell negative regulator of cytokine signaling pathway. To further understand CISH functions in NK cells, we developed a conditional Cish-deficient mouse model in NK cells (Cishfl/flNcr1Ki/+). We detected no developmental or homeostatic difference in NK cells.

However, global gene expression of Cishfl/flNcr1Ki/+ NK cells compared to Cish+/+Ncr1Ki/+ NK cells revealed upregulation of pathways and genes associated with NK cell cycling and activation. We show that CISH does not only regulate interleukin-15 (IL-15) signaling pathways but also natural cytotoxicity receptors (NCR) pathways. Indeed, CISH protein expression level increases upon NCR triggering. Primed Cishfl/flNcr1Ki/+ NK cells display increased activation upon NCR stimulation. Cishfl/flNcr1Ki/+ NK cells display lower activation thresholds and Cishfl/flNcr1Ki/+ mice are more resistant to tumor metastasis. Remarkably, we found that Cishfl/flNcr1Ki/+ mice were also more resistant to primary breast cancer growth in addition to superior control of spontaneous tumor metastasis. CISH deletion favors NK cell accumulation to the primary tumor, optimizes NK cell killing properties and decreases TIGIT immune checkpoint receptor expression, limiting NK cell exhaustion. Finally, we argue that specifically enhancing NK cell function is sufficient to boost anti-tumor response to both primary and secondary tumor models. Using CRISPRi, we then targeted CISH in human NK-92 or primary NK cells. According to the results in our mouse model, CISH deletion favors NCR signaling and anti-tumor functions in human NK cells. Our results validate CISH as an emerging therapeutic target to enhance NK cell immunotherapy. Force correction analysis method for derivation of multidimensional free energy landscapes from adaptively-biased replica simulations A methodology is proposed for the calculation of multidimensional free-energy landscapes of molecular systems, based on analysis of multiple Molecular Dynamics trajectories wherein adaptive biases have been applied to enhance the sampling of different collective variables. In this approach, which we refer to as Force Correction Analysis Method (FCAM), local averages of the total and biasing forces are evaluated post-hoc, and the latter are subtracted from the former to obtain unbiased estimates of the mean force across collective-variable space. Multidimensional free-energy surfaces and minimum free-energy pathways are then derived from integration of the mean force landscape through kinetic Monte Carlo algorithm. To evaluate the proposed method, a series of numerical tests and comparisons with existing approaches were carried out for small molecules, peptides, and proteins, based on all-atom trajectories generated with standard, concurrent and replica-exchange Metadynamics in collective-variable spaces ranging from one- to six-dimensional. The tests confirm the correctness of the FCAM formulation and demonstrate that calculated mean forces and free energies converge rapidly and accurately, outperforming other methods used to unbias this kind of simulation data.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=102 SRC="FIGDIR/small/431654v5_ufig1.gif" ALT="Figure 1">
View larger version (45K):
org.highwire.dtl.DTLVardef@1412289org.highwire.dtl.DTLVardef@1425d53org.highwire.dtl.DTLVardef@1fcf5d0org.highwire.dtl.DTLVardef@6e6941_HPS_FORMAT_FIGEXP  M_FIG TOC/Abstract Graphic

C_FIG Spontaneous differentiation leads to emergence of hybrid cell states relate to poor prognosis in oral cancer PurposeCellular dynamics between phenotypically heterogeneous subpopulations of cancer cells within individual tumor is shown to be responsible for drug tolerance and overall poor prognosis; however, evidences were largely missing in oral cancer. Therefore, this study was undertaken to describe the dynamic phenotypic states among oral cancer cells, its influence on transcriptomic heterogeneity as well as its clinical significance.

Experimental DesignWe multiplexed phenotypic markers of putative oral-stem-like cancer cells (SLCCs) and characterized diversity among CD44-positive oral cancer cell subpopulations with respect to distinct expression of CD24 and aldehyde dehydrogenase (ALDH)-activity in multiple cell lines. Population trajectories were characterized by Markov model and cell states were defined based on the population specific RNA sequencing (RNAseq). ssGSEA based gene expression signatures were explore for prognostic significance.

ResultsOral cancer cells followed two distinct patterns of spontaneous repopulation dynamics with stochastic inter-conversions on  ALDH-axis, however a strict non-interconvertible transition on  CD24-axis. Interestingly, plastic  ALDH-axis was harnessed to enrich ALDHHigh subpopulations in response to Cisplatin treatment, to adapt a drug tolerant state. Phenotype-specific RNAseq results suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate states of stemness within the differentiating oral cancer cells. Further, survival analysis with each subpopulation-specific gene signature strongly suggested that the cell-state dynamics may act as possible mechanism to drive ITH, resulting in poor prognosis in patient.

ConclusionsOur results emphasized the prognostic power of the population dynamics in oral cancer. Importantly, we have described the phenotypic-composition of heterogeneous subpopulations critical for global tumor behaviour in oral cancer; which is a prerequisite knowledge important for precision treatment, however largely lacking for most solid tumors.

Graphical AbstractWe have characterized diversity among CD44-positive oral cancer cells lines with respect to distinct expression of CD24 and ALDH-activity. Subpopulations showed stochastic inter-conversions on ALDH-axis but a strict non-interconvertible transition of CD24Low to CD24High phenotype, even in response to chemotherapy-induced stress. RNAseq study suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate alternate states of stemness within the differentiating oral cancer cells. The described population dynamics demonstrtaed influence tumor behaviour possibly by increasing intratumoral heterogeneity in aggressive oral tumors.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=112 SRC="FIGDIR/small/457509v3_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@1dc6a7org.highwire.dtl.DTLVardef@dd0948org.highwire.dtl.DTLVardef@18c9603org.highwire.dtl.DTLVardef@cca8fb_HPS_FORMAT_FIGEXP  M_FIG C_FIG Translational relevanceIntratumoral heterogeneity (ITH) has been the clinically important factor, impacting aggressive cancer behaviour, drug tolerance and overall poor prognosis. Recent high-throughput studies have provided better cellular and molecular resolution of ITH; however, the prerequisite knowledge which defines the composition of subpopulations critical for global tumor behaviour is majorly lacking for most of the solid tumors. By combining phenotypic markers, we have defined four subpopulations of oral cancer cells. These subpopulations showed stochastic inter-conversions as well as a strict non-interconvertible transition among them to acheive heterogeneity. Importantly, transcriptional states of each subpopulations indicated a clinically relevant signatures for patient prognosis. Also, we observed interconversions of these subpopulations in response to Cisplatin to accumulate drug-tolerant cell state, as rapid and reversible strategy to respond to chemotherapy induced stress. Thus, the characteristics of described phenotypic subgroups may be translated to the clinic for estimating the extent of intratumoral heterogeneity in oral cancer patients.