Sex determination gene transformer regulates the male-female difference in Drosophila fat storage via the adipokinetic hormone pathway Sex differences in whole-body fat storage exist in many species. For example, Drosophila females store more fat than males. Yet, the mechanisms underlying this sex difference in fat storage remain incompletely understood. Here, we identify a key role for sex determination gene transformer (tra) in regulating the male-female difference in fat storage. Normally, a functional Tra protein is present only in females, where it promotes female sexual development. We show that loss of Tra in females reduced whole-body fat storage, whereas gain of Tra in males augmented fat storage. Tras role in promoting fat storage was largely due to its function in neurons, specifically the Adipokinetic hormone (Akh)-producing cells (APCs). Our analysis of Akh pathway regulation revealed a male bias in APC activity and Akh pathway function, where this sex-biased regulation influenced the sex difference in fat storage by limiting triglyceride accumulation in males. Importantly, Tra loss in females increased Akh pathway activity, and genetically manipulating the Akh pathway rescued Tra-dependent effects on fat storage. This identifies sex-specific regulation of Akh as one mechanism underlying the male-female difference in whole-body triglyceride levels, and provides important insight into the conserved mechanisms underlying sexual dimorphism in whole-body fat storage.
The Skeletal Muscle Circadian Clock Regulates Titin Splicing Through RBM20 Circadian rhythms are maintained by a cell autonomous, transcriptional-translational feedback loop known as the molecular clock. While previous research suggests a role of the molecular clock in regulating skeletal muscle structure and function, no mechanisms have connected the molecular clock to sarcomere filaments. Utilizing inducible, skeletal muscle specific, Bmal1 knockout (iMSBmal1-/-) mice, we showed that knocking out skeletal muscle clock function alters titin isoform expression through both RNAseq and LC-MS. This alteration in titins spring length resulted in sarcomere length heterogeneity. We demonstrate the direct link between altered titin splicing and sarcomere length in vitro using U7 snRNPs that truncate the region of titin altered in iMSBmal1-/- muscle. Finally, we identified a mechanism whereby the skeletal muscle clock regulates titin isoform expression through transcriptional regulation of Rbm20, a potent splicing regulator of titin. Our findings demonstrate the importance of the skeletal muscle circadian clock in maintaining titin isoform through regulation of RBM20 expression. Because circadian rhythm disruption is a feature of many chronic diseases, our results highlight a novel pathway that could be targeted to maintain skeletal muscle structure and function in a range of pathologies.
Mining divergent and conserved signals in 5' splicing site sequences across fungi, metazoa and plants The main steps of the splicing process are similar across eukaryotes. However differences in splicing factors, gene architecture and sequence divergences suggest clade-specific features of splicing and its regulation. In each organism the ensemble of 5 splicing sequences reflects the balance between natural nucleotidic variability and minimal molecular constraints to assure splicing fidelity. This compromise shapes the underlying statistical patterns in donor sequences composition.

In this work we aimed to mine conserved and divergent signals in splicing donor sequences. As 5 donor sequences are a major cue for proper recognition of splicing sites we reasoned that statistical regularities of their sequence composition might reflect biological functionality and evolutionary history associated to splicing mechanisms.

We considered a regularized maximum entropy modeling framework to mine for non-trivial two-site correlations in donor sequences of 30 different eukaryote organisms. Our approach allowed us to accommodate and extend within a unified framework many of the regularities observed in previous works, like the negative epistatic effects between exonic and intronic consensus sites. In addition, for each analyzed organism, we could identify minimal sets of two-site coupling patterns that could generate, at a given regularization level, observed one-site and two-site frequencies in donor sequences. Noticeably, performing a systematic and comparative analysis of 5ss we showed that lineage information could be traced from joint di-nucleotide probabilities. Specifically, we could identify characteristic two-site coupling patterns for plants and for animals and argue that they could echo differences in splicing regulation previously reported between these groups.

Author summaryThe sequence composition of 5 splicing sites of eukaryote organisms reflect a complex scenario. Nucleotide variability has to coexist with the need to correctly define exon/intron boundaries and the fidelity of splicing seems to depend on a pattern of trade-offs between substitutions at different positions of the splicing site. Better understanding of these patterns may help to gain insight into the details underlying the splicing process and its evolution.

In this work we aimed to study conserved and divergent signatures embedded in the sequence composition of eukaryotic 5 splicing sites. We developed generative probabilistic models that allowed us to analyze sequence composition of donor sequences for several eukaryote organisms. Our regularized models served to incrementally disentangle the minimal set of coupling parameters needed to accurately reproduced observed 1-site and 2-site nucleotide frequencies. Focusing our study on di-nucelotide probabilities we found that they actually carried phylogenetic signal. In particular, our comparative analysis allowed us to identify differential two-site coupling patterns for animal and plants that might be related to specific differences in splicing regulation. Structure-guided synthesis of FK506 and FK520 analogs with increased selectivity exhibit in vivo therapeutic efficacy against Cryptococcus Calcineurin is an essential virulence factor that is conserved across human fungal pathogens including Cryptococcus neoformans, Aspergillus fumigatus, and Candida albicans. Although an excellent target for antifungal drug development, the serine-threonine phosphatase activity of calcineurin is conserved in mammals and inhibition of this activity results in immunosuppression. FK506 (tacrolimus) is a naturally produced macrocyclic compound that inhibits calcineurin by binding to the immunophilin FKBP12. Previously, our fungal calcineurin-FK506-FKBP12 structure-based approaches identified a non-conserved region of FKBP12 that can be exploited for fungal-specific targeting. These studies led to the design of an FK506 analog, APX879, modified at the C22 position that was less immunosuppressive yet maintained antifungal activity. We now report high resolution protein crystal structures of fungal FKBP12 and a human truncated calcineurin-FKBP12 bound to a natural FK506 analog, FK520 (ascomycin). Based on information from these structures and the success of APX879, we synthesized and screened a novel panel of C22-modified compounds derived from both FK506 and FK520. One compound, JH-FK-05 demonstrates broad-spectrum antifungal activity in vitro and is non-immunosuppressive in vivo. In murine models of pulmonary and disseminated C. neoformans infection, JH-FK-05 treatment significantly reduced fungal burden and extended animal survival alone and in combination with fluconazole. Furthermore, molecular dynamic simulations performed with JH-FK-05 binding to fungal and human FKBP12 identified additional residues outside of the C22 and C21 positions that could be modified to generate novel FK506 analogs with improved antifungal activity.

SignificanceDue to rising rates of antifungal drug resistance and a limited armamentarium of antifungal treatments, there is a paramount need for novel antifungal drugs to treat systemic fungal infections. Calcineurin has been established as an essential and conserved virulence factor in several fungi, making it an attractive antifungal target. However, due to the immunosuppressive action of calcineurin inhibitors, they have not been successfully utilized clinically for antifungal treatment in humans. Recent availability of crystal structures of fungal calcineurin bound inhibitor complexes have enabled the structure-guided design of FK506 analogs and led to a breakthrough in the development of a compound increased for fungal specificity. The development of a calcineurin inhibitor with reduced immunosuppressive activity and therapeutic antifungal activity would add a significant tool to the treatment options for these invasive fungal infections with exceedingly high rates of mortality. Chromatin-contact atlas reveals disorder-mediated protein interactions and moonlighting chromatin-associated RBPs RNA-binding proteins (RBPs) play diverse roles in regulating co-transcriptional RNA-processing and chromatin functions, but our knowledge of the repertoire of chromatin-associated RBPs (caRBPs) and their interactions with chromatin remains limited. Here, we developed SPACE (Silica Particle Assisted Chromatin Enrichment) to isolate global and regional chromatin components with high specificity and sensitivity, and SPACEmap to identify the chromatin-contact regions in proteins. Applied to mouse embryonic stem cells, SPACE identified 1,459 chromatin-associated proteins, [~]48% of which are annotated as RBPs, indicating their dual roles in chromatin and RNA-binding. Additionally, SPACEmap stringently verified chromatin-binding of 404 RBPs and identified their chromatin-contact regions. Notably, SPACEmap showed that about half of the caRBPs bind to chromatin by intrinsically disordered regions (IDRs). Studying SPACE and total proteome dynamics from mES cells grown in 2iL and serum medium indicates significant correlation (R = 0.62). One of the most dynamic caRBPs is Dazl, which we find co-localized with PRC2 at transcription start sites of genes that are distinct from Dazl mRNA binding. Dazl and other PRC2-colocalised caRBPs are rich in intrinsically disordered regions (IDRs), which could contribute to the formation and regulation of phase-separated PRC condensates. Together, our approach provides an unprecedented insight into IDR-mediated interactions and caRBPs with moonlighting functions in native chromatin. Water levels primarily drive variation in photosynthesis and nutrient use of scrub Red Mangroves in the southeastern Florida Everglades We investigated how mangrove-island micro-elevation (i.e., habitat: center vs. edge) affects tree physiology in a scrub mangrove forest of the southeastern Everglades. We measured leaf gas exchange rates of scrub Rhizophora mangle L. trees monthly during 2019, hypothesizing that CO2 assimilation (Anet) and stomatal conductance (gsw) would decline with increasing water levels and salinity, expecting more-considerable differences at mangrove-island edges than centers, where physiological stress is greatest. Water levels varied between 0 and 60 cm from the soil surface, rising during the wet season (May-October) relative to the dry season (November-April). Porewater salinity ranged from 15 to 30 ppt, being higher at mangrove-island edges than centers. Anet maximized at 15.1 {micro}mol m-2 s-1, and gsw was typically <0.2 mol m-2 s-1, both of which were greater in the dry than the wet season and greater at island centers than edges, with seasonal variability being roughly equal to variation between habitats. After accounting for season and habitat, water level positively affected Anet in both seasons but did not affect gsw. Our findings suggest that inundation stress (i.e., water level) is the primary driver of variation in leaf gas exchange rates of scrub mangroves in the Florida Everglades, while also constraining Anet more than gsw. The interaction between inundation stress due to permanent flooding and habitat varies with season as physiological stress is alleviated at higher-elevation mangrove-island center habitats during the dry season. Freshwater inflows during the wet season, increase water levels and inundation stress at higher-elevation mangrove-island centers, but also potentially alleviate salt and sulfide stress in soils. Thus, habitat heterogeneity leads to differences in nutrient and water acquisition and use between trees growing in island centers versus edges, creating distinct physiological controls on photosynthesis, which likely affect carbon flux dynamics of scrub mangroves in the Everglades. Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Where there's smoke, there's fuel: dynamic vegetation data improve predictions of wildfire hazard in the Great Basin Wildfires are a growing management concern in western US rangelands, where invasive annual grasses have altered fire regimes and contributed to an increased incidence of catastrophic large wildfires. Fire activity in arid, non-forested ecosystems is thought to be largely controlled by interannual variation in fuel amount, which in turn is controlled by antecedent weather. Thus, long-range forecasting of fire activity in rangelands should be feasible given annual estimates of fuel quantity. Using a 32 yr time series of spatial data, we employed machine learning algorithms to predict the relative probability of large (>405 ha) wildfire in the Great Basin based on fine-scale annual and 16-day estimates of cover and production of vegetation functional groups, weather, and multitemporal scale drought indices. We evaluated the predictive utility of these models with a leave-one-year-out cross-validation, building spatial hindcasts of fire probability for each year that we compared against actual footprints of large wildfires. Herbaceous aboveground biomass production, bare ground cover, and long-term drought indices were the most important predictors of burning. Across 32 fire seasons, 88% of the area burned in large wildfires coincided with the upper 3 deciles of predicted fire probabilities. At the scale of the Great Basin, several metrics of fire activity were moderately to strongly correlated with average fire probability, including total area burned in large wildfires, number of large wildfires, and maximum fire size. Our findings show that recent years of exceptional fire activity in the Great Basin were predictable based on antecedent weather-driven growth of fine fuels and reveal a significant increasing trend in fire probability over the last three decades driven by widespread changes in fine fuel characteristics. Allopatric divergence of cooperators limits cheating range and effects of a defector mutation Social defectors may meet diverse cooperators. Genotype-by-genotype interactions may constrain the ranges of cooperators upon which particular defectors can cheat, limiting cheater spread. The bacterium Myxococcus xanthus undergoes cooperative multicellular development, but some developmental defectors cheat on cooperators during this process. In this study, interactions between a cheater disrupted at the signaling gene csgA and allopatrically diversified cooperators reveal a very small cheating range. Expectedly, the cheater failed to cheat on all natural-isolate cooperators owing to non-cheater-specific antagonisms. Surprisingly, lab-evolved cooperators that diverged from their cheating-susceptible ancestor by fewer than 20 mutations and without experiencing cheating had already exited the csgA mutants cheating range. Cooperators might also diversify in the potential for a mutation to reduce expression of cooperative trait or generate a cheating phenotype. A new csgA mutation constructed in several highly diverged cooperators generated diverse sporulation phenotypes, ranging from a complete defect to no defect, indicating that genetic backgrounds can limit the set of genomes in which a mutation creates a defector. Our results suggest that natural populations feature geographic mosaics of cooperators that have diversified in their susceptibility to particular cheaters and in the phenotypes generated by any given cooperation-gene mutation.

Significance statementSelection on cooperators exploited by obligate cheaters can induce evolution of resistance to cheating. Here we show that cooperators can also rapidly evolve immunity to cheating simply as a byproduct of evolutionary divergence in environments in which cooperation and cheating at the focal trait do not occur because the trait is not expressed. We also find that differences in the genomic context in which a cooperation-gene mutation arises can profoundly alter its phenotypic effect and determine whether the mutation generates a social defect at all - a pre-requisite for obligate cheating. These findings suggest that general divergence of social populations under a broad range of environmental conditions can restrict both the set of mutations that might generate social defectors in the first place and the host range of such defectors once they arise. Outlier detection in multimodal MRI identifies rare individual phenotypes among 20,000 brains Outliers in neuroimaging represent spurious data or the data of unusual phenotypes that deserve special attention such as clinical follow-up. Outliers have usually been detected in a supervised or semi-supervised manner for labeled neuroimaging cohorts. There has been much less work using unsupervised outlier detection on large unlabeled cohorts like the UK Biobank brain imaging dataset. Given its large sample size, rare imaging phenotypes within this unique cohort are of interest, as they are often clinically relevant and could be informative for discovering new processes. Here we developed a two-level outlier detection and screening methodology to characterize individual outliers of multiple different brain imaging phenotypes from 20,000 UK Biobank subjects. In primary screening, every subject was parameterized with an outlier score per imaging phenotype to quantitate the degree of outlierness. This approach enabled the assessment of test-retest reliability via outlier scores, which ranged from excellent reliability for ventricular volume, white matter lesion volume, and fractional anisotropy, to good reliability for mean diffusivity and cortical thickness. Resting-state functional connectivity was eliminated for individual-level outlier screening due to its low test-retest reliability of outlier scores. In secondary screening, the extreme outliers (1110 subjects) were examined individually, and those arising from data collection/processing errors were eliminated. A subgroup (108 subjects) of the remaining non-artifactual outliers were radiologically reviewed, and radiological findings were identified in 98%. This study establishes an unsupervised framework for investigating rare individual imaging phenotypes within a large neuroimaging cohort. Rescue of stalled clathrin-mediated endocytosis by asymmetric Arp2/3-mediated actin assembly Actin assembly facilitates vesicle formation in several trafficking pathways. Clathrin-mediated endocytosis (CME) shows elevated actin assembly dependence under high membrane tension. Why actin assembly at CME sites occurs heterogeneously even within the same cell, and how assembly forces are harnessed, are not fully understood. Here, endocytic dynamics, actin presence, and geometry of CME proteins from three different functional modules, were analyzed using three-dimensional (3D) super-resolution microscopy, live-cell imaging, and machine-learning-based computation. When hundreds of CME events were compared, sites with actin assembly showed a distinct signature, a delay between completion of coat expansion and vesicle scission, indicating that actin assembly occurs preferentially at stalled CME sites. N-WASP is recruited to one side of CME sites where it is positioned to stimulate asymmetric actin assembly. We propose that asymmetric actin assembly rescues stalled CME sites by pulling vesicles into the cell much like a bottle opener pulls off a bottle cap. 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. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. EMBED: a low dimensional reconstruction of subject-specific gut microbiome dynamics based on ecological normal modes Dimensionality reduction can offer unique insights into high dimensional microbiome dynamics by leveraging collective abundance fluctuations of multiple bacteria driven by similar ecological perturbations. However, methods providing lower-dimensional representations of microbiome dynamics both at the community and individual taxa level are not currently available. To that end, we present EMBED: Essential MicroBiomE Dynamics, a probabilistic non-linear tensor factorization approach. Similar to normal mode analysis in structural biophysics, EMBED infers ecological normal modes (ECNs), which represent the unique orthogonal modes capturing the collective behavior of microbial communities. A very small number of ECNs can accurately approximate microbiome dynamics across multiple data sets. Inferred ECNs reflect specific ecological behaviors, providing natural templates along which the dynamics of individual bacteria may be partitioned. Moreover, the multi-subject treatment in EMBED systematically identifies subject-specific and universal abundance dynamics that are not detected by traditional approaches. Collectively, these results highlight the utility of EMBED as a versatile dimensionality reduction tool for studies of microbiome dynamics. SARS-CoV-2 Can Infect Human Embryos The spread of SARS-CoV-2 has led to a devastating pandemic, with infections resulting in a range of symptoms collectively known as COVID-19. The full repertoire of human tissues and organs susceptible to infection is an area of active investigation, and some studies have implicated the reproductive system. The effects of COVID-19 on human reproduction remain poorly understood, and particularly the impact on early embryogenesis and establishment of a pregnancy are not known. In this work, we explore the susceptibility of early human embryos to SARS-CoV-2 infection. We note that ACE2 and TMPRSS2, two canonical cell entry factors for SARS-CoV-2, are co-expressed in cells of the trophectoderm in blastocyst-stage preimplantation embryos. Using fluorescent reporter virions pseudotyped with Spike (S) glycoprotein from SARS-CoV-2, we observe robust infection of trophectoderm cells, and this permissiveness could be attenuated with blocking antibodies targeting S or ACE2. When exposing human blastocysts to the live, fully infectious SARS-CoV-2, we detected cases of infection that compromised embryo health. Therefore, we identify a new human target tissue for SARS-CoV-2 with potential medical implications for reproductive health during the COVID-19 pandemic and its aftermath. Transmembrane helices are an overlooked and evolutionarily conserved source of major histocompatibility complex class I and II epitopes Cytolytic T cell responses are predicted to be biased towards membrane proteins. The peptide-binding grooves of most haplotypes of histocompatibility complex class I (MHC-I) are relatively hydrophobic, therefore peptide fragments derived from human transmembrane helices (TMHs) are predicted to be presented more often as would be expected based on their abundance in the proteome. However, the physiological reason of why membrane proteins might be over-presented is unclear. In this study, we show that the over-presentation of TMH-derived peptides is general, as it is predicted for bacteria and viruses and for both MHCI and MHC-II. Moreover, we show that TMHs are evolutionarily more conserved, because single nucleotide polymorphisms (SNPs) are present relatively less frequently in TMH-coding chromosomal regions compared to regions coding for extracellular and cytoplasmic protein regions. Thus, our findings suggest that both cytolytic and helper T cells respond more to membrane proteins, because these are evolutionary more conserved. We speculate that TMHs therefore are less prone to escape mutations that enable pathogens to evade T cell responses. Spatial organization of nuclear pores in Xenopus laevis oocytes Nuclear pores are protein assemblies inserted in the nuclear envelope of eukaryotic cells, acting as main gates for communication between nucleus and cytoplasm. So far, nuclear pores have been extensively studied to determine their structure and composition, yet their spatial organization and geometric arrangement on the nuclear surface are still poorly understood. Here, we analyze super-resolution images of the surface of Xenopus laevis oocyte nuclei during development, and characterize the arrangement of nuclear pores using tools commonly employed to study the atomic structural and topological features of soft matter. To interpret the experimental results, we hypothesize an effective interaction among nuclear pores and implemented it in extensive numerical simulations of octagonal clusters mimicking typical pore shapes. Thanks to our simple model, we find simulated spatial distributions of nuclear pores that are in excellent agreement with experiments, suggesting that an effective interaction among nuclear pores exists and could explain their geometrical arrangement. Furthermore, our results show that the statistical features of the geometric arrangement of nuclear pores do not depend on the type of pore-pore interaction, attractive or repulsive, but are mainly determined by the octagonal symmetry of each single pore. These results pave the way to further studies needed to determine the biological nature of pore-pore interactions. CRCM5484: A BET- BDII Selective Compound With Differential Anti-Leukemic Drug Modulation Differentially screening the Fr-PPIChem chemical library on the BET BRD4-BDII versus -BDI bromodomains led to the discovery of a BDII selective tetrahydropyridothienopyrimidinone (THPTP)-based compound. Structure-activity relationship (SAR) and hit-to-lead approaches allowed us to develop CRCM5484, a potent inhibitor of BET proteins with a preferential and 475-fold selectivity for the second bromodomain of the BRD3 protein (BRD3-BDII) over its first bromodomain (BRD3-BDI). Its very low activity was demonstrated in various cell-based assays, corresponding with recent data describing other selective BDII compounds. However, screening on a drug sensitivity and resistance-profiling platform revealed its ability to modulate the anti-leukemic activity in combination with various FDA-approved and/or in-development drugs in a cell- and context-dependent differential manner. Altogether, the results confirm the originality of the THPTP molecular mode of action in the BD cavity and its potential as starting scaffold for the development of potent and selective bromodomain inhibitors. The genome of the cereal pest Sitophilus oryzae: a transposable element haven BackgroundThe rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

ResultsWe sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 Kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino-acids required for insect development and cuticle biosynthesis.

ConclusionsHere we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes. Key features of the genetic architecture and evolution of host-microbe interactions revealed by high-resolution genetic mapping of the mucosa-associated gut microbiome in hybrid mice Determining the forces that shape diversity in host-associated bacterial communities is critical to understanding the evolution and maintenance of metaorganisms. To gain deeper understanding of the role of host genetics in shaping gut microbial traits, we employed a powerful genetic mapping approach using inbred lines derived from the hybrid zone of two incipient house mouse species. Further, we uniquely performed our analysis on microbial traits measured at the gut mucosal interface, which is in more direct contact with host cells and the immune system. A high number of mucosa-associated bacterial taxa have significant heritability estimates; heritabilities are greater for 16S rRNA transcript-compared to gene copy-based traits, and interestingly, are positively correlated with cospeciation rate estimates. Genomewide association mapping identifies 443 loci influencing 123 taxa, with narrow genomic intervals pinpointing promising candidate genes and pathways. Importantly, we identified an enrichment of candidate genes associated with several human diseases, including inflammatory bowel disease, and functional categories including innate immunity and G-protein-coupled receptors. These results highlight key features of the genetic architecture of mammalian host-microbe interactions and how they diverge as new species form. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. 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 Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. NUDT18 catalyzes the hydrolysis of active metabolites of the antivirals Remdesivir, Ribavirin and Molnupiravir Remdesivir and Molnupiravir have gained considerable interest due to their activity against SARS-CoV-2. Cellular hydrolysis of their active triphosphate forms, Remdesivir-TP and Molnupiravir-TP, would decrease drug efficiency. We therefore tested Remdesivir-TP as a substrate against a panel of human hydrolases and found that NUDT18 catalyzes the hydrolysis of Remdesivir-TP. The kcat value of NUDT18 for Remdesivir-TP was determined to 2.6 s-1 and the Km value was 156 M, suggesting that NUDT18 catalyzed hydrolysis of Remdesivir-TP occurs in cells. We demonstrate that the triphosphates of the antivirals Ribavirin and Molnupiravir are hydrolyzed by NUDT18, albeit with a lower efficiency compared to Remdesivir-TP. NUDT18 also hydrolyses the triphosphates of Sofosbuvir and Aciclovir although with significantly lower activity. These results suggest that NUDT18 can act as a cellular sanitizer of modified nucleotides and may influence the antiviral efficacy of Remdesivir, Molnupiravir and Ribavirin. NUDT18 is expressed in respiratory epithelial cells and may limit the antiviral efficacy of Remdesivir and Molnupiravir against SARS-CoV2 replication by decreasing the intracellular concentration of their active metabolites at their intended site of action.