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.

Background: While there is ample evidence for health risks associated with heat and other extreme weather events today, little is known about the impact of weather patterns on population health in preindustrial societies.

Objective: To investigate the impact of weather patterns on population health in Sweden before and during industrialization.

Methods: We obtained records of monthly mortality and of monthly mean temperatures and precipitation for Skellefteå parish, northern Sweden, for the period 1800-1950. The associations between monthly total mortality, as well as monthly mortality due to infectious and cardiovascular diseases, and monthly mean temperature and cumulative precipitation were modelled using a time series approach for three separate periods, 1800−1859, 1860-1909, and 1910-1950.

Results: We found higher temperatures and higher amounts of precipitation to be associated with lower mortality both in the medium term (same month and two-months lag) and in the long run (lag of six months up to a year). Similar patterns were found for mortality due to infectious and cardiovascular diseases. Furthermore, the effect of temperature and precipitation decreased over time.

Conclusions: Higher temperature and precipitation amounts were associated with reduced death counts with a lag of up to 12 months. The decreased effect over time may be due to improvements in nutritional status, decreased infant deaths, and other changes in society that occurred in the course of the demographic and epidemiological transition.

Contribution: The study contributes to a better understanding of the complex relationship between weather and mortality and, in particular, historical weather-related mortality.

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.

Background: Most excess deaths that occur during extreme hot weather events do not have natural heat recorded as an underlying or contributing cause. This study aims to identify the specific individuals who died because of hot weather using only secondary data. A novel approach was developed in which the expected number of deaths was repeatedly sampled from all deaths that occurred during a hot weather event, and compared with deaths during a control period. The deaths were compared with respect to five factors known to be associated with hot weather mortality. Individuals were ranked by their presence in significant models over 100 trials of 10,000 repetitions. Those with the highest rankings were identified as probable excess deaths. Sensitivity analyses were performed on a range of model combinations. These methods were applied to a 2009 hot weather event in greater Vancouver, Canada.

Results: The excess deaths identified were sensitive to differences in model combinations, particularly between univariate and multivariate approaches. One multivariate and one univariate combination were chosen as the best models for further analyses. The individuals identified by multiple combinations suggest that marginalized populations in greater Vancouver are at higher risk of death during hot weather.

Conclusions: This study proposes novel methods for classifying specific deaths as expected or excess during a hot weather event. Further work is needed to evaluate performance of the methods in simulation studies and against clinically identified cases. If confirmed, these methods could be applied to a wide range of populations and events of interest.

The Quadrangles Av-11 and Av-12 on Vesta are located at the northern rim of the giant Rheasilvia south polar impact basin. The primary geologic units in Av-11 and Av-12 include material from the Rheasilvia impact basin formation, smooth material and different types of impact crater structures (such as bimodal craters, dark and bright crater ray material and dark ejecta material). Av-11 and Av-12 exhibit almost the full range of mass wasting features observed on Vesta, such as slump blocks, spur-and-gully morphologies and landslides within craters. Processes of collapse, slope instability and seismically triggered events force material to slump down crater walls or scarps and produce landslides or rotational slump blocks. The spur-and-gully morphology that is known to form on Mars is also observed on Vesta; however, on Vesta this morphology formed under dry conditions.

A variety of geologic landforms and features are observed within quadrangle Av-13 Tuccia in the southern hemisphere of Vesta. The quadrangle covers parts of the highland Vestalia Terra as well as the floors of the large Rheasilvia and Veneneia impact basins, which results in a substantial elevation difference of more than 40 km between the northern and the southern portions of the quadrangle. Measurements of crater size–frequency distributions within and surrounding the Rheasilvia basin indicate that gravity-driven mass wasting in the interior of the basin has been important, and that the basin has a more ancient formation age than would be expected from the crater density on the basin floor alone. Subsequent to its formation, Rheasilvia was superimposed by several mid-sized impact craters. The most prominent craters are Tuccia, Eusebia, Vibidia, Galeria, and Antonia, whose geology and formation ages are investigated in detail in this work. These impact structures provide a variety of morphologies indicating different sorts of subsequent impact-related or gravity-driven mass wasting processes. Understanding the geologic history of the relatively young craters in the Rheasilvia basin is important in order to understand the even more degraded craters in other regions of Vesta.

In this paper we present a time-stratigraphic scheme and geologic time scale for the protoplanet Vesta, based on global geologic mapping and other analyses of NASA Dawn spacecraft data, complemented by insights gained from laboratory studies of howardite–eucrite–diogenite (HED) meteorites and geophysical modeling. On the basis of prominent impact structures and their associated deposits, we propose a time scale for Vesta that consists of four geologic time periods: Pre-Veneneian, Veneneian, Rheasilvian, and Marcian. The Pre-Veneneian Period covers the time from the formation of Vesta up to the Veneneia impact event, from 4.6 Ga to >2.1 Ga (using the asteroid flux-derived chronology system) or from 4.6 Ga to 3.7 Ga (under the lunar-derived chronology system). The Veneneian Period covers the time span between the Veneneia and Rheasilvia impact events, from >2.1 to 1 Ga (asteroid flux-derived chronology) or from 3.7 to 3.5 Ga (lunar-derived chronology), respectively. The Rheasilvian Period covers the time span between the Rheasilvia and Marcia impact events, and the Marcian Period covers the time between the Marcia impact event until the present. The age of the Marcia impact is still uncertain, but our current best estimates from crater counts of the ejecta blanket suggest an age between ∼120 and 390 Ma, depending upon choice of chronology system used. Regardless, the Marcia impact represents the youngest major geologic event on Vesta. Our proposed four-period geologic time scale for Vesta is, to a first order, comparable to those developed for other airless terrestrial bodies.

We report on a preliminary global geologic map of Vesta, based on data from the Dawn spacecraft’s High-Altitude Mapping Orbit (HAMO) and informed by Low-Altitude Mapping Orbit (LAMO) data. This map is part of an iterative mapping effort; the geologic map has been refined with each improvement in resolution. Vesta has a heavily-cratered surface, with large craters evident in numerous locations. The south pole is dominated by an impact structure identified before Dawn’s arrival. Two large impact structures have been resolved: the younger, larger Rheasilvia structure, and the older, more degraded Veneneia structure. The surface is also characterized by a system of deep, globe-girdling equatorial troughs and ridges, as well as an older system of troughs and ridges to the north. Troughs and ridges are also evident cutting across, and spiraling arcuately from, the Rheasilvia central mound.

However, no volcanic features have been unequivocally identified. Vesta can be divided very broadly into three terrains: heavily-cratered terrain; ridge-and-trough terrain (equatorial and northern); and terrain associated with the Rheasilvia crater. Localized features include bright and dark material and ejecta (some defined specifically by color); lobate deposits; and mass-wasting materials. No obvious volcanic features are evident. Stratigraphy of Vesta’s geologic units suggests a history in which formation of a primary crust was followed by the formation of impact craters, including Veneneia and the associated Saturnalia Fossae unit. Formation of Rheasilvia followed, along with associated structural deformation that shaped the Divalia Fossae ridge-and-trough unit at the equator. Subsequent impacts and mass wasting events subdued impact craters, rims and portions of ridge-and-trough sets, and formed slumps and landslides, especially within crater floors and along crater rims and scarps. Subsequent to the formation of Rheasilvia, discontinuous low-albedo deposits formed or were emplaced; these lie stratigraphically above the equatorial ridges that likely were formed by Rheasilvia. The last features to be formed were craters with bright rays and other surface mantling deposits.

Executed progressively throughout data acquisition, the iterative mapping process provided the team with geologic proto-units in a timely manner. However, interpretation of the resulting map was hampered by the necessity to provide the team with a standard nomenclature and symbology early in the process. With regard to mapping and interpreting units, the mapping process was hindered by the lack of calibrated mineralogic information. Topography and shadow played an important role in discriminating features and terrains, especially in the early stages of data acquisition.