title: Examining Batch Effect in Histopathology as a Distributionally Robust Optimization Problem

abstract: Computer vision (CV) approaches applied to digital pathology have informed biological discovery and development of tools to help inform clinical decision-making. However, batch effects in the images have the potential to introduce spurious confounders and represent a major challenge to effective analysis and interpretation of these data. Standard methods to circumvent learning such confounders include (i) application of image augmentation techniques and (ii) examination of the learning process by evaluating through external validation (e.g., unseen data coming from a comparable dataset collected at another hospital). Here, we show that the source site of a histopathology slide can be learned from the image using CV algorithms in spite of image augmentation, and we explore these source site predictions using interpretability tools. A CV model trained using Empirical Risk Minimization (ERM) risks learning this source-site signal as a spurious correlate in the weak-label regime, which we abate by using a training method with abstention. We find that a patch based classifier trained using abstention outperformed a model trained using ERM by 9.9, 10 and 19.4% F1 in the binary classification tasks of identifying tumor versus normal tissue in lung adenocarcinoma, Gleason score in prostate adenocarcinoma, and tumor tissue grade in clear cell renal cell carcinoma, respectively, at the expense of up to 80% coverage (defined as the percent of tiles not abstained on by the model). Further, by examining the areas abstained by the model, we find that the model trained using abstention is more robust to heterogeneity, artifacts and spurious correlates in the tissue. Thus, a method trained with abstention may offer novel insights into relevant areas of the tissue contributing to a particular phenotype. Together, we suggest using data augmentation methods that help mitigate a digital pathology models reliance on potentially spurious visual features, as well as selecting models that can identify features truly relevant for translational discovery and clinical decision support.

category: cancer biology


title: A data compendium of Mycobacterium tuberculosis antibiotic resistance

abstract: The Comprehensive Resistance Prediction for Tuberculosis: an International Consortium (CRyPTIC) presents here a compendium of 15,211 Mycobacterium tuberculosis global clinical isolates, all of which have undergone whole genome sequencing (WGS) and have had their minimum inhibitory concentrations to 13 antitubercular drugs measured in a single assay. It is the largest matched phenotypic and genotypic dataset for M. tuberculosis to date. Here, we provide a summary detailing the breadth of data collected, along with a description of how the isolates were collected and uniformly processed in CRyPTIC partner laboratories across 23 countries. The compendium contains 6,814 isolates resistant to at least one drug, including 2,129 samples that fully satisfy the clinical definitions of rifampicin resistant (RR), multi-drug resistant (MDR), pre-extensively drug resistant (pre-XDR) or extensively drug resistant (XDR). Accurate prediction of resistance status (sensitive/resistant) to eight antitubercular drugs by using a genetic mutation catalogue is presented along with the presence of suspected resistance-conferring mutations for isolates resistant to the newly introduced drugs bedaquiline, clofazimine, delamanid and linezolid. Finally, a case study of rifampicin mono-resistance demonstrates how this compendium could be used to advance our genetic understanding of rare resistance phenotypes. The compendium is fully open-source and it is hoped that the dataset will facilitate and inspire future research for years to come.

category: microbiology


title: In vivo serial passaging of human-simian immunodeficiency virus clones identifies viral characteristics that enhance persistent viral replication

abstract: We previously reported that a human immunodeficiency virus type 1 with a simian immunodeficiency virus vif substitution (HSIV-vifNL4-3) could replicate in pigtailed macaques (PTMs), demonstrating that Vif is a species-specific tropism factor of primate lentiviruses. However, infections did not result in high peak viremia or setpoint plasma viral loads, as observed during SIV infection of PTMs. Here, we characterized variants isolated from one of the original infected animals with CD4 depletion after nearly four years of infection to identify determinants of increased replication fitness. In our studies, we found that the HSIV-vif clones did not express the HIV-1 Vpr protein due to interference from the vpx open reading frame in singly spliced vpr mRNA. To examine whether these viral genes contribute to persistent viral replication, we generated infectious HSIV-vif clones expressing either the HIV-1 Vpr or SIV Vpx protein. And then to determine viral fitness determinants of HSIV-vif, we conducted three rounds of serial in vivo passaging in PTMs, starting with an initial inoculum containing a mixture of CXCR4-tropic (Vpr- HSIV-vifNL4-3 isolated at 196 (C/196) and 200 (C/200) weeks post-infection from a PTM with depressed CD4 counts) and CCR5-tropic HSIV (Vpr+ HSIV-vif derivatives based NL-AD8 and Bru-Yu2 and a Vpx expressing HSIV-vifYu2). Interestingly, all infected PTMs showed peak plasma viremia close to or above 105 copies/ml and persistent viral replication for more than 20 weeks. The passage 3 PTM showed peak viral loads greater than 105 viral RNA copies/ml. Infectious molecular clones (IMCs) recovered from the passage 3 PTM (HSIV-P3 IMCs) included mutations required for HIV-1 Vpr expression and those mutations encoded by the CXCR4-tropic HSIV-vifNL4-3 isolates C/196 and C/200. The data indicate that the biological isolates selected during long-term infection acquired HIV-1 Vpr expression to enhance their replication fitness in PTMs. Further passaging of HSIV-P3 IMCs in vivo may generate pathogenic variants with higher replication capacity, which will be a valuable resource as challenge virus in vaccine and cure studies.


title: Effect of sodium silicate on drinking water biofilm development

abstract:
Sodium silicates have been studied for sequestration of iron, coagulation, and corrosion control, but their impact on biofilm formation has not been documented in detail. This study investigated the impact of sodium silicate corrosion control on biomass accumulation in drinking water systems in comparison to orthophosphate, a common corrosion inhibitor. Biofilm growth was measured by determining ATP concentrations, and the bacterial community was characterized using 16S ribosomal RNA (rRNA) sequencing. A pilot-scale study with cast-iron pipe loops, annular reactors (ARs), and polycarbonate coupons demonstrated significantly lower biofilm ATP concentrations in the sodium silicate-treated AR than the orthophosphate-treated AR when the water temperature exceeded 20{degrees}C. However, an elevated sodium silicate dose (48 mg L-1 of SiO2) disturbed and dispersed the biofilm formed inside the AR, resulting in elevated effluent ATP concentrations. Two separate experiments confirmed that biomass accumulation was higher in the presence of orthophosphate at high water temperatures (20{degrees}C) only. No significant differences were identified in biofilm ATP concentrations at lower water temperatures (below 20{degrees}C). Differences in bacterial communities between the orthophosphate- and sodium silicate-treated systems were not statistically significant, even though orthophosphate promoted higher biofilm growth. However, the genera Halomonas and Mycobacterium--which include opportunistic pathogens--were present at greater relative abundances in the orthophosphate-treated system compared to the sodium silicate system.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=111 SRC="FIGDIR/small/460503v1_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@14ac86corg.highwire.dtl.DTLVardef@ac500borg.highwire.dtl.DTLVardef@bb2980org.highwire.dtl.DTLVardef@1a1b5c0_HPS_FORMAT_FIGEXP  M_FIG C_FIG Orthophosphate promotes more biofilm growth in comparison to sodium silicates at water temperatures above 20{degrees}C.

Water impact statementSodium silicates have been used in drinking water treatment for decades, both as sequestrants and as corrosion inhibitors. However, their impact on biofilm formation is poorly understood, and this risks drinking water quality. This study aims to further clarify the effects of corrosion inhibitors on biofilm development, including inhibitors that are not phosphate-based.


title: Phylodynamic Pattern of Genetic Clusters, Paradigm Shift on Spatio-temporal Distribution of Clades, and Impact of Spike Glycoprotein Mutations of SARS-CoV-2 Isolates from India

abstract: Background: The COVID-19 pandemic is associated with high morbidity and mortality, with the emergence of numerous variants. The dynamics of SARS-CoV-2 with respect to clade distribution is uneven, unpredictable and fast changing. Aims: Retrieving the complete genomes of SARS-CoV-2 from India and subjecting them to analysis on phylogenetic clade diversity, Spike (S) protein mutations and their functional consequences such as immune escape features and impact on infectivity. Methods: Whole genome of SARS-CoV-2 isolates (n=4,326) deposited from India during the period from January 2020 to December 2020 is retrieved from GISAID and various analyses performed using in silico tools. Results: Notable clade dynamicity is observed indicating the emergence of diverse SARS-CoV-2 variants across the country. GR clade is predominant over the other clades and the distribution pattern of clades is uneven. D614G is the commonest and predominant mutation found among the S-protein followed by L54F. Mutation score prediction analyses reveal that there are several mutations in S-protein including the RBD and NTD regions that can influence the virulence of virus. Besides, mutations having immune escape features as well as impacting the immunogenicity and virulence through changes in the glycosylation patterns are identified. Conclusions: The study has revealed emergence of variants with shifting of clade dynamics within a year in India. It is shown uneven distribution of clades across the nation requiring timely deposition of SARS-CoV-2 sequences. Functional evaluation of mutations in S-protein reveals their significance in virulence, immune escape features and disease severity besides impacting therapeutics and prophylaxis.

category: evolutionary biology


title: Mechanistic insights into dynamic mutual regulation of USP14 and proteasome

abstract: Proteasomal degradation of ubiquitylated proteins is sophisticatedly regulated at multiple levels1-3. A primary regulatory checkpoint is the removal of ubiquitin chains from substrates by the deubiquitylating enzyme USP14 that associates reversibly with the proteasome. How USP14 is activated and regulates the proteasome function remains unknown4-7. Here we report cryo-electron microscopy (cryo-EM) structures of human USP14 in complex with the 26S proteasome in nine conformational states at 3.0-3.6 [A] resolution, captured during polyubiquitylated protein degradation. Time-resolved cryo-EM analysis of the conformational continuum revealed two parallel pathways of proteasome state transitions induced by USP14 and captured transient conversion of substrate-engaged intermediates into substrate-inhibited intermediates. On the substrate-engaged pathway, USP14 activation allosterically reprograms conformational landscape of the AAA-ATPase motor and stimulates opening of the core particle gate8-10, enabling observation of a near-complete cycle of asymmetric ATP hydrolysis around the ATPase ring during processive substrate unfolding. Dynamic USP14-ATPase interactions decouple the ATPase activity from RPN11-catalysed deubiquitylation11-13 and kinetically introduce three regulatory checkpoints on the proteasome, at the steps of ubiquitin recognition, substrate translocation initiation and ubiquitin chain recycling. These findings provide unprecedented insights into the complete functional cycles of USP14-regulated proteasome and of USP14 activation-deubiquitylation-disassembly and establish mechanistic foundations for USP14-targeted therapeutic discovery.

category: biochemistry


title: Minos: variant adjudication and joint genotyping of cohorts of bacterial genomes

abstract:
Short-read variant calling for bacterial genomics is a mature field, and there are many widely-used software tools. Different underlying approaches (eg pileup, local or global assembly, paired-read use, haplotype use) lend each tool different strengths, especially when considering non-SNP (single nucleotide polymorphism) variation or potentially distant reference genomes. It would therefore be valuable to be able to integrate the results from multiple variant callers, using a robust statistical approach to "adjudicate" at loci where there is disagreement between callers. To this end, we present a tool, Minos, for variant adjudication by mapping reads to a genome graph of variant calls. Minos allows users to combine output from multiple variant callers without loss of precision. Minos also addresses a second problem of joint genotyping SNPs and indels in bacterial cohorts, which can also be framed as an adjudication problem.

We benchmark on 62 samples from 3 species (Mycobacterium tuberculosis, Staphylococcus aureus, Klebsiella pneumoniae) and an outbreak of 385 M. tuberculosis samples. Finally, we joint genotype a large M. tuberculosis cohort (N{approx}15k) for which the rifampicin phenotype is known. We build a map of non-synonymous variants in the RRDR (rifampicin resistance determining region) of the rpoB gene and extend current knowledge relating RRDR SNPs to heterogeneity in rifampicin resistance levels. We replicate this finding in a second M. tuberculosis cohort (N{approx}13k).

Minos is released under the MIT license, available at https://github.com/iqbal-lab-org/minos.

category: bioinformatics


title: The value of genotype-specific reference for transcriptome analyses

abstract: High quality transcriptome assembly using short reads from RNA-seq data still heavily relies upon reference-based approaches, of which the primary step is to align RNA-seq reads to a single reference genome of haploid sequence. However, it is increasingly apparent that while different genotypes within a species share "core" genes, they also contain variable numbers of "specific" genes that are only present a subset of individuals. Using a common reference may thus lead to a loss of genotype-specific information in the assembled transcript dataset and the generation of erroneous, incomplete or misleading transcriptomics analysis results. With the recent development of pan-genome information in many species, it is important that we understand the limitations of single genotype references for transcriptomics analysis. In this study, we quantitively evaluated the advantages of using genotype-specific reference genomes for transcriptome assembly and analysis using cultivated barley as a model. We mapped barley cultivar Barke RNA-seq reads to the Barke genome and to the cultivar Morex genome (common barley genome reference) to construct a genotype specific Reference Transcript Dataset (sRTD) and a common Reference Transcript Datasets (cRTD), respectively. We compared the two RTDs according to their transcript diversity, transcript sequence and structure similarity and the accuracy they provided for transcript quantification and differential expression analysis. Our evaluation shows that the sRTD has a significantly higher diversity of transcripts and alternative splicing events. Despite using a high-quality reference genome for assembly of the cRTD, we miss ca. 40% transcripts present in the sRTD and cRTD only has ca. 70% true assemblies. We found that the sRTD is more accurate for transcript quantification as well as differential expression and differential alternative splicing analysis. However, gene level quantification and comparative expression analysis are less affected by the source RTD, which indicates that analysing transcriptomic data at the gene level may be a reasonable compromise when a high-quality genotype-specific reference is not available.


title: Mechanism of broad-spectrum Cas9 inhibition by AcrIIA11

abstract: Mobile genetic elements evade CRISPR-Cas adaptive immunity by encoding anti-CRISPR proteins (Acrs). Acrs inactivate CRISPR-Cas systems via diverse mechanisms but are generally specific for a narrow subset of Cas nucleases that share high sequence similarity. Here, we demonstrate that AcrIIA11 inhibits diverse Cas9 sub-types in vitro and human cells. Single-molecule fluorescence imaging reveals that AcrIIA11 interferes with the first steps of target search by reducing S. aureus Cas9s diffusion on non-specific DNA. DNA cleavage is inhibited because the AcrIIA11:Cas9 complex is kinetically trapped at PAM-rich decoy sites, preventing Cas9 from reaching its target. This work establishes that DNA trapping can be used to inhibit a broad spectrum of Cas9 orthologs in vitro and during mammalian genome editing.

category: biophysics


title: Transitioning to confined spaces impacts bacterial swimming and escape response

abstract:
Symbiotic bacteria often navigate complex environments before colonizing privileged sites in their host organism. Chemical gradients are known to facilitate directional taxis of these bacteria, guiding them towards their eventual destination. However, less is known about the role of physical features in shaping the path the bacteria take and defining how they traverse a given space. The flagellated marine bacterium Vibrio fischeri,which forms a binary symbiosis with the Hawaiian bobtail squid, Euprymna scolopes, must navigate tight physical confinement, squeezing through a bottleneck constricting to ~2 m in width on the way to its eventual home. Using microfluidic in vitro experiments, we discovered that V. fischeri cells alter their behavior upon entry into confined space, straightening their swimming paths and promoting escape from confinement. Using a computational model, we attributed this escape response to two factors: reduced directional fluctuation and a refractory period between reversals. Additional experiments in asymmetric capillary tubes confirmed that V. fischeri quickly escape from tapered ends, even when drawn into the ends by chemoattraction. This avoidance was apparent down to a limit of confinement approaching the diameter of the cell itself, resulting in a balance between chemoattraction and evasion of physical confinement. Our findings demonstrate that non-trivial distributions of swimming bacteria can emerge from simple physical gradients in the level of confinement. Tight spaces may serve as an additional, crucial cue for bacteria while they navigate complex environments to enter specific habitats.

Significance StatementSymbiotic bacteria that navigate to and through specific host tissues often face tight physical confinement. This work reveals that confinement-associated changes in swimming can dramatically alter taxis, shaping bacterial localization in conjuncture with other motility-directing cues. This work helps explain how bacteria can avoid getting stuck in confined areas while transiting to privileged spaces, adding confinement as an environmental cue that symbiotic bacteria use to shape their motility behavior.


title: Molecular architecture of SARS-CoV-2 envelope by integrative modeling

abstract: Despite tremendous efforts by research community during the COVID-19 pandemic, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M (membrane) protein. The molecular dynamics simulations allowed us to test the envelop stability under different configurations and revealed that M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns formed by Ms transmembrane and intra-virion (endo) domains. These results were in agreement with the experimental data, demonstrating a generic and versatile integrative approach to model the structure of a virus de novo, providing insights into critical roles of structural proteins in the viral assembly and integration, and proposing new targets for the antiviral therapies.


title: Hydrogen sulfide release via the ACE inhibitor Zofenopril prevents intimal hyperplasia in human vein segments and in a mouse model of carotid artery stenosis

abstract:
ObjectivesHypertension is a major risk factor for intimal hyperplasia (IH) and restenosis following vascular and endovascular interventions. Pre-clinical studies suggest that hydrogen sulfide (H2S), an endogenous gasotransmitter, limits restenosis. While there is no clinically available pure H2S releasing compound, the sulfhydryl-containing angiotensin-converting enzyme inhibitor Zofenopril is a source of H2S. Here, we hypothesized that Zofenopril, due to H2S release, would be superior to other non-sulfhydryl containing angiotensin converting enzyme inhibitor (ACEi), in reducing intimal hyperplasia.

MaterialsSpontaneously hypertensive male Cx40 deleted mice (Cx40-/-) or WT littermates were randomly treated with Enalapril 20 mg (Mepha Pharma) or Zofenopril 30 mg (Mylan SA). Discarded human vein segments and primary human smooth muscle cells (SMC) were treated with the active compound Enalaprilat or Zofenoprilat.

MethodsIH was evaluated in mice 28 days after focal carotid artery stenosis surgery and in human vein segments cultured for 7 days ex vivo. Human primary smooth muscle cell (SMC) proliferation and migration were studied in vitro.

ResultsCompared to control animals (intima/media thickness=2.3{+/-}0.33), Enalapril reduced IH in Cx40-/- hypertensive mice by 30% (1.7{+/-}0.35; p=0.037), while Zofenopril abrogated IH (0.4{+/-}0.16; p<.0015 vs. Ctrl and p>0.99 vs. sham-operated Cx40-/- mice). In WT normotensive mice, enalapril had no effect (0.9665{+/-}0.2 in control vs 1.140{+/-}0.27; p>.99), while Zofenopril also abrogated IH (0.1623{+/-}0.07, p<.008 vs. Ctrl and p>0.99 vs. sham-operated WT mice). Zofenoprilat, but not Enalaprilat, also prevented intimal hyperplasia in human veins segments ex vivo. The effect of Zofenopril on carotid and SMC correlated with reduced SMC proliferation and migration. Zofenoprilat inhibited the MAPK and mTOR pathways in SMC and human vein segments.

ConclusionZofenopril provides extra beneficial effects compared to non-sulfhydryl ACEi to reduce SMC proliferation and restenosis, even in normotensive animals. These findings may hold broad clinical implications for patients suffering from vascular occlusive diseases and hypertension.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=100 SRC="FIGDIR/small/460108v1_ufig1.gif" ALT="Figure 1">
View larger version (48K):
org.highwire.dtl.DTLVardef@a37b62org.highwire.dtl.DTLVardef@f977e0org.highwire.dtl.DTLVardef@2d3266org.highwire.dtl.DTLVardef@147d604_HPS_FORMAT_FIGEXP  M_FIG C_FIG What this paper addsThe current strategies to reduce intimal hyperplasia (IH) principally rely on local drug delivery, in endovascular approach. The oral angiotensin converting enzyme inhibitor (ACEi) Zofenopril has additional effects compared to other non-sulfyhydrated ACEi to prevent intimal hyperplasia and restenosis. Given the number of patients treated with ACEi worldwide, these findings call for further prospective clinical trials to test the benefits of sulfhydrated ACEi over classic ACEi for the prevention of restenosis in hypertensive patients.

category: physiology


title: The conventional dendritic cell 1 subset primes CD8+ T cells and traffics tumor antigen to drive anti-tumor immunity in the brain

abstract:
The central nervous system (CNS) antigen presenting cell (APC) which primes anti-tumor CD8+ T cell responses remains undefined. Elsewhere, the conventional dendritic cell 1 (cDC1) performs this role. However, steady-state brain cDC1 are rare; cDC localize to choroid plexus and dura. Using preclinical glioblastoma models and cDC1-deficient mice, we explored the role of cDC1 in CNS anti-tumor immunity. We determined that cDC1 mediate checkpoint blockade-induced survival benefit and prime neoantigen-specific CD8+ T cells against brain tumors. We observed that cDC, including cDC1, isolated from the tumor, the dura, and the CNS-draining cervical lymph nodes harbored a traceable fluorescent tumor-antigen. In patient samples, we observed several APC subsets (including the CD141+ cDC1-equivalent) infiltrating glioblastomas, meningiomas, and dura. In these same subsets, we identified a tumor-specific fluorescent metabolite of 5- aminolevulinic acid, which labels tumor cells during fluorescence-guided glioblastoma resection. Together, these data elucidate the specialized behavior of cDC1 and suggest cDC1 play a significant role in CNS anti-tumor immunity.

One Sentence SummaryCNS cDC1 engage in previously undefined behavior to establish immune responses against brain tumors.

category: immunology


title: Identification of CRF66_BF, a new HIV-1 circulating recombinant form of South American origin

abstract: Circulating recombinant forms (CRFs) are important components of the HIV-1 pandemic. Among 108 reported in the literature, 17 are BF1 intersubtype recombinant, most of which are of South American origin. Among these, all 5 identified in the Southern Cone and neighboring countries, except Brazil, derive from a common recombinant ancestor related to CRF12_BF, which circulates widely in Argentina, as deduced from coincident breakpoints and clustering in phylogenetic trees. In a HIV-1 molecular epidemiological study in Spain, we identified a phylogenetic cluster of 20 samples from 3 separate regions which were of F1 subsubtype, related to the Brazilian strain, in protease-reverse transcriptase (Pr-RT) and of subtype B in integrase. Remarkably, 14 individuals from this cluster (designated BF9) were Paraguayans and only 4 were native Spaniards. HIV-1 transmission was predominantly heterosexual, except for a subcluster of 6 individuals, 5 of which were men who have sex with men. Ten additional database sequences, from Argentina (n=4), Spain (n=3), Paraguay (n=1), Brazil (n=1), and Italy (n=1), branched within the BF9 cluster. To determine whether it represents a new CRF, near full-length genome (NFLG) sequences were obtained for 6 viruses from 3 Spanish regions. Bootscan analyses showed a coincident BF1 recombinant structure, with 5 breakpoints, located in p17gag, integrase, gp120, gp41-rev overlap, and nef, which was identical to that of two BF1 recombinant viruses from Paraguay previously sequenced in NFLGs. Interestingly, none of the breakpoints coincided with those of CRF12_BF. In a maximum likelihood phylogenetic tree, all 8 NFLG sequences grouped in a strongly supported clade segregating from previously identified CRFs and from the CRF12_BF  family clade. These results allow us to identify a new HIV-1 CRF, designated CRF66_BF. Through a Bayesian coalescent analysis, the most recent common ancestor of CRF66_BF was estimated around 1984 in South America, either in Paraguay or Argentina. Among Pr-RT sequences obtained by us from HIV-1-infected Paraguayans living in Spain, 14 (20.9%) of 67 were of CRF66_BF, suggesting that CRF66_BF may be one of the major HIV-1 genetic forms circulating in Paraguay. CRF66_BF is the first reported non-Brazilian South American HIV-1 CRF_BF unrelated to CRF12_BF.


title: Characterization of the microbiota associated with 12-week-old bovine fetuses exposed to divergent in utero nutrition

abstract: A recent study reported the existence of a diverse microbiota in 5-to-7-month-old calf fetuses, suggesting that colonization of the bovine gut with so-called "pioneer" microbiota may begin during mid-gestation. In the present study, we investigated 1) the presence of microbiota in bovine fetuses at early gestation (12 weeks), and 2) whether the fetal microbiota is influenced by the maternal rate of gain or dietary supplementation with vitamins and minerals (VTM) during early gestation. Amniotic and allantoic fluids, and intestinal and placental (cotyledon) tissue samples obtained from fetuses (n = 33) on day 83 of gestation were processed for the assessment of fetal microbiota using 16S rRNA gene sequencing. The sequencing results revealed that a diverse and complex microbial community was present in each of these fetal compartments evaluated. Allantoic and amniotic fluids, and fetal intestinal and placenta microbiota each had distinctly different (0.047 [≥] R2 [≥] 0.019, P [≤] 0.031) microbial community structures. Allantoic fluid had a greater (P < 0.05) microbial richness (number of OTUs) (Mean 122) compared to amniotic fluid (84), intestine (63) and placenta (66). Microbial diversity (Shannon index) was similar for the intestinal and placental samples, and both were less diverse compared with fetal fluid microbiota (P < 0.05). Thirty-nine different archaeal and bacterial phyla were detected across all fetal samples, with Proteobacteria (55%), Firmicutes (16.2%), Acidobacteriota (13.6%) and Bacteroidota (5%) predominating. Among the 20 most relatively abundant bacterial genera, Acidovorax, Acinetobacter, Brucella, Corynebacterium, Enterococcus, Exiguobacterium and Stenotrophomonas differed by fetal sample type (P < 0.05). A total of 55 taxa were shared among the four different microbial communities. qPCR of bacteria in the intestine and placenta samples as well as scanning electron microscopy imaging of fetal fluids provided additional evidence for the presence of a microbiota in these samples. Minor effects of maternal rate of gain and VTM supplementation, and their interactions on microbial richness and composition were detected. Overall, the results of this study indicate that colonization with pioneer microbiota may occur during early gestation in bovine fetuses, and that that the maternal nutritional regime during gestation may influence the early fetal microbiota.


title: Poly (acetyl arginyl) glucosamine attenuates Pseudomonas aeruginosa in a rat lung infection model

abstract: Pseudomonas aeruginosa is a common opportunistic pathogen that can cause chronic infections in multiple disease states, including respiratory infections in patients with cystic fibrosis (CF) and non-CF bronchiectasis. Like many opportunists, P. aeruginosa forms multicellular biofilm communities that are widely thought to be an important determinant of bacterial persistence and resistance to antimicrobials and host immune effectors during chronic/recurrent infections. Poly (acetyl, arginyl) glucosamine (PAAG) is a glycopolymer which has antimicrobial activity against a broad range of bacterial species, and also has mucolytic activity which can normalize rheologic properties of cystic fibrosis mucus. In this study, we sought to evaluate the effect of PAAG on P. aeruginosa bacteria within biofilms in vitro, and in the context of experimental pulmonary infection in a rodent infection model. PAAG treatment caused significant bactericidal activity against P. aeruginosa biofilms, and a reduction in the total biomass of preformed P. aeruginosa biofilms on abiotic surfaces, as well as on the surface of immortalized cystic fibrosis human bronchial epithelial cells. Studies of membrane integrity indicated that PAAG causes changes to P. aeruginosa cell morphology and dysregulates membrane polarity. PAAG treatment reduced infection and consequent tissue inflammation in experimental P. aeruginosa rat infections. Based on these findings we conclude that PAAG represents a novel means to combat P. aeruginosa infection, which may warrant further evaluation as a therapeutic.


title: Revealing spatio-temporal dynamics with long-term trypanosomatid live-cell imaging

abstract:
Trypanosoma brucei, the causative agent of human African trypanosomiasis, employs a flagellum for dissemination within the parasites mammalian and insect hosts. T. brucei cells are highly motile in culture and must be able to move in all three dimensions for reliable cell division. These characteristics have made long-term microscopic imaging of live T. brucei cells challenging, which has limited our understanding of a variety of important cell-cycle events. To address this issue, we have devised an imaging approach that confines cells to small volumes that can be imaged continuously for up to 24 h. This system employs cast agarose microwells generated using a PDMS stamp that can be made with different dimensions to maximize cell viability and imaging quality. Using this approach, we have imaged individual T. brucei through multiple rounds of cell division with high spatial and temporal resolution. We have employed this method to study the differential rate of T. brucei daughter cell division and show that the approach is compatible with loss-of-function experiments such as small molecule inhibition and RNAi. We have also developed a strategy that employs in-well "sentinel" cells to monitor potential toxicity due to imaging. This live-cell imaging method will provide a novel avenue for studying a wide variety of cellular events in trypanosomatids that have previously been inaccessible.

ImportanceTrypanosoma brucei causes severe diseases that affect humans and livestock in Sub-Saharan Africa. Efficient strategies for manipulating the T. brucei genome have provided a wealth of information about protein localization and function in diverse cellular processes. However, employing live-cell imaging for phenotypic analysis in T. brucei remains a significant challenge because immobilization of this highly motile parasite rapidly leads to morphologic defects and cell death. While fixed-cell imaging can provide snapshots of cellular events, it cannot provide the direct causal link or precise timing of events that comes from watching a living cell change over time. Our strategy using agarose microwells now allows long-term live cell T. brucei imaging with a simple apparatus that can be adapted for a wide range of experimental conditions.


title: Single cell infection with influenza A virus using drop-based microfluidics

abstract: Influenza A virus (IAV) is an RNA virus with high genetic diversity which necessitates the development of new vaccines targeting emerging mutations each year. As IAV exists in genetically heterogeneous populations, current studies focus on understanding population dynamics at the single cell level. These studies include novel methodology that can be used for probing populations at the single cell level, such as single cell sequencing and microfluidics. Here, we introduce a drop-based microfluidics method to study IAV infection at a single cell level by isolating infected host cells in microscale drops. Single human alveolar basal epithelial (A549), Madin-Darby Canine Kidney cells (MDCK) and MDCK + human siat7e gene (Siat7e) cells infected with the pandemic A/California/07/2009 (H1N1) strain were encapsulated within 50 m radii drops and incubated at 37{degrees}C. We demonstrate that drops remain stable over 24 hours, that 75% of cells remain viable, and that IAV virus can propagate within the drops. Drop-based microfluidics therefore enables single cell analysis of viral populations produced from individually infected cells.


title: Dynamic behavior of cell-free mitochondrial DNA in human saliva

abstract: Mitochondria contain their own genome that can be released in multiple biofluids such as blood and cerebrospinal fluid, as cell-free mitochondrial DNA (cf-mtDNA). In clinical studies, single measures of blood cf-mtDNA predict mortality, and higher cf-mtDNA levels are associated with mental and physical stress. However, the dynamics of cf-mtDNA has not been defined, and whether it can be measured non-invasively like other neuroendocrine markers in saliva has not been examined. Here we report cf-mtDNA in human saliva and establish its natural within-person dynamic behavior across multiple weeks. In a small proof-of-principle cohort of healthy adults, we first develop an approach to rapidly quantify salivary cf-mtDNA without DNA extraction, and demonstrate the existence of saliva cf-mtDNA. We then deploy this approach to perform an intensive repeated-measures analysis of two healthy men studied at 4 daily timepoints over 53-60 consecutive days (n=212-220 observations each) with parallel measures of steroid hormones, self-reported daily mood, and health-related behaviors. Salivary cf-mtDNA exhibited a robust awakening response reaching up to two orders of magnitude 30-45 minutes after awakening, varied from day-to-day, and moderately correlated with the cortisol awakening response. No consistent association with self-reported daily mood/health-related behaviors were found, although this requires further examination in more extensive studies. Dynamic variation in cf-mtDNA was inversely related with salivary interleukin 6 (IL6), inconsistent with a pro-inflammatory effect of salivary cf-mtDNA. The highly dynamic behavior of salivary cf-mtDNA opens the door to non-invasive studies examining the relevance of mtDNA signaling in relation to human health.

category: molecular biology


title: Fast and reliable sgRNA efficiency testing using HIReff

abstract: CRISPR/Cas9 is the method of choice for gene editing like the endogenous knock-in of sequences in order to investigate protein function, abundance or intracellular localization. One of the crucial steps in the preparation of CRISPR/Cas9-mediated knock-ins is the design of sgRNAs, which need to be tested carefully in order to minimize off-target binding and reach highest cleavage efficiency. Usually, sgRNA is evaluated via mismatch cleavage assays, like Surveyor or T7 endonuclease 1 assay. We demonstrate that these methods are often highly cost- and time-intensive with a low sensitivity and high fail rate. As an alternative, we present a new HITI-based sgRNA efficiency (HIReff) test to precisely evaluate sgRNA efficiency. HIReff is based on a sophisticated integration vector with on-site generation of a linear donor fragment that allows a comparably easy read-out via fluorescence signal and integrates several internal controls. Next to a quantifiable sgRNA assessment, HIReff provides additional information on the  gene/protein to be studied abundance, subcellular localization and promoter activity and allows derivation of fluorescence protein labeled clonal lines. We highlight benefits of HIReff in comparison to commonly used enzyme-based assays and demonstrate improved practicability and high sensitivity, while being less time-, labor- and costintensive at the same time. Our results suggest HIReff as a fast and easy-to-use alternative for sgRNA efficiency testing.


title: Human - murine concordance of molecular signatures in nerve-sparing murine partial bladder outlet obstruction (NeMO)

abstract: Recently we demonstrated the utility of a nerve-sparing mid-urethra model of partial outlet obstruction (NeMO) that has high consistency and minimal mortalities, unlike the traditional model proximal to the bladder neck. Our goal was to uncover potential therapeutic targets by investigating the genome wide transcriptional changes and pathways altered in NeMO to compare with published human bladder obstruction data. We performed RNAseq and analysed the differentially upregulated and downregulated genes for associated pathways, transcription factor binding site analysis (TFBS), upstream regulators and Gene Set Enrichment Analysis (GSEA). NeMO increased bladder mass, relative bladder mass and hyperactivity, and decreased voiding efficiency. In NeMO vs. sham, 831 genes were differentially expressed (adjusted p<0.05) and correlated significantly with at least one physiologic parameter. Gene ontology revealed an enrichment for matrix pathways in the upregulated genes, and for cardiac contraction, oxidative phosphorylation and pyruvate metabolism in downregulated genes. TFBS analysis revealed a differential regulation of up vs downregulated genes, with KLF4 strongly associated with the downregulated genes. Downregulated genes of Human bladder obstruction were also associated with the TFBS of KLF4. GSEA of the NeMO gene set confirmed the DAVID results, but also showed a cluster of cytokine activation genes. In human bladder underactive obstruction, cytokines were also highly upregulated. The common cytokine pathway upregulation provided an example of the use of RNAseq for uncovering potential new therapeutic targets. As TNF and the innate immune pathways were strongly implicated in both human and mouse, and TNF is produced by macrophages, we depletion macrophages with clodronate (CL) during NeMO. Although CL did not block hypertrophy, it significantly decreased NeMO-induced hyperactive voiding (p<0.01) and increased voiding efficiency (p<0.05). The expression of several cytokines/chemokines correlated significantly with bladder functional parameters such as residual volumes, and hyperactivity. Conclusions: Gene expression signatures of NeMO were consistent with human bladder obstruction, supporting the use of the nerve-sparing mouse obstruction model for therapeutic exploration.


title: Evaluation of microbiome enrichment and host DNA depletion in human vaginal samples using Oxford Nanopore's adaptive sequencing

abstract: Metagenomic sequencing is promising for clinical applications to study microbial composition concerning disease or patient outcomes. Alterations of the vaginal microbiome are associated with adverse pregnancy outcomes, like preterm premature rupture of membranes and preterm birth. Methodologically these samples often have to deal with low relative amounts of prokaryotic DNA and high amounts of host DNA (> 90%), decreasing the overall microbial resolution. Nanopores adaptive sampling method offers selective DNA depletion or target enrichment to directly reject or accept DNA molecules during sequencing without specialized sample preparation. Here, we demonstrate how selective  human host depletion resulted in a 1.70 fold ({+/-} 0.27 fold) increase in total sequencing depth, providing higher taxonomic profiling sensitivity. At the same time, the microbial composition remains consistent with the control experiments. The complete removal of all human host sequences is not yet possible and should be considered as an ethical approval statement might still be necessary. Adaptive sampling increased microbial sequencing yield in all 15 sequenced clinical routine vaginal samples, making it a valuable tool for clinical surveillance and medical-based research, which can be used in addition to other host depletion methods before sequencing.


title: CD56 at the human NK cell lytic immunological synapse

abstract: CD56 is the main identifying cell surface molecule of NK cells and has been recently identified as a regulator of cytotoxic function in NK cell lines. Despite its newly defined role in lytic granule polarization and exocytosis, biological questions remain involving the localization and function of CD56 at the immunological synapse. Here we use confocal and structured illumination microscopy to demonstrate recruitment of CD56 to the peripheral supramolecular activating cluster (pSMAC) of the immunological synapse of lytic effector cells. We provide additional data demonstrating that cell lines that are less dependent on CD56 for function are not utilizing alternative pathways of cytotoxicity, and that those that are dependent on CD56 have normal expression of activating and adhesion receptors. Finally, we use actin reporter (LifeAct) expressing NK92 cell lines and live cell confocal microscopy to visualize live cell killing events with WT and CD56-KO cells. This work further characterizes the novel role for CD56 in cytotoxic function of NK cells and provides deeper insight into the role of CD56 at the NK cell immunological synapse.

category: cell biology


title: Intronization enhances expression of S-protein and other transgenes challenged by cryptic splicing

abstract: The natural habitat of SARS-CoV-2 is the cytoplasm of a mammalian cell where it replicates its genome and expresses its proteins. While SARS-CoV-2 genes and hence its codons are presumably well optimized for mammalian protein translation, they have not been sequence optimized for nuclear expression. The cDNA of the Spike protein harbors over a hundred predicted splice sites and produces mostly aberrant mRNA transcripts when expressed in the nucleus. While different codon optimization strategies increase the proportion of full-length mRNA, they do not directly address the underlying splicing issue with commonly detected cryptic splicing events hindering the full expression potential. Similar splicing characteristics were also observed in other transgenes. By inserting multiple short introns throughout different transgenes, significant improvement in expression was achieved, including >7-fold increase for Spike transgene. Provision of a more natural genomic landscape offers a novel way to achieve multi-fold improvement in transgene expression.

category: genetics


title: Non-enzymatic assembly of active chimeric ribozymes from aminoacylated RNA oligonucleotides

abstract:
Aminoacylated tRNAs, which harbor a covalent linkage between amino acids and RNA, are a universally conserved feature of life. Because they are essential substrates for ribosomal translation, aminoacylated oligonucleotides must have been present in the RNA World prior to the evolution of the ribosome. One possibility we are exploring is that the aminoacyl ester linkage served another function before being recruited for ribosomal protein synthesis. The nonenzymatic assembly of ribozymes from short RNA oligomers under realistic conditions remains a key challenge in demonstrating a plausible pathway from prebiotic chemistry to the RNA World. Here, we show that aminoacylated RNAs can undergo template-directed assembly into chimeric amino acid-RNA polymers that are active ribozymes. We demonstrate that such chimeric polymers can retain the enzymatic function of their all-RNA counterparts by generating chimeric hammerhead, RNA ligase, and aminoacyl transferase ribozymes. Amino acids with diverse side chains form linkages that are well tolerated within the RNA backbone, potentially bringing novel functionalities to ribozyme catalysis. Our work suggests that aminoacylation chemistry may have played a role in primordial ribozyme assembly. Increasing the efficiency of this process provides an evolutionary rationale for the emergence of sequence and amino acid specific aminoacyl-RNA synthetase ribozymes, which could then have generated the substrates for ribosomal protein synthesis.

Significance StatementThe emergence of the primordial ribosome from the RNA World would have required access to aminoacylated RNA substrates. The spontaneous generation of such substrates without enzymes is inefficient, and it remains unclear how they could be selected for in a prebiotic milieu. In our study we identify a role for aminoacylated RNA in ribozyme assembly, a longstanding problem in the origin of life research. We show that aminoacylated RNAs, but not unmodified RNAs, rapidly assemble into chimeric amino acid-bridged ribozymes that retain their native enzymatic activity. Our work therefore addresses two key challenges within the origin-of-life field: we demonstrate assembly of functional ribozymes and we identify a potential evolutionary benefit for RNA aminoacylation that is independent of coded peptide translation.


title: Decoding the Mechanism of Specific RNA Targeting by Ribosomal Methytransferases

abstract: Methylation of specific nucleotides is integral for ribosomal biogenesis and serves as a common way to confer antibiotic resistance by pathogenic bacteria. Here, by determining the high-resolution structure of 30S-KsgA by cryo-EM, a state was captured, where KsgA juxtaposes between helices h44 and h45, separating them, thereby enabling remodeling of the surrounded rRNA and allowing the cognate site to enter the methylation pocket. With the structure as a guide, factors that direct the enzyme to its cognate site with high fidelity were unearthed by creating several mutant versions of the ribosomes, where interacting bases in the catalytic helix h45 and surrounding helices h44, h24, and h27 were mutated and evaluated for their methylation efficiency. The biochemical studies delineated specificity hotspots that enable KsgA to achieve an induced fit. This study enables the identification of distal exclusive allosteric pocket and other divergent structural elements in each rMTase, which can be exploited to develop strategies to reverse methylation, mediated drug resistance.


title: Stable breast cancer prognosis

abstract:
Predicting breast cancer prognosis helps improve the treatment and management of the disease. In the last decades, many prediction models have been developed for breast cancer prognosis based on transcriptomic data. A common assumption made by these models is that the test and training data follow the same distribution. However, in practice, due to the heterogeneity of breast cancer and the different environments (e.g. hospitals) where data are collected, the distribution of the test data may shift from that of the training data. For example, new patients likely have different breast cancer stage distribution from those in the training dataset. Thus these existing methods may not provide stable prediction performance for breast cancer prognosis in situations with the shift of data distribution. In this paper, we present a novel stable prediction method for reliable breast cancer prognosis under data distribution shift. Our model, known as Deep Global Balancing Cox regression (DGBCox), is based on the causal inference theory. In DGBCox, firstly high-dimensional gene expression data is transferred to latent network-based representations by a deep auto-encoder neural network. Then after balancing the latent representations using a proposed causality-based approach, causal latent features are selected for breast cancer prognosis. Causal features have persistent relationships with survival outcomes even under distribution shift across different environments according to the causal inference theory. Therefore, the proposed DGBCox method is robust and stable for breast cancer prognosis. We apply DGBCox to 12 test datasets from different breast cancer studies. The results show that DGBCox outperforms benchmark methods in terms of both prediction accuracy and stability. We also propose a permutation importance algorithm to rank the genes in the DGBCox model. The top 50 ranked genes suggest that the cell cycle and the organelle organisation could be the most relevant biological processes for stable breast cancer prognosis.

Author summaryVarious prediction models have been proposed for breast cancer prognosis. The prediction models usually train on a dataset and predict the survival outcomes of patients in new test datasets. The majority of these models share a common assumption that the test and training data follow the same distribution. However, as breast cancer is a heterogeneous disease, the assumption may be violated in practice. In this study, we propose a novel method for reliable breast cancer prognosis when the test data distribution shifts from that of the training data. The proposed model has been trained on one dataset and applied to twelve test datasets from different breast cancer studies. In comparison with the benchmark methods in breast cancer prognosis, our model shows better prediction accuracy and stability. The top 50 important genes in our model provide clues to the relationship between several biological mechanisms and clinical outcomes of breast cancer. Our proposed method in breast cancer can potentially be adapted to apply to other cancer types.


title: MetaTrass: High-quality Metagenomic Taxonomic Read Assembly of Single-Species based on co-barcoding sequencing data and references

abstract: With the development of sequencing technologies and computational analysis in metagenomics, the genetic diversity of non-conserved regions has been receiving intensive attention to unravel the human gut microbial community. However, it remains a challenge to obtain enough microbial draft genomes at a high resolution from a single sample. In this work, we presented MetaTrass with a strategy of binning first and assembling later to assemble high-quality draft genomes based on metagenomics co-barcoding reads and the public reference genomes. We applied the tool to the single tube long fragment reads datasets for four human faecal samples, and generated more high-quality draft genomes with longer contiguity and higher resolution than the common combination strategies of genome assembling and binning. A total of 178 high-quality genomes was successfully assembled by MetaTrass, but the maximum of 58 was generated by the optimal common combination strategy in our tests. These high-quality genomes paved the way for genetic diversity and lineage analysis among different samples. With the high capability of assembling high-quality genomes of metagenomics datasets, MetaTrass will facilitate the study of spatial characters and dynamics of complex microbial communities at high resolution. The open-source code of MetaTrass is available at https://github.com/BGI-Qingdao/MetaTrass.


title: MetaTrass: High-quality metagenome assembly on the human gut microbiome by co-barcoding sequencing reads


title: DeepREAL: A Deep Learning Powered Multi-scale Modeling Framework Towards Predicting Out-of-distribution Receptor Activity of Ligand Binding

abstract: Drug discovery has witnessed intensive exploration of the problem of drug-target physical interactions over two decades, however, a strong drug binding affinity to a single target often fails to translate into desired clinical outcomes. A critical knowledge gap needs to be filled for correlating drug-target interactions with phenotypic responses: predicting the receptor activities or function selectivity upon the ligand binding (i.e., agonist vs. antagonist) on a genome-scale and for novel chemicals. Two major obstacles compound the difficulty on this direction: known data of receptor activity is far too scarce to train a robust model in light of genome-scale applications, and real-world applications need to deploy a model on data from various shifted distributions. To address these challenges, we have developed an end-to-end deep learning framework, DeepREAL, for multi-scale modeling of genome-wide receptor activities of ligand binding. DeepREAL utilizes self-supervised learning on tens of millions of protein sequences and pre-trained binary interaction classification to solve the data distribution shift and data scarcity problems. Extensive benchmark studies that simulate real-world scenarios demonstrate that DeepREAL achieves state-of-the-art performance in out-of-distribution settings.


title: Epitope scaffolding using alpha-synuclein cyclic peptides to generate oligomer-selective antibodies for Parkinson's disease

abstract: Effectively scaffolding epitopes on immunogens, in order to raise conformationally selective antibodies through active immunization, is a central problem in treating protein misfolding diseases, particularly neurodegenerative diseases such as Alzheimers disease or Parkinsons disease. We seek to selectively target conformations enriched in toxic, oligomeric propagating species while sparing the healthy forms of the protein that are often more abundant. To this end, we scaffolded epitopes in cyclic peptides by varying the number of flanking glycines, to best mimic a misfolding-specific conformation of an epitope of -synuclein enriched in the oligomer ensemble, as characterized by a region most readily disordered and solvent-exposed in a stressed, partially denatured protofibril. We screen and rank the cyclic peptide scaffolds of -synuclein in silico based on their ensemble overlap properties with the fibril, oligomer-model, and isolated monomer ensembles. We introduce a method for screening against structured off-pathway targets in the human proteome, by selecting scaffolds with minimal conformational similarity between their epitope and the same solvent-exposed primary sequence in structured human proteins. Different cyclic peptide scaffolds with variable numbers of glycines can have markedly different conformational ensembles. Ensemble comparison and overlap was quantified by the Jensen-Shannon Divergence, and a new measure introduced here--the embedding depth, which determines the extent to which a given ensemble is subsumed by another ensemble, and which may be a more useful measure in sculpting the conformational-selectivity of an antibody.


title: Dysregulation of the Tweak/Fn14 pathway in skeletal muscle of spinal muscular atrophy mice

abstract: Spinal muscular atrophy (SMA) is a childhood neuromuscular disorder caused by depletion of the survival motor neuron (SMN) protein. SMA is characterized by the selective death of spinal cord motor neurons, leading to progressive muscle wasting. Loss of skeletal muscle in SMA is a combination of denervation-induced muscle atrophy and intrinsic muscle pathologies. Elucidation of the pathways involved is essential to identify the key molecules that contribute to and sustain muscle pathology. The tumor necrosis factor-like weak inducer of apoptosis (TWEAK)/TNF receptor superfamily member fibroblast growth factor inducible 14 (Fn14) pathway has been shown to play a critical role in the regulation of denervation-induced muscle atrophy as well as muscle proliferation, differentiation and metabolism in adults. However, it is not clear whether this pathway would be important in highly dynamic and developing muscle. We thus investigated the potential role of the TWEAK/Fn14 pathway in SMA muscle pathology, using the severe Taiwanese Smn-/-;SMN2 and the less severe Smn2B/- SMA mice, which undergo a progressive neuromuscular decline in the first three post-natal weeks. Here, we report significantly dysregulated expression of the TWEAK/Fn14 pathway during disease progression in skeletal muscle of the two SMA mouse models. In addition, siRNA-mediated Smn knockdown in C2C12 myoblasts suggests a genetic interaction between Smn and the TWEAK/Fn14 pathway. Further analyses of SMA, Tweak-/- and Fn14-/- mice revealed dysregulated myopathy, myogenesis and glucose metabolism pathways as a common skeletal muscle feature, and providing further evidence in support of a relationship between the TWEAK/Fn14 pathway and Smn. Finally, a pharmacological intervention (Fc-TWEAK) to upregulate the activity of the TWEAK/Fn14 pathway improved disease phenotypes in the two SMA mouse models. Our study provides novel mechanistic insights into the molecular players that contribute to muscle pathology in SMA and into the role of the TWEAK/Fn14 pathway in developing muscle.

category: pathology


title: Identifying locations susceptible to micro-anatomical reentry using a spatial network representation of atrial fibre maps

abstract:
Micro-anatomical reentry has been identified as a potential driver of atrial fibrillation (AF). In this paper, we introduce a novel computational method which aims to identify which atrial regions are most susceptible to micro-reentry. The approach, which considers the structural basis for micro-reentry only, is based on the premise that the accumulation of electrically insulating interstitial fibrosis can be modelled by simulating percolation-like phenomena on spatial networks. Our results suggest that at high coupling, where micro-reentry is rare, the micro-reentrant substrate is highly clustered in areas where the atrial walls are thin and have convex wall morphology. However, as transverse connections between fibres are removed, mimicking the accumulation of interstitial fibrosis, the substrate becomes less spatially clustered, and the bias to forming in thin, convex regions of the atria is reduced. Comparing our algorithm on image-based models with and without atrial fibre structure, we find that strong longitudinal fibre coupling can suppress the micro-reentrant substrate, whereas regions with disordered fibre orientations have an enhanced risk of micro-reentry. We suggest that with further development, these methods may have future potential for patient-specific risk stratification, taking a longitudinal view of the development of the micro-reentrant substrate.

Author summaryAtrial fibrillation (AF) is the most common abnormal heart rhythm, yet, despite extensive research, treatment success rates remain poor. In part, this is because there is an incomplete understanding of the mechanistic origin of AF. In this paper, we investigate one proposed mechanism of AF, the formation of "micro-reentrant circuits", which can be thought of as a "short circuit", forming when electrically insulating fibrosis (structural repair tissue) infiltrates the space between heart muscle cells. Previously, such circuits have been found in experimental hearts, but identifying these circuits clinically is difficult. Here, we aim to take a small step towards developing computational methods for identifying where in the atria these circuits are most likely to form, drawing on techniques from network science. Our approach indicates that a number of factors are key to determining where circuits form, most notably the thickness of the heart muscle, and the alignment of muscle fibres.


title: Network- and Enrichment-based Inference of Phenotypes and Targets from large-scale Disease Maps

abstract: Disease maps have emerged as computational knowledge bases for exploring and modeling diseasespecific molecular processes. By capturing molecular interactions, disease-associated processes, and phenotypes in standardized representations, disease maps provide a platform for applying bioinformatics and systems biology approaches. Applications range from simple map exploration to algorithm-driven target discovery and network perturbation. The web-based MINERVA environment for disease maps provides a platform to develop tools not only for mapping experimental data but also to identify, analyze and simulate disease-specific regulatory networks. We have developed a MINERVA plugin suite based on network topology and enrichment analyses that facilitate multi-omics data integration and enable in silico perturbation experiments on disease maps. We demonstrate workflows by analyzing two RNA-seq datasets on the Atlas of Inflammation Resolution (AIR). Our approach improves usability and increases the functionality of disease maps by providing easy access to available data and integration of selfgenerated data. It supports efficient and intuitive analysis of omics data, with a focus on disease maps.

category: systems biology


title: Synthetic lethality-based prediction of anti-SARS-CoV-2 targets

abstract: Novel strategies are needed to identify drug targets and treatments for the COVID-19 pandemic. The altered gene expression of virus-infected host cells provides an opportunity to specifically inhibit viral propagation via targeting the synthetic lethal (SL) partners of such altered host genes. Pursuing this antiviral strategy, here we comprehensively analyzed multiple in vitro and in vivo bulk and single-cell RNA-sequencing datasets of SARS-CoV-2 infection to predict clinically relevant candidate antiviral targets that are SL with altered host genes. The predicted SL-based targets are highly enriched for infected cell inhibiting genes reported in four SARS-CoV-2 CRISPR-Cas9 genome-wide genetic screens. Integrating our predictions with the results of these screens, we further selected a focused subset of 26 genes that we experimentally tested in a targeted siRNA screen using human Caco-2 cells. Notably, as predicted, knocking down these targets reduced viral replication and cell viability only under the infected condition without harming non-infected cells. Our results are made publicly available, to facilitate their in vivo testing and further validation.


title: Blood-based untargeted metabolomics in Relapsing-Remitting Multiple Sclerosis revealed the testable therapeutic target

abstract: Metabolic aberrations impact the pathogenesis of multiple sclerosis (MS) and possibly can provide clues for new treatment strategies. Using untargeted metabolomics, we measured serum metabolites from 35 relapsing-remitting patients and 14 healthy age-matched controls. Out of 632 known metabolites detected, 60 were significantly altered in relapsing-remitting MS (RRMS). Bioinformatics analysis identified an altered "metabotype" in RRMS patients, represented by 4 changed metabolic pathways of glycerophospholipid, citrate cycle, sphingolipid, and pyruvate metabolism. Interestingly, the common upstream metabolic pathway feeding these 4 pathways is the glycolysis pathway. Real-time bioenergetic analysis of the patient derived peripheral blood mononuclear cells, showed enhanced glycolysis, supporting the altered metabolic state of immune cells. Experimental autoimmune encephalomyelitis mice treated with the glycolytic inhibitor, 2-deoxy-D-glucose ameliorated the disease progression and inhibited the disease pathology significantly by promoting the anti-inflammatory phenotype of monocytes/macrophage in the central nervous system. Our study suggests that targeting glycolysis offers a potential target for MS.

category: neuroscience


title: Adaptation and serial choice bias are unaltered in autism

abstract: Autism Spectrum Disorder (ASD) or autism is characterized by social and non-social symptoms, including sensory hyper- and hyposensitivities. A suggestion has been put forward that some of these symptoms could be explained by differences in how sensory information is integrated with its context, including a lower tendency to leverage the past in the processing of new perceptual input. At least two history-dependent effects of opposite directions have been described in the visual perception literature: a repulsive adaptation effect, where perception of a stimulus is biased away from an adaptor stimulus, and an attractive serial choice bias, where perceptual choices are biased towards the previous choice. In this study, we investigated whether autistic participants differed in either bias from typically developing controls (TD). Sixty-four adolescent participants (31 with ASD, 33 TD) were asked to categorize oriented line stimuli in two tasks which were designed so that we would induce either adaptation or serial choice bias. Although our tasks successfully induced both biases, in comparing the two groups, we found no differences in the magnitude of adaptation nor in the modulation of perceptual choices by the previous choice. In conclusion, we find no evidence of a decreased integration of the past in visual perception of autistic individuals.


title: Structural Disconnection of the Posterior Medial Frontal Cortex Reduces Speech Error Monitoring

abstract: Optimal performance in any task relies on the ability to detect and repair errors. The anterior cingulate cortex and the broader posterior medial frontal cortex (pMFC) are active during error processing. However, it is unclear whether damage to the pMFC impairs error monitoring. We hypothesized that successful error monitoring critically relies on connections between the pMFC and broader cortical networks involved in executive functions and the task being monitored. We tested this hypothesis in the context of speech error monitoring in people with post-stroke aphasia. Diffusion weighted images were collected in 51 adults with chronic left-hemisphere stroke and 37 age-matched control participants. Whole-brain connectomes were derived using constrained spherical deconvolution and anatomically-constrained probabilistic tractography. Support vector regressions identified white matter connections in which lost integrity in stroke survivors related to reduced error detection during confrontation naming. Lesioned connections to the bilateral pMFC were related to reduce error monitoring, including many connections to regions associated with speech production and executive function. We conclude that connections to the pMFC support error monitoring. Error monitoring in speech production is supported by the structural connectivity between the pMFC and regions involved in speech production and executive function. Interactions between pMFC and other taskrelevant processors may similarly be critical for error monitoring in other task contexts.


title: Bodily ownership of an independent supernumerary limb: An exploratory study

abstract: Can our brain perceive a sense of ownership towards an independent supernumerary limb; one that can be moved independently of any other limb and provides its own independent movement feedback? Following the rubberhand illusion experiment, a plethora of studies has shown that the human representation of  self is very plastic. But previous studies have almost exclusively investigated ownership towards  substitute artificial limbs, which are controlled by the movements of a real limb and/or limbs from which non-visual sensory feedback is provided on an existing limb. Here, to investigate whether the human brain can own an independent artificial limb, we first developed a novel independent robotic  sixth finger. We allowed participants to train using the finger and examined whether it induced changes in the body representation using behavioral as well as cognitive measures. Our results suggest that unlike a substituted artificial limb (like in the rubber hand experiment), it is more difficult for humans to perceive a sense of ownership towards an independent limb. However, ownership does seem possible, as we observed clear tendencies of changes in the body representation that correlated with the cognitive reports of the sense of ownership. Our results provide the first evidence to show that an independent supernumerary limb can be embodied by the human brain.


title: Permutational immune analysis reveals architectural similarities between inflammaging, Down syndrome and autoimmunity

abstract:
People with Down syndrome show cellular and clinical features of dysregulated aging of the immune system, including naive-memory shift in the T cell compartment and increased incidence of autoimmunity. However, a quantitative understanding of how various immune compartments change with age in Down syndrome remains lacking. Here we performed deep immunophenotyping of a cohort of individuals with Down syndrome across the lifespan, selecting for individuals not affected by autoimmunity. We simultaneously interrogated age- and sex-matched healthy neurotypical controls and people with type 1 diabetes, as a representative autoimmune disease. We built a new analytical software, IMPACD, that enabled us to rapidly identify many features of immune dysregulation in Down syndrome that are recapitulated in other autoimmune diseases. We found significant quantitative and qualitative dysregulation of naive CD4+ and CD8+ T cells in Down syndrome and identified IL-6 as a candidate driver of some of these changes, thus extending the consideration of immunopathologic cytokines in Down syndrome beyond interferons. Notably, we successfully used immune cellular composition to generate three quantitative models of aging (i.e. immune clocks) trained on control subjects. All three immune clocks demonstrated significantly advanced immune aging in people with Down syndrome. Notably, one of these clocks, informed by Down syndrome-relevant biology, also showed advanced immune aging in people with type 1 diabetes. Together, our findings demonstrate a novel approach to studying immune aging in Down syndrome which may have implications in the context of other autoimmune diseases.

One Sentence SummaryPermutational analysis of immune landscape reveals advanced immune aging in people with Down syndrome and in people with type 1 diabetes.


title: Insertive vaginal sex has a profound impact on the penile immune correlates of HIV susceptibility

abstract:
The penis is the primary site of HIV acquisition in heterosexual men. Elevated penile inflammatory cytokines increase sexual acquisition risk, and topically applied cytokines enhance foreskin HIV susceptibility in an explant model. However, the impact of penile-vaginal sex on these immune parameters is undefined. Heterosexual couples were recruited to the Sex, Couples and Science (SECS) Study, with the collection of penile swabs, semen, cervico-vaginal secretions, and blood after a period of abstinence, and repeated sampling up to 72 hours after either condomless (n=30) or condom-protected (n=8) penile-vaginal sex. Soluble immune parameters were quantified by multiplex immunoassay. Co-primary immune endpoints were penile levels of IL-8 and MIG, cytokines previously linked to penile HIV acquisition. One hour after sex there were dramatic increases in penile IL-8 and MIG levels, regardless of condom use, with a gradual return to baseline by 72 hours; similar patterns were observed for other chemoattractant chemokines. Penile cytokine changes were similar in circumcised and uncircumcised men, and repeated measures ANOVA and ANCOVA models demonstrated that the degree of change after condomless sex was explained by cytokine levels in their partners cervico-vaginal secretions. This may have important implications for the biology of penile HIV acquisition.

Author summaryIn heterosexual men, the penis is the primary site of Human Immunodeficiency Virus (HIV) acquisition. Levels of inflammatory cytokines in the coronal sulcus are associated with an increased HIV risk, and we hypothesized that these may be altered after insertive penile sex. Therefore, we designed the Sex, Couples and Science Study (SECS study) to define the impact of penile-vaginal sex on the penile immune correlates of HIV susceptibility. We found that multiple coronal sulcus cytokines increased dramatically and rapidly after sex, regardless of condom use, with a return to baseline levels by 72 hours. The changes observed after condomless sex were strongly predicted by cytokine concentrations in the vaginal secretions of the female partner, and were similar in circumcised and uncircumcised men. We believe that these findings have important implications for understanding the immunopathogenesis of penile HIV acquisition; in addition, they have important implications for the design of clinical studies of penile HIV acquisition and prevention.


title: In Silico Epitope Prediction And VP1 Modelling For Foot-And-Mouth Serotype SAT 2 For Vaccine Design In East Africa

abstract: Foot and Mouth Disease Virus has seven distinct, geographically localized, serotypes and a vaccination targeting one serotype does not confer immunity against another serotype. The use of inactivated vaccines is not safe and confers an immunity with a relatively shorter time. Using the VP1 sequences isolated in East Africa, we have predicted epitopes able to induce humoral and cell-mediated immunity in cattle. The Wu-Kabat variability index calculated in this study reflects the variable, including the known GH loop, and conserved regions, with the latter being good candidates for region-tailored vaccine design. Furthermore, we modelled the identified epitopes on a 3D model (PDB ID:5aca) to represent the epitopes structurally. This study can be used for in vitro and in vivo experiments.


title: EROS-mediated control of NOX2 and P2X7 biosynthesis

abstract: EROS (Essential for Reactive Oxygen Species) protein is indispensable for expression of the gp91phox-p22phox heterodimer of the phagocyte NADPH oxidase. EROS deficiency in humans causes the severe immunodeficiency, chronic granulomatous disease (CGD), but its mechanism of action remains unknown. We elucidate the role of EROS, showing it acts at the earliest stages of gp91phox maturation. It binds the immature 58kDa gp91phox directly, interacting with the OST glycosylation machinery and prevents gp91phox degradation. EROS also regulates the purine receptors P2X7 and P2X1 through direct interactions and P2X7 is almost absent in EROS deficiency. Accordingly, lack of EROS results in markedly abnormal P2X7 signalling, inflammasome activation and T cell responses. The loss of both ROS and P2X7 signalling leads to resistance to influenza infection. Our work identifies EROS as a highly selective chaperone for key proteins in innate and adaptive immunity and a rheostat for immunity to infection. It has profound implications for our understanding of immune physiology, immunodeficiency and gene therapy.


title: Source-to-Target Automatic Rotating Estimation (STARE) -- A publicly-available, blood-free quantification approach for PET tracers with irreversible kinetics: Theoretical framework and validation for FDG

abstract:
IntroductionFull quantification of positron emission tomography (PET) data requires an input function. This generally means arterial blood sampling, which is invasive, labor-intensive and burdensome. There is no current, standardized method to fully quantify PET radiotracers with irreversible kinetics in the absence of blood data. Here, we present Source-to-Target Automatic Rotating Estimation (STARE), a novel, data-driven approach to quantify the net influx rate (Ki) of irreversible PET radiotracers, that requires only individual-level PET data and no blood data. We validate STARE with [18F]FDG PET and assess its performance using simulations.

MethodsSTARE builds upon a source-to-target tissue model, where the tracer time activity curves (TACs) in multiple "target" regions are expressed at once as a function of a "source" region, based on the two-tissue irreversible compartment model, and separates target region Ki from source Ki by fitting the source-to-target model across all target regions simultaneously. To ensure identifiability, data-driven, subject-specific anchoring is used in the STARE minimization, which takes advantage of the PET signal in a vasculature cluster in the FOV that is automatically extracted and partial volume-corrected. To avoid the need for any a priori determination of a single source region, each of the considered regions acts in turn as the source, and a final Ki is estimated in each region by averaging the estimates obtained in each source rotation.

ResultsIn a large dataset of [18F]FDG human scans (N=69), STARE Ki estimates were in good agreement with corresponding arterial blood-based estimates (regression slope=0.88, r=0.80), and were precisely estimated, as assessed by comparing STARE Ki estimates across several runs of the algorithm (coefficient of variation across runs=6.74 {+/-} 2.48%). In simulations, STARE Ki estimates were largely robust to factors that influence the individualized anchoring used within its algorithm.

ConclusionThrough simulations and application to [18F]FDG PET data, feasibility is demonstrated for STARE blood-free, data-driven quantification of Ki. Future work will include applying STARE to PET data obtained with a portable PET camera and to other irreversible radiotracers.


title: Droplet-based Single-cell Total RNA-seq Reveals Differential Non-Coding Expression and Splicing Patterns during Mouse Development

abstract: In recent years, single-cell transcriptome sequencing has revolutionized biology, allowing for the unbiased characterization of cellular subpopulations. However, most methods amplify the termini of polyadenylated transcripts capturing only a small fraction of the total cellular transcriptome. This precludes the detection of many long non-coding, short non-coding and non-polyadenylated protein-coding transcripts. Additionally, most workflows do not sequence the full transcript hindering the analysis of alternative splicing. We therefore developed VASA-seq to detect the total transcriptome in single cells. VASA-seq is compatible with both plate-based formats and droplet microfluidics. We applied VASA-seq to over 30,000 single cells in the developing mouse embryo during gastrulation and early organogenesis. The dynamics of the total single-cell transcriptome result in the discovery of novel cell type markers many based on non-coding RNA, an in vivo cell cycle analysis and an improved RNA velocity characterization. Moreover, it provides the first comprehensive analysis of alternative splicing during mammalian development.


title: IPMK physically binds to the SWI/SNF complex and modulates BRG1 occupancy

abstract: Inositol polyphosphate multikinase (IPMK), a key enzyme in the inositol polyphosphate (IP) metabolism, is a pleiotropic signaling factor involved in major biological events including transcriptional control. In yeasts, IPMK and its IP products were known to promote the activity of SWI/SNF chromatin remodeling complex, which plays a critical role in gene expression by regulating chromatin accessibility. However, the direct linkage between IPMK and chromatin remodelers remains unclear, raising a question on how IPMK contributes to the transcriptional regulation in mammals. By employing unbiased screenings and in vivo/in vitro immunoprecipitations, here we demonstrated that IPMK physically associates with native mammalian SWI/SNF complexes by directly binding to SMARCB1, BRG1, and SMARCC1. Furthermore, we identified the specific domains required for the IPMK-SMARCB1 binding. Notably, using CUT&RUN and ATAC-seq assays, we discovered that IPMK co-localizes with BRG1 and regulates BRG1 localization as well as BRG1-mediated chromatin accessibility in a genome-wide manner (including promoter-TSS) in mouse embryonic stem cells. Finally, our mRNA-seq analyses revealed that IPMK and SMARCB1 regulate common gene sets, validating a functional link between IPMK and SWI/SNF complex. Together, these findings establish an importance of IPMK in promoter targeting of the SWI/SNF complex, thereby contributing to SWI/SNF-meditated chromatin accessibility and transcription.


title: Predictors of human-infective RNA virus discovery in the United States, China and Africa, an ecological study

abstract:
BackgroundThe variation in the pathogen type as well as the spatial heterogeneity of predictors make the generality of any associations with pathogen discovery debatable. Our previous work confirmed that the association of a group of predictors differed across different types of RNA viruses, yet there have been no previous comparisons of the specific predictors for RNA virus discovery in different regions. The aim of the current study was to close the gap by investigating whether predictors of discovery rates within three regions--the United States, China and Africa--differ from one another and from those at the global level.

MethodsBased on a comprehensive list of human-infective RNA viruses, we collated published data on first discovery of each species in each region. We used a Poisson boosted regression tree (BRT) model to examine the relationship between virus discovery and 33 predictors representing climate, socio-economics, land use, and biodiversity across each region separately. The discovery probability in three regions in 2010-2019 was mapped using the fitted models and historical predictors.

ResultsThe numbers of human-infective virus species discovered in the United States, China and Africa up to 2019 were 95, 80 and 107 respectively, with China lagging behind the other two regions. In each region, discoveries were clustered in hotspots. BRT modelling suggested that in all three regions RNA virus discovery was best predicted by land use and socio- economic variables, followed by climatic variables and biodiversity, though the relative importance of these predictors varied by region. Map of virus discovery probability in 2010- 2019 indicated several new hotspots outside historical high-risk areas. Most new virus species since 2010 in each region (6/6 in the United States, 19/19 in China, 12/19 in Africa) were discovered in high risk areas as predicted by our model.

ConclusionsThe drivers of spatiotemporal variation in virus discovery rates vary in different regions of the world. Within regions virus discovery is driven mainly by land-use and socio- economic variables; climate and biodiversity variables are consistently less important predictors than at a global scale. Potential new discovery hotspots in 2010-2019 are identified. Results from the study could guide active surveillance for new human-infective viruses in local high risk areas.

FundingDarwin Trust of Edinburgh; European Union.

category: ecology


title: DNA metabarcoding reveals host-specific communities of arthropods residing in fungal fruit bodies

abstract: Biological communities within living organisms are structured by their hosts traits. How host traits affect biodiversity and community composition is poorly explored for some associations, such as arthropods within fungal fruit bodies. Using DNA metabarcoding, we revealed the arthropod communities in living fruit bodies of eleven wood-decay fungi from boreal forests and investigated how they were affected by different fungal traits. Arthropod diversity was higher in fruit bodies with a larger surface area-to-volume ratio, suggesting that colonisation is crucial to maintain arthropod populations. Diversity was not higher in long-lived fruit bodies, most likely because these fungi invest in physical or chemical defences against arthropods. Arthropod community composition was structured by all measured host traits, namely fruit body size, thickness, surface area, morphology and toughness. Notably, we identified a community gradient where soft and short-lived fruit bodies harboured more true flies, while tougher and long-lived fruit bodies had more oribatid mites and beetles, which might reflect different development times of the arthropods. Ultimately, close to 75% of the arthropods were specific to one or two fungal hosts. Besides revealing surprisingly diverse and host-specific arthropod communities within fungal fruit bodies, our study provided insight on how host traits structure communities.


title: Planktonic protist diversity across contrasting Subtropical and Subantarctic waters of the southwest Pacific

abstract:
Planktonic protists are an essential component of marine pelagic ecosystems where they mediate important trophic and biogeochemical functions. Although these functions are largely influenced by their taxonomic affiliation, the composition and spatial variability of planktonic protist communities remain poorly characterized in vast areas of the ocean. Here, we investigated the diversity of these communities in contrasting oceanographic conditions of the southwest Pacific sector (33-58{degrees}S) using DNA metabarcoding of the 18S rRNA gene. Seawater samples collected during twelve cruises (n = 482, 0-2000 m) conducted east of New Zealand were used to characterize protist communities in Subtropical (STW) and Subantarctic (SAW) water masses and the Subtropical Front (STF) that separates them. Diversity decreased with latitude and temperature but tended to be lowest in the STF. Sample ordination resulting from the abundance of amplicon single variants (ASVs) corresponded to the different water masses. Overall, Dinophyceae (34% of standardized total number of reads) and Chlorophyta (27%) co-dominated the euphotic zone, but their relative abundance and composition at class and lower taxonomic levels varied consistently between water masses. Among Chlorophyta, several picoplanktonic algae species of the Mamiellophyceae class including Ostreococcus lucimarinus dominated in STW, while the Chloropicophyceae species Chloroparvula pacifica was most abundant in SAW. Bacillariophyta (7%), Prymnesiophyceae (5%), and Pelagophyceae (3%) classes were less abundant but showed analogous water mass specificity at class and finer taxonomic levels. Protist community composition in the STF had mixed characteristics and showed regional differences with the southern STF (50{degrees}S) having more resemblance with subantarctic communities than the STF over the Chatham Rise region (42-44{degrees}S). Below the euphotic zone, Radiolaria sequences dominated the dataset (52%) followed by Dinophyceae (27%) and other heterotrophic groups like Marine Stramenopiles and ciliates (3%). Among Radiolaria, several unidentified ASVs assigned to Spumellarida were most abundant, but showed significantly different distribution between STW and SAW highlighting the need to further investigate the taxonomy and ecology of this group. This study represents a significant step forward towards characterizing protistan communities composition in relation to major water masses and fronts in the South Pacific providing new insights about the biogeography and ecological preferences of different taxa from class to species and genotypic level.

HighlightsO_LIWater-mass preference of different taxa emerged at class, species and genotypic level.
C_LIO_LIMamiellophyceae green algae dominated in subtropical waters.
C_LIO_LIDinophyceae and Chloropicophyceae green algae dominated in subantarctic waters.
C_LIO_LIA diverse assemblage of Radiolaria dominated the mesopelagic zone.
C_LIO_LISmall rather than large taxa dominated phytoplankton blooms in subtropical waters.
C_LI


title: Genome-wide association mapping within a single Arabidopsis thaliana population reveals a richer genetic architecture for defensive metabolite diversity

abstract: A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genic heterogeneity and epistasis, would be expected to increase in severity with the geographic range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (< 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits.


title: Convergent transcriptomic targets of propranolol and primidone identify potential biomarkers for essential tremor

abstract: Essential tremor (ET) is one of the most common movement disorders, affecting nearly 5% of individuals over 65 years old. Despite its high heritability, few genetic risk loci for ET have been identified. Recent advances in pharmacogenomics have generated a wealth of data that led to the identification of molecular signatures in response to hundreds of chemical compounds. Among the different forms of data, gene expression has proven to be quite successful for the inference of drug response in cell models. We sought to leverage this approach in the context of ET where many patients are responsive two drugs: propranolol and primidone. Propranolol- and primidone-specific transcriptomic drug targets, as well as convergent gene targets across both drugs, could provide insights into the pathogenesis of ET and identify possible targets of interest for future treatments. In this study, cerebellar DAOY and neural progenitor cells were treated for 5 days with clinical concentrations of propranolol and primidone, after which RNA-sequencing was used to identify differentially expressed genes. The expression of genes previously implicated in genetic and transcriptomic studies of ET and other movement disorders, such as TRAPPC11, were significantly upregulated by propranolol. Pathway enrichment analysis identified multiple terms related to calcium signalling, endosomal sorting, axon guidance, and neuronal morphology. Convergent differentially expressed genes across all treatments and cell types were also found to be significantly more mutationally constrained, implying that they might harbour rare deleterious variants implicated in disease. Furthermore, these genes were enriched within cell types having high expression of ET related genes in both cortical and cerebellar tissues. Altogether, our results highlight potential cellular and molecular mechanisms associated with tremor reduction and identify relevant genetic biomarkers for drug-responsiveness in ET.

category: genomics


title: Behavioural measurements of hearing-sensitivity in swarming Anopheles coluzzii mosquitoes reveal an unprecedentedly low threshold

abstract: Mosquitoes of many species mate in station-keeping swarms. Mating chases ensue as soon as a male detects the flight tones of a female with his auditory organs. Previous studies of hearing thresholds have mainly used electrophysiological methods that prevent the mosquito from flying naturally. The main aim of this study was to quantify behaviourally the sound-level threshold at which males can hear females. Free-flying male Anopheles coluzzii were released in a large arena (~2 m high x 2 m x 1 m) with a conspicuous object on the ground that stimulates swarming behaviour. Males were exposed to a range of natural and synthetic played-back sounds of female flight. We monitored the responses of males and their distance to the speaker by recording changes in their wingbeat frequency and angular speed. We show that the mean male behavioural threshold of particle-velocity hearing lies between 13-20 dB SVL (95%-CI). A conservative estimate of 20 dB SVL (i.e., < 0.5 m/s particle velocity) is already 12 to 26 dB lower than most of the published electrophysiological measurements from the Johnstons organ. In addition, we suggest that 1) the first harmonic of female flight-sound is sufficient for males to detect her presence, 2) males respond with a greater amplitude to single-female sounds than to the sound of a group of females and 3) the response of males to the playback of the flight sound of a live female is the same as that of a recorded sound of constant frequency and amplitude.

category: animal behavior and cognition


title: Behavioural analysis of swarming mosquitoes reveals higher hearing sensitivity than previously measured with electrophysiology methods


title: Behavior in mice subjected to a Token Slot-machine: effect of Unpredictable Chronic Mild Stress

abstract:
Several studies have succeeded in teaching animals (primates, pigeons, rats, but not mice) the value of tokens by having them executing a task using a vending-machine apparatus, where in order to receive the primary reinforcement (food), the animals had to perform a specific action that allowed them to obtain the secondary reinforcement (tokens: metal balls). We tried to assess this kind of behavior in mice that had previously been trained to use some tokens, with the aim of rewarding them not with food, but with other tokens, as a result of a token economy task. We found that mice exhibit economic behavior. Further on, our research tried to investigate the effect of stress on their operant decision-making. Therefore, the mice were divided into two groups: a Control group (n=10) and a group subjected to an Unpredictable Chronic Mild Stress (UCMS) treatment (n=8). We found that chronic stress increases some aspects of sub-optimal economic activity.

Summary statementWe designed an original model enabling to assess behavior in mice that had previously been trained to use some tokens, with the aim of rewarding them not with food, but with other tokens. Further on, our research investigated the effect of stress on their operant decision-making.