Zhang, Yuji

ABSTRACT While BAG3 has been identified as a causative gene for dilated cardiomyopathy, the major pathological events in BAG3‐ related cardiomyopathy that could be targeted for therapeutic benefit remain to be discovered. Here, we aim to uncover novel pathological events through genetic studies in a zebrafish bag3 cardiomyopathy model. Given the known cardioprotective effects of mtor inhibition and the fact that transcription factor EB (tfeb ) encodes a direct downstream phosphorylation target of mTOR signaling, we generated a cardiomyocyte‐specific transgenic line overexpressing tfeb ( Tg[cmlc2:tfeb] ). This overexpression was sufficient to restore defective proteostasis and rescue cardiac dysfunction in the bag3 cardiomyopathy model. Importantly, we detected accelerated cardiac senescence in the bag3 cardiomyopathy model, which can be mitigated by Tg(cmlc2:tfeb) . We compared cardiac transcriptomes between the Tg(cmlc2:tfeb) transgenic fish and the mtor xu015/+ mutant and found that inhibition of the fatty acid binding protein a ( fabp7a ) gene exerts therapeutic effects. Consistent with this genetic evidence, we detected elevated fabp7a expression in the bag3 cardiomyopathy model, whereas cardiomyocyte‐specific overexpression of fabp7a induced dysregulated proteostasis, accelerated cardiac senescence, and cardiac dysfunction. To elucidate the functions of Fabp7a in normative cardiac aging, we turned to the African Turquoise Killifish. We noted elevated Fabp7a expression in the hearts of aged killifish, and pharmacological inhibition of Fabp7a mitigated the cardiac aging process. Together, this study uncovered accelerated cardiac senescence as a key pathological event in bag3 cardiomyopathy and reveals that manipulating the mTOR‐Tfeb‐Fabp7a axis can mitigate this pathology and confer cardioprotective effects.

Abstract To test the hypothesis that complex immune phenotypes can be defined by stable gene expression profiles we are using a RNAseq-based transcriptome analysis to relate phenotypes of activated peripheral blood CD4+ cells from healthy subjects and stage IV melanoma patients to their freshly-isolated counterparts. We have defined molecular immune response signatures (MIRS) containing genes that consistently differentiate individuals from one another. Two different phenotypic groups are defined using this approach, both of which are mathematically distinct from randomly-selected groups. First, a MIRS generated using healthy donor T cells stimulated with α-CD3/CD28 beads defines a distinct response phenotype and successfully predicts 17 independent healthy subject and melanoma patient responses, demonstrating the predictive capability of the approach. Second, our comparison of melanoma patient signatures with those of healthy subjects revealed a group of healthy individuals with distinct expression signatures that did not cluster with those of any of the melanoma patients. These results indicate that our RNAseq-based transcriptome analysis has the potential to identify phenotypic groups as well as gene sets that have response-predictive capabilities. Our work has both discovery and translational implications, as our approach can be used to investigate complex genetic contributions to immune responses and predict how individual patients will respond to therapeutic interventions.

Previously we showed the generation of a protein trap library made with the gene-break transposon (GBT) in zebrafish ( Danio rerio ) that could be used to facilitate novel functional genome annotation towards understanding molecular underpinnings of human diseases (Ichino et al, 2020). Here, we report a significant application of this library for discovering essential genes for heart rhythm disorders such as sick sinus syndrome (SSS). SSS is a group of heart rhythm disorders caused by malfunction of the sinus node, the heart’s primary pacemaker. Partially owing to its aging-associated phenotypic manifestation and low expressivity, molecular mechanisms of SSS remain difficult to decipher. From 609 GBT lines screened, we generated a collection of 35 zebrafish insertional cardiac (ZIC) mutants in which each mutant traps a gene with cardiac expression. We further employed electrocardiographic measurements to screen these 35 ZIC lines and identified three GBT mutants with SSS-like phenotypes. More detailed functional studies on one of the arrhythmogenic mutants, GBT411 , in both zebrafish and mouse models unveiled Dnajb6 as a novel SSS causative gene with a unique expression pattern within the subpopulation of sinus node pacemaker cells that partially overlaps with the expression of hyperpolarization activated cyclic nucleotide gated channel 4 (HCN4), supporting heterogeneity of the cardiac pacemaker cells.

Myosin heavy chain 7 (MYH7) is a major causative gene for hypertrophic cardiomyopathy, but the affected signaling pathways and therapeutics remain elusive. In this research, we identified ventricle myosin heavy chain like ( vmhcl ) as a zebrafish homolog of human MYH7 , and we generated vmhcl frameshift mutants. We noted vmhcl -based embryonic cardiac dysfunction (VEC) in the vmhcl homozygous mutants and vmhcl -based adult cardiomyopathy (VAC) phenotypes in the vmhcl heterozygous mutants. Using the VEC model, we assessed 7 known cardiomyopathy signaling pathways pharmacologically and 11 candidate genes genetically via CRISPR/Cas9 genome editing technology based on microhomology-mediated end joining (MMEJ). Both studies converged on therapeutic benefits of mTOR or mitogen-activated protein kinase (MAPK) inhibition of VEC. While mTOR inhibition rescued the enlarged nuclear size of cardiomyocytes, MAPK inhibition restored the prolonged cell shape in the VEC model. The therapeutic effects of mTOR and MAPK inhibition were later validated in the VAC model. Together, vmhcl/myh7 loss of function is sufficient to induce cardiomyopathy in zebrafish. The VEC and VAC models in zebrafish are amenable to both efficient genetic and chemical genetic tools, offering a rapid in vivo platform for discovering candidate signaling pathways of MYH7 cardiomyopathy.

Background: Rapid advance of genome technologies accelerate the discovery of genetic basis of cardiomyopathy and heart failure and also enable system biology studies to pinpoint underlying mechanism. However, the limited throughput of mammalian models restricted the number of genes that can be studied in a particular lab. Adult zebrafish has been recently pursued as a new model with higher throughput. However, as a non-mammalian model, its conservation is not tested. Objective: To assess the conservation of zebrafish for genetic studies of human dilated cardiomyopathy (DCM) via transcriptome analysis of 51 known DCM-causative genes. Methods and Results: By conducting RNA-sequencing (RNA-seq) analysis of larva and adult zebrafish, we identified genes with high expression level in the heart and fetal gene program using differential expression between embryonic and adult stages. We then searched zebrafish orthologues for 51 reported human DCM-causative or associated genes and identified zebrafish orthologues for 49 of them. While 30 genes have a single orthologue, 14 genes have two homologues and the remaining 5 genes have more than three. We then applied the transcriptome data to prioritize these homologues for the 19 DCM causative genes with more than one homologue. Based on the cardiac abundance and cardiac enrichment hypothesis, we are able to recommend a single zebrafish homologue of high priority for 12 out of 19 DCM genes, 2 zebrafish homologues of high priority for ACTC1. Interestingly, our expressional data suggested zebrafish othologues for human MYH6 and MYH7 , respectively. Similar to that in mammals, these two zebrafish othologues are oppositely expressed during zebrafish embryonic and adult stage. Conclusions: Orthologues for the majority of DCM causative genes can be found in Zebrafish, supporting its usage as a conserved vertebrate model for studying DCM. The definition of cardiac transcriptome in zebrafish will facilitate the future system biology studies. This vertebrate model with higher throughput can be further leveraged to validate the novel variants identified from human patients, to understand underlying signaling pathways and to develop novel therapeutics.

Abstract Introduction We have previously developed a primary human bone marrow model in CLL which allows us to study marrow derived mesenchymal stroma cells (MSC) and mimic in vivoMSC-leukemic interactions in CLL (Ding, Blood, 2010) These interactions proved to be bidirectional with functional implications for both cell types. However the biologic and clinical implications of these bilateral cellular interactions in relation to MSC function in CLL patients have not been completely studied. MSC are clearly established as an essential component in forming supportive marrow niches for hematopoiesis. Our focus in this study thus was to further investigate on the status of CLL MSC function in comparison with normal MSC and the impact of the leukemic CLL B cell on MSC function and its role in forming marrow niches in CLL. Methods Bone marrow (BM) aspirates of CLL patients or normal age-matched subjects were used to generate MSC and these MSC were then assessed for various functions. Specifically we assessed the growth, differentiation capacity, gene expression, cytokine secretion profile and transcriptome analysis of CLL and normal MSC using colony forming unit (CFU) assay, MSC differentiation assay, Affymetrix U133 2.0 plus array and RNA sequencing respectively. To determine MSC CFU frequency, fixed number of nucleated cells (5 or 10 million) from BM aspirates were added in MSC growth medium for 10 days to form CFU. CFU frequency was then calculated by dividing the number of MSC colonies by the fixed number of inputted nucleated cells. In order to test the specific gene changes that occurred in MSC after their interaction with CLL leukemic cells, normal MSC were co-cultured with CLL cells for 1 week and subsequently followed by cell growth, gene expression analysis with RNA sequencing and cytokine profiles using a multiplex assay. Results CLL MSC generated significantly lower CFU colonies compared to normal subjects (CLL [n=24] vs normal [n=6] mean: 2.2 vs. 22 CFU per million nucleated cells, p=0.01). In order to test if the lower number of CLL MSC CFU was due to their reduced growth capacity or lower frequency, we then performed the CFU assay using identical numbers of CLL and normal MSC at the start of culture. CLL MSC was found to have significantly lower CFU capacity compared to normal MSC (mean: 2.1 [CLL] vs 4.5 [normal] CFU/10 input cell, n=7 in each group, p=0.03). To investigate the mechanism(s) for the reduced growth capacity of CLL MSC, we compared the gene expression profiles of non-stimulated CLL vs. normal MSC (n=10 for each group) with the following outcomes: CLL MSC had increased expression of cell cycle inhibitor: p16(Ink4a) and p57, both key markers of cell senescence (Baker, Nature, 2011); CLL MSC overexpressed homeobox B cluster of transcription factors, key regulator for cell development and differentiation. Wnt inhibitors DKK1/DKK2 and non-canonical Wnt 5b known to be associated with aging were found to be 10 and 2 fold more elevated in CLL MSC. In addition, MSC abundantly produced a milieu of the cytokines MCP-1/IL-8/IL-6/IL-1Ra after interaction with CLL leukemic cells. These cytokines are associated with a senescent associated secretory phenotype (SASP). Co-culture experiments of normal MSC with CLL Leukemic cells revealed significant impact of leukemic cells on normal MSC functions. Co-culture with a low density of CLL leukemic cells (CLL to MSC ratio: ≤10,000/1) increased MSC growth (1.2 fold). However, high density co-culture of CLL cells with MSC (CLL to MSC ratio of 100,000/1) inhibited MSC growth (0.8 fold). Transcriptome analysis of normal MSC before and after one week of co-culture with CLL leukemic cells showed that normal MSC had alteration of ~1500 genes. Many of the modified genes (CDKN2B, DKK2, LIF, HGF, FOXQ1 etc) are known to be involved in regulation of cell growth and senescence. Conclusion Collectively, these results demonstrated that bone marrow derived CLL MSC have decreased proliferative potential and possess senescent phenotypes. Importantly these senescent features were inducible in normal MSC after their interaction with CLL leukemic cells implying that leukemic B cells can educate MSC to become senescent. The implications of the CLL MSC senescence status are multiple and include: defective hematopoietic niches due to dysfunctional MSC, immune dysregulation and the enhanced nurturing of CLL B cell by these MSC. We are now investigating several of these latter features in our current work. Disclosures Shanafelt: Genentech: Research Funding; GlaxoSmithKline: Research Funding; Celegene: Research Funding; Cephalon: Research Funding; Pharmacyclics/Jannsen: Research Funding; Hospira: Research Funding; Polyphenon E Int'l: Research Funding. Kay:Genetech: Research Funding; Pharmacyclics: Research Funding; Hospira: Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees; Gilead: Membership on an entity's Board of Directors or advisory committees.

ABSTRACT The adult zebrafish is an emerging vertebrate model for studying human cardiomyopathies; however, whether the simple zebrafish heart can model different subtypes of cardiomyopathies, such as dilated cardiomyopathy (DCM), remains elusive. Here, we generated and characterized an inherited DCM model in adult zebrafish and used this model to search for therapeutic strategies. We employed transcription activator-like effector nuclease (TALEN) genome editing technology to generate frame-shift mutants for the zebrafish ortholog of human BCL2-associated athanogene 3 (BAG3), an established DCM-causative gene. As in mammals, the zebrafish bag3 homozygous mutant (bag3e2/e2) exhibited aberrant proteostasis, as indicated by impaired autophagy flux and elevated ubiquitinated protein aggregation. Through comprehensive phenotyping analysis of the mutant, we identified phenotypic traits that resembled DCM phenotypes in mammals, including cardiac chamber enlargement, reduced ejection fraction characterized by increased end-systolic volume/body weight (ESV/BW), and reduced contractile myofibril activation kinetics. Nonbiased transcriptome analysis identified the hyperactivation of the mechanistic target of rapamycin (mTOR) signaling in bag3e2/e2 mutant hearts. Further genetic studies showed that mtorxu015/+, an mTOR haploinsufficiency mutant, repaired abnormal proteostasis, improved cardiac function and rescued the survival of the bag3e2/e2 mutant. This study established the bag3e2/e2 mutant as a DCM model in adult zebrafish and suggested mtor as a candidate therapeutic target gene for BAG3 cardiomyopathy.

Abstract Genome-wide identification and characterization of long noncoding RNAs (lncRNA) in individual immune cell lineages helps us better understand the driving mechanisms behind melanoma and advance personalized patient treatment. To elucidate the transcriptional landscape in diverse immune cell types of peripheral blood cells (PBC) in stage IV melanoma, we used whole transcriptome RNA sequencing to profile lncRNAs in CD4+, CD8+, and CD14+ PBC from 132 patient samples. Our integrative computational approach identified 27,625 expressed lncRNAs, 2,744 of which were novel. Both T cells (i.e., CD4+ and CD8+ PBC) and monocytes (i.e., CD14+ PBC) exhibited differential transcriptional expression profiles between patients with melanoma and healthy subjects. Cis- and trans-level coexpression analysis suggested that lncRNAs are potentially involved in many important immune-related pathways and the programmed cell death receptor 1 checkpoint pathways. We also identified nine gene coexpression modules significantly associated with melanoma status, all of which were significantly enriched for three mRNA translation processes. Age and melanoma traits closely correlated with each other, implying that melanoma contains age-associated immune changes. Our computational prediction analysis suggests that many cis- and trans-regulatory lncRNAs could interact with multiple transcriptional and posttranscriptional regulatory elements in CD4+, CD8+, and CD14+ PBC, respectively. These results provide novel insights into the regulatory mechanisms involving lncRNAs in individual immune cell types in melanoma and can help expedite cell type-specific immunotherapy treatments for such diseases. Significance: These findings elucidate melanoma-associated changes to the noncoding transcriptional landscape of distinct immune cell classes, thus providing cell type-specific guidance to targeted immunotherapy regimens. Cancer Res; 78(15); 4411–23. ©2018 AACR.

Abstract Human immunity exhibits remarkable heterogeneity among individuals, which engenders variable responses to immune perturbations in human populations. Population studies reveal that, in addition to interindividual heterogeneity, systemic immune signatures display longitudinal stability within individuals, and these signatures may reliably dictate how given individuals respond to immune perturbations. We hypothesize that analyzing relationships among these signatures at the population level may uncover baseline immune phenotypes that correspond with response outcomes to immune stimuli. To test this, we quantified global gene expression in peripheral blood CD4+ cells from healthy individuals at baseline and following CD3/CD28 stimulation at two time points 1 mo apart. Systemic CD4+ cell baseline and poststimulation molecular immune response signatures (MIRS) were defined by identifying genes expressed at levels that were stable between time points within individuals and differential among individuals in each state. Iterative differential gene expression analyses between all possible phenotypic groupings of at least three individuals using the baseline and stimulated MIRS gene sets revealed shared baseline and response phenotypic groupings, indicating the baseline MIRS contained determinants of immune responsiveness. Furthermore, significant numbers of shared phenotype-defining sets of determinants were identified in baseline data across independent healthy cohorts. Combining the cohorts and repeating the analyses resulted in identification of over 6000 baseline immune phenotypic groups, implying that the MIRS concept may be useful in many immune perturbation contexts. These findings demonstrate that patterns in complex gene expression variability can be used to define immune phenotypes and discover determinants of immune responsiveness.