Comparative analysis of cattle breeds as satellite cell donors for cultured beef BackgroundCultured meat is a promising new field with the potential for considerable environmental and animal welfare benefits. One technological approach to cultured meat production utilises the proliferative and differentiative capacity of muscle-derived satellite cells (SCs) to produce large volumes of cultured muscle tissue from small biopsies of donor animals. Differing genotypes between cattle breeds lead to predictable phenotypic traits, resulting in breeds being favoured for their respective meat or milk production characteristics in the livestock industry. However, whilst these breeds show significant differences in muscle growth, it is unclear whether the physiological differences observed between them in vivo are reflected in differences in SC behaviour in vitro, particularly with respect to proliferation, differentiation and cellular longevity, and hence whether particular breeds might represent preferred SC donors for a cultured beef bioprocess.

ResultsComparing SCs isolated from five breeds (Belgian Blue, Holstein Friesian, Galloway, Limousin and Simmental), we found that the proliferation rates were largely unaffected by the donor breed. In contrast, potentially meaningful differences were observed in the kinetics and extent of myogenic differentiation. Furthermore, whilst differentiation dropped for all breeds with increasing population doublings (PDs), SCs from Belgian Blue and Limousin cattle showed significantly longer retention of differentiation capacity over long-term passaging.

ConclusionSCs from all breeds were able to proliferate and differentiate, although Limousin and (particularly) Belgian Blue cattle, both breeds commonly used for traditional meat production, may represent preferred donors for cultured beef production.
The H4K20 demethylase DPY-21 regulates thedynamics of condensin DC binding Condensin is a multi-subunit SMC complex that binds to and compacts chromosomes. Here we addressed the regulation of condensin binding dynamics using C. elegans condensin DC, which represses X chromosomes in hermaphrodites for dosage compensation. We established fluorescence recovery after photobleaching (FRAP) using the SMC4 homolog DPY-27 and showed that a well-characterized ATPase mutation abolishes its binding. Next, we performed FRAP in the background of several chromatin modifier mutants that cause varying degrees of X-chromosome derepression. The greatest effect was in a null mutant of the H4K20me2 demethylase DPY-21, where the mobile fraction of condensin DC reduced from [~]30% to 10%. In contrast, a catalytic mutant of dpy-21 did not regulate condensin DC mobility. Hi-C data in the dpy-21 null mutant showed little change compared to wild type, uncoupling Hi-C measured long-range DNA contacts from transcriptional repression of the X chromosomes. Together, our results indicate that DPY-21 has a non-catalytic role in regulating the dynamics of condensin DC binding, which is important for transcription repression.
Common anti-haemostatic medications increase the severity of uropathogenic Escherichia coli sepsis Uropathogenic Escherichia coli (UPEC) causes urinary tract infections that can result in sepsis. The haemostatic system is protective in the pyelonephritis stage of ascending UPEC infection, but the role of the haemostatic system has not been investigated during sepsis. Here we utilize a zebrafish-UPEC systemic infection model to visualize infection-induced coagulation and examine the effects of commonly prescribed anti-haemostatic medications on the infection severity. Treatment of systemically infected zebrafish with warfarin, aspirin, or ticagrelor reduced host survival, while stabilization of clots with aminocaproic acid increased host survival. Anti-haemostatic drug treatment increased UPEC burden. Our findings provide evidence that commonly prescribed anti-haemostatic medications may worsen the outcome of severe UPEC infection. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Stress-dependent cell stiffening by tardigrade tolerance proteins through reversible formation of cytoskeleton-like filamentous network and gel-transition Tardigrades are able to tolerate almost complete dehydration by entering a reversible ametabolic state called anhydrobiosis and resume their animation upon rehydration. Dehydrated tardigrades are exceptionally stable and withstand various physical extremes. Although trehalose and late embryogenesis abundant (LEA) proteins have been extensively studied as potent protectants against dehydration in other anhydrobiotic organisms, tardigrades produce high amounts of tardigrade-unique protective proteins. Cytoplasmic-abundant heat-soluble (CAHS) proteins are uniquely invented in the lineage of eutardigrades, a major class of the phylum Tardigrada and are essential for their anhydrobiotic survival. However, the precise mechanisms of their action in this protective role are not fully understood. In the present study, we first postulated the presence of tolerance proteins that form protective condensates via phase separation in a stress-dependent manner and searched for tardigrade proteins that reversibly form condensates upon dehydration-like stress. Through comprehensive analysis, we identified 336 such proteins, collectively dubbed "Desolvataion-induced ReversiblY condensing Proteins (DRYPs)". Unexpectedly, we rediscovered CAHS proteins as highly enriched in DRYPs, 3 of which were major components of DRYPs. We revealed that these CAHS proteins reversibly polymerize into many cytoskeleton-like filaments depending on hyperosmotic stress in cultured cells and undergo reversible gel-transition in vitro. CAHS filamentation increases cell stiffness to resist deformation and improves resistance to dehydration-like stress. The conserved putative helical C-terminal region is necessary and sufficient for filament formation by CAHS proteins, and mutations disrupting the secondary structure of this region impaired both the filament formation and the gel transition. On the basis of these results, we propose that CAHS proteins are novel cytoskeletal proteins that form filamentous networks and undergo gel-transition in a stress-dependent manner to provide on-demand physical stabilization of cell integrity against deformative forces during dehydration and also contribute to the exceptional physical stability in a dehydrated state. Drug screening at single-organoid resolution via bioprinting and interferometry High-throughput drug screening is an established approach to investigate tumor biology and identify therapeutic leads. Traditional platforms for high-throughput screening use two-dimensional cultures of immortalized cell lines which do not accurately reflect the biology of human tumors. More clinically relevant model systems, such as three-dimensional tumor organoids, can be difficult to screen and scale. For example, manually seeded organoids coupled to destructive endpoint assays allow for the characterization of response to treatment, but do not capture the transitory changes and intra-sample heterogeneity underlying clinically observed resistance to therapy. We therefore developed a pipeline to generate bioprinted tumor organoids linked to label-free, real-time imaging via high-speed live cell interferometry (HSLCI) and machine learning-based quantitation of individual organoids. Bioprinting cells gives rise to 3D organoid structures that preserve tumor histology and gene expression. HSLCI imaging in tandem with machine learning-based image segmentation and organoid classification tools enables accurate, label-free parallel mass measurements for thousands of bioprinted organoids. We demonstrate that our method quantitatively identifies individual organoids as insensitive, transiently sensitive, or persistently sensitive to specific treatments. This opens new avenues for rapid, actionable therapeutic selection using automated tumor organoid screening. 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. Optimized detection of allelic imbalances specific for homologous recombination deficiency improves the prediction of clinical outcomes in cancer Homologous recombination DNA-repair deficiency (HRD) is a common driver of genomic instability and confers a therapeutic vulnerability in human cancer. Using a systematic pan- cancer characterization, we revealed distinct patterns of somatic allelic imbalances (AIs) in high-grade serous ovarian cancer (HGSC). Using machine learning on a multi-omics dataset, we generated an optimized algorithm to detect HRD in HGSC (ovaHRDscar). ovaHRDscar improved the prediction of clinical outcomes in three independent validation cohorts. Characterization of 98 spatiotemporally distinct tumor samples indicated ovary/adnex as the preferred site to assess HRD. Cancer-specific optimization improved the prediction of clinical outcomes also in triple-negative breast cancer (tnbcHRDscar). In conclusion, our systematic approach based on multi-omics data improves the detection of HRD with the premise to improve patient selection for HR-targeted therapies. Determining significant correlation between pairs of extant characters in a small parsimony framework When studying the evolutionary relationships among a set of species, the principle of parsimony states that a relationship involving the fewest number of evolutionary events is likely the correct one. Due to its simplicity, this principle was formalized in the context of computational evolutionary biology decades ago by, e.g., Fitch and Sankoff. Because the parsimony framework does not require a model of evolution, unlike maximum likelihood or Bayesian approaches, it is often a good starting point when no reasonable estimate of such a model is available.

In this work, we devise a method for determining if pairs of discrete characters are significantly correlated across all most parsimonious reconstructions, given a set of species on these characters, and an evolutionary tree. The first step of this method is to use Sankoffs algorithm to compute all most parsimonious assignments of ancestral states (of each character) to the internal nodes of the phylogeny. Correlation between a pair of evolutionary events (e.g., absent to present) for a pair of characters is then determined by the (co-) occurrence patterns between the sets of their respective ancestral assignments. The probability of obtaining a correlation this extreme (or more) under a null hypothesis where the events happen randomly on the evolutionary tree is then used to assess the significance of this correlation. We implement this method: parcours (PARsimonious CO-occURrenceS) and use it to identify significantly correlated evolution among vocalizations and morphological characters in the Felidae family.

The parcours tool is freely available at https://github.com/murraypatterson/parcours Neurological mechanism of sensory deficits after developmental exposure to non-dioxin-like polychlorinated biphenyls (PCBs) in zebrafish larvae The most abundant polychlorinated biphenyl (PCB) congeners found in the environment and in humans are neurotoxic. This is of particular concern for early life stages because the exposure of the more vulnerable developing nervous system to chemicals can result in neurobehavioral disorders. To uncover currently unknown links between PCB target mechanisms and neurobehavioral deficits, we investigated the effects of the non-dioxin-like (NDL) congener PCB153 on neuronal morphology and synaptic transmission linked to the proper execution of a sensorimotor response using zebrafish as a vertebrate model. Zebrafish that were exposed during development to concentrations similar to those found in human cord blood and PCB contaminated sites showed a delay in startle response. A similar delay was observed for other NDL congeners but not for the potent dioxin-like congener PCB126. Morphological and biochemical data demonstrate that while exposure to PCB153 induced swelling of afferent sensory neurons, the disruption of dopaminergic and GABAergic signaling is the dominant mechanism associated with motor movement. The effects on important and broadly conserved signaling mechanisms in vertebrates suggest that NDL PCBs may contribute to neurodevelopmental abnormalities in humans and, with the startle response being critical for the survival of fish, to evolutionary adaptation in wildlife. Estimating the fundamental niche: accounting for the uneven availability of existing environmental conditions in the calibration area Studies that question important conceptual and methodological aspects of the field of ecological niche modeling (and species distribution modeling) have cast doubts on whether it is possible to estimate the fundamental niche of a species using presence-only data. The main limitation in niche estimation is that presence data come from the realized niche, which is only a subset of the fundamental niche. Most of the existing methods lack the ability to overcome this limitation and therefore fit niches that resemble the realized niche. To obtain a more accurate estimate of the fundamental niche, we propose using the geographic region that is accessible to a species (based on its dispersal ability) to determine a sampling distribution, in environmental space, from which we can quantify the likelihood of observing a particular environmental combination in a sample of presence points. We incorporate this sampling distribution into a multivariate normal model (i.e., a Mahalanobis distance model) by creating a weight function that modifies the probabilities of observing different environmental combinations in a sample of presences. This modification accounts for the uneven availability of environmental conditions. We show that the parameters of this weighted-normal model can be approximated with a maximum likelihood estimation approach; and then used to draw ellipsoids (confidence regions) that represent the fundamental niche of the species. We illustrate the application of our model with two worked examples. First, we use presence data for an invasive species and an accessible area that includes only its native range to evaluate whether the fitted model predicts confirmed establishments of the species outside its native range. Second, we use presence data for closely related species with known accessible areas to demonstrate how the different dispersal abilities of the species constrain a classic Mahalanobis distance model. Overall, we show that accounting for the distribution of environmental conditions that are accessible to a species indeed affects the estimation of the ellipsoids used to model its fundamental niche. Sperm Histone H3 Lysine 4 tri-methylation serves as a metabolic sensor of paternal obesity and is associated with the inheritance of metabolic dysfunction ObjectiveParental environmental exposures can strongly influence descendant risks for adult disease. How paternal obesity changes the sperm chromatin leading to the acquisition of metabolic disease in offspring remains controversial and ill-defined. The objective of this study was to assess: (1) whether obesity induced by a high-fat diet alters sperm histone methylation; (2) whether paternal obesity can induce metabolic disturbances across generations; (3) whether there could be cumulative damage to the sperm epigenome leading to enhanced metabolic dysfunction in descendants; and (4) whether obesity-sensitive regions associate with embryonic epigenetic and transcriptomic profiles. Using a genetic mouse model of epigenetic inheritance, we investigated the role of histone H3 lysine 4 methylation (H3K4me3) in the paternal transmission of metabolic dysfunction. This transgenic mouse overexpresses the histone demethylase enzyme KDM1A in the developing germline and has an altered sperm epigenome at the level of histone H3K4 methylation. We hypothesized that challenging transgenic sires with a high-fat diet would further erode the sperm epigenome and lead to enhanced metabolic disturbances in the next generations.

MethodsTo assess whether paternal obesity can have inter- or transgenerational impacts, and if so, to identify potential mechanisms of this non-genetic inheritance, we used wildtype C57BL/6NCrl and transgenic males with a pre-existing altered sperm epigenome. To induce obesity, sires were fed either a control or high-fat diet (10% or 60% kcal fat, respectively) for 10-12 weeks, then bred to wildtype C57BL/6NCrl female fed a regular diet. F1 and F2 descendants were characterized for metabolic phenotypes by examining the effects of paternal obesity by sex, on body weight, fat mass distribution, the liver transcriptome, intraperitoneal glucose and insulin tolerance tests. To determine whether obesity altered the F0 sperm chromatin, native chromatin immunoprecipitation-sequencing targeting H3K4me3 was performed. To gain insight into mechanisms of paternal transmission, we compared our sperm H3K4me3 profiles with embryonic and placental chromatin states, histone modification and gene expression profiles.

ResultsObesity-induced alterations in H3K4me3 occurred at genes implicated in metabolic, inflammatory, and developmental processes. These processes were associated with offspring metabolic dysfunction and corresponded to genes enriched for H3K4me3 in embryos, and overlapped embryonic and placenta gene expression profiles. Transgenerational susceptibility to metabolic disease was only observed when obese F0 had a pre-existing modified sperm epigenome. This coincided with increased H3K4me3 alterations in sperm and more severe phenotypes affecting their offspring.

ConclusionsOur data suggest sperm H3K4me3 might serve as a metabolic sensor that connects paternal diet with offspring phenotypes via the placenta. This non-DNA based knowledge of inheritance has the potential to improve our understanding of how environment shapes heritability and may lead to novel routes for the prevention of disease. This study highlights the need to further study the connection between the sperm epigenome, placental development and childrens health. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. SufB intein splicing in Mycobacterium tuberculosis is influenced by two remote conserved N-extein Histidines Inteins are auto-processing domains that implement a multi-step biochemical reaction termed protein splicing, marked by cleavage and formation of peptide bonds. They excise from a precursor protein, generating a functional protein via covalent bonding of flanking exteins. We report the kinetic study of splicing and cleavage reaction in [Fe-S] cluster assembly protein SufB from Mycobacterium tuberculosis. Although it follows a canonical intein splicing pathway, distinct features are added by extein residues present in the active site. Sequence analysis identified two conserved histidines in the N-extein region; His-5 and His-38. Kinetic analyses of His-5Ala and His-38Ala SufB mutants exhibited significant reductions in splicing and cleavage rates relative to the SufB wild-type precursor protein. Structural analysis and molecular dynamics simulations suggested that Mtu SufB displays a unique mechanism where two remote histidines work concurrently to facilitate N-terminal cleavage reaction. His-38 is stabilized by the solvent-exposed His-5, and can impact N-S acyl shift by direct interaction with the catalytic Cys1. Development of inteins as biotechnological tools or as pathogen specific novel antimicrobial targets requires a more complete understanding of such unexpected roles of conserved extein residues in protein splicing. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. An evolutionary ecological mechanism of CyanoHABs: unavoidable consequences when low nutrient requirements meet eutrophication The mechanism of cyanobacterial harmful algal blooms (CyanoHABs) is complicated and confusing. One major reason is they are studied primarily from an ecological perspective and on bloom-forming species only. This narrow angle loses a broader evolutionary and ecological context in which CyanoHABs occur and fails to provide information on relevant components to achieve a wholistic understanding.

To derive a comprehensive mechanism of CyanoHABs, we examine CyanoHABs through the overlooked evolutionary and ecological lenses: evolutionary radiation, ecological comparison with co-living algae, and recently identified genomic functional repertoire between blooming and non-blooming species. We found four key factors: (1) elaborate diverse functional repertoire and low nutrient requirement in cyanobacteria molded by early adaptive evolution, (2) cyanobacteria having lower nutrient requirements than co-living algae, (3) there is no directed evolution in biological functions toward water eutrophication in cyanobacteria, (4) the CyanoHAB-associated functional repertoire are more abundant and complete in blooming than non-blooming species. These factors lead us to postulate a preliminary mechanism of CyanoHABs as a synergistic quad: superior functional repertoire, established with long adaptive radiation under nutrient-deficient conditions and not evolved toward eutrophic conditions, enables cyanobacteria to efficiently utilize elevated nutrients under current eutrophic levels for excess growth and CyanoHABs thereof, due to their lower nutrient requirements than co-living algae. This working comprehensive mechanism without doubt needs further empirical testing, which can be undertaken with more comparative studies of multiple species using integrated systems biology approaches. Brain expressed FKBP5 delineates a therapeutic subtype of severe mental illness Deducing genes capable of classifying biologically-distinct psychiatric subtypes, and their targets for treatment, is a priority approach in psychiatry. FKBP5 is one such gene with strong evidence of utility to delineate a trans-diagnostic psychiatric subtype. Yet how brain-expressed FKBP5 is affected in psychiatric disorders in humans is not fully understood and critical for propelling FKBP5-targeting treatment development. We performed a large-scale postmortem study (n=895) of FKBP5 using dorsolateral prefrontal cortex samples derived from individuals with severe psychiatric disorders with a comprehensive battery of bulk/single-cell omics and histological analyses. We observed consistently heightened FKBP5 mRNA and protein in psychopathology, moderated by genotype and age, and accompanied by DNA methylation changes in key enhancers. These effects were most prominent in superficial-layer pyramidal cells. Heightened FKBP5 was also differentially associated with downstream pathways according to age, being specifically associated with synaptic transmission in early adulthood, and neuroinflammation and neurodegeneration in later life. The renal inflammatory network of nephronophthisis STRUCTURED ABSTRACTO_ST_ABSBACKGROUNDC_ST_ABSThe majority of genetic kidney disease leading to kidney failure is caused by mutations in ciliary genes. How cilia malfunction leads to progressive kidney damage is poorly understood, but recent evidence links ciliopathy genes to CCL2 dependent macrophage recruitment in autosomal dominant polycystic kidney disease (ADPKD), the most studied renal ciliopathy. Whether or not renal inflammation is involved in other renal ciliopathies is unclear.

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

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

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

SIGNIFICANCE STATEMENTMutations in genes encoding primary cilia proteins are the leading cause of genetic kidney failure. In autosomal dominant polycystic kidney disease (ADPKD), deregulated cilia signaling leads to kidney infiltration by macrophages through the chemokine CCL2. Little is known about renal inflammation in nephronophthisis (NPH), the most frequent pediatric renal ciliopathy. Using NPH mice models, tissues and cells from NPH patients, we unveil renal inflammation as preeminent feature of NPH. Remarkably, the renal inflammatory evoked by ciliary gene mutations in NPH does not overlap with ADPKD: it is CCL2 independent, involves a prominent recruitment of neutrophils and T cells and a specific cytokine signature. This unforeseen findings strengthen the link between primary cilia and renal inflammation. Simpati: patient classifier identifies signature pathways based on similarity networks for the disease prediction BACKGROUNDPathway-based patient classification is a supervised learning task which implies a model learning pathways as features to predict the classes of patients. The counterpart of enrichment tools for the pathway analysis are fundamental methods for clinicians and biomedical scientists. They allow to find signature cellular functions which help to define and annotate a disease phenotype. They provide results which lead human experts to manually classify patients. It is a paradox that pathwaybased classifiers which natively resolve this objective are not strongly developed. They could simulate the human way of thinking, decipher hidden multivariate relationships between the deregulated pathways and the disease phenotype, and provide more information than a probability value. Instead, there are currently only two classifiers of such kind, they require a nontrivial hyperparameter tuning, are difficult to interpret and lack in providing new insights. There is the need of new classifiers which can provide novel perspectives about pathways, be easy to apply with different biological omics and produce new data enabling a further analysis of the patients.

RESULTSWe propose Simpati, an innovative and interpretable patient classifier based on pathway-specific patient similarity networks. The first classifier to adopt ad-hoc novel algorithms for such graph type. It standardizes the biological high-throughput dataset of patients profiles with a propagation algorithm that considers the interconnected nature of the cells molecules for inferring a new activity score. This allows Simpati to classify with dense, sparse, and non-homogenous omic data. Simpati organizes patients molecules in pathways represented by patient similarity networks for being interpretable, handling missing data and preserving the patient privacy. A network represents patients as nodes and a novel similarity measure determines how much every pair act co-ordinately in a pathway. Simpati detects signature biological processes based on how much the topological properties of the related networks separate the patient classes. In this step, it includes a new cohesive subgroup detection algorithm to handle patients not showing the same pathway activity as the other class members. An unknown patient is then classified by a unique recommender system which considers how much is similar to known patients and distant from being an outlier. Simpati outperforms previously published classifiers on five cancer datasets described with two biological omics, classifies well with sparse data, identifies more relevant pathways associated to the patients disease than the competitors and has the lowest computational requirements.

CONCLUSIONSimpati can serve as generic-purpose pathway-based classifier of patient classes. It provides signature pathways to unveil the altered biological mechanisms of a disease phenotype and to classify patients according to the learnt pathway-specific similarities. The signature condition and patient prediction can be deciphered considering the patient similarity networks which must reveal the members of a patient class more cohesive and similar than the non-members. Simpati divides the pathways in up and downinvolved. Upinvolved when the signaling cascades generated by the altered molecules of the disease patients impact stronger the pathway than the ones of the control class. We provide an R implementation, a graphical user interface and a visualization function for the patient similarity networks. The software is available at: https://github.com/LucaGiudice/Simpati Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

MethodsWe propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F -penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

ResultsWe evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F -penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach. Once-daily feeding is associated with better health in companion dogs: Results from the Dog Aging Project A variety of diets have been studied for possible anti-aging effects. In particular, studies of intermittent fasting and time-restricted feeding in laboratory rodents have found evidence of beneficial health outcomes. Companion dogs represent a unique opportunity to study diet in a large mammal that shares human environments. The Dog Aging Project has been collecting data on thousands of companion dogs of all different ages, sizes, and breeds since 2019. We leveraged this diverse cross-sectional dataset to investigate associations between feeding frequency and cognitive function (n = 10,474) as well as nine broad categories of health conditions (n = 24,238). Controlling for sex, age, breed, and other potential confounders, we found that dogs fed once daily rather than more frequently had lower mean scores on a cognitive dysfunction scale, and lower odds of having gastrointestinal, dental, orthopedic, kidney/urinary, and liver/pancreas disorders. Therefore, we find that once-daily feeding is associated with better health in multiple domains. Future research with longitudinal data can provide stronger evidence for a possible causal effect of feeding frequency on health in companion dogs. 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. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Genetic variation at mouse and human ribosomal DNA influences associated epigenetic states BackgroundRibosomal DNA (rDNA) displays substantial inter-individual genetic variation in human and mouse. A systematic analysis of how this variation impacts epigenetic states and expression of the rDNA has thus far not been performed.

ResultsUsing a combination of long- and short-read sequencing, we establish that 45S rDNA units in the C57BL/6J mouse strain exist as distinct genetic haplotypes that influence the epigenetic state and transcriptional output of any given unit. DNA methylation dynamics at these haplotypes are dichotomous and life-stage specific: at one haplotype, the DNA methylation state is sensitive to the in utero environment, but refractory to post-weaning influences, whereas other haplotypes entropically gain DNA methylation during ageing only. On the other hand, individual rDNA units in human show limited evidence of genetic haplotypes, and hence little discernible correlation between genetic and epigenetic states. However, in both species, adjacent units show similar epigenetic profiles, and the overall epigenetic state at rDNA is strongly positively correlated with total rDNA copy number. Analysis of different mouse inbred strains reveals that in some strains, such as 129S1/SvImJ, rDNA copy number is only approximately 150 copies per diploid genome and DNA methylation levels are <5%.

ConclusionsOur work demonstrates that rDNA-associated genetic variation has a considerable influence on rDNA epigenetic state and consequently rRNA expression outcomes. In the future, it will be important to consider the impact of inter-individual rDNA (epi)genetic variation on mammalian phenotypes and diseases. Heat-Dependent Hairpin Melting Drives Thermo-gated TRPV1 Opening The capsaicin receptor TRPV1 is a bio-thermometer in a primary afferent sensory neuron for noxious heat detection. However, the underlying molecular mechanism is unclear. Here, a hairpin thermodynamic model, together with graph theory, was to examine a role of temperature-dependent hairpin sizes in controlling the temperature threshold and sensitivity and gating states of TRPV1. The results showed that heat-dependent melting of the largest hairpin initiates a temperature threshold for channel opening while the smaller hairpins stabilize the open state and increase the heat efficacy or sensitivity. Therefore, this model may provide a broad structural basis for the thermo-gated TRP channels. The OPAQUE1/DISCORDIA2 myosin XI is required for phragmoplast guidance during asymmetric cell division in maize Asymmetric divisions produce cells with different fates and are critical for development. Here we show that a maize myosin XI protein, OPAQUE1 (O1), is necessary for asymmetric divisions during maize stomatal development. We analyzed stomatal precursor cells prior to and during asymmetric division to determine why o1 mutants have abnormal division planes. Cell polarization and nuclear positioning occur normally in the o1 mutant, and the future site of division is correctly specified. The defect in o1 occurs during late cytokinesis, when the plant-specific phragmoplast -- made of microtubules, actin and other proteins - forms the nascent cell plate. The phragmoplast becomes misguided and does not meet the previously established division site. Initial phragmoplast guidance is correct in o1. However, as phragmoplast expansion continues, phragmoplasts in o1 stomatal precursor cells become misguided and do not meet the cortex at the established division site. To understand how this myosin protein contributes to phragmoplast guidance, we identified O1-interacting proteins. Other myosins, specific actin-binding proteins, and maize kinesins related to the Arabidopsis thaliana division site markers PHRAGMOPLAST ORIENTING KINESINs (POKs) interact with O1. We propose that different myosins are important at multiple steps of phragmoplast expansion, and the O1 actin motor and POK-like microtubule motors work together to ensure correct late-stage phragmoplast guidance.