Publications

Assays detecting blood transcriptome changes are studied for infectious disease diagnosis. Blood-based RNA alternative splicing (AS) events, which have not been well characterized in pathogen infection, have potential normalization and assay platform stability advantages over gene expression for diagnosis. Here, we present a computational framework for developing AS diagnostic biomarkers. Leveraging a large prospective cohort of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and whole-blood RNA sequencing (RNA-seq) data, we identify a major functional AS program switch upon viral infection. Using an independent cohort, we demonstrate the improved accuracy of AS biomarkers for SARS-CoV-2 diagnosis compared with six reported transcriptome signatures. We then optimize a subset of AS-based biomarkers to develop microfluidic PCR diagnostic assays. This assay achieves nearly perfect test accuracy (61/62 = 98.4 percent) using a naive principal component classifier, significantly more accurate than a gene expression PCR assay in the same cohort. Therefore, our RNA splicing computational framework enables a promising avenue for host-response diagnosis of infection.

From insects to mammals, oocytes and sperm develop within germline cysts comprising cells connected by intercellular bridges (ICBs). In numerous insects, formation of the cyst is accompanied by growth of the fusome—a membranous organelle that permeates the cyst. Fusome composition and function are best understood in Drosophila melanogaster: during oogenesis, the fusome dictates cyst topology and size and facilitates oocyte selection, while during spermatogenesis, the fusome synchronizes the cyst’s response to DNA damage. Despite its distinct and sex-specific roles during insect gametogenesis, elucidating fusome growth and inheritance in females and its structure and connectivity in males has remained challenging. Here, we take advantage of advances in three-dimensional (3D) confocal microscopy and computational image processing tools to reconstruct the topology, growth, and distribution of the fusome in both sexes. In females, our experimental findings inform a theoretical model for fusome assembly and inheritance and suggest that oocyte selection proceeds through an ‘equivalency with a bias’ mechanism. In males, we find that cell divisions can deviate from the maximally branched pattern observed in females, leading to greater topological variability. Our work consolidates existing disjointed experimental observations and contributes a readily generalizable computational approach for quantitative studies of gametogenesis within and across species.

Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens.

The aggregation of epitopes that are also able to bind major histocompatibility complex (MHC) alleles raises questions around the potential connection between the formation of epitope aggregates and their affinities to MHC receptors. We first performed a general bioinformatic assessment over a public dataset of MHC class II epitopes, finding that higher experimental binding correlates with higher aggregation-propensity predictors. We then focused on the case of P10, an epitope used as a vaccine candidate against Paracoccidioides brasiliensis that aggregates into amyloid fibrils. We used a computational protocol to design variants of the P10 epitope to study the connection between the binding stabilities towards human MHC class II alleles and their aggregation propensities. The binding of the designed variants was tested experimentally, as well as their aggregation capacity. High-affinity MHC class II binders in vitro were more disposed to aggregate forming amyloid fibrils capable of binding Thioflavin T and congo red, while low affinity MHC class II binders remained soluble or formed rare amorphous aggregates. This study shows a possible connection between the aggregation propensity of an epitope and its affinity for the MHC class II cleft.

Omicron BA.1 is a highly infectious variant of SARS-CoV-2 that carries more than thirty mutations on the spike protein in comparison to the Wuhan wild type (WT). Some of the Omicron mutations, located on the receptor-binding domain (RBD), are exposed to the surrounding solvent and are known to help evade immunity. However, the impact of buried mutations on the RBD conformations and on the mechanics of the spike opening is less evident. Here, we use all-atom molecular dynamics (MD) simulations with metadynamics to characterize the thermodynamic RBD-opening ensemble, identifying significant differences between WT and Omicron. Specifically, the Omicron mutations S371L, S373P, and S375F make more RBD interdomain contacts during the spike's opening. Moreover, Omicron takes longer to reach the transition state than WT. It stabilizes up-state conformations with fewer RBD epitopes exposed to the solvent, potentially favoring immune or antibody evasion.

Multicellular organisms contain diverse tissues built from multiple cell types. It remains unclear how large numbers of interacting cells can precisely coordinate their gene expression during tissue self-organization. We develop a generalized model of multicellular gene expression that includes intracellular and intercellular gene interactions in tissue-like collectives. Motivated by modern transcriptomics, we represent multistable cellular phenotypes by mapping the binarized transcriptional patterns of individual cells onto Hopfield networks. We incorporate spatial cell-cell signaling by coupling transcriptional states of adjacent cells on a square lattice. We show that tuning the intercellular signaling strength results in a cascade of transitions toward different collective states with emergent single-cell phenotypes. Despite an enormous number of possible tissue states, we find that intercellular signaling tends to stabilize a small number of compositionally and spatially simple tissue types. These results establish a theoretical framework to investigate how cell collectives self-organize into distinct stable patterns.

The F-BAR protein Cdc15 is essential for cytokinesis in Schizosaccharomyces pombe and plays a key role in attaching the cytokinetic ring (CR) to the plasma membrane (PM). Cdc15’s abilities to bind to the membrane and oligomerize via its F-BAR domain are inhibited by phosphorylation of its intrinsically disordered region (IDR). Multiple cell polarity kinases regulate Cdc15 IDR phosphostate, and of these the DYRK kinase Pom1 phosphorylation sites on Cdc15 have been shown in vivo to prevent CR formation at cell tips. Here, we compared the ability of Pom1 to control Cdc15 phosphostate and cortical localization to that of other Cdc15 kinases: Kin1, Pck1, and Shk1. We identified distinct but overlapping cohorts of Cdc15 phosphorylation sites targeted by each kinase, and the number of sites correlated with each kinases’ abilities to influence Cdc15 PM localization. Coarse-grained simulations predicted that cumulative IDR phosphorylation moves the IDRs of a dimer apart and toward the F-BAR tips. Further, simulations indicated that the overall negative charge of phosphorylation masks positively charged amino acids necessary for F-BAR oligomerization and membrane interaction. Finally, simulations suggested that dephosphorylated Cdc15 undergoes phase separation driven by IDR interactions. Indeed, dephosphorylated but not phosphorylated Cdc15 undergoes liquid–liquid phase separation to form droplets in vitro that recruit Cdc15 binding partners. In cells, Cdc15 phosphomutants also formed PM-bound condensates that recruit other CR components. Together, we propose that a threshold of Cdc15 phosphorylation by assorted kinases prevents Cdc15 condensation on the PM and antagonizes CR assembly.

Segmentation of the Drosophila embryo is initiated by localized maternal signals. In this context, anteriorly localized Bicoid activates the gap genes in the anterior half of the embryo while posteriorly localized Nanos represses the translation of maternal hunchback mRNA to pattern the posterior half. The non-segmented termini are patterned by the localized activation of mitogen-activated protein kinase. Yet, the spatial extent of the terminal patterning system in regulating gap genes beyond poles remains unknown. We investigated the patterning potential of the terminal system using mutagenized embryos that lack both the anterior and the posterior maternal signaling systems. Using a combination of quantitative imaging and mathematical modeling, we analyzed the spatial patterns of gap genes in the early Drosophila embryo. We found that this mutant embryo develops symmetric cuticle patterns along the anteroposterior axis with two segments on each side. Notably, the terminal system can affect the expression of Krüppel in the torso region. Our mathematical model recapitulates the experimental data and reveals the potential bistability in the terminal patterning system. Collectively, our study suggests that the terminal system can act as a long-range inductive signal and establish multiple gene expression boundaries along the anteroposterior axis of the developing embryo.

The previously reported Q is a thermoresponsive coiled-coil protein capable of higher-order supramolecular assembly into fibers and hydrogels with upper critical solution temperature (UCST) behavior. Here, we introduce a new coiled-coil protein that is redesigned to disfavor lateral growth of its fibers and thus achieve a higher crosslinking density within the formed hydrogel. We also introduce a favorable hydrophobic mutation to the pore of the coiled-coil domain for increased thermostability of the protein. We note that an increase in storage modulus of the hydrogel and crosslinking density is coupled with a decrease in fiber diameter. We further fully characterize our α-helical coiled-coil (Q2) hydrogel for its structure, nano-assembly, and rheology relative to our previous single domain protein, Q, over the time of its gelation demonstrating the nature of our hydrogel self-assembly system. In this vein, we also characterize the ability of Q2 to encapsulate the small hydrophobic small molecule, curcumin, and its impact on the mechanical properties of Q2. The design parameters here not only show the importance of electrostatic potential in self-assembly but also provide a step towards predictable design of electrostatic protein interactions.

Drosophila oogenesis is a powerful and tractable model for studies of cell and developmental biology due to the multitude of well-characterized events in both germline and somatic cells, the ease of genetic manipulation in fruit flies, and the large number of egg chambers produced by each fly. Recent improvements in live imaging and ex vivo culturing protocols have enabled researchers to conduct more detailed, longer-term studies of egg chamber development, enabling insights into fundamental biological processes. Here, we present a protocol for dissection, culturing, and imaging of late-stage egg chambers to study intercellular and directional cytoplasmic flow during “nurse cell dumping.” This critical developmental process towards the latter stages of oogenesis (stages 10b/11) results in rapid growth of the oocyte and shrinkage of the nurse cells and is accompanied by dynamic changes in cell shape.