Micro topographical instruction of bacterial attachment, biofilm formation and in vivo host response Bio-instructive materials that prevent bacterial biofilm formation and drive an appropriate host immune response have the potential to significantly reduce the burden of medical device-associated infections. As bacterial surface attachment is sensitive to surface topography, we experimentally surveyed 2,176 combinatorially generated shapes using an unbiased high-throughput micro topographical polystyrene polymer chip screen. This identified topographies that reduced colonization in vitro by up to 15-fold compared with a flat surface for both motile and non-motile bacterial pathogens. Equivalent reductions were achieved on polyurethane, a polymer commonly used in medical devices. Using machine learning methods, a set of design rules based on generalisable topo-descriptors was established for predicting bacteria-resistant micro topographies. In a murine foreign body infection model, anti-attachment topographies were shown to be refractory to P. aeruginosa colonization and serendipitously, the fibrotic response to the implant was reduced, offering control of foreign body response using simple topographical patterning of non-eluting implants.
Optical regulation of endogenous RhoA reveals selection of cellular responses by signal amplitude How protein signaling networks respond to different input strengths is an important but poorly understood problem in cell biology. For example, the small GTPase RhoA regulates both focal adhesion (FA) growth or disassembly, but whether RhoA serves as a switch selecting between cellular outcomes, or if outcomes are simply modulated by additional factors in the cell, is not clear. Here, we develop a photoswitchable RhoA guanine exchange factor, psRhoGEF, to precisely control endogenous RhoA activity. We also develop a FRET-based biosensor to allow visualization of RhoA activity together with psRhoGEF control. Using these new optical tools, we discover that low levels of RhoA activation preferentially induce FA disassembly in a Src-dependent manner, while high levels induce both FA growth and disassembly in a ROCK-dependent manner. Thus, rheostatic control of RhoA activation with photoswitchable RhoGEF reveals that cells can use signal amplitude to produce multiple responses to a single biochemical signal.
Host species identity shapes the diversity and structure of insect microbiota As for most of the life that inhabits our planet, microorganisms play an essential role in the fitness of insects, including nutrition, reproduction, defence, and many other functions. More recently, we assisted in an exponential growth of studies describing the taxonomical composition of bacterial communities across insects phylogeny. However, there is still an outstanding question that needs to be answered: which factors contribute most to shape insects microbiomes? This study tries to find an answer to this question by taking advantage of publicly available sequencing data and reanalysing over 4,000 samples of insect-associated bacterial communities under a common framework. Results suggest that insect taxonomy has a wider impact on the structure and diversity of their associated microbial communities than the other factors considered (diet, sex, life stage, sample origin, and treatment). However, when specifically testing for signatures of co-diversification of insect species and their microbiota, analyses found weak support for this, suggesting that while insect species strongly drive the structure and diversity of insect microbiota, the diversification of those microbial communities did not follow their hosts phylogeny. Furthermore, a parallel survey of the literature highlights several methodological limitations that need to be considered in future research endeavours. Higher-order epistasis and phenotypic prediction Contemporary high-throughput mutagenesis experiments are providing an increasingly detailed view of the complex patterns of genetic interaction that occur between multiple mutations within a single protein or regulatory element. By simultaneously measuring the effects of thousands of combinations of mutations, these experiments have revealed that the genotype-phenotype relationship typically reflects genetic interactions not only between pairs of sites, but also higher-order interactions among larger numbers of sites. However, modeling and understanding these higher-order interactions remains challenging. Here, we present a method for reconstructing sequence-to-function mappings from partially observed data that can accommodate all orders of genetic interaction. The main idea is to make predictions for unobserved genotypes that match the type and extent of epistasis found in the observed data. This information on the type and extent of epistasis can be extracted by considering how phenotypic correlations change as a function of mutational distance, which is equivalent to estimating the fraction of phenotypic variance due to each order of genetic interaction (additive, pairwise, three-way, etc.). Using these estimated variance components, we then define an empirical Bayes prior that in expectation matches the observed pattern of epistasis, and reconstruct the genotype-phenotype mapping by conducting Gaussian process regression under this prior. To demonstrate the power of this approach, we present an application to the antibody-binding domain GB1 and also provide a detailed exploration of a dataset consisting of high-throughput measurements for the splicing efficiency of human pre-mRNA 5' splice sites, for which we also validate our model predictions via additional low-throughput experiments. Distribution and molecular evolution of the anti-CRISPR family AcrIF7 Anti-CRISPRs are proteins capable of blocking CRISPR-Cas systems and typically encoded in mobile genetic elements. Since their discovery, numerous anti-CRISPR families have been identified. However, little is known about the distribution and sequence diversity of members within a family, nor how these traits influence the anti-CRISPRs function and evolution. Here we use AcrIF7 to explore the dissemination and molecular evolution of an anti-CRISPR family. We uncovered five sub-clusters and prevalent anti-CRISPR variants within the group. Remarkably, AcrIF7 homologs display high similarity despite their broad geographical, ecological and temporal distribution. Although mainly associated with Pseudomonas aeruginosa, AcrIF7 was identified in distinct genetic backgrounds indicating horizontal dissemination, primarily by phages. Using mutagenesis, we recreated variation observed in databases but also extended the sequence diversity of the group. Characterisation of the variants identified residues key for the anti-CRISPR function and other contributing to its mutational tolerance. Moreover, molecular docking revealed that variants with affected function lose key interactions with its CRISPR-Cas target. Analysis of publicly available data and the generated variants suggests that the dominant AcrIF7 variant corresponds to the minimal and optimal anti-CRISPR selected in the family. Our study provides a blueprint to investigate the molecular evolution of anti-CRISPR families. Correlation between speed and turning naturally arises for sparsely sampled cell movements Mechanisms regulating cell movement are not fully understood. One feature of cell movement that determines how far cells displace from an initial position is persistence, the ability to perform movements in a direction similar to the previous movement direction. Persistence is thus determined by turning angles between two sequential displacements and be characterized by an average turning angle or persistence time. Recent studies found that a cells average speed and turning are negatively correlated, suggesting a fundamental cell-intrinsic program whereby cells with a lower turning ability (i.e., larger persistence time) are intrinsically faster (or faster cells turn less). By simulating correlated or persistent random walks (PRWs) using two different frameworks (one based on von Mises-Fisher (vMF) distribution and another based on Ornstein-Uhlenbeck (OU) process) we show that the negative correlation between speed and turning naturally arises when cell trajectories are sub-sampled, i.e., when the frequency of sampling is lower than frequency at which cells make movements. This effect is strongest when the sampling frequency is on the order of magnitude with the typical cell persistence time and when cells vary in persistence time. Both conditions are observed for datasets of T cell movements in vivo that we have analyzed. In simulations the correlation arises due to randomness of cell movements resulting in highly variable persistence times for individual cells that with sub-sampling leads to large variability of average cell speeds. Interestingly, previously suggested methodology of calculating displacement of cell cohorts with different speeds resulted in similar results whether or not there is a cell-intrinsic correlation between cell speed and persistence. For both vMF- and OU-based simulations of PRWs we could find parameter values (distribution of persistence times, speeds, and sampling frequency) that matched experimentally measured correlations between speed and turning for two datasets of T cell movement in vivo suggesting that such simple correlations are not fully informative on the intrinsic link between speed and persistence. Our results thus suggest that sub-sampling may contribute to (and perhaps fully explains) the observed correlation between speed and turning at least for some cell trajectory data and emphasize the role of sampling frequency in inference of critical cellular parameters of cell motility such as speeds.

Secondary AbstractMeasurement of cell movements often results in a negative correlation between average speed and average turning angle suggesting an existence of a universal, cell-intrinsic movement program. We show that such a negative correlation may arise if cells in the population differ in their ability for persistent movement when the movement data are sub-sampled. We show that sub-sampling of cell trajectories generated using two different frameworks of persistent random walk can match the experimentally observed correlation between speed and turning for T cell movements in vivo. Determinants of Murine Skin Microbiota Composition in Homeostasis and Wound Healing Animal microbiota have complex interactions with hosts and environment that determines its composition. Yet the ability of hosts to determine their microbiota composition is less well studied. In this study, to investigate the role host genetics in determining skin microbiota, we used 30 different mouse strains from the recombinant inbred panel, the Collaborative Cross. Murine skin microbiota composition was strongly dependent on murine strain with > 50% of the variation explained by murine strain. In particular, a quantitative trait locus on chromosome 4 associates both with Staphylococcus abundance and principal-component multi-trait analyses. Additionally, excisional wound associated changes in microbiota composition were not uniform across mouse strains and were host-specific, the genetic background accounting for about 40% of the variation in microbiota. Genetic background also had the highest effect on the healing speed of wounds accounting for over 50% of the variation while mouse age and microbiota composition change accounted only for 20% and 5% of the healing speed despite reaching statistical significance. In conclusion, host genetics has a significant impact on the skin microbiota composition during both homeostasis and wound healing. These findings have long reaching implications in our understanding of associations between microbiota dysbiosis and disease. Long term analysis of social structure: evidence of age-based consistent associations in male Alpine ibex Despite its recognized importance for understanding the evolution of animal sociality as well as for conservation, long term analysis of social networks of animal populations is still relatively uncommon. We investigated social network dynamics in males of a gregarious mountain ungulate (Alpine ibex, Capra ibex) over ten years focusing on groups, sub-groups and individuals, exploring the dynamics of sociality over different scales. Despite the social structure changing between seasons, the Alpine ibex population was highly cohesive: fission-fusion dynamics lead almost every male in the population to associate with each other male at least once. Nevertheless, we found that male Alpine ibex showed preferential associations that were maintained across seasons and years. Age seemed to be the most important factor driving preferential associations while other characteristics, such as social status, appeared less crucial. We also found that centrality measures were influenced by age and were also related to individual physical condition. The multi-scale and long-term frame of our study helped us show that ecological constrains, such as resource availability, may play a role in shaping associations in a gregarious species, but they cannot solely explain sociality and preferential association that are likely also to be driven by life-history linked physiological and social needs. Our results highlight the importance of long-term studies based on individually recognizable subjects to help us build on our understanding of the evolution of animal sociality. Selective Src-Family B Kinase Inhibition Promotes Pulmonary Artery Endothelial Cell Dysfunction. Protein tyrosine kinases (PTKs) are essential for eukaryotic signaling. By targeting select PTKs, the group of drugs known as Tyrosine kinase inhibitors (TKIs) have proven to be effective for treating multiple diseases ranging from cancer to pulmonary fibrosis. However, some TKIs also paradoxically lead to the development of adverse conditions such as pulmonary arterial hypertension (PAH) by promoting endothelial cell dysfunction (ECD). We hypothesize that (1) subsets of PTKs may disproportionately modulate signaling pathways critical for endothelial homeostasis, such as BMPR2 signaling, and (2) inhibiting those pro-endothelial PTKs can promote the development of ECD. Herein we use an agnostic high-throughput siRNA screen to investigate how PTKs affect the canonical BMPR2 signaling pathway. Our major finding is that within the Src-family of non-receptor PTKs, the Src-B family promotes canonical BMPR2 signaling while the Src-A family suppresses it. We focus on two representative members of each family, Lck (for Src-B) and Fyn (for Src-A) that are the strongest activators or inhibitors of BMPR2 signaling in the screen. We confirm that Lck is expressed in the endothelium of pulmonary arteries and show that Lck knockout (termed si-Lck) in pulmonary artery endothelial cells (PAECs) suppresses canonical BMPR2 signaling while Fyn knockout (termed si-Fyn) promotes canonical BMPR2 signaling. Furthermore, Lck and Fyn are responsible for opposing functional behaviors in PAECs: si-Lck promotes apoptosis and interferes with tube formation while si-Fyn suppresses apoptosis and promotes tube formation. After analyzing the whole-transcriptome signature of si-Lck and si-Fyn PAECs we find that in addition to BMPR2 signaling suppression, si-Lck (and not si-Fyn) increases a broad number of ECD markers and increases canonical NF-{kappa}{beta} signaling. In summary, for the first time we show that Src-A and B Family of PTKs exert differential control over key endothelial cell signaling pathways resulting in direct phenotypic consequences. This knowledge may help to guide the design of more precise TKIs which avoid adverse drug reactions brought about through endothelial cell dysfunction. caliPER: A software for blood-free parametric Patlak mapping using PET/MRI input function. Routine clinical use of absolute PET quantification techniques is limited by the need for serial arterial blood sampling for input function and more importantly by the lack of automated pharmacokinetic analysis tools that can be readily implemented in clinic with minimal effort. PET/MRI provides the ability for absolute quantification of PET probes without the need for serial arterial blood sampling using image-derived input functions (IDIFs). Here we introduce caliPER, a modular and scalable software for simplified pharmacokinetic modelling of PET probes with irreversible uptake or binding based on PET/MR IDIFs and Patlak Plot analysis. caliPER generates regional values or parametric maps of net influx rate (Ki) using reconstructed dynamic PET images and anatomical MRI aligned to PET for IDIF vessel delineation. We evaluated the performance of caliPER for blood-free region-based and pixel-wise Patlak analyses of [18F] FDG by comparing caliPER IDIF to serial arterial blood input functions and its application in imaging brain glucose hypometabolism in Frontotemporal dementia. IDIFs corrected for partial volume errors including spill-out and spill-in effects were similar to arterial blood input functions with a general bias of around 6-8%, even for arteries <5 mm. The Ki and cerebral metabolic rate of glucose estimated using caliPER IDIF were similar to estimates using arterial blood sampling (<2%) and within limits of whole brain values reported in literature. Overall, caliPER is a promising tool for irreversible PET tracer quantification and can simplify the ability to perform parametric analysis in clinical settings without the need for blood sampling.

HighlightsO_LIcaliPER is an adaptable image processing software for extracting image-derived input functions and generating parametric maps of irreversible PET tracer uptake.
C_LIO_LIAnatomical (T1-weighted/T2-weighted/time-of-flight angiography) MRI carefully aligned to PET provides a robust approach for delineation of vessels on PET, eliminating the need for serial blood sampling for input functions.
C_LIO_LIApplication of caliPER in modelling glucose uptake in patients with Frontotemporal dementia, demonstrates the feasibility of absolute quantification of cerebral metabolic rate of glucose in clinical populations.
C_LI Glutathione regulates transcriptional activation of iron transporters via S-nitrosylation of bHLH factors to modulate subcellular iron homeostasis While glutathione (GSH) is known to regulate iron (Fe) deficiency response in plants, its involvement in modulating subcellular Fe homeostasis remained elusive. In this study, we report that the GSH depleted mutants, cad2-1 and pad2-1 displayed increased sensitivity to Fe deficiency with significant down-regulation of the vacuolar Fe exporters, AtNRAMP3 and AtNRAMP4 and the chloroplast Fe importer, AtPIC1. Moreover, the pad2-1 mutant accumulated higher Fe content in vacuoles and lower in chloroplasts compared with Col-0 under Fe limited condition. Exogenous GSH treatment could enhance the chloroplast Fe content in Col-0 plants but failed to do so in the nramp3nramp4 double mutant demonstrating the role of GSH in modulating subcellular Fe homeostasis. Pharmacological experiments, mutant analysis and promoter assay revealed that this regulation involved transcriptional activation of the transporter genes by a GSH-GSNO module. The Fe responsive bHLH transcription factors (TFs), AtbHLH29, AtbHLH38 and AtbHLH101 were identified to interact with the promoters of these genes which were in turn activated via S-nitrosylation. Together, the present study delineates the role of GSH-GSNO module in regulating subcellular Fe homeostasis by transcriptional activation of the Fe transporters AtNRAMP3, AtNRAMP4 and AtPIC1 via S-nitrosylation of the bHLH TFs during Fe deficiency.

Summary statementGlutathione regulates subcellular iron homeostasis under iron deficiency via GSNO dependent transcriptional activation of AtNRAMP3, AtNRAMP4 and AtPIC1 genes by S-nitrosylation of the iron responsive bHLH transcription factors, AtbHLH29, AtbHLH38 and AtbHLH101. Liver Endothelium Promotes HER3-mediated Cell Survival in KRAS Wild-type and Mutant Colorectal Cancer Cells We previously showed that liver endothelial cells (ECs) secreted soluble factors in a paracrine fashion and activated human epidermal growth factor receptor 3 (HER3, also known as ERBB3) for promoting colorectal cancer (CRC) growth and chemoresistance. However, RAS proteins play a critical role in receptor tyrosine kinase signaling pathways, and KRAS mutations mediate CRC resistance to therapies targeting EGFR, another HER protein. Therefore, the role of KRAS mutation status in EC-induced HER3 activation and CRC survival was investigated as it has therapeutic implications. We used CRC cell lines and patient-derived xenografts harboring KRAS wild-type or mutant genes and demonstrated that liver EC-secreted factors promoted HER3-mediated CRC cell growth independent of KRAS mutation status. Also, blocking HER3 in CRC cells by siRNAs or a HER3 antibody seribantumab blocked EC-induced CRC survival. Our findings highlight the potential of utilizing HER3 targeted therapies for treating patients with mCRC independent of RAS mutational status. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Directed evolution of angiotensin-converting enzyme 2 (ACE2) peptidase activity profiles for therapeutic applications Angiotensin Converting-Enzyme 2 (ACE2) is currently being investigated for its ability to beneficially modulate the Angiotensin receptor (ATR) therapeutic axis to treat multiple human diseases, but its broad substrate scope and diverse physiological roles limit its potential as a therapeutic agent. In this work we engineer ACE2 variants with enhanced hydrolytic activity and specificity toward Angiotensin-II (Ang-II) that may lead to new therapeutic candidates with improved efficacy and reduced off-target side effects. We established a yeast display-based liquid chromatography screen that enabled use of directed evolution to discover ACE2 variants with improved Ang-II activity and specificity relative to the off-target peptide substrates Apelin-13 and Angiotensin-I (Ang-I). We screened ACE2 active site libraries to reveal three substitution-tolerant positions (M360, T371 and Y510) that can be mutated to enhance ACE2s activity profile and followed up on these hits with focused double mutant libraries to further improve the enzyme. Relative to wildtype ACE2, our top variant (T371L/Y510Ile) displayed a sevenfold increase in Ang-II turnover number (kcat), a sixfold diminished catalytic efficiency (kcat/Km) on Apelin-13, and an overall decreased activity on other ACE2 substrates that were not directly assayed in directed evolution screen. At physiologically relevant substrate concentrations, T371L/Y510Ile hydrolyzes more Ang-II than wildtype ACE2 with concomitant Ang-II:Apelin-13 specificity improvements reaching 30-fold. Our efforts have delivered ATR axis-acting therapeutic candidates with relevance to both established and unexplored ACE2 therapeutic applications and provide a foundation for further ACE2 engineering efforts. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. First detection of herpesvirus and frequency of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Multiparametric profiling of cellular lipid metabolism in hypercholesterolemia Systematic insight into cellular dysfunctions can improve understanding of disease etiology, risk assessment and patient stratification. We present a multiparametric high-content imaging platform enabling quantification of low-density lipoprotein (LDL) uptake and lipid storage in cytoplasmic droplets of primary leukocyte subpopulations. We validated this platform with samples from 65 individuals with variable blood LDL-cholesterol (LDL-c) levels, including familial hypercholesterolemia (FH) and non-FH subjects. We integrated lipid storage data into a novel readout, lipid mobilization, measuring the efficiency with which cells deplete lipid reservoirs. Lipid mobilization correlated positively with LDL uptake and negatively with hypercholesterolemia and age, improving differentiation of individuals with normal and elevated LDL-c. Moreover, combination of cell-based readouts with a polygenic risk score for LDL-c explained hypercholesterolemia better than the genetic risk score alone. This platform provides functional insights into cellular lipid trafficking from a few mls of blood and is applicable to dissect metabolic disorders, such as hypercholesterolemia.

MotivationWe have limited information on how cellular lipid uptake and processing differ between individuals and influence the development of metabolic diseases, such as hypercholesterolemia. Available assays are labor intensive, require skilled personnel and are difficult to scale to higher throughput, making it challenging to obtain systematic functional cell-based data from individuals. To overcome this problem, we established a scalable automated analysis pipeline enabling reliable quantification of multiple cellular readouts, including lipid uptake, storage and mobilization, from different white blood cell populations. This approach provides new personalized insights into the cellular basis of hypercholesterolemia and obesity.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/440471v4_ufig1.gif" ALT="Figure 1">
View larger version (84K):
org.highwire.dtl.DTLVardef@9222c1org.highwire.dtl.DTLVardef@27e289org.highwire.dtl.DTLVardef@89ba82org.highwire.dtl.DTLVardef@33c9cd_HPS_FORMAT_FIGEXP  M_FIG C_FIG The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here, we assembled 26 haplotype-resolved genome assemblies from 13 genetically diverse sheep using PacBio HiFi sequencing. We constructed a graph-based ovine pan-genome and discovered 142,422 biallelic insertions and deletions, 7,028 divergent alleles and 13,419 multiallelic variations. We then used a graph-based approach to genotype the biallelic SVs in 684 individuals from 45 domestic breeds and two wild species. Integration with RNA-seq data allows to identify candidate expression-associated SVs. We demonstrate a direct link of SVs and phenotypes by localizing the putative causative insertion in HOXB13 gene responsible for the long-tail trait and identifying multiple large SVs associated with the fat-tail. Beyond generating a benchmark resource for ovine structural variants, our study highlights that animal genetic research will greatly benefit from using a pan-genome graph rather than a single reference genome. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. Cell-type specific circadian bioluminescence rhythms in Dbp reporter mice Circadian rhythms are endogenously generated physiological and molecular rhythms with a cycle length of about 24 h. Bioluminescent reporters have been exceptionally useful for studying circadian rhythms in numerous species. Here, we report development of a reporter mouse generated by modification of a widely expressed and highly rhythmic gene encoding D-site albumin promoter binding protein (Dbp). In this line of mice, firefly luciferase is expressed from the Dbp locus in a Cre-recombinase-dependent manner, allowing assessment of bioluminescence rhythms in specific cellular populations. A mouse line in which luciferase expression was Cre-independent was also generated. The Dbp reporter alleles do not alter Dbp gene expression rhythms in liver or circadian locomotor activity rhythms. In vivo and ex vivo studies show the utility of the reporter alleles for monitoring rhythmicity. Our studies reveal cell-type specific characteristics of rhythms among neuronal populations within the suprachiasmatic nuclei ex vivo. In vivo studies show Dbp-driven bioluminescence rhythms in the liver of Albumin-Cre;DbpKI/+ "liver reporter" mice. After a shift of the lighting schedule, locomotor activity achieved the proper phase relationship with the new lighting cycle more rapidly than hepatic bioluminescence did. As previously shown, restricting food access to the daytime altered the phase of hepatic rhythmicity. Our model allowed assessment of the rate of recovery from misalignment once animals were provided with food ad libitum. These studies confirm the previously demonstrated circadian misalignment following environmental perturbations and reveal the utility of this model for minimally invasive, longitudinal monitoring of rhythmicity from specific mouse tissues. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. Strength of interactions in the Notch gene regulatory network determines patterning and fate in the notochord Development of multicellular organisms requires the generation of gene expression patterns that determines cell fate and organ shape. Groups of genetic interactions known as Gene Regulatory Networks (GRNs) play a key role in the generation of such patterns. However, how the topology and parameters of GRNs determine patterning in vivo remains unclear due to the complexity of most experimental systems. To address this, we use the zebrafish notochord, an organ where coin-shaped precursor cells are initially arranged in a simple unidimensional geometry. These cells then differentiate into vacuolated and sheath cells. Using newly developed transgenic tools together with in vivo imaging, we identify jag1a and her6/her9 as the main components of a Notch GRN that generates a lateral inhibition pattern and determines cell fate. Making use of this experimental system and mathematical modeling we show that lateral inhibition patterning requires that ligand-receptor interactions are stronger within the same cell than in neighboring cells. Altogether, we establish the zebrafish notochord as an experimental system to study pattern generation, and identify and characterize how the properties of GRNs determine self-organization of gene patterning and cell fate. 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. wnt16 exerts pleiotropic effects on bone and lean mass in zebrafish Osteoporosis, a disease of bone fragility, and sarcopenia, a condition of reduced muscle mass and strength, commonly occur in the same individual--a condition termed osteosarcopenia. Genetic variants at the CPED1-WNT16 locus have pleiotropic effects on human bone mineral density (BMD) and lean tissue mass. However, our current understanding of Wnt16 and its influence on bone fails to explain effects on lean mass at this locus. In this study, we sought to better understand the biological mechanism underlying pleiotropy at the CPED1-WNT16 locus on BMD and lean mass. Using single cell analysis, microCT imaging, and genetic approaches in zebrafish, we show that wnt16 exerts dual influence on muscle and bone, and provide evidence that this occurs through a function in somite development, a period in which precursors for both of these tissues are specified. We show that mutations in wnt16 alter vertebral size and lean mass during axial skeletogenesis. We show that wnt16+ cells transcriptionally resemble cells within the amniote dermomyotome, and provide evidence they function as myogenic precursors during embryonic development. We also provide evidence that WNT16 is a gene of major effect on lean mass at the CPED-WNT16 locus. Our studies support the potential for wnt16 to exert pleiotropic effects on bone and lean mass through a function in myogenic precursors. 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. Fear memory-associated synaptic and mitochondrial changes revealed by deep learning-based processing of electron microscopy data Serial section electron microscopy (ssEM) can provide comprehensive 3D ultrastructural information of the brain with exceptional computational cost. Targeted reconstruction of subcellular structures from ssEM datasets is less computationally demanding but still highly informative. We thus developed a Region-CNN-based deep learning method to identify, segment, and reconstruct synapses and mitochondria to explore the structural plasticity of synapses and mitochondria in the auditory cortex of mice subjected to fear conditioning. Upon reconstructing over 135,000 mitochondria and 160,000 synapses, we found that fear conditioning significantly increases the number of mitochondria but decreases their size, and promotes the formation of multi-contact synapses comprising a single axonal bouton and multiple postsynaptic sites from different dendrites. Modeling indicates that such multi-contact configuration increases the information storage capacity of new synapses by over 50%. With high accuracy and speed in reconstruction, our method yields structural and functional insight into cellular plasticity associated with fear learning.