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.

High-resolution, global quantification of fossil fuel CO[subscript 2] emissions is emerging as a critical need in carbon cycle science and climate policy. We build upon a previously developed fossil fuel data assimilation system (FFDAS) for estimating global high-resolution fossil fuel CO[subscript 2] emissions. We have improved the underlying observationally based data sources, expanded the approach through treatment of separate emitting sectors including a new pointwise database of global power plants, and extended the results to cover a 1997 to 2010 time series at a spatial resolution of 0.1°. Long-term trend analysis of the resulting global emissions shows subnational spatial structure in large active economies such as the United States, China, and India. These three countries, in particular, show different long-term trends and exploration of the trends in nighttime lights, and population reveal a decoupling of population and emissions at the subnational level. Analysis of shorter-term variations reveals the impact of the 2008–2009 global financial crisis with widespread negative emission anomalies across the U.S. and Europe. We have used a center of mass (CM) calculation as a compact metric to express the time evolution of spatial patterns in fossil fuel CO[subscript 2] emissions. The global emission CM has moved toward the east and somewhat south between 1997 and 2010, driven by the increase in emissions in China and South Asia over this time period. Analysis at the level of individual countries reveals per capita CO[subscript 2] emission migration in both Russia and India. The per capita emission CM holds potential as a way to succinctly analyze subnational shifts in carbon intensity over time. Uncertainties are generally lower than the previous version of FFDAS due mainly to an improved nightlight data set.

Large urban emissions of greenhouse gases result in large atmospheric enhancements relative to background that are easily measured. Using CO₂ mole fractions and Δ¹⁴C and δ¹³C values of CO₂ in the Los Angeles megacity observed in inland Pasadena (2006–2013) and coastal Palos Verdes peninsula (autumn 2009–2013), we have determined time series for CO₂ contributions from fossil fuel combustion (C_ff) for both sites and broken those down into contributions from petroleum and/or gasoline and natural gas burning for Pasadena. We find a 10 % reduction in Pasadena C_ff during the Great Recession of 2008–2010, which is consistent with the bottom-up inventory determined by the California Air Resources Board. The isotopic variations and total atmospheric CO₂ from our observations are used to infer seasonality of natural gas and petroleum combustion. The trend of CO₂ contributions to the atmosphere from natural gas combustion is out of phase with the seasonal cycle of total natural gas combustion seasonal patterns in bottom-up inventories but is consistent with the seasonality of natural gas usage by the area's electricity generating power plants. For petroleum, the inferred seasonality of CO₂ contributions from burning petroleum is delayed by several months relative to usage indicated by statewide gasoline taxes. Using the high-resolution Hestia-LA data product to compare C_ff from parts of the basin sampled by winds at different times of year, we find that variations in observed fossil fuel CO₂ reflect seasonal variations in wind direction. The seasonality of the local CO₂ excess from fossil fuel combustion along the coast, on Palos Verdes peninsula, is higher in autumn and winter than spring and summer, almost completely out of phase with that from Pasadena, also because of the annual variations of winds in the region. Variations in fossil fuel CO₂ signals are consistent with sampling the bottom-up Hestia-LA fossil CO₂ emissions product for sub-city source regions in the LA megacity domain when wind directions are considered.

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.

Atmospheric CO₂ inversions estimate surface carbon fluxes from an optimal fit to atmospheric CO₂ measurements, usually including prior constraints on the flux estimates. Eleven sets of carbon flux estimates are compared, generated by different inversions systems that vary in their inversions methods, choice of atmospheric data, transport model and prior information. The inversions were run for at least 5 yr in the period between 1990 and 2010. Mean fluxes for 2001–2004, seasonal cycles, interannual variability and trends are compared for the tropics and northern and southern extra-tropics, and separately for land and ocean. Some continental/basin-scale subdivisions are also considered where the atmospheric network is denser. Four-year mean fluxes are reasonably consistent across inversions at global/latitudinal scale, with a large total (land plus ocean) carbon uptake in the north (−3.4 Pg C yr^-1 (±0.5 Pg C yr^-1 standard deviation), with slightly more uptake over land than over ocean), a significant although more variable source over the tropics (1.6 ± 0.9 Pg C yr^-1) and a compensatory sink of similar magnitude in the south (−1.4 ± 0.5 Pg C yr^-1) corresponding mainly to an ocean sink. Largest differences across inversions occur in the balance between tropical land sources and southern land sinks. Interannual variability (IAV) in carbon fluxes is larger for land than ocean regions (standard deviation around 1.06 versus 0.33 Pg C yr[superscript −1] for the 1996–2007 period), with much higher consistency among the inversions for the land. While the tropical land explains most of the IAV (standard deviation ~ 0.65 Pg C yr^-1), the northern and southern land also contribute (standard deviation ~ 0.39 Pg C yr^-1). Most inversions tend to indicate an increase of the northern land carbon uptake from late 1990s to 2008 (around 0.1 Pg C yr^-1, predominantly in North Asia. The mean seasonal cycle appears to be well constrained by the atmospheric data over the northern land (at the continental scale), but still highly dependent on the prior flux seasonality over the ocean. Finally we provide recommendations to interpret the regional fluxes, along with the uncertainty estimates.

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.

This paper presents an analysis of methane emissions from the Los Angeles Basin at monthly timescales across a 4-year time period – from September 2011 to August 2015. Using observations acquired by a ground-based near-infrared remote sensing instrument on Mount Wilson, California, combined with atmospheric CH₄–CO₂ tracer–tracer correlations, we observed −18 to +22 % monthly variability in CH₄ : CO₂ from the annual mean in the Los Angeles Basin. Top-down estimates of methane emissions for the basin also exhibit significant monthly variability (−19 to +31 % from annual mean and a maximum month-to-month change of 47 %). During this period, methane emissions consistently peaked in the late summer/early fall and winter. The estimated annual methane emissions did not show a statistically significant trend over the 2011 to 2015 time period.

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.

Attributing observed CO2 variations to human or natural cause is critical to deducing and tracking emissions from observations. We have used in situ CO2, CO, and planetary boundary layer height (PBLH) measurements recorded during the CalNex-LA (CARB et al., 2008) ground campaign of 15 May-15 June 2010, in Pasadena, CA, to deduce the diurnally varying anthropogenic component of observed CO2 in the megacity of Los Angeles (LA). This affordable and simple technique, validated by carbon isotope observations and WRF-STILT (Weather Research and Forecasting model - Stochastic Time-Inverted Lagrangian Transport model) predictions, is shown to robustly attribute observed CO2 variation to anthropogenic or biogenic origin over the entire diurnal cycle. During CalNex-LA, local fossil fuel combustion contributed up to similar to 50% of the observed CO2 enhancement overnight, and similar to 100% of the enhancement near midday. This suggests that sufficiently accurate total column CO2 observations recorded near midday, such as those from the GOSAT or OCO-2 satellites, can potentially be used to track anthropogenic emissions from the LA megacity.

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.