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.

Contemporary human populations conform to ecogeographic predictions that animals will become more compact in cooler climates and less compact in warmer ones. However, it remains unclear to what extent this pattern reflects plastic responses to current environments or genetic differences among populations. Analyzing anthropometric surveys of 232,684 children and adults from across 80 ethnolinguistic groups in sub-Saharan Africa, Asia and the Americas, we confirm that body surface-to-volume correlates with contemporary temperature at magnitudes found in more latitudinally diverse samples (Adj. R² = 0.14-0.28). However, far more variation in body surface-to-volume is attributable to genetic population structure (Adj. R² = 0.50-0.74). Moreover, genetic population structure accounts for nearly all of the observed relationship between contemporary temperature and body surface-to-volume among children and adults. Indeed, after controlling for population structure, contemporary temperature accounts for no more than 4% of the variance in body form in these groups. This effect of genetic affinity on body form is also independent of other ecological variables, such as dominant mode of subsistence and household wealth per capita. These findings suggest that the observed fit of human body surface-to-volume with current climate in this sample reflects relatively large effects of existing genetic population structure of contemporary humans compared to plastic response to current environments.

Background: The transition from the home to college is a phase in which emerging adults shift toward more unhealthy eating and physical activity patterns, higher body mass indices, thus increasing risk of overweight/obesity. Currently, little is understood about how changing friendship networks shape weight gain behaviors. This paper describes the recruitment, data collection, and data analytic protocols for the SPARC (Social impact of Physical Activity and nutRition in College) study, a longitudinal examination of the mechanisms by which friends and friendship networks influence nutrition and physical activity behaviors and weight gain in the transition to college life.

Methods: The SPARC study aims to follow 1450 university freshmen from a large university over an academic year, collecting data on multiple aspects of friends and friendship networks. Integrating multiple types of data related to student lives, ecological momentary assessments (EMAs) are administered via a cell phone application, devilSPARC. EMAs collected in four 1-week periods (a total of 4 EMA waves) are integrated with linked data from web-based surveys and anthropometric measurements conducted at four times points (for a total of eight data collection periods including EMAs, separated by ~1 month). University databases will provide student card data, allowing integration of both time-dated data on food purchasing, use of physical activity venues, and geographical information system (GIS) locations of these activities relative to other students in their social networks.

Discussion: Findings are intended to guide the development of more effective interventions to enhance behaviors among college students that protect against weight gain during college.

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

The heterocyclic indole-alkaloid scytonemin is a sunscreen found exclusively among cyanobacteria. An 18-gene cluster is responsible for scytonemin production in Nostoc punctiforme ATCC 29133. The upstream genes scyABCDEF in the cluster are proposed to be responsible for scytonemin biosynthesis from aromatic amino acid substrates. In vitro studies of ScyA, ScyB, and ScyC proved that these enzymes indeed catalyze initial pathway reactions. Here we characterize the role of ScyD, ScyE, and ScyF, which were logically predicted to be responsible for late biosynthetic steps, in the biological context of N. punctiforme. In-frame deletion mutants of each were constructed (ΔscyD, ΔscyE, and ΔscyF) and their phenotypes studied. Expectedly, ΔscyE presents a scytoneminless phenotype, but no accumulation of the predicted intermediaries. Surprisingly, ΔscyD retains scytonemin production, implying that it is not required for biosynthesis. Indeed, scyD presents an interesting evolutionary paradox: it likely originated in a duplication event from scyE, and unlike other genes in the operon, it has not been subjected to purifying selection. This would suggest that it is a pseudogene, and yet scyD is highly conserved in the scytonemin operon of cyanobacteria. ΔscyF also retains scytonemin production, albeit exhibiting a reduction of the production yield compared with the wild-type. This indicates that ScyF is not essential but may play an adjuvant role for scytonemin synthesis. Altogether, our findings suggest that these downstream genes are not responsible, as expected, for the late steps of scytonemin synthesis and we must look for those functions elsewhere. These findings are particularly important for biotechnological production of this sunscreen through heterologous expression of its genes in more tractable organisms.

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: Prior studies have shown that using uterotonics to augment or induce labor before arrival at comprehensive Emergency Obstetric and Neonatal Care (CEmONC) settings (henceforth, “outside uterotonics”) may contribute to perinatal mortality in low- and middle-income countries. We estimate its effect on perinatal mortality in rural Bangladesh.

Methods: Using hospital records (23986 singleton term births, Jan 1, 2009-Dec 31, 2015) from rural Bangladesh, we use a logistic regression model to estimate the increased risk of perinatal death from uterotonics administered outside a CEmONC facility.

Results: Among term births (≥37 weeks gestation), the risk of perinatal death adjusted for key confounders is significantly increased among women reporting uterotonic use outside of CEmONC (OR = 3 · 0, 95 % CI = 2 · 4,3 · 7). This increased risk is particularly high for fresh stillbirths (OR = 4 · 0, 95 % CI = 3 · 0,5 · 3) and intrapartum-related causes of early neonatal deaths (birth asphyxia) (OR = 3 · 1, 95 % CI = 2 · 2,4 · 5).

Conclusions: In this sample, outside uterotonic use was associated with substantially increased risk of fresh stillbirths, deaths due to birth asphyxia, and all perinatal deaths. In settings of high uterotonic use outside of controlled settings, substantial improvement in both stillbirth and early neonatal mortality may be made by reducing such use.

N₂ fixation and ammonia oxidation (AO) are the two most important processes in the nitrogen (N) cycle of biological soil crusts (BSCs). We studied the short-term response of acetylene reduction assay (ARA) rates, an indicator of potential N₂ fixation, and AO rates to temperature (T, -5°C to 35°C) in BSC of different successional stages along the BSC ecological succession and geographic origin (hot Chihuahuan and cooler Great Basin deserts). ARA in all BSCs increased with T until saturation occurred between 15 and 20°C, and declined at 30–35°C. Culture studies using cyanobacteria isolated from these crusts indicated that the saturating effect was traceable to their inability to grow well diazotrophically within the high temperature range. Below saturation, temperature response was exponential, with Q₁₀ significantly different in the two areas (~ 5 for Great Basin BSCs; 2–3 for Chihuahuan BSCs), but similar between the two successional stages. However, in contrast to ARA, AO showed a steady increase to 30–35°C in Great Basin, and Chihuhuan BSCs showed no inhibition at any tested temperature. The T response of AO also differed significantly between Great Basin (Q₁₀ of 4.5–4.8) and Chihuahuan (Q₁₀ of 2.4–2.6) BSCs, but not between successional stages. Response of ARA rates to T did not differ from that of AO in either desert. Thus, while both processes scaled to T in unison until 20°C, they separated to an increasing degree at higher temperature. As future warming is likely to occur in the regions where BSCs are often the dominant living cover, this predicted decoupling is expected to result in higher proportion of nitrates in soil relative to ammonium. As nitrate is more easily lost as leachate or to be reduced to gaseous forms, this could mean a depletion of soil N over large landscapes globally.