MTD: a unique pipeline for host and meta-transcriptome joint and integrative analyses of RNA-seq data RNA-seq data contains not only host transcriptomes but also non-host information that comprises transcripts from active microbiota in the host cells. Therefore, joint and integrative analyses of both host and meta-transcriptome can reveal gene expression of microbial community in a given sample as well as the correlative and interactive dynamics of host response to the microbiome. However, there are no convenient tools that can systemically analyze host-microbiota interactions by simultaneously quantifying host and meta-transcriptome in the same sample at the tissue and the single-cell level, which poses a challenge for interested researchers with a limited expertise in bioinformatics. Here, we developed a software pipeline that can comprehensively and synergistically analyze and correlate the host and meta-transcriptome in a single sample using bulk and single-cell RNA-seq data. This pipeline, named MTD, can extensively identify and quantify microbiome, including viruses, bacteria, protozoa, fungi, plasmids, and vectors in the host cells and correlate the microbiome with the host transcriptome. MTD is easy to install and run, involving only a few lines of simple commands. It empowers researchers with unique genomics insights into host immune responses to microorganisms.
RATTLE: Reference-free reconstruction and quantification of transcriptomes from Nanopore sequencing Nanopore sequencing enables the efficient and unbiased measurement of transcriptomes from any sample. However, current methods for transcript identification and quantification rely of mapping reads to a reference genome, which precludes the study of species with a partial or missing reference or the identification of disease-specific transcripts not readily identifiable from a reference. Here we present RATTLE, a tool to perform reference-free reconstruction and quantification of transcripts using only Nanopore reads. Using simulated data and experimental data from isoform spike-ins, human tissues, and cell lines, we show that RATTLE accurately determines transcript sequences and their abundances, and shows good scalability with the number of transcripts. RATTLE provides unprecedented access to transcriptomes from any sample and species without relying on a reference or additional technologies.
Movement planning and feedback are sufficient for sensory prediction error computation and motor adaptation Prediction errors guide many forms of learning, providing teaching signals that help us improve our performance. Implicit motor adaptation, for instance, is driven by sensory prediction errors (SPEs), which occur when the expected and observed consequences of a movement differ. Traditionally, SPE computation is thought to require movement execution. However, recent work suggesting that the brain generates and accounts for sensory predictions based on motor imagery or planning alone calls this assumption into question. Here, by measuring implicit adaptation during a visuomotor task, we tested whether motor planning and well-timed sensory feedback are sufficient for SPE computation. Human participants were cued to reach to a target and were, on a subset of trials, rapidly cued to withhold these movements. Errors displayed both on trials with and without movements induced single-trial implicit learning. Learning following trials without movements persisted even when movement trials had never been paired with errors, and when the direction of movement and sensory feedback trajectories were decoupled. These observations demonstrate that the brain can compute SPEs without generating overt movements, leading to the adaptation of planned movements even when they are not performed. Harnessing changes in open chromatin determined by ATAC-seq to generate insulin-responsive reporter constructs. BackgroundGene regulation is critical for proper cellular function. Next-generation sequencing technology has revealed the presence of regulatory networks that regulate gene expression and essential cellular functions. Studies investigating the epigenome have begun to uncover the complex mechanisms regulating transcription. Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is quickly becoming the assay of choice for many epigenomic investigations. However, whether intervention-mediated changes in accessible chromatin determined by ATAC-seq can be harnessed to generate intervention-inducible reporter constructs has not been systematically assayed.

ResultsWe used the insulin signaling pathway as a model to investigate chromatin regions and gene expression changes using ATAC- and RNA-seq in insulin-treated Drosophila S2 cells. We found correlations between ATAC- and RNA-seq data, especially when stratifying differentially-accessible chromatin regions by annotated feature type. In particular, our data demonstrated a strong correlation between chromatin regions annotated to distal promoters (1-2 kb from the transcription start site) and downstream gene expression. We cloned candidate distal promoter regions upstream of luciferase and demonstrate insulin-inducibility of several of these reporters.

ConclusionsInsulin-induced chromatin accessibility determined by ATAC-seq reveals enhancer regions that drive insulin-inducible reporter gene expression. Stimulated oxygen transport in tissue by magnetic needle acupuncture AimTo understand physiology of magnetic needle acupuncture by investigating O2 transport in tissue during needle intervention.

MethodsO2 transport in tissue is modeled by utilizing COMSOL with magnetic needle inserted into muscle tissue in a 2D porous media. The damaged tissue has been mimicked by an extracted tissue block with 1st order O2 consumption rate. The convection-diffusion O2 transport in the damaged tissue has been further evaluated by varying magnetic flux density B0 of the needle (0-1 T), myoglobin concentration (0-1 mM), O2 tension (5-100 Torr), O2 consumption rates, tissue permeability (10-12-10-6 m2) and porosity (0.1-0.9).

Results1) Active O2 transport carried by interstitial flow is enhanced with the intervention of a magnetic needle by generating a high gradient magnetic field around the tip, which exerts a strong force (104 N/m3) on the diamagnetic interstitial fluid to accelerate the flow. 2) This interstitial flow can reach 30 m/s at B0 = 1T and strongly correlates to the needle B0 and tissue permeability. 3) The needle stirs the interstitial flow can pump O2 flux by 1-2 orders of magnitude compared to that without magnetic field. 4) The enhancement of active O2 transport by magnetic needle is site-specific to the tissue in the vicinity of the tip. This enhancement is more effective in edema condition with a high tissue permeability (>10-9 m2).

ConclusionsThe dramatic enhancement of O2 transport to restore the O2 mitochondria metabolism for dysfunctional muscle tissues is the fundamental physiological mechanism of magnetic needle acupuncture. In vitro evolution of uracil glycosylase towards DnaKJ and GroEL binding evolves different misfolded states Natural evolution is driven by random mutations that improve fitness. In vitro evolution mimics this process, however, on a short time-scale and is driven by the given bait. Here, we used directed in vitro evolution of a random library of Uracil glycosylase (eUNG) displayed on yeast surface to select for binding to chaperones GroEL, DnaK+DnaJ+ATP (DnaKJ) or E.coli cell extract (CE). Using binding to the eUNG inhibitor Ugi as probe for native foldedess, the CE selected population was further divided to Ugi binders (+U) (native) or non-binders (-U). We found that GroEL, DnaKJ and CE-U select and enrich for mutations causing eUNG to misfold, with the three being enriched in mutations in buried and conserved positions, with a tendency to increase positive charge. Still, each selection, as well as CE+U and natural evolution of eUNG has its own trajectory. While GroEL and CE-U selected for mutants highly sensitive to protease cleavage, DnaKJ selected for partially structured misfolded species with a tendency to refold, making them less sensitive to proteases. CE+U selected for more neutral mutations than natural evolution. In a more general context, our results show that GroEL has a higher tendency to purge promiscuous misfolded protein mutants from the system, while DnaKJ binds mutants misfolding-prone species that are, upon chaperone release, more likely to natively refold. CE-U shares some of the properties of GroEL and DnaKJ selected populations, while harboring also unique properties, explained by the existence of additional chaperones in CE, such as Tig, HtpG and ClpB. Foraging with friends, moving with anyone: effects of the social environment on movement-integrated habitat selection O_LIMovement links the distribution of habitats with the social environment of animals using those habitats; yet integrating movement, habitat selection, and socioecology remains an opportunity for further study.
C_LIO_LIHere, our objective was to disentangle the roles of habitat selection and social association as drivers of collective movement in a gregarious ungulate. To accomplish this objective, we (1) assessed whether socially familiar individuals form discrete social communities and whether social communities have high spatial, but not necessarily temporal, overlap; and (2) we modelled the relationship between collective movement and selection of foraging habitats using socially informed integrated step selection analysis.
C_LIO_LIWe used social network analysis to assign individuals to social communities and determine short and long-term social preference among individuals. Using integrated step selection functions (iSSF), we then modelled the effect of social processes, i.e., nearest neighbour distance and social preference, and movement behaviour on patterns of habitat selection.
C_LIO_LIBased on assignment of individuals to social communities and home range overlap analyses, individuals assorted into discrete social communities, and these communities had high spatial overlap. By unifying social network analysis with iSSF, we identified movement-dependent social association, where individuals foraged with more familiar individuals, but moved collectively with any between foraging patches.
C_LIO_LIOur study demonstrates that social behaviour and space use are inter-related based on spatial overlap of social communities and movement-dependent habitat selection. Movement, habitat selection, and social behaviour are linked in theory. Here, we put these concepts into practice to demonstrate that movement is the glue connecting individual habitat selection to the social environment.
C_LI Proteasome complexes experience profound structural and functional rearrangements throughout mammalian spermatogenesis During spermatogenesis, spermatogonia undergo a series of mitotic and meiotic divisions on their path to spermatozoa. To achieve this, a succession of processes requiring high proteolytic activity are in part orchestrated by the proteasome. The spermatoproteasome (s20S) is specific to the developing gametes, in which the gamete-specific 4s subunit replaces the 4 isoform found in the constitutive proteasome (c20S). Although the s20S is conserved across species and was shown to be crucial for germ cell development, its mechanism, function and structure remain incompletely characterized. Here, we used advanced mass spectrometry (MS) methods to map the composition of proteasome complexes and their interactomes throughout spermatogenesis. We observed that the s20S becomes highly activated as germ cells enter meiosis, mainly through a particularly extensive 19S activation and, to a lesser extent, PA200 binding. Additionally, the proteasome population shifts from c20S (98%) to s20S (>82-92%) during differentiation, presumably due to the shift from 4 to 4s expression. We demonstrated that s20S, but not c20S, interacts with components of the meiotic synaptonemal complex, where it may localize via association with the PI31 adaptor protein. In vitro, s20S preferentially binds to 19S, and displays higher trypsin- and chymotrypsin-like activities, both with and without PA200 activation. Moreover, using MS methods to monitor protein dynamics, we identified significant differences in domain flexibility between 4 and 4s. We propose that these differences induced by 4s incorporation result in significant changes in the way the s20S interacts with its partners, and dictate its role in germ cell differentiation. HOXB7 induces STAT3-mediated transformation and lung metastasis in an immortalized mammary gland NMuMG cells Genes from the homeobox family are often dysregulated in cancer. Particularly, homeobox B7 (HOXB7) amplification and overexpression correlate with poor prognosis in various cancers, such as gastric, pancreatic, and lung cancers. Moreover, HOXB7 is known to contribute to cancer progression by promoting epithelial to mesenchymal transition, drug resistance, and angiogenesis. In this study, we show that HOXB7 overexpression correlates with a poor prognosis in patients with breast cancer. Coamplification of HOXB7 with ERBB2 enhanced the expression of ERBB2 and the malignancy of HER2-positive breast cancer cells. Additionally, HOXB7 overexpression in an immortalized murine mammary gland epithelial cell line, NMuMG, induced cellular transformation through an activation of the JAK-STAT signaling pathway, which resulted in lung metastasis upon orthotropic transplantation in mice. These findings reveal that HOXB7 alone is sufficient to induce tumorigenesis and lung metastasis using HOXB7-expressing NMuMG cells. 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. Ineffectual AEC1 differentiation from KRT8hi transitional state without fibrosis is associated with fatal COVID-19 ARDS ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. If Research is Not the Evidence, What is it? The Lack of Research Citations in Health Vlogs Doctor-influencers have the knowledge needed to understand scientific research, and they have the social media popularity to share it with large audiences. In this article, I explore the treatment of evidence in vlogs by Egyptian doctor-influencers: how they present it and what they use instead of it. I answer the following questions: 1) How do doctor-influencers present evidence on social media? 2) If evidence is not the research, then what is it? The data was collected through in-depth interviews with 12 Egyptian doctor-influencers from different specializations; two focus groups; and a critical discourse analysis of 48 of the influencers most popular and most engaging videos. I found that the doctor-influencers cited academic research only in the videos about medical controversies, new treatments, misinformation, or common medical mistakes. In most of their videos, the Egyptian doctor-influencers scarcely referred to medical research. This is because they believe that their audiences would not understand evidence due to their low educational attainment, weak research skills, and the language barrier (as most medical research is published in English rather than Arabic). These findings highlight the challenges of communicating science in societies with a low level of education. In these settings, communicators may need to move beyond mediating research studies on social media and instead focus on simplifying, translating, and connecting evidence-based information to the cultural context. Thirteen independent genetic loci associated with preserved processing speed in a study of cognitive resilience in 330,097 individuals in the UK Biobank Cognitive resilience is the ability to withstand the negative effects of stress on cognitive functioning and is important for maintaining quality of life while aging. UK Biobank (UKB) does not have direct measurements of the same cognitive phenotype at distal timepoints. Therefore, we used number of education years (EY) as a proxy phenotype for past cognitive performance. Current cognitive performance was determined based on processing speed. This approach captured an average time span of 40 years between past and current cognitive performance in 330,097 individuals. A confounding factor was that EY is highly polygenic and masked the genetics of resilience. To overcome this, we employed Genomics Structural Equation Modelling (GenomicSEM) to perform a GWAS-by-subtraction using two GWAS, one GWAS of EY and resilience and a second GWAS of EY but not resilience. Subtracting one from the other generated a GWAS of Resilience. Replication of this approach was shown using independent discovery and replication samples within UKB and the full UKB GWAS results show significant genetic correlation with a GWAS of cognitive change in the independent Health and Retirement Study (N=9,526; P=1.5x10-3). We found 13 independent genetic loci for Resilience. Functional analyses showed enrichment in several brain regions and involvement of specific cell types, including GABAergic neurons (P=6.59x10-8) and glutamatergic neurons (P=6.98x10-6) in the cortex. Gene-set analyses implicated the biological process "neuron differentiation" (P=9.7x10-7) and the cellular component "synaptic part" (P=2.14x10-6). Mendelian randomization analysis showed a causative effect of white matter volume on cognitive resilience. These results enhance neurobiological understanding of resilience. 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. Herbivorous fish of Sinai marine protected areas (Gulf of Aqaba): Structure baseline for potential thermal stress impact management Herbivorous reef fishes are critically important for maintaining the health and resilience of coral reefs. In the Egyptian Gulf of Aqaba (GoA), many of these herbivores are important fishery targets, however, little is known about their biomass and size structure patterns. Here, we conduct a large-scale ([~]250 km of coastline) assessment of the status of herbivorous fish populations on 30 reef sites from eight coastal regions in the GoA. These regions are subjected to three levels of fishing and protection, ranging from almost no fishing (No-take, NT), and moderately fished (via Gear-restriction, GR) to heavily fished (Open-access, OA). We found that the NT fishery reserve was the most effective in maintaining herbivorous fish richness, body-size, and biomass across all functional groups, though numerical density was not. Total herbivore biomass was 4.3 and 2.8 times higher on NT reefs and GR reefs respectively than on OA reefs. Among GoA regions, only Ras Mohammed and Sharm El-Sheikh (unfished), and Nabq (moderately fished) met the global mean herbivore biomass target of [~]30 kg/500 m2. Our work also highlights three important results regarding the impacts of fishing and fisheries management on herbivorous fish populations: (i) size structure was heavily skewed toward smaller individuals, with fishes less than (20 cm) accounting for 81.1% of the population in heavily fished OA reefs vs. 52.8% at unfished NT reefs, (ii) biomass of larger herbivore individuals (>35 cm) accounted for less than 1% of the total biomass in OA reefs vs. 37% at NT reefs, and (iii) large-bodied target species, Cetoscarus bicolor, Naso unicorn, and Kyphosus spp., accounting for 15.9% of the total biomass at NT reefs on average, while they were virtually absent from OA reefs. Collectively, these findings suggest that many principal fishery species in the central-northern regions of the GoA have been overexploited, which in turn can lead to reef degradation. Therefore, there is an urgent need to develop a participatory management approach and enforce restrictions on destructive fishing gear in order to promote recovery of herbivore biomass.

HighlightsO_LILarge-scale assessment of the status of herbivorous fishes in the northern Red Sea.
C_LIO_LIComparison of herbivorous fish assemblages across three levels of fishing pressure.
C_LIO_LILarge-bodied herbivore species were either rare or absent at heavily fished reefs.
C_LIO_LILoss of these key herbivores may boost algal dominance, thus diminish reef functionality.
C_LIO_LIPromoting local community participation in the fishery management process.
C_LI Onset and kinetics of B cell polarization depend on substrate mechanics and antigen concentration The immune synapse is the tight contact zone between a lymphocyte and a cell presenting its cognate antigen. This structure serves as a signaling platform and entails a polarization of intra-cellular components, necessary to the immunological function of the cell. While the surface properties of the presenting cell are known to control the formation of the synapse, their impact on polarization has not yet been studied.

Using functional lipid droplets as tunable artificial presenting cells combined with a microfluidic pairing device, we simultaneously observe synchronized synapses and dynamically quantify polarization patterns of individual B cells. By assessing how the ligand concentration, the surface fluidity and the substrate rigidity impact this polarization, we show that its onset depends on the antigen concentration at the synapse, and that the substrate rigidity controls both its onset and its kinetics. Our experimental system enables a fine phenotyping of monoclonal cell populations based on their synaptic readout.

Significance statementTo fulfill their immune function, B cells acquire antigens presented on the surface of antigen-presenting cells and organize their cytoskeleton and lysosomal traffic towards the contact. The polarization mechanism is evolutionary conserved, making the system interesting both from the immunology and the cell biology point of view. We developed and finely characterized antigen-coated droplets as surrogate presenting cells with a physiological-flow microfluidic cell-pairing chip to explore the B cell lysosome polarization in response to changes in droplet mechanics and antigen concentration. We observed that polarization is sensitive to mechanical properties (responding more strongly but less fast on stiffer substrates) and to the local antigen concentration at the contact. Hence, surprisingly, a late cell response is sensitive to the mechanics of the contact. XAV-19, a swine glyco-humanized polyclonal antibody against SARS-CoV-2 Spike receptor-binding domain, targets multiple epitopes and broadly neutralizes variants. Amino acid substitutions and deletions in Spike protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants can reduce the effectiveness of monoclonal antibodies (mAbs). In contrast, heterologous polyclonal antibodies raised against S protein, through the recognition of multiple target epitopes, have the potential to maintain neutralization capacities. XAV-19 is a swine glyco-humanized polyclonal neutralizing antibody raised against the receptor binding domain (RBD) of the Wuhan-Hu-1 Spike protein of SARS-CoV-2. XAV-19 target epitopes were found distributed all over the RBD and particularly cover the receptor binding motives (RBM), in direct contact sites with the Angiotensin Converting Enzyme-2 (ACE-2). Therefore, in Spike/ACE2 interaction assays, XAV-19 showed potent neutralization capacities of the original Wuhan Spike and of the United Kingdom (Alpha/B.1.1.7) and South African (Beta/B.1.351) variants. These results were confirmed by cytopathogenic assays using Vero E6 and live virus variants including the Brazil (Gamma/P.1) and the Indian (Delta/B.1.617.2) variants. In a selective pressure study with the Beta strain on Vero E6 cells conducted over 1 month, no mutation was associated with addition of increasing doses XAV-19. The potential to reduce viral load in lungs was confirmed in a human ACE2 transduced mouse model. XAV-19 is currently evaluated in patients hospitalized for COVID-19-induced moderate pneumonia in a phase 2a-2b (NCT04453384) where safety was already demonstrated and in an ongoing 2/3 trial (NCT04928430) to evaluate the efficacy and safety of XAV-19 in patients with moderate-to-severe COVID-19. Owing to its polyclonal nature and its glyco-humanization, XAV-19 may provide a novel safe and effective therapeutic tool to mitigate the severity of coronavirus disease 2019 (Covid-19) including the different variants of concern identified so far. Sperm Histone H3 Lysine 4 tri-methylation serves as a metabolic sensor of paternal obesity and is associated with the inheritance of metabolic dysfunction ObjectiveParental environmental exposures can strongly influence descendant risks for adult disease. How paternal obesity changes the sperm chromatin leading to the acquisition of metabolic disease in offspring remains controversial and ill-defined. The objective of this study was to assess: (1) whether obesity induced by a high-fat diet alters sperm histone methylation; (2) whether paternal obesity can induce metabolic disturbances across generations; (3) whether there could be cumulative damage to the sperm epigenome leading to enhanced metabolic dysfunction in descendants; and (4) whether obesity-sensitive regions associate with embryonic epigenetic and transcriptomic profiles. Using a genetic mouse model of epigenetic inheritance, we investigated the role of histone H3 lysine 4 methylation (H3K4me3) in the paternal transmission of metabolic dysfunction. This transgenic mouse overexpresses the histone demethylase enzyme KDM1A in the developing germline and has an altered sperm epigenome at the level of histone H3K4 methylation. We hypothesized that challenging transgenic sires with a high-fat diet would further erode the sperm epigenome and lead to enhanced metabolic disturbances in the next generations.

MethodsTo assess whether paternal obesity can have inter- or transgenerational impacts, and if so, to identify potential mechanisms of this non-genetic inheritance, we used wildtype C57BL/6NCrl and transgenic males with a pre-existing altered sperm epigenome. To induce obesity, sires were fed either a control or high-fat diet (10% or 60% kcal fat, respectively) for 10-12 weeks, then bred to wildtype C57BL/6NCrl female fed a regular diet. F1 and F2 descendants were characterized for metabolic phenotypes by examining the effects of paternal obesity by sex, on body weight, fat mass distribution, the liver transcriptome, intraperitoneal glucose and insulin tolerance tests. To determine whether obesity altered the F0 sperm chromatin, native chromatin immunoprecipitation-sequencing targeting H3K4me3 was performed. To gain insight into mechanisms of paternal transmission, we compared our sperm H3K4me3 profiles with embryonic and placental chromatin states, histone modification and gene expression profiles.

ResultsObesity-induced alterations in H3K4me3 occurred at genes implicated in metabolic, inflammatory, and developmental processes. These processes were associated with offspring metabolic dysfunction and corresponded to genes enriched for H3K4me3 in embryos, and overlapped embryonic and placenta gene expression profiles. Transgenerational susceptibility to metabolic disease was only observed when obese F0 had a pre-existing modified sperm epigenome. This coincided with increased H3K4me3 alterations in sperm and more severe phenotypes affecting their offspring.

ConclusionsOur data suggest sperm H3K4me3 might serve as a metabolic sensor that connects paternal diet with offspring phenotypes via the placenta. This non-DNA based knowledge of inheritance has the potential to improve our understanding of how environment shapes heritability and may lead to novel routes for the prevention of disease. This study highlights the need to further study the connection between the sperm epigenome, placental development and childrens health. 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. Discovery of highly potent pancoronavirus fusion inhibitors that also effectively inhibit COVID-19 variants from the UK (Alpha), South Africa (Beta), and India (Delta) We report here the discovery of several highly potent small molecules that showed low nM potency against SARS-CoV (IC50: as low as 13 nM), SARS-CoV-2 (IC50: as low as 23 nM), and MERS-CoV (IC50: as low as 76 nM) in pseudovirus based assays with excellent selectivity indices (SI: as high as > 5000) demonstrating their pancoronavirus inhibition. Some compounds also show 100% inhibition of CPE (IC100) at 1.25 {micro}M against an authentic SARS-CoV-2 (US_WA-1/2020). Furthermore, the most active inhibitors also potently inhibited variants of concerns (VOCs), such as the UK (B.1.1.7), South Africa (B.1.351), and Delta variant (B.1.617.2), originated in India. We confirmed that one of the potent inhibitors binds to the prefusion spike protein trimer of SARS-CoV-2 by SPR. Besides, we showed that they inhibit virus-mediated cell-cell fusion. The ADME data of one of the most active inhibitors, NBCoV1, show drug-like properties. In vivo PK of NBCoV1 in rats demonstrated excellent half-life (t1/2) of 11.3 h, mean resident time (MRT) of 14.2 h, and oral bioavailability. We expect the lead inhibitors to pave the way for further development to preclinical and clinical candidates. Sequence-dependent model of genes with dual σ factor preference Escherichia coli uses {sigma} factors to quickly control large gene cohorts during stress conditions. While most of its genes respond to a single {sigma} factor, approximately 5% of them have dual {sigma} factor preference. The most common are those responsive to both {sigma}70, which controls housekeeping genes, and {sigma}38, which activates genes during stationary growth and stresses. Using RNA-seq and flow-cytometry measurements, we show that  {sigma}70+38 genes are nearly as upregulated in stationary growth as  {sigma}38 genes. Moreover, we find a clear quantitative relationship between their promoter sequence and their response strength to changes in {sigma}38 levels. We then propose and validate a sequence dependent model of {sigma}70+38 genes, with dual sensitivity to {sigma}38 and {sigma}70, that is applicable in the exponential and stationary growth phases, as well in the transient period in between. We further propose a general model, applicable to other stresses and {sigma} factor combinations. Given this, promoters controlling {sigma}70+38 genes (and variants) could become important building blocks of synthetic circuits with predictable, sequence-dependent sensitivity to transitions between the exponential and stationary growth phases. Usutu virus escapes langerin-induced restriction to productively infect human Langerhans cells, unlike West Nile virus. Usutu virus (USUV) and West Nile virus (WNV) are phylogenetically close emerging arboviruses transmitted by mosquitoes, and constitute a global public health threat. Since USUV and WNV enter the body through the skin, the first immune cells they encounter are skin-resident dendritic cells, the most peripheral outpost of immune defense. This unique network is composed of Langerhans cells (LCs) and dermal DCs, which reside in the epidermis and the dermis, respectively.

Using human skin explants, we show that while both viruses can replicate in keratinocytes, they can also infect resident DCs with distinct tropism, since WNV preferentially infects DCs in the dermis, whereas USUV has a greater propensity to infect LCs. Using both purified human epidermal LCs (eLCs) and monocyte derived LCs (MoLCs), we confirm that LCs sustain a faster and more efficient replication of USUV compared with WNV and that this correlates with a more intense innate immune response to USUV compared with WNV. Next, we show that ectopic expression of the LC-specific C-type lectin receptor (CLR), langerin, in HEK293T cells allows WNV and USUV to bind and enter, but supports the subsequent replication of USUV only. Conversely, blocking or silencing langerin in MoLCs or eLCs made them resistant to USUV infection, thus demonstrating that USUV uses langerin to enter and replicate in LCs. Altogether, our results demonstrate that LCs constitute privileged target cells for USUV in human skin, because langerin favors its entry and replication. Intriguingly, this suggests that USUV efficiently escapes the antiviral functions of langerin, which normally safeguards LCs from most viral infections. 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 hexose transporter SWEET5 confers galactose sensitivity to Arabidopsis pollen germination via the galactokinase GALK Galactose is an abundant and essential sugar used for the biosynthesis of many macromolecules in different organisms, including plants. Galactose metabolism is tightly and finely controlled since excess galactose and derivatives are inhibitory. In Arabidopsis, root growth and pollen germination were strongly inhibited by excess galactose. However, the mechanism of galactose induced inhibition during pollen germination remains obscure. In this study, we characterized a plasma-membrane localized transporter, AtSWEET5, that transports glucose and galactose. SWEET5 protein level started to accumulate at the tricellular stage of pollen development and peaked in mature pollen before rapidly declining after pollen was germinated. SWEET5 levels are responsible for the dosage-dependent sensitivity of galactose and galactokinase (GALK) is essential for the inhibitory effects of galactose during pollen germination. Overall, SWEET5 and GALK contribute to the maintenance of galactose metabolic homeostasis during pollen germination. SWEET5 serves as a major low-affinity hexose transporter at the early stage of pollen germination.

One-sentence summarySWEET5 mediates pollen galactose sensitivity via GALK that is required for efficient galactose uptake in pollen during pollen germination. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. A combined human in silico and CRISPR/Cas9-mediated in vivo zebrafish based approach for supporting gene target validation in early drug discovery The clinical heterogeneity of heart failure has challenged our understanding of the underlying genetic mechanisms of this disease. In this respect, large-scale patient DNA sequencing studies have become an invaluable strategy for identifying potential genetic contributing factors. The complex aetiology of heart failure, however, also means that in vivo models are vital to understand the links between genetic perturbations and functional impacts. Traditional approaches (e.g. genetically-modified mice) are optimal for assessing small numbers of proposed target genes, but less practical when multiple targets are identified. The zebrafish, in contrast, offers great potential for higher throughput in vivo gene functional assessment to aid target prioritisation and support definitive studies undertaken in mice. Here we used whole-exome sequencing and bioinformatics on human patient data to identify 3 genes (API5, HSPB7, and LMO2) suggestively associated with heart failure that were also predicted to play a broader role in disease aetiology. The role of these genes in cardiovascular system development and function was then further investigated using in vivo CRISPR/Cas9-mediated gene mutation analysis in zebrafish. We observed multiple impacts in F0 knockout zebrafish embryos (crispants) following effective somatic mutation, including reductions in ventricle size, pericardial oedema, and chamber malformation. In the case of lmo2, there was also a significant impact on cardiovascular function as well as an expected reduction in erythropoiesis. The data generated from both the human in silico and zebrafish in vivo assessments undertaken supports roles for API5, HSPB7, and LMO2 in human cardiovascular disease and identifies them as potential drug targets for further investigation. The data presented also supports the use of human in silico genetic variant analysis, in combination with zebrafish crispant phenotyping, as a powerful approach for assessing gene function as part of an integrated multi-level drug target validation strategy.