Optimal transport analysis reveals trajectories in steady-state systems

Understanding how cells change their identity and behaviour in living systems is an important question in many fields of biology. The problem of inferring cell trajectories from single-cell measurements has been a major topic in the single-cell analysis community, with different methods developed for equilibrium and non-equilibrium systems (e.g. haematopoeisis vs. embryonic development). We show that optimal transport analysis, a technique originally designed for analysing time-courses, may also be applied to infer cellular trajectories from a single snapshot of a population in equilibrium. Therefore, optimal transport provides a unified approach to inferring trajectories that is applicable to both stationary and non-stationary systems. Our method, StationaryOT, is mathematically motivated in a natural way from the hypothesis of a Waddington's epigenetic landscape. We implement StationaryOT as a software package and demonstrate its efficacy in applications to simulated data as well as single-cell data from Arabidopsis thaliana root development.

Models of heterogeneous dopamine signaling in an insect learning and memory center

The Drosophila mushroom body exhibits dopamine dependent synaptic plasticity that underlies the acquisition of associative memories. Recordings of dopamine neurons in this system have identified signals related to external reinforcement such as reward and punishment. However, other factors including locomotion, novelty, reward expectation, and internal state have also recently been shown to modulate dopamine neurons. This heterogeneity is at odds with typical modeling approaches in which these neurons are assumed to encode a global, scalar error signal. How is dopamine dependent plasticity coordinated in the presence of such heterogeneity? We develop a modeling approach that infers a pattern of dopamine activity sufficient to solve defined behavioral tasks, given architectural constraints informed by knowledge of mushroom body circuitry. Model dopamine neurons exhibit diverse tuning to task parameters while nonetheless producing coherent learned behaviors. Notably, reward prediction error emerges as a mode of population activity distributed across these neurons. Our results provide a mechanistic framework that accounts for the heterogeneity of dopamine activity during learning and behavior.

Dopamine neurons across the animal kingdom are involved in the formation of associative memories. While numerous studies have recorded activity in these neurons related to external and predicted rewards, the diversity of these neurons’ activity and their tuning to non-reward-related quantities such as novelty, movement, and internal state have proved challenging to account for in traditional modeling approaches. Using a well-characterized model system for learning and memory, the mushroom body of Drosophila fruit flies, Jiang and Litwin-Kumar provide an account of the diversity of signals across dopamine neurons. They show that models optimized to solve tasks like those encountered by flies exhibit heterogeneous activity across dopamine neurons, but nonetheless this activity is sufficient for the system to solve the tasks. The models will be useful to generate testable hypotheses about dopamine neuron activity across different experimental conditions.

Targeting fibroblast growth factor receptors to combat aggressive ependymoma

Ependymomas (EPN) are central nervous system tumors comprising both aggressive and more benign molecular subtypes. However, therapy of the high-risk subtypes posterior fossa group A (PF-A) and supratentorial RELA-fusion positive (ST-RELA) is limited to gross total resection and radiotherapy, as effective systemic treatment concepts are still lacking. We have recently described fibroblast growth factor receptors 1 and 3 (FGFR1/FGFR3) as oncogenic drivers of EPN. However, the underlying molecular mechanisms and their potential as therapeutic targets have not yet been investigated in detail. Making use of transcriptomic data across 467 EPN tissues, we found that FGFR1 and FGFR3 were both widely expressed across all molecular groups. FGFR3 mRNA levels were enriched in ST-RELA showing the highest expression among EPN as well as other brain tumors. We further identified high expression levels of fibroblast growth factor 1 and 2 (FGF1, FGF2) across all EPN subtypes while FGF9 was elevated in ST-EPN. Interrogation of our EPN single-cell RNA-sequencing data revealed that FGFR3 was further enriched in cycling and progenitor-like cell populations. Corroboratively, we found FGFR3 to be predominantly expressed in radial glia cells in both mouse embryonal and human brain datasets. Moreover, we detected alternative splicing of the FGFR1/3-IIIc variant, which is known to enhance ligand affinity and FGFR signaling. Dominant-negative interruption of FGFR1/3 activation in PF-A and ST-RELA cell models demonstrated inhibition of key oncogenic pathways leading to reduced cell growth and stem cell characteristics. To explore the feasibility of therapeutically targeting FGFR, we tested a panel of FGFR inhibitors in 12 patient-derived EPN cell models revealing sensitivity in the low-micromolar to nano-molar range. Finally, we gain the first clinical evidence for the activity of the FGFR inhibitor nintedanib in the treatment of a patient with recurrent ST-RELA. Together, these preclinical and clinical data suggest FGFR inhibition as a novel and feasible approach to combat aggressive EPN.

A recurrent, homozygous EMC10 frameshift variant is associated with a syndrome of developmental delay with variable seizures and dysmorphic features

PURPOSE

The endoplasmic reticulum membrane complex (EMC) is a highly conserved, multifunctional 10-protein complex related to membrane protein biology. In seven families, we identified 13 individuals with highly overlapping phenotypes who harbor a single identical homozygous frameshift variant in EMC10.

METHODS

Using exome, genome, and Sanger sequencing, a recurrent frameshift EMC10 variant was identified in affected individuals in an international cohort of consanguineous families. Multiple families were independently identified and connected via Matchmaker Exchange and internal databases. We assessed the effect of the frameshift variant on EMC10 RNA and protein expression and evaluated EMC10 expression in normal human brain tissue using immunohistochemistry.

RESULTS

A homozygous variant EMC10 c.287delG (Refseq NM_206538.3, p.Gly96Alafs*9) segregated with affected individuals in each family, who exhibited a phenotypic spectrum of intellectual disability (ID) and global developmental delay (GDD), variable seizures and variable dysmorphic features (elongated face, curly hair, cubitus valgus, and arachnodactyly). The variant arose on two founder haplotypes and results in significantly reduced EMC10 RNA expression and an unstable truncated EMC10 protein.

CONCLUSION

We propose that a homozygous loss-of-function variant in EMC10 causes a novel syndromic neurodevelopmental phenotype. Remarkably, the recurrent variant is likely the result of a hypermutable site and arose on distinct founder haplotypes.

CD300lf is the primary physiologic receptor of murine norovirus but not human norovirus

Murine norovirus (MNoV) is an important model of human norovirus (HNoV) and mucosal virus infection more broadly. Viral receptor utilization is a major determinant of cell tropism, host range, and pathogenesis. The bona fide receptor for HNoV is unknown. Recently, we identified CD300lf as a proteinaceous receptor for MNoV. Interestingly, its paralogue CD300ld was also sufficient for MNoV infection in vitro. Here we explored whether CD300lf is the sole physiologic receptor in vivo and whether HNoV can use a CD300 ortholog as an entry receptor. We report that both CD300ld and CD300lf are sufficient for infection by diverse MNoV strains in vitro. We further demonstrate that CD300lf is essential for both oral and parenteral MNoV infection and to elicit anti-MNoV humoral responses in vivo. In mice deficient in STAT1 signaling, CD300lf is required for MNoV-induced lethality. Finally, we demonstrate that human CD300lf (huCD300lf) is not essential for HNoV infection, nor does huCD300lf inhibit binding of HNoV virus-like particles to glycans. Thus, we report huCD300lf is not a receptor for HNoV.

Relaxed constraint and functional divergence of the progesterone receptor (PGR) in the human stem-lineage

The steroid hormone progesterone, acting through the progesterone receptor (PR), a ligand-activated DNA-binding transcription factor, plays an essential role in regulating nearly every aspect of female reproductive biology. While many reproductive traits regulated by PR are conserved in mammals, Catarrhine primates evolved several derived traits including spontaneous decidualization, menstruation, and a divergent (and unknown) parturition signal, suggesting that PR may also have evolved divergent functions in Catarrhines. There is conflicting evidence, however, whether the progesterone receptor gene (PGR) was positively selected in the human lineage. Here we show that PGR evolved rapidly in the human stem-lineage (as well as other Catarrhine primates), which likely reflects an episode of relaxed selection intensity rather than positive selection. Coincident with the episode of relaxed selection intensity, ancestral sequence resurrection and functional tests indicate that the major human PR isoforms (PR-A and PR-B) evolved divergent functions in the human stem-lineage. These results suggest that the regulation of progesterone signaling by PR-A and PR-B may also have diverged in the human lineage and that non-human animal models of progesterone signaling may not faithfully recapitulate human biology.

Surfactant-enhanced DNA accessibility to nuclease accelerates phenotypic β-lactam antibiotic susceptibility testing of Neisseria gonorrhoeae

Rapid antibiotic susceptibility testing (AST) for Neisseria gonorrhoeae (Ng) is critically needed to counter widespread antibiotic resistance. Detection of nucleic acids in genotypic AST can be rapid, but it has not been successful for β-lactams (the largest antibiotic class used to treat Ng). Rapid phenotypic AST for Ng is challenged by the pathogen's slow doubling time and the lack of methods to quickly quantify the pathogen's response to β-lactams. Here, we asked two questions: (1) Is it possible to use nucleic acid quantification to measure the β-lactam susceptibility phenotype of Ng very rapidly, using antibiotic-exposure times much shorter than the 1- to 2-h doubling time of Ng? (2) Would such short-term antibiotic exposures predict the antibiotic resistance profile of Ng measured by plate growth assays over multiple days? To answer these questions, we devised an innovative approach for performing a rapid phenotypic AST that measures DNA accessibility to exogenous nucleases after exposure to β-lactams (termed nuclease-accessibility AST [nuc-aAST]). We showed that DNA in antibiotic-susceptible cells has increased accessibility upon exposure to β-lactams and that a judiciously chosen surfactant permeabilized the outer membrane and enhanced this effect. We tested penicillin, cefixime, and ceftriaxone and found good agreement between the results of the nuc-aAST after 15-30 min of antibiotic exposure and the results of the gold-standard culture-based AST measured over days. These results provide a new pathway toward developing a critically needed phenotypic AST for Ng and additional global-health threats.

Modulation of flagellar rotation in surface-attached bacteria: A pathway for rapid surface-sensing after flagellar attachment

Attachment is a necessary first step in bacterial commitment to surface-associated behaviors that include colonization, biofilm formation, and host-directed virulence. The Gram-negative opportunistic pathogen Pseudomonas aeruginosa can initially attach to surfaces via its single polar flagellum. Although many bacteria quickly detach, some become irreversibly attached and express surface-associated structures, such as Type IV pili, and behaviors, including twitching motility and biofilm initiation. P. aeruginosa that lack the GTPase FlhF assemble a randomly placed flagellum that is motile; however, we observed that these mutant bacteria show defects in biofilm formation comparable to those seen for non-motile, aflagellate bacteria. This phenotype was associated with altered behavior of ΔflhF bacteria immediately following surface-attachment. Forward and reverse genetic screens led to the discovery that FlhF interacts with FimV to control flagellar rotation at a surface, and implicated cAMP signaling in this pathway. Although cAMP controls many transcriptional programs in P. aeruginosa, known targets of this second messenger were not required to modulate flagellar rotation in surface-attached bacteria. Instead, alterations in switching behavior of the motor appeared to result from direct or indirect effects of cAMP on switch complex proteins and/or the stators associated with them.

Registration of published randomized trials: a systematic review and meta-analysis

BACKGROUND

Prospective trial registration is a powerful tool to prevent reporting bias. We aimed to determine the extent to which published randomized controlled trials (RCTs) were registered and registered prospectively.

METHODS

We searched MEDLINE and EMBASE from January 2005 to October 2017; we also screened all articles cited by or citing included and excluded studies, and the reference lists of related reviews. We included studies that examined published RCTs and evaluated their registration status, regardless of medical specialty or language. We excluded studies that assessed RCT registration status only through mention of registration in the published RCT, without searching registries or contacting the trial investigators. Two independent reviewers blinded to the other's work performed the selection. Following PRISMA guidelines, two investigators independently extracted data, with discrepancies resolved by consensus. We calculated pooled proportions and 95% confidence intervals using random-effects models.

RESULTS

We analyzed 40 studies examining 8773 RCTs across a wide range of clinical specialties. The pooled proportion of registered RCTs was 53% (95% confidence interval 44% to 58%), with considerable between-study heterogeneity. A subset of 24 studies reported data on prospective registration across 5529 RCTs. The pooled proportion of prospectively registered RCTs was 20% (95% confidence interval 15% to 25%). Subgroup analyses showed that registration was higher for industry-supported and larger RCTs. A meta-regression analysis across 19 studies (5144 RCTs) showed that the proportion of registered trials significantly increased over time, with a mean proportion increase of 27%, from 25 to 52%, between 2005 and 2015.

CONCLUSIONS

The prevalence of trial registration has increased over time, but only one in five published RCTs is prospectively registered, undermining the validity and integrity of biomedical research.

A phase II study of the dual mTOR inhibitor MLN0128 in patients with metastatic castration resistant prostate cancer

Background MLN0128 is a first-in-class, dual mTOR inhibitor with potential to outperform standard rapalogs through inhibition of TORC1 and TORC2. This phase II study was designed to assess antitumor activity of MLN0128 in metastatic castration-resistant prostate cancer (mCRPC). Methods Eligible patients had mCRPC previously treated with abiraterone acetate and/or enzalutamide. Five patients started MLN0128 at 5 mg once daily, subsequently dose reduced to 4 mg because of toxicity. Four subsequent patients started MLN0128 at 4 mg daily. Primary endpoint was progression-free survival at 6 months. Results Nine patients were enrolled and median time on treatment was 11 weeks (range: 3-30). Best response was stable disease. All patients had a rise in PSA on treatment, with a median 159% increase from baseline (range: 12-620%). Median baseline circulating tumor cell count was 1 cell/mL (range: 0-40); none had a decrease in cell count posttreatment. Grade ≤ 2 adverse events included fatigue, anorexia, and rash. The most common serious adverse events were grade 3 dyspnea and maculopapular rash. Eight patients discontinued treatment early because of radiographic progression (n = 1), grade 3 toxicity (n = 5), or investigator discretion (n = 2). Four patients had immediate PSA decline following drug discontinuation, suggesting MLN0128 could cause compensatory increase of androgen receptor (AR) activity. Correlative studies of pretreatment and posttreatment biopsy specimens revealed limited inhibition of AKT phosphorylation, 4EBP1 phosphorylation, and eIF4E activity. Conclusions Clinical efficacy of MLN0128 in mCRPC was limited likely due to dose reductions secondary to toxicity, PSA kinetics suggesting AR activation resulting from mTOR inhibition, and poor inhibition of mTOR signaling targets.

Rift valley fever viral load correlates with the human inflammatory response and coagulation pathway abnormalities in humans with hemorrhagic manifestations

Rift Valley fever virus is an arbovirus that affects both livestock and humans throughout Africa and in the Middle East. Despite its endemicity throughout Africa, it is a rare event to identify an infected individual during the acute phase of the disease and an even rarer event to collect serial blood samples from the affected patient. Severely affected patients can present with hemorrhagic manifestations of disease. In this study we identified three Ugandan men with RVFV disease that was accompanied by hemorrhagic manifestations. Serial blood samples from these men were analyzed for a series of biomarkers specific for various aspects of human pathophysiology including inflammation, endothelial function and coagulopathy. There were significant differences between biomarker levels in controls and cases both early during the illness and after clearance of viremia. Positive correlation of viral load with markers of inflammation (IP-10, CRP, Eotaxin, MCP-2 and Granzyme B), markers of fibrinolysis (tPA and D-dimer), and markers of endothelial function (sICAM-1) were all noted. However, and perhaps most interesting given the fact that these individuals exhibited hemorrhagic manifestations of disease, was the finding of a negative correlation between viral load and P-selectin, ADAMTS13, and fibrinogen all of which are associated with coagulation pathways occurring on the endothelial surface.

Phenotypically silent Cre recombination within the postnatal ventricular conduction system

The cardiac conduction system (CCS) is composed of specialized cardiomyocytes that initiate and maintain cardiac rhythm. Any perturbation to the normal sequence of electrical events within the heart can result in cardiac arrhythmias. To understand how cardiac rhythm is established at the molecular level, several genetically modified mouse lines expressing Cre recombinase within specific CCS compartments have been created. In general, Cre driver lines have been generated either by homologous recombination of Cre into an endogenous locus or Cre expression driven by a randomly inserted transgene. However, haploinsufficiency of the endogenous gene compromises the former approach, while position effects negatively impact the latter. To address these limitations, we generated a Cre driver line for the ventricular conduction system (VCS) that preserves endogenous gene expression by targeting the Contactin2 (Cntn2) 3’ untranslated region (3’UTR). Here we show that Cntn2^{3’UTR-IRES-Cre-EGFP/+} mice recombine floxed alleles within the VCS and that Cre expression faithfully recapitulates the spatial distribution of Cntn2 within the heart. We further demonstrate that Cre expression initiates after birth with preservation of native Cntn2 protein. Finally, we show that Cntn2^{3’UTR-IRES-Cre-EGFP/+} mice maintain normal cardiac mechanical and electrical function. Taken together, our results establish a novel VCS-specific Cre driver line without the adverse consequences of haploinsufficiency or position effects. We expect that our new mouse line will add to the accumulating toolkit of CCS-specific mouse reagents and aid characterization of the cell-autonomous molecular circuitry that drives VCS maintenance and function.

High-Throughput Accurate Single-Cell Screening of Euglena gracilis with Fluorescence-Assisted Optofluidic Time-Stretch Microscopy

The development of reliable, sustainable, and economical sources of alternative fuels is an important, but challenging goal for the world. As an alternative to liquid fossil fuels, algal biofuel is expected to play a key role in alleviating global warming since algae absorb atmospheric CO2 via photosynthesis. Among various algae for fuel production, Euglena gracilis is an attractive microalgal species as it is known to produce wax ester (good for biodiesel and aviation fuel) within lipid droplets. To date, while there exist many techniques for inducing microalgal cells to produce and accumulate lipid with high efficiency, few analytical methods are available for characterizing a population of such lipid-accumulated microalgae including E. gracilis with high throughout, high accuracy, and single-cell resolution simultaneously. Here we demonstrate high-throughput, high-accuracy, single-cell screening of E. gracilis with fluorescence-assisted optofluidic time-stretch microscopy-a method that combines the strengths of microfluidic cell focusing, optical time-stretch microscopy, and fluorescence detection used in conventional flow cytometry. Specifically, our fluorescence-assisted optofluidic time-stretch microscope consists of an optical time-stretch microscope and a fluorescence analyzer on top of a hydrodynamically focusing microfluidic device and can detect fluorescence from every E. gracilis cell in a population and simultaneously obtain its image with a high throughput of 10,000 cells/s. With the multi-dimensional information acquired by the system, we classify nitrogen-sufficient (ordinary) and nitrogen-deficient (lipid-accumulated) E. gracilis cells with a low false positive rate of 1.0%. This method holds promise for evaluating cultivation techniques and selective breeding for microalgae-based biofuel production.

Dynamic ASXL1 Exon Skipping and Alternative Circular Splicing in Single Human Cells

Circular RNAs comprise a poorly understood new class of noncoding RNA. In this study, we used a combination of targeted deletion, high-resolution splicing detection, and single-cell sequencing to deeply probe ASXL1 circular splicing. We found that efficient circular splicing required the canonical transcriptional start site and inverted AluSx elements. Sequencing-based interrogation of isoforms after ASXL1 overexpression identified promiscuous linear splicing between all exons, with the two most abundant non-canonical linear products skipping the exons that produced the circular isoforms. Single-cell sequencing revealed a strong preference for either the linear or circular ASXL1 isoforms in each cell, and found the predominant exon skipping product is frequently co-expressed with its reciprocal circular isoform. Finally, absolute quantification of ASXL1 isoforms confirmed our findings and suggests that standard methods overestimate circRNA abundance. Taken together, these data reveal a dynamic new view of circRNA genesis, providing additional framework for studying their roles in cellular biology.

Computational Investigation of the Interplay of Substrate Positioning and Reactivity in Catechol O-Methyltransferase

Catechol O-methyltransferase (COMT) is a SAM- and Mg2+-dependent methyltransferase that regulates neurotransmitters through methylation. Simulations and experiments have identified divergent catecholamine substrate orientations in the COMT active site: molecular dynamics simulations have favored a monodentate coordination of catecholate substrates to the active site Mg2+, and crystal structures instead preserve bidentate coordination along with short (2.65 Å) methyl donor-acceptor distances. We carry out longer dynamics (up to 350 ns) to quantify interconversion between bidentate and monodentate binding poses. We provide a systematic determination of the relative free energy of the monodentate and bidentate structures in order to identify whether structural differences alter the nature of the methyl transfer mechanism and source of enzymatic rate enhancement. We demonstrate that the bidentate and monodentate binding modes are close in energy but separated by a 7 kcal/mol free energy barrier. Analysis of interactions in the two binding modes reveals that the driving force for monodentate catecholate orientations in classical molecular dynamics simulations is derived from stronger electrostatic stabilization afforded by alternate Mg2+ coordination with strongly charged active site carboxylates. Mixed semi-empirical-classical (SQM/MM) substrate C-O distances (2.7 Å) for the bidentate case are in excellent agreement with COMT X-ray crystal structures, as long as charge transfer between the substrates, Mg2+, and surrounding ligands is permitted. SQM/MM free energy barriers for methyl transfer from bidentate and monodentate catecholate configurations are comparable at around 21-22 kcal/mol, in good agreement with experiment (18-19 kcal/mol). Overall, the work suggests that both binding poses are viable for methyl transfer, and accurate descriptions of charge transfer and electrostatics are needed to provide balanced relative barriers when multiple binding poses are accessible, for example in other transferases.

The Role of MAPT Haplotype H2 and Isoform 1N/4R in Parkinsonism of Older Adults

Background and Objective

Recently, we have shown that the Parkinson’s disease (PD) susceptibility locus MAPT (microtubule associated protein tau) is associated with parkinsonism in older adults without a clinical diagnosis of PD. In this study, we investigated the relationship between parkinsonian signs and MAPT transcripts by assessing the effect of MAPT haplotypes on alternative splicing and expression levels of the most common isoforms in two prospective clinicopathologic studies of aging.

Materials and Methods

using regression analysis, controlling for age, sex, study and neuropathology, we evaluated 976 subjects with clinical, genotyping and brain pathology data for haplotype analysis. For transcript analysis, we obtained MAPT gene and isoform-level expression from the dorsolateral prefrontal cortex for 505 of these subjects.

Results

The MAPT H2 haplotype was associated with lower total MAPT expression (p = 1.2x10-14) and global parkinsonism at both study entry (p = 0.001) and proximate to death (p = 0.050). Specifically, haplotype H2 was primarily associated with bradykinesia in both assessments (p<0.001 and p = 0.008). MAPT total expression was associated with age and decreases linearly with advancing age (p<0.001). Analysing MAPT alternative splicing, the expression of 1N/4R isoform was inversely associated with global parkinsonism (p = 0.008) and bradykinesia (p = 0.008). Diminished 1N/4R isoform expression was also associated with H2 (p = 0.001).

Conclusions

Overall, our results suggest that age and H2 are associated with higher parkinsonism score and decreased total MAPT RNA expression. Additionally, we found that H2 and parkinsonism are associated with altered expression levels of specific isoforms. These findings may contribute to the understanding of the association between MAPT locus and parkinsonism in elderly subjects and in some extent to age-related neurodegenerative diseases.

Human Epithelial Cells Discriminate between Commensal and Pathogenic Interactions with Candida albicans

The commensal fungus, Candida albicans, can cause life-threatening infections in at risk individuals. C. albicans colonizes mucosal surfaces of most people, adhering to and interacting with epithelial cells. At low concentrations, C. albicans is not pathogenic nor does it cause epithelial cell damage in vitro; at high concentrations, C. albicans causes mucosal infections and kills epithelial cells in vitro. Here we show that while there are quantitative dose-dependent differences in exposed epithelial cell populations, these reflect a fundamental qualitative difference in host cell response to C. albicans. Using transcriptional profiling experiments and real time PCR, we found that wild-type C. albicans induce dose-dependent responses from a FaDu epithelial cell line. However, real time PCR and Western blot analysis using a high dose of various C. albicans strains demonstrated that these dose-dependent responses are associated with ability to promote host cell damage. Our studies support the idea that epithelial cells play a key role in the immune system by monitoring the microbial community at mucosal surfaces and initiating defensive responses when this community is dysfunctional. This places epithelial cells at a pivotal position in the interaction with C. albicans as epithelial cells themselves promote C. albicans stimulated damage.