The Heterochromatin protein 1 is a master regulator in RNA splicing precision deficient in ulcerative colitis Defects in RNA splicing have been linked to numerous human disorders, but remain poorly explored in inflammatory bowel disease (IBD). Here, we report that, in the gut epithelium of patients with ulcerative colitis (UC), the expression of the chromatin and alternative splicing regulator HP1{gamma} is strongly reduced. Accordingly, inactivation of the HP1{gamma} gene in the mouse gut triggered several IBD-like traits, including inflammation and dysbiosis. In parallel, we discovered that its loss of function broadly increased splicing noise, reducing requirement for canonical splicing consensus sequences, and favoring the usage of cryptic splice sites at numerous genes with key functions in gut biology. This notably resulted in the production of progerin, a noncanonical toxic splice variant of prelamin A mRNA, responsible for the Hutchinson Gilford Progeria Syndrome (HGPS) of premature aging. Likewise, production of progerin transcript was found to be a signature of colonic cells from UC patients. Thus, our study identifies HP1{gamma} as a regulator of RNA metabolism in vivo, providing a unique mechanism linking anti-inflammation and accuracy of RNA splicing in the gut epithelium. HP1 defect may confer a general disturbance in RNA splicing precision to scrutinize in IBD and more generally in accelerating aging diseases.
The OCT2 Transporter Regulates Dopamine D1 Receptor Signaling at the GolgiApparatus Dopamine is the main catecholamine in the brain and the kidney, where it controls a number of physiological functions such as locomotion, cognition, emotion, endocrine regulation and renal function. As a membrane impermeant hormone and neurotransmitter, dopamine is thought to signal by binding and activating dopamine receptors, members of the G protein couple receptor (GPCR) family, only on the plasma membrane. Here, using novel nanobody-based biosensors, we demonstrate for the first time that the dopamine D1 receptor (D1DR), the primary mediator of dopaminergic signaling in the brain and kidney, not only functions on the plasma membrane but becomes activated at the Golgi apparatus in the presence of its ligand. We present evidence that activation of the Golgi pool of D1DR is dependent on Organic Cation Transporter 2 (OCT2), a dopamine transporter, providing an explanation for how the membrane impermeant dopamine accesses subcellular pools of D1DR. We further demonstrate that dopamine activates Golgi-D1DR in the striatal medium spiny neurons (MSN) and this activity depends on OCT2 function. We also introduce a new approach to selectively interrogate compartmentalized D1DR signaling by inhibiting Gs coupling, using a nanobody-based chemical recruitment system. Using this strategy, we show that Golgi-localized D1DRs regulate cAMP production and mediate local protein kinase A activation. Together, our data suggest that spatially compartmentalized signaling hubs are previously unappreciated regulatory aspects of D1DR signaling. Our data also provide further evidence for the role of transporters in regulating subcellular GPCR activity.
Geometric machine learning informed by ground truth: Recovery of conformational continuum from single-particle cryo-EM data of biomolecules This work is based on the manifold-embedding approach to study biological molecules exhibiting continuous conformational changes. Previous work established a method capable of reconstructing 3D movies and accompanying energetics of atomic-level structures from single-particle cryo-EM images of macromolecules displaying multiple conformational degrees of freedom. Here, we introduce an unsupervised geometric machine learning approach that is informed by detailed heuristic analysis of manifolds formed by simulated heterogeneous cryo-EM datasets generated from an atomic structure. These simulated data were generated with increasing complexity to account for multiple conformational motions, state occupancies and typical microscope parameters in a wide range of signal-to-noise ratios. Using these datasets as ground-truth, we provide detailed exposition of our findings using several conformational motions while exploring the available parameter space. Guided by these insights, we build a framework to leverage the high-dimensional geometric information obtained towards reconstituting a quasi-continuum of conformational states in the form of a free-energy landscape and respective 3D density maps for all states therein. As shown by a direct comparison of results, this framework offers substantial improvements relative to the previous work. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Hepatic MIR20B promotes nonalcoholic fatty liver disease by suppressing PPARA BackgroundNon-alcoholic fatty liver disease (NAFLD) is characterized by excessive lipid accumulation and imbalances in lipid metabolism in the liver. Although nuclear receptors (NRs) play a crucial role in hepatic lipid metabolism, the underlying mechanisms of NR regulation in NAFLD remain largely unclear.

MethodsUsing network analysis and RNA-seq to determine the correlation between NRs and microRNA in human NAFLD patients, we revealed that MIR20B specifically targets PPARA. MIR20B mimic and anti-MIR20B were administered to human HepG2 and Huh-7 cells and mouse primary hepatocytes as well as high fat diet (HFD)- or methionine-deficient diet (MCD)-fed mice to verify the specific function of MIR20B in NAFLD. We tested the inhibition of the therapeutic effect of a PPAR agonist, fenofibrate, by MIR20B and the synergic effect of combination of fenofibrate with anti-MIR20B.

ResultsWe revealed that MIR20B specifically targets PPARA through miRNA regulatory network analysis of nuclear receptor genes in NAFLD. The expression of MIR20B was upregulated in free fatty acid (FA)- treated hepatocytes and the livers of both obesity-induced mice and NAFLD patients. Overexpression of MIR20B significantly increased hepatic lipid accumulation and triglyceride levels. Furthermore, MIR20B significantly reduced FA oxidation and mitochondrial biogenesis by targeting PPARA. In MIR20B-introduced mice, the effect of fenofibrate to ameliorate hepatic steatosis was significantly suppressed. Finally, inhibition of MIR20B significantly increased FA oxidation and uptake, resulting in improved insulin sensitivity and a decrease in NAFLD progression. Moreover, combination of fenofibrate and anti-MIR20B exhibited the synergic effect on improvement of NAFLD in MCD-fed mice.

ConclusionsTaken together, our results demonstrate that the novel MIR20B targets PPARA, plays a significant role in hepatic lipid metabolism, and present an opportunity for the development of novel therapeutics for NAFLD.

FundingThis research was funded by Korea Mouse Phenotyping Project (2016M3A9D5A01952411), the National Research Foundation of Korea (NRF) grant funded by the Korea government (2020R1F1A1061267, 2018R1A5A1024340, NRF-2021R1I1A2041463, 2020R1I1A1A01074940), and the Future-leading Project Research Fund (1.210034.01) of UNIST. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. Predicting gene regulatory networks from multi-omics to link genetic risk variants and neuroimmunology to Alzheimers disease phenotypes BackgroundGenome-wide association studies have found many genetic risk variants associated with Alzheimers disease (AD). However, how these risk variants affect deeper phenotypes such as disease progression and immune response remains elusive. Also, our understanding of cellular and molecular mechanisms from disease risk variants to various phenotypes is still limited. To address these problems, we performed an integrative multi-omics analysis of genotype, transcriptomics, and epigenomics for revealing gene regulatory mechanisms from disease variants to AD phenotypes.

MethodFirst, given the population gene expression data of a cohort, we construct and cluster its gene co-expression network to identify gene co-expression modules for various AD phenotypes. Next, we predict transcription factors (TFs) regulating co-expressed genes and AD risk SNPs that interrupt TF binding sites on regulatory elements. Finally, we construct a gene regulatory network (GRN) linking SNPs, interrupted TFs, and regulatory elements to target genes and gene modules for each phenotype in the cohort. This network thus provides systematic insights into gene regulatory mechanisms from risk variants to AD phenotypes.

ResultsOur analysis predicted GRNs in three major AD-relevant regions: Hippocampus, Dorsolateral Prefrontal Cortex (DLPFC), Lateral Temporal Lobe (LTL). Comparative analyses revealed cross-region-conserved and region-specific GRNs, in which many immunological genes are present. For instance, SNPs rs13404184 and rs61068452 disrupt SPI1 binding and regulation of AD gene INPP5D in the Hippocampus and LTL. However, SNP rs117863556 interrupts bindings of REST to regulate GAB2 in DLPFC only. Driven by emerging neuroinflammation in AD, we used Covid-19 as a proxy to identify possible regulatory mechanisms for neuroimmunology in AD. To this end, we looked at the GRN subnetworks relating to genes from shared AD-Covid pathways. From those subnetworks, our machine learning analysis prioritized the AD-Covid genes for predicting Covid-19 severity. Decision Curve Analysis also validated our AD-Covid genes outperform known Covid-19 genes for classifying severe Covid-19 patients. This suggests AD-Covid genes along with linked SNPs can be potential novel biomarkers for neuroimmunology in AD. Finally, our results are open-source available as a comprehensive functional genomic map for AD, providing a deeper mechanistic understanding of the interplay among multi-omics, brain regions, gene functions like neuroimmunology, and phenotypes. Heterogeneity of circulating epithelial cells in breast cancer at single-cell resolution: identifying tumor and hybrid cells Circulating tumor cells and hybrid cells formed by the fusion of tumor cells with normal cells are leading players in metastasis and have prognostic relevance. Circulating tumor cells and hybrid cells are identified as CD45-negative and CD45-positive epithelial cells. However, such an approach is challenging because epithelial cells are observed in the blood of healthy individuals. In this study, we applied single-cell RNA sequencing to profile CD45-negative and CD45-positive circulating epithelial cells (CECs) in 20 breast cancer patients and one healthy donor. DNA ploidy analysis was used to identify the tumor and hybrid cells among CD45[boxh] and CD45+ CECs in patients, respectively. Functional enrichment analysis was applied to characterize aneuploid and diploid cells. Diploid cells were also annotated to generate cell-type candidates and analyzed for copy-number aberrations (CNAs) to confirm or refute their tumor origin. CD45[boxh] and CD45+ CECs were found in cancer patients (25.5 (range 0-404) and median 6.5 (0-147)) and the healthy donor (8 and 11 cells) and divided into three clusters. Two CD45[boxh] CEC clusters were predominantly aneuploid (97% and 98%), but one cluster contained more diploid (59%) than aneuploid cells. CD45+ CECs were mostly diploid: only clusters 1 and 2 had aneuploid cells (16% and 2%). Diploid CD45[boxh] and CD45+ CECs were annotated as different immune cells and surprisingly harbored many CNAs. Cancer-associated signaling pathways were found only in aneuploid cells of CD45[boxh] CEC cluster 1 and diploid cells of CD45+ CEC cluster 1. Thus, our findings suggest that CECs in breast cancer patients are a highly heterogeneous population comprising aneuploid (tumor and hybrid) and diploid (normal) cells. DNA ploidy analysis is an effective instrument for identifying tumor and hybrid cells among CD45[boxh] and CD45+ CECs, respectively. GWAS-associated variants, non-genetic factors, and transient transcriptome in Multiple Sclerosis etiopathogenesis: a colocalization analysis A clinically actionable understanding of multiple sclerosis (MS) etiology goes through GWAS interpretation, prompting research on new gene regulatory models. We previously suggested a stochastic etiologic model where small-scale random perturbations could reach a threshold for MS development. The recently described mapping of the transient transcriptome (TT), including intergenic and intronic RNAs, seems appropriate to verify this model through a rigorous colocalization analysis. We show that genomic regions coding for the TT were significantly enriched for MS-associated GWAS variants and DNA binding sites for molecular transducers mediating putative, non-genetic, etiopathogenetic factors for MS (e.g., vitamin D deficiency, Epstein Barr virus latent infection, B cell dysfunction). These results suggest a model whereby TT-coding regions are hotspots of convergence between genetic ad non-genetic factors of risk/protection for MS, and plausibly for other complex disorders. Our colocalization analysis also provides a freely available data resource (www.mscoloc.com) for future research on MS transcriptional regulation. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. 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 compound of baicalein, wogonin, and oroxylin-A inhibiting EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. 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. SARS-CoV-2 Beta Variant Substitutions Alter Spike Glycoprotein Receptor Binding Domain Structure and Stability The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) and the subsequent COVID-19 pandemic has visited a terrible cost on the world in the forms of disease, death, and economic turmoil. The rapid development and deployment of extremely effective vaccines against SARS-CoV-2 have seemingly brought within reach the end of the pandemic. However, the virus has acquired mutations; and emerging variants of concern (VOC) are more infectious and reduce the efficacy of existing vaccines. While promising efforts to combat these variants are underway, the evolutionary pressures leading to these variants are poorly understood. To that end, here we have studied the effects on the structure and function of the SARS-CoV-2 spike glycoprotein receptor-binding domain of three amino-acid substitutions found in several variants of concern, including alpha (B.1.1.7), beta (B.1.351), and gamma (P.1). We found that these substitutions alter the RBD structure, stability, and ability to bind to ACE2, in such a way as to possibly have opposing and compensatory effects. These findings provide new insights into how these VOC may have been selected for infectivity while maintaining the structure and stability of the receptor binding domain. Adaptive Wild Introgression Mediate Temporal Wheat Root Water Influxes under Drought Stress through Plant-Atmospheric Interrelations Drought intensity as experienced by plants depends upon soil moisture status and atmospheric variables such as temperature, radiation, and air vapour pressure deficit (VPD). Although the role of shoot architecture with these edaphic and atmospheric factors is well-characterized, the extent to which shoot and root dynamic interactions as a continuum are controlled by genotypic variation is less known. Here, we targeted these interactions using a wild emmer introgression line (IL20) with a distinct drought-induced shift in the shoot-to-root ratio and its drought-sensitive recurrent parent Svevo. Using a gravimetric platform, we show that IL20 maintained higher root water influx and gas exchange under drought stress, which supported a greater growth. Interestingly, the advantage of IL20 in root water influx and transpiration was expressed earlier during the daily diurnal cycle under lower VPD and therefore supported higher transpiration efficiency. Application of structural equation model indicates that under drought, VPD and radiation are antagonistic to transpiration rate, whereas the root water influx operates as feedback for the higher atmospheric responsiveness of leaves. Collectively, our results suggest that a drought-induced shift in root-to-shoot ratio can improve plant water uptake potential in a short preferable time window determined by both water and atmospheric parameters. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. HIGHER VIRAL VIRULENCE SELECTS FOR FASTER EVOLUTION OF HOST RESISTANCE Parasites and pathogens vary strikingly in their virulence and the resulting selection they impose on their hosts. While the evolution of different virulence levels is well studied, the evolution of host resistance in response to different virulence levels is less understood and as of now mainly based on observations and theoretical predictions with few experimental tests. Increased virulence can increase selection for host resistance evolution if resistance costs are outweighed by the benefits of avoiding infection. To test this, we experimentally evolved the bacterium Vibrio alginolyticus in the presence of two variants of the filamentous phage, VALG{Phi}8, that differ in their virulence. The bacterial host exhibited two alternative defence strategies against future phage infection: (1) super infection exclusion (SIE) whereby phage-infected cells were immune to subsequent infection at a cost of reduced growth, and (2) surface receptor mutations (SRM) in genes encoding the MSHA type-IV pilus providing resistance to infection by preventing phage attachment. While SIE emerged rapidly against both phages, SRM evolved faster against the high virulence compared to the low virulence phage. Using a mathematical model of our system we show that increasing virulence strengthens selection for SRM due to the higher costs of infection suffered by SIE immune hosts. In both the experiments and the model, higher levels of SRM in the host population drove more rapid phage extinction. Thus, by accelerating the evolution of host resistance, more virulent phages caused shorter epidemics. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. Polycomb suppresses a female gene regulatory network to ensure testicular differentiation Gonadal sex determination is controlled by the support cells of testes and ovaries. In testes, the epigenetic mechanism that maintains cellular memory to suppress female sexual differentiation remains unknown. Here, we show that Polycomb suppresses a female gene regulatory network in Sertoli cells, the specific support cells for postnatal testes. Through genetic ablation, we removed Polycomb repressive complex 1 (PRC1) from embryonic Sertoli cells after sex determination. PRC1-depleted postnatal Sertoli cells exhibited defective proliferation and cell death, leading to the degeneration of adult testes. In adult Sertoli cells, PRC1 suppressed the specific, critical genes required for granulosa cells, the support cells of ovaries, thereby inactivating the female gene regulatory network. The underlying chromatin of female genes was coated with Polycomb-mediated repressive modifications: PRC1-mediated H2AK119ub and PRC2-mediated H3K27me3. Taken together, we identify a critical mechanism centered on Polycomb that maintains the male fate in adult testes. kmtricks: Efficient and flexible construction of Bloom filters for large sequencing data collections When indexing large collections of short-read sequencing data, a common operation that has now been implemented in several tools (Sequence Bloom Trees and variants, BIGSI, ..) is to construct a collection of Bloom filters, one per sample. Each Bloom filter is used to represent a set of k-mers which approximates the desired set of all the non-erroneous k-mers present in the sample. However, this approximation is imperfect, especially in the case of metagenomics data. Erroneous but abundant k-mers are wrongly included, and non-erroneous but low-abundant ones are wrongly discarded. We propose kmtricks, a novel approach for generating Bloom filters from terabase-sized collections of sequencing data. Our main contributions are 1/ an efficient method for jointly counting k-mers across multiple samples, including a streamlined Bloom filter construction by directly counting, partitioning and sorting hashes instead of k-mers, which is approximately four times faster than state-of-the-art tools; 2/ a novel technique that takes advantage of joint counting to preserve low-abundant k-mers present in several samples, improving the recovery of non-erroneous k-mers. Our experiments highlight that this technique preserves around 8x more k-mers than the usual yet crude filtering of low-abundance k-mers in a large metagenomics dataset.

Availabilityhttps://github.com/tlemane/kmtricks

FundingThe work was funded by IPL Inria Neuromarkers, ANR Inception (ANR-16-CONV-0005), ANR Prairie (ANR-19-P3IA-0001), ANR SeqDigger (ANR-19-CE45-0008). β-glucan imprinting remodels macrophage function in response to environmental cues In vitro, exposure of human primary monocytes to the fungal {beta}-glucan enhances their pro-inflammatory responsiveness towards several pathogens. Yet, the role of environmental condition of this process remains unclear. Here we found that {beta}-glucan-induced innate immune memory counteract the anti-inflammatory status of the macrophages. In response to {beta}-glucan imprinting, M-CSF-(M2-like-) macrophages increase their pro-inflammatory responsiveness to secondary stimuli associated with decrease of the M-CSF differentiation hallmarks. In contrast, in GM-CSF-(M1-like-) environment, {beta}-glucan imprinting reduced the pro-inflammatory canonical feature of the macrophages, together with their terminal differentiation marks. Comparing M-CSF and GM-CSF environment, we observed that {beta}-glucan-imprinted macrophages present comparable functions in terms of cytokine responses, phagocytosis, oxidative burst, and angiogenesis. This effect is mediated through Dectin-1 and associated with altered expression of the master regulators of macrophage terminal differentiation, IRF5 and IRF3. {beta}-glucan-induced innate immune memory skews the commitment of the macrophages in complex environment towards similar and less terminally differentiated cells. Together, these observations suggest a potential therapeutic role for {beta}-glucan-induced modulation of innate memory in different pathological contexts. Accuracy and Power Analysis of Social Networks Built From Count Data O_LIPower analysis is used to estimate the probability of correctly rejecting a null hypothesis for a given statistical model and dataset. Conventional power analyses assume complete information, but the stochastic nature of behavioural sampling can mean that true and estimated networks are poorly correlated. Power analyses do not currently take the effect of sampling into account. This could lead to inaccurate estimates of statistical power, potentially yielding misleading results.
C_LIO_LIHere we develop a method for computing network correlation: the correlation between an estimated social network and its true network, using a Gamma-Poisson model of social event rates for networks constructed from count data. We use simulations to assess how the level of network correlation affects the power of nodal regression analyses. We also develop a generic method of power analysis applicable to any statistical test, based on the concept of diminishing returns.
C_LIO_LIWe demonstrate that our network correlation estimator is both accurate and moderately robust to its assumptions being broken. We show that social differentiation, mean social event rate, and the harmonic mean of sampling times positively impacts the strength of network correlation. We also show that the required level of network correlation to achieve a given power level depends on many factors, but that 0.80 network correlation usually corresponds to around 0.80 power for nodal regression in ideal circumstances.
C_LIO_LIWe provide guidelines for using our network correlation estimator to verify the accuracy of networks built from count data, and to conduct power analysis. This can be used prior to data collection, in post hoc analyses, or even for subsetting networks in dynamic network analysis. The network correlation estimator and custom power analysis methods have been made available as an R package.
C_LI Beta-cell specific Insr deletion promotes insulin hypersecretion and improves glucose tolerance prior to global insulin resistance Abstract

Insulin receptor (Insr) protein can be found at higher levels in pancreatic {beta}-cells than in most other tissues, but the consequences of {beta}-cell insulin resistance remain enigmatic. Ins1cre allele was used to delete Insr specifically in {beta}-cells of both female and male mice. Experimental mice were compared to Ins1cre-containing littermate controls at multiple ages and on multiple diets. RNA-seq of purified recombined {beta}-cells revealed transcriptomic consequences of Insr loss, which differed between female and male mice. Action potential and calcium oscillation frequencies were increased in Insr knockout {beta}- cells from female, but not male mice, whereas only male {beta}InsrKO mice had reduced ATP-coupled oxygen consumption rate and reduced expression of genes involved in ATP synthesis. Female {beta}InsrKO and {beta}InsrHET mice exhibited elevated insulin release in perifusion experiments, during hyperglycemic clamps, and following i.p. glucose challenge. Deletion of Insr did not alter {beta}-cell area up to 9 months of age, nor did it impair hyperglycemia-induced proliferation. Based on our data, we adapted a mathematical model to include {beta}-cell insulin resistance, which predicted that {beta}-cell Insr knockout would improve glucose tolerance depending on the degree of whole-body insulin resistance. Indeed, glucose tolerance was significantly improved in female {beta}InsrKO and {beta}InsrHET mice when compared to controls at 9, 21 and 39 weeks, and also in insulin-sensitive 4-week old males. We did not observe improved glucose tolerance in older male mice or in high fat diet-fed mice, corroborating the prediction that global insulin resistance obscures the effects of {beta}-cell specific insulin resistance. The propensity for hyperinsulinemia was associated with mildly reduced fasting glucose and increased body weight. We further validated our main in vivo findings using the Ins1-CreERT transgenic line and found that male mice had improved glucose tolerance 4 weeks after tamoxifen-mediated Insr deletion. Collectively, our data show that loss of {beta}-cell Insr contributes to glucose-induced hyperinsulinemia, thereby improving glucose homeostasis in otherwise insulin sensitive dietary and age contexts. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. The Inflammatory Role of Milk Fat Globule Epidermal Growth Factor VIII in Angiotensin II Induced Arterial Remodeling BackgroundAngiotensin II (Ang II) and milk fat globule-epidermal growth factor VIII (MFG-E8) are involved in age-associated arterial remodeling; however, the inflammatory role of MFG-E8 in Ang II associated arterial remodeling with aging remains unknown.

Methods and ResultsIn this study, 30-week-old MFG-E8 knock out (KO) and age-matched wild-type (WT) mice were infused with Ang II or saline. After infusion the with Ang II, the aortic molecular, cellular, and structural remodeling were observed in mice and compared to those infused with saline, but these effects were dependent on the expression of MFG-E8: (1) In the WT mice, Ang II infusion substantially increased intimal-medial thickness, elastic lamina degradation, collagen deposition, and the proliferation of VSMCs; in contrast, in the KO mice, these effects were significantly reduced; (2) In the WT mice, Ang II treatment significantly increased the activation and expression of MMP2, TGF-{beta}1, and its downstream signaling molecule p-SMAD2, and collagen type I production, however, in the KO mice, these molecular effects were significantly reduced; (3) In the WT mice, Ang II treatment increased inflammatory p-NF-{kappa}B p65, MCP1, TNF-, ICAM1, and VCAM1 molecular expression, while conversely, in the KO mice, no significant inflammatory changes were found; (4) Importantly, compared to untreated control mice with a wide range of age from 4-96 weeks, Ang II infused "younger" mice produced an "older" arterial inflammatory phenotype, which was alleviated by MFG-E8 deficiency.

ConclusionsMFG-E8 mediates Ang II associated arterial inflammatory remodeling. Targeting MFG-E8 is a novel molecular approach to curb adverse arterial remodeling during aging and hypertension.

CLINICAL PERSPECTIVEO_ST_ABSWhat Is New?C_ST_ABS* Both Ang II and MFG-E8 increases are involved in proinflammatory arterial remodeling mediating the molecular, cellular and tissue events in aging and hypertension.
* MFG-E8 is essential for Ang II induced and age-associated adverse arterial remodeling via the increase of proinflammation, intimal medial thickening, elastin fragmentation, collagen deposition, and VSMC proliferation.


What Are the Clinical Implications?Since MFG-E8 mediates Ang II induced proinflammation in arterial wall remodeling in aging and hypertension, targeting MFG-E8 is a potential molecular approach to curb inflammatory arterial remodeling, maintaining the health of the vascular system during aging and hypertension. Dramatic consequences of reducing erythrocyte membrane cholesterol on Plasmodium falciparum Cholesterol is the most abundant lipid in the erythrocyte. During its blood stage development, the malaria parasite establishes an active cholesterol gradient across the various membrane systems within the infected erythrocyte. Interestingly, some antimalarial compounds have recently been shown to disrupt cholesterol homeostasis in intraerythrocytic stages of Plasmodium falciparum. These studies point to the importance of cholesterol for parasite growth. Previously, reduction of cholesterol from the erythrocyte membrane by treatment with methyl-{beta}-cyclodextrin (M{beta}CD) was shown to inhibit parasite invasion and growth. In addition, M{beta}CD treatment of trophozoite stage P. falciparum was shown to result in parasite expulsion from the host cell. We have revisited these phenomena by using live video microscopy, ultrastructural analysis, and response to antimalarial compounds. By using time-lapse video microscopy of fluorescently tagged parasites, we show that M{beta}CD treatment for just 30 min causes dramatic expulsion of the trophozoite stage parasites. This forceful expulsion occurs within 10 sec. Remarkably, the plasma membrane of the host cell from which the parasite has been expelled does not appear to be compromised. The parasitophorous vacuolar membrane (PVM) continued to surround the extruded parasite, but the PVM appeared damaged. Treatment with antimalarial compounds targeting PfATP4 or PfNCR1 prevented M{beta}CD-mediated extrusion of the parasites, pointing to a potential role of cholesterol dynamics underlying the expulsion phenomena. We also confirmed the essential role of erythrocyte plasma membrane cholesterol for invasion and growth of P. falciparum. This defect can be partially complemented by cholesterol and desmosterol but not with epicholesterol, revealing stereospecificity underlying cholesterol function. Overall, our studies advance previous observations and reveal unusual cell biological features underlying cholesterol depletion of infected erythrocyte plasma membrane.

ImportanceMalaria remains a major challenge in much of the world. Symptoms of malaria are caused by the growth of parasites belonging to Plasmodium spp. inside the red blood cells (RBC), leading to their destruction. The parasite depends upon its host for much of its nutritional need. Cholesterol is a major lipid in the RBC plasma membrane, which is the only source of this lipid for malarial parasites. We have previously shown that certain new antimalarial compounds disrupt cholesterol homeostasis in P. falciparum. Here we use live time-lapse video microscopy to show dramatic expulsion of the parasite from the host RBC when the cholesterol content of the RBC is reduced. Remarkably, this expulsion is inhibited by the antimalarials that disrupt lipid homeostasis. We also show stereospecificity of cholesterol in supporting parasite growth inside RBC. Overall, these results point to a critical role of cholesterol in physiology of malaria parasites.