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.

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

We present a phylogeographic study of at least six reproductively isolated lineages of new world harvester ants within the Pogonomyrmex barbatus and P. rugosus species group. The genetic and geographic relationships within this clade are complex: Four of the identified lineages show genetic caste determination (GCD) and are divided into two pairs. Each pair has evolved under a mutualistic system that necessitates sympatry. These paired lineages are dependent upon one another because their GCD requires interlineage matings for the production of F1 hybrid workers, and intralineage matings are required to produce queens. This GCD system maintains genetic isolation among these interdependent lineages, while simultaneously requiring co-expansion and emigration as their distributions have changed over time. It has also been demonstrated that three of these four GCD lineages have undergone historical hybridization, but the narrower sampling range of previous studies has left questions on the hybrid parentage, breadth, and age of these groups. Thus, reconstructing the phylogenetic and geographic history of this group allows us to evaluate past insights and hypotheses and to plan future inquiries in a more complete historical biogeographic context. Using mitochondrial DNA sequences sampled across most of the morphospecies’ ranges in the U.S.A. and Mexico, we conducted a detailed phylogeographic study. Remarkably, our results indicate that one of the GCD lineage pairs has experienced a dramatic range expansion, despite the genetic load and fitness costs of the GCD system. Our analyses also reveal a complex pattern of vicariance and dispersal in Pogonomyrmex harvester ants that is largely concordant with models of late Miocene, Pliocene, and Pleistocene range shifts among various arid-adapted taxa in North America.

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.

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.

The electronic band structure of MoS₂, MoSe₂, WS₂, and WSe₂, crystals has been studied at various hydrostatic pressures experimentally by photoreflectance (PR) spectroscopy and theoretically within the density functional theory (DFT). In the PR spectra direct optical transitions (A and B) have been clearly observed and pressure coefficients have been determined for these transitions to be: α_A = 2.0 ± 0.1 and α_B = 3.6 ± 0.1 meV/kbar for MoS₂, α_A = 2.3 ± 0.1 and α_B = 4.0 ± 0.1 meV/kbar for MoSe₂, α_A = 2.6 ± 0.1 and α_B = 4.1 ± 0.1 meV/kbar for WS₂, α_A = 3.4 ± 0.1 and α_B = 5.0 ± 0.5 meV/kbar for WSe₂. It has been found that these coefficients are in an excellent agreement with theoretical predictions. In addition, a comparative study of different computational DFT approaches has been performed and analyzed. For indirect gap the pressure coefficient have been determined theoretically to be −7.9, −5.51, −6.11, and −3.79, meV/kbar for MoS₂, MoSe₂, WS₂, and WSe₂, respectively. The negative values of this coefficients imply a narrowing of the fundamental band gap with the increase in hydrostatic pressure and a semiconductor to metal transition for MoS₂, MoSe₂, WS₂, and WSe₂, crystals at around 140, 180, 190, and 240 kbar, respectively.

Binary transition metal dichalcogenide monolayers share common properties such as a direct optical bandgap, spin-orbit splittings of hundreds of meV, light–matter interaction dominated by robust excitons and coupled spin-valley states. Here we demonstrate spin-orbit-engineering in Mo[_(1-x)]W_xSe₂ alloy monolayers for optoelectronics and applications based on spin- and valley-control. We probe the impact of the tuning of the conduction band spin-orbit spin-splitting on the bright versus dark exciton population. For MoSe₂ monolayers, the photoluminescence intensity decreases as a function of temperature by an order of magnitude (4–300 K), whereas for WSe₂ we measure surprisingly an order of magnitude increase. The ternary material shows a trend between these two extreme behaviors. We also show a non-linear increase of the valley polarization as a function of tungsten concentration, where 40% tungsten incorporation is sufficient to achieve valley polarization as high as in binary WSe₂.

We present two-dimensional Mg(OH)₂ sheets and their vertical heterojunctions with CVD-MoS₂ for the first time as flexible 2D insulators with anomalous lattice vibration and chemical and physical properties. New hydrothermal crystal growth technique enabled isolation of environmentally stable monolayer Mg(OH)₂ sheets. Raman spectroscopy and vibrational calculations reveal that the lattice vibrations of Mg(OH)₂ have fundamentally different signature peaks and dimensionality effects compared to other 2D material systems known to date. Sub-wavelength electron energy-loss spectroscopy measurements and theoretical calculations show that Mg(OH)₂ is a 6 eV direct-gap insulator in 2D, and its optical band gap displays strong band renormalization effects from monolayer to bulk, marking the first experimental confirmation of confinement effects in 2D insulators. Interestingly, 2D-Mg(OH)₂ sheets possess rather strong surface polarization (charge) effects which is in contrast to electrically neutral h-BN materials. Using 2D-Mg(OH)₂ sheets together with CVD-MoS₂ in the vertical stacking shows that a strong change transfer occurs from n-doped CVD-MoS₂ sheets to Mg(OH)₂, naturally depleting the semiconductor, pushing towards intrinsic doping limit and enhancing overall optical performance of 2D semiconductors. Results not only establish unusual confinement effects in 2D-Mg(OH)₂, but also offer novel 2D-insulating material with unique physical, vibrational, and chemical properties for potential applications in flexible optoelectronics.

Transition metal trichalcogenides form a class of layered materials with strong in-plane anisotropy. For example, titanium trisulfide (TiS₃) whiskers are made out of weakly interacting TiS₃ layers, where each layer is made of weakly interacting quasi-one-dimensional chains extending along the b axis. Here we establish the unusual vibrational properties of TiS₃ both experimentally and theoretically. Unlike other two-dimensional systems, the Raman active peaks of TiS₃ have only out-of-plane vibrational modes, and interestingly some of these vibrations involve unique rigid-chain vibrations and S–S molecular oscillations. High-pressure Raman studies further reveal that the A_g^S-S S-S molecular mode has an unconventional negative pressure dependence, whereas other peaks stiffen as anticipated. Various vibrational modes are doubly degenerate at ambient pressure, but the degeneracy is lifted at high pressures. These results establish the unusual vibrational properties of TiS₃ with strong in-plane anisotropy, and may have relevance to understanding of vibrational properties in other anisotropic two-dimensional material systems.