Coronavirus-associated molecular mimicry through homology to a SARS-CoV-2 peptide could be leading to susceptibility in patients with HLA-A*02:01 and HLA-A*24:02 serotypes This study aims to predict autoimmunity-related pathological mechanisms that possess risk for individuals with specific human leukocyte antigen (HLA) serotypes and shared by certain coronaviruses including SARS-CoV-2, based on homology to a SARS-CoV-2 peptide. With the given aim, 1-) coronavirus-associated sequences, which are homologous to the 15mer SARS-CoV-2 peptide CFLGYFCTCYFGLFC, are obtained. 2-) Human peptides that have at least 7 residue matches with those coronavirus sequences, and the SARS-CoV-2 15mer, are found. 3-) Epitope pairs, which are sourced by those aligned coronavirus and human sequences are identified. 4-) Epitope pairs that are predicted to bind strongly not only to the same HLA allele with each other but also to the same HLA allele as those of the respective alignment of the SARS-CoV-2 peptide are selected. Following are the identified proteins or peptides (with HLA-A*02:01 or HLA-A*24:02 epitopes), as described in 1-to-4: Immunoglobulin heavy chain junction regions, CRB1 isoform I precursor, slit homolog 2 protein, hCG1995581, hCG2028737, phospholipid phosphatase-related protein type 2. Among those, CRB1 isoform I precursor sequence with the predicted HLA-A*24:02 epitope aligns with the highest number of different sequences. Results imply autoimmunity risk in COVID-19 patients with HLA-A*02:01 and HLA-A*24:02 serotypes, through molecular mimicry, as a shared pathogenicity risk that can be prevalent upon getting infected with certain coronaviruses. These can pave way to improved risk groups assessment and autoimmunity treatment options, for COVID-19 and its associated diseases. Also, the approach in this study can be used to predict prospective pathologies of the transmissible variants in susceptible humans.
An Immune Cell Atlas Reveals Dynamic COVID-19 Specific Neutrophil Programming Amenable to Dexamethasone Therapy SARS-CoV-2 is a novel coronavirus that causes acute respiratory distress syndrome (ARDS), death and long-term sequelae. Innate immune cells are critical for host defense but are also the primary drivers of ARDS. The relationships between innate cellular responses in ARDS resulting from COVID-19 compared to other causes of ARDS, such as bacterial sepsis is unclear. Moreover, the beneficial effects of dexamethasone therapy during severe COVID-19 remain speculative, but understanding the mechanistic effects could improve evidence-based therapeutic interventions. To interrogate these relationships, we developed an scRNA-Seq and plasma proteomics atlas (biernaskielab.ca/COVID_neutrophil). We discovered that compared to bacterial ARDS, COVID-19 was associated with distinct neutrophil polarization characterized by either interferon (IFN) or prostaglandin (PG) active states. Neutrophils from bacterial ARDS had higher expression of antibacterial molecules such as PLAC8 and CD83. Dexamethasone therapy in COVID patients rapidly altered the IFNactive state, downregulated interferon responsive genes, and activated IL1R2+ve neutrophils. Dexamethasone also induced the emergence of immature neutrophils expressing immunosuppressive molecules ARG1 and ANXA1, which were not present in healthy controls. Moreover, dexamethasone remodeled global cellular interactions by changing neutrophils from information receivers into information providers. Importantly, male patients had higher proportions of IFNactive neutrophils, a greater degree of steroid-induced immature neutrophil expansion, and increased mortality benefit compared to females in the dexamethasone era. Indeed, the highest proportion of IFNactive neutrophils was associated with mortality. These results define neutrophil states unique to COVID-19 when contextualized to other life-threatening infections, thereby enhancing the relevance of our findings at the bedside. Furthermore, the molecular benefits of dexamethasone therapy are also defined, and the identified pathways and plasma proteins can now be targeted to develop improved therapeutics.
A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea The importance of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae, Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including: (i) phylogenomics, (ii) the distribution of shared derived alleles, and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics, D-foil, and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggests that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species gene pools. Disc and Actin-Associated Protein 1 Influence Attachment in the Intestinal Parasite Giardia lamblia The deep-branching eukaryote Giardia lamblia is an extracellular parasite that attaches to the host intestine via a microtubule-based structure called the ventral disc. Control of attachment is mediated in part by the movement of two regions of the ventral disc that either permit or exclude the passage of fluid under the disc. Several known disc-associated proteins (DAPs) contribute to disc structure and function, but no force-generating protein has been identified among them. We recently identified several Giardia actin (GlActin) interacting proteins at the ventral disc, which could potentially employ actin polymerization for force generation and disc conformational changes. One of these proteins, Disc and Actin Associated Protein 1 (DAAP1), is highly enriched at the two regions of the disc previously shown to be important for fluid flow during attachment. In this study, we investigate the role of both GlActin and DAAP1 in ventral disc morphology and function. We confirmed interaction between GlActin and DAAP1 through coimmunoprecipitation, and used immunofluorescence to localize both proteins throughout the cell cycle and during trophozoite attachment. Similar to other DAPs, the association of DAAP1 with the disc is stable, except during cell division when the disc disassembles. Depletion of GlActin by translation-blocking antisense morpholinos resulted in both impaired attachment and defects in the ventral disc, indicating that GlActin contributes to disc-mediated attachment. Depletion of DAAP1 through CRISPR interference resulted in intact discs but impaired attachment, gating, and flow under the disc. As attachment is essential for infection, elucidation of these and other molecular mediators is a promising area for development of new therapeutics against a ubiquitous parasite.

Author SummaryGiardia lamblia is a single-celled organism and one of the most common gastrointestinal parasites worldwide. In developing countries, recurrent Giardia infections are common, due to lack of access to clean water. Giardia infections can lead to diarrhea, vomiting, dehydration, disruption of the intestinal microbiome, and chronic infections can lead to irritable bowel syndrome. Because existing drug treatments have side effects and Giardias resistance to drugs is increasing, new treatment strategies are needed. The parasites attachment to the hosts intestine is mediated by a Giardia-specific structure that resembles a suction cup and is called the ventral adhesive disc. We previously identified DAAP1, a protein which interacts with Giardia actin and localizes to the ventral disc. Here, we explore the relationship between these two proteins and investigate their role in disc-based attachment. Most disc proteins, including DAAP1, are unrelated to any human proteins, making them appealing drug targets to inhibit parasite attachment and infection. Snrpb, the gene mutated in CCMS, is haploinsufficient and required in neural crest cells for proper splicing of developmentally important genes essential for craniofacial morphogenesis Heterozygous mutations in SNRPB, an essential core component of the five small ribonucleoprotein particles of the spliceosome, are responsible for Cerebrocostomandibular Syndrome (CCMS). We show that Snrpb heterozygous embryos arrest shortly after implantation. Additionally, heterozygous deletion of Snrpb in the developing brain and neural crest cells models craniofacial malformations found in CCMS, and results in death shortly after birth. RNAseq analysis of mutant heads prior to morphological defects revealed increased exon-skipping and intron-retention in association with increased 5 splice site strength. We found increased exon-skipping in negative regulators of the P53-pathway, increased levels of nuclear P53, P53-target genes. However, removing TrP53 in Snrpb heterozygous mutant neural crest cells did not completely rescue craniofacial development. We also found a small but significant increase in exon-skipping of several transcripts required for head and midface development, including Smad2 and Rere. Furthermore, mutant embryos exhibited ectopic or missing expression of Fgf8 and Shh, which are required to coordinate face and brain development. Thus, we propose that mis-splicing of transcripts that regulate P53-activity and craniofacial-specific genes both contribute to craniofacial malformations. A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. DNA supercoiling-mediated collective behavior of co-transcribing RNA polymerases Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or "push" between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate via the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the torsional response of supercoiled DNA. Using stochastic simulations, we demonstrate long-range cooperation between co-transcribing RNAPs mediated by DNA supercoiling. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern a key cellular process. Inhibitory proteins block substrate access by occupying the active site cleft of Bacillus subtilis intramembrane metalloprotease SpoIVFB Intramembrane proteases function in numerous signaling pathways that impact health, but knowledge about regulation of intramembrane proteolysis is limited. We examined inhibition of intramembrane metalloprotease SpoIVFB by proteins BofA and SpoIVFA. We found that BofA residues in and near a predicted transmembrane segment are required for SpoIVFB inhibition, and cross-linking experiments indicated that this transmembrane segment occupies the SpoIVFB active site cleft. SpoIVFA is also required for SpoIVFB inhibition. The inhibitory proteins block access of the substrate N-terminal Proregion to the membrane-embedded SpoIVFB active site, based on additional cross-linking experiments; however, the inhibitory proteins did not prevent interaction between the substrate C-terminal region and the SpoIVFB soluble domain. A structural model was built of SpoIVFB in complex with BofA and parts of SpoIVFA and substrate, using partial homology and constraints from cross-linking and co-evolutionary analyses. The model predicts that conserved BofA residues interact to stabilize a transmembrane segment and a membrane-embedded C-terminal region. SpoIVFA is predicted to bridge the BofA C-terminal region and SpoIVFB, forming a membrane-embedded inhibition complex. Our results reveal a novel mechanism of intramembrane protease inhibition with clear implications for relief from inhibition in vivo and design of inhibitors as potential therapeutics. A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice Flexible reprogramming of Pristionchus pacificus motivation for attacking Caenorhabditis elegans in predator-prey competition Animals with diverse diets must adapt their food priorities to a wide variety of environmental conditions. This diet optimization problem is especially complex for predators that compete with prey for food. Although predator-prey competition is widespread and ecologically critical, it remains difficult to disentangle predatory and competitive motivations for attacking competing prey. Here, we dissect the foraging decisions of the omnivorous nematode Pristionchus pacificus to reveal that its seemingly failed predatory attempts against Caenorhabditis elegans are actually motivated acts of efficacious territorial aggression. While P. pacificus easily kills and eats larval C. elegans with a single bite, adult C. elegans typically survives and escapes from bites. However, nonfatal biting can provide competitive benefits by reducing access of adult C. elegans and its progeny to bacterial food that P. pacificus also eats. We show that P. pacificus considers costs and benefits of both predatory and territorial outcomes to decide which food goal, prey or bacteria, should guide its motivation for biting. These predatory and territorial motivations impose different sets of rules for adjusting willingness to bite in response to changes in bacterial abundance. In addition to biting, predatory and territorial motivations also influence which search tactic P. pacificus uses to increase encounters with C. elegans. When treated with an octopamine receptor antagonist, P. pacificus switches from territorial to predatory motivation for both biting and search. Overall, we demonstrate that P. pacificus assesses alternate outcomes of attacking C. elegans and flexibly reprograms its foraging strategy to prioritize either prey or bacterial food. 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. Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. Discovery and functional assessment of a novel adipocyte population driven by intracellular Wnt/β-catenin signaling in mammals Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. Allopurinol Blocks The Formation And Progression Of Aortic Aneurysm In A Mouse Model Of Marfan Syndrome BackgroundThe pathogenesis and progression of aortic aneurysm in Marfan syndrome (MFS) involves dysregulated TGF-{beta} and nitric oxide signaling, altered hemodynamics, and biomechanical forces. Increasing evidence indicates that redox stress participates in MFS aortopathy, though its contribution is not well established. We reported elevated reactive oxygen species (ROS) formation and NADPH oxidase NOX4 upregulation in MFS patients and mice aortic samples. Here we address the contribution of xanthine oxidoreductase (XOR), which is transcribed as dehydrogenase (XDH) and then post-translationally modified to oxidase (XO), catabolizing mainly purines into uric acid and ROS.

Methods and ResultsWe report that in MFS patients, XOR protein levels increased in aortic samples. In MFS mice (Fbn1C1041G/+), XOR mRNA transcripts augmented in the aorta of young but not older MFS mice. Moreover, the enzymatic activity of XO form increased with respect to XDH. The administration of the XOR inhibitor allopurinol (ALO) blocked the progression of aortic root aneurysm in MFS mice. ALO was also protective when administrated before the appearance of aneurysm. MFS-induced elastic fiber fragmentation, fibrotic remodeling, nuclear translocation of pNRF2, and increased 3-nitrotyrosine levels in the aortic tunica media as well as endothelial dysfunction were all prevented by allopurinol. The large overproduction of H2O2 in MFS aorta did not occur in allopurinol-treated mice whose uric acid plasma levels remained unaltered.

ConclusionsThis study strengthens the concept that redox stress is an important determinant of aortic aneurysm formation and progression in MFS and supports a clinical trial for ALO in the pharmacological treatment of MFS aortopathy.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/464182v3_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@eab02aorg.highwire.dtl.DTLVardef@1d3b39eorg.highwire.dtl.DTLVardef@1b54d0forg.highwire.dtl.DTLVardef@1d61ad0_HPS_FORMAT_FIGEXP  M_FIG C_FIG Hierachical Resampling for Bagging in Multi-Study Prediction with Applications to Human Neurochemical Sensing We propose the "study strap ensemble," which combines advantages of two common approaches to fitting prediction models when multiple training datasets ("studies") are available: pooling studies and fitting one model versus averaging predictions from multiple models each fit to individual studies. The study strap ensemble fits models to bootstrapped datasets, or "pseudo-studies." These are generated by resampling from multiple studies with a hierarchical resampling scheme that generalizes the randomized cluster bootstrap. The study strap is controlled by a tuning parameter that determines the proportion of observations to draw from each study. When the parameter is set to its lowest value, each pseudo-study is resampled from only a single study. When it is high, the study strap ignores the multi-study structure and generates pseudo-studies by merging the datasets and drawing observations like a standard bootstrap. We empirically show the optimal tuning value often lies in between, and prove that special cases of the study strap draw the merged dataset and the set of original studies as pseudo-studies. We extend the study strap approach with an ensemble weighting scheme that utilizes information in the distribution of the covariates of the test dataset.

Our work is motivated by neuroscience experiments using real-time neurochemical sensing during awake behavior in humans. Current techniques to perform this kind of research require measurements from an electrode placed in the brain during awake neurosurgery and rely on prediction models to estimate neurotransmitter concentrations from the electrical measurements recorded by the electrode. These models are trained by combining multiple datasets that are collected in vitro under heterogeneous conditions in order to promote accuracy of the models when applied to data collected in the brain. A prevailing challenge is deciding how to combine studies or ensemble models trained on different studies to enhance model generalizability.

Our methods produce marked improvements in simulations and in this application. All methods are available in the studyStrap CRAN package. 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. Suppression of sost/sclerostin and dickkopf-1 promote intervertebral disc structure in mice Intervertebral disc (IVD) degeneration is a leading cause of low back pain and characterized by accelerated extracellular matrix breakdown and IVD height loss but there is no approved pharmacological therapeutic. Deletion of Wnt signaling receptor Lrp5 induces IVD degeneration and suggests that Wnt signaling in the IVD may be responsive to inhibition of Wnt signaling inhibitors sost(gene)/sclerostin(protein) or dickkopf-1 (dkk1). Anti-sclerostin antibody (Scl-Ab) is an FDA-approved bone therapeutic that activates Wnt signaling. We (1) determined if pharmacological neutralization of sclerostin, dkk1 or their combination stimulate Wnt signaling and promote IVD structure and (2) determined the extent of the response of the IVD to global, persistent deletion of sost. Nine-week-old C57Bl/6J female mice (n=6-7/grp) were subcutaneously injected 2x/wk for 5.5 wk with scl-Ab (25 mg/kg), dkk1-Ab (25 mg/kg), 3:1 scl-Ab/dkk1-Ab (18.75:6.25 mg/kg) or vehicle (Veh). Separately, IVD of sost KO and WT (wildtype) mice (n=8, grp) were harvested at 16 weeks of age. First, compared to vehicle, scl-Ab, dkk1-Ab and 3:1 scl-Ab/dkk1-Ab similarly increased lumbar IVD height and {beta}-catenin gene expression. Despite these similarities, scl-Ab decreased cellular stress-related heat shock protein gene expressions while neither dkk1-Ab nor scl-Ab/dkk1-Ab altered the same. Genetically and compared to WT, sost KO increased MRI-determined hydration and proteoglycan staining in the IVD. Notably, persistent deletion of sost was compensated by upregulation of dkk1, which consequently reduced the cell nuclear fraction for Wnt signaling transcription factor {beta}-catenin in whole IVD. Lastly, RNA-sequencing pathway analysis confirmed the compensatory suppression of Wnt signaling and determined a reduction of cellular stress pathways. Together, suppression of sost/sclerostin or dkk1 each promote IVD structure by stimulating Wnt signaling, but sclerostin and dkk1 may differentially regulate cellular stress pathways. Ultimately, postmenopausal women prescribed scl-Ab injections to prevent vertebral fracture may also benefit from a restoration of IVD height and health.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=72 SRC="FIGDIR/small/449486v2_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@db2e17org.highwire.dtl.DTLVardef@1af1bbborg.highwire.dtl.DTLVardef@145cc4corg.highwire.dtl.DTLVardef@167a19e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Suppression of Wnt signaling inhibitors by genetic or pharmacological approaches promotes intervertebral disc structure and hydration by Wnt signaling. However, persistent activation of Wnt signaling induces a compensatory reduction of Wnt signaling that shifts IVD cells toward a chondrocyte-like (CLC) phenotype. AF: annulus fibrosus, NC: notochordal cell, NP: nucleus pulposus, PG: proteoglycan 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. First Report of Lactoria cornuta from Karnataka, West coast of India, Role of tropical cyclones in Range Extension India is known for its rich biodiversity and is fortunate to have several endemic species from different classes of vertebrates. India is home to 7.5% of the global fish diversity, with 91 endemic species of ray finned fishes, the actinopterygians. Some fish species have never been reported until recently, one such example is the long-horned Cowfish (Ostracidae), best known for inhabiting only certain areas of the Indo-Pacific oceans. However, this fish has been reported from Bay of Bengal, Costs of Goa and Lakshadweep. This is the first time this fish has been found on the Coasts of Karnataka, west coast of India. Given its morphology, the migration seems highly improbable since it can only endure lethargic swimming. In this article, we have explored the role of other external forces that could have contributed to its range extension in the east and west coast of India. In this paper we are exploring the role of Super Cyclones as a vital force in determining the expansion of its range. The fact of collection longhorn cowfish for the first time from the coast of Karwar coincides with the passage of Cyclone Ockhi on the west coast of India this suggests the possible role of super cyclone in its range extension. A year later, another specimen was recovered from a location very close to the first, which indicates the successful establishment of Lacturia cornuta in its new environment. Related morphometric and meristic studies of our specimens are consistent with previous studies reported in the Bay of Bengal. Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems. Hierarchical Bayesian Integrated Model for Estimating Migratory Bird Harvest in Canada The Canadian Wildlife Service (CWS) requires reliable estimates of the harvest of migratory game birds, including waterfowl, to effectively manage populations of these hunted species. The National Harvest Survey is an annual survey of hunters who purchase Canadas mandatory migratory game bird hunting permit, integrating information from a survey of hunting activity with information from a separate survey of species composition in the harvest. We use these survey data to estimate the number of birds harvested for each species, as well as hunting activity metrics such as the number of active hunters and days spent hunting. The analytical methods used to generate these estimates have not changed since the survey was first designed in the early 1970s. Here we describe a new hierarchical Bayesian integrated model, which replaces the series of ratio estimators that comprised the old model. We are now using this new model to generate estimates for migratory bird harvests as of the 2019-2020 hunting season, and to generate updated estimates for all earlier years. The hierarchical Bayesian model uses over-dispersed Poisson distributions to model mean hunter activity and harvest (zero inflated Poisson and zero truncated Poisson, respectively). It also includes multinomial distributions to model some key components including, variation in total harvest across periods of the hunting season, the species composition of the harvest within each of those periods, and the age and sex composition in the harvests of a given species. We estimated the parameters of the Poisson and the multinomial distributions for each year as random effects using first-difference time-series. This time-series component allows the model to share information across years and reduces the sensitivity of the estimates to annual sampling noise. The new model estimates are generally very similar to those from the old model, particularly for the species that occur most commonly in the harvest, and so the results do not suggest any major changes to harvest management decisions and regulations. However, estimates for all species from the new model are more precise and less susceptible to annual sampling error, particularly for species that occur less commonly in the harvest (e.g., sea ducks and other species of conservation concern). This new model, with its hierarchical Bayesian framework, will also facilitate future improvements and elaborations, allowing the incorporation of prior information from the rich literature and knowledge in game bird management and biology. Construction of a new chromosome-scale, long-read reference genome assembly for the Syrian hamster, Mesocricetus auratus BackgroundThe Syrian hamster (Mesocricetus auratus) has been suggested as a useful mammalian model for a variety of diseases and infections, including infection with respiratory viruses such as SARS-CoV-2. The MesAur1.0 genome assembly was generated in 2013 using whole-genome shotgun sequencing with short-read sequence data. Current more advanced sequencing technologies and assembly methods now permit the generation of near-complete genome assemblies with higher quality and greater continuity.

FindingsHere, we report an improved assembly of the M. auratus genome (BCM_Maur_2.0) using Oxford Nanopore Technologies long-read sequencing to produce a chromosome-scale assembly. The total length of the new assembly is 2.46 Gbp, similar to the 2.50 Gbp length of a previous assembly of this genome, MesAur1.0. BCM_Maur_2.0 exhibits significantly improved continuity with a scaffold N50 that is 6.7 times greater than MesAur1.0. Furthermore, 21,616 protein coding genes and 10,459 noncoding genes are annotated in BCM_Maur_2.0 compared to 20,495 protein coding genes and 4,168 noncoding genes in MesAur1.0. This new assembly also improves the unresolved regions as measured by nucleotide ambiguities, where approximately 17.11% of bases in MesAur1.0 were unresolved compared to BCM_Maur_2.0 in which the number of unresolved bases is reduced to 3.00%.

ConclusionsAccess to a more complete reference genome with improved accuracy and continuity will facilitate more detailed, comprehensive, and meaningful research results for a wide variety of future studies using Syrian hamsters as models. Synergy of chemotherapy and macrophage depletion leads to B cell-mediated T cell memory activation and durable triple negative breast cancer regression Immunosuppressive elements within the tumor microenvironment such as Tumor Associated Macrophages (TAMs) can present a barrier to successful anti-tumor responses by cytolytic T cells. We employed preclinical syngeneic p53 null mouse models of triple negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was used to successfully treat several highly aggressive TNBC murine mammary tumors and lung metastasis. Using this regimen and single cell RNA sequencing we characterized tumor infiltrating lymphocytes (TILs) including helper T cells and antigen-presenting B cells that were highly enriched in good responders to combination therapy. Using high dimensional imaging techniques, we identified the close spatial localization of B220+ CD86+ activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks post-treatment in one model that also exhibited long-term tumor regression post-treatment. We also characterized the transcriptional and metabolic heterogeneity of TAMS in these two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model is highly expressed and conserved in human claudin-low breast cancers, and high expression of the T12 signature correlated with reduced overall survival. This T12 tumor TAM signature may help identify human claudin-low breast cancer patients that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rationale combinations of immunotherapy.

SignificanceA treatment regimen that harnessed the immunostimulatory effects of cyclophosphamide coupled with the inhibition of CSF1R was used to successfully treat several highly aggressive claudin-low TNBC murine mammary tumors and lung metastasis. sgcocaller and comapr: personalised haplotype assembly and comparative crossover map analysis using single-gamete sequencing data Profiling gametes of an individual enables the construction of personalised haplotypes and meiotic crossover landscapes, now achievable at larger scale than ever through the availability of high-throughput single-cell sequencing technologies. However, high-throughput single-gamete data commonly have low depth of coverage per gamete, which challenges existing gametebased haplotype phasing methods. In addition, haplotyping a large number of single gametes from high-throughput singlecell DNA sequencing data and constructing meiotic crossover profiles using existing methods requires intensive processing. Here, we introduce efficient software tools for the essential tasks of generating personalised haplotypes and calling crossovers in gametes from single-gamete DNA sequencing data (sgcocaller), and constructing, visualising, and comparing individualised crossover landscapes from single gametes (comapr). With additional data pre-possessing, the tools can also be applied to bulk-sequenced samples. We demonstrate that sgcocaller is able to generate impeccable phasing results for high-coverage datasets, on which it is more accurate and stable than existing methods, and also performs well on low-coverage single-gamete sequencing datasets for which current methods fail. Our tools achieve highly accurate results with user-friendly installation, comprehensive documentation, efficient computation times and minimal memory usage. NUCLEAR RETENTION OF LEPTIN/UPD2 REGULATES ORGANISMAL RESILIENCE TO NUTRIENT EXTREMES Adipokines released from the adipocytes function as a systemic adipometer; they impinge on neural circuits to signal nutrient status. On starvation, adipokines must be retained to signal energy deficit; else, it significantly reduces starvation survival. But how fat cells retain adipokines is unclear. Here, we demonstrate that Atg8, a cell-intrinsic autophagy factor, regulates the starvation-induced acute retention of the Leptin Drosophila ortholog Upd2. We show that on starvation, as a direct consequence of Atg8s lipidation, Upd2 accumulates in the nucleus. We illustrate that Upd2s nuclear retention is critical to fat mobilization and increased starvation resilience. Furthermore, nuclear Upd2 promotes the expression of a secreted innate immune gene signature. This hints at an unanticipated connection between adipokine nuclear retention and increased innate immunity. In conclusion, we propose that, during starvation, Atg8s role is not just limited to autophagy but is critical for withholding adipokines in the nucleus to promote starvation resilience.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=178 HEIGHT=200 SRC="FIGDIR/small/428913v5_ufig1.gif" ALT="Figure 1">
View larger version (33K):
org.highwire.dtl.DTLVardef@1826425org.highwire.dtl.DTLVardef@17cca52org.highwire.dtl.DTLVardef@1c70ac7org.highwire.dtl.DTLVardef@94a024_HPS_FORMAT_FIGEXP  M_FIG C_FIG O_LIIn fed state Upd2 requires Atg8 for nuclear exit and cytosolic localization.
C_LIO_LIAtg8s lipidation on starvation results in Upd2s nuclear accumulation.
C_LIO_LIUpd2 nuclear retention on starvation increases fat mobilization and post-starvation hunger.
C_LIO_LIOn starvation Upd2 nuclear retention increases expression of a secreted innate immune signature.
C_LI