Protein sequence sampling and prediction from structural data The increasing integration between protein engineering and machine learning has led to many interesting results. A problem still to solve is to evaluate the likelihood that a sequence will fold into a target structure. This problem can be also viewed as sequence prediction from a known structure.

In the current work, we propose improvements in the recent architecture of Geometric Vector Perceptrons [1] in order to optimize the sampling of sequences from a known backbone structure. The proposed model differs from the original in that there is: (i) no updating in the vectorial embedding, only in the scalar one, (ii) only one layer of decoding. The first aspect improves the accuracy of the model and reduces the use of memory, the second allows for training of the model with several tasks without incurring data leakage.

We treat the trained classifier as an Energy-Based Model and sample sequences by sampling amino acids in a non-autoreggresive manner in the empty positions of the sequence using energy-guided criteria and followed by random mutation optimization. We improve the median identity of samples from 40.2% to 44.7%.

An additional question worth investigating is whether sampled and original sequences fold into similar structures independent of their identity. We chose proteins in our test set whose sampled sequences show low identity (under 30%) but for which our model predicted favorable energies. We used AlphaFold [2, 3] and observed that the predicted structures for sampled sequences highly resemble the predicted structures for original sequences, with an average TM-score of 0.84.
Pipeline to detect the positional and directional relationship between transposable elements and adjacent genes in host genome AO_SCPLOWBSTRACTC_SCPLOWUnderstanding the relationship between transposable elements (TEs) and their closest positional genes in the host genome is a key point to explore their potential role in genome evolution. Transposable elements can regulate and affect gene expression not only because of their mobility within the genome but also because of their transcriptional activity. A comprehensive knowledge of structural organization between transposable elements and neighboring genes is important to study TE functional role in gene regulation. We implemented a pipeline to display the relationship between the distribution of TEs and adjacent genes in host genomes. Our tool is freely available here: https://github.com/marieBvr/TEs_genes_relationship_pipeline.
First detection of herpesvirus and prevalence 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. Library-free BoxCarDIA solves the missing value problem in label-free quantitative proteomics The last decade has seen significant advances in the application of quantitative mass spectrometry-based proteomics technologies to tackle important questions in plant biology. The current standard for quantitative proteomics in plants is the use of data-dependent acquisition (DDA) analysis with or without the use of chemical labels. However, the DDA approach preferentially measures higher abundant proteins, and often requires data imputation due to quantification inconsistency between samples. In this study we systematically benchmarked a recently developed library-free data-independent acquisition (directDIA) method against a state-of-the-art DDA label-free quantitative proteomics workflow for plants. We next developed a novel acquisition approach combining MS1-level BoxCar acquisition with MS2-level directDIA analysis that we call BoxCarDIA. DirectDIA achieves a 33% increase in protein quantification over traditional DDA, and BoxCarDIA a further 8%, without any changes in instrumentation, offline fractionation, or increases in mass-spectrometer acquisition time. BoxCarDIA, especially, offers wholly reproducible quantification of proteins between replicate injections, thereby addressing the long-standing missing-value problem in label-free quantitative proteomics. Further, we find that the gains in dynamic range sampling by directDIA and BoxCarDIA translate to deeper quantification of key, low abundant, functional protein classes (e.g., protein kinases and transcription factors) that are underrepresented in data acquired using DDA. We applied these methods to perform a quantitative proteomic comparison of dark and light grown Arabidopsis cell cultures, providing a critical resource for future plant interactome studies. Our results establish BoxCarDIA as the new method of choice in quantitative proteomics using Orbitrap-type mass-spectrometers, particularly for proteomes with large dynamic range such as that of plants. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. 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. Loss of Acta2 in cardiac fibroblasts does not prevent myofibroblast differentiation or affect cardiac repair after myocardial infarction In response to myocardial infarction (MI), quiescent cardiac fibroblasts differentiate into myofibroblasts mediating tissue repair in the infarcted area. One of the most widely accepted markers of myofibroblast differentiation is the expression of Acta2 which encodes smooth muscle alpha-actin (SMA) that is assembled into stress fibers. However, the requirement of Acta2 in the myofibroblast differentiation of cardiac fibroblasts and its role in post-MI cardiac repair were still not known. To answer these questions, we generated a tamoxifen-inducible cardiac fibroblast-specific Acta2 knockout mouse line. Surprisingly, mice that lacked Acta2 in cardiac fibroblasts had a normal survival rate after MI. Moreover, Acta2 deletion did not affect the function or overall histology of infarcted hearts. No difference was detected in the proliferation, migration, or contractility between WT cardiac fibroblasts and Acta2-null cardiac myofibroblasts. Additional analysis identified that Acta2-null cardiac myofibroblasts had a normal total filamentous actin level and total actin level. Acta2 deletion caused a unique compensatory increase in the transcription level of Actg2 and a possible increase in the protein abundance of cytoplasmic actin isoforms. In conclusion, the deletion of Acta2 does not prevent the myofibroblast differentiation of cardiac fibroblasts or affect the post-MI cardiac repair, and the increase in the expression of non-SMA actin isoforms and the functional redundancy between actin isoforms are likely able to compensate for the loss of Acta2 in cardiac myofibroblasts. Fine-scale adaptations to environmental variation and growth strategies drive phyllosphere Methylobacterium diversity. Methylobacterium is a prevalent bacterial genus of the phyllosphere. Despite its ubiquity, little is known about the extent to which its diversity reflects neutral processes like migration and drift, versus environmental filtering of life history strategies and adaptations. In two temperate forests, we investigated how phylogenetic diversity within Methylobacterium was structured by biogeography, seasonality, and growth strategies. Using deep, culture-independent barcoded marker gene sequencing coupled with culture-based approaches, we uncovered a considerable diversity of Methylobacterium in the phyllosphere. We cultured different subsets of Methylobacterium lineages depending upon the temperature of isolation and growth (20 {degrees}C or 30 {degrees}C), suggesting long-term adaptation to temperature. To a lesser extent than temperature adaptation, Methylobacterium diversity was also structured across large (>100km; between forests) and small geographical scales (<1.2km within forests), among host tree species, and was dynamic over seasons. By measuring growth of 79 isolates at different temperature treatments, we observed contrasting growth performances, with strong lineage- and season-dependent variations in growth strategies. Finally, we documented a progressive replacement of lineages with a high-yield growth strategy typical of cooperative, structured communities, in favor of those characterized by rapid growth, resulting in convergence and homogenization of community structure at the end of the growing season. Together our results show how Methylobacterium is phylogenetically structured into lineages with distinct growth strategies, which helps explain their differential abundance across regions, host tree species, and time. This works paves the way for further investigation of adaptive strategies and traits within a ubiquitous phyllosphere genus.

ImportanceMethylobacterium is a bacterial group tied to plants. Despite its ubiquity and importance to their hosts, little is known about the processes driving Methylobacterium community dynamics. By combining traditional culture-dependent and independent (metabarcoding) approaches, we monitored Methylobacterium diversity in two temperate forests over a growing season. On the surface of tree leaves, we discovered remarkably diverse and dynamic Methylobacterium communities over short temporal (from June to October) and spatial scales (within 1.2 km). Because we cultured different subsets of Methylobacterium diversity depending on the temperature of incubation, we suspected that these dynamics partly reflected climatic adaptation. By culturing strains in lab conditions mimicking seasonal variations, we found that diversity and environmental variations were indeed good predictors of Methylobacterium growth performances. Our findings suggest that Methylobacterium community dynamics at the surface of tree leaves results from the succession of strains with contrasted growth strategies in response to environmental variations. Thermal modulation of Zebrafish exploratory statistics reveals constraints on individual behavioral variability BackgroundVariability is a hallmark of animal behavior. It contributes to survival by endowing individuals and populations with the capacity to adapt to ever-changing environmental conditions. Intra-individual variability is thought to reflect both endogenous and exogenous modulations of the neural dynamics of the central nervous system. However, how variability is internally regulated and modulated by external cues remains elusive. Here we address this question by analyzing the statistics of spontaneous exploration of freely swimming zebrafish larvae, and by probing how these locomotor patterns are impacted when changing the water temperatures within an ethologically relevant range.

ResultsWe show that, for this simple animal model, five short-term kinematic parameters - interbout interval, turn amplitude, travelled distance, turn probability and orientational flipping rate - together control the long-term exploratory dynamics. We establish that the bath temperature consistently impacts the means of these parameters, but leave their pairwise covariance unchanged. These results indicate that the temperature merely controls the sampling statistics within a well-defined kinematic space delineated by this robust statistical structure. At a given temperature, individual animals explore the behavioral space over a timescale of tens of minutes, suggestive of a slow internal state modulation that could be externally biased through the bath temperature. By combining these various observations into a minimal stochastic model of navigation, we show that this thermal modulation of locomotor kinematics results in a thermophobic behavior, complementing direct gradient-sensing mechanisms.

ConclusionsThis study establishes the existence of a well-defined locomotor space accessible to zebrafish larvae during spontaneous exploration, and quantifies self-generated modulation of locomotor patterns. Intra-individual variability reflects a slow diffusive-like probing of this space by the animal. The bath temperature in turn restricts the sampling statistics to sub-regions, endowing the animal with basic thermophobicity. This study suggests that in Zebrafish, as well as in other ectothermic animals, ambient temperature could be used to efficiently manipulate internal states in a simple and ethological way. Topology and adenocarcinoma cell localization dataset on the labyrinthin biomarker. The discovery and characterization of tumor associated antigens is increasingly important to advance the field of immuno-oncology. Data on the, 1) topology, 2) amino acid (a.a.) homology analyses and 3) cell surface localization are presented as evidence in support of labyrinthin as a pan-adenocarcinoma marker. Multiple bioinformatics tools were used to collectively predict calcium binding domain(s), N-myristoylation sites, kinase II phosphorylation sites, and transmembrane domain(s) with predicted N- and C-termini orientations. Sequence homology comparisons were found and alignments made for labyrinthin (255 a.a.) vs. the intracellular aspartyl/asparaginyl beta-hydroxylase (ASPH; 758 a.a.) and the ASPH-gene related protein junctate (299 a.a.), which are both type II proteins. The epitopes for commonly used MCA 44-3A6 anti-labyrinthin and FB-50 anti-ASPH/junctate antibodies were present in all three proteins. The selective expression of labyrinthin was investigated by third party fluorescent activated cell sorter (FACS) analyses of MCA 44-3A6 binding to non-permeabilized adenocarcinoma and normal cell line models. Labyrinthin was detected on non-permeabilized A549 human lung adenocarcinoma cells, but not on normal WI-38 human lung fibroblasts nor primary cultures of normal human glandular-related cells. Microscopic images of immunofluorescent labelled MCA 44-3A6 binding to A549 cells at random cell cycle stages complement the FACS results by further showing that labyrinthin persisted on the cell surfaces along with some cell internalization within 20 minutes. The dataset presented herein is related to research article: "Labyrinthin: A distinct pan-adenocarcinoma diagnostic and immunotherapeutic tumor specific antigen" by Babich et al., Heliyon, 2022 [1].

Specifications Table

O_TBL View this table:
org.highwire.dtl.DTLVardef@19222bborg.highwire.dtl.DTLVardef@1264bcdorg.highwire.dtl.DTLVardef@13a8158org.highwire.dtl.DTLVardef@1dcf0org.highwire.dtl.DTLVardef@503a9f_HPS_FORMAT_FIGEXP  M_TBL C_TBL Value of the data O_LIBecause identification and characterization of neoantigens [3, 4, 5] are rare yet increasingly important for diagnostics and therapeutics, the identification of labyrinthin as a pan-adenocarcinoma target is significant for further study by clinicians and cancer researchers.
C_LIO_LIThe data provide a basis to determine to what extent labyrinthin is a distinguishing marker for adenocarcinomas vs. conventional methods.
C_LIO_LIThe results indicate that therapeutic strategies can be employed that take advantage of labyrinthins convenient presence on adenocarcinoma cell surfaces.
C_LIO_LIThe discovery of a calcium binding sequence in labyrinthin provides a basis to research if it is key to explaining the long-known phenomena of elevated calcium and irregular signaling in cancer cells [8, 9].
C_LIO_LIThe data reveal a need to explore distinctions between ASPH and labyrinthin, in particular because ASPH [10, 11] and labyrinthin-based treatments have already been FDA approved for clinical trials.
C_LI Common variants in breast cancer risk loci predispose to distinct tumor subtypes BackgroundGenome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER), but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

MethodsAmong 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

ResultsEighty-five of 173 variants were associated with at least one tumor feature (false discovery rate <5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at P<0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

ConclusionThis report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. Modular Dynamic Biomolecular Modelling with Bond Graphs: The Unification of Stoichiometry, Thermodynamics, Kinetics and Data. Renewed interest in dynamic simulation models of biomolecular systems has arisen from advances in genome-wide measurement and applications of such models in biotechnology and synthetic biology. In particular, genome-scale models of cellular metabolism beyond the steady state are required in order to represent transient and dynamic regulatory properties of the system. Development of such whole-cell models requires new modelling approaches. Here we propose the energy-based bond graph methodology, which integrates stoichiometric models with thermo-dynamic principles and kinetic modelling. We demonstrate how the bond graph approach intrinsically enforces thermodynamic constraints, provides a modular approach to modelling, and gives a basis for estimation of model parameters leading to dynamic models of biomolecular systems. The approach is illustrated using a well-established stoichiometric model of Escherichia coli (E. coli) and published experimental data. Diffusion and distal linkages govern interchromosomal dynamics during meiotic prophase The pairing of homologous chromosomes (homologs) in meiosis is essential for distributing the correct numbers of chromosomes into haploid gametes. In budding yeast, pairing depends on the formation of 150-200 Spo11-mediated double-strand breaks (DSBs) that are distributed among 16 homolog pairs, but it is not known if all, or only a subset of these DSBs, contribute to the close juxtaposition of homologs. Having established a system to measure the position of fluorescently tagged chromosomal loci in 3D space over time, we analyzed locus trajectories to determine how frequently, and how long, loci spend colocalized or apart. Continuous imaging revealed highly heterogeneous cell-to-cell behavior of foci, with the majority of cells exhibiting a "mixed" phenotype where foci move into and out of proximity, even at late stages of prophase, suggesting that the axial structures of the synaptonemal complex may be more dynamic than anticipated. The observed plateaus of the mean-squared change in distance (MSCD) between foci informed the development of a biophysical model of two diffusing polymers that captures the loss of centromere linkages as cells enter meiosis, nuclear confinement, and the formation of Spo11-dependent linkages. The predicted number of linkages per chromosome in our theoretical model closely approximates the small number (~2-4) of estimated synapsis-initiation sites, suggesting that excess DSBs have negligible effects on the overall juxtaposition of homologs. These insights into the dynamic in-terchromosomal behavior displayed during homolog pairing demonstrate the power of combining time-resolved in vivo analysis with modeling at the granular level.

Significance StatementEssential for sexual reproduction, meiosis is a specialized cell division required for the production of haploid gametes. Critical to this process is the pairing, recombination, and segregation of homologous chromosomes (homologs). While pairing and recombination are linked, it is not known how many linkages are sufficient to hold homologs in proximity. Here, we reveal that random diffusion and the placement of a small number of linkages are sufficient to establish the apparent "pairing" of homologs. We also show that colocalization between any two loci is more dynamic than anticipated. Our study is the first to provide observations of live interchromosomal dynamics during meiosis and illustrates the power of combining single-cell measurements with theoretical polymer modeling. A framework for assessing the skill and value of operational recruitment forecasts Forecasting variation in the recruitment to fish stocks is one of the most challenging and long-running problems in fisheries science and essentially remains unsolved today. Traditionally recruitment forecasts are developed and evaluated based on explanatory and goodness-of-fit approaches that do not reflect their ability to predict beyond the data on which they were developed. Here we propose a new generic framework that allows the skill and value of recruitment forecasts to be assessed in a manner that is relevant to their potential use in an operational setting. We assess forecast skill based on predictive power using a retrospective forecasting approach inspired by meterology, and emphasise the importance of assessing these forecasts relative to a baseline. We quantify the value of these forecasts using an economic cost-loss decision model that is directly relevant to many forecast users. We demonstrate this framework using four stocks of lesser sandeel (Ammodytes marinus) in the North Sea, showing for the first time in an operationally realistic setting that skilful and valuable forecasts are feasible in two of these areas. This result shows the ability to produce valuable short-term recruitment forecasts, and highlights the need to revisit our approach to and understanding of recruitment forecasting. A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. 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. Omicron booster in ancestral strain vaccinated mice augments protective immunities against both the Delta and Omicron variants A booster vaccination is called for constraining the evolving epidemic of SARS-CoV-2. However, the necessity of a new COVID-19 vaccine is currently unclear. To compare the effect of an Omicron-matched S DNA vaccine and an ancestral S DNA vaccine in boosting cross-reactive immunities, we firstly immunized mice with two-dose of a DNA vaccine encoding the spike protein of the ancestral Wuhan strain. Then the mice were boosted with DNA vaccines encoding spike proteins of either the Wuhan strain or the Omicron variant. Specific antibody and T cell responses were measured at 4 weeks post boost. Our data showed that the Omicron-matched vaccine efficiently boosted RBD binding antibody and neutralizing antibody responses against both the Delta and the Omicron variants. Of note, antibody responses against the Omicron variant elicited by the Omicron-matched vaccine were much stronger than those induced by the ancestral S DNA vaccine. Meanwhile, CD8+ T cell responses against both the ancestral Wuhan strain and the Omicron strain also tended to be higher in mice boosted by the Omicron-matched vaccine than those in mice boosted with the ancestral S DNA vaccine, albeit no significant difference was observed. Our findings suggest that an Omicron-matched vaccine is preferred for boosting cross-reactive immunities. The Impact of Neuron Morphology on Cortical Network Architecture It has become increasingly clear that the neurons in the cerebral cortex are not randomly interconnected. This wiring specificity can result from synapse formation mechanisms that interconnect neurons depending on their activity or genetically defined identity. Here we report that in addition to these synapse formation mechanisms, the structural composition of the neuropil provides a third prominent source by which wiring specificity emerges in cortical networks. This structurally determined wiring specificity reflects the packing density, morphological diversity and similarity of the dendritic and axonal processes. The higher these three factors are, the more recurrent the topology of the networks. Conversely, low density, diversity and similarity yield feedforward networks. These principles predict connectivity patterns from subcellular to network scales that are remarkably consistent with empirical observations from a rich body of literature. Thus, cortical network architectures reflect the specific morphological properties of their constituents to much larger degrees than previously thought. Behavioral and postural analyses establish sleep-like states for mosquitoes that can impact host landing and blood feeding Sleep is an evolutionarily conserved process that has been described in different animal systems. For insects, sleep characterization has been primarily achieved using behavioral and electrophysiological correlates in a few systems. Sleep in mosquitoes, which are important vectors of disease-causing pathogens, has not been directly examined. This is surprising as circadian rhythms, which have been well studied in mosquitoes, influence sleep in other systems. In this study, we characterized sleep in mosquitoes using body posture analysis and behavioral correlates, and quantified the effect of sleep deprivation on sleep rebound, host landing and blood-feeding propensity. Body and appendage position metrics revealed a clear distinction between the posture of mosquitoes in their putative sleep and awake states for multiple species, which correlate with a reduction in responsiveness to host cues. Sleep assessment informed by these posture analyses indicated significantly more sleep during periods of low activity. Nighttime and daytime sleep deprivation resulting from the delivery of vibration stimuli induced sleep rebound in the subsequent phase in day and night active mosquitoes, respectively. Lastly, sleep deprivation suppressed host landing in both laboratory and field settings and also impaired blood feeding of a human host when mosquitoes would normally be active. These results suggest that quantifiable sleep states occur in mosquitoes, and highlight the potential epidemiological importance of mosquito sleep. Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent global expansion A novel coronavirus, SARS-CoV-2, has caused over 190 million cases and over 4 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e. M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic one year later. By analyzing 4,013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of initial phase of global ignition over a timespan of 14-week. 11 major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted a half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implying that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261,323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e. including all the publicly available viral genomes) in order to recapitulate our findings over one-year timespan. By 25 December 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the four variants of concern (VOC) are the descendants of WE1 type. This study demonstrates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time. Highly efficient synthetic CRISPR RNA/Cas9-based mutagenesis for rapid cardiovascular phenotypic screening in F0 zebrafish The zebrafish is a valuable vertebrate model to study cardiovascular formation and function due to the facile visualization and rapid development of the circulatory system in its externally growing embryos. Despite having distinct advantages, zebrafish have paralogs of many important genes, making reverse genetics approaches inefficient since generating animals bearing multiple gene mutations requires substantial efforts. Here, we present a simple and robust synthetic CRISPR RNA/Cas9-based mutagenesis approach for generating biallelic F0 zebrafish knockouts. Using a dual-guide synthetic CRISPR RNA/Cas9 ribonucleoprotein (dgRNP) system, we compared the efficiency of biallelic gene disruptions following the injections of one, two, and three dgRNPs per gene into the cytoplasm or yolk. We show that simultaneous cytoplasmic injections of three distinct dgRNPs per gene into one-cell stage embryos resulted in the most efficient and consistent biallelic gene disruptions. Importantly, this triple dgRNP approach enables efficient inactivation of cell autonomous and cell non-autonomous gene function, likely due to the low mosaicism of biallelic disruptions. In support of this finding, we provide evidence that the F0 animals generated by this method fully phenocopied the endothelial and peri-vascular defects observed in corresponding stable mutant homozygotes. Moreover, this approach faithfully recapitulated the trunk vessel phenotypes resulting from the genetic interaction between two vegfr2 zebrafish paralogs. Mechanistically, investigation of genome editing and mRNA decay indicates that the combined mutagenic actions of three dgRNPs per gene lead to an increased probability of frameshift mutations, enabling efficient biallelic gene disruptions. Therefore, our approach offers a highly robust genetic platform to quickly assess novel and redundant gene function in F0 zebrafish. A novel Fiji/ImageJ plugin for the rapid analysis of blebbing cells When confined, cells have recently been shown to undergo a phenotypic switch to what has been termed, fast amoeboid (leader bleb-based) migration. However, as this is a nascent area of research, few tools are available for the rapid analysis of cell behavior. Here, we demonstrate that a novel Fiji/ImageJ-based plugin, Analyze_Blebs, can be used to quickly obtain cell migration parameters and morphometrics from time lapse images. As validation, we show that Analyze_Blebs can detect significant differences in cell migration and morphometrics, such as the largest bleb size, upon introducing different live markers of F-actin, including F-tractin and LifeAct tagged with green and red fluorescent proteins. We also demonstrate, using flow cytometry, that live markers increase total levels of F-actin. Furthermore, that F-tractin increases cell stiffness, which was found to correlate with a decrease in migration, thus reaffirming the importance of cell mechanics as a determinant of Leader Bleb-Based Migration (LBBM).

Highlight summaryA new plugin, Analyze_Blebs, enables the rapid analysis of cell migration and morphometrics of fast amoeboid cells. Morphometrics combined with cell stiffness are found to be predictive of confined, fast amoeboid (leader bleb-based), migration. Gene loss during a transition to multicellularity Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being highly similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited a high degree of novel interactions, suggesting that gene loss plays a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest developmental complexity may be driven by gene loss, and not just by gene gain.

Significance StatementIncreased developmental complexity is thought to evolve mainly through genetic innovation and co-option. Comparison of the genomes of the closely related volvocines revealed that despite large changes in developmental complexity their genomes are exhibiting significant gene loss. Further, a burst of gene loss is shown to occur at the transition to undifferentiated multicellularity by gene inactivation and decay. This likely results in changes to protein-protein interactions within the cell, suggesting evolutionary novelty does not always require gene gain. Using empirical and modelling approaches we demonstrate gene loss can more easily produce protein network novelty than gene gain. We propose that gene loss can be a driver of biological innovation, allowing for reconfiguration of gene networks and differential use of existing functional repertoires. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics.