Trait-based responses to land use and canopy dynamics modify long-term diversity changes in forest understories AimLand use is the foremost cause of global biodiversity decline, but species do not respond equally to land-use practices. Instead, it is suggested that responses vary with species traits, but long-term data on the trait-mediated effects of land-use on communities is scarce. Here we study how forest understory communities have been affected by two land-use practices during 4-5 decades, and whether changes in plant diversity are related to changes in functional composition.

LocationFinland

Time period1968-2019

Major taxa studiedVascular plants

MethodsWe resurveyed 245 vegetation plots in boreal herb-rich forest understories, and used hierarchical Bayesian linear models to relate changes in diversity, species composition, average plant size, and leaf economic traits to reindeer abundance, forest management intensity, and changes in climate, canopy cover and composition.

ResultsIntensively managed forests decreased in species richness and had increased turnover, but management did not affect functional composition. Increased reindeer densities corresponded with increased leaf dry matter content, evenness and diversity, and decreased height and specific leaf area. Successional development in the canopy was associated with increased specific leaf area and decreased leaf dry matter content and height in the understory over the study period. Effects of reindeer abundance and canopy density on diversity were partially mediated by vegetation height, which had a negative relationship with evenness. Observed changes in climate had no discernible effect on any variable.

Main conclusionsFunctional traits are useful in connecting vegetation changes to the mechanisms that drive them, and provide unique information compared to turnover and diversity metrics. These trait-dependent selection effects could inform whether species benefit or suffer from land use changes and explain observed biodiversity changes under global change.
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.
Feed efficiency of lactating Holstein cows was not as repeatable across diets as within diet over subsequent lactation stages BackgroundImproving feed efficiency has become a common target for dairy farmers to meet the requirement of producing more milk with fewer resources. To improve feed efficiency, a prerequisite is to ensure that the cows identified as most or least efficient will remain as such, independently of diet composition. Therefore, the current research analysed the ability of lactating dairy cows to maintain their feed efficiency while changing the energy density of the diet by changing its concentration in starch and fibre. A total of 60 lactating Holstein cows, including 33 primiparous cows, were first fed a high starch diet-low fibre (diet S+F-), then switched over to a low starch diet-high fibre (diet S-F+). To know if diet affect feed efficiency, we compared the ability of feed efficiency to be maintained within a diet over subsequent lactation stages, known as repeatability of feed efficiency, with its ability to be maintained across diets, known as reproducibility of feed efficiency. To do so we used two indicators: the estimation of the error of repeatability/reproducibility, which is commonly used in metrology, and the coefficient of correlation of concordance (CCC), which is used in biology. The effect of diet change could also lead to a change in cows sorting behaviour which could potentially affect feed efficiency if for example the most efficient cows select more concentrate than the least efficient. We therefore analysed the relationship between the differences in individual feed refusals composition and the differences in feed efficiency. To do so, the composition of each feed refusal was described with its near infra-red (NIR) spectroscopy and was performed on each individual feed ingredient, diet and refusals and used as composition variable. The variability of the NIR spectra of the refusals was described with its principal components thanks to a principal component analysis (PCA). The Pearson correlation was estimated to check the relationship between feed efficiency and refusals composition, i.e. sorting behaviour.

ResultsThe error of reproducibility of feed efficiency across diets was 2.95 MJ/d. This error was significantly larger than the errors of repeatability estimated within diet, which were 2.01 MJ/d within diet S-F+ and 2.40 MJ/d within diet S+F-. The CCC was 0.64 between feed efficiency estimated within diet S+F- and feed efficiency estimated within diet S-F+. This CCC was smaller than the one observed for feed efficiency estimated within diet between two subsequent lactation stages (CCC = 0.72 within diet S+F- and 0.85 within diet S-F+). Feed efficiency was poorly correlated to the first two principal components, which explained 90% of the total variability of the NIR spectra of the individual refusals. This suggests that feed sorting behaviour did not explain differences in feed efficiency.

ConclusionsFeed efficiency was significantly less reproducible across diets than repeatable within the same diet over subsequent lactation stages, but cows ranking for feed efficiency was not significantly affected by diet change. This loss in repeatability across diets could be due to a more pronounced feed sorting subsequent to the change in diet composition. However, the differences in sorting behaviour between cows were not associated to feed efficiency differences in this trial neither with the S+F- diet nor with the S-F+ diet. Those results have to be confirmed on diets having different forage to concentrate ratios to ensure that the least and most efficient cows will not change. Brain-to-brain concordance in child-parent dyads underlies children's psychological wellbeing Negative family environments can be harmful to childrens psychological wellbeing via altered child-parent interactions. Single-brain paradigms have uncovered changes in brain systems related to negative family environments, but how neurobiological reciprocity between child and parent brains contribute to childrens psychological wellbeing remains unknown. Across two studies examining 395 child-parent dyads, we show that children from negative family environments experienced more psychopathological symptoms. Dyad-based analysis of movie-watching functional magnetic resonance imaging data revealed higher inter-subject correlations in dorsal- and ventral-MPFC, and greater concordance of activity with regions critical for socioemotional skills in child-parent than child-stranger dyads. Negative family environment reduced inter-subject correlation of ventral-MPFC with the hippocampus in child-parent dyads, which further mediated childrens internalizing symptoms. The ventral-MPFC and hippocampus exhibited reciprocal activity, which helped child-parent dyads form shared meaning during movie watching. Our findings suggest child-parent brain-to-brain concordance may be a neurobiological mechanism for how family environments affect psychological wellbeing. P2X2 receptor subunit interfaces are missense variant hotspots where mutations tend to increase apparent ATP affinity Background and PurposeP2X receptors (P2XRs) are trimeric ligand-gated ion channels (LGICs) that open a cation-selective pore in response to ATP binding to their large extracellular domain (ECD). The seven known P2XR subtypes can assemble as homo- or heterotrimeric complexes and contribute to numerous physiological functions, including nociception, inflammation and hearing. The overall structure of P2XRs is well established, but little is known about the spectrum and prevalence of human genetic variations and the functional implications in specific domains.

Experimental ApproachHere we examine the impact of P2X2 receptor (P2X2R) inter-subunit interface missense variants identified in the human population or through structural predictions. We test both single and double mutants through electrophysiological and biochemical approaches.

Key resultsWe demonstrate that predicted ECD inter-subunit interfaces display a higher-than-expected density of missense variations and that the majority of mutations that disrupt putative inter-subunit interactions result in channels with higher apparent ATP affinity. Lastly, we show that double mutants at the subunit interface show significant energetic coupling, especially if located in close proximity.

Conclusions and ImplicationsWe provide the first structural mapping of the mutational burden across the human population in a LGIC and show that the density of missense mutations is constrained between protein domains, indicating evolutionary selection at the domain level. Our data may indicate that, unlike other LGICs, P2X2Rs have evolved an intrinsically high threshold for activation, possibly to allow for additional modulation or as a cellular protection mechanism against overstimulation.

Bullet point summary What is already known: O_LIP2X2 receptors are ATP-activated ion channels implicated in hearing and nociceptice pathways
C_LI
 What this study adds: O_LIA structural mapping of missense variants observed in the human population
C_LIO_LIWe identify the intersubunit-interface as a variant hotspot and decipher functional impact of mutations
C_LI
 Clinical significance: O_LIThe development of both inhibitors and activators of P2X2 receptor function may be required
C_LI Intramolecular 13C signals inform about carbon allocation in plants Stable isotopes at natural abundance are key tools to study physiological processes occurring outside the temporal scope of manipulation and monitoring experiments. Whole-molecule carbon isotope ratios (13C/12C) enable assessments of plant carbon uptake yet conceal information about carbon allocation. Here, we identify an intramolecular 13C/12C signal at treering glucose C-5 and C-6 and develop experimentally testable theories on its origin. More specifically, we assess the potential of processes within C3 metabolism for signal introduction based (inter alia) on constraints on signal propagation posed by metabolic networks. We propose that the intramolecular signal reports carbon allocation into major metabolic pathways in actively photosynthesising leaf cells including the anaplerotic, shikimate, and non-mevalonate pathway. We support our theoretical framework by linking it to previously reported whole-molecule 13C/12C increases in cellulose of ozone-treated Betula pendula and a highly significant relationship between the intramolecular signal and tropospheric ozone concentration. Our theory postulates a pronounced preference of leaf-cytosolic triose-phosphate isomerase to catalyse the forward reaction in vivo (dihydroxyacetone phosphate to glyceraldehyde 3-phosphate). In conclusion, intramolecular 13C/12C analysis resolves information about carbon uptake and allocation enabling more comprehensive assessments of carbon metabolism than whole-molecule 13C/12C analysis.

HighlightIntramolecular 13C/12C analysis resolves information about carbon uptake and allocation (and associated environmental controls) enabling more comprehensive assessments of carbon metabolism, plant-environment interactions, and environmental variability than whole-molecule 13C/12C analysis. Hidden dynamics of human proteasome autoregulation during protein degradation The cellular functions are executed by macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving conformational continuum of important complexes at atomic level remains a key challenge in the methodology of computational biophysics. Here we introduce a deep-learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy reconstructions of nonequilibrium conformational continuum. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. Applications of this methodology to six experimental datasets demonstrated its generality in breaking resolution barrier of visualizing dynamic components of functional complexes, in choreographing continuous inter-subunit motions and in exploring their  hidden conformational space. Importantly, AlphaCryo4D identified 64 conformers of the substrate-bound human 26S proteasome and reconstituted single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which together illuminates a grand hierarchical allostery for proteasome autoregulation. 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. Tunable Cellular Localization and Extensive Cytoskeleton-Interplay of Reflectins Reflectin proteins are nature copolymers consist of repeated canonical domains, locating in a biophotonic system called Bragg lamellae and manipulating the dynamic structural coloration of iridocytes. Their biological functions are intriguing, but the underlying mechanism is not fully understood. In this study, reflectin A1, A2, B1 and C were found to present distinguished cyto-/nucleoplasmic localization preferences. Comparable intracellular localization can be reproduced by truncated reflectin variants, suggesting a conceivable evolutionary order among reflectin proteins. Secondly, the size-dependent access of reflectin variants into nucleus demonstrate a potential model of how reflectins get into Bragg lamellae. Thirdly, RfA1 was found to extensively interact with cytoskeleton, including its binding to actin and microtubule minus-end-directed movement. This implies cytoskeleton system plays fundamental role during the organization and transportation of reflectin proteins. Findings presented here recommend reflectins as programmable biomaterials which can be used to decipher their evolution processes, to delineate their biological mechanism, and to achieve tunable intracellular targeting as editable tags. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. Genomic insights into longan evolution from a chromosome-level genome assembly and population genomics of longan accessions Longan (Dimocarpus longan) is a subtropical fruit best known for its nutritious fruit and regarded as a precious tonic and traditional medicine since ancient times. High-quality chromosome-scale genome assembly is valuable for functional genomic study and genetic improvement of longan. Here, we report a chromosome-level reference genome sequence for longan cultivar JDB with an assembled genome of 455.5 Mb in size anchored to fifteen chromosomes, representing a significant improvement of contiguity (contig N50=12.1 Mb, scaffold N50= 29.5 Mb) over a previous draft assembly. A total of 40,420 protein-coding genes were predicted in D. longan genome. Synteny analysis suggests longan shares the widespread gamma event with core eudicots, but has no other whole genome duplications. Comparative genomics showed that D. longan genome experienced significant expansions of gene families related to phenylpropanoid biosynthesis and UDP-glucosyltransferase. Deep genome sequencing analysis of longan cultivars identified longan biogeography as a major contributing factor for genetic diversity, and revealed a clear population admixture and introgression among cultivars of different geographic origins, postulating a likely migration trajectory of longan overall confirmed by existing historical records. Finally, genome-wide association studies (GWAS) of longan cultivars identified quantitative trait loci (QTL) for six different fruit quality traits and revealed a shared QTL containing three genes for total soluble solid and seed weight. The chromosome-level reference genome assembly, annotation and population genetic resource for D. longan will facilitate the molecular studies and breeding of desirable longan cultivars in the future. Use of lignocellulosic residue from second-generation ethanol production to enhance methane production through co-digestion This is a pioneer study evaluating the methane (CH4) production potential from residues of integrated 1st (vinasse and filter cake) and 2nd (deacetylation pretreatment liquor from straw) generation (1G2G) sugarcane biorefinery, providing a fully chemical characterization of them and their relation with the anaerobic digestion (AD) process. Small-scale assays provided fundamentals for basing the co-digestion optimization by assessing the optimal co-substrates synergistic conditions. Biochemical Methane Potential (BMP) tests showed co-digestion enhanced CH4 yield of isolated substrates, reaching up to 605 NmLCH4 gVS-1. The association of vinasse and deacetylation liquor as co-substrates increased the BMP by ~38% mostly by nutritionally benefiting the methanogenic activity. The kinetic analysis confirmed that the deacetylation liquor was the co-substrate responsible for improving the CH4 production in the co-digestion systems due to the highest CH4 conversion rate. The alkaline characteristic of the liquor (pH~12) also prevented alkalizing from being added to the co-digestion, an input that normally makes the process economically unfeasible to implement on an industrial scale due to the large quantities required for buffering the reactor. The filter cake had the lowest BMP (262 NmLCH4 gVS-1) and digestibility ([&le;] 40%), further limited by the required stirring to improve the mass transfer of biochemical reactions. The present study drives towards more sustainable use of vinasse, the most voluminous waste from the sugarcane industry, and lignin-rich residues derived from pre-treatment alkaline methods, aiming at an energy-efficient utilization, by at least 16% when compared to the traditional vinasse AD. The experimental and modeling elements from this work indicated the lignin-rich liquor is the main responsible for putting the co-digestion as a disruptive technological arrangement within the 1G2G sugarcane biorefineries, reinforcing the biogas production as the hub of the bioeconomy in the agroindustrial sector. Structure of the IL-27 quaternary receptor signaling complex Interleukin 27 (IL-27) is a heterodimeric cytokine that functions to constrain T cell-mediated inflammation and plays an important role in immune homeostasis. Binding of IL-27 to cell surface receptors IL-27R and gp130 results in activation of receptor-associated Janus Kinases and nuclear translocation of STAT1 and STAT3 transcription factors. Despite the emerging therapeutic importance of this cytokine axis in cancer and autoimmunity, a molecular blueprint of the IL-27 receptor signaling complex, and its relation to other gp130/IL-12 family cytokines, is currently unclear. We used cryogenic-electron microscopy (cryo-EM) to determine the quaternary structure of IL-27 (p28/Ebi3) bound to receptor subunits, IL-27R and gp130. The resulting 3.47 [A] resolution structure revealed a three-site assembly mechanism nucleated by the central p28 subunit of the cytokine. The overall topology and molecular details of this binding are reminiscent of IL-6 but distinct from related heterodimeric cytokines IL-12 and IL-23. These results indicate distinct receptor assembly mechanisms used by heterodimeric cytokines with important consequences for targeted agonism and antagonism of IL-27 signaling. Modulations in motor unit discharge are related to changes in fascicle length during isometric contractions The integration of electromyography (EMG) and ultrasound imaging has provided important information about the mechanisms of muscle activation and contraction. Unfortunately, EMG does not allow an accurate assessment of the interplay between the neural drive received by muscles, changes in fascicle length (FL) and the force/torque produced. We aimed to assess the relationship between modulations in tibialis anterior (TA) motor unit (MU) firing rate, FL and dorsiflexion torque (DT) using ultrasound-transparent high-density EMG electrodes. EMG and ultrasound images were recorded simultaneously from TA, using a 32-electrode silicon matrix, while performing isometric dorsiflexion, at diverse ankle joint positions (0{degrees} and 30{degrees} plantar flexion) and torques (20% and 40% of maximum). EMG signals were decomposed into individual MUs and changes in FL were assessed with a fascicle-tracking algorithm. MU firings were converted into a cumulative spike train (CST) that was cross-correlated with DT (CST-DT) and FL (CST-FL). High cross-correlations were found for CST-FL, 0.60 (range: 0.31-0.85) and CST-DT 0.71 (range: 0.31-0.88). Cross-correlation lags revealed that the delay between CST-FL (~75ms) was significantly smaller than CST-DT (~150ms, p<0.001). These delays affected the interpretation of MU recruitment/de-recruitment thresholds, with FL showing similar lengths for both recruitment and de-recruitment. This study is the first to demonstrate that changes in TA FL are closely related to both modulations in MU firing frequency and DT. These relationships allow assessment of the interplay between neural drive, muscle contraction and resultant torque, thereby providing a better understanding of the mechanisms responsible for the generation of muscle force.

NEW AND NOTEWORTHYBy employing ultrasound-transparent high-density surface EMG electrodes, we show that modulations in tibialis anterior motor unit discharge rate were closely related to both changes in its fascicle length and resultant torque. These relationships allowed quantifying delays between neural drive and muscle shortening as well as muscle shortening and torque during submaximal isometric contractions, providing an accurate estimation of the time required to generate muscle force and subsequent production of torque via the tendon. Naturally occurring combinations of receptors from single cell transcriptomics in endothelial cells VEGF inhibitor drugs have been successful, especially in ophthalmology, but not all patients respond to them. Combinations of drugs are likely to be needed for more effective therapies of angiogenesis-related diseases. In this paper we describe naturally occurring combinations of receptors in endothelial cells that might help to understand how cells communicate and to identify targets for drug combinations. We also develop and share a new software tool called DECNEO to identify them.

Single-cell gene expression data are used to identify a set of co-expressed endothelial cell receptors, conserved among species (mice and humans) and enriched, within a network, of connections to up-regulated genes. This set includes several receptors previously shown to play a role in angiogenesis. Multiple statistical tests from large datasets, including an independent validation set, support the reproducibility, evolutionary conservation and role in angiogenesis of these naturally occurring combinations of receptors. We also show tissue-specific combinations and, in the case of choroid endothelial cells, consistency with both well-established and recent experimental findings, presented in a separate paper.

The results and methods presented here advance the understanding of signaling to endothelial cells. The methods are generally applicable to the decoding of intercellular combinations of signals. Tackling polypharmacology of kinase inhibitors and selecting chemical kinase probes by transcription factor activity profiling. Protein kinase inhibitors (PKI) evolved as promising drugs for multiple diseases. However, PKIs promiscuity causes polypharmacological effects challenging therapeutic development. We show that the polypharmacology of PKIs can be inferred from their impact on transcription factors (TF) lying at the apexes of signaling pathways. Using a high-content reporter assay, we evaluated TF activity profiles (TFAP) in cells treated with Akt, CDK, Aurora, Raf, MEK, and ERK inhibitors. For each kinase family, we found kinase-specific consensus TFAP signatures signifying the on-target PKI activity. Remarkably, proximal kinases had the highest similarity of their PKIs consensus signatures. Furthermore, we show that individual PKIs exhibited the consensus, "on-target signatures" within certain "specificity windows" and distinct "off-target signatures" at other concentrations. We also show that the TFAP signatures permit pinpointing the off-target PKI activity at irrelevant kinases and non-kinase effectors. Therefore, the TFAP approach provides clear-cut quantitative metrics for assessing PKI activity and concentration ranges where the on-target activity dominates cell response. By illuminating PKI biology invisible to target-based techniques, this effect-based approach streamlines the selection of kinase chemical probes and PKI drug prioritization. Combined Lineage Tracing and scRNA-seq Reveals Unexpected First Heart Field Predominance of Human iPSC Differentiation During mammalian development, the left and right ventricles arise from early populations of cardiac progenitors known as the first and second heart fields, respectively. While these populations have been extensively studied in non-human model systems, their identification and study in vivo human tissues have been limited due to the ethical and technical limitations of accessing gastrulation stage human embryos. Human induced pluripotent stem cells (hiPSCs) present an exciting alternative for modeling early human embryogenesis due to their well-established ability to differentiate into all embryonic germ layers. Here, we describe the development of a TBX5/MYL2 lineage tracing reporter system that allows for the identification of FHF-progenitors and their descendants including left ventricular cardiomyocytes. Furthermore, using single cell RNA sequencing (scRNA-seq) with oligonucleotide-based sample multiplexing, we extensively profiled differentiating hiPSCs across 12 timepoints in two independent iPSC lines. Surprisingly, our reporter system and scRNA-seq analysis revealed a predominance of FHF differentiation using the small molecule Wnt-based 2D differentiation protocol. We compared this data with existing murine and 3D cardiac organoid scRNA-seq data and confirmed the dominance of left ventricular cardiomyocytes (>90%) in our hiPSC-derived progeny. Together, our work provides the scientific community with a powerful new genetic lineage tracing approach as well as a single cell transcriptomic atlas of hiPSCs undergoing cardiac differentiation. cvlr: Finding heterogeneously methylated genomicregions using ONT reads SummaryNanopore reads encode information on the methylation status of cytosines in CpG dinucleotides. The length of the reads makes it comparatively easy to look at patterns consisting of multiple loci; here we exploit this property to look for regions where one can define subpopulations of cells based on methylation patterns. As a benchmark we run our clustering algorithm on known imprinted genes and show that the clustering based on methylation is consistent with the phasing of the genome; we then scan chromosome 15 looking for windows corresponding to heterogeneous methylation. We can also compute the covariance of methylation across these regions while keeping into account the mixture of different types of reads.

Availabilityhttps://github.com/EmanueleRaineri/releases

Contactemanuele.raineri@cnag.crg.eu, simon.heath@cnag.crg.eu

Supplementary informationTables, figures, and some further explanations of the algorithms are available as online supplementary information. Structure-function analysis of the nsp14 N7-guanine methyltransferase reveals an essential role in Betacoronavirus replication As coronaviruses (CoVs) replicate in the host cell cytoplasm, they rely on their own capping machinery to ensure the efficient translation of their mRNAs, protect them from degradation by cellular 5 exoribonucleases, and escape innate immune sensing. The CoV nonstructural protein 14 (nsp14) is a bi-functional replicase subunit harboring an N-terminal 3'-to-5' exoribonuclease (ExoN) domain and a C-terminal (N7-guanine)-methyltransferase (N7-MTase) domain that is presumably involved in viral mRNA capping. Here, we aimed to integrate structural, biochemical, and virological data to assess the importance of conserved N7-MTase residues for nsp14s enzymatic activities and virus viability. We revisited the crystal structure of severe acute respiratory syndrome (SARS)-CoV nsp14 to perform an in silico comparative analysis between betacoronaviruses. We identified several residues likely involved in the formation of the N7-MTase catalytic pocket, which presents a fold distinct from the Rossmann fold observed in most known MTases. Next, for SARS-CoV and Middle East respiratory syndrome-CoV, site-directed mutagenesis of selected residues was used to assess their importance for in vitro enzymatic activity. Most of the engineered mutations abolished N7-MTase activity, while not affecting nsp14-ExoN activity. Upon reverse engineering of these mutations into different betacoronavirus genomes, we identified two substitutions (R310A and F426A in SARS-CoV nsp14) abrogating virus viability and one mutation (H424A) yielding a crippled phenotype across all viruses tested. Our results identify the N7-MTase as a critical enzyme for betacoronavirus replication and define key residues of its catalytic pocket that can be targeted to design inhibitors with a potential pan-coronaviral activity spectrum.

Significance StatementThe ongoing SARS-CoV-2 pandemic emphasizes the urgent need to develop efficient broad-spectrum anti-CoV drugs. The structure-function characterization of conserved CoV replicative enzymes is key to identifying the most suitable drug targets. Using a multidisciplinary comparative approach and different betacoronaviruses, we characterized the key conserved residues of the nsp14 (N7-guanine)-methyltransferase, a poorly defined subunit of the CoV mRNA-synthesizing machinery. Our study highlights the unique structural features of this enzyme and establishes its essential role in betacoronavirus replication, while identifying two residues that are critical for the replication of the four betacoronaviruses tested, including SARS-CoV-2. Batch experiments demonstrating a two-stage bacterial process coupling methanotrophic and heterotrophic bacteria for 1-alkene production from methane Methane (CH4) is a sustainable carbon feedstock source for aerobic CH4-oxidizing bacteria (methanotrophs) to produce value-added chemicals. Under substrate-limited (e.g., CH4, oxygen and nitrogen) conditions, CH4 oxidation results in the production of various short-chain organic acids and platform chemicals. These CH4-derived products could be broadened by utilizing them as a feedstock for heterotrophic bacteria. As a proof of concept, a two-stage system for CH4 abatement and 1-alkene production was developed in this study. Types I and II methanotrophs, i.e., Methylobacter tundripaludum SV96 and Methylocystis rosea SV97, respectively, were investigated in batch tests under different CH4 and air supplementation schemes. CH4 oxidation under either microaerobic or aerobic conditions induced the production of formate, acetate, succinate, and malate in M. tundripaludum SV96, accounting for 4.8-7.0% of consumed CH4-carbon while M. rosea SV97 produced the same compounds except for malate, and with lower efficiency than M. tundripaludum SV96, accounting for 0.7-1.8% of consumed CH4-carbon For the first time, the organic acids-rich spent media of methanotrophs were successfully used for 1-alkene production using engineered Acinetobacter baylyi ADP1  tesA-undA cells. The highest yield of 1-undecene was obtained from spent medium of M. tundripaludum SV96 at 68.9 {+/-} 11.6 mol mol Csubstrate-1. LPS-TLR4 pathway exaggerates alcohol-induced acute-on-chronic liver failure via provoking NETs formation Backgrounds: Intrahepatic infiltration of neutrophils is a character of alcoholic acute-on-chronic liver failure (AACLF) and neutrophil extracellular traps (NETs) are an important strategy for neutrophils to fix and kill invading microorganisms. Intestinal bacteria and the gut-liver axis have been thought to play a key role in many liver diseases also including AACLF. However, whether NETs appear in AACLF and play a role in AACLF is still unsure. Methods: WT, NE KO, and TLR4 KO mice were used to build the AACLF model, and the intestinal bacteria were eliminated at the same time and LPS was given. Then the formation of NETs and AACLF related markers were detected. Results: The serum MPO-DNA and LPS concentration was increased in AACLF patients and a correlation was revealed between these two indexes. More intrahepatic NETs formed in AACLF mice by testing MPO-DNA, Cit H3, and NE. These markers decreased with gut detergent and restored markers with gut detergent plus LPS supplement. While NETs formation failed to change with gut microbiome or combine LPS supplement in TLR4 KO mice. As we tested AACLF related characters, liver injury, intrahepatic fat deposition, inflammation, and fibrosis alleviated with depletion of NE. These related marks were also attenuated with gut sterilization by antibiotics and recovered with combined treatment with antibiotics plus LPS. But the liver injury, intrahepatic fat, fibro deposition, and liver inflammation-related markers did show a significant difference in TLR4 KO mice when they received the same treatment. Conclusion: Intestinal-derived LPS promotes NETs formation in AACLF through the TLR4 pathway and further accelerates the AACLF process by NETs. A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Spatial phylogenetics of Japanese ferns: Patterns, processes, and implications for conservation PremiseBiodiversity is often only measured with species richness. However, this ignores evolutionary history and is not sufficient for making conservation decisions. Here, we characterize multiple facets and drivers of biodiversity to understand how these relate to bioregions and conservation status in the ferns of Japan.

MethodsWe compiled a community dataset of 1,239 20 km x 20 km grid-cells including 672 taxa based on > 300,000 specimen records. We combined this with a phylogeny and functional traits to analyze taxonomic, phylogenetic, and functional diversity, and modeled biodiversity metrics in response to environmental factors and reproductive mode. Hierarchical clustering was used to delimit bioregions. Conservation status and threats were assessed by comparing the overlap of significantly diverse grid-cells with conservation zones and range maps of native Japanese deer.

ResultsTaxonomic richness was highest at mid-latitudes. Phylogenetic and functional diversity and phylogenetic endemism were highest in small southern islands. Relative phylogenetic and functional diversity were high at high and low latitudes, and low at mid-latitudes. Grid-cells were grouped into three (phylogenetic) or four (taxonomic) major bioregions. Temperature and apomixis were identified as drivers of biodiversity patterns. Conservation status was generally high for grid-cells with significantly high biodiversity, but the threat due to herbivory by deer was greater for taxonomic richness than other metrics.

ConclusionsOur integrative approach reveals previously undetected patterns and drivers of biodiversity in the ferns of Japan. Future conservation efforts should recognize that threats can vary by biodiversity metric and consider multiple multiple metrics when establishing conservation priorities. Outsourced hearing in an orb-weaving spider that uses its web as an auditory sensor Hearing is a fundamental sense of many animals, including all mammals, birds, some reptiles, amphibians, fish, and arthropods. The auditory organs of these animals are extremely diverse in anatomy after hundreds of millions of years of evolution, yet all are made up of cellular tissues and are morphologically part of bodies of animals. Here we show hearing in the orb-weaving spider, Larinioides sclopetarius is not constrained by the organisms body but is extended through outsourcing hearing to its extended phenotype, the proteinaceous, self-manufactured orb-web. We find the wispy, wheel-shaped orb-web acts as a hyperacute acoustic "antenna" to capture the sound-induced air particle movements that approach the maximum physical efficiency, better than the acoustic responsivity of all previously known eardrums. By sensing the motion of web threads, the spider remotely detects and localizes the source of an incoming airborne acoustic wave such as those emitted by approaching prey or predators. By outsourcing its acoustic sensors to its web, the spider is released from body size constraints and permits the araneid spider to increase its sound-sensitive surface area enormously, up to 10,000 times greater than the spider itself. The spider also enables the flexibility to functionally adjust and regularly regenerate its "external eardrum" according to its needs. The "outsourcing" and "supersizing" of auditory function in spiders provides unique features for studying extended and regenerative sensing, and designing novel acoustic flow detectors for precise fluid dynamic measurement and manipulation. A platform of patient-derived microtumors identifies therapeutic vulnerabilities in ovarian cancer. BackgroundIn light of the frequent development of therapeutic resistance in cancer treatment, there is a strong need for personalized model systems representing patient tumor heterogeneity, while enabling parallel drug testing and identification of appropriate treatment responses in individual patients. Using ovarian cancer as a prime example of a heterogeneous tumor disease with complex microenvironment and high recurrence rates, we developed a 3D preclinical ovarian cancer model comprised of patient-derived microtumors (PDM) and autologous tumor-infiltrating lymphocytes (TILs) for identification of treatment vulnerabilities and validation of anti-cancer drug efficacy using immunohistochemistry, immune cell phenotyping, functional assays and protein profiling analyses.

MethodsPDM and TILs were isolated from fresh primary ovarian cancer tissue specimen using mechanical disruption and limited enzymatic digestion and were subsequently cultured in suspension in defined media in the absence of serum. The heterogeneous cellular composition of isolated PDM as well as autologous TILs was analyzed by FFPE immunohistochemistry and multi-color flow cytometry, respectively. For in-depth protein profiling of PDM we established Reverse Phase Protein Array (RPPA) analyses of >110 total and phospho-proteins. Treatment efficacy in response to chemotherapeutics as well as immunotherapeutic compounds was assessed in PDM and PDM-TIL co-cultures using a functional viability assay in microplate format.

ResultsThe enzymatic digestion of primary ovarian cancer tissue and suspension culture in defined serum-free media allowed fast and efficient recovery of patient-derived microtumors (PDM). Immunohistochemical analyses demonstrated histopathological comparability of ovarian cancer PDM with corresponding patient tumor tissue. Reverse Phase Protein Array (RPPA)-based analyses of >110 total and phospho-proteins enabled the identification of patient-specific sensitivities to standard, platinum-based therapy and thereby the prediction of potential treatment-responders. Finally, combining PDM and autologous TILs for individual efficacy testing of immune checkpoint inhibitor treatment demonstrated the potential for patient-specific enhancement of cytotoxic TIL activity by this therapeutic approach.

Conclusion3D patient-derived ovarian cancer microtumors represent a preclinical, ex vivo tumor model that reflects intertumoral heterogeneity and represent the cellular complexity of individual patient tumors. Combining protein pathway analysis and anti-cancer drug efficacy testing of PDM enables drug mode-of-action analyses and therapeutic sensitivity prediction within a clinically relevant time frame after surgery. Follow-up studies in larger cohorts are currently under way to further evaluate the applicability of this platform to support clinical decision-making and personalizing cancer treatment.