Genetic diversity and connectivity of the Ostreid herpesvirus 1 populations in France: a first attempt to phylogeographic inference for a marine mollusc disease The genetic diversity of viral populations is a key driver of the spatial and temporal diffusion of viruses; yet, studying the diversity of whole genomes from natural populations still remains a challenge. Phylodynamic approaches are commonly used for RNA viruses harboring small genomes, but have only rarely been applied to DNA viruses with larger genomes. Here, we used the Pacific Oyster Mortality Syndrome (POMS, a disease that affects oyster farms around the world) as a model to study the genetic diversity of its causative agent, the Ostreid herpesvirus 1 (OsHV-1) in the three main French oyster-farming areas. Using ultra-deep sequencing on individual moribund oysters and an innovative combination of bioinformatics tools, we de novo assembled 21 OsHV-1 new genomes. Combining quantification of major and minor genetic variations, phylogenetic analysis and ancestral state reconstruction of discrete traits approaches; we assessed the connectivity of OsHV-1 viral populations between the three oyster-farming areas. Our results suggest that the Marennes- Oleron Bay represents the main source of OsHV-1 diversity, from where the virus has dispersed to other farming areas, a scenario consistent with current practices of oyster transfers in France. We demonstrate that phylodynamic approaches can be applied to aquatic DNA viruses to determine how epidemiological, immunological, and evolutionary processes act and potentially interact to shape their diversity patterns.

ImportancePhylogeography is a field of research that attempts to reconstruct the relationships between individual genotypes within a species and then correlate these genealogical relationships with their geographic and temporal origin. This field of research has become an essential step in the understanding of pandemics, in particular to determine the origin, spread and evolution of a pathogen as currently illustrated in studies on viral pandemics. However, because phylogeographic analyses are based on genome variation, stable genomes yield less information than labile genomes. Accordingly, viruses with double-stranded DNA (dsDNA) genomes generally have lower nucleotide diversity than RNA viruses. In this study, by combining the use of both major and minor genetic variations with phylogeographic analyses of the oyster herpesvirus OsHV-1, we highlight genealogical relationships that are not depicted in phylogenetic trees based on consensus viral genomes only. These data offer a plausible scenario reflecting the origin and spread of OsHV-1 populations between oyster- farming sites.
The activity and functions of soil microbial communities in the Finnish sub-Arctic vary across vegetation types Increased microbial activity in high-latitude soils due to climate change might lead to higher greenhouse gas (GHG) emissions. However, mechanisms of microbial GHG production and consumption in tundra soils are not thoroughly understood. We analyzed 116 soil metatranscriptomes from 73 sites in the Finnish sub-Arctic to investigate how the diversity and functional potential of bacterial and archaeal communities vary across vegetation types and soil layers. Soils differed in physicochemical conditions, with meadow soils being characterized by higher pH and low soil organic matter (SOM) and carbon/nitrogen ratio whereas dwarf shrub-dominated ecosystems with high SOM and low pH. Actinobacteria, Acidobacteria, Alphaproteobacteria, and Planctomycetes predominated all communities but there were significant differences on genus level between vegetation types, as plant polymer degrading groups were more active in shrub-dominated soils compared to meadows. Given that climate change scenarios predict expansion in dwarf shrubs at high latitudes, our results indicate that the rate of carbon turnover in tundra soils may increase in the future. Additionally, transcripts of methanotrophs were detected in the mineral layer of all soils, potentially moderating methane fluxes from deeper layers. In all, this study provides new insights into possible shifts in tundra microbial diversity and activity with climate change.
14-3-3ζ constrains insulin secretion by regulating mitochondrial function in pancreatic β-cells While critical for neurotransmitter synthesis in the brain, members of the 14-3-3 protein family are often assumed to have redundant, over-lapping roles due to their high sequence homology and ubiquitous expression. Despite this assumption, various mammalian 14-3-3 isoforms have now been implicated in regulating cellular and organismal metabolism; however, these functions were primarily observed in cell lines or from systemic knockout mouse models. To date, we have begun to define the contributions of 14-3-3{zeta} in adipocytes, but whether 14-3-3{zeta} has additional metabolic roles in other cell types, such as the pancreatic {beta}-cell, is unclear. We previously documented a pro-survival role of 14-3-3{zeta} in MIN6 insulinoma cells, as depletion of 14-3-3{zeta} induced cell death, but paradoxically, whole-body deletion of 14-3-3{zeta} in mice resulted in significantly enlarged {beta}-cell area with no effects on insulin secretion. To better understand the role of 14-3-3{zeta} in {beta}-cells, we generated {beta}-cell-specific 14-3-3{zeta} knockout ({beta}14-3-3{zeta}KO) mice, and while no differences in {beta}-cell mass were observed, {beta}14-3-3{zeta}KO mice displayed potentiated insulin secretion due to enhanced mitochondrial function and ATP synthesis. Deletion of 14-3-3{zeta} led to profound changes to the {beta}-cell transcriptome, where pathways associated with mitochondrial respiration and oxidative phosphorylation were upregulated. Acute treatment of mouse islets and human islets with pan-14-3-3 inhibitors recapitulated the potentiation in glucose-stimulated insulin secretion (GSIS) and mitochondrial function, suggesting that 14-3-3{zeta} is a critical isoform in {beta}-cells that regulates GSIS. In dysfunctional db/db islets and islets from type 2 diabetic donors, expression of Ywhaz/YWHAZ, the gene encoding 14-3-3{zeta}, was inversely associated with insulin secretory capacity, and pan-14-3-3 protein inhibition was capable of enhancing GSIS and mitochondrial function. Taken together, this study demonstrates important regulatory functions of 14-3-3{zeta} and its related isoforms in insulin secretion and mitochondrial function in {beta}-cells. A deeper understanding of how 14-3-3{zeta} influences {beta}-cell function will further advance our knowledge of how insulin secretion from {beta}-cells is regulated. ScRNA-seq Expression of IFI27 and APOC2 Identifies Four Alveolar Macrophage Superclusters in Healthy BALF Alveolar macrophages (AMs) reside on the luminal surface of the airways and alveoli, ensuring proper gas exchange by ingesting cellular debris and pathogens, and regulating inflammatory responses. Therefore, understanding the heterogeneity and diverse roles played by AMs and recruited monocytes is critical for treating airway diseases. We performed single-cell RNA sequencing on 113,213 bronchoalveolar lavage cells from four healthy and three uninflamed cystic fibrosis subjects and identified four AM superclusters (families) based on the expression of IFI27 and APOC2 genes. These 4 AM families have at least eight distinct functional members (subclusters) named after their differentially expressed gene(s): IGF1, CCL18, CXCL5, Cholesterol, Chemokine, Metallothionein, Interferon and small-cluster AMs. Interestingly, the Chemokine cluster further divides with each subcluster selectively expressing a unique combination of chemokines. One of the most striking observations, besides the heterogeneity, is the conservation of AM family members in relatively equal ratio across all AM superclusters and individuals. Finally, transcriptional data and TotalSeq technology were used to investigate cell surface markers that distinguish resident AMs from recruited monocytes. Overall, functional specializations of the AM subclusters suggest that there are AM niches with defined programming states, highlighting a clear division of labor. Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

MethodsAnti-GD2 CAR T cells were generated using CRISPR/Cas9 within nine days using recombinant Cas9 protein and nucleic acids, without any viral vectors or animal components. The CAR was specifically targeted to the T Cell Receptor Alpha Constant gene (TRAC). T cell products were characterized at the level of the genome, transcriptome, proteome, and secretome using CHANGE-seq, scRNA-seq, spectral cytometry, and ELISA assays. Functionality was evaluated in vivo in an NSG xenograft neuroblastoma model.

ResultsIn comparison to traditional retroviral CAR T cells, virus-free CRISPR CAR (VFC-CAR) T cells exhibit TRAC-targeted genomic integration of the CAR transgene, elevation of transcriptional and protein characteristics associated with a memory phenotype, and low tonic signaling prior to infusion arising in part from the the knockout of the TCR. Upon exposure to the GD2 target antigen, anti-GD2 VFC-CAR T cells exhibited specific cytotoxicity against GD2+ cells in vitro and induced solid tumor regression in vivo, with robust homing, persistence, and low exhaustion against a human neuroblastoma xenograft model.

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors. Non-canonical Wnt/Ror2 signaling regulates cell-matrix crosstalk to prompt directional tumor cell invasion and dissemination in breast cancer Bidirectional cell-extracellular matrix (ECM) interactions represent fundamental exchanges during tumor progression. We demonstrate the noncanonical Wnt receptor, Ror2, regulates tumor cell-driven matrix remodeling in models of breast cancer. Wnt/Ror2 loss-of-function triggers tumor cell invasion, accompanied by changes in actin cytoskeleton, adhesion, and collagen crosslinking gene expression programs. E-cadherin levels decline upon Ror2 depletion, and spatially, we pinpoint the upregulation and redistribution of 5 and {beta}3 integrins together with the production of fibronectin in areas of invasion. Wnt/{beta}-catenin-dependent and Wnt/Ror2 alternative Wnt signaling appear to regulate distinct functions for tumor cells regarding their ability to modify cell-ECM exchanges during invasion. Furthermore, blocking either integrin or focal adhesion kinase (FAK), a downstream mediator of integrin-mediated signal transduction, abrogates the enhanced migration observed upon Ror2 loss. These results reveal a critical function for the alternative Wnt receptor, Ror2, as a determinant of reciprocal communication between tumor cells and ECM during cancer invasion and metastasis. X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles, and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here we use wholegenome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 11% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution. Native proline-rich motifs exploit sequence context to target actin-remodeling Ena/VASP proteins The human proteome is replete with short linear motifs (SLiMs) of 4-6 residues that are critical for protein-protein interactions, yet the importance of the sequence surrounding such motifs is underexplored. We devised a proteomic screen to systematically examine the influence of SLiM sequence context on protein-protein interactions. Focusing on the EVH1 domain of ENAH, an actin regulator that is upregulated in invasive cancers, we screened 36-residue proteome-derived peptides for binding. We discovered a pocket on the ENAH EVH1 domain that diverged from its orthologs to recognize extended SLiMs, and we found that proteins with two EVH1-binding SLiMs can wrap around a single domain. We also found that the ciliary protein PCARE uses an extended 23-residue region to obtain higher affinity than any known ENAH EVH1-binding motif. Our screen provides a way to uncover the effects of broader proteomic context on motif-mediated interactions, revealing diverse mechanisms of contextual control over EVH1 interactions and establishing that SLiMs cant be fully understood outside of their native context. Game of theories: an agent-based model of how variation in contest cost-benefits determines the evolution of contest resolution rules Game-theory based models are used to understand rules that animals use to settle contests over indivisible resources. However, the empirical literature of contests indicates controversial support to models, with some species supporting different models and other species showing no support to any model. Since strategies used to resolve contests may have different associated costs, it is possible that different conditions have determined the evolution of distinct assessment strategies used by animals. We used an agent-based model to explore the importance of the following conditions: resource availability, probability of reproduction with resource, and damage costs on evolution of assessment strategies. We used self- and mutual-assessment models as a heurist framework to build agents with different assessment strategies. In our model, agents competed for resources in scenarios with different combinations of resource availability, probability of reproduction with resource, and damage costs. We found that agents following the self-assessment with damage strategy were prevalent in scenarios with no probability of reproduction without the resource, independently of other variables. We also found that agents following the non-aggressive strategy occurred in all scenarios. However, agents using the non-aggressive strategy were prevalent only in scenarios with probability of reproduction with the resource. Finally, we observed that agents using mutual-assessment occurred only in a scenario with high risk of damage, low availability of resources, and with probability of reproduction without the resource. These results indicate that agents following the self-assessment with damage and non-aggressive strategies may be able to stay at most scenarios. The effect of recent habitat change on genetic diversity at putatively adaptive and neutral loci in Primula veris in semi-natural grasslands Recent habitat change in semi-natural grasslands due to a lack of management has been shown to affect the genetic diversity of grassland plants. However, it is unknown how a change in local environment affects genetic diversity at adaptive loci. We applied RADseq (restriction-site associated DNA sequencing) to extract > 3,000 SNPs across 568 individuals from 32 Estonian populations of Primula veris, a plant species common to semi-natural grasslands. We evaluated the effect of recent grassland overgrowth following management abandonment on the genetic diversity at putatively neutral and adaptive loci, which we distinguished by applying three methods, i.e., linear and categorical environmental association analyses, and an FST outlier test. For validation, we randomised the genotype to sample assignments. Effects of recent habitat change on genetic diversity differed between neutral and adaptive SNP sets. Genetic diversity assessed at putatively neutral loci was similar in open and overgrown habitats but showed a significant difference between these habitat types at putatively adaptive loci: overgrown (i.e. newly established) habitats exhibited higher genetic diversity at putatively adaptive loci than open (i.e. old) habitats, likely due to the exertion of novel selection pressures imposed by new habitat conditions. This increase in genetic diversity at putatively adaptive loci in the new environment points to currently ongoing selection processes where genetic adaptation to the old habitat is potentially lost through altered allele frequencies. Our study suggests that a recent change in local habitat conditions may not be reflected in neutral loci whereas putatively adaptive loci can inform about potential selection processes. Trial and error: a hierarchical modeling approach to test-retest assessment The concept of test-retest reliability indexes the consistency of a measurement across time. High reliability is critical for any scientific study, but specifically for the study of individual differences. Evidence of poor reliability of commonly used behavioral and functional neuroimaging tasks is mounting. Reports on low reliability of task-based fMRI have called into question the adequacy of using even the most common, well-characterized cognitive tasks with robust population-level effects, to measure individual differences. Here, we lay out a hierarchical framework that estimates reliability as a correlation divorced from trial-level variability, and show that reliability estimates tend to be higher compared to the conventional framework that adopts condition-level modeling and ignores across-trial variability. We examine how estimates from the two frameworks diverge and assess how different factors (e.g., trial and subject sample sizes, relative magnitude of cross-trial variability) impact reliability estimates. We also show that, under specific circumstances, the two statistical frameworks converge. Results from the two approaches are approximately equivalent if (a) the trial sample size is sufficiently large, or (b) cross-trial variability is in the same order of magnitude as, or less than, cross-subject variability. As empirical data indicate that cross-trial variability is large in most tasks, this work highlights that a large number of trials (e.g., greater than 100) may be required to achieve precise reliability estimates. We reference the tools TRR and 3dLMEr for the community to apply trial-level models to behavior and neuroimaging data and discuss how to make these new measurements most useful for current studies. Estimation of differential cell cycle kinetics in higher plant root meristem with cellular fate and positional resolution. Plant root development is a complex spatial-temporal process that originates in the root apical meristem (RAM). To shape the organs structure signaling between the different cells and cell files must be highly coordinated. Thereby, diverging kinetics of chromatin remodeling and cell growth in these files need to be integrated and balanced by differential cell growth and local differences in cell proliferation frequency. Understanding the local differences in cell cycle duration in the RAM and its correlation with chromatin organization is crucial to build a holistic view on the different regulatory processes and requires a quantitative estimation of the chromatin geometry and underlying mitotic cell cycle phases timing at every cell file and every position. Unfortunately, so far precise methods for such analysis are missing.

This study presents a robust and straightforward pipeline to determine in parallel the duration of cell cycles key stages in all cell layers of a plants root and their nuclei organization. The methods combine marker-free techniques based on the detection of the nucleus, deep analysis of the chromatin phase transition, incorporation of 5-ethynyl-2'-deoxyuridine (EdU), and mitosis with a deep-resolution plant phenotyping platform to analyze all key cell cycle events kinetics.

In the Arabidopsis thaliana L. RAM S-phase duration was found to be as short as 20-30 minutes in all cell files. The subsequent G2-phase duration however depends on the cell type/position and varies from 3.5 hours in the pericycle to more than 4.5 hours in the epidermis. Overall, S+G2+M duration in Arabidopsis under our condition is 4 hours in the pericycle and up to 5.5 hours in the epidermis.

Endocycle duration was determined as the time required to achieve 100% EdU index in the transition zone and estimated to be in the range of 3-4 hours.

Besides Arabidopsis, we show that the presented technique is applicable also to root tips of other dicot and monocot plants (tobacco (Nicotiana tabacum L.), tomato (Lycopersicon esculentum L.) and wheat (Triticum aestivum L.). 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. Desalination and temperature increase will shift seasonal grazing patterns of invasive Gammarus tigrinus on charophytes Charophytes are a refuge for zooplankton and stabilize sediments, but they are also a food source for various animal species (water birds, fishes, invertebrates). Especially the introduction of new species, like Gammarus tigrinus, into the Baltic Sea led to yet not understood changes in the food web. Furthermore, future projections point to increased water temperatures at lowered salinity levels affecting species capacity to acclimatize to changing abiotic factors. In this study we investigated the influence of temperature and salinity on the grazing pressure of Gammarus tigrinus on two charophyte species: Chara aspera and Chara tomentosa. The grazing experiments were conducted in a full factorial design with the factors salinity (3 - 13 g kg-1), temperature (5 - 30 {degrees}C), and charophyte species. Grazing rates were determined as mass deviation within 48 hours considering biomass changes in the presence and absence of gammarids. Grazing rate were further used to calculate charophyte losses in two coastal lagoons with different salinity concentrations for recent and future time periods. The potential grazing peak of about 24 {degrees}C is not yet reached in these coastal waters but may be reached in the near future as shown by our future projection results. However, the temperature increase, and desalination will cause a shift in seasonal individual grazing patterns from summer to spring and autumn. Desalination and temperature increase can lead to a shift in optimal habitats for G. tigrinus in the future. A dynamical systems treatment of transcriptomic trajectories in hematopoiesis Inspired by Waddingtons illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein the dynamics of an exogenous signal couples to a cells enormously complex signaling and transcriptional machinery, eliciting a qualitative transition in the collective state of a cell - its fate. It remains unclear, however, whether the dynamical systems framework can go beyond a word-based caricature of the system and provide sharp quantitative insights that further our understanding of differentiation. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the individual concentration variations of a myriad of genes. Here, starting from formal mathematical derivations, we challenge these transcriptomic trajectories to a rigorous statistical evaluation of whether they display signatures consistent with bifurcations. After pinpointing bifurcations along transcriptomic trajectories of the neutrophil branch of hematopoeitic differentiation we are able to further leverage the primitive features of a linear instability to identify the single-direction in gene expression space along which the bifurcation unfolds and identify possible gene contributors. This scheme identifies transcription factors long viewed to play a crucial role in the process of neutrophil differentiation in addition to identifying a host of other novel genetic players. Most broadly speaking, we provide evidence that, though very high-dimensional, a bifurcating dynamical systems formalism might be appropriate for the process of cellular differentiation and that it can be leveraged to provide insights. Ambitiously, our work attempts to take a step beyond data-analysis and towards the construction of falsifiable mathematical models that describe the dynamics of the entire transcriptome. Distinguishing excess mutations and increased cell death based on variant allele frequencies Tumors often harbor orders of magnitude more mutations than healthy tissues. The increased number of mutations may be due to an elevated mutation rate or frequent cell death and correspondingly rapid cell turnover, or a combination of the two. It is difficult to disentangle these two mechanisms based on widely available bulk sequencing data, where sequences from individual cells are intermixed and, thus, the cell lineage tree of the tumor cannot be resolved. Here we present a method that can simultaneously estimate the cell turnover rate and the rate of mutations from bulk sequencing data. Our method works by simulating tumor growth and finding the parameters with which the observed data can be reproduced with maximum likelihood. Applying this method to a real tumor sample, we find that both the mutation rate and the frequency of death may be high.

Author SummaryTumors frequently harbor an elevated number of mutations, compared to healthy tissue. These extra mutations may be generated either by an increased mutation rate or the presence of cell death resulting in increased cellular turnover and additional cell divisions for tumor growth. Separating the effects of these two factors is a nontrivial problem. Here we present a method which can simultaneously estimate cell turnover rate and genomic mutation rate from bulk sequencing data. Our method is based on the estimation of the parameters of a generative model of tumor growth and mutations. Applying our method to a human hepatocellular carcinoma sample reveals an elevated per cell division mutation rate and high cell turnover. High-quality Arabidopsis thaliana Genome Assembly with Nanopore and HiFi Long Reads Arabidopsis thaliana is an important and long-established model species for plant molecular biology, genetics, epigenetics, and genomics. However, the latest version of reference genome still contains significant number of missing segments. Here, we report a high-quality and almost complete Col-0 genome assembly with two gaps (Col-XJTU) using combination of Oxford Nanopore Technology ultra-long reads, PacBio high-fidelity long reads, and Hi-C data. The total genome assembly size is 133,725,193 bp, introducing 14.6 Mb of novel sequences compared to the TAIR10.1 reference genome. All five chromosomes of Col-XJTU assembly are highly accurate with consensus quality (QV) scores > 60 (ranging from 62 to 68), which are higher than those of TAIR10.1 reference (QV scores ranging from 45 to 52). We have completely resolved chromosome (Chr) 3 and Chr5 in a telomere-to-telomere manner. Chr4 has been completely resolved except the nucleolar organizing regions, which comprise long repetitive DNA fragments. The Chr1 centromere (CEN1), reportedly around 9 Mb in length, is particularly challenging to assemble due to the presence of tens of thousands of CEN180 satellite repeats. Using the cutting-edge sequencing data and novel computational approaches, we assembled about 4 Mb of sequence for CEN1 and a 3.5-Mb-long CEN2. We investigated the structure and epigenetics of centromeres. We detected four clusters of CEN180 monomers, and found that the centromere-specific histone H3-like protein (CENH3) exhibits a strong preference for CEN180 cluster 3. Moreover, we observed hypomethylation patterns in CENH3-enriched regions. We believe that this high-quality genome assembly, Col-XJTU, would serve as a valuable reference to better understand the global pattern of centromeric polymorphisms, as well as genetic and epigenetic features in plants. Mitochondrial double-stranded RNAs as a pivotal mediator in the pathogenesis of Sjogrens syndrome ObjectiveSjogrens syndrome (SS) is a systemic autoimmune disease that targets the exocrine glands, resulting in impaired saliva and tear secretion. To date, type I interferons (I-IFNs) are increasingly recognized as pivotal mediators in SS, but their endogenous drivers have not been elucidated. This study investigates the role of mitochondrial double-stranded RNAs (mt-dsRNAs) in regulating I-IFN response in SS.

MethodsSaliva and tear from SS patients and controls (n=73 for saliva and n=16 for tear), the salivary glands of the SS-prone-non-obese-diabetic mouse, and primary human salivary glandular cells were screened for mt-dsRNAs by RT-qPCR. The human salivary cell line (NS-SV-AC) grown as three-dimensional spheroids were subject to dsRNA stress to measure mt-dsRNA induction and recapitulation of SS glandular inflammatory features. Acetylcholine, SS-IgG, upadacitinib (JAK1 inhibitor), or 2-C'-methyladenosine (mitochondrial transcription inhibitor) were applied to characterize the roles of mt-dsRNAs. To identify endogenous dsRNA-sensor and confirm the mitochondrial origin of cytoplasmic dsRNAs, the immunoprecipitation of dsRNAs was performed.

Resultsmt-dsRNAs were elevated in the SS specimens with salivary ND5 and tear CYTB1 being statistically associated with secretory dysfunction/inflammation and corneal/conjunctival damage, respectively. Stimulation of the spheroids with dsRNA stress of poly I:C induced mt-dsRNAs, p-PKR, and I-IFNS via the JAK1/STAT pathway whereas the inhibition of mt-RNA synthesis or JAK1 attenuated the glandular signature. The inhibitory effect of acetylcholine on mt-dsRNAs and I-IFNS induction was reversed by SS-IgG.

Conclusionmt-dsRNAs amplify the impact of dsRNA stress on SS glandular signatures in vitro, potentially propagating a pseudo-viral signal in the SS target tissue.

SummaryMitochondrial double-stranded RNA levels were elevated in the tear and saliva of SS patients, which was associated with secretory dysfunction and tissue inflammation. These RNAs amplified type I interferon signature as well as glandular phenotypes reported in SS. Inhibitors of mitochondrial RNA transcription or JAK1 in salivary gland acinar cell spheroids attenuated the mitochondrial RNA-mediated changes. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic re-organization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa and Pseudomonas putida. We combined growth dynamics and dynamic changes of intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis and growth and uncovered fundamental limitations of the regulatory network: after nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Graphical synopsis O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Standfirst textMicrobes face a series of fundamental trade-offs that limit their ability to optimize simultaneously for both glycolytic and gluconeogenic growth.

Bullet pointsO_LILag times between glycolysis and gluconeogenesis show asymmetry in many microbes: A long lag in one direction, but a short lag in the other.
C_LIO_LILong lag times are caused by an inability to sense fluxes after nutrient shifts.
C_LIO_LIWith existing regulation, lag time asymmetry can only be overcome by reducing either growth rate or increasing futile cycling in metabolism.
C_LI Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. A natural mutator allele shapes mutation spectrum variation in mice Summary ParagraphAlthough germline mutation rates and spectra can vary within and between species, genetic modifiers of these traits have long eluded detection. In this study, we searched for loci that influence germline mutagenesis using a uniquely powerful resource: a panel of recombinant inbred mouse lines known as the BXD, descended from the laboratory mouse strains C57BL/6J (B) and DBA/2J (D). Each BXD lineage has been maintained by brother-sister mating in the near absence of natural selection, accumulating de novo mutations for up to 50 years on a known genetic background that is a unique linear mosaic of B and D haplotypes. We show that mice inheriting D haplotypes at a quantitative trait locus (QTL) on chromosome 4 accumulate C>A germline mutations at a 50% higher rate than those inheriting B haplotypes, primarily due to the activity of a C>A-dominated mutational signature known as SBS18. The B and D QTL haplotypes encode different alleles of the DNA repair gene Mutyh, which underlies the heritable colorectal cancer syndrome in which SBS18 was first identified. The B and D Mutyh alleles are present in wild populations of Mus musculus domesticus, providing evidence that common genetic variation modulates germline mutagenesis in a model mammalian species.