Mapping the planet's critical natural assets Sustaining the organisms, ecosystems, and processes that underpin human well-being is necessary to achieve sustainable development. Here we identify critical natural assets, natural and semi-natural ecosystems that provide 90% of the total current magnitude of 14 types of natures contributions to people (NCP). Critical natural assets for maintaining local-scale NCP (12 of the 14 NCP mapped) comprise 30% of total global land area and 24% of national territorial waters, while 44% of land area is required for maintaining all NCP (including those that accrue at the global scale, carbon storage and moisture recycling). At least 87% of the worlds population lives in the areas benefiting from critical natural assets for local-scale NCP, while only 16% lives on the lands containing these assets. Critical natural assets also overlap substantially with areas important for biodiversity (covering area requirements for 73% of birds and 66% of mammals) and cultural diversity (representing 96% of global Indigenous and non-migrant languages). Many of the NCP mapped here are left out of international agreements focused on conserving species or mitigating climate change, yet this analysis shows that explicitly prioritizing critical natural assets for NCP could simultaneously advance development, climate, and conservation goals. Crafting policy and investment strategies that protect critical natural assets is essential for sustaining human well-being and securing Earths life support systems.
Host species identity shapes the diversity and structure of insect microbiota As for most of the life that inhabits our planet, microorganisms play an essential role in the fitness of insects, including nutrition, reproduction, defence, and many other functions. More recently, we assisted in an exponential growth of studies describing the taxonomical composition of bacterial communities across insects phylogeny. However, there is still an outstanding question that needs to be answered: which factors contribute most to shape insects microbiomes? This study tries to find an answer to this question by taking advantage of publicly available sequencing data and reanalysing over 4,000 samples of insect-associated bacterial communities under a common framework. Results suggest that insect taxonomy has a wider impact on the structure and diversity of their associated microbial communities than the other factors considered (diet, sex, life stage, sample origin, and treatment). However, when specifically testing for signatures of co-diversification of insect species and their microbiota, analyses found weak support for this, suggesting that while insect species strongly drive the structure and diversity of insect microbiota, the diversification of those microbial communities did not follow their hosts phylogeny. Furthermore, a parallel survey of the literature highlights several methodological limitations that need to be considered in future research endeavours.
A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Deficiency of gluconeogenic enzyme PCK1 promotes murine NASH progression and fibrosis by activation PI3K/AKT/PDGF axis Background and AimsNonalcoholic steatohepatitis (NASH) is a chronic liver disease characterized by hepatic lipid accumulation, inflammation, and progressive fibrosis. However, the pathomechanisms underlying NASH are incompletely explored. Phosphoenolpyruvate carboxykinase 1 (PCK1) catalyzes the first rate-limiting step of gluconeogenesis. This study was designed to determine the role of PCK1 in regulating NASH progression.

MethodsLiver metabolism, hepatic steatosis, and fibrosis were evaluated at 24 weeks in liver-specific Pck1-knockout mice fed with NASH diet or chow diet. Gain- or loss-of-function approaches were used to explore the underlying mechanism in vitro. AKT and RhoA inhibitors were evaluated for NASH treatment in vivo.

ResultsHepatic PCK1 was downregulated in patients with NASH and mouse models of NASH. Mice with liver Pck1 deficiency displayed hepatic lipid disorder and liver injury fed with normal diet, while showed aggravated fibrosis and inflammation when fed NASH diet. Mechanistically, PCK1 deficiency upregulated genes involved in fatty acid transport and lipid droplet formation. Moreover, metabolomics analysis showed the accumulation of glycerol 3-phosphate, the substrate of triglyceride synthesis. Furthermore, the loss of hepatic PCK1 could activate the RhoA/PI3K/AKT pathway, which leads to increased secretion of PDGF-AA and promotes the activation of hepatic stellate cells. Accordingly, treatment with RhoA and AKT inhibitors alleviated NASH progression in the presence of Pck1 deletion in vivo.

ConclusionsPCK1 deficiency plays a key role in the development of hepatic steatosis and fibrosis by facilitating the RhoA/PI3K/AKT/PDGF-AA axis. These findings provide a novel insight into therapeutic approaches for the treatment of NASH. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. Absence of cardiolipin from the outer-leaflet of a mitochondrial inner-membrane mimic restricts Opa1-mediated fusion Cardiolipin is a tetra-acylated di-phosphatidylglycerol lipid enriched in the matrix-facing (inner) leaflet of the mitochondrial inner membrane. Cardiolipin plays an important role in regulating mitochondria function and dynamics. Yet, the mechanisms connecting cardiolipin distribution and mitochondrial protein function remain indirect. In our previous work, we established an in vitro system reconstituting mitochondrial inner membrane fusion mediated by Opa1. We found that the long form of Opa1 (l-Opa1) works together with the proteolytically processed short form (s-Opa1) to mediate fast and efficient membrane fusion. Here, we extend our reconstitution system to generate supported lipid bilayers with asymmetric cardiolipin distribution. Using this system, we find the presence of cardiolipin on the inter-membrane space-facing (outer) leaflet is important for membrane tethering and fusion. We discuss how the presence of cardiolipin in this leaflet may influence protein and membrane properties, and future applications for this approach. 5S-IGS rDNA in wind-pollinated trees ( Fagus  L.) encapsulates 55 million years of reticulate evolution and hybrid origins of modern species Standard models of plant speciation assume strictly dichotomous genealogies in which a species, the ancestor, is replaced by two offspring species. The reality in wind-pollinated trees with long evolutionary histories is more complex: species evolve from other species through isolation when genetic drift exceeds gene flow; lineage mixing can give rise to new species (hybrid taxa such as nothospecies and allopolyploids). The multi-copy, potentially multi-locus 5S rDNA is one of few gene regions conserving signal from dichotomous and reticulate evolutionary processes down to the level of intra-genomic recombination. Therefore, it can provide unique insights into the dynamic speciation processes of lineages that diversified tens of millions of years ago. Here, we provide the first high-throughput sequencing (HTS) of the 5S intergenic spacers (5S-IGS) for a lineage of wind-pollinated subtropical to temperate trees, the Fagus crenata - F. sylvatica s.l. lineage, and its distant relative F. japonica. The observed 4,963 unique 5S-IGS variants reflect a complex history of hybrid origins, lineage sorting, mixing via secondary gene flow, and intra-genomic competition between two or more paralogous-homoeologous 5S rDNA lineages. We show that modern species are genetic mosaics and represent a striking case of ongoing reticulate evolution during the past 55 million years.

Significance statementThe evolution of extra-tropical wind-pollinated tree genera involves dynamic speciation processes. High-throughput sequencing of the multi-copy, potentially multi-locus 5S rDNA reveals a complex history of hybrid origins, lineage sorting and mixing, and intra-genomic competition between paralogous-homeologous loci in the core group of Eurasian beech trees (genus Fagus) and their distant relative, F. japonica. The modern species are genetic mosaics and represent a striking case of at least 55 million years of ongoing reticulate evolution. De-noising Spatial Transcriptomics Data Based on Position and Image Information Spatial transcriptomics (ST) technologies provide gene expression close to or even superior to single-cell resolution while retaining the physical locations of sequencing and often also providing matched pathology images. However, the expression data captured by ST technologies suffer from high noise levels, as a result of the shallow coverage in each sequencing unit. The extra experimental steps for preserving the spatial locations of sequencing could result in even more severe noises, compared to regular single-cell RNA-sequencing (scRNA-seq). Fortunately, such noises could be largely removed by leveraging information from the physical locations of sequencing, and the tissue and cellular organization reflected by corresponding pathology images. In this work, we demonstrated the extensive levels of noise in ST data. We developed a mathematical model, named Sprod, to impute accurate ST gene expression based on latent space and graph learning of matched location and imaging data. We comprehensively validated Sprod and demonstrated its advantages over prior methods for removing drop-outs in scRNA-seq data. We further showed that, after adequate imputation by Sprod, differential expression analyses, pseudotime analyses, and cell-to-cell interaction inferences yield significantly more informative results. Overall, we envision denoising by Sprod to become a key first step to empower ST technologies for biomedical discoveries and innovations. Community Evaluation of Glycoproteomics Informatics Solutions Reveals High-Performance Search Strategies of Serum N- and O-Glycopeptide Data Glycoproteome profiling (glycoproteomics) is a powerful yet analytically challenging research tool. The complex tandem mass spectra generated from glycopeptide mixtures require sophisticated analysis pipelines for structural determination. Diverse software aiding the process have appeared, but their relative performance remains untested. Conducted through the HUPO Human Proteome Project - Human Glycoproteomics Initiative, this community study, comprising both developers and users of glycoproteomics software, evaluates the performance of informatics solutions for system-wide glycopeptide analysis. Mass spectrometry-based glycoproteomics datasets from human serum were shared with all teams. The relative team performance for N- and O-glycopeptide data analysis was comprehensively established and validated through orthogonal performance tests. Excitingly, several high-performance glycoproteomics informatics solutions were identified. While the study illustrated that significant informatics challenges remain, as indicated by a high discordance between annotated glycopeptides, lists of high-confidence (consensus) glycopeptides were compiled from the standardised team reports. Deep analysis of the performance data revealed key performance-associated search variables and led to recommendations for improved "high coverage" and "high accuracy" glycoproteomics search strategies. This study concludes that diverse software for comprehensive glycopeptide data analysis exist, points to several high-performance search strategies, and specifies key variables that may guide future software developments and assist informatics decision-making in glycoproteomics. Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

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

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

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach. Patient-specific network connectivity combined with a next generation neural mass model to test clinical hypothesis of seizurepropagation Dynamics underlying epileptic seizures span multiple scales in space and time, therefore, understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. In this view, mathematical models have been developed, ranging from single neuron to neural population.

In this study we consider a neural mass model able to exactly reproduce the dynamics of heterogeneous spiking neural networks. We combine the mathematical modelling with structural information from non-invasive brain imaging, thus building large-scale brain network models to explore emergent dynamics and test clinical hypothesis. We provide a comprehensive study on the effect of external drives on neuronal networks exhibiting multistability, in order to investigate the role played by the neuroanatomical connectivity matrices in shaping the emergent dynamics. In particular we systematically investigate the conditions under which the network displays a transition from a low activity regime to a high activity state, which we identify with a seizure-like event. This approach allows us to study the biophysical parameters and variables leading to multiple recruitment events at the network level. We further exploit topological network measures in order to explain the differences and the analogies among the subjects and their brain regions, in showing recruitment events at different parameter values.

We demonstrate, along the example of diffusion-weighted magnetic resonance imaging (MRI) connectomes of 20 healthy subjects and 15 epileptic patients, that individual variations in structural connectivity, when linked with mathematical dynamic models, have the capacity to explain changes in spatiotemporal organization of brain dynamics, as observed in network-based brain disorders. In particular, for epileptic patients, by means of the integration of the clinical hypotheses on the epileptogenic zone (EZ), i.e. the local network where highly synchronous seizures originate, we have identified the sequence of recruitment events and discussed their links with the topological properties of the specific connectomes. The predictions made on the basis of the implemented set of exact mean-field equations turn out to be in line with the clinical pre-surgical evaluation on recruited secondary networks. Integrative Modelling of Signalling Network Dynamics Identifies Cell Type-selective Therapeutic Strategies for FGFR4-driven Cancers Oncogenic FGFR4 signalling represents a potential therapeutic target in many cancer types, including triple negative breast cancer (TNBC) and hepatocellular carcinoma (HCC). However, resistance to single-agent therapy directed at FGFR4 remains a major challenge, prompting the need to identify more effective combinatorial therapeutic strategies. Here, we integrated computational network modelling and experimental validation to characterise dynamic reprogramming of the FGFR4 signalling network in TNBC following FGFR4 kinase inhibition. We found that AKT, which signals downstream of FGFR4, displayed a rapid and potent reactivation following FGFR4 targeting. Through model-based simulation and systematic prediction of the effect of co-targeting specific network nodes, we predicted, and validated experimentally, strong synergism of co-targeting FGFR4 and particular ErbB kinases or AKT, but not the upstream kinase PI3K. Further, incorporation of protein expression data from hundreds of cancer cell types enabled us to adapt our model to other diverse cellular contexts, leading to the prediction that while AKT rebound occurs frequently, it is not a general phenomenon. Instead, ERK is reactivated in a subset of cell types, including the FGFR4-driven HCC cell line Hep3B. This was subsequently corroborated, and moreover, co-targeting FGFR4 and MEK in Hep3B cells markedly enhanced inhibition of cell proliferation. Overall, these findings provide novel insights into the dynamics of drug-induced network remodelling in cancer cells, highlight the impact of protein expression heterogeneity on network response to targeted therapy and identify candidate cell type-selective combination treatments for FGFR4-driven cancer. Atrazine induced In vivo immunotoxicity in Bivalves Atrazine is ubiquitously used broad-spectrum herbicide to control the weeds in agriculture. The present study aimed to evaluate the acute toxicity and immunotoxicity of Atrazine in two ecologically and economically important bivalves. Acute toxicity of atrazine evaluated in triplicates by taking control and six experimental groups each comprising of 30 animals and treated with a range of atrazine from 2 PPM to 12 PPM for 96 hours. Mortalities were recorded for every 24 hours until 96 hours and data analyzed by one-way ANOVA and Dunnett T-test. The results indicated a significant increase in mortalities with increase in dose and time of exposure in both species. The values of LC50 were determined as 6.10 PPM and 4.90 PPM respectively for Perna viridis and Paphia malabarica. Furthermore, the immunotoxic potential of atrazine assessed by treating mussels and clams with the five sub-lethal doses of atrazine for 14 days and quantifying the viability of hemocytes by using simple yet reliable Tryphan blue exclusion assay. The results of the present study suggest atrazine-induced immunotoxicity by decreasing the number of viable hemocytes in bivalves. Hemocytes with phagocytic function are indispensable to confer innate immunity in bivalves, decreased viability of these cells leads to compromised immunity. This study is first of its kind to implicate atrazine with the immunotoxicity in bivalves. The DEAD box RNA helicase DDX42 is an intrinsic inhibitor of positive-strand RNA viruses Genome-wide screens are powerful approaches to unravel new regulators of viral infections. Here, we used a CRISPR/Cas9 screen to reveal new HIV-1 inhibitors. This approach led us to identify the RNA helicase DDX42 as an intrinsic antiviral inhibitor. DDX42 was previously described as a non-processive helicase, able to bind RNA secondary structures such as G-quadruplexes, with no clearly defined function ascribed. Our data show that depletion of endogenous DDX42 significantly increased HIV-1 DNA accumulation and infection in cell lines and primary cells. DDX42 overexpression inhibited HIV-1, whereas a dominant-negative mutant increased infection. Importantly, DDX42 also restricted retrotransposition of LINE-1, infection with other retroviruses and positive-strand RNA viruses, including CHIKV and SARS-CoV-2. However, DDX42 did not inhibit infection with three negative-strand RNA viruses, arguing against a general, unspecific effect on target cells, which was confirmed by RNA-seq analysis. DDX42 was found in the vicinity of viral elements by proximity ligation assays, and cross-linking RNA immunoprecipitation confirmed a specific interaction of DDX42 with RNAs from sensitive viruses. This strongly suggested a direct mode of action of DDX42 on viral ribonucleoprotein complexes. Taken together, our results show for the first time a new and important role of DDX42 in intrinsic antiviral immunity. Optimized detection of allelic imbalances specific for homologous recombination deficiency improves the prediction of clinical outcomes in cancer Homologous recombination DNA-repair deficiency (HRD) is a common driver of genomic instability and confers a therapeutic vulnerability in human cancer. Using a systematic pan- cancer characterization, we revealed distinct patterns of somatic allelic imbalances (AIs) in high-grade serous ovarian cancer (HGSC). Using machine learning on a multi-omics dataset, we generated an optimized algorithm to detect HRD in HGSC (ovaHRDscar). ovaHRDscar improved the prediction of clinical outcomes in three independent validation cohorts. Characterization of 98 spatiotemporally distinct tumor samples indicated ovary/adnex as the preferred site to assess HRD. Cancer-specific optimization improved the prediction of clinical outcomes also in triple-negative breast cancer (tnbcHRDscar). In conclusion, our systematic approach based on multi-omics data improves the detection of HRD with the premise to improve patient selection for HR-targeted therapies. Mitochondrial genomes of two parasitic Cuscuta species lack clear evidence of horizontal gene transfer and retain unusually fragmented ccmFC genes BackgroundThe intimate association between parasitic plants and their hosts favours the exchange of genetic material, potentially leading to horizontal gene transfer (HGT) between plants. With the recent publication of several parasitic plant nuclear genomes, there has been considerable focus on such non-sexual exchange of genes. To enhance the picture on HGT events in a widely distributed parasitic genus, Cuscuta (dodders), we assembled and analyzed the organellar genomes of two recently sequenced species, C. australis and C. campestris, making this the first account of complete mitochondrial genomes (mitogenomes) for this genus.

ResultsThe mitogenomes are 265,696 and 275,898 bp in length and contain a typical set of mitochondrial genes, with ten missing or pseudogenized genes often lost from angiosperm mitogenomes. Each mitogenome also possesses a structurally unusual ccmFC gene, which exhibits splitting of one exon and a shift to trans-splicing of its intron. Based on phylogenetic analysis of mitochondrial genes from across angiosperms and similarity-based searches, there is little to no indication of HGT into the Cuscuta mitogenomes. A few candidate regions for plastome-to-mitogenome transfer were identified, with one suggestive of possible HGT.

ConclusionsThe lack of HGT is surprising given examples from the nuclear genomes, and may be due in part to the relatively small size of the Cuscuta mitogenomes, limiting the capacity to integrate foreign sequences. 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. A dysmorphic mouse model reveals developmental interactions of chondrocranium and dermatocranium The cranial endo- and dermal skeletons, which comprise the vertebrate skull, evolved independently over 470 million years ago and form separately during embryogenesis. In mammals, much of the cartilaginous chondrocranium is transient, undergoing endochondral ossification or disappearing, so its role in skull morphogenesis is not well studied and it remains an enigmatic structure. We provide complete three-dimensional (3D) reconstructions of the laboratory mouse chondrocranium from embryonic day 13.5 through 17.5 using a novel methodology of uncertainty-guided segmentation of phosphotungstic enhanced 3D microcomputed tomography images with sparse annotation. We evaluate the embryonic mouse chondrocranium and dermatocranium in 3D and delineate the effects of a Fgfr2 variant on embryonic chondrocranial cartilages and on the association with forming dermal bones using the Fgfr2cC342Y/+ Crouzon syndrome mouse. We show that the dermatocranium develops outside of and in shapes that conform to the chondrocranium. Results reveal direct effects of the Fgfr2 variant on embryonic cartilage, on chondrocranium morphology, and on the association between chondrocranium and dermatocranium development. Histologically we observe a trend of relatively more chondrocytes, larger chondrocytes, and/or more matrix in the Fgfr2cC342Y/+ embryos at all timepoints before the chondrocranium begins to disintegrate at E16.5. The chondrocrania and forming dermatocrania of Fgfr2cC342Y/+ embryos are relatively large, but a contrasting trend begins at E16.5 and continues into early postnatal (P0 and P2) timepoints, with the skulls of older Fgfr2cC342Y/+ mice reduced in most dimensions compared to Fgfr2c+/+ littermates. Our findings have implications for the study and treatment of human craniofacial disease, for understanding the impact of chondrocranial morphology on skull growth, and potentially on the evolution of skull morphology. In vitro pectate lyase activity and carbon uptake assays and whole genome sequencing of Bacillus amyloliquefaciens subsp. plantarum strains for a pectin defective pathway The pectin lyase activity of 59 Bacillus amyloliquefaciens subsp. plantarum (Bap) strains was tested in vitro on Pectate Agar (PA) and Tris-Spizizen Salts (TSS) medium. Bap strains were cultured on TSA medium and washed three times with sterile water before the inoculation on PA media. Higher and lower pectate lyase activity were observed in six (AP193, AP203, AP299, AP80, AP102, and AP52) and four (AP 194, AP214, AP215, and AP305) Bap strains compared to other Bap strains. A total of 12 Bap strains (AP67, AP71, AP77, AP78, AP85, AP102, AP108, AP135, AP143, AP189, AP192, and AP193) grew vigorously on TSS medium. A total of six Bap strains (AP194, AP204, AP214, AP216, AP219, and HD73) had lower growth compared to other Bap strains. Pectin (1%) were used for in vitro PA and TSS medium. Pectate lyase and utilization activity were not found in Bacillus thuringiensis subsp. kurstaki strain HD73 compared to Bap strains. A draft genome sequence for strains AP194 and AP214 that were negative for pectin utilization were generated using an Illumina MiSeq. RAST analysis revealed that the pectin-associated gene altronate hydrolase (uxaA) absent in AP 214 strain. Multiple amino acid alignments of exuT and uxuB gene sequence showed dissimilarities among AP194, AP214, and reference Bap strains. Genomically Mined Acoustic Reporter Genes Enable In Vivo Monitoring of Tumors and Tumor-Homing Bacteria A major outstanding challenge in the fields of biological research, synthetic biology and cell-based medicine is the difficulty of visualizing the function of natural and engineered cells noninvasively inside opaque organisms. Ultrasound imaging has the potential to address this challenge as a widely available technique with a tissue penetration of several centimeters and spatial resolution below 100 {micro}m. Recently, the first genetically encoded acoustic reporters were developed based on bacterial gas vesicles to link ultrasound signals to molecular and cellular function. However, the properties of these first-generation acoustic reporter genes (ARGs) resulted in limited sensitivity and specificity for imaging in the in vivo context. Here, we describe second-generation ARGs with greatly improved acoustic properties and expression characteristics, identified through a phylogenetic screen of candidate gene clusters from diverse bacteria and archaea. The resulting constructs offer major qualitative and quantitative improvements, including much stronger ultrasound contrast, the ability to produce nonlinear signals distinguishable from background tissue, and stable long-term expression. We demonstrate the utility of these next-generation ARGs by imaging in situ gene expression in mouse models of breast cancer and tumor-homing probiotic bacteria, revealing the unique spatial distributions of tumor growth and colonization by therapeutic cells noninvasively in living subjects. Propelling and perturbing appendages together facilitate strenuous ground self-righting Terrestrial animals must self-right when overturned on the ground, but this locomotor task is strenuous. To do so, the discoid cockroach often pushes its wings against the ground to begin a somersault which rarely succeeds. As it repeatedly attempts this, the animal probabilistically rolls to the side to self- right. During winged self-righting, the animal flails its legs vigorously. Here, we studied whether wing opening and leg flailing together facilitate strenuous ground self-righting. Adding mass to increase hind leg flailing kinetic energy increased the animals self-righting probability. We then developed a robot with similar strenuous self-righting behavior and used it as a physical model for systematic experiments. The robots self-righting probability increased with wing opening and leg flailing amplitudes. A potential energy landscape model revealed that, although wing opening did not generate sufficient kinetic energy to overcome the high pitch potential energy barrier to somersault, it reduced the barrier for rolling, facilitating the small kinetic energy from leg flailing to probabilistically overcome it to self-right. The model also revealed that the stereotyped body motion during self-righting emerged from physical interaction of the body and appendages with the ground. Our work demonstrated the usefulness of potential energy landscape for modeling self-righting transitions.

Impact statementAnimals can, and robots should, use different types of appendages together to propel and perturb themselves to self-right when overturned, a strenuous yet crucial locomotor task. 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. The Mitochondrial RNA Granule Modulates Manganese-Dependent Cell Toxicity Prolonged manganese exposure causes manganism, a neurodegenerative movement disorder. The identity of adaptive and non-adaptive cellular processes targeted by manganese remains mostly unexplored. Here we study mechanisms engaged by manganese in genetic cellular models known to increase susceptibility to manganese exposure, the plasma membrane manganese efflux transporter SLC30A10 and the mitochondrial Parkinsons gene PARK2. We found that SLC30A10 and PARK2 mutations as well as manganese exposure compromised the mitochondrial RNA granule as well as mitochondrial transcript processing. These RNA granule defects led to impaired assembly and function of the mitochondrial respiratory chain. Notably, cells that survived a cytotoxic manganese challenge had impaired RNA granule function, thus suggesting that this granule phenotype was adaptive. CRISPR gene editing of subunits of the mitochondrial RNA granule, FASTKD2 or DHX30, as well as pharmacological inhibition of mitochondrial transcription-translation, were protective rather than deleterious for survival of cells acutely exposed to manganese. Similarly, adult Drosophila mutants with defects in the mitochondrial RNA granule component scully were safeguarded from manganese-induced mortality. We conclude that the downregulation of the mitochondrial RNA granule function is a protective mechanism for acute metal toxicity.

Significance StatementMutations in the manganese efflux transporter SLC30A10 and the mitochondrial Parkinsons gene PARK2, cause neurodegeneration and increased susceptibility to toxic manganese exposure. Thus, molecular processes affected in both mutants could offer insight into fundamental mechanisms conferring susceptibility or resilience to environmental and genetic factors associated with neurodegeneration. Here we report that SLC30A10 and PARK2 mutations compromise an understudied structure, the mitochondrial RNA granule, which is required for processing polycistronic mitochondrial RNAs. Cells and Drosophila lacking mitochondrial RNA granule components were resistant to manganese exposure. We conclude that the downregulation of the mitochondrial RNA granule function is an adaptive mechanism for cells exposed to manganese. The renal inflammatory network of nephronophthisis STRUCTURED ABSTRACTO_ST_ABSBACKGROUNDC_ST_ABSThe majority of genetic kidney disease leading to kidney failure is caused by mutations in ciliary genes. How cilia malfunction leads to progressive kidney damage is poorly understood, but recent evidence links ciliopathy genes to CCL2 dependent macrophage recruitment in autosomal dominant polycystic kidney disease (ADPKD), the most studied renal ciliopathy. Whether or not renal inflammation is involved in other renal ciliopathies is unclear.

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

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

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

SIGNIFICANCE STATEMENTMutations in genes encoding primary cilia proteins are the leading cause of genetic kidney failure. In autosomal dominant polycystic kidney disease (ADPKD), deregulated cilia signaling leads to kidney infiltration by macrophages through the chemokine CCL2. Little is known about renal inflammation in nephronophthisis (NPH), the most frequent pediatric renal ciliopathy. Using NPH mice models, tissues and cells from NPH patients, we unveil renal inflammation as preeminent feature of NPH. Remarkably, the renal inflammatory evoked by ciliary gene mutations in NPH does not overlap with ADPKD: it is CCL2 independent, involves a prominent recruitment of neutrophils and T cells and a specific cytokine signature. This unforeseen findings strengthen the link between primary cilia and renal inflammation. Inflammasome activation in infected macrophages drives COVID-19 pathology Severe COVID-19 is characterized by persistent lung inflammation, inflammatory cytokine production, viral RNA, and sustained interferon (IFN) response all of which are recapitulated and required for pathology in the SARS-CoV-2 infected MISTRG6-hACE2 humanized mouse model of COVID-19 with a human immune system1-20. Blocking either viral replication with Remdesivir21-23 or the downstream IFN stimulated cascade with anti-IFNAR2 in vivo in the chronic stages of disease attenuated the overactive immune-inflammatory response, especially inflammatory macrophages. Here, we show SARS-CoV-2 infection and replication in lung-resident human macrophages is a critical driver of disease. In response to infection mediated by CD16 and ACE2 receptors, human macrophages activate inflammasomes, release IL-1 and IL-18 and undergo pyroptosis thereby contributing to the hyperinflammatory state of the lungs. Inflammasome activation and its accompanying inflammatory response is necessary for lung inflammation, as inhibition of the NLRP3 inflammasome pathway reverses chronic lung pathology. Remarkably, this same blockade of inflammasome activation leads to the release of infectious virus by the infected macrophages. Thus, inflammasomes oppose host infection by SARS-CoV-2 by production of inflammatory cytokines and suicide by pyroptosis to prevent a productive viral cycle.