ASU Scholarship Showcase

Linnorm is a novel normalization and transformation method for the analysis of single cell RNA sequencing (scRNA-seq) data. Linnorm is developed to remove technical noises and simultaneously preserve biological variations in scRNA-seq data, such that existing statistical methods can be improved. Using real scRNA-seq data, we compared Linnorm with existing normalization methods, including NODES, SAMstrt, SCnorm, scran, DESeq and TMM. Linnorm shows advantages in speed, technical noise removal and preservation of cell heterogeneity, which can improve existing methods in the discovery of novel subtypes, pseudo-temporal ordering of cells, clustering analysis, etc. Linnorm also performs better than existing DEG analysis methods, including BASiCS, NODES, SAMstrt, Seurat and DESeq2, in false positive rate control and accuracy.

Two distinct monocyte (Mo)/macrophage (Mp) subsets (Ly6C^low and Ly6C^hi) orchestrate cardiac recovery process following myocardial infarction (MI). Prostaglandin (PG) E₂ is involved in the Mo/Mp-mediated inflammatory response, however, the role of its receptors in Mos/Mps in cardiac healing remains to be determined. Here we show that pharmacological inhibition or gene ablation of the Ep3 receptor in mice suppresses accumulation of Ly6C^low Mos/Mps in infarcted hearts. Ep3 deletion in Mos/Mps markedly attenuates healing after MI by reducing neovascularization in peri-infarct zones. Ep3 deficiency diminishes CX3C chemokine receptor 1 (CX3CR1) expression and vascular endothelial growth factor (VEGF) secretion in Mos/Mps by suppressing TGFβ1 signaling and subsequently inhibits Ly6C^low Mos/Mps migration and angiogenesis. Targeted overexpression of Ep3 receptors in Mos/Mps improves wound healing by enhancing angiogenesis. Thus, the PGE₂/Ep3 axis promotes cardiac healing after MI by activating reparative Ly6C^low Mos/Mps, indicating that Ep3 receptor activation may be a promising therapeutic target for acute MI.

Modeling of transcriptional regulatory networks (TRNs) has been increasingly used to dissect the nature of gene regulation. Inference of regulatory relationships among transcription factors (TFs) and genes, especially among multiple TFs, is still challenging. In this study, we introduced an integrative method, LogicTRN, to decode TF–TF interactions that form TF logics in regulating target genes. By combining cis-regulatory logics and transcriptional kinetics into one single model framework, LogicTRN can naturally integrate dynamic gene expression data and TF-DNA-binding signals in order to identify the TF logics and to reconstruct the underlying TRNs. We evaluated the newly developed methodology using simulation, comparison and application studies, and the results not only show their consistence with existing knowledge, but also demonstrate its ability to accurately reconstruct TRNs in biological complex systems.

Strategies are needed to improve repopulation of decellularized lung scaffolds with stromal and functional epithelial cells. We demonstrate that decellularized mouse lungs recellularized in a dynamic low fluid shear suspension bioreactor, termed the rotating wall vessel (RWV), contained more cells with decreased apoptosis, increased proliferation and enhanced levels of total RNA compared to static recellularization conditions. These results were observed with two relevant mouse cell types: bone marrow-derived mesenchymal stromal (stem) cells (MSCs) and alveolar type II cells (C10). In addition, MSCs cultured in decellularized lungs under static but not bioreactor conditions formed multilayered aggregates. Gene expression and immunohistochemical analyses suggested differentiation of MSCs into collagen I-producing fibroblast-like cells in the bioreactor, indicating enhanced potential for remodeling of the decellularized scaffold matrix. In conclusion, dynamic suspension culture is promising for enhancing repopulation of decellularized lungs, and could contribute to remodeling the extracellular matrix of the scaffolds with subsequent effects on differentiation and functionality of inoculated cells.

We know very little about how soil-borne pollutants such as selenium (Se) can impact pollinators, even though Se has contaminated soils and plants in areas where insect pollination can be critical to the functioning of both agricultural and natural ecosystems. Se can be biotransferred throughout the food web, but few studies have examined its effects on the insects that feed on Se-accumulating plants, particularly pollinators. In laboratory bioassays, we used proboscis extension reflex (PER) and taste perception to determine if the presence of Se affected the gustatory response of honey bee (Apis mellifera L., Hymenoptera: Apidae) foragers. Antennae and proboscises were stimulated with both organic (selenomethionine) and inorganic (selenate) forms of Se that commonly occur in Se-accumulating plants. Methionine was also tested. Each compound was dissolved in 1 M sucrose at 5 concentrations, with sucrose alone as a control. Antennal stimulation with selenomethionine and methionine reduced PER at higher concentrations. Selenate did not reduce gustatory behaviors. Two hours after being fed the treatments, bees were tested for sucrose response threshold. Bees fed selenate responded less to sucrose stimulation. Mortality was higher in bees chronically dosed with selenate compared with a single dose. Selenomethionine did not increase mortality except at the highest concentration. Methionine did not significantly impact survival. Our study has shown that bees fed selenate were less responsive to sucrose, which may lead to a reduction in incoming floral resources needed to support coworkers and larvae in the field. If honey bees forage on nectar containing Se (particularly selenate), reductions in population numbers may occur due to direct toxicity. Given that honey bees are willing to consume food resources containing Se and may not avoid Se compounds in the plant tissues on which they are foraging, they may suffer similar adverse effects as seen in other insect guilds.

Octopamine plays an important role in many behaviors in invertebrates. It acts via binding to G protein coupled receptors located on the plasma membrane of responsive cells. Several distinct subtypes of octopamine receptors have been found in invertebrates, yet little is known about the expression pattern of these different receptor subtypes and how each subtype may contribute to different behaviors. One honey bee (Apis mellifera) octopamine receptor, AmOA1, was recently cloned and characterized. Here we continue to characterize the AmOA1 receptor by investigating its distribution in the honey bee brain. We used two independent antibodies produced against two distinct peptides in the carboxyl-terminus to study the distribution of the AmOA1 receptor in the honey bee brain. We found that both anti-AmOA1 antibodies revealed labeling of cell body clusters throughout the brain and within the following brain neuropils: the antennal lobes; the calyces, pedunculus, vertical (alpha, gamma) and medial (beta) lobes of the mushroom body; the optic lobes; the subesophageal ganglion; and the central complex. Double immunofluorescence staining using anti-GABA and anti-AmOA1 receptor antibodies revealed that a population of inhibitory GABAergic local interneurons in the antennal lobes express the AmOA1 receptor in the cell bodies, axons and their endings in the glomeruli. In the mushroom bodies, AmOA1 receptors are expressed in a subpopulation of inhibitory GABAergic feedback neurons that ends in the visual (outer half of basal ring and collar regions) and olfactory (lip and inner basal ring region) calyx neuropils, as well as in the collar and lip zones of the vertical and medial lobes. The data suggest that one effect of octopamine via AmOA1 in the antennal lobe and mushroom body is to modulate inhibitory neurons.

Neural representations of odors are subject to computations that involve sequentially convergent and divergent anatomical connections across different areas of the brains in both mammals and insects. Furthermore, in both mammals and insects higher order brain areas are connected via feedback connections. In order to understand the transformations and interactions that this connectivity make possible, an ideal experiment would compare neural responses across different, sequential processing levels. Here we present results of recordings from a first order olfactory neuropile – the antennal lobe (AL) – and a higher order multimodal integration and learning center – the mushroom body (MB) – in the honey bee brain. We recorded projection neurons (PN) of the AL and extrinsic neurons (EN) of the MB, which provide the outputs from the two neuropils. Recordings at each level were made in different animals in some experiments and simultaneously in the same animal in others. We presented two odors and their mixture to compare odor response dynamics as well as classification speed and accuracy at each neural processing level. Surprisingly, the EN ensemble significantly starts separating odor stimuli rapidly and before the PN ensemble has reached significant separation. Furthermore the EN ensemble at the MB output reaches a maximum separation of odors between 84–120 ms after odor onset, which is 26 to 133 ms faster than the maximum separation at the AL output ensemble two synapses earlier in processing. It is likely that a subset of very fast PNs, which respond before the ENs, may initiate the rapid EN ensemble response. We suggest therefore that the timing of the EN ensemble activity would allow retroactive integration of its signal into the ongoing computation of the AL via centrifugal feedback.

Investigation into the causes underlying the rapid, global amphibian decline provides critical insight into the effects of changing ecosystems. Hypothesized and confirmed links between amphibian declines, disease, and environmental changes are increasingly represented in published literature. However, there are few long-term amphibian studies that include data on population size, abnormality/injury rates, disease, and habitat variables to adequately assess changes through time. We cultured and identified microorganisms isolated from abnormal/injured and repressed tissue regeneration sites of the endangered Ozark Hellbender, Cryptobranchus alleganiensis bishopi, to discover potential causative agents responsible for their significant decline in health and population. This organism and our study site were chosen because the population and habitat of C. a. bishopi have been intensively studied from 1969–2009, and the abnormality/injury rate and apparent lack of regeneration were established.

Although many bacterial and fungal isolates recovered were common environmental organisms, several opportunistic pathogens were identified in association with only the injured tissues of C.a. bishopi. Bacterial isolates included Aeromonas hydrophila, a known amphibian pathogen, Granulicetella adiacens, Gordonai terrae, Stenotrophomonas maltophilia, Aerococcus viridans, Streptococcus pneumoniae and a variety of Pseudomonads, including Pseudomonas aeruginosa, P. stutzeri, and P. alcaligenes. Fungal isolates included species in the genera Penicillium, Acremonium, Cladosporium, Curvularia, Fusarium, Streptomycetes, and the Class Hyphomycetes. Many of the opportunistic pathogens identified are known to form biofilms. Lack of isolation of the same organism from all wounds suggests that the etiological agent responsible for the damage to C. a. bishopi may not be a single organism. To our knowledge, this is the first study to profile the external microbial consortia cultured from a Cryptobranchid salamander. The incidence of abnormalities/injury and retarded regeneration in C. a. bishopi may have many contributing factors including disease and habitat degradation. Results from this study may provide insight into other amphibian population declines.

It remains challenging to predict regulatory variants in particular tissues or cell types due to highly context-specific gene regulation. By connecting large-scale epigenomic profiles to expression quantitative trait loci (eQTLs) in a wide range of human tissues/cell types, we identify critical chromatin features that predict variant regulatory potential. We present cepip, a joint likelihood framework, for estimating a variant’s regulatory probability in a context-dependent manner. Our method exhibits significant GWAS signal enrichment and is superior to existing cell type-specific methods. Furthermore, using phenotypically relevant epigenomes to weight the GWAS single-nucleotide polymorphisms, we improve the statistical power of the gene-based association test.

A distinct pathovar of Salmonella enterica serovar Typhimurium, ST313, has emerged in sub-Saharan Africa as a major cause of fatal bacteremia in young children and HIV-infected adults. D23580, a multidrug resistant clinical isolate of ST313, was previously shown to have undergone genome reduction in a manner that resembles that of the more human-restricted pathogen, Salmonella enterica serovar Typhi. It has since been shown through tissue distribution studies that D23580 is able to establish an invasive infection in chickens. However, it remains unclear whether ST313 can cause lethal disease in a non-human host following a natural course of infection. Herein we report that D23580 causes lethal and invasive disease in a murine model of infection following peroral challenge. The LD⁵⁰ of D23580 in female BALB/c mice was 4.7 x 10⁵ CFU. Tissue distribution studies performed 3 and 5 days post-infection confirmed that D23580 was able to more rapidly colonize the spleen, mesenteric lymph nodes and gall bladder in mice when compared to the well-characterized S. Typhimurium strain SL1344. D23580 exhibited enhanced resistance to acid stress relative to SL1344, which may lend towards increased capability to survive passage through the gastrointestinal tract as well as during its intracellular lifecycle. Interestingly, D23580 also displayed higher swimming motility relative to SL1344, S. Typhi strain Ty2, and the ST313 strain A130. Biochemical tests revealed that D23580 shares many similar metabolic features with SL1344, with several notable differences in the Voges-Proskauer and catalase tests, as well alterations in melibiose, and inositol utilization. These results represent the first full duration infection study using an ST313 strain following the entire natural course of disease progression, and serve as a benchmark for ongoing and future studies into the pathogenesis of D23580.