sgcocaller and comapr: personalised haplotype assembly and comparative crossover map analysis using single-gamete sequencing data Profiling gametes of an individual enables the construction of personalised haplotypes and meiotic crossover landscapes, now achievable at larger scale than ever through the availability of high-throughput single-cell sequencing technologies. However, high-throughput single-gamete data commonly have low depth of coverage per gamete, which challenges existing gametebased haplotype phasing methods. In addition, haplotyping a large number of single gametes from high-throughput singlecell DNA sequencing data and constructing meiotic crossover profiles using existing methods requires intensive processing. Here, we introduce efficient software tools for the essential tasks of generating personalised haplotypes and calling crossovers in gametes from single-gamete DNA sequencing data (sgcocaller), and constructing, visualising, and comparing individualised crossover landscapes from single gametes (comapr). With additional data pre-possessing, the tools can also be applied to bulk-sequenced samples. We demonstrate that sgcocaller is able to generate impeccable phasing results for high-coverage datasets, on which it is more accurate and stable than existing methods, and also performs well on low-coverage single-gamete sequencing datasets for which current methods fail. Our tools achieve highly accurate results with user-friendly installation, comprehensive documentation, efficient computation times and minimal memory usage.
Accommodating site variation in neuroimaging data using normative and hierarchical Bayesian models AO_SCPLOWBSTRACTC_SCPLOWThe potential of normative modeling to make individualized predictions from neuroimaging data has enabled inferences that go beyond the case-control approach. However, site effects are often confounded with variables of interest in a complex manner and can bias estimates of normative models, which has impeded the application of normative models to large multi-site neuroimaging data sets. In this study, we suggest accommodating for these site effects by including them as random effects in a hierarchical Bayesian model. We compared the performance of a linear and a non-linear hierarchical Bayesian model in modeling the effect of age on cortical thickness. We used data of 570 healthy individuals from the ABIDE (autism brain imaging data exchange) data set in our experiments. In addition, we used data from individuals with autism to test whether our models are able to retain clinically useful information while removing site effects. We compared the proposed single stage hierarchical Bayesian method to several harmonization techniques commonly used to deal with additive and multiplicative site effects using a two stage regression, including regressing out site and harmonizing for site with ComBat, both with and without explicitly preserving variance related to age and sex as biological variation of interest. In addition, we made predictions from raw data, in which site has not been accommodated for. The proposed hierarchical Bayesian method showed the best predictive performance according to multiple metrics. Beyond that, the resulting z-scores showed little to no residual site effects, yet still retained clinically useful information. In contrast, performance was particularly poor for the regression model and the ComBat model in which age and sex were not explicitly modeled. In all two stage harmonization models, predictions were poorly scaled, suffering from a loss of more than 90 % of the original variance. Our results show the value of hierarchical Bayesian regression methods for accommodating site variation in neuroimaging data, which provides an alternative to harmonization techniques. While the approach we propose may have broad utility, our approach is particularly well suited to normative modelling where the primary interest is in accurate modelling of inter-subject variation and statistical quantification of deviations from a reference model.

1 HighlightsO_LIDevelopment and presentation of normative modeling approach based on hierarchical Bayesian modeling that can be applied to large multi-site neuroimaging data sets.
C_LIO_LIComparison of performance of Hierarchical Bayesian model including site as predictor to several common ways to harmonize for multi-site effects.
C_LIO_LIPresentation of normative modeling as site correction tool.
C_LI
Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. Once-daily feeding is associated with better health in companion dogs: Results from the Dog Aging Project A variety of diets have been studied for possible anti-aging effects. In particular, studies of intermittent fasting and time-restricted feeding in laboratory rodents have found evidence of beneficial health outcomes. Companion dogs represent a unique opportunity to study diet in a large mammal that shares human environments. The Dog Aging Project has been collecting data on thousands of companion dogs of all different ages, sizes, and breeds since 2019. We leveraged this diverse cross-sectional dataset to investigate associations between feeding frequency and cognitive function (n = 10,474) as well as nine broad categories of health conditions (n = 24,238). Controlling for sex, age, breed, and other potential confounders, we found that dogs fed once daily rather than more frequently had lower mean scores on a cognitive dysfunction scale, and lower odds of having gastrointestinal, dental, orthopedic, kidney/urinary, and liver/pancreas disorders. Therefore, we find that once-daily feeding is associated with better health in multiple domains. Future research with longitudinal data can provide stronger evidence for a possible causal effect of feeding frequency on health in companion dogs. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. Spatio-temporal, optogenetic control of gene expression in organoids Organoids derived from stem cells become increasingly important to study human development and to model disease. However, methods are needed to control and study spatio-temporal patterns of gene expression in organoids. To this aim, we combined optogenetics and gene perturbation technologies to activate or knock-down RNA of target genes, at single-cell resolution and in programmable spatio-temporal patterns. To illustrate the usefulness of our approach, we locally activated Sonic Hedgehog (SHH) signaling in an organoid model for human neurodevelopment. High-resolution spatial transcriptomic and single-cell analyses showed that this local induction was sufficient to generate stereotypically patterned organoids in three dimensions and revealed new insights into SHHs contribution to gene regulation in neurodevelopment.

With this study, we propose optogenetic perturbations in combination with spatial transcriptomics as a powerful technology to reprogram and study cell fates and tissue patterning in organoids. Embryonic exposure to cannabidiol disrupts active neural circuits, an effect increased by Δ9-tetrahydrocannabinol and involving CB1R and CB2R in zebrafish In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Human cortical areas are sparsely connected: Combining histology with diffusion MRI to estimate the absolute number of axons The tracts between cortical areas are conceived as playing a central role in cortical information processing, but their actual numbers have never been determined in humans. Here we estimate the absolute number of axons linking cortical areas from a whole-cortex diffusion-MRI (dMRI) connectome, calibrated using the histologically-measured callosal fiber density. Median connectivity is estimated as ~6200 axons between cortical areas within-hemisphere and ~1300 axons inter-hemispherically, with axons connecting functionally-related areas surprisingly sparse. For example, we estimate that <5% of the axons in the trunk of the arcuate and superior longitudinal fasciculi connect Wernickes and Brocas areas. These results suggest that detailed information is transmitted between cortical areas either via linkage of the dense local connections or via rare, extraordinarily privileged long-range connections. Manipulating the microbiome alters regenerative outcomes in Xenopus laevis tadpoles via lipopolysaccharide signalling Xenopus laevis tadpoles can regenerate functional tails, containing spinal cord, notochord, muscle, fin, blood vessels and nerves, except for a brief refractory period at around one week of age. At this stage, amputation of the tadpoles tail may either result in scarless wound healing, or the activation of a regeneration programme, which replaces the lost tissues. We recently demonstrated a link between bacterial lipopolysaccharides and successful tail regeneration in refractory stage tadpoles, and proposed that this could result from lipopolysaccharides binding to Toll-like receptor 4 (TLR4). Here, we have used 16S rRNA sequencing to show that the tadpole skin microbiome is highly variable between sibships and that the community can be altered by raising embryos in the antibiotic gentamicin. Six gram-negative genera, including Delftia and Chryseobacterium, were over-represented in tadpoles that underwent tail regeneration. Lipopolysaccharides purified from a commensal Chryseobacterium spp. XDS4, an exogenous Delftia spp. or Escherichia coli could significantly increase the number of antibiotic-raised tadpoles that attempted regeneration. Conversely, the quality of regeneration was impaired in native-raised tadpoles exposed to the antagonistic lipopolysaccharide of Rhodobacter sphaeroides. Editing TLR4 using CRISPR/Cas9 also reduced regeneration quality, but not quantity, at the level of the cohort. However, we found that the editing level of individual tadpoles was a poor predictor of regenerative outcome. In conclusion, our results suggest that variable regeneration in refractory stage tadpoles depends at least in part on the skin microbiome and lipopolysaccharide signalling, but that signalling via TLR4 cannot account for all of this effect. Variation in white spruce needle respiration at the species range limits: a potential impediment to northern expansion White spruce (Picea glauca) spans a massive range from arctic treeline to temperate forests, yet the variability in respiratory physiology and related implications for tree carbon balance at the extremes of this distribution remain as enigmas. Working at both the most northern and southern extents of the white spruce distribution range more than 5000 km apart, we measured the short- term temperature response of dark respiration (R/T) at upper and lower canopy positions. R/T curves were fit to both polynomial and thermodynamic models so that model parameters could be compared among locations, canopy positions, and with previously published data. Respiration measured at 25{degrees}C (R25) was 68% lower at the southern location than at the northern location (0.73{+/-}0.15 vs. 2.27{+/-}0.02 mol m-2 s-1), resulting in a significantly lower (p< 0.01) intercept in R/T response in temperate trees. Only at the southern location did upper canopy leaves have a steeper temperature response than lower canopy leaves, likely reflecting steeper canopy gradients in light. No differences were observed in the maximum temperature of respiration. At the northern range limit respiration is nearly twice that of the average R25 reported in a global leaf respiration database. This large carbon cost likely challenges tree survival and contributes to restricting the location of the northern treeline. We predict that without significant thermal acclimation, foliage respiration will increase with projected end-of-the-century warming and will likely constrain the future range limits of this important boreal species.

Summary StatementWhite spruce (Picea glauca) needle respiration at the northern limit of the species range is three times higher than at the southern range limit (when measured at 25 {degrees}C). This high carbon cost likely challenges tree survival and contributes to the location of the northern treeline. Evaluation of US state pollinator plans using 3 evidence-based policymaking frameworks After the US federal government created a national pollinator protection plan in 2015, many states followed with their own. Since their goal is to promote pollinating insect conservation, we wanted to know whether the state plans are using best practices for evidence-based science policy. In early 2019 we found and downloaded every existing, publicly available US state pollinator protection plan. We then used content analysis to assess the goals, scope, and implementation of state-level pollinator protection plans across the US. This analysis was conducted using three distinct frameworks for evidence-based policymaking: US Department of Interior Adaptive Resources Management (ARM), US Environmental Protection Agency management pollinator protection plan (MP3) guidance, and Pew Trusts Pew-MacAthur Results First Project elements of evidence-based state policymaking (PEW) framework. Then we scored them using the framework criteria, to assess whether the plans were using known best practices for evidence based policymaking. Of the 31 states with a state pollinator plan, Connecticut was the state with the lowest total score across the three evaluation frameworks. The state with the highest overall scores, across the three frameworks, was Missouri. Most states did not score highly on the majority of the frameworks. Overall, many state plans were lacking policy elements that address monitoring, evaluation, and adjustment. These missing elements impact states ability to achieve their conservation goals. Our results indicate that states can improve their pollinator conservation policies to better match evidence-based science policy guidance, regardless of which framework is used. Cargo-mediated mechanisms reduce inter-motor mechanical interference, promote load-sharing and enhance processivity in teams of molecular motors In cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of cytoskeletal filaments. How motors that typically mechanically interfere with each other, work as teams is unclear. Here we explored the possibility that purely physical mechanisms may potentially enhance teamwork, both at the single motor and cargo level. To explore these mechanisms, we developed a 3D dynamical simulation of cargo transport along microtubules by teams of canonically non-cooperative kinesin-1 motors. We accounted for cargo membrane fluidity by explicitly simulating the Brownian dynamics of motors on the cargo surface and considered both the load and ATP dependence of single motor functioning. We showed explicitly, for the first time, that surface fluidity leads to the reduction of negative mechanical interference between kinesins, enhancing load sharing thereby decreasing single motor off-rates and increasing processivity. Remarkably, we also showed that, independent of fluidity, increasing numbers of bound motors pull the cargo closer to the microtubule, increasing the on-rates of individual unbound motors, resulting in a cooperative increase in bound motor numbers that depends on 3D cargo geometry. At the cargo level, surface fluidity makes more motors available for binding, though this effect is significant only at low ATP or high motor density. Interestingly, we find that the fluidity induced reduction in mechanical interference dominates over the increased availability of motors for typical physiological and in vitro conditions, this allowing us to reconcile different experimental results in different regimes. Finally, we show that these effects can altogether result in enhanced mechanical efficiency and tunable cargo run-lengths for teams of molecular motors over physiological ranges of fluidity with implications for new experimental validation efforts, transport in vivo, artificial cargo design and multi-motor-driven mechanical processes in general.

Author summaryIn cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of protein filaments. How do intrinsically non-cooperative motors function as teams within cells? In this paper, we show, using computer simulations, that the fluid surfaces of cellular cargo reduce the mechanical interference between motors allowing better load sharing and decreasing their detachment rates, thereby increasing the distance over which they can carry cargo. We additionally found that motors pull cargo closer to the filament increasing the attachment rates of other unbound motors. These effects act synergistically with an increased availability of motors due to fluidity to further increase travel distances. Our simulation model accounted for the dynamics of the motors on the cargo surface and their load and ATP dependent kinetics, allowing us to connect single motor properties to overall transport. Our work on understanding how teamwork arises in mechanically coupled motors sheds new light on cellular processes, reconciles existing observations, encourages new experimental validation efforts and can also suggest new ways of improving mechanical efficiency and transport of artificial cargo powered by motor teams. Evolutionary rates of testes-expressed genes differ in monogamous and promiscuous Peromyscus species Reproductive proteins, including those expressed in the testes, are among the fastest evolving proteins across the tree of life. Sexual selection on traits involved in sperm competition is thought to be a primary driver of testes gene evolution and is expected to differ between promiscuous and monogamous species due to intense competition between males to fertilize females in promiscuous lineages and lack thereof in monogamous ones. Here, we employ the rodent genus Peromyscus as a model to explore differences in evolutionary rates between testes-expressed genes in monogamous and promiscuous species. We find candidate genes that may be associated with increased sperm production in promiscuous species and gene ontology categories that show patterns of molecular convergence associated with phenotypic convergence in independently evolved monogamous species. Overall, our results highlight the possible molecular consequences of differences in mating system, likely due to differences in selective pressures. Improved AlphaFold modeling with implicit experimental information Machine learning prediction algorithms such as AlphaFold1 and RoseTTAFold2 can create remarkably accurate protein models, but these models usually have some regions that are predicted with low confidence or poor accuracy3-6. We hypothesized that by implicitly including experimental information, a greater portion of a model could be predicted accurately, and that this might synergistically improve parts of the model that were not fully addressed by either machine learning or experiment alone. An iterative procedure was developed in which AlphaFold models are automatically rebuilt based on experimental density maps and the rebuilt models are used as templates in new AlphaFold predictions. We find that including experimental information improves prediction beyond the improvement obtained with simple rebuilding guided by the experimental data. This procedure for AlphaFold modeling with density has been incorporated into an automated procedure for crystallographic and electron cryo-microscopy map interpretation. Robust and Remote Measurement of Heart Rate Based on A Surveillance Camera Arrhythmia is a marked symptom of many cardiovascular diseases (CVDs). The accurate and timely detection of heart rate can greatly reduce the harm of arrhythmia to people. However, it is still a challenge to robustly and remotely measure heart rate in daily life due to the changing environmental conditions during measurement, such as the varying light intensity, the movement of people, and the uncertain distance between the sensor to people. In this study, we propose a method to accurately measure human heart rate within a distance of 4.5 meters under different light intensities by simply using a surveillance camera. After a 20-second color video of a persons hand is captured by the camera, a method based on the Fast Fourier Transform (FFT) algorithm is designed to extract the blood volume pulse wave to calculate the heart rate. According to the comparison between the real heart rate and results measured by electrocardiography (ECG), the proposed method achieves an accuracy of 95.8% when the measurement is performed within a distance of 4.5 meters and 90% when within 5.0 meters. Our experiments show that when the illuminance varies between 100-1000 lux (lighting level indoor), we still get the correct results. Our experiments also demonstrate that the proposed method accurately obtain heart rate even when the light intensity is below 32 lux (300-500 lux in a workplace environment). The methods strong adaptability to changing environmental conditions makes it applicable to many scenarios, such as homes of the elderly, classrooms, and other public spaces. SARS-CoV-2 B.1.617.2 Delta variant emergence, replication and immune evasion The SARS-CoV-2 B.1.617.2 (Delta) variant was first identified in the state of Maharashtra in late 2020 and spread throughout India, outcompeting pre-existing lineages including B.1.617.1 (Kappa) and B.1.1.7 (Alpha). In vitro, B.1.617.2 is 6-fold less sensitive to serum neutralising antibodies from recovered individuals, and 8-fold less sensitive to vaccine-elicited antibodies as compared to wild type Wuhan-1 bearing D614G. Serum neutralising titres against B.1.617.2 were lower in ChAdOx-1 versus BNT162b2 vaccinees. B.1.617.2 spike pseudotyped viruses exhibited compromised sensitivity to monoclonal antibodies against the receptor binding domain (RBD) and N-terminal domain (NTD), in particular to the clinically approved bamlavinimab and imdevimab monoclonal antibodies. B.1.617.2 demonstrated higher replication efficiency in both airway organoid and human airway epithelial systems as compared to B.1.1.7, associated with B.1.617.2 spike being in a predominantly cleaved state compared to B.1.1.7. Additionally we observed that B.1.617.2 had higher replication and spike mediated entry as compared to B.1.617.1, potentially explaining B.1.617.2 dominance. In an analysis of over 130 SARS-CoV-2 infected healthcare workers across three centres in India during a period of mixed lineage circulation, we observed substantially reduced ChAdOx-1 vaccine efficacy against B.1.617.2 relative to non-B.1.617.2. Compromised vaccine efficacy against the highly fit and immune evasive B.1.617.2 Delta variant warrants continued infection control measures in the post-vaccination era. Synergy of chemotherapy and macrophage depletion leads to B cell-mediated T cell memory activation and durable triple negative breast cancer regression Immunosuppressive elements within the tumor microenvironment such as Tumor Associated Macrophages (TAMs) can present a barrier to successful anti-tumor responses by cytolytic T cells. We employed preclinical syngeneic p53 null mouse models of triple negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was used to successfully treat several highly aggressive TNBC murine mammary tumors and lung metastasis. Using this regimen and single cell RNA sequencing we characterized tumor infiltrating lymphocytes (TILs) including helper T cells and antigen-presenting B cells that were highly enriched in good responders to combination therapy. Using high dimensional imaging techniques, we identified the close spatial localization of B220+ CD86+ activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks post-treatment in one model that also exhibited long-term tumor regression post-treatment. We also characterized the transcriptional and metabolic heterogeneity of TAMS in these two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model is highly expressed and conserved in human claudin-low breast cancers, and high expression of the T12 signature correlated with reduced overall survival. This T12 tumor TAM signature may help identify human claudin-low breast cancer patients that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rationale combinations of immunotherapy.

SignificanceA treatment regimen that harnessed the immunostimulatory effects of cyclophosphamide coupled with the inhibition of CSF1R was used to successfully treat several highly aggressive claudin-low TNBC murine mammary tumors and lung metastasis. Bias correction for inverse variance weighting Mendelian randomization BackgroundInverse-variance weighted two-sample Mendelian randomization (IVW-MR) is the most widely used approach that utilizes genome-wide association studies (GWAS) summary statistics to infer the existence and the strength of the causal effect between an exposure and an outcome. Estimates from this approach can be subject to different biases due to the use of weak instruments and winners curse, which can change as a function of the overlap between the exposure and outcome samples.

MethodsWe developed a method (MRlap) that simultaneously considers weak instrument bias and winners curse, while accounting for potential sample overlap. Assuming spike-and-slab genomic architecture and leveraging LD-score regression and other techniques, we could analytically derive, reliably estimate, and hence correct for the bias of IVW-MR using association summary statistics only.

ResultsWe tested our approach using simulated data for a wide range of realistic settings. In all the explored scenarios, our correction reduced the bias, in some situations by as much as 30 folds. When applied to real data on obesity-related exposures, we observed significant differences between IVW-based and corrected effects, both for non-overlapping and fully overlapping samples. Additionally, our results are consistent with the fact that the strength of the biases will decrease as the sample size increases. Finally, we also showed that the overall bias is also dependent on the genetic architecture of the exposure, and traits with low heritability and/or high polygenicity are more strongly affected.

ConclusionsOur method not only reduces bias in causal effect estimation but also enables the use of much larger GWAS sample sizes, by allowing for potentially overlapping samples. Gain-of-function cardiomyopathic mutations in RBM20 rewire splicing regulation and re-distribute ribonucleoprotein granules within processing bodies RNA binding motif protein 20 (RBM20) is a key regulator of alternative splicing in the heart, and its mutation leads to malignant dilated cardiomyopathy (DCM). To understand the mechanism of RBM20-associated DCM, we engineered isogenic human induced pluripotent stem cells (iPSCs) with heterozygous or homozygous DCM-associated missense mutations in RBM20 (R636S) as well as RBM20 knockout (KO) iPSCs. iPSC-derived engineered heart tissues made from these cell lines recapitulated contractile dysfunction of RBM20-associated DCM and revealed greater dysfunction with missense mutations than KO. Analysis of RBM20 RNA binding by eCLIP revealed a gain-of-function preference of mutant RBM20 for 3' UTR sequences that are shared with amyotrophic lateral sclerosis (ALS) and processing-body associated RNA binding proteins (FUS, DDX6). Deep RNA sequencing revealed that the RBM20 R636S mutant has unique gene, splicing, polyadenylation and circular RNA defects that differ from RBM20 KO, impacting distinct cardiac signaling pathways. Splicing defects specific to KO or R636S mutations were supported by data from R636S gene-edited pig hearts and eCLIP. Super-resolution microscopy verified that mutant RBM20 maintains limited nuclear localization potential; rather, the mutant protein associates with cytoplasmic processing bodies (DDX6) under basal conditions, and with stress granules (G3BP1) following acute stress. Taken together, our results highlight a novel pathogenic mechanism in cardiac disease through splicing-dependent and -independent pathways that are likely to mediate differential contractile phenotypes and stress-associated heart pathology. EMBED: Essential Microbiome Dynamics, a dimensionality reduction approach for longitudinal microbiome studies 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. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Assessment of Sex Bias in Housekeeping Gene Expression in Adipose Tissue Through the Massive Analysis of Transcriptomics Data As the housekeeping genes (HKGs) generally involved in maintaining essential cell functions are typically assumed to exhibit constant expression levels across cell types, they are commonly employed as internal controls in gene expression studies. Nevertheless, multiple studies indicate that not all HKGs display stable expression across cells and tissues and under various healthy and diseased conditions, which can introduce systematic errors into experimental results. The selection and validation of HKGs as controls for each studied condition represent crucial steps in ensuring the validity of obtained results; however, up till now, sex has not been typically considered as a biological variable.

In this study, we evaluate the expression profiles of six classical HKGs (four metabolic: GAPDH, HPRT, PPIA, and UBC, and two ribosomal: 18S and RPL19) to determine expression stability in adipose tissues of Homo sapiens and Mus musculus and check sex bias and their overall suitability as internal controls. We also assess the stabile expression of genes included in distinct whole-transcriptome microarrays available from the Gene Expression Omnibus database (GEO) to identify sex-unbiased HKGs suitable for use as internal controls. We perform a sex-based analysis to describe any sexual dimorphisms in mRNA expression stability.

A novel computational strategy based on meta-analysis techniques proves that certain classical HKGs fail to function adequately as controls when analyzing human adipose tissue (HAT) considering sex as a variable. The extensively used 18S gene displays sex-based variability in adipose tissue, although PPIA and RPL19 do not, and hence, represent robust HKGs. We propose new sex-unbiased human and mouse HKGs (suHKG) derived from sex-specific expression profiles, such as RPS8 and UBB. All results generated during this study are readily available by accessing an open web resource (https://bioinfo.cipf.es/metafun-HKG) for consultation and reuse in further studies. N-VEGF, the forgotten isoform of VEGF-A, activate cellular stress survival circuits Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. Heterologous vaccination regimens with self-amplifying RNA and Adenoviral COVID vaccines induce robust immune responses in mice Several vaccines have demonstrated efficacy against SARS-CoV-2 mediated disease, yet there is limited data on the immune response induced by heterologous vaccination regimens using alternate vaccine modalities. Here, we present a detailed description of the immune response, in mice, following vaccination with a self-amplifying RNA (saRNA) vaccine and an adenoviral vectored vaccine (ChAdOx1 nCoV-19/AZD1222) against SARS-CoV-2. We demonstrate that antibody responses are higher in two dose heterologous vaccination regimens than single dose regimens. Neutralising titres after heterologous prime-boost were at least comparable or higher than the titres measured after homologous prime boost vaccination with viral vectors. Importantly, the cellular immune response after a heterologous regimen is dominated by cytotoxic T cells and Th1+ CD4 T cells which is superior to the response induced in homologous vaccination regimens in mice. These results underpin the need for clinical trials to investigate the immunogenicity of heterologous regimens with alternate vaccine technologies. 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.