Parameterized syncmer schemes improve long-read mapping MotivationSequencing long reads presents novel challenges to mapping. One such challenge is low sequence similarity between the reads and the reference, due to high sequencing error and mutation rates. This occurs, e.g., in a cancer tumor, or due to differences between strains of viruses or bacteria. A key idea in mapping algorithms is to sketch sequences with their minimizers. Recently, syncmers were introduced as an alternative sketching method that is more robust to mutations and sequencing errors.

ResultsWe introduce parameterized syncmer schemes, a generalization of syncmers, and provide a theoretical analysis for multi-parameter schemes. By combining these schemes with downsampling or minimizers we can achieve any desired compression and window guarantee. We implemented the use of parameterized syncmer schemes in the popular minimap2 and Winnowmap2 mappers. In tests on simulated and real long read data from a variety of genomes, the syncmer-based algorithms, with scheme parameters selected on the basis of the theoretical analysis, reduced unmapped reads by 20-60% at high compression while usually using less memory. The advantage was more pronounced at low sequence identity. At sequence identity of 75% and medium compression, syncmer-minimap had only 37% as many unmapped reads, and 8% fewer of the reads that did map were incorrectly mapped. Even at lower compression and error rates, parameterized syncmer based mapping mapped more reads than the original minimizer-based mappers as well as mappers using the original syncmer schemes. We conclude that using parameterized syncmer schemes can improve mapping of long reads in a wide range of settings.

Availabilityhttps://github.com/Shamir-Lab/syncmer_mapping

Supplementary informationSupplementary data are available at https://github.com/Shamir-Lab/syncmer_mapping.

Author summaryPopular long read mappers use minimizers, the minimal hashed k-mers from overlapping windows, as alignment seeds. Recent work showed that syncmers, which select a fixed set of k-mers as seeds, are more likely to be conserved under errors or mutations than minimizers, making them potentially useful for mapping error-prone long reads. We introduce a framework for creating syncmers, that we call parameterized syncmer schemes, which generalize those introduced so far, and provide a theoretical analysis of their properties. We implemented parameterized syncmer schemes in the minimap2 and Winnowmap2 long read mappers. Using parameters selected on the basis of our theoretical analysis we demonstrate improved mapping performance, with fewer unmapped and incorrectly mapped reads on a variety of simulated and real datasets. The improvements are consistent across a broad range of compression rates and sequence identities, with the most significant improvements for lower sequence identity (high error or mutation rates) and high compression.
MEDICC2: whole-genome doubling aware copy-number phylogenies for cancer evolution Chromosomal instability (CIN) and somatic copy-number alterations (SCNA) play a key role in the evolutionary process that shapes cancer genomes. SC-NAs comprise many classes of clinically relevant events, such as localised amplifications, gains, losses, loss-of-heterozygosity (LOH) events, and recently discovered parallel evolutionary events revealed by multi-sample phasing. These events frequently appear jointly with whole genome doubling (WGD), a transformative event in tumour evolution involving tetraploidization of genomes preceded or followed by individual chromosomal copy-number changes and associated with an overall increase in structural CIN.

While SCNAs have been leveraged for phylogeny reconstruction in the past, existing methods do not take WGD events into account and cannot model parallel evolution. They frequently make use of the infinite sites assumption, do not model horizontal dependencies between adjacent genomic loci and can not infer ancestral genomes. Here we present MEDICC2, a new phylogeny inference algorithm for allele-specific SCNA data that addresses these shortcomings. MEDICC2 dispenses with the infinite sites assumption, models parallel evolution and accurately identifies clonal and subclonal WGD events. It times SCNAs relative to each other, quantifies SCNA burden in single-sample studies and infers phylogenetic trees and ancestral genomes in multi-sample or single-cell sequencing scenarios with thousands of cells.

We demonstrate MEDICC2s ability on simulated data, real-world data of 2,778 single sample tumours from the Pan-cancer analysis of whole genomes (PCAWG), 10 bulk multi-region prostate cancer patients and two recent single-cell datasets of triple-negative breast cancer comprising several thousands of single cells.
A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. Engineered bacteria detect tumor DNA Advances in bacterial engineering have catalysed the development of living cell diagnostics and therapeutics1-3, including microbes that respond to diseases such as gut inflammation4, intestinal bleeding5, pathogens6 and hypoxic tumors7. Bacteria can easily access the entire gastrointestinal tract via oral administration8, and they can produce outputs that can be noninvasively measured in stool4 or urine7. Cellular memory, such as bistable switches4,9,10 or genomic rearrangement11, has been used to allow bacteria to store information over time. However, living biosensors have not yet been engineered to detect specific DNA sequences or mutations from outside the cell. Here, we engineer naturally competent Acinetobacter baylyi to detect donor DNA from the genomes of colorectal cancer (CRC) cells and organoids. We characterize the functionality of the biosensors in vitro with co-culture assays and then validate in vivo with sensor bacteria delivered orally or rectally into mice injected with orthotopic donor CRC organoids. We observe horizontal gene transfer from the tumor to the sensor bacteria in vivo, allowing their detection in stool. The sensor bacteria achieved 100% discrimination between mice with and without CRC using both delivery methods. Our findings establish a framework for biosensing applications that require the detection of mutations or organisms within environments that are difficult to sample. In addition, the platform can be readily expanded to include in situ production and delivery of therapeutic payloads at the detection site. Synthetic antigen presenting cells reveal the diversity and functional specialisation of extracellular vesicles composing the fourth signal of T cell immunological synapses. The T cell Immunological Synapse (IS) is a pivotal hub for the regulation of adaptive immunity by endowing the exchange of information between cells engaged in physical contacts. Beyond the integration of antigen (signal one), co-stimulation (signal two), and cytokines (signal three), the IS facilitates the delivery of T-cell effector assemblies including supramolecular attack particles (SMAPs) and extracellular vesicles (EVs). How these particulate outputs differ among T -cell subsets and how subcellular compartments and signals exchanged at the synapse contribute to their composition is not fully understood. Here we harnessed bead-supported lipid bilayers (BSLBs) as a tailorable and versatile technology for the study of synaptic particle biogenesis and composition in different T-cell subsets, including CART. These synthetic antigen-presenting cells (APCs) facilitated the characterisation of trans-synaptic vesicles (tSV) as a heterogeneous population of EVs comprising among others PM-derived synaptic ectosomes and CD63+ exosomes. We harnessed BSLB to unveil the factors influencing the vesicular release of CD40L, as a model effector, identifying CD40 trans presentation, T-cell activation, ESCRT upregulation/recruitment, antigen density/potency, co-repression by PD-1 ligands, and its processing by ADAM10 as major determinants. Further, BSLB made possible the comparison of microRNA (miR) species associated with tSV and steadily released EVs. Altogether, our data provide evidence for a higher specialisation of tSV which are enriched not only in effector immune receptors but also in miR and RNA-binding proteins. Considering the molecular uniqueness and functional complexity of the tSV output, which is also accompanied by SMAPs, we propose their classification as signal four.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=198 HEIGHT=200 SRC="FIGDIR/small/445691v3_ufig1.gif" ALT="Figure 1">
View larger version (67K):
org.highwire.dtl.DTLVardef@100f3aorg.highwire.dtl.DTLVardef@57adfforg.highwire.dtl.DTLVardef@6052fcorg.highwire.dtl.DTLVardef@1e85c21_HPS_FORMAT_FIGEXP  M_FIG C_FIG HighlightsO_LIBead Supported Lipid Bilayers (BSLB) reconstituting antigen-presenting cells support synapse assembly by T cells and the release of effector particles.
C_LIO_LIBSLB facilitate the dissection of the cellular machineries and synapse composition shaping the released tSV.
C_LIO_LItSV and their steadily released counterparts have a different composition. TSV show a higher enrichment of effectors including immune receptors, miR, RNA- and other nucleic acid-binding proteins, than EVs.
C_LI 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. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. 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. The DC1 domain protein BINUCLEATE POLLEN is required for pollen development in Arabidopsis thaliana Development of the male gametophyte is a tightly regulated process that requires precise control of cell division and gene expression. A relevant aspect to understand the events underlying pollen development regulation constitutes the identification and characterization of the genes required for this process. In this work we showed that the DC1 domain protein BINUCLEATE POLLEN (BNP) is essential for pollen development and germination. Pollen grains carrying the defective BNP allele failed to complete mitosis II and are impaired in pollen germination. By yeast two hybrid screening and bimolecular fluorescence complementation assays, we identified a set of BNP-interacting proteins. Among confirmed interactors we found NAC family transcriptional regulators Vascular plant One-Zinc finger 1 (VOZ1) and VOZ2. VOZ1 localization changes during pollen development, moving to the vegetative nucleus at the tricellular stage. We observed that this relocalization requires BNP, as in their absence in pollen from bnp/BNP plants, VOZ1 nuclear localization resulted impaired. As voz1voz2 double mutants showed the same developmental defect observed in bnp pollen grains, we propose that BNP requirement to complete microgametogenesis could be linked to its interaction with VOZ1/2 proteins. BNP could have a role of scaffold protein, recruiting VOZ1/2 to the endosomal system into assemblies that are required for their further translocation to the nucleus, where they act as transcriptional regulators. 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. Structural and Functional Analyses Explain Pea KAI2 Receptor Diversity and Reveal Stereoselective Catalysis During Signal Perception KAI2 are plant /{beta} hydrolase receptors, which perceive smoke-derived butenolide signals (karrikins) and putative endogenous, yet unidentified phytohormones (KAI2-ligands, KLs). The number of functional KAI2 receptors varies among plant species. It has been suggested that KAI2 gene duplication and sub-functionalization plays an adaptative role for diverse environments or ligand diversification by altering the receptor responsiveness to specific KLs. Legumes represent one of the largest families of flowering plants and contain many essential agronomic crops. Prior to legume diversification, KAI2 underwent duplication, resulting in KAI2A and KAI2B. Integrating plant genetics, ligand perception and enzymatic assays, and protein crystallography, we demonstrate that Pisum sativum KAI2A and KAI2B act as receptors and enzymes with divergent ligand stereoselectivity. KAI2B has a stronger affinity than KAI2A towards the KAI2-ligand (-)-GR24 and remarkably hydrolyses a broader range of substrates including the strigolactone-like isomer (+)-GR24. We determine the crystal structures of PsKAI2B in apo and butenolide-bound states. The biochemical and structural analyses as well as recorded mass spectra of KAI2s reveal a transient intermediate on the catalytic serine and a stable adduct on the catalytic histidine, further illuminating the role of KAI2 not only as receptors but also as bona fide enzymes. Our work uncovers the stereoselectivity of ligand perception and catalysis by evolutionarily diverged KAI2 receptors in KAR/KL signaling pathways and proposes adaptive sensitivity to KAR/KL and strigolactone phytohormones by KAI2B. Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin BackgroundEvolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation at many enhancers, distal regulatory elements that control spatio-temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific, code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

ResultsWe present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that is able to predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Then, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity and liver enhancers associated with liver regeneration have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

ConclusionThe framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements. Ultrahigh-Throughput Ambient MS: Direct Analysis at 22 Samples per Second by Infrared Matrix-Assisted Laser Desorption Electrospray Ionization Mass Spectrometry Infrared Matrix-Assisted Laser Desorption Electrospray Ionization (IR-MALDESI) mass spectrometry is an ambient-direct sampling method being developed for high-throughput, label-free, biochemical screening of large-scale compound libraries. Here, we report the development of an ultrahigh-throughput continuous motion IR-MALDESI sampling approach capable of acquiring data at rates up to 22.7 samples per second in a 384-well microtiter plate. At top speed, less than 1% analyte carryover is observed from well-to-well and signal intensity relative standard deviations (RSD) of 11.5% and 20.9% for 3 M 1-hydroxymidazolam and 12 M dextrorphan, respectively, are achieved. The ability to perform parallel kinetics studies on 384 samples with ~30s time resolution using an isocitrate dehydrogenase 1 (IDH1) enzyme assay is shown. Finally, we demonstrate the repeatability and throughput of our approach by measuring 115,200 samples from 300 microtiter plate reads consecutively over 5.54 hours with RSDs under 8.14% for each freshly introduced plate. Taken together, these results demonstrate the use of IR-MALDESI at sample acquisition rates that surpass other currently reported direct sampling mass spectrometry approaches used for high throughput compound screening.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/465730v3_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@5f916corg.highwire.dtl.DTLVardef@4473b4org.highwire.dtl.DTLVardef@82e8forg.highwire.dtl.DTLVardef@128b57b_HPS_FORMAT_FIGEXP  M_FIG C_FIG An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. Sexual signals persist over deep time: ancient co-option of bioluminescence for courtship displays in cypridinid ostracods Although the the diversity, beauty, and intricacy of sexually selected courtship displays command the attention of evolutionists, the longevity of these traits in deep time is poorly understood. Population-based theory suggests sexual selection could either lower or raise extinction risk, resulting in high or low persistence of lineages with sexually selected traits. Furthermore, empirical studies that directly estimate longevity of sexually selected traits are uncommon. Sexually selected signals - including bioluminescent courtship - originated multiple times during evolution, allowing empirical study of their longevity after careful phylogenetic and divergence time analyses. Here, we estimate the first transcriptome-based molecular phylogeny and divergence times of Cypridinidae. We report extreme longevity of bioluminescent courtship, a trait important in mate choice and probably under sexual selection. Our relaxed-clock estimates of divergence times coupled with stochastic character mapping show luminous courtship evolved only once in Cypridinidae in a Sub-Tribe we name Luxorina at least 151 Million Years Ago (Ma) from cypridinid ancestors that used bioluminescence only in anti-predator displays, defining a Tribe we name Luminini. This time-calibrated molecular phylogeny of cypridinids will serve as a foundation for integrative and comparative studies on the biochemistry, molecular evolution, courtship, diversification, and ecology of cypridinid bioluminescence. The persistence of luminous courtship for hundreds of millions of years indicates that rates of speciation within the group exceeded extinction risk, allowing for the persistence of a diverse clade of signalling species and that sexual selection did not cause rapid loss of associated traits. Ion Permeation, Selectivity, and Electronic Polarization in Fluoride Channels Fluoride channels (Fluc) export toxic F- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the permeation mechanism at the molecular level has remained elusive. Herein we have applied constant-pH molecular dynamics and free energy sampling methods to investigate fluoride permeation through a Fluc protein from Escherichia coli. We find that fluoride is facile to permeate in its charged form, i.e., F-, by traversing through a non-bonded network. The extraordinary F- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F- permeation rate calculated using an electronically polarizable force field is significantly more accurate compared to the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F- - Fluc interactions.

O_TEXTBOXSIGNIFICANCE

A comprehensive atomistic-level computational study is presented of the mechanism of fluoride permeation through fluoride channels. The mechanism of fluoride permeation of the F- anion is established and the microscopic determinants of F- selectivity revealed. The essential nature of electronic polarization during F- permeation is also demonstrated through the computational modeling.

C_TEXTBOX Co-option of immune effectors by the hormonal signalling system triggering metamorphosis in Drosophila melanogaster Insect metamorphosis is regulated by the production, secretion and degradation of two peripheral hormones: 20-hydroxyecdysone (ecdysone) and juvenile hormone (JH). In addition to their roles in developmental regulation, increasing evidence suggests that these hormones are involved in innate immunity processes, such as phagocytosis and the induction of antimicrobial peptide (AMP) production. AMP regulation includes systemic responses as well as local responses at surface epithelia that contact with the external environment. At pupariation, Drosophila melanogaster increases dramatically the expression of three AMP genes, drosomycin (drs), drosomycin-like 2 (drsl2) and drosomycin-like 5 (drsl5). We show that the systemic action of drs at pupariation is dependent on ecdysone signalling in the fat body and operates via the ecdysone downstream target, Broad-Z4. In parallel, ecdysone also regulates local responses, specifically through the activation of drsl2 expression in the gut. Finally, we confirm the relevance of this ecdysone dependent AMP expression for the control of bacterial persistence by showing that flies lacking drs expression in the fat body have higher bacterial persistence over metamorphosis. Together, our data emphasize the importance of the association between ecdysone signalling and immunity using in vivo studies and establish a new role for ecdysone at pupariation, which impacts developmental success by regulating the immune system in a stage-dependent manner. We speculate that this co-option of immune effectors by the hormonal system may constitute a pre-emptive mechanism to control bacterial numbers in the pupa at the core of metamorphosis evolution. Brain deletion of Serf2 shifts amyloid conformation in a mouse model Loss of protein homeostasis accelerates ageing and contributes to age-related brain diseases like Parkinsons and Alzheimers diseases. Disease-specific proteins become aggregation-prone and form amyloid deposits that characterize disease, but the biological mechanisms driving their transition in cells and tissues are poorly understood. MOAG-4/SERF accelerates aggregation and toxicity of disease-related proteins in worm- and cell models. Whether MOAG-4/SERF affects amyloid pathology in brain, however, has remained unknown. Here we show that brain-specific knockout of SERF2 is viable and shifts aggregation in a mouse model for amyloid-beta aggregation. Without affecting amyloid-beta levels, knockout of SERF2 redirects its assembly into structurally and morphologically distinct amyloid fibrils and deposits. Our results show that SERF modifies the structural outcome of the aggregation process in mammalian brain. With structural polymorphisms recently being associated with disease severity and progression, our results may help to understand the molecular basis for structural diversity and find possibilities for structure-based interventions.

HighlightsO_LILoss of SERF2 slows embryonic development and causes perinatal lethality
C_LIO_LISERF2 affects proliferation in a cell-autonomous fashion
C_LIO_LIBrain-specific Serf2 knockout does not affect viability or A{beta} production
C_LIO_LIBrain deletion of Serf2 shifts the amyloid conformation of A{beta}
C_LI Chromosome breaks in breast cancers occur near herpes tumor virus sequences and explain why cancer drugs stop working One explanation for why breast cancer has a relentless tendency to come back after treatment is that breast cancer phenomes are fluid, so it is difficult for therapy to target new phenotypes that arise. In breast cancer, viral interactions cause chromosome breakage and fusion, producing anaphase bridges containing two centromeres. Cell division then leads to more chromosome breaks, more mutations, and new cancer phenotypes. Analyses of public data from about 2100 breast cancers implicate Epstein-Barr virus (EBV or Human Herpes Virus 4) in the cycles that break chromosomes and produce anaphase bridges. EBV generates breakage cycles with highly fragmented chromosomes (chromothripsis) as an inherent feature. These episodes continually produce new cancer driver mutations. DNA at chromothripsis boundaries in breast cancers consistently matches EBV tumor variants. Breakpoints are often near known EBV binding sites, and breast cancer shares an EBV methylation signature with other known EBV cancers. Breast cancer breakpoints on all chromosomes cluster around the same positions as in cancers definitively associated with EBV infection (nasopharyngeal cancer and Burkitts lymphoma). The status of EBV in the microbiome and the presence of anaphase bridges are phenotypic markers of breast cancer that are not explicitly targeted by chemotherapy. The breast cancer phenome continually supplies new cancer driver genes that foster breast cancer recurrence.

SummaryBreaks in breast cancer chromosomes and in fragments of anaphase bridges occur mainly near regions that resemble Epstein-Barr virus (EBV) tumor variants

Extensive breaks are an inherent outcome of anaphase bridges and the reassembly of highly fragmented chromosomes continually generates new cancer driver mutations and new cancer phenotypes.

Treatment does not explicitly target the underlying breakage and reassembly cycles so breast cancer often comes back.

The breast cancer phenome is not fixed and can continuously generate new phenotypes

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=154 HEIGHT=200 SRC="FIGDIR/small/467751v4_ufig1.gif" ALT="Figure 1">
View larger version (38K):
org.highwire.dtl.DTLVardef@6537cforg.highwire.dtl.DTLVardef@132fa50org.highwire.dtl.DTLVardef@1cfcb04org.highwire.dtl.DTLVardef@11979aa_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. A novel Boolean network inference strategy to model early hematopoiesis aging Hematopoietic stem cell (HSC) aging is a multifactorial event that leads to changes in HSC properties and function. These changes are intrinsically coordinated and affect the early hematopoiesis, involving hematopoietic stem and progenitor cells (HSPCs). The objective of this work is to better understand the mechanisms and factors controlling these changes. We have therefore developed an original strategy to construct a Boolean network of genes explaining the priming and homeostasis of HSCs (graphical abstract). Based on our previous scRNA-seq data, we performed an exhaustive analysis of the transcriptional network and identified active transcription modules or regulons along the differentiation trajectory of selected HSPC states. This global view of transcriptional regulation led us to focus on 15 components, 13 selected TFs (Tal1, Fli1, Gata2, Gata1, Zfpm1, Egr1, Junb, Ikzf1, Myc, Cebpa, Bclaf1, Klf1, Spi1) and 2 complexes regulating the ability of HSC to cycle (CDK4/6 - Cyclin D and CIP/KIP). We then defined the connections controlling the differentiation dynamics of HSC states and constructed an influence graph between the TFs involved in the dynamics by mixing observations from our scRNA-seq data and knowledge from the literature. Then, using answer set programming (ASP) and in silico perturbation analysis, we obtained a Boolean model which is the solution of a Boolean satisfiability problem. Finally, perturbation of the model based on age-related changes revealed important regulations, such as the overactivation of Egr1 and Junb or the loss of Cebpa activation by Gata2, which were found to be relevant for the myeloid bias of aged HSC. Our work shows the efficiency of the combination of manual and systematic methods to elaborate a Boolean model. The developed strategy led to the proposal of new regulatory mechanisms underlying the differentiation bias of aged HSCs, explaining the decreased transcriptional priming of HSCs to all mature cell types except megakaryocytes. In vitro susceptibility to β-lactam antibiotics and viability of Neisseria gonorrhoeae strains producing plasmid-mediated broad- and extended-spectrum β-lactamases N. gonorrhoeae strains producing plasmid-mediated broad- and extended-spectrum {beta}-lactamases (ESBLs) have been obtained in vitro and studied for their viability and susceptibility to {beta}-lactam antibiotics. The artificial pblaTEM-1 and pblaTEM-20 plasmids were constructed by site-directed mutagenesis from a pblaTEM-135 plasmid of the Toronto/Rio type, which was extracted from the clinical isolate. MIC values were determined for a series of {beta}-lactam antibiotics, including benzylpenicillin, ampicillin, cefuroxime, ceftriaxone, cefixime, cefotaxime, cefepime, meropenem, imipenem, and doripenem. The N. gonorrhoeae strain carrying the pblaTEM-20 plasmid exhibited a high level of resistance to penicillins and II-IV generation cephalosporins (MIC [&ge;] 2 mg/L) but not to carbapenems (MIC [&le;] 0.008 mg/L). However, this strain stopped growing after 6 hours of cultivation. The reduced viability was not associated with the plasmid loss but can be explained by the following factors: the presence of the plasmid itself, which requires additional costs for its reproduction, and the expression of ESBL, which can affect the structure of the peptidoglycan layer in the cell membrane. The cell growth was mathematically modeled using the generalized Verhulst equation. The parameter, characterizing the maximum possible number of cells grown under given conditions, decreased within the wild type (without plasmids) - pblaTEM-135 - pblaTEM-1 - pblaTEM-20 series, i.e., the plasmid-bearing strains had reduced viability compared to that of the wild-type strain. The kinetics for the cell death of N. gonorrhoeae strains without the pblaTEM-20 plasmid in the presence of ceftriaxone can be described by a modified Chick-Watson law. For the N. gonorrhoeae strain that harbors the pblaTEM-20 plasmid, the kinetics for cell number changes with time reflected several processes: the hydrolysis of ceftriaxone by TEM-20 {beta}-lactamase, growth and gradual death of cells. The demonstrated reduced viability of N. gonorrhoeae strains with the pblaTEM-20 plasmid probably explains the absence of clinical isolates of N. gonorrhoeae that produce ESBL. Accounting for stimulations that do not elicit motor-evoked potentials when mapping cortical representations of multiple muscles The representation of muscles in the cortex can be mapped using navigated transcranial magnetic stimulation. The mapping can be quantified via measures like the centre of gravity or the size of the region of excitability. Determining these measures typically relies only on stimulation points that yield motor-evoked potentials (MEPs); stimulations that do not elicit an MEP, i.e., non-MEP points, are ignored entirely. In this study, we show how incorporating non-MEP points may affect the estimates of the size and centroid of the excitable area in eight hand and forearm muscles after mono-phasic single-pulse TMS. We performed test-retest assessments in twenty participants and estimated the reliability of centroids and sizes of the corresponding areas using inter-class correlation coefficients. For most muscles, the reliability turned out good, if not excellent. As expected, removing the non-MEP points significantly decreased area sizes and area weights, suggesting that conventional approaches that do not account for non-MEP points are likely to overestimate the regions of excitability. Dynamic Phosphorylation of RGS Provides Spatial Regulation of MAP Kinase And Promotes Completion of Cytokinesis during the Yeast Pheromone Response Yeast use a G-protein coupled receptor (GPCR) signaling pathway to detect mating pheromone, arrest in G1, and direct polarized growth towards the potential mating partner. The primary negative regulator of this pathway is the regulator of G-protein signaling (RGS), Sst2, which induces G GTPase activity and subsequent inactivation of all downstream signaling including a MAPK cascade. The MAPK Fus3 phosphorylates the RGS in response to pheromone, but the role of this modification is unknown. We set out to examine the role of RGS phosphorylation during the pheromone response. We found that RGS phosphorylation peaks early in the pheromone response and diminishes RGS localization to the polarity site and focuses MAPK complexes there. At later time points, RGS is predominantly unphosphorylated, which promotes RGS localization to the polar cap and broadens the distribution of MAPK complexes relative to the Cdc42 polarity machinery. Surprisingly, we found that phosphorylation of the RGS is required for the completion of cytokinesis prior to pheromone induced growth. The completion of cytokinesis in the presence of pheromone is promoted by the formin Bnr1 and the kelch-repeat protein, Kel1, both proteins previously found to interact with the RGS. Artificial reefs geographical location matters more than shape, age and depth for sessile invertebrate colonization in the Gulf of Lion (NorthWestern Mediterranean Sea) AO_SCPLOWBSTRACTC_SCPLOWArtificial reefs (ARs) have been used to support fishing activities. Sessile invertebrates are essential components of trophic networks within ARs, supporting fish productivity. However, colonization by sessile invertebrates is possible only after effective larval dispersal from source populations, usually in natural habitat. While most studies focused on short term colonization by pioneer species, we propose to test the relevance of geographic location, shape, age and depth of immersion on the ARs long term colonization by species found in natural stable communities in the Gulf of Lion. We recorded the presence of five sessile invertebrates species, with contrasting life history traits and regional distribution in the natural rocky habitat, on ARs with different shapes deployed during two immersion time periods (1985 and the 2000s) and in two depth ranges (<20m and >20m). At the local level ([~]5kms), neither shape, depth nor immersion duration differentiated ARs assemblages. At the regional scale (>30kms), colonization patterns differed between species, resulting in diverse assemblages. This study highlights the primacy of geographical positioning over shape, immersion duration and depth in ARs colonization, suggesting it should be accounted for in maritime spatial planning. A chromosome-level assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and improves time of divergence estimates for the Cannabaceae family O_LIWe present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly resulted in fragmented assemblies. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur.
C_LIO_LIThe chromosome-level assembly of Cascade was developed by scaffolding the previously-reported Cascade assembly generated with PacBio long-read sequencing, and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations, and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assess molecular evolution in gene sequences, gene family expansion and contraction, and time divergence using Bayesian inference.
C_LIO_LIWe identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families.
C_LIO_LIThe Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.
C_LI