A dynamical systems treatment of transcriptomic trajectories in hematopoiesis Inspired by Waddingtons illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein the dynamics of an exogenous signal couples to a cells enormously complex signaling and transcriptional machinery, eliciting a qualitative transition in the collective state of a cell - its fate. It remains unclear, however, whether the dynamical systems framework can go beyond a word-based caricature of the system and provide sharp quantitative insights that further our understanding of differentiation. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the individual concentration variations of a myriad of genes. Here, starting from formal mathematical derivations, we challenge these transcriptomic trajectories to a rigorous statistical evaluation of whether they display signatures consistent with bifurcations. After pinpointing bifurcations along transcriptomic trajectories of the neutrophil branch of hematopoeitic differentiation we are able to further leverage the primitive features of a linear instability to identify the single-direction in gene expression space along which the bifurcation unfolds and identify possible gene contributors. This scheme identifies transcription factors long viewed to play a crucial role in the process of neutrophil differentiation in addition to identifying a host of other novel genetic players. Most broadly speaking, we provide evidence that, though very high-dimensional, a bifurcating dynamical systems formalism might be appropriate for the process of cellular differentiation and that it can be leveraged to provide insights. Ambitiously, our work attempts to take a step beyond data-analysis and towards the construction of falsifiable mathematical models that describe the dynamics of the entire transcriptome.
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
New Segmentation and Feature Extraction Algorithm for Classification of White Blood Cells in Peripheral Smear Images This article addresses a new method for the classification of white blood cells (WBCs) using image processing techniques and machine learning methods. The proposed method consists of three steps: detecting the nucleus and cytoplasm, extracting features, and classification. At first, a new algorithm is designed to segment the nucleus. For the cytoplasm to be detected, only a part of it located inside the convex hull of the nucleus is involved in the process. This attitude helps us overcome the difficulties of segmenting the cytoplasm. In the second phase, three shapes and four novel color features are devised and extracted. Finally, by using an SVM model, the WBCs are classified. The segmentation algorithm can detect the nucleus with a dice similarity coefficient of 0.9675. The proposed method can categorize WBCs in Raabin-WBC, LISC, and BCCD datasets with accuracies of 94.65 %, 92.21 %, and 94.20 %, respectively. It is worth mentioning that the hyperparameters of the classifier are fixed only with the Raabin-WBC dataset, and these parameters are not readjusted for LISC and BCCD datasets. The obtained results demonstrate that the proposed method is robust, fast, and accurate. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Do males know? Evidence-driven rainmaking mating strategy by male bonobos to meet the fertile window of females Female bonobos exhibit prolonged receptivity, one suggested function of which is to lower male mating competition. However, it is questionable whether easier access to receptive females can reduce male-male competition, given the exclusive nature of male reproductive success. We tested whether males could determine the fertile phase of females. We found that ovulation probability predicted male mating effort. High-ranking males copulated with females with higher fertility, and male-male agonistic interactions increased when there were fertile females in the party. When there were multiple females with maximal swelling, males concentrated their mating effort on females with an older infant whose maximal swelling started earlier, and they continued mating efforts until detumescence (rainmaking). These findings suggest that male bonobos distinguish between fertile and non-fertile phases of females and that having more receptive females in the party does not reduce male-male competition for fertile females.

TeaserMales use the rainmakers rule to meet the periovulatory phase of the female bonobo for better reproductive success. An Essential Role of UBXN3B in B Lymphopoiesis Hematopoiesis is finely regulated to enable timely production of the right numbers and types of mature immune cells to maintain tissue homeostasis. Dysregulated hematopoiesis may compromise antiviral immunity and/or exacerbate immunopathogenesis. Herein, we report an essential role of UBXN3B in maintenance of hematopoietic homeostasis and restriction of immunopathogenesis during respiratory viral infection. Ubxn3b deficient (Ubxn3b-/-) mice are highly vulnerable to SARS-CoV-2 and influenza A infection, characterized by more severe lung immunopathology, lower virus-specific IgG, significantly fewer B cells, but more myeloid cells than Ubxn3b+/+ littermates. This aberrant immune compartmentalization is recapitulated in uninfected Ubxn3b-/- mice. Mechanistically, UBXN3B controls precursor B-I (pre-BI) transition to pre-BII and subsequent proliferation in a cell-intrinsic manner, by maintaining BLNK protein stability and pre-BCR signaling. These results reveal an essential role of UBXN3B for the early stage of B cell development. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Fixed or random? On the reliability of mixed-effects models for a small number of levels in grouping variables O_LIBiological data are often intrinsically hierarchical. Due to their ability to account for such dependencies, mixed-effects models have become a common analysis technique in ecology and evolution. While many questions around their theoretical foundations and practical applications are solved, one fundamental question is still highly debated: When facing a low number of levels should we model a grouping (blocking, clustering) variable as a random or fixed effect? In such situation, the variance of the random effect is imprecise, but whether this affects the statistical properties of the population effect is unclear.
C_LIO_LIHere, we analyzed the consequences of including a grouping variable as fixed or random effect in the correctly specified and other possible miss-specified models (too simple or too complex models) for data with small number of levels (2 - 8). For all these options, we calculated type I error rates and power. Moreover, we show how these statistical properties change with the study design.
C_LIO_LIWe found that the model choice does not influence the statistical properties of the population effect when the effect is the same at all levels However, if an ecological effect differs among levels, using a random slope and intercept model, and switching to a fixed-effect model only in case of a singular fit, avoids overconfidence in the results. Additionally, power and type I error are strongly influenced by the number of and difference between levels.
C_LIO_LIWe conclude that inferring the correct random effect structure is of high importance to get correct statistical properties. When in doubt, we recommend starting with the simpler model and using model diagnostics to identify missing components. When having identified the correct structure, we encourage to start with a mixed-effects model independent of the number of levels and switch to a fixed-effect model only in case of a singular fit. With these recommendations, we allow for more informative choices about study design and data analysis and thus make ecological inference with mixed-effects models more robust for small number of levels.
C_LI N-VEGF, the forgotten isoform of VEGF-A, activate cellular stress survival circuits Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. The potential of dynamic physiological traits in young tomato plants to predict field-yield performance To address the challenge of predicting tomato yields in the field, we used whole-plant functional phenotyping to evaluate water relations under well-irrigated and drought conditions. The genotypes tested are known to exhibit variability in their yields in wet and dry fields. The examined lines included two lines with recessive mutations that affect carotenoid biosynthesis, zeta z2083 and tangerine t3406, both isogenic to the processing tomato variety M82. The two mutant lines were reciprocally grafted onto M82, and multiple physiological characteristics were measured continuously, before, during and after drought treatment in the greenhouse. A comparative analysis of greenhouse and field yields showed that the whole-canopy stomatal conductance (gsc) in the morning and cumulative transpiration (CT) were strongly correlated with field measurements of total yield (TY: r2 = 0.9 and 0.77, respectively) and plant vegetative weight (PW: r2 = 0.6 and 0.94, respectively). Furthermore, the minimum CT during drought and the rate of recovery when irrigation was resumed were both found to predict resilience. Emergent dynamics of a three-node regulatory network explain phenotypic switching and heterogeneity: a case study of Th1/Th2/Th17 cell differentiation Naive helper (CD4+) T-cells can differentiate into distinct functional subsets including Th1, Th2, and Th17 phenotypes. Each of these phenotypes has a  master regulator - T-bet (Th1), GATA3 (Th2) and ROR{gamma}T (Th17) - that inhibits the other two master regulators. Such mutual repression among them at a transcriptional level can enable multistability, giving rise to six experimentally observed phenotypes - Th1, Th2, Th17, hybrid Th/Th2, hybrid Th2/Th17 and hybrid Th1/Th17. However, the dynamics of switching among these phenotypes, particularly in the case of epigenetic influence, remains unclear. Here, through mathematical modeling, we investigated the coupled transcription-epigenetic dynamics in a three-node mutually repressing network to elucidate how epigenetic changes mediated by any  master regulator can influence the transition rates among different cellular phenotypes. We show that the degree of plasticity exhibited by one phenotype depends on relative strength and duration of mutual epigenetic repression mediated among the master regulators in a three-node network. Further, our model predictions can offer putative mechanisms underlying relatively higher plasticity of Th17 phenotype as observed in vitro and in vivo. Together, our modeling framework characterizes phenotypic plasticity and heterogeneity as an outcome of emergent dynamics of a three-node regulatory network, such as the one mediated by T-bet/GATA3/ROR{gamma}T. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Evaluation of methods for estimating coalescence times using ancestral recombination graphs The ancestral recombination graph (ARG) is a structure that describes the joint genealogies of sampled DNA sequences along the genome. Recent computational methods have made impressive progress towards scalably estimating whole-genome genealogies. In addition to inferring the ARG, some of these methods can also provide ARGs sampled from a defined posterior distribution. Obtaining good samples of ARGs is crucial for quantifying statistical uncertainty and for estimating population genetic parameters such as effective population size, mutation rate, and allele age. Here, we use standard neutral coalescent simulations to benchmark the estimates of pairwise coalescence times from three popular ARG inference programs: ARGweaver, Relate, and tsinfer+tsdate. We compare 1) the true coalescence times to the inferred times at each locus; 2) the distribution of coalescence times across all loci to the expected exponential distribution; 3) whether the sampled coalescence times have the properties expected of a valid posterior distribution. We find that inferred coalescence times at each locus are most accurate in ARGweaver, and often more accurate in Relate than in tsinfer+tsdate. However, all three methods tend to overestimate small coalescence times and underestimate large ones. Lastly, the posterior distribution of ARGweaver is closer to the expected posterior distribution than Relates, but this higher accuracy comes at a substantial trade-off in scalability. The best choice of method will depend on the number and length of input sequences and on the goal of downstream analyses, and we provide guidelines for the best practices. Paternal knockdown of tRNA (cytosine-5-)-methyltransferase (Dnmt2) increases offspring susceptibility to infection in red flour beetles Intergenerational effects from fathers to offspring are increasingly reported from diverse organisms, but the underlying mechanisms remain speculative. Paternal trans-generational immune priming (TGIP) was demonstrated in the red flour beetle Tribolium castaneum: non-infectious bacterial exposure of fathers protects their offspring against an infectious challenge for at least two generations. Epigenetic processes, such as cytosine methylation of nucleic acids, have been proposed to enable transfer of information from fathers to offspring. Here we studied a potential role in TGIP of the Dnmt2 gene (renamed as Trdmt1 in humans), which encodes a highly conserved enzyme that methylates different RNAs, including specific cytosines of a set of tRNAs. Dnmt2 has previously been reported to be involved in intergenerational epigenetic inheritance in mice and protection against viruses in fruit flies. We first studied gene expression and found that Dnmt2 is expressed in various life history stages and tissues of T. castaneum, with high expression in the reproductive organs. RNAi-mediated knockdown of Dnmt2 in fathers was systemic, slowed down offspring larval development and increased mortality of the adult offspring upon bacterial infection. However, these effects were independent of bacterial exposure of the fathers. In conclusion, our results point towards a role of Dnmt2 for paternal effects, while elucidation of the mechanisms behind paternal TGIP needs further studies. CHRM1 is a Druggable Melanoma Target Whose Endogenous Activity is Determined by Inherited Genetic Variation in DOPA Production Melanoma risk is 30 times higher in people with lightly pigmented skin compared to those with darkly pigmented skin. Here we show that this difference results from more than melanin pigment and its ultraviolet radiation (UVR) shielding effect. Using primary human melanocytes representing the full human skin pigment continuum and several preclinical melanoma models, we show that cell-intrinsic differences between dark and light melanocytes regulate melanocyte proliferative capacity, overall cellular differentiation state, and susceptibility to malignant transformation, independently of melanin and UV exposure. We determined that these differences result from dihydroxyphenylalanine (DOPA), a melanin precursor synthesized at higher levels in melanocytes from dark skin. Although DOPA was not previously known to have specific signaling activity, we used both high throughput pharmacologic and genetic in vivo CRISPR screens to determine that DOPA limits melanocyte and melanoma cell proliferation by directly inhibiting the muscarinic acetylcholine receptor M1 (CHRM1), a G Protein-Coupled Receptor (GPCR) not previously known to bind DOPA, nor to affect melanoma pathobiology. Pharmacologic CHRM1 antagonism in melanoma leads to depletion of c-Myc and FOXM1, both of which are proliferation drivers associated with aggressive melanoma. In preclinical mouse melanoma models using both immune deficient and syngeneic immune competent mice, pharmacologic inhibition of CHRM1 or FOXM1 inhibited tumor growth. CHRM1 and FOXM1 may be new therapeutic targets for melanoma. 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. DOMINANT AWN INHIBITOR encodes the ALOG protein originating from gene duplication and inhibits awn elongation by suppressing cell proliferation and elongation in sorghum The awn, a needle-like structure extending from the tip of the lemma in grass species, plays a role in environmental adaptation and fitness. In some crops, awns appear to have been eliminated during domestication. Although numerous genes involved in awn development have been identified, several dominant genes that eliminate awns are also known to exist. For example, in sorghum (Sorghum bicolor), the dominant awn inhibiting gene has been known since 1921; however, its molecular features remain uncharacterized. In this study, we conducted quantitative trait locus analysis and a genome-wide association study of awn-related traits in sorghum and identified DOMINANT AWN INHIBITOR (DAI), which encodes the ALOG family protein on chromosome 3. DAI appeared to be present in most awnless sorghum cultivars, likely because of its effectiveness. Detailed analysis of the ALOG protein family in cereals revealed that DAI originated from duplication of its twin paralog (DAIori) on chromosome 10. Observations of immature awns in near-isogenic lines revealed that DAI inhibits awn elongation by suppressing both cell proliferation and elongation. We also found that only DAI gained a novel function to inhibit awn elongation through an awn-specific expression pattern distinct from that of DAIori. Interestingly, heterologous expression of DAI with its own promoter in rice inhibited awn elongation in the awned cultivar Kasalath. We found that DAI originated from gene duplication, providing an interesting example of gain-of-function that occurs only in sorghum but shares its functionality with rice and sorghum. Chromosome counting in the mouse zygote using low-invasive super-resolution live-cell imaging In preimplantation embryos, an abnormal chromosome number causes developmental failure and a reduction in the pregnancy rate. Conventional chromosome testing methods requiring biopsy reduce the risk of associated genetic diseases; nevertheless, the reduction in cell number also reduces the pregnancy rate. Therefore, we attempted to count the chromosomes in mouse (Slc:ICR) embryos using super-resolution live-cell imaging as a new method of chromosome counting that does not reduce the cell number or viability. We counted the forty chromosomes at the first mitosis by injecting embryos with histone H2B-mCherry mRNA under conditions by which pups could be obtained; however, the results were often an underestimation of chromosome number and varied by embryo and time point. Therefore, we developed a method to count the chromosomes via CRISPR/dCas-mediated live-cell fluorescence in situ hybridization targeting the sequence of the centromere region, enabling us to count the chromosomes more accurately in mouse embryos. The methodology presented here may provide useful information for assisted reproductive technologies, such as those used in livestock animals/humans, as a technique for assessing the chromosomal integrity of embryos prior to transfer.

Impact StatementLow-toxic super-resolution observation enables chromosome counting in preimplantation embryos without cell collection. Spatially localized fMRI metrics as predictive andhighly distinct state-independent fingerprints Individual subjects can be accurately identified in a database based on their functional connectome (FC), i.e. the whole-brain patterns of synchrony in fluctuations of brain activation as measured with functional MRI (fMRI). Accurate identification requires a representation of brain function to be simultaneously reliable (i.e. low intraindividual variability) and idiosyncratic (i.e. high interindividual variability), hence calling it a "fingerprint". Importantly, several studies have shown lower test-retest reliability (i.e. higher intraindividual variability) of the global FC compared to local representations of brain function such as regional homogeneity (ReHo) and (fractional) amplitude of low-frequency fluctuations ((f)ALFF). Therefore, here, with resting-state and task fMRI-data from the Human Connectome Project and the enhanced Nathan Kline Institute, we demonstrate that the local functional fingerprint, and especially regional homogeneity (ReHo), achieves near-perfect identification accuracies using four different ways to quantify identification. This finding is replicable across various parcellations as well as resilient to confounding effects. Moreover, using a machine-learning setup, we show that the small intraindividual ReHo fingerprint variability across sessions is meaningful for explaining individual-level intelligence. Further analyses reveal that the attention networks and the Default Mode Network contributed the most to the individual "uniqueness" of the ReHo fingerprint. Last, with another publicly available dataset, we show that ReHo is also stable in individuals with schizophrenia and that its (in)stability across sessions relates to intelligence scores. Altogether, our findings suggest that the ReHo fingerprint is a good candidate for further exploration of applicability in precision medicine. Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states. 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. Multi-omics investigation of tacrolimus nephrotoxicity BackgroundTacrolimus, an immunosuppressive drug prescribed to a majority of transplanted patients is nephrotoxic, through still unclear mechanisms. This study aims to evaluate the impact of tacrolimus on a lineage of proximal tubular cells using a multi-omics approach.

MethodsLLC-PK1 cells were exposed to 5 M of tacrolimus for 24h. Intracellular proteins and metabolites, and extracellular metabolites were extracted and analysed by LC-MS/MS. The transcriptional expression of the dysregulated proteins PCK-1, FBP1 and FBP2 was measured using RT-qPCR.

ResultsIn our cell model, tacrolimus impacted different metabolic pathways including those of arginine (e.g., citrulline, ornithine) (p < 0.0001), amino acids (e.g., valine, isoleucine, aspartic acid) (p < 0.0001) and pyrimidine (p<0.01). In addition, it induced oxidative stress (p < 0.01) as shown by a decrease in total cell glutathione quantity. It impacted cell energy through an increase in Krebs cycle intermediates (e.g., citrate, aconitate, fumarate) (p < 0.01) and down-regulation of PCK-1 (p < 0.05) and FPB1 (p < 0.01), which are key enzymes in gluconeogenesis. Apart from glucose synthesis, gluconeogenesis is an important process in kidney mediated acid-base balance control.

ConclusionThe variations found using this multi-omics approach clearly point towards a dysregulation of energy production in epithelial cells of the renal tubule, and potentially of their functions, that may be implicated in tacrolimus nephrotoxicity in the clinics. Microbial loop of a Proterozoic ocean analogue Meromictic Lake Cadagno, an ancient ocean analogue, is known for its permanent stratification and persistent anoxygenic microbial bloom within the chemocline. Although the anaerobic microbial ecology of the lake has been extensively studied for at least 25 years, a comprehensive picture of the microbial food web linking the bacterial layer to phytoplankton and viruses, with explicit measures of primary and secondary production, is still missing. This study sought to better understand the abundances and productivity of microbes in the context of nutrient biogeochemical cycling across the stratified zones of Lake Cadagno. Photosynthetic pigments and chloroplast 16S rRNA gene phylogenies suggested the presence of eukaryotic phytoplankton through the water column. Evidence supported high abundances of Ankyra judayi, a high-alpine adapted chlorophyte, in the oxic mixolimnion where oxygenic-primary production peaked. Through the no- and low-oxygen chemocline and monimolimnion, chlorophytes related to Closteriopsis acicularis, a known genus of meromictic lakes, and Parachlorella kessleri were observed. Chromatium, anoxygenic phototrophic sulfur bacteria, dominated the chemocline along with Lentimicrobium, a genus of known fermenters whose abundance was newly reported in Lake Cadagno. Secondary production peaked in the chemocline suggesting primary producers depend on hetrotrophs for nutrient remineralization. As previously observed, sulfur reducing bacteria (SRBs), especially Desulfocapsa and Desulfobulbus, were present in the chemocline and anoxic monimolimnion. Virus-to-microbe ratios (VMR) peaked in the zone of phytoplankton, yet were at a minimum at the peak of Chromatium. These dynamic trends suggest viruses may play a role in the modulation of oxygenic and anoxygenic photo- and chemosynthesis in Lake Cadagno and other permanently stratified systems.

ImportanceAs a window to the past, the study offers insights into the role of microbial guilds of Proterozoic ocean chemoclines in the production and recycling of organic matter of sulfur- and ammonia-containing ancient oceans. The new observations described here suggest that eukaryotic algae were persistent in the low oxygen upper-chemocline in association with purple and green sulfur bacteria in the lower half of the chemocline. Further, this study provides the first insights into Lake Cadagno viral ecology. High viral abundances suggested viruses may be essential components of the chemocline where their activity may result in the release and recycling of organic matter. The framework developed in this study through the integration of diverse geochemical and biological data types lays the foundation for future studies to quantitatively resolve the processes performed by discrete populations comprising the microbial loop in this early anoxic ocean analogue. 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. Measuring cystic fibrosis drug responses in organoids derived from 2D differentiated nasal epithelia Cystic Fibrosis (CF) is caused by genetic defects that impair the cystic fibrosis transmembrane conductance regulator (CFTR) channel in airway epithelial cells. These defects may be overcome by specific CFTR modulating drugs, for which the efficacy can be predicted in a personalized manner using 3D nasal-brushing-derived airway organoids in a forskolin-induced swelling assay. Despite of this, previously described application of 3D airway organoids in CFTR function assays have not been fully optimal. In this report we therefore describe an alternative method of culturing nasal brushing-derived airway organoids, which are created from an equally differentiated airway epithelial monolayer of a 2D air-liquid interface culture. In addition, we have defined organoid culture conditions, with the growth factor/cytokine combination neuregulin-1{beta} and interleukine-1{beta}, which enabled consistent detection of CFTR modulator responses in nasal airway organoids cultures from subjects with CF. Alpha satellite insertion close to an ancestral centromeric region Human centromeres are mainly composed of alpha satellite DNA hierarchically organized as higher-order repeats (HORs). Alpha satellite dynamics is shown by sequence homogenization in centromeric arrays and by its transfer to other centromeric locations, for example during the maturation of new centromeres. We identified during prenatal aneuploidy diagnosis by FISH a de novo insertion of alpha satellite DNA from the centromere of chromosome 18 (D18Z1) into cytoband 15q26. Although bound by CENP-B, this locus did not acquire centromeric functionality as demonstrated by lack of constriction and absence of CENP-A binding. The insertion was associated with a 2.8 kbp deletion and likely occurred in the paternal germline. The site was enriched in long terminal repeats (LTRs) and located ~10 Mbp from the location where a centromere was ancestrally seeded and became inactive in the common ancestor of humans and apes 20-25 million years ago. Long read mapping to the T2T-CHM13 human genome assembly revealed that the insertion derives from a specific region of chromosome 18 centromeric 12-mer HOR array in which the monomer size follows a regular pattern. The rearrangement did not directly disrupt any gene or predicted regulatory element and did not alter the methylation status of the surrounding region, consistent with the absence of phenotypic consequences in the carrier. This case demonstrates a likely rare but new class of structural variation that we name  alpha satellite insertion. It also expands our knowledge on alphoid DNA dynamics and conveys the possibility that alphoid arrays can relocate near vestigial centromeric sites.