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.
A chromosome-level genome assembly of a maize elite breeding line Dan340 The maize cultivar Dan340 is an excellent backbone inbred line of Luda Red Cob Group with several desirable characters, such as disease resistance, lodging resistance, high combining ability, wide adaptability and so on. In this study, we constructed a high-quality chromosome-level reference genome for Dan340 by combining PacBio long HiFi sequencing, Illumina short reads and chromosomal conformational capture (Hi-C) sequencing reads. The final assembly of Dan340 genome was 2,348.72 Mb, including 2,738 contigs and 2,315 scaffolds with a N50 of 41.49 Mb and 215.35 Mb, respectively. The percent of high quality Illumina reads mapped to the reference genome was up to 97.48%. The assembly of this genome will not only facilitate our understanding about intraspecific genome diversity in maize, but also provides a novel resource for maize breeding.

Key findingsO_LIThe final assembly of Dan340 genome was 2,348.72 Mb, including 2,738 contigs and 2,315 scaffolds with a N50 of 41.49 Mb and 215.35 Mb, respectively.
C_LIO_LIThe percent of reads mapped to the reference genome was up to 97.48%.
C_LIO_LIThe results showed that 96.84% of the plant single-copy orthologues were complete. Complete single-copy and multi copy genes accounted for 87.36% and 9.48% of the genes, respectively. Taken together, these results indicated that our Dan340 genome assembly presented high quality and coverage.
C_LI
Cross-dataset reproducibility of human retinotopic maps Population receptive field (pRF) models fit to fMRI data are used to non-invasively measure retinotopic maps in human visual cortex, and these maps are a fundamental component of visual neuroscience experiments. Here, we examined the reproducibility of retinotopic maps across two datasets: a newly acquired retinotopy dataset from New York University (NYU) (n=44) and a public dataset from the Human Connectome Project (HCP) (n=181). Our goal was to assess the degree to which pRF properties are similar across datasets, despite substantial differences in their experimental protocols. The two datasets simultaneously differ in their stimulus apertures, participant pool, fMRI protocol, MRI field strength, and preprocessing pipeline. We assessed the cross-dataset reproducibility of the two datasets in terms of the similarity of vertex-wise pRF estimates and in terms of large-scale polar angle asymmetries in cortical magnification. Within V1, V2, V3, and hV4, the group-median NYU and HCP vertex-wise polar angle estimates were nearly identical. Both eccentricity and pRF size estimates were also strongly correlated between the two datasets, but with a slope different from 1; the eccentricity and pRF size estimates were systematically greater in the NYU data. Next, to compare large-scale map properties, we quantified two polar angle asymmetries in V1 cortical magnification previously identified in the HCP data. The NYU dataset confirms earlier reports that more cortical surface area represents horizontal than vertical visual field meridian, and lower than upper vertical visual field meridian. Together, our findings show that the retinotopic properties of V1, V2, V3, and hV4 can be reliably measured across two datasets, despite numerous differences in their experimental design. fMRI-derived retinotopic maps are reproducible because they rely on an explicit computational model of the fMRI response. In the case of pRF mapping, the model is grounded in physiological evidence of how visual receptive fields are organized, allowing one to quantitatively characterize the BOLD signal in terms of stimulus properties (i.e., location and size). The new NYU Retinotopy Dataset will serve as a useful benchmark for testing hypotheses about the organization of visual areas and for comparison to the HCP 7T Retinotopy Dataset. If Research is Not the Evidence, What is it? The Lack of Research Citations in Health Vlogs Doctor-influencers have the knowledge needed to understand scientific research, and they have the social media popularity to share it with large audiences. In this article, I explore the treatment of evidence in vlogs by Egyptian doctor-influencers: how they present it and what they use instead of it. I answer the following questions: 1) How do doctor-influencers present evidence on social media? 2) If evidence is not the research, then what is it? The data was collected through in-depth interviews with 12 Egyptian doctor-influencers from different specializations; two focus groups; and a critical discourse analysis of 48 of the influencers most popular and most engaging videos. I found that the doctor-influencers cited academic research only in the videos about medical controversies, new treatments, misinformation, or common medical mistakes. In most of their videos, the Egyptian doctor-influencers scarcely referred to medical research. This is because they believe that their audiences would not understand evidence due to their low educational attainment, weak research skills, and the language barrier (as most medical research is published in English rather than Arabic). These findings highlight the challenges of communicating science in societies with a low level of education. In these settings, communicators may need to move beyond mediating research studies on social media and instead focus on simplifying, translating, and connecting evidence-based information to the cultural context. Three membrane fusion pore families determine the pathway to pore dilation During exocytosis secretory vesicles fuse with a target membrane and release neurotransmitters, hormones or other bioactive molecules through a membrane fusion pore. The initially small pore may subsequently dilate for full contents release, as commonly observed in amperometric traces. The size, shape and evolution of the pore is critical to the course of contents release, but exact fusion pore solutions accounting for membrane tension and bending energy constraints have not been available. Here we obtained exact solutions for fusion pores between two membranes. We find three families: a narrow pore, a wide pore and an intermediate tether-like pore. For high tensions these are close to the catenoidal and tether solutions recently reported for freely hinged membrane boundaries. We suggest membrane fusion initially generates a stable narrow pore, and the dilation pathway is a transition to the stable wide pore family. The unstable intermediate pore is the transition state that sets the energy barrier for this dilation pathway. Pore dilation is mechanosensitive, as the energy barrier is lowered by increased membrane tension. Finally, we study fusion pores in nanodiscs, powerful systems for the study of individual pores. We show that nanodiscs stabilize fusion pores by locking them into the narrow pore family.

SignificanceDuring neurotransmission, hormone release and other fundamental processes, secretory vesicles fuse their membranes with target membranes to release contents through an initially small membrane fusion pore that subsequently dilates. Dilation is assisted by proteins such as SNAREs and synaptotagmin. While macroscopic soap film shapes are well characterized, finding exact solutions for microscopic cellular membrane surfaces is made more complex by bending energy constraints. Here, computational analysis revealed three families of fusion pores between two membranes. Our work suggests membrane fusion generates a member of the narrow pore family, and pore dilation is a transition to the wide pore family. The energy barrier that SNAREs or synaptotagmin must surmount to achieve dilation is set by a third unstable intermediate pore family. Barbary ground squirrels do not have a sentinel system but instead synchronise vigilance Coordinated behavior, such as hunting in lions and coordinated vigilance as antipredator behavior, are examples of benefits of group-living. Instead of asynchronous vigilance, some social species synchronize their vigilance bouts or take turns acting as sentinels. To increase our knowledge on the evolution of vigilance behavior, we studied whether vigilance is coordinated in Barbary ground squirrels, Atlantoxerus getulus. We show that vigilance was synchronized instead of taking turns. Multiple non-mutually exclusive hypotheses could explain synchronization: Barbary ground squirrels may perch because 1) neighbors are perched (copying effect), 2) perch synchrony may be an emergent property of the ecology as all squirrels may be satiated at the same time (collective behavior), or 3) the benefits are large in terms of evading ambush predators and scanning effectiveness (watch each others back). Particularly, in habitats where the field of view is obstructed by man-made structures and multiple individuals may be necessary to watch for terrestrial predators, synchronized vigilance may have greater fitness benefits than sentinel behavior. We conclude that it is essential to test assumptions of coordination and, thus, to analyze coordination to describe sentinel systems.

Significance StatementVigilance behavior can be vital to an animals survival. Taking turns acting as sentinels or synchronizing vigilance bouts reduces the cost of the trade-off between feeding and predation risk. A sentinel system assumes that sentinels are vigilant from raised positions, warn group members of danger, and alternate vigilance bouts. However, the assumption of alternating vigilance bouts remains poorly tested. We tested this assumption in invasive Barbary ground squirrels. We found that instead of alternating, individuals synchronized their vigilance bouts. Perch synchrony may be 1) a response to perching group members (copying effect), 2) an emergent property of the species ecology, and 3) an adaptation to anthropogenically altered habitats (watch each others back). 2,3-cAMP treatment mimics stress molecular response in Arabidopsis thaliana The role of the RNA degradation product 2,3-cyclic adenosine monophosphate (2,3-cAMP) is poorly understood. Recent studies have identified 2,3-cAMP in plant material and determined its role in stress signaling. The level of 2,3-cAMP increases upon wounding, dark, and heat, and 2,3-cAMP by binding to an RNA-binding protein, Rbp47b, promotes stress granule (SG) assembly. To gain further mechanistic insight into 2,3-cAMP function, we used a multi-omics approach combining transcriptomics, metabolomics, and proteomics to dissect Arabidopsis response to 2,3-cAMP treatment. We demonstrated that 2,3-cAMP is metabolized into adenosine, suggesting that the well-known cyclic nucleotide-adenosine pathway from human cells might also exist in plants. Transcriptomic analysis revealed only minor overlap between 2,3-cAMP-and adenosine-treated plants, suggesting that these molecules act through independent mechanisms. Treatment with 2,3-cAMP changed the levels of hundreds of transcripts, proteins, and metabolites, many previously associated with plant stress responses including protein and RNA degradation products, glucosinolates, chaperones and SG components. Finally, we demonstrated that 2,3-cAMP treatment influences the movement of processing bodies, supporting the role of 2,3-cAMP in the formation and motility of membraneless organelles. Identifying and correcting multiple sources of misspecification in GWAS summary statistics for 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. Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems. 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. Experimental evolution reveals post-transcriptional regulation as a novel driver of Leishmania fitness gain The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network driving parasite fitness, with genome instability causing highly reproducible, gene dosage-dependent changes in protein abundance linked to post-transcriptional regulation. These in turn were associated with a gene dosage-independent reduction in abundance of flagellar transcripts and a coordinated increase in abundance of coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of fitness-adapted ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium. Mechanisms of splicing substrate definition by modulation of pre-mRNA 3D structural scaffolds Coordination of different serine-arginine-rich (SR) proteins - a class of critical splicing activators - facilitates recognition of the highly degenerate mammalian splice signal sequences against the background sequences. Yet, the mechanistic details of their actions remain unclear. Here we show that cooperative binding of SR proteins to exonic and intronic motifs remodels the pre-mRNA 3D structural scaffold. The scaffold generated by the pre-mRNA-specific combinations of different SR proteins in an appropriate stoichiometry is recognized by U1 snRNP - the 5' splice site recognizing factor. A large excess of U1 snRNP particles displaces majority of the bound SR protein molecules from the remodeled pre-mRNA. A higher than optimal stoichiometry of SR proteins occludes the binding sites on the pre-mRNA, raising the U1 snRNP levels required for SR protein displacement and potentially impeding spliceosome assembly. This novel splicing-related function of U1 snRNP is important for distinguishing the substrate and the non-substrate by U2AF65 - the primary 3' splice site-recognizing factor. Overall, this work elucidates early regulatory steps of mammalian splicing substrate definition by SR proteins. 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. Low mammalian species richness is associated with Kyasanur Forest disease outbreak risk in deforested landscapes in the Western Ghats, India Kyasanur forest disease virus (KFDV) is a rapidly expanding tick-borne zoonotic virus with natural foci in the forested region of the Western Ghats of South India. The Western Ghats is one of the worlds most important biodiversity hotspots and, like many such areas of high biodiversity, is under significant pressure from anthropogenic landscape change. The current study sought to quantify mammalian species richness using ensemble models of the distributions of a sample of species extant in the Western Ghats and to explore its association with KFDV outbreaks, as well as the modifying effects of deforestation on this association. Species richness was quantified as a composite of individual species distributions, as derived from ensembles of boosted regression tree, random forest, and generalised additive models. Species richness was further adjusted for the potential biotic constraints of sympatric species. Both species richness and forest loss demonstrated strong positive associations with KFDV outbreaks, however forest loss substantially modified the association between species richness and outbreaks. High species richness was associated with increased KFDV risk but only in areas of low forest loss. In contrast, lower species richness was associated with increased KFDV risk in areas of greater forest loss. This relationship persisted when species richness was adjusted for biotic constraints at the taluk-level. In addition, the taluk-level species abundances of three monkey species (Macaca radiata, Semnopithecus hypoleucus, and Semnopithecus priam) were also associated with outbreaks. These results suggest that increased monitoring of wildlife in areas of significant habitat fragmentation may add considerably to critical knowledge gaps in KFDV epidemiology and infection ecology and should be incorporated into novel One Health surveillance development for the region. In addition, the inclusion of some primate species as sentinels of KFDV circulation into general wildlife surveillance architecture may add further value.

HighlightsO_LILocal mammalian species richness is estimated across the entire Western Ghats region
C_LIO_LILow species richness is associated with high KFDV risk in deforested landscapes
C_LIO_LIThis work identifies key landscapes for wildlife disease surveillance development
C_LI 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. Microvesicles Transfer Mitochondria and Increase Mitochondrial Function in Brain Endothelial Cells We have demonstrated, for the first time that microvesicles, a sub-type of extracellular vesicles (EVs) derived from hCMEC/D3: a human brain endothelial cell (BEC) line transfer polarized mitochondria to recipient BECs in culture and to neurons in mice acute brain cortical and hippocampal slices. This mitochondrial transfer increased ATP levels by 100 to 200-fold (relative to untreated cells) in the recipient BECs exposed to oxygen-glucose deprivation, an in vitro model of cerebral ischemia. We have also demonstrated that transfer of microvesicles, the larger EV fraction, but not exosomes resulted in increased mitochondrial function in hypoxic endothelial cultures. Gene ontology and pathway enrichment analysis of EVs revealed a very high association to glycolysis-related processes. In comparison to heterotypic macrophage- derived EVs, BEC-derived EVs demonstrated a greater selectivity to transfer mitochondria and increase endothelial cell survival under ischemic conditions.

HighlightsO_LIMicrovesicles transfer mitochondria to endothelial cells and brain slice neurons
C_LIO_LIMitochondrial transfer increased ATP in ischemic brain endothelial cells (BECs)
C_LIO_LITransfer of microvesicles increased mitochondrial function in hypoxic BECs
C_LIO_LITransfer of exosomes did not affect mitochondrial function in hypoxic BECs
C_LIO_LIHomotypic BEC-derived EVs result in greater ATP levels in the recipient BECs
C_LI Crystal structures of alphavirus nonstructural protein 4 (nsP4) reveal an intrinsically dynamic RNA-dependent RNA polymerase fold. Alphaviruses such as Ross River virus (RRV), chikungunya virus, Sindbis virus (SINV), and Venezuelan equine encephalitis virus are mosquito-borne pathogens that can cause arthritis or encephalitis diseases. Nonstructural protein 4 (nsP4) of alphaviruses possesses RNA-dependent RNA polymerase (RdRp) activity essential for viral RNA replication. No 3D structure has been available for nsP4 of any alphaviruses despite its importance for understanding alphaviral RNA replication and for the design of antiviral drugs. Here, we report crystal structures of the RdRp domain of nsP4 from both RRV and SINV determined at resolutions of 2.6 and 1.9 [A]. The structure of the alphavirus RdRp domain appears most closely related to RdRps from pestiviruses, noroviruses, and picornaviruses. Hydrogendeuterium exchange mass spectrometry (HDX-MS) and nuclear magnetic resonance (NMR) methods, showed that in solution, nsP4 is highly dynamic with an intrinsically disordered N-terminal domain. Both full-length nsP4 and the RdRp domain were capable to catalyze RNA polymerization. Structure-guided mutagenesis using a trans-replicase system identified nsP4 regions critical for viral RNA replication.

Key PointsO_LICrystal structures of alphavirus nsP4 RNA polymerase domain from RRV and SINV.
C_LIO_LInsP4 protein is highly dynamic with an intrinsically disordered N-terminal domain.
C_LIO_LIOptimized RNA elongation activity assay to facilitate antiviral discovery.
C_LI SARS-CoV-2 variants resist antibody neutralization and broaden host ACE2 usage New SARS-CoV-2 variants continue to emerge from the current global pandemic, some of which can replicate faster and with greater transmissibility and pathogenicity. In particular, UK501Y.V1 identified in UK, SA501Y.V2 in South Africa, and BR501Y.V3 in Brazil are raising serious concerns as they spread quickly and contain spike protein mutations that may facilitate escape from current antibody therapies and vaccine protection. Here, we constructed a panel of 28 SARS-CoV-2 pseudoviruses bearing single or combined mutations found in the spike protein of these three variants, as well as additional nine mutations that within or close by the major antigenic sites in the spike protein identified in the GISAID database. These pseudoviruses were tested against a panel of monoclonal antibodies (mAbs), including some approved for emergency use to treat SARS-CoV-2 infection, and convalescent patient plasma collected early in the pandemic. SA501Y.V2 pseudovirus was the most resistant, in magnitude and breadth, against mAbs and convalescent plasma, followed by BR501Y.V3, and then UK501Y.V1. This resistance hierarchy corresponds with Y144del and 242-244del mutations in the N-terminal domain as well as K417N/T, E484K and N501Y mutations in the receptor binding domain (RBD). Crystal structural analysis of RBD carrying triple K417N-E484K-N501Y mutations found in SA501Y.V2 bound with mAb P2C-1F11 revealed a molecular basis for antibody neutralization and escape. SA501Y.V2 and BR501Y.V3 also acquired substantial ability to use mouse and mink ACE2 for entry. Taken together, our results clearly demonstrate major antigenic shifts and potentially broadening the host range of SA501Y.V2 and BR501Y.V3, which pose serious challenges to our current antibody therapies and vaccine protection. JACKIE: Fast enumeration of genomic single- and multi-copy target sites and their off-targets for CRISPR and other engineered nuclease systems ZFP-, TALE-, and CRISPR-based methods for genome, epigenome editing and imaging have provided powerful tools to interrogate functions of genomes. Targeting sequence design is vital to the success of these experiments. While existing design software mainly focus on designing target sequence for specific elements, we report here the implementation of JACKIE (Jackie and Alberts Comprehensive K-mer Instances Enumerator), a suite of software for enumerating all single- and multi-copy sites in the genome that can be incorporated for genome-scale designs as well as loaded onto genome browsers alongside other tracks for convenient web-based graphic-user-interface (GUI)-enabled design. We also implement fast algorithms to identify sequence neighborhoods or off-target counts of targeting sequences so that designs with low probability of off-target can be identified among millions of design sequences in reasonable time. We demonstrate the application of JACKIE-designed CRISPR site clusters for genome imaging. The limits of evolutionary convergence in sympatry: reproductive interference and historical constraints leading to local diversity in warning traits Mutualistic interactions between defended species represent a striking case of evolutionary convergence in sympatry, driven by the increased protection against predators brought by mimicry in warning traits. However, such convergence is often limited: sympatric defended species frequently display different or imperfectly similar warning traits. The phylogenetic distance between sympatric species may indeed prevent evolution towards the exact same signal. Moreover, warning traits are also involved in mate recognition, so that trait convergence might result in heterospecific courtship and mating. Here, we develop a mathematical model to investigate the strength and direction of evolution of warning trait in defended species with different ancestral traits. Specifically, we determine the effect of phenotypic distances between ancestral trait states of sympatric defended species and of costs of heterospecific sexual interactions on imperfect mimicry and trait divergence. Our analytical results confirm that reproductive interference and historical constraints limit the convergence of warning traits, leading to either complete divergence, or imperfect mimicry. Our model reveals that imperfect mimicry evolves only when ancestral trait values differ between species, due to historical constraints and highlights the importance of female and predator discriminations in the evolution of such imperfect mimicry. Our study thus provides new predictions on how reproductive interference interacts with historical constraints and may promote the emergence of novel warning traits, enhancing mimetic diversity. Precartilage condensation during limb skeletogenesis occurs by tissue phase separation controlled by a bistable cell-state switch with suppressed oscillatory dynamics The tetrapod limb skeleton is initiated in unpatterned limb bud mesenchyme by the formation of precartilage condensations. Here, based on time-lapse videographic analysis of a forming condensation in a high-density culture of chicken limb bud mesenchyme, we observe a phase transition to a more fluidized state for cells within spatial compacted foci (protocondensations that will progress to condensations), as reflected in their spatial confinement, cell-substratum interaction and speed of motion. Previous work showed that galectin-8 and galectin-1A, two proteins of the galactoside-binding galectin family, are the earliest determinants of this process in the chicken limb bud, and that their interactions in forming skeletogenic patterns of condensations can be interpreted mathematically through a reaction-diffusion-adhesion framework. Based on this framework, we use an ordinary differential equation-based approach to analyze the core switching modality of the galectin reaction network and characterize the states of the network independent of the diffusive and adhesive arms of the patterning mechanism. We identify two steady states where the concentrations of both galectins are respectively, negligible, and very high. An explicit Lyapunov function shows that there are no periodic solutions. For sigmoidal galectin production terms, the model exhibits a bistable switch that arises from a monostable state via saddle-node bifurcation. Our model therefore predicts that the galectin network exists in low and high expression states separated in space or time without any intermediate states. This provides a causal basis for the observed outside vs. inside transition observed in the in vitro video data. We performed a quantitative analysis of the distribution of galectin-1A in cultures of condensing chick limb mesenchymal cells and found that the interior of the protocondensations had concentrations of this protein (compared to the immediate exterior) over and above that expected from its higher cell density, consistent with the models predictions. The galectin-based patterning network is thus suggested, on theoretical grounds, to incorporate a core switch independent of any spatial or temporal dynamics, that drives the chondrogenic cell state transition in limb skeletogenesis. Prostaglandin EP3 Receptor signaling is required to prevent insulin hypersecretion and metabolic dysfunction in a non-obese mouse model of insulin resistance When homozygous for the LeptinOb mutation (Ob), Black-and-Tan Brachyury (BTBR) mice become morbidly obese and severely insulin resistant, and by 10 weeks of age, frankly diabetic. Previous work has shown Prostaglandin EP3 Receptor (EP3) expression and activity is up-regulated in islets from BTBR-Ob mice as compared to lean controls, actively contributing to their beta-cell dysfunction. In this work, we aimed to test the impact of beta-cell-specific EP3 loss on the BTBR-Ob phenotype by crossing Ptger3 floxed mice with the Rat insulin promoter (RIP)-CreHerr driver strain. Instead, germline recombination of the floxed allele in the founder mouse - an event whose prevalence we identified as directly associated with underlying insulin resistance of the background strain - generated a full-body knockout. Full-body EP3 loss provided no diabetes protection to BTBR-Ob mice, but, unexpectedly, significantly worsened BTBR-lean insulin resistance and glucose tolerance. This in vivo phenotype was not associated with changes in beta-cell fractional area or markers of beta-cell replication ex vivo. Instead, EP3-null BTBR-lean islets had essentially uncontrolled insulin hypersecretion. The selective up-regulation of constitutively-active EP3 splice variants in islets from young, lean BTBR mice as compared to C57BL/6J, where no phenotype of EP3 loss has been observed, provides a potential explanation for the hypersecretion phenotype. In support of this, high islet EP3 expression in Balb/c females vs. Balb/c males was fully consistent with their sexually-dimorphic metabolic phenotype after loss of EP3-coupled Gz protein. Taken together, our findings provide a new dimension to the understanding of EP3 as a critical brake on insulin secretion.

New and NoteworthyIslet Prostaglandin EP3 receptor (EP3) signaling is well-known as up-regulated in the pathophysiological conditions of type 2 diabetes, contributing to beta-cell dysfunction. Unexpected findings in mouse models of non-obese insulin sensitivity and resistance provide a new dimension to our understanding of EP3 as a key modulator of insulin secretion. A previously-unknown relationship between mouse insulin resistance and the penetrance of Rat insulin promoter-driven germline floxed allele recombination is critical to consider when creating beta-cell-specific knockouts.

For Table of Contents Use Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=121 SRC="FIGDIR/small/671289v3_ufig1.gif" ALT="Figure 1">
View larger version (37K):
org.highwire.dtl.DTLVardef@15cf1d5org.highwire.dtl.DTLVardef@108f565org.highwire.dtl.DTLVardef@126e479org.highwire.dtl.DTLVardef@61fdf8_HPS_FORMAT_FIGEXP  M_FIG C_FIG Spatial regulation of the anaphase-promoting complex/cyclosome by a microtubule adaptor The anaphase-promoting complex/cyclosome (APC/C) coordinates advancement through mitosis via temporally controlled polyubiquitination of effector proteins. Despite the long-appreciated spatial organization of key events in mitosis mediated largely by cytoskeletal networks, the spatial regulation of APC/C, the major mitotic E3 ligase, is poorly understood. Here, we describe a microtubule-resident protein, PLEKHA5, as an interactor of APC/C and spatial regulator of its activity in mitosis. PLEKHA5 knockdown delayed mitotic progression, causing accumulation of APC/C substrates dependent upon the PLEKHA5-APC/C interaction. To assess the spatial importance of this interaction, we developed microtubule-localized proximity biotinylation tools, which revealed that depletion of PLEKHA5 decreased the extent of APC/C association with the microtubule cytoskeleton. This decreased APC/C microtubule-localization in turn prevented efficient loading of APC/C with its co-activator CDC20, leading to defects in E3 ligase catalytic activity. We propose that PLEKHA5 functions as an adaptor of APC/C that promotes its subcellular localization to microtubules and facilitates its activation by CDC20, thus ensuring the timely turnover of key mitotic APC/C substrates and proper progression through mitosis. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation.