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.

Introduction: Urbanization can considerably impact animal ecology, evolution, and behavior. Among the new conditions that animals experience in cities is anthropogenic noise, which can limit the sound space available for animals to communicate using acoustic signals. Some urban bird species increase their song frequencies so that they can be heard above low-frequency background city noise. However, the ability to make such song modifications may be constrained by several morphological factors, including bill gape, size, and shape, thereby limiting the degree to which certain species can vocally adapt to urban settings. We examined the relationship between song characteristics and bill morphology in a species (the house finch, Haemorhous mexicanus) where both vocal performance and bill size are known to differ between city and rural animals.

Results: We found that bills were longer and narrower in more disturbed, urban areas. We observed an increase in minimum song frequency of urban birds, and we also found that the upper frequency limit of songs decreased in direct relation to bill morphology.

Conclusions: These findings are consistent with the hypothesis that birds with longer beaks and therefore longer vocal tracts sing songs with lower maximum frequencies because longer tubes have lower-frequency resonances. Thus, for the first time, we reveal dual constraints (one biotic, one abiotic) on the song frequency range of urban animals. Urban foraging pressures may additionally interact with the acoustic environment to shape bill traits and vocal performance.

The termites evolved eusociality and complex societies before the ants, but have been studied much less. The recent publication of the first two termite genomes provides a unique comparative opportunity, particularly because the sequenced termites represent opposite ends of the social complexity spectrum. Zootermopsis nevadensis has simple colonies with totipotent workers that can develop into all castes (dispersing reproductives, nest-inheriting replacement reproductives, and soldiers). In contrast, the fungus-growing termite Macrotermes natalensis belongs to the higher termites and has very large and complex societies with morphologically distinct castes that are life-time sterile. Here we compare key characteristics of genomic architecture, focusing on genes involved in communication, immune defenses, mating biology and symbiosis that were likely important in termite social evolution. We discuss these in relation to what is known about these genes in the ants and outline hypothesis for further testing.

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.

Color vision in birds is mediated by four types of cone photoreceptors whose maximal sensitivities (λ_max) are evenly spaced across the light spectrum. In the course of avian evolution, the λ_max of the most shortwave-sensitive cone, SWS1, has switched between violet (λ_max > 400 nm) and ultraviolet (λ_max < 380 nm) multiple times. This shift of the SWS1 opsin is accompanied by a corresponding short-wavelength shift in the spectrally adjacent SWS2 cone. Here, we show that SWS2 cone spectral tuning is mediated by modulating the ratio of two apocarotenoids, galloxanthin and 11’,12’-dihydrogalloxanthin, which act as intracellular spectral filters in this cell type. We propose an enzymatic pathway that mediates the differential production of these apocarotenoids in the avian retina, and we use color vision modeling to demonstrate how correlated evolution of spectral tuning is necessary to achieve even sampling of the light spectrum and thereby maintain near-optimal color discrimination.

Eusocial insects, mostly Hymenoptera, have evolved unique colonial lifestyles that rely on the perception of social context mainly through pheromones, and chemoreceptors are hypothesized to have played important adaptive roles in the evolution of sociality. However, because chemoreceptor repertoires have been characterized in few social insects and their solitary relatives, a comprehensive examination of this hypothesis has not been possible. Here, we annotate ∼3,000 odorant and gustatory receptors in recently sequenced Hymenoptera genomes and systematically compare >4,000 chemoreceptors from 13 hymenopterans, representing one solitary lineage (wasps) and three independently evolved eusocial lineages (ants and two bees). We observe a strong general tendency for chemoreceptors to expand in Hymenoptera, whereas the specifics of gene gains/losses are highly diverse between lineages. We also find more frequent positive selection on chemoreceptors in a facultative eusocial bee and in the common ancestor of ants compared with solitary wasps. Our results suggest that the frequent expansions of chemoreceptors have facilitated the transition to eusociality. Divergent expression patterns of odorant receptors between honeybee and ants further indicate differential roles of chemoreceptors in parallel trajectories of social evolution.

Epigenetic inheritance plays an important role in mediating alternative phenotype in highly social species. In order to gain a greater understanding of epigenetic effects in societies, we investigated DNA methylation in the termite Zootermopsis nevadensis. Termites are the most ancient social insects, and developmentally distinct from highly-studied, hymenopteran social insects. We used replicated bisulfite-sequencing to investigate patterns of DNA methylation in both sexes and among castes of Z. nevadensis. We discovered that Z. nevadensis displayed some of the highest levels of DNA methylation found in insects. We also found strong differences in methylation between castes. Methylated genes tended to be uniformly and highly expressed demonstrating the antiquity of associations between intragenic methylation and gene expression. Differentially methylated genes were more likely to be alternatively spliced than not differentially methylated genes, and possessed considerable enrichment for development-associated functions. We further observed strong overrepresentation of multiple transcription factor binding sites and miRNA profiles associated with differential methylation, providing new insights into the possible function of DNA methylation. Overall, our results show that DNA methylation is widespread and associated with caste differences in termites. More generally, this study provides insights into the function of DNA methylation and the success of insect societies.

Gut-associated microbiota of ants include Rhizobiales bacteria with affiliation to the genus Bartonella. These bacteria may enable the ants to fix atmospheric nitrogen, but no genomes have been sequenced yet to test the hypothesis. Sequence reads from a member of the Rhizobiales were identified in the data collected in a genome project of the ant Harpegnathos saltator. We present an analysis of the closed 1.86 Mb genome of the ant-associated bacterium, for which we suggest the species name Candidatus Tokpelaia hoelldoblerii. A phylogenetic analysis reveals a relationship to Bartonella and Brucella, which infect mammals. Novel gene acquisitions include a gene for a putative extracellular protein of more than 6,000 amino acids secreted by the type I secretion system, which may be involved in attachment to the gut epithelium. No genes for nitrogen fixation could be identified, but genes for a multi-subunit urease protein complex are present in the genome. The urease genes are also present in Brucella, which has a fecal-oral transmission pathway, but not in Bartonella, which use blood-borne transmission pathways. We hypothesize that the gain and loss of the urease function is related to transmission strategies and lifestyle changes in the host-associated members of the Rhizobiales.

Accumulating data from genome-wide association studies (GWAS) have provided a collection of novel candidate genes associated with complex diseases, such as atherosclerosis. We identified an atherosclerosis-associated single-nucleotide polymorphism (SNP) located in the intron of the long noncoding RNA (lncRNA) LINC00305 by searching the GWAS database. Although the function of LINC00305 is unknown, we found that LINC00305 expression is enriched in atherosclerotic plaques and monocytes. Overexpression of LINC00305 promoted the expression of inflammation-associated genes in THP-1 cells and reduced the expression of contractile markers in co-cultured human aortic smooth muscle cells (HASMCs). We showed that overexpression of LINC00305 activated nuclear factor-kappa beta (NF-κB) and that inhibition of NF-κB abolished LINC00305-mediated activation of cytokine expression. Mechanistically, LINC00305 interacted with lipocalin-1 interacting membrane receptor (LIMR), enhanced the interaction of LIMR and aryl-hydrocarbon receptor repressor (AHRR), and promoted protein expression as well as nuclear localization of AHRR. Moreover, LINC00305 activated NF-κB exclusively in the presence of LIMR and AHRR. In light of these findings, we propose that LINC00305 promotes monocyte inflammation by facilitating LIMR and AHRR cooperation and the AHRR activation, which eventually activates NF-κB, thereby inducing HASMC phenotype switching.