Modelling and assessing additional transmission routes for porcine reproductive and respiratory syndrome virus: vehicle movements and feed ingredients Accounting for multiple modes of livestock disease dissemination in epidemiological models remains a challenge. We developed and calibrated a mathematical model for transmission of porcine reproductive and respiratory syndrome virus (PRRSV), tailored to fit nine modes of between-farm transmission pathways including: farm-to-farm proximity (local transmission), contact network of batches of pigs transferred between farms (pig movements), re-break probabilities for farms with previous PRRSV outbreaks, with the addition of four different contact networks of transportation vehicles (vehicles to transport pigs to farms, pigs to markets, feed and crew) and the amount of animal by-products within feed ingredients (e.g. animal fat or meat and bone meal). The model was calibrated on weekly PRRSV outbreaks data. We assessed the role of each transmission pathway considering the dynamics of specific types of production (i.e., sow, nursery). Although our results estimated that the networks formed by transportation vehicles were more densely connected than the network of pigs transported between-farms, pig movements and farm proximity were the main PRRSV transmission routes regardless of farm types. Among the four vehicle networks, vehicles transporting pigs to farms explained a large proportion of infections, sow = 20.9%; nursery = 15%; and finisher = 20.6%. The animal by-products showed a limited association with PRRSV outbreaks through descriptive analysis, and our model results showed that the contribution of animal fat contributed only 2.5% and meat and bone meal only 0.03% of the infected sow farms. Our work demonstrated the contribution of multiple routes of PRRSV dissemination, which has not been deeply explored before. It also provides strong evidence to support the need for cautious, measured PRRSV control strategies for transportation vehicles and further research for feed by-products modeling. Finally, this study provides valuable information and opportunities for the swine industry to focus effort on the most relevant modes of PRRSV between-farm transmission.
Confronting population models with experimental microcosm data: from trajectory matching to state-space models O_LIPopulation and community ecology traditionally has a very strong theoretical foundation with well-known models, such as the logistic and its many variations, and many modification of the classical Lotka-Volterra predator-prey and interspecific competition models. More and more, these classical models are confronted to data via fitting to empirical time-series for purposes of projections or for estimating model parameters of interest. However, the interface between mathematical population or community models and data, provided by a statistical model, is far from trivial.
C_LIO_LIIn order to help empiricists, especially researchers working with experimental laboratory populations in micro- and mesocosms, make informed decisions, we here ask which error structure one should use when fitting classical deterministic ODE models to empirical data, from single species to community dynamics and trophic interactions. We use both realistically simulated data and empirical data from microcosms to answer this question in a Bayesian framework.
C_LIO_LIWe find that many model parameters can be estimated precisely with an appropriate choice of error structure using pure observation error or state-space models, if observation errors are not too high. However, Allee effect models are typically hard to identify and state-space models should be preferred with when model complexity increases.
C_LIO_LIOur work shows that, at least in the context of experimental laboratory populations, deterministic models can be used to describe stochastic population dynamics that include process variability and observation error. Also, many scenarios do not require a complex state-space model formulation and simpler trajectory matching is sufficient for accurate parameter estimates. Finally, we provide a comprehensive tutorial for fitting these models in R.
C_LI
Dramatic differences in male and female mortality trends for selected European cohorts over ~20 years. Longevity is of considerable interest. Collation of recent data after World War II by the Human Mortality Database allowed analyses, previously unattainable, of modal death-ages for sufficient numbers of selected European cohorts. The aim was to track modes and medians/means ([&ge;]60 years old (y)) of all-cause mortality for both sexes. The only highest-quality, large-number Lexis data available were analyzed: from nine countries: Denmark, Finland, France, Iceland, Italy, Netherlands, Norway, Sweden and Switzerland; raw-data modes (and medians/means [&ge;]60y, plus thin-plate-spline averages), were analyzed, plus pooled data. Here we show that for cohorts 1880-[~]1900 dramatic sex differences existed between death-age changes with all countries except Iceland showing male modal negative trends lasting [~]10-20 years and medians in all countries near-constant or negative lasting [~]10-20 years; whereas females from most countries showed fairly constant positive trends (except Finnish modes and Norwegian medians). For cohorts [~]1900-1919 male and female modal trends were positive (except Dutch and Icelandic cohorts and Finnish females). The net results were that male mortality modes for Danish, Icelandic, Italian, Dutch, Swedish and Norwegian 1919 cohorts were roughly the same as for 1880 cohorts, whereas female death-age modes increased. Results clarify previously knowledge concerning sex differences during this period. Despite improved environment in late adulthood over this period, this did not translate into increased male longevity and earlier events might have sealed their fate, especially in Denmark, Italy, Netherlands, Norway, aand Sweden (and, later, Iceland). RNAi screen of RING/U-box domain ubiquitin ligases identifies critical regulators of tissue regeneration in planarians Regenerative processes depend on the interpretation of signals to coordinate cell behaviors. The role of ubiquitin-mediated signaling is known to be important in many cellular and biological contexts, but its role in regeneration is not well understood. To investigate how ubiquitylation impacts tissue regeneration in vivo, we are studying planarians that are capable of regenerating after nearly any injury using a population of stem cells. Here we used RNAi to screen RING/U-box E3 ubiquitin ligases that are highly expressed in planarian stem cells and stem cell progeny. RNAi screening identified nine genes with functions in regeneration, including the spliceosomal factor prpf19 and histone modifier rnf2; based on their known roles in developmental processes, we further investigated these two genes. We found that prpf19 was required for animal survival but not for stem cell maintenance, suggesting a role in promoting cell differentiation. Because RNF2 is the catalytic subunit of the Polycomb Repressive Complex 1 (PRC1), we also examined other putative members of this complex (CBX and PHC). We observed a striking phenotype of regional tissue misspecification in cbx and phc RNAi planarians. To identify genes regulated by PRC1, we performed RNA-seq after knocking down rnf2 or phc. Although these proteins are predicted to function in the same complex, we found that the set of genes differentially expressed in rnf2 versus phc RNAi were largely non-overlapping. Using in situ hybridization, we showed that rnf2 regulates gene expression levels within a tissue type, whereas phc is necessary for the spatial restriction of gene expression, findings consistent with their respective in vivo phenotypes. This work not only uncovered roles for RING/U-box E3 ligases in stem cell regulation and regeneration, but also identified differential gene targets for two putative PRC1 factors required for maintaining cell-type-specific gene expression in planarians. A rationalized definition of tumor suppressor microRNAsexcludes miR-34a While several microRNAs (miRNAs) have been proposed to act as tumor suppressors, a consensual definition of tumor suppressing miRNAs is still missing. Similarly to coding genes, we propose that tumor suppressor miRNAs must show evidence of genetic or epigenetic inactivation in cancers, and exhibit an anti-tumorigenic (e.g., anti-proliferative) activity under endogenous expression levels. Here we observe that this definition excludes the most extensively studied tumor suppressor candidate miRNA, miR-34a. In analyzable cancer types, miR-34a does not appear to be down-regulated in primary tumors relatively to normal adjacent tissues. Deletion of miR-34a is occasionally found in human cancers, but it does not seem to be driven by an anti-tumorigenic activity of the miRNA, since it is not observed upon smaller, miR-34a-specific alterations. Its anti-proliferative action was observed upon large, supra-physiological transfection of synthetic miR-34a in cultured cells, and our data indicates that endogenous miR-34a levels do not have such an effect. Our results therefore argue against a general tumor suppressive function for miR-34a, providing an explanation to the lack of efficiency of synthetic miR-34a administration against solid tumors. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. LPS-TLR4 pathway exaggerates alcohol-induced acute-on-chronic liver failure via provoking NETs formation Backgrounds: Intrahepatic infiltration of neutrophils is a character of alcoholic acute-on-chronic liver failure (AACLF) and neutrophil extracellular traps (NETs) are an important strategy for neutrophils to fix and kill invading microorganisms. Intestinal bacteria and the gut-liver axis have been thought to play a key role in many liver diseases also including AACLF. However, whether NETs appear in AACLF and play a role in AACLF is still unsure. Methods: WT, NE KO, and TLR4 KO mice were used to build the AACLF model, and the intestinal bacteria were eliminated at the same time and LPS was given. Then the formation of NETs and AACLF related markers were detected. Results: The serum MPO-DNA and LPS concentration was increased in AACLF patients and a correlation was revealed between these two indexes. More intrahepatic NETs formed in AACLF mice by testing MPO-DNA, Cit H3, and NE. These markers decreased with gut detergent and restored markers with gut detergent plus LPS supplement. While NETs formation failed to change with gut microbiome or combine LPS supplement in TLR4 KO mice. As we tested AACLF related characters, liver injury, intrahepatic fat deposition, inflammation, and fibrosis alleviated with depletion of NE. These related marks were also attenuated with gut sterilization by antibiotics and recovered with combined treatment with antibiotics plus LPS. But the liver injury, intrahepatic fat, fibro deposition, and liver inflammation-related markers did show a significant difference in TLR4 KO mice when they received the same treatment. Conclusion: Intestinal-derived LPS promotes NETs formation in AACLF through the TLR4 pathway and further accelerates the AACLF process by NETs. Unifying Transcriptome Regulation and Cellular Operational Efficiency: Application of the Locality and Caching Principles via a Cell-to-Computer Analogy Gene expression is time-consuming, and the delay from transcription activation to produced proteins is sequentially longer from bacteria to yeast and to humans. How human cells overcome the delay and attain operational efficiency, i.e., quick proteomic response to signals, is not well understood. The computer has endured the same system latency issue due to much slower information retrieval (hard drive (HD) to memory and to central processing unit (CPU)) than CPU execution, and mitigated it via efficient memory management, namely, the spatiotemporal locality principles that control specialized user programs and the permanent caching of core system functions - the operating system (OS) kernel. In this study, we unified gene expression and HD-memory-CPU information flow as instances of the Shannon information theory, that support the respective system operations and both consist of three components: information storage, the execution/decoding step, and the channel for the dynamic storage-to-execution information flow; the gene expression machinery and their regulators, and the OS kernel, were deemed as the respective channels. This abstraction prompted a multi-omic comparative analysis, leading to three observations. First, the temporal locality principle explains the mRNA stabilization-by-translation regulatory mechanism and controls specialized cellular functions. Second, the caching principle explains cytoplasmic mRNA sequestration and the defiance of the locality principle by highly sequestered mRNAs. Third, strikingly, in both systems, the caching principle controls the information channels; similar to permanent caching of OS kernel, basic gene expression machinery and their regulators are the top most sequestered mRNAs. In short, the locality and the caching principles act in the cells to control specialized and core cellular functions, respectively. The results close the knowledge gap between, and thus unify, transcriptome regulation and cellular operational latency mitigation. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Inability to switch from ARID1A-BAF to ARID1B-BAF impairs exit from pluripotency and commitment towards neural crest formation in ARID1B-related neurodevelopmental disorders The BAF complex modulates chromatin accessibility. Specific BAF configurations have functional consequences, and subunit switches are essential for cell differentiation. ARID1B and its paralog ARID1A encode for mutually exclusive BAF subunits. De novo ARID1B haploinsufficient mutations cause a neurodevelopmental disorder spectrum, including Coffin-Siris syndrome, which is characterized by neurological and craniofacial features. Here, we reprogrammed ARID1B+/- Coffin-Siris patient-derived skin fibroblasts into iPSCs and modeled cranial neural crest cell (CNCC) formation. We discovered that ARID1B is active only during the first stage of this process, coinciding with neuroectoderm specification, where it is part of a lineage-specific BAF configuration (ARID1B-BAF), which includes SMARCA4 and nine additional subunits. ARID1B-BAF acts as a gatekeeper, ensuring exit from pluripotency and lineage commitment, by attenuating NANOG, SOX2 and thousands of enhancers directly regulated by these two pluripotency factors at the iPSC stage.

In iPSCs, these enhancers are maintained active by an ARID1A-containing BAF. At the onset of differentiation, cells transition from ARID1A-BAF to ARID1B-BAF, eliciting attenuation of the NANOG/SOX2 networks, and triggering pluripotency exit. Coffin-Siris patient cells fail to perform the ARID1A/ARID1B switch and maintain ARID1A-BAF at pluripotency enhancers throughout all stages of CNCC formation. This leads to a persistent and aberrant SOX2 and NANOG activity, which impairs CNCC formation. In fact, despite showing the typical neural crest signature (TFAP2A+, SOX9+), ARID1B-haploinsufficient CNCCs are also NANOG/OCT4-positive, in stark contrast with the ARID1B-wt CNCCs, which are NANOG/OCT4-negative.

These findings suggest a connection between ARID1B mutations, neuroectoderm formation, and a pathogenic mechanism for Coffin-Siris syndrome. Mitofusin 1 and 2 regulation of mitochondrial DNA content is a critical determinant of glucose homeostasis The dynamin-like GTPases Mitofusin 1 and 2 (Mfn1 and Mfn2) are essential for mitochondrial function, which has been principally attributed to their regulation of fission/fusion dynamics. Here, we report that Mfn1 and 2 are critical for glucose-stimulated insulin secretion (GSIS) primarily through control of mitochondrial DNA (mtDNA) content. Whereas Mfn1 and Mfn2 individually were dispensable for glucose homeostasis, combined Mfn1/2 deletion in {beta}-cells reduced mtDNA content, impaired mitochondrial morphology and networking, and decreased respiratory function, ultimately resulting in severe glucose intolerance. Importantly, gene dosage studies unexpectedly revealed that Mfn1/2 control of glucose homeostasis was dependent on maintenance of mtDNA content, rather than mitochondrial structure. Mfn1/2 maintain mtDNA content by regulating the expression of the crucial mitochondrial transcription factor Tfam, as Tfam overexpression ameliorated the reduction in mtDNA content and GSIS in Mfn1/2-deficient {beta}-cells. Thus, the primary physiologic role of Mfn1 and 2 in {beta}-cells is coupled to the preservation of mtDNA content rather than mitochondrial architecture, and Mfn1 and 2 may be promising targets to overcome mitochondrial dysfunction and restore glucose control in diabetes. Targeted isolation of panels of diverse human protective broadly neutralizing antibodies against SARS-like viruses The emergence of current SARS-CoV-2 variants of concern (VOCs) and potential future spillovers of SARS-like coronaviruses into humans pose a major threat to human health and the global economy 1-7. Development of broadly effective coronavirus vaccines that can mitigate these threats is needed 8, 9. Notably, several recent studies have revealed that vaccination of recovered COVID-19 donors results in enhanced nAb responses compared to SARS-CoV-2 infection or vaccination alone 10-13. Here, we utilized a targeted donor selection strategy to isolate a large panel of broadly neutralizing antibodies (bnAbs) to sarbecoviruses from two such donors. Many of the bnAbs are remarkably effective in neutralization against sarbecoviruses that use ACE2 for viral entry and a substantial fraction also show notable binding to non-ACE2-using sarbecoviruses. The bnAbs are equally effective against most SARS-CoV-2 VOCs and many neutralize the Omicron variant. Neutralization breadth is achieved by bnAb binding to epitopes on a relatively conserved face of the receptor binding domain (RBD) as opposed to strain-specific nAbs to the receptor binding site that are commonly elicited in SARS-CoV-2 infection and vaccination 14-18. Consistent with targeting of conserved sites, select RBD bnAbs exhibited in vivo protective efficacy against diverse SARS-like coronaviruses in a prophylaxis challenge model. The generation of a large panel of potent bnAbs provides new opportunities and choices for next-generation antibody prophylactic and therapeutic applications and, importantly, provides a molecular basis for effective design of pan-sarbecovirus vaccines. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. Transcriptional regulation of HSCs in Aging and MDS reveals DDIT3 as a Potential Driver of Dyserythropoiesis Myelodysplastic syndromes (MDS) are hematopoietic stem cell (HSC) malignancies characterized by ineffective hematopoiesis, with increased incidence in elderly individuals. In this work, we analyzed the transcriptome of human HSCs purified from young and elderly healthy donors, as well as MDS patients, identifying transcriptional alterations following eight different patterns of expression. While aging-associated lesions seemed to predispose HSCs to myeloid transformation, disease-specific alterations may trigger MDS development. Among MDS-specific lesions, we detected the upregulation of the transcription factor DDIT3. Overexpression of DDIT3 in human healthy HSCs induced an MDS-like transcriptional state, and a delay in erythropoiesis. Such effect was associated with downregulation of transcription factors required for normal erythropoiesis, and with a failure in the activation of their transcriptional programs. Moreover, DDIT3 knockdown in CD34+ cells from MDS patients with anemia was able to restore erythropoiesis. These results identify DDIT3 as a driver of dyserythropoiesis, and a potential therapeutic target to restore the inefficient erythropoiesis characterizing MDS patients.

STATEMENT OF SIGNIFICANCEThis study defines how human aging and MDS development are characterized by transcriptional alterations in HSCs that follow different patterns, some of which may contribute to myeloid transformation. Among them, we demonstrate how MDS-specific upregulation of DDIT3 in HSCs induces dyserythropoiesis, while its knockdown in HSPCs from MDS patients restores proper erythroid differentiation. Bayesian Additive Regression Trees for Genotype by Environment Interaction Models We propose a new class of models for the estimation of genotype by environment (GxE) interactions in plant-based genetics. Our approach, named AMBARTI, uses semi-parametric Bayesian additive regression trees to accurately capture marginal genotypic and environment effects along with their interaction in a fully Bayesian model. We demonstrate that our approach is competitive or superior to similar models widely used in the literature via both simulation and a real world data set. Furthermore, we introduce new types of visualisation to properly assess both the marginal and interactive predictions from the model. An R package that implements our approach is available at https://github.com/ebprado/ambarti. Collective polymerase dynamics emerge from DNA supercoiling during transcription All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts and transcription stops. DNA supercoiling can also mediate non-local interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton) whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.

PACS numbers: Loci for prediction of penicillin and tetracycline susceptibility in Neisseria gonorrhoeae: a genome wide association study BackgroundWhile Neisseria gonorrhoeae poses an urgent public health threat because of increasing antimicrobial resistance, much of the circulating population remains susceptible to historical treatment regimens. Point-of-care diagnostics that report susceptibility could allow for reintroduction of these regimens, but development of such diagnostics has been limited to ciprofloxacin, for which susceptibility can be predicted from a single locus.

MethodsWe assembled a dataset of 12,045 N. gonorrhoeae genomes with phenotypic resistance data for tetracycline (n = 3,611) and penicillin (n = 6,935). Using conditional genome wide association studies (GWAS), we sought to define genetic variants associated with susceptibility to penicillin and tetracycline. We evaluated the sensitivity and specificity of these variants for predicting susceptibility and non-resistance in our collection of gonococcal genomes.

FindingsIn our conditional penicillin GWAS, the presence of a genetic variant defined by a non-mosaic penA allele without an insertion at codon 345 was significantly associated with penicillin susceptibility and had the highest negative effect size of significant variants (p = 5.0 x 10-14, {beta} = -2.5). In combination with the absence of blaTEM, this variant predicted penicillin susceptibility with high specificity (99.8%) and modest sensitivity (36.7%). For tetracycline, the wild type allele at rpsJ codon 57, encoding valine, was significantly associated with tetracycline susceptibility (p = 5.6 x 10-16, {beta} = -1.6) after conditioning on the presence of tetM. The combination of rpsJ codon 57 allele and tetM absence predicted tetracycline susceptibility with high specificity (97.2%) and sensitivity (88.7%).

InterpretationAs few as two genetic loci can predict susceptibility to penicillin and tetracycline in N. gonorrhoeae with high specificity. Molecular point-of-care diagnostics targeting these loci have the potential to increase available treatments for gonorrhea.

FundingNational Institute of Allergy and Infectious Diseases, the National Science Foundation, and the Smith Family Foundation

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed with the terms "Neisseria gonorrhoeae" and "diagnostic" or "assay" plus "penicillin" or "tetracycline" for reports in any language published up to July 1, 2021. We additionally searched for "Neisseria gonorrhoeae" and "genome wide association study". We found that previously proposed molecular diagnostics for penicillin and tetracycline susceptibility either exclusively focused on plasmid-mediated resistance (i.e., targeting blaTEM or tetM) or did not include variants in genes encoding antibiotic targets (e.g., did not include penA or rpsJ). Targets for molecular surveillance have focused on resistance-associated alleles rather than susceptibility-associated alleles. We did not find any previous penicillin or tetracycline conditional genome wide association studies (GWAS) in N. gonorrhoeae.

Added value of this studyTo identify targets for molecular diagnostics that predict penicillin and tetracycline susceptibility, we conducted GWAS conditioning on the presence of plasmid-mediated resistance determinants to detect chromosomal loci with the highest association with susceptibility. We discovered a sequence (penA_01) that differentiates susceptible isolates from those with a resistance-associated insertion at codon 345 and from those with mosaic penA alleles, which is associated with penicillin susceptibility. We also found that rpsJ codon 57 was the chromosomal locus contributing the most to tetracycline susceptibility. The combination of these chromosomal loci and the absence of plasmid encoded determinants predicts penicillin and tetracycline susceptibility with high specificity in both a large global collection of N. gonorrhoeae and a validation dataset consisting of recently published genomes from CDCs Gonococcal Isolate Surveillance Program (GISP) surveillance collected in 2018.

Implications of all the available evidenceThe chromosomal loci penA_01 and rpsJ codon 57 in combination with plasmid loci blaTEM and tetM are candidates for the development of point-of-care molecular diagnostics for penicillin and tetracycline susceptibility. The loci may be combined with the currently available ciprofloxacin susceptibility diagnostics to predict susceptibility to multiple antibiotics. Additionally, our study suggests that conditional GWAS focused on variants associated with susceptibility may be a promising approach to identify minimal sets of loci for molecular diagnostics and surveillance. DNA-controlled Spatio-temporal Patterning of a Cytoskeletal Active Gel Living cells move and change their shape because signaling chemical reactions modify the state of their cytoskeleton; an active gel that converts chemical energy into mechanical forces. To create life-like materials, it is thus key to engineer chemical pathways that drive active gels. Here, we describe the preparation of DNA-responsive surfaces that control the activity of a cytoskeletal active gel comprised of microtubules: a DNA signal triggers the release of molecular motors from the surface into the gel bulk, generating forces that structure the gel. Depending on the DNA sequence and concentration, the gel forms a periodic band pattern or contracts globally. Finally, we show that the structuration of the active gel can be spatially controlled in the presence of a gradient of DNA concentration. We anticipate that such DNA-controlled active matter will contribute to the development of life-like materials with self-shaping properties. Arabidopsis AtPME2 has a pH-dependent processivity and control cell wall mechanical properties Pectin methylesterases (PMEs) modify homogalacturonans (HG) chemistry and play a key role in regulating primary cell wall mechanical properties. How PME activity can fine-tune pectin structure in the growing plant has remained elusive. Here we report on the Arabidopsis AtPME2, which we found to be highly expressed during lateral root emergence and dark-grown hypocotyl elongation. We produced the mature active enzyme using heterologous expression in Pichia pastoris and characterized it through the use of a generic plant PME antiserum suitable for detecting recombinant and native enzyme independent of species source. At neutral pH AtPME2 is preferentially active on pectins with a degree of 55-70% methylesterification and can be inhibited by PME inhibitor protein (PMEI). We show that the mode of action for AtPME2 can switch from full processivity (at pH 8), creating large blocks of unmethylated galacturonic acid, to low processivity (at pH 5) and relate these observations to the differences in electrostatic potential of the protein at acidic and alkaline pH. To assess the role of AtPME2 in development, we characterized two knock-out lines. We show that in the context of acidified apoplast, low-processive demethylesterification by AtPME2 can loosen the cell wall, with consequent increase in cell elongation and etiolated hypocotyl length. Our study brings insights into how the pH-dependent regulation by PME activity could affect pectin structure and associated cell wall mechanical properties in expansion.

One sentence summaryThe processivity of AtPME2, a pectin methylesterase that fine-tunes cell wall pectins is modulated by pH in vitro and impacts the mechanical properties of the wall, affecting development in planta. Social integration predicts survival in female white-faced capuchin monkeys Across multiple species of social mammals, a growing number of studies have found that that individual sociality is associated with survival. In long-lived species, like primates, lifespan is one of the main components of fitness. We used 18 years of data from the Lomas Barbudal Monkey project to quantify social integration in 11 capuchin (Cebus capucinus) groups and tested whether female survivorship was associated with females tendencies to interact with three types of partners: (1) all group members, (2) adult females, (3) adult males. We found strong evidence that females who engaged more with other females in affiliative interactions and foraged in close proximity experienced increased survivorship. We found some weak evidence that females might also benefit from engaging in more support in agonistic contexts with other females. These benefits were evident in models that account for the females rank and group size. Female interactions with all group members also increased survival, but the estimates of the effects were more uncertain. In interactions with adult males, only females who provided more grooming to males survived longer. The results presented here suggest that social integration may result in survival-related benefits. Females might enjoy these benefits through exchanging grooming for other currencies, such as coalitionary support or tolerance. Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae The endosymbiotic origin of plastids from cyanobacteria gave eukaryotes photosynthetic capabilities and launched the diversification of countless forms of algae. These primary plastids are found in members of the eukaryotic supergroup Archaeplastida. All known archaeplastids still retain some form of primary plastids, which are widely assumed to have a single origin. Here, we used single-cell genomics from natural samples combined with phylogenomics to infer the evolutionary origin of the phylum Picozoa, a globally distributed but seemingly rare group of marine microbial heterotrophic eukaryotes. Strikingly, the analysis of 43 single-cell genomes shows that Picozoa belong to Archaeplastida, specifically related to red algae and the phagotrophic rhodelphids. These picozoan genomes support the hypothesis that Picozoa lack a plastid, and further reveal no evidence of an early cryptic endosymbiosis with cyanobacteria. These findings change our understanding of plastid evolution as they either represent the first complete plastid loss in a free-living taxon, or indicate that red algae and rhodelphids obtained their plastids independently of other archaeplastids. 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. Epitope-tagging of the endogenous murine BiP/GRP78/Hspa5 locus allows direct analysis of the BiP interactome and protein misfolding in vivo BiP/GRP78, encoded by the Hspa5 gene, is the major HSP70 family member in the endoplasmic reticulum (ER) lumen, and controls ER protein folding. The essential functions of BiP in facilitating proper protein folding are mainly mediated through its dynamic interaction with unfolded or misfolded client proteins, and by serving as a negative regulator of the Unfolded Protein Response. A mechanistic understanding of the dynamics of BiP interaction with its protein partners is essential to understand ER biology, and therefore, we have sought to develop a tractable model to study misfolded protein interaction with BiP. For this purpose, we have used homologous recombination to insert a 3xFLAG epitope tag into the endogenous murine Hspa5 gene, just upstream from the essential KDEL signal necessary for ER localization of BiP. Tagging BiP in this way did not alter Hspa5 expression under basal or ER-stress induced conditions in hepatocytes ex vivo or in fibroblasts. Furthermore, the tag did not alter the cellular localization of BiP or its functionality. All of these findings in primary tissue culture were also confirmed in vivo in livers of heterozygous mice harboring one WT and one FLAG-tagged Hspa5 allele. Hepatocyte-specific BiP-FLAG modification did not alter liver function or UPR signaling. Importantly, immunoprecipitation with anti-FLAG antibody completely pulled down FLAG-tagged BiP from lysates of BiP-FLAG expressing livers. Affinity purification-mass spectrometry (AP-MS) of BiP-FLAG protein complexes isolated from the BiP-FLAG-expressing livers of tunicamycin (TM)- and vehicle-treated mice revealed a marked increase in interaction of glycoproteins with BiP-FLAG in response to inhibition of N-glycosylation due to TM-treatment, validating utility of our BiP-Flag mice as a tool to identify ER misfolded proteins in vivo. Significantly, our AP-MS analysis also provided in vivo evidence demonstrating that BiP-FLAG binds to UPR transducers IRE1 and PERK under basal conditions but is released upon TM-treatment to activate UPR. We have also employed this mouse model to demonstrate that proinsulin in pancreatic {beta} cells misfolds and interacts with BiP-FLAG in healthy mice. In summary, we generated a novel model that can be used to investigate the BiP interactome in vivo under physiological and pathophysiological conditions in a cell type-specific manner. This tool provides, for the first time, an unbiased approach to identify unfolded and misfolded BiP-client proteins, and a new approach to study ER protein misfolding in a cell-type and temporal manner to uncover the role of BiP in many essential ER processes. Hippocampal sharp wave-ripples and the associated sequence replay emerge from structured synaptic interactions in a network model of area CA3 Hippocampal place cells are activated sequentially as an animal explores its environment. These activity sequences are internally recreated ("replayed"), either in the same or reversed order, during bursts of activity (sharp wave-ripples; SWRs) that occur in sleep and awake rest. SWR-associated replay is thought to be critical for the creation and maintenance of long-term memory. In order to identify the cellular and network mechanisms of SWRs and replay, we constructed and simulated a data-driven model of area CA3 of the hippocampus. Our results show that the chain-like structure of recurrent excitatory interactions established during learning not only determines the content of replay, but is essential for the generation of the SWRs as well. We find that bidirectional replay requires the interplay of the experimentally confirmed, temporally symmetric plasticity rule, and cellular adaptation. Our model provides a unifying framework for diverse phenomena involving hippocampal plasticity, representations, and dynamics, and suggests that the structured neural codes induced by learning may have greater influence over cortical network states than previously appreciated. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important.