Faculty and Staff

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.

The emergence of the disease chytridiomycosis caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd) has been implicated in dramatic global amphibian declines. Although many species have undergone catastrophic declines and/or extinctions, others appear to be unaffected or persist at reduced frequencies after Bd outbreaks. The reasons behind this variance in disease outcomes are poorly understood: differences in host immune responses have been proposed, yet previous studies suggest a lack of robust immune responses to Bd in susceptible species. Here, we sequenced transcriptomes from clutch-mates of a highly susceptible amphibian, Atelopus zeteki, with different infection histories. We found significant changes in expression of numerous genes involved in innate and inflammatory responses in infected frogs despite high susceptibility to chytridiomycosis.

We show evidence of acquired immune responses generated against Bd, including increased expression of immunoglobulins and major histocompatibility complex genes. In addition, fungal-killing genes had significantly greater expression in frogs previously exposed to Bd compared with Bd-naïve frogs, including chitinase and serine-type proteases. However, our results appear to confirm recent in vitro evidence of immune suppression by Bd, demonstrated by decreased expression of lymphocyte genes in the spleen of infected compared with control frogs. We propose susceptibility to chytridiomycosis is not due to lack of Bd-specific immune responses but instead is caused by failure of those responses to be effective. Ineffective immune pathway activation and timing of antibody production are discussed as potential mechanisms. However, in light of our findings, suppression of key immune responses by Bd is likely an important factor in the lethality of this fungus.

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.

Amphibians vary in their response to infection by the amphibian-killing chytrid fungus, Batrachochytrium dendrobatidis (Bd). Highly susceptible species are the first to decline and/or disappear once Bd arrives at a site. These competent hosts likely facilitate Bd proliferation because of ineffective innate and/or acquired immune defenses. We show that Atelopus zeteki, a highly susceptible species that has undergone substantial population declines throughout its range, rapidly and exponentially increases skin Bd infection intensity, achieving intensities that are several orders of magnitude greater than most other species reported. We experimentally infected individuals that were never exposed to Bd (n = 5) or previously exposed to an attenuated Bd strain (JEL427-P39; n = 3). Within seven days post-inoculation, the average Bd infection intensity was 18,213 zoospores (SE: 9,010; range: 0 to 66,928).

Both average Bd infection intensity and zoospore output (i.e., the number of zoospores released per minute by an infected individual) increased exponentially until time of death (t₅₀ = 7.018, p<0.001, t₄₆ = 3.164, p = 0.001, respectively). Mean Bd infection intensity and zoospore output at death were 4,334,422 zoospores (SE: 1,236,431) and 23.55 zoospores per minute (SE: 22.78), respectively, with as many as 9,584,158 zoospores on a single individual. The daily percent increases in Bd infection intensity and zoospore output were 35.4% (SE: 0.05) and 13.1% (SE: 0.04), respectively. We also found that Bd infection intensity and zoospore output were positively correlated (t₄₃ = 3.926, p<0.001). All animals died between 22 and 33 days post-inoculation (mean: 28.88; SE: 1.58). Prior Bd infection had no effect on survival, Bd infection intensity, or zoospore output. We conclude that A. zeteki, a highly susceptible amphibian species, may be an acute supershedder. Our results can inform epidemiological models to estimate Bd outbreak probability, especially as they relate to reintroduction programs.

Laboratory investigations into the amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), have accelerated recently, given the pathogen’s role in causing the global decline and extinction of amphibians. Studies in which host animals were exposed to Bd have largely assumed that lab-maintained pathogen cultures retained the infective and pathogenic properties of wild isolates. Attenuated pathogenicity is common in artificially maintained cultures of other pathogenic fungi, but to date, it is unknown whether, and to what degree, Bd might change in culture. We compared zoospore production over time in two samples of a single Bd isolate having different passage histories: one maintained in artificial media for more than six years (JEL427-P39), and one recently thawed from cryopreserved stock (JEL427-P9). In a common garden experiment, we then exposed two different amphibian species, Eleutherodactylus coqui and Atelopus zeteki, to both cultures to test whether Bd attenuates in pathogenicity with in vitro passages. The culture with the shorter passage history, JEL427-P9, had significantly greater zoospore densities over time compared to JEL427-P39. This difference in zoospore production was associated with a difference in pathogenicity for a susceptible amphibian species, indicating that fecundity may be an important virulence factor for Bd. In the 130-day experiment, Atelopus zeteki frogs exposed to the JEL427-P9 culture experienced higher average infection intensity and 100% mortality, compared with 60% mortality for frogs exposed to JEL427-P39. This effect was not observed with Eleutherodactylus coqui, which was able to clear infection. We hypothesize that the differences in phenotypic performance observed with Atelopus zeteki are rooted in changes of the Bd genome. Future investigations enabled by this study will focus on the underlying mechanisms of Bd pathogenicity.

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.

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.

Infection after renal transplantation remains a major cause of morbidity and death, especially infection from the extensively drug-resistant bacteria, A. baumannii. A total of fourteen A. baumannii isolates were isolated from the donors’ preserved fluid from DCD (donation after cardiac death) renal transplantation and four isolates in the recipients’ draining liquid at the Kidney Disease Center, The First Affiliated Hospital, College of Medicine, Zhejiang University, from March 2013 to November 2014. An outbreak of A. baumannii emerging after DCD renal transplantation was tracked to understand the transmission of the pathogen. PFGE displayed similar DNA patterns between isolates from the same hospital. Antimicrobial susceptibility tests against thirteen antimicrobial agents were determined using the K-B diffusion method and eTest. Whole-genome sequencing was applied to investigate the genetic relationship of the isolates. With the clinical data and research results, we concluded that the A. baumannii isolates 3R1 and 3R2 was probably transmitted from the donor who acquired the bacteria during his stay in the ICU, while isolate 4R1 was transmitted from 3R1 and 3R2 via medical manipulation. This study demonstrated the value of integration of clinical profiles with molecular methods in outbreak investigation and their importance in controlling infection and preventing serious complications after DCD transplantation.

Competing endogenous RNAs (ceRNAs) are RNA molecules that sequester shared microRNAs (miRNAs) thereby affecting the expression of other targets of the miRNAs. Whether genetic variants in ceRNA can affect its biological function and disease development is still an open question. Here we identified a large number of genetic variants that are associated with ceRNA's function using Geuvaids RNA-seq data for 462 individuals from the 1000 Genomes Project. We call these loci competing endogenous RNA expression quantitative trait loci or ‘cerQTL’, and found that a large number of them were unexplored in conventional eQTL mapping. We identified many cerQTLs that have undergone recent positive selection in different human populations, and showed that single nucleotide polymorphisms in gene 3΄UTRs at the miRNA seed binding regions can simultaneously regulate gene expression changes in both cis and trans by the ceRNA mechanism. We also discovered that cerQTLs are significantly enriched in traits/diseases associated variants reported from genome-wide association studies in the miRNA binding sites, suggesting that disease susceptibilities could be attributed to ceRNA regulation. Further in vitro functional experiments demonstrated that a cerQTL rs11540855 can regulate ceRNA function. These results provide a comprehensive catalog of functional non-coding regulatory variants that may be responsible for ceRNA crosstalk at the post-transcriptional level.