Presence of a SARS-COV-2 protein enhances Amyloid Formation of Serum Amyloid A A marker for the severeness and disease progress of COVID-19 is overexpression of serum amyloid A (SAA) to levels that in other diseases are associated with a risk for SAA amyloidosis. In order to understand whether SAA amyloidosis could also be a long-term risk of SARS-COV-2 infections we have used long all-atom molecular dynamic simulations to study the effect of a SARS-COV-2 protein segment on SAA amyloid formation. Sampling over 40 {micro}s we find that presence of the nine-residue segment SK9, located at the C-terminus of the Envelope protein, increases the propensity for SAA fibril formation by three mechanisms: it reduces the stability of the lipid-transporting hexamer shifting the equilibrium toward monomers, it increases the frequency of aggregation-prone configurations in the resulting chains, and it raises the stability of SAA fibrils. Our results therefore suggest that SAA amyloidosis and related pathologies may be a long-term risk of SARS-COV-2 infections.
A quartz crystal microbalance method to quantify the size of hyaluronan and other glycosaminoglycans on surfaces Hyaluronan (HA) is a major component of peri- and extra-cellular matrices and plays important roles in many biological processes such as cell adhesion, proliferation and migration. The abundance, size distribution and presentation of HA dictate its biological effects and are also useful indicators of pathologies and disease progression. Methods to assess the molecular mass of free-floating HA and other glycosaminoglycans (GAGs) are well established. In many biological and technological settings, however, GAGs are displayed on surfaces, and methods to obtain the size of surface-attached GAGs are lacking. Here, we present a method to size HA that is end-attached to surfaces. The method is based on the quartz crystal microbalance with dissipation monitoring (QCM-D) and exploits that the softness and thickness of films of grafted HA increase with HA size. These two quantities are sensitively reflected by the ratio of the dissipation shift ({Delta}D) and the negative frequency shift (-{Delta}f) measured by QCM-D upon the formation of HA films. Using a series of size-defined HA preparations, ranging in size from [~]2 kDa tetrasaccharides to [~]1 MDa polysaccharides, we establish a monotonic yet non-linear standard curve of the {Delta}D/-{Delta}f ratio as a function of HA size, which reflects the distinct conformations adopted by grafted HA chains depending on their size and surface coverage. We demonstrate that the standard curve can be used to determine the mean size of HA, as well as other GAGs, such as chondroitin sulfate and heparan sulfate, of preparations of previously unknown size in the range from 1 to 500 kDa, with a resolution of better than 10%. For polydisperse samples, our analysis shows that the process of surface-grafting preferentially selects smaller GAG chains, and thus reduces the average size of GAGs that are immobilised on surfaces comparative to the original solution sample. Our results establish a quantitative method to size HA and other GAGs grafted on surfaces, and also highlight the importance of sizing GAGs directly on surfaces. The method should be useful for the development and quality control of GAG-based surface coatings in a wide range of research areas, from molecular interaction analysis to biomaterials coatings.
Discovery of highly potent pancoronavirus fusion inhibitors that also effectively inhibit COVID-19 variants from the UK (Alpha), South Africa (Beta), and India (Delta) We report here the discovery of several highly potent small molecules that showed low nM potency against SARS-CoV (IC50: as low as 13 nM), SARS-CoV-2 (IC50: as low as 23 nM), and MERS-CoV (IC50: as low as 76 nM) in pseudovirus based assays with excellent selectivity indices (SI: as high as > 5000) demonstrating their pancoronavirus inhibition. Some compounds also show 100% inhibition of CPE (IC100) at 1.25 {micro}M against an authentic SARS-CoV-2 (US_WA-1/2020). Furthermore, the most active inhibitors also potently inhibited variants of concerns (VOCs), such as the UK (B.1.1.7), South Africa (B.1.351), and Delta variant (B.1.617.2), originated in India. We confirmed that one of the potent inhibitors binds to the prefusion spike protein trimer of SARS-CoV-2 by SPR. Besides, we showed that they inhibit virus-mediated cell-cell fusion. The ADME data of one of the most active inhibitors, NBCoV1, show drug-like properties. In vivo PK of NBCoV1 in rats demonstrated excellent half-life (t1/2) of 11.3 h, mean resident time (MRT) of 14.2 h, and oral bioavailability. We expect the lead inhibitors to pave the way for further development to preclinical and clinical candidates. C. elegans PEZO-1 is a Mechanosensitive Ion Channel Involved in Food Sensation PIEZO channels are force sensors essential for physiological processes, including baroreception and proprioception. The Caenorhabditis elegans genome encodes an ortholog gene of the Piezo family, pezo-1, which is expressed in several tissues, including the pharynx. This myogenic pump is an essential component of the C. elegans alimentary canal, whose contraction and relaxation are modulated by mechanical stimulation elicited by food content. Whether pezo-1 encodes a mechanosensitive ion channel and contributes to pharyngeal function remain unknown. Here, we leverage genome editing, genetics, microfluidics, and electropharyngeogram recording to establish that pezo-1 is expressed in the pharynx, including in a proprioceptive-like neuron, and regulates pharyngeal function. Knockout (KO) and gain-of-function (GOF) mutants reveal that pezo-1 is involved in fine-tuning pharyngeal pumping frequency, sensing osmolarity, and food mechanical properties. Using pressure-clamp experiments in primary C. elegans embryo cultures, we determine that pezo-1 KO cells do not display mechanosensitive currents, whereas cells expressing wild-type or GOF PEZO-1 exhibit mechanosensitivity. Moreover, infecting the Spodoptera frugiperda cell line with a baculovirus containing the G-isoform of pezo-1 (among the longest isoforms) demonstrates that pezo-1 encodes a mechanosensitive channel. Our findings reveal that pezo-1 is a mechanosensitive ion channel that regulates food sensation in worms. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. An alternative UPF1 isoform drives conditional remodeling of nonsense-mediated mRNA decay The nonsense-mediated mRNA decay (NMD) pathway monitors translation termination to degrade transcripts with premature stop codons and regulate thousands of human genes. Here we show that an alternative mammalian-specific isoform of the core NMD factor UPF1, termed UPF1LL, enables condition-dependent remodeling of NMD specificity. Previous studies indicate that the extension of a conserved regulatory loop in the UPF1LL helicase core confers a decreased propensity to dissociate from RNA upon ATP hydrolysis relative to the major UPF1 isoform, designated UPF1 SL. Using biochemical and transcriptome-wide approaches, we find that UPF1LL overcomes the protective RNA binding proteins PTBP1 and hnRNP L to preferentially bind and down-regulate long 3UTRs normally shielded from NMD. Unexpectedly, UPF1LL supports induction of NMD on new populations of substrate mRNAs in response to activation of the integrated stress response and impaired translation efficiency. Thus, while canonical NMD is abolished by moderate translational repression, UPF1LL activity is enhanced, providing a mechanism to rapidly rewire NMD specificity in response to cellular stress. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. Identifying and testing marker-trait associations for growth and phenology in three pine species: implications for genomic prediction In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 SNPs from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments. The coiled-coil domain of E. coli FtsLB is a structurally detuned element critical for modulating its activation in bacterial cell division The FtsLB complex is a key regulator of bacterial cell division, existing in either an off or on state which supports the activation of septal peptidoglycan synthesis. In Escherichia coli, residues known to be critical for this activation are located in a region near the C-terminal end of the periplasmic coiled-coil domain of FtsLB, raising questions about the precise role of this conserved domain in the activation mechanism. Here, we investigate an unusual cluster of polar amino acids found within the core of the FtsLB coiled coil. We hypothesized that these amino acids likely reduce the structural stability of the domain and thus may be important for governing conformational changes. We found that mutating these positions to hydrophobic residues increased the thermal stability of FtsLB but caused cell division defects, suggesting that the coiled-coil domain is a "detuned" structural element. In addition, we identified suppressor mutations within the polar cluster, indicating that the precise identity of the polar amino acids is important for fine-tuning the structural balance between the off and on states. We propose a revised structural model of the tetrameric FtsLB (named the "Y-model") in which the periplasmic domain splits into a pair of coiled-coil branches. In this configuration, the hydrophilic terminal moieties of the polar amino acids remain more favorably exposed to water than in the original four-helix bundle model ("I-model"). We propose that a shift in this architecture, dependent on its marginal stability, is involved in activating the FtsLB complex and triggering septal cell wall reconstruction. VeloViz: RNA-velocity informed embeddings for visualizing cellular trajectories 0Single cell transcriptomic technologies enable genome-wide gene expression measurements in individual cells but can only provide a static snapshot of cell states. RNA velocity analysis can infer cell state changes from single cell transcriptomics data. To interpret these cell state changes as part of underlying cellular trajectories, current approaches rely on visualization with principal components, t-distributed stochastic neighbor embedding, and other 2D embeddings derived from the observed single cell transcriptional states. However, these 2D embeddings can yield different representations of the underlying cellular trajectories, hindering the interpretation of cell state changes. We developed VeloViz to create RNA-velocity-informed 2D and 3D embeddings from single cell transcriptomics data. Using both real and simulated data, we demonstrate that VeloViz embeddings are able to consistently capture underlying cellular trajectories across diverse trajectory topologies, even when intermediate cell states may be missing. By taking into consideration the predicted future transcriptional states from RNA velocity analysis, VeloViz can help visualize a more reliable representation of underlying cellular trajectories. VeloViz is available as an R package on GitHub (https://github.com/JEFworks-Lab/veloviz) with additional tutorials at https://JEF.works/veloviz/. LabGym: quantification of user-defined animal behaviors using learning-based holistic assessment Quantifying animal behavior is important for many branches of biological research. Current computational tools for quantifying behaviors typically rely on a few pre-defined, simplified features to identify a behavior. However, such an approach underuses the behavioral information and restricts the tools applicability to limited behavior types or animal species. Here we report a new tool, LabGym, for quantifying animal behaviors without such limitations. Combining a novel approach for effective evaluation of animal motion with customizable convolutional recurrent networks for capturing spatiotemporal details, LabGym provides holistic behavioral assessment and accurately identify user-defined animal behaviors without restrictions on behavior types or animal species. It then calculates various parameters to provide quantitative measurements of the intensity and dynamics of each behavior. LabGym requires neither any intermediate step that causes information loss nor programming knowledge from users for post-hoc analysis. It tracks multiple animals simultaneously in various experimental settings for high-throughput and versatile analysis, and provides users a way to generate visualizable behavioral datasets that are valuable resources for research in neuroscience, behavioral science, and computer science. We demonstrate its efficacy in capturing subtle behavioral changes in animals ranging from soft-bodied invertebrates to mammals. Poly(rC)-binding protein 1 limits hepatitis C virus assembly and secretion The hepatitis C virus (HCV) co-opts numerous cellular elements - including proteins, lipids, and microRNAs - to complete its viral life cycle. The cellular RNA-binding protein, poly(rC)-binding protein 1 (PCBP1), was previously reported to bind to the 5 untranslated region (UTR) of the HCV genome; however, its importance in the viral life cycle has remained unclear. Herein, we sought to clarify the role of PCBP1 in the HCV life cycle. Using the HCV cell culture (HCVcc) system, we found that knockdown of endogenous PCBP1 resulted in an overall decrease in viral RNA accumulation, yet resulted in an increase in extracellular viral titers. To dissect PCBP1s specific role in the HCV life cycle, we carried out assays for viral entry, translation, genome stability, RNA replication, as well as virion assembly and secretion. We found that PCBP1 knockdown did not directly affect viral entry, translation, RNA stability, or RNA replication, but resulted in an overall increase in infectious particle secretion. This increase in virion secretion was evident even when viral RNA synthesis was inhibited, and blocking virus secretion could partially restore the viral RNA accumulation decreased by PCBP1 knockdown. We therefore propose a model where endogenous PCBP1 normally limits virion assembly and secretion, which increases viral RNA accumulation in infected cells by preventing the departure of viral genomes packaged into virions. Overall, our findings improve our understanding of how cellular RNA-binding proteins influence viral genomic RNA utilization during the HCV life cycle. Gradients in cellular metabolism regulate development of tonotopy along the cochlea. In vertebrates with elongated auditory organs, mechanosensory hair cells (HCs) are organised such that complex sounds are broken down into their component frequencies along the basal-to-apical long (tonotopic) axis. To generate the frequency-specific characteristics required at the appropriate positions, nascent HCs must determine their tonotopic properties during development. This relies on complex signalling between the developing HC and its local niche along the axis within the auditory organ. Here we apply NAD(P)H fluorescence lifetime imaging (FLIM) and live imaging of mitochondria to reveal metabolic gradients along the tonotopic axis of the developing cochlea. We further show that re-shaping these gradients during development, by inhibiting cytosolic glucose metabolism, alters normal Bmp7 and Chordin like-1 signalling leading to flattening of tonotopic patterning along the axis. Our work supports a causal link between morphogen signalling and metabolic reprogramming in specifying tonotopic morphologies in auditory HCs. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Assessment of Sex Bias in Housekeeping Gene Expression in Adipose Tissue Through the Massive Analysis of Transcriptomics Data As the housekeeping genes (HKGs) generally involved in maintaining essential cell functions are typically assumed to exhibit constant expression levels across cell types, they are commonly employed as internal controls in gene expression studies. Nevertheless, multiple studies indicate that not all HKGs display stable expression across cells and tissues and under various healthy and diseased conditions, which can introduce systematic errors into experimental results. The selection and validation of HKGs as controls for each studied condition represent crucial steps in ensuring the validity of obtained results; however, up till now, sex has not been typically considered as a biological variable.

In this study, we evaluate the expression profiles of six classical HKGs (four metabolic: GAPDH, HPRT, PPIA, and UBC, and two ribosomal: 18S and RPL19) to determine expression stability in adipose tissues of Homo sapiens and Mus musculus and check sex bias and their overall suitability as internal controls. We also assess the stabile expression of genes included in distinct whole-transcriptome microarrays available from the Gene Expression Omnibus database (GEO) to identify sex-unbiased HKGs suitable for use as internal controls. We perform a sex-based analysis to describe any sexual dimorphisms in mRNA expression stability.

A novel computational strategy based on meta-analysis techniques proves that certain classical HKGs fail to function adequately as controls when analyzing human adipose tissue (HAT) considering sex as a variable. The extensively used 18S gene displays sex-based variability in adipose tissue, although PPIA and RPL19 do not, and hence, represent robust HKGs. We propose new sex-unbiased human and mouse HKGs (suHKG) derived from sex-specific expression profiles, such as RPS8 and UBB. All results generated during this study are readily available by accessing an open web resource (https://bioinfo.cipf.es/metafun-HKG) for consultation and reuse in further studies. Embedded 3D printing in self-healing annealable composites for precise patterning of functionally mature human neural constructs Human in vitro models of neural tissue with controllable cellular identity, tunable microenvironment, and defined spatial arrangement are needed to facilitate studies of brain development and disease. Towards this end, embedded printing in jammed microgel supports (i.e., granular gels) holds great promise as it allows precise and programmable patterning of extremely soft and compliant tissue constructs. However, in contrast to the vast material landscape available for bulk hydrogels, granular printing support formulations are restricted to a handful of materials without the ability for facile adjustment of biofunctional properties of the cellular microenvironment. Therefore, there has been a need for novel materials that take advantage of versatile biomimicry of bulk hydrogels while providing high-fidelity support for embedded printing akin to granular gels. To address this need, we present a modular platform for bioengineering of neuronal networks via direct embedded 3D printing of human stem cells inside Self-Healing Annealable Particle-Extracellular matrix (SHAPE) composites. SHAPE composites consist of soft microgels immersed in viscous extracellular-matrix solution to enable precise freeform patterning of human stem cells and consequent generation and long-term maintenance of mature subtype-specific neurons that extend projections within the volume of the annealed support. The developed approach further allows multi-ink deposition, live spatial and temporal monitoring of oxygen levels, as well as creation of vascular channels. Due to its modularity, SHAPE biomanufacturing toolbox not only offers a solution for functional modeling of mechanically sensitive neural constructs, but also has potential to be applied to a wide range of biomaterials with different crosslinking mechanisms to model tissues and diseases where recapitulation of complex architectural features and topological cues is essential. A non-canonical convergence of carbohydrate and glutamine metabolism is required after metabolic rewiring in a solid environment The fate of pyruvate, which is modulated by the activity of the mitochondrial pyruvate carrier (MPC), is a defining metabolic feature in many cancers. Diffuse large B-cell lymphomas (DLBCLs) are a genetically and metabolically heterogeneous cancer. Although MPC expression and activity differed between DLBCL subgroups, mitochondrial pyruvate was uniformly consumed by glutamate pyruvate transaminase 2 (GPT2) to support -ketoglutarate production as part of glutaminolysis. This led us to discover that glutamine exceeds pyruvate as a carbon source for the tricarboxylic acid (TCA) cycle in DLBCLs. Furthermore, we found that MPC inhibition unexpectedly leads to decreased glutaminolysis, which is contrary to previous observations in other cell types. We also discovered that MPC inhibition and depletion only decreased DLBCL proliferation in an extracellular matrix (ECM) environment and in vivo xenografts, but not in the typical DLBCL suspension environment. We also have found that the metabolic profile of DLBCL cells in ECM is markedly different from cells in suspension environment. Thus, we report that besides the canonical glutamate dehydrogenase (GDH)-mediated glutaminolysis, the non-canonical GPT2 mediated consumption and assimilation of glutamine and pyruvate in DLBCLs enables their proliferation in an extracellular environment-dependent manner.

HIGHLIGHTSO_LIGlutamine, but not glucose, is a major carbon source for the tricarboxylic acid cycle in diffuse large B-cell lymphomas.
C_LIO_LIMitochondrial pyruvate supports glutaminolysis in diffuse large B-cell lymphomas by supplying pyruvate for glutamate pyruvate transaminase 2-mediated  -ketoglutarate production.
C_LIO_LIMitochondrial pyruvate carrier inhibition leads to decreased glutaminolysis in diffuse large B-cell lymphomas.
C_LIO_LI -ketoglutarate production is important for diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment.
C_LIO_LIMitochondrial pyruvate carrier activity supports diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment and in mouse xenografts.
C_LI Mathematical Modeling and Analysis of Mitochondrial Retrograde Signaling Dynamics Mitochondria, semi-autonomous eukaryotic organelles, participate in energy production and metabolism, making mitochondrial quality control crucial. As most mitochondrial proteins are encoded by nuclear genes, quality control depends on proper mitochondria-nucleus communication, designated mitochondrial retrograde signaling. Early studies focused on retrograde signaling participants and specific gene knockouts. However, mitochondrial signal modulation remains elusive. Using yeast, we simulated signal propagation following mitochondrial damage and proposed a mathematical model based ordinary differential equations. Mitochondrial retrograde signaling decisions were described using a Boolean model. Dynamics were analyzed through an ordinary differential equation-based model and extended to evaluate the model response to noisy damage signals. Simulation revealed localized protein concentration dynamics, including waveforms, frequency response, and robustness under noise. Retrograde signaling is bistable with three localized steady states, and increased damage compromises robustness. We elucidated mitochondrial retrograde signaling, thus providing a basis for drug design against yeast and fungi.

Author SummaryThe yeast RTG pathway regulates mitochondrial metabolism and mitochondrial quality through passing mitochondrial signal to the nucleus to modulate gene expressions. Using microscopic data of RTG proteins translocation from cytosol to nucleus, the parameters were found by fitting 16 knockout conditions simulated from the proposed differential equation-based model. Further, through dose response, ultrasensitivity, frequency response and noise interference, we demonstrated a switch-like property of RTG pathway activation, a capacity charging property of RTG protein translocation, a low pass filter property of signals in different frequencies, and noise amplification due to mitochondrial signal. The proposed mathematical model enables us to understand the dynamics and mechanisms in mitochondrial retrograde signaling in yeast and provides potential antifungal treatment strategies. Leveraging a surrogate outcome to improve inference on a partially missing target outcome Sample sizes vary substantially across tissues in the Genotype-Tissue Expression (GTEx) project, where considerably fewer samples are available from certain inaccessible tissues, such as the substantia nigra (SSN), than from accessible tissues, such as blood. This severely limits power for identifying tissue-specific expression quantitative trait loci (eQTL) in undersampled tissues. Here we propose Surrogate Phenotype Regression Analysis (SO_SCPLOWPRAYC_SCPLOW) for leveraging information from a correlated surrogate outcome (e.g. expression in blood) to improve inference on a partially missing target outcome (e.g. expression in SSN). Rather than regarding the surrogate outcome as a proxy for the target outcome, SO_SCPLOWPRAYC_SCPLOW jointly models the target and surrogate outcomes within a bivariate regression framework. Unobserved values of either outcome are treated as missing data. We describe and implement an expectation conditional maximization algorithm for performing estimation in the presence of bilateral outcome missingness. SO_SCPLOWPRAYC_SCPLOW estimates the same association parameter estimated by standard eQTL mapping and controls the type I error even when the target and surrogate outcomes are truly uncorrelated. We demonstrate analytically and empirically, using simulations and GTEx data, that in comparison with marginally modeling the target outcome, jointly modeling the target and surrogate outcomes increases estimation precision and improves power. Ineffectual AEC1 differentiation from KRT8hi transitional state without fibrosis is associated with fatal COVID-19 ARDS ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. How much does it take to be old? Modelling the time since the last harvesting to infer the distribution of overmature forests in France AimThe distribution of overmature forests in metropolitan France is poorly known, with only a few well-studied prominent sites, and has never been evaluated countrywide. Here, we modelled French forest reserves time since the last harvesting operation - a proxy for forest maturity - then inferred the current statistical distribution of overmature forests (i.e. forests over 50 years without harvesting) in France.

LocationMetropolitan France

MethodsWe used inventories from forest reserves and managed forests to calibrate a generalised linear mixed model explaining the time since the last harvesting with selected structural attributes and environmental variables. We then projected this model on the independent National Forest Inventory dataset. We thus obtained an updated estimation of the proportion and a rough distribution of overmature forest stands in metropolitan France.

ResultsWe found that high basal area of very large trees, high volumes of standing and downed deadwood, high diversity of tree-related microhabitats and more marginally diversity of decay stages best characterized the time since the last harvesting. Volumes of stumps and high density of coppices translating legacy of past forest management also distinguished more overmature plots. Our projection yielded an estimated 3% of French forests over 50 years without harvesting mostly located in more inaccessible areas (i.e. mountainous areas) and a promising proportion of future overmature forests if left unharvested.

Main conclusionsOur study showed that the time since the last harvesting is a good proxy for a combination of stand structure attributes key in characterising overmature temperate forests. It gives the first robust statistical estimate of the proportion of overmature forests and may serve to report on their status in metropolitan France. Our method could be implemented at a larger spatial scale, notably in countries with accessible National Forest Inventory and calibration data, to produce indicators at international level. Omicron BA.1 and BA.2 Variants Increase the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE2 SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. Currently, the most common Omicron variants are the original BA.1 Omicron strain and the BA.2 variant, which became more prevalent since it first appeared. To investigate how these mutations affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ~34% and ~51% stronger binding free energies for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells. 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. Disc and Actin-Associated Protein 1 Influence Attachment in the Intestinal Parasite Giardia lamblia The deep-branching eukaryote Giardia lamblia is an extracellular parasite that attaches to the host intestine via a microtubule-based structure called the ventral disc. Control of attachment is mediated in part by the movement of two regions of the ventral disc that either permit or exclude the passage of fluid under the disc. Several known disc-associated proteins (DAPs) contribute to disc structure and function, but no force-generating protein has been identified among them. We recently identified several Giardia actin (GlActin) interacting proteins at the ventral disc, which could potentially employ actin polymerization for force generation and disc conformational changes. One of these proteins, Disc and Actin Associated Protein 1 (DAAP1), is highly enriched at the two regions of the disc previously shown to be important for fluid flow during attachment. In this study, we investigate the role of both GlActin and DAAP1 in ventral disc morphology and function. We confirmed interaction between GlActin and DAAP1 through coimmunoprecipitation, and used immunofluorescence to localize both proteins throughout the cell cycle and during trophozoite attachment. Similar to other DAPs, the association of DAAP1 with the disc is stable, except during cell division when the disc disassembles. Depletion of GlActin by translation-blocking antisense morpholinos resulted in both impaired attachment and defects in the ventral disc, indicating that GlActin contributes to disc-mediated attachment. Depletion of DAAP1 through CRISPR interference resulted in intact discs but impaired attachment, gating, and flow under the disc. As attachment is essential for infection, elucidation of these and other molecular mediators is a promising area for development of new therapeutics against a ubiquitous parasite.

Author SummaryGiardia lamblia is a single-celled organism and one of the most common gastrointestinal parasites worldwide. In developing countries, recurrent Giardia infections are common, due to lack of access to clean water. Giardia infections can lead to diarrhea, vomiting, dehydration, disruption of the intestinal microbiome, and chronic infections can lead to irritable bowel syndrome. Because existing drug treatments have side effects and Giardias resistance to drugs is increasing, new treatment strategies are needed. The parasites attachment to the hosts intestine is mediated by a Giardia-specific structure that resembles a suction cup and is called the ventral adhesive disc. We previously identified DAAP1, a protein which interacts with Giardia actin and localizes to the ventral disc. Here, we explore the relationship between these two proteins and investigate their role in disc-based attachment. Most disc proteins, including DAAP1, are unrelated to any human proteins, making them appealing drug targets to inhibit parasite attachment and infection. Human readaptation outpaces loss of antibiotic resistance in livestock-associated MRSA Mobile genetic elements (MGEs) are agents of horizontal gene transfer in bacteria, but can also be vertically inherited by daughter cells. Establishing the dynamics that led to contemporary patterns of MGEs in bacterial genomes is central to predicting the emergence and evolution of novel and resistant pathogens. Here we develop existing methods to fully characterise the evolutionary history of MGEs in methicillin-resistant Staphylococcus aureus (MRSA) clonal-complex (CC) 398. CC398 is the dominant MRSA in European livestock, and a growing cause of human infections. We reconstruct MGE dynamics using a collection of 1,180 CC398 genomes, sampled from livestock and humans, over 27 years. We find that the emergence of livestock-associated CC398 coincided with the acquisition of a Tn916 transposon carrying a tetracycline resistance gene, which has been stably vertically inherited for 57 years. This was followed by the acquisition of a type V SCCmec that carries methicillin, tetracycline and heavy metal resistance genes, which has been maintained by vertical and horizontal transmission for 35 years, with occasional truncations and replacements with type IV SCCmec. In contrast, a class of prophages that carry a human immune evasion gene cluster, that are largely absent from livestock-associated CC398, have been repeatedly gained and lost in both human- and livestock-associated CC398. These contrasting dynamics mean that when livestock-associated MRSA infects humans, adaptation to the human host outpaces loss of antibiotic resistance. In addition, stable inheritance means that antibiotic resistance in livestock-associated MRSA may persist despite ongoing reductions in antibiotic and zinc oxide use in farming. Lateralized phenotypic differences after intrahippocampal kainic acid injection in female mice Clinical evidence indicates that patients with temporal lobe epilepsy (TLE) often show differential outcomes of comorbid conditions in relation to the lateralization of the seizure focus. A particularly strong relationship exists between the side of seizure focus and the propensity for distinct reproductive endocrine comorbidities in women with TLE. Therefore, here we evaluated whether targeting of left or right dorsal hippocampus for intrahippocampal kainic acid (IHKA) injection, a model of TLE, produces different outcomes in hippocampal granule cell dispersion, body weight gain, and multiple measures of reproductive endocrine dysfunction in female mice. One, two, and four months after IHKA or saline injection, in vivo measurements of estrous cycles and weight were followed by ex vivo examination of hippocampal dentate granule cell dispersion, circulating ovarian hormone and corticosterone levels, ovarian morphology, and pituitary gene expression. IHKA mice with right-targeted injection (IHKA-R) showed greater granule cell dispersion and pituitary Fshb expression compared to mice with left-targeted injection (IHKA-L). By contrast, pituitary expression of Lhb and Gnrhr were higher in IHKA-L mice compared to IHKA-R, but these values were not different from respective saline-injected controls. IHKA-L mice also showed an increased rate of weight gain compared to IHKA-R mice. Increases in estrous cycle length, however, were similar in both IHKA-L and IHKA-R mice. These findings indicate that although major reproductive endocrine dysfunction phenotypes present similarly after targeting left or right dorsal hippocampus for IHKA injection, distinct underlying mechanisms based on lateralization of epileptogenic insult may contribute to produce similar emergent reproductive endocrine outcomes.

HighlightsO_LIGreater granule cell dispersion with right-sided IHKA injection
C_LIO_LIChanges in pituitary gene expression vary with side of injection
C_LIO_LIIncreased weight gain after left-sided injection
C_LIO_LISimilar estrous cycle disruption after injection of left or right hippocampus
C_LI